Quiero cargar archivos .adf en R. Los datos provienen de esta página: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057
Intenté el siguiente código que encontré después de algunas investigaciones en Internet. El problema es que en la clase RasterLayer obtengo valores negativos que no deberían estar allí. ¡No sé por qué sucede esto, así que espero que alguien pueda ayudarme!
Código:
library(rgdal)
library(RColorBrewer)
dpath<- path...
x <- new("GDALReadOnlyDataset", dpath)
getDriver(x)
getDriverLongName(getDriver(x))
xx<-asSGDF_GROD(x)
r <- raster(xx)La salida para 'r' es:
Clase r: Dimensiones RasterLayer: 2160, 4320, 9331200 (nrow, ncol, ncell) resolución: 0.08333333, 0.08333333 (x, y) extensión: -180, 180, -90, 90 (xmin, xmax, ymin, ymax) coord. árbitro. : + proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + fuente de datos no_defs: en nombres de memoria: valores de banda1: -997, 16 (min, max)
El '16' en los valores se refiere a las 16 clases de duración del período de crecimiento. Pero me pregunto de dónde provienen esos '-997'. Tal vez algo anda mal con el coord. ¿árbitro?
Aquí también hay un resumen de datos de 'xx':
Resumen de datos: mín. 1st Qu. Mediana media 3er Qu. Max. NA -997 3 5 -9 8 16 7123158
Y si miramos los datos en xx más de cerca:
tabla (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
Realmente hay solo esta cosa '-997' en ella. Creo que las NA son océanos, ¿hay algo mal con la carga de datos o simplemente no entiendo los datos?
fuente
Respuestas:
Tienes casi razón:
NODATA se establece en -32768 para los océanos. Además, -997 está configurado para grandes lagos que no están excluidos de la costa.
Dado que el contenido de píxeles (período de crecimiento) no tiene sentido en los lagos, también puede tratar con seguridad -997 como NODATA.
fuente