组合栅格堆栈

时间:2016-04-08 08:29:56

标签: r raster r-raster

我有多个光栅堆栈,其中包含来自CRU TS 3.23的温度数据

我在文件中读到如下:

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")

有什么方法可以将所有这些栅格组合起来创建一个大栅格?我正在研究降水趋势,当它们在一个物体中时更容易做到!感谢

1 个答案:

答案 0 :(得分:2)

raster::stack(nc, nc2, nc3, nc4, nc5, nc6, nc7)应该做到这一点。已经描述了非常类似的方法,例如, here

但是,请注意,您希望stack的对象需要具有相同的坐标参考系统(CRS),分辨率(res)和{{1} }。否则,extent很可能会失败。如果上述条件中的一个(或几个)不同,您可以部署stackraster::projectRaster来创建具有统一空间属性的对象。