重复相同的栅格图层以创建栅格堆栈

Question

我正在尝试从栅格层创建栅格堆栈，其中栅格堆栈就是重复一定次数的同一栅格层。

我可以这样做：

library(raster)
rasterstack <- addLayer(rasterlayer, rasterLayer, rasterLayer) 

这是可行的。但是，我希望堆栈大约有1000层。我想我可以循环通过，但我想知道是否有更复杂的方法来做到这一点。

我尝试这样做的原因是为了计算栅格堆叠的加权平均值，其中每个层都是不同的时间段，并且权重位于不同的栅格层对象中。我希望，如果我从权重栅格层创建一个与数据具有相同层数的栅格堆栈，我将能够执行以下操作：

  weightedmean <- weighted.mean( data.RasterStack, weights.RasterStack )

Answer 1

示例数据

library(raster)
r <- raster(ncol=10, nrow=10, vals=1:100)

解决方案

n <- 10  # number of copies
s <- stack(lapply(1:n, function(i) r)) 

或

s <- stack(replicate(n, r))