堆积光栅在循环中

Question

我需要在循环中堆叠一些栅格，如：

for(month in 1:12){
.
.
.
"some algorithm spiting out a raster called 'sm_esa'"
sm_esa_stack<-stack(sm_esa)
}

最后，我想创建一个包含12层（每层一个月）的堆栈。但是我的最后一行显然会覆盖每一个新的栅格而不是叠加。任何提示？

Answer 1

在循环外实例化一个空堆栈，在循环的每次迭代之后，通过堆叠当前堆栈和新的Rasterlayer将新的rasterLayer添加到堆栈。

x <- stack()
for(month in 1:12){
.
.
.
"some algorithm spiting out a raster called 'sm_esa'"
x <- stack( x , sm_esa )
}