Я хочу загрузить файлы .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 (мин, макс)
«16» в значениях относится к 16 классам длины вегетационного периода. Но мне интересно, откуда взялись эти «-997». Может быть, что-то не так с координа. исх?
Вот также сводка данных 'xx':
Сводка данных: мин. 1 кв. Медиана Среднее 3 кв. Максимум. АН -997 3 5 -9 8 16 7123158
И если мы посмотрим на данные в хх более внимательно:
таблица (хе $ 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. Я думаю, что АН - это ozeans, так что-то не так с загрузкой данных, или я просто не понимаю данные?
источник
Ответы:
Вы почти правы
NODATA установлен на -32768 для океанов. Кроме того, -997 установлен для великих озер, которые не исключены береговой линией.
Поскольку пиксельное содержание (период роста) не имеет смысла на озерах, вы можете спокойно относиться к -997 как к NODATA.
источник