Я хочу завантажити .adf файли в R. Дані з цієї сторінки: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057
Я спробував наступний код, який знайшов після деяких досліджень в Інтернеті. Проблема полягає в тому, що в класі RasterLayer я отримую негативні значення, яких там не повинно бути. Я не знаю, чому це відбувається, тому, сподіваюся, хтось може мені допомогти !?
Код:
library(rgdal)
library(RColorBrewer)
dpath<- path...
x <- new("GDALReadOnlyDataset", dpath)
getDriver(x)
getDriverLongName(getDriver(x))
xx<-asSGDF_GROD(x)
r <- raster(xx)
Вихід для 'r':
r клас: Розміри RasterLayer: 2160, 4320, 9331200 (nrow, ncol, ncell) роздільна здатність: 0,08333333, 0,08333333 (x, y) ступінь: -180, 180, -90, 90 (xmin, xmax, ymin, ymax) координата. реф. : + proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + no_defs Джерело даних: в іменах пам'яті: band1 значення: -997, 16 (min, max)
Значення "16" відноситься до 16 класів тривалості періоду вегетації. Але мені цікаво, звідки беруться ті '-997'. Можливо, щось не так з координатою. ref?
Ось також резюме даних "xx":
Підсумок даних: Хв. 1-е Кв. Середнє середнє 3-е ква. Макс. NA -997 3 5 -9 8 16 7123158
А якщо ми уважніше розглянемо дані в xx:
таблиця (xx $ band1)
-997 1 2 3 4 5 6 7 8 9 10 11 12 31711 429643 83011 166674 207228 270161 240958 183342 118608 98795 88473 73743 56022 13 14 15 16 30104 45521 52216 31832
У цьому дійсно є лише ця річ «-997». Я думаю, що НА є озеанами, то чи щось не так з завантаженням даних чи я просто не розумію дані?