堆積光柵在循環中

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 )
}