在堆栈中复制栅格

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)