合并栅格堆栈

Question

合并栅格堆栈

4

我有多个栅格堆栈，包含来自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")

有没有办法将所有这些光栅图像合并为一个大图像？我正在研究降水趋势，如果它们都在一个对象中，那么做起来会更容易！谢谢。

- A_J_M

1个回答

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- fdetsch · Accepted Answer

raster::stack(nc, nc2, nc3, nc4, nc5, nc6, nc7)可以实现此功能。类似的方法已经在这里进行了描述。

但是，请注意，您想要stack的对象需要具有相同的坐标参考系统（CRS）、分辨率（res）和extent。否则，stack很可能会失败。如果上述任何一个条件不同，您可以使用raster::projectRaster和raster::resample来创建具有统一空间属性的对象。