如何用Terra批量合并单波段栅格为多波段栅格并解决循环报错？

错误原因 #

你原来的代码之所以失败，是因为用普通向量存储terra::rast()生成的栅格对象——普通向量没法容纳这类复杂的空间对象，赋值时自然会报错，最后向量里全是NA。

正确实现方法 #

方法一：列表存储+批量合并 #

先把所有栅格读到列表里，再合并成多波段栅格：

# 加载terra包即可，无需加载raster包
library(terra)

# 获取工作目录下所有tif文件，full.names确保拿到完整路径
raster_paths <- list.files(pattern = "\\.tif$", full.names = TRUE)

# 批量读取栅格到列表
raster_list <- lapply(raster_paths, rast)

# 合并列表里的栅格为多波段栅格
multi_band_raster <- do.call(c, raster_list)

方法二：一步到位直接读取 #

terra::rast()支持直接传入路径向量，一步生成多波段栅格，更省事：

library(terra)

raster_paths <- list.files(pattern = "\\.tif$", full.names = TRUE)
multi_band_raster <- rast(raster_paths)

生成GIF动画 #

合并完成后，直接用animate()处理多波段栅格就能生成动画，还能指定保存路径：

# 生成并保存GIF，可调整宽高、帧率等参数
animate(multi_band_raster, filename = "栅格动画.gif", width = 800, height = 600)

注意点 #

• 所有待合并的栅格必须空间范围、分辨率、投影完全一致，不然合并会失败。
• pattern = "\\.tif$"是精准匹配后缀为tif的文件，避免误匹配其他带“tif”字符的文件。
• 一定要加full.names = TRUE，不然可能因路径问题找不到文件。

内容的提问来源于stack exchange，提问作者89_Simple