当 Rasterstacks 在一个未知长度的列表中时求和 RasterStacks

Question

我有一个 Rasterstacks 列表，每个都有 2 层。我想将列表中的所有 Rasterstack 添加到具有 2 层的单个 Rasterstack 中。假设我有 2 个 Rasterstacks 的列表，每个都有 2 层，我将使用以下命令

raster_stack_sum = raster_list[[1]] + raster_list[[2]]

这将返回具有 2 层的单个 Rasterstack。我的问题是当列表的长度（即 RasterStacks 的数量）未知时，我将如何做到这一点。当列表的长度（在本例中为 5）未知时，是否有等效的更短方法来运行以下代码行？

raster_stack_sum = raster_list[[1]] + raster_list[[2]] + 
raster_list[[3]]+ raster_list[[4]] + raster_list[[5]]

Answer 1

您可以使用stack 从Raster* 对象列表中创建一个新的RasterStack。

示例数据：

library(raster)
s <- stack(system.file("external/rlogo.grd", package="raster")) 
x <- list(s, s, s)

解决方案：

sum(stack(x))

是什么让您说“不推荐存储 Rasterstacks 列表”？

编辑

要逐层汇总 RasterStack 对象列表，我会使用循环

y <- x[[1]]
for (i in 2:length(x)) { y = y + x[[i]] }