Wie lese ich * .adf-Dateien in R?

9

Ich möchte ADF-Dateien in R laden. Die Daten stammen von dieser Seite: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057

Ich habe den folgenden Code ausprobiert, den ich nach einigen Recherchen im Internet gefunden habe. Das Problem ist, dass ich in der Klasse RasterLayer negative Werte bekomme, die nicht vorhanden sein sollten. Ich weiß nicht, warum das passiert, also kann mir hoffentlich jemand helfen!?

Code:

library(rgdal)
library(RColorBrewer)
dpath<- path...

x <- new("GDALReadOnlyDataset", dpath)
getDriver(x)
getDriverLongName(getDriver(x))
xx<-asSGDF_GROD(x)
r <- raster(xx)

Die Ausgabe für 'r' lautet:

r Klasse: RasterLayer Abmessungen: 2160, 4320, 9331200 (nrow, ncol, ncell) Auflösung: 0,08333333, 0,08333333 (x, y) Ausdehnung: -180, 180, -90, 90 (xmin, xmax, ymin, ymax) Koordinate. ref. : + proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + no_defs Datenquelle: in Speichernamen: band1 Werte: -997, 16 (min, max)

Die '16' in den Werten bezieht sich auf die 16 Klassen der Länge der Wachstumsperiode. Aber ich frage mich, woher diese '-997' kommen. Vielleicht stimmt etwas nicht mit der Koordinate. ref?

Hier ist auch eine Datenzusammenfassung von 'xx':

Datenübersicht: min. 1. Qu. Median Mean 3rd Qu. Max. NA's -997 3 5 -9 8 16 7123158

Und wenn wir uns die Daten in xx genauer ansehen:

Tabelle (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

Da ist wirklich nur dieses '-997'-Ding drin. Ich denke, die NAs sind Ozeane. Stimmt etwas mit dem Laden der Daten nicht oder verstehe ich die Daten einfach nicht?

Franz Krah
quelle
1
Bitte ändern Sie den Titel in etwas Intelligentes und Genaues.
Mdsumner
fwiw, es ist nicht nötig, rgdal direkt zu verwenden, nur Raster (dpath)
mdsumner

Antworten:

3

Sie haben fast recht:

NODATA ist für Ozeane auf -32768 eingestellt. Darüber hinaus ist -997 für große Seen festgelegt, die von der Küste nicht ausgeschlossen sind.

Da der Pixelgehalt (Wachstumsphase) auf Seen keinen Sinn ergibt, können Sie -997 auch sicher als NODATA behandeln.

AndreJ
quelle