在堆栈中复制栅格

Question

我想将每3天可用的土壤湿度栅格与每日可用的降雨栅格进行比较。我制作每个堆栈并重新采样到适当的分辨率。现在，为了轻松地比较堆栈，能够复制土壤湿度堆栈中的每一层并将其插入自身旁边两次是很好的。这与Stacking an existing RasterStack multiple times基本相同 除了我需要对大堆进行排序，以便所有栅格都按时间顺序排列。有没有办法做到这一点？

（我知道我可以在第一次堆叠文件之前复制文件，但这需要重新采样3倍的堆栈。由于重新采样是我脚本中最慢的一部分，应该有更好的方法。）

Answer 1

这样的东西？

# example data
r <- raster(ncol=10, nrow=10)
r[]=1:ncell(r)
x <- brick(r,r,r,r,r,r)
x <- x * 1:6

y <- list()
for (i in 1:nlayers(x)) {
    r <- raster(x, i)
    y <- c(y, r, r, r)
}
s <- stack(y)