Як читати * .adf файли в R?


9

Я хочу завантажити .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». Я думаю, що НА є озеанами, то чи щось не так з завантаженням даних чи я просто не розумію дані?


1
Будь ласка, змініть заголовок на щось розумне та точне.
mdsumner

fwiw, взагалі не потрібно використовувати rgdal, просто растр (dpath)
mdsumner

Відповіді:


3

Ви майже маєте рацію:

Для океанів NODATA встановлено значення -32768. Крім того, -997 встановлений для великих озер, які не виключаються береговою лінією.

Оскільки вміст пікселів (період вирощування) не має сенсу для озер, ви можете сміливо ставитися до -997 як до NODATA.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.