在RStudio中使用do.call合并多波段卫星影像Raster Brick出错求助

这个错误的根源很明确：你给存储影像的列表list.data设置了文件名作为元素名称，而do.call()会把这些名称当作merge()函数的参数名传递，但raster::merge()并不识别这些自定义的参数名，导致函数认为必填的x参数缺失。

快速修复方案 #

只需要移除列表的命名即可，有两种方式：

1. 直接删除names(list.data) <- list.filenames这一行代码；
2. 如果需要保留文件名做其他用途，在调用do.call()前清空列表名称：

names(list.data) <- NULL
merged_scene <- do.call(merge, list.data)

更简洁的代码优化 #

另外，你可以用lapply()替代for循环来更高效地生成影像列表，同时直接获取完整文件路径避免手动拼接，让代码更简洁易读（还能适配不同系统的路径分隔符）：

library(raster)
folder <- file.path(getwd(), "Images-folder")
# full.names=TRUE直接返回完整文件路径，无需手动拼接
list.filenames <- list.files(folder, pattern = "\\.tif$", full.names = TRUE)
# 用lapply生成无命名的影像列表
list.data <- lapply(list.filenames, brick)
# 执行合并
merged_scene <- do.call(merge, list.data)

替代方案：用purrr::reduce()合并 #

如果你习惯tidyverse风格，也可以用purrr::reduce()来逐步合并列表中的影像，这种方式不需要担心列表命名的问题，代码可读性也更强：

library(raster)
library(purrr)
folder <- file.path(getwd(), "Images-folder")
list.filenames <- list.files(folder, pattern = "\\.tif$", full.names = TRUE)
list.data <- lapply(list.filenames, brick)
merged_scene <- reduce(list.data, merge)

为什么命名会导致错误？ #

当你给list.data的元素命名后，do.call(merge, list.data)实际上等价于执行：

merge("raster1.tif" = raster1, "raster2.tif" = raster2, ...)

而raster::merge()的参数定义是merge(x, y, ..., tol=0, filename="", overwrite=FALSE, ...)，它只识别x、y、tol这些官方参数，完全不知道你传入的文件名参数是什么，因此会判定必填的x参数缺失，抛出你看到的错误。

内容的提问来源于stack exchange，提问作者James Ball