2016-04-08 10 views
2

ich mehrere Raster-Stacks enthalten Temperaturdaten von CRU TS 3.23Kombination Raster stapelt

Ich habe habe in den Dateien wie folgt lauten:

nc = stack("cru_ts3.23.1951.1960.pre.dat.nc", varname = "pre") 
nc2 = stack("cru_ts3.23.1961.1970.pre.dat.nc", varname = "pre") 
... 
nc7 = stack("cru_ts3.23.2011.2014.pre.dat.nc", varname = "pre") 

Gibt es eine Möglichkeit, all diese Rastern kombinieren zu schaffen ein großer? Ich schaue auf Niederschlagstrends, die viel einfacher sind, wenn sie alle in einem Objekt sind! Danke

Antwort

2

raster::stack(nc, nc2, nc3, nc4, nc5, nc6, nc7) sollte den Trick tun. Ein sehr ähnlicher Ansatz wurde z. here.

Beachten Sie jedoch, dass die Objekte, die Sie ein identisches Koordinaten Referenzen System (CRS) zu stack benötigt haben möchten, Auflösung (res) und extent. Andernfalls wird wahrscheinlich stack fehlschlagen. Sie könnten raster::projectRaster und raster::resample bereitstellen, um Objekte mit einheitlichen räumlichen Eigenschaften zu erstellen, wenn eines (oder mehrere) der obigen Kriterien unterschiedlich ist.

+1

Sie können das Stapeln auch in einer Zeile durchführen mit: 'stack (" cru_ts3.23.1951.1960.pre.dat.nc "," cru_ts3.23.1951.1970.pre.dat.nc "," cru_ts3.23.1951 .2014.pre.dat.nc ", varname =" pre ")' – maRtin

+1

Wahre Geschichte @maRtin, ich wollte nur bei den oben erstellten Objekten bleiben. – fdetsch