Voglio caricare i file .adf in R. I dati provengono da questa pagina: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057
Ho provato il seguente codice che ho trovato dopo alcune ricerche su Internet. Il problema è che nella classe RasterLayer ottengo valori negativi che non dovrebbero esserci. Non so perché questo accada, quindi spero che qualcuno mi possa aiutare !?
Codice:
library(rgdal)
library(RColorBrewer)
dpath<- path...
x <- new("GDALReadOnlyDataset", dpath)
getDriver(x)
getDriverLongName(getDriver(x))
xx<-asSGDF_GROD(x)
r <- raster(xx)L'output per 'r' è:
classe r: dimensioni dello strato raster: 2160, 4320, 9331200 (nrow, ncol, ncell) risoluzione: 0,08333333, 0,08333333 (x, y) estensione: -180, 180, -90, 90 (xmin, xmax, ymin, ymax) coord. rif. : + proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + origine dati no_defs: in nomi di memoria: valori band1: -997, 16 (min, max)
Il '16' nei valori si riferisce alle 16 classi di lunghezza del periodo di crescita. Ma mi chiedo da dove vengano quei '-997'. Forse qualcosa che non va nel coord. ref?
Ecco anche un riepilogo dei dati di 'xx':
Riepilogo dati: min. 1 ° Qu. Media mediana 3a Qu. Max. NA -997 3 5-9 8 16 7123158
E se guardiamo più da vicino i dati in xx:
tavolo (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
C'è davvero solo questa cosa '-997'. Penso che le NA siano ozean, quindi c'è qualcosa che non va nel caricamento dei dati o non capisco i dati?