R语言批量读取tif栅格文件 两种读取方法均报错问题咨询

问题背景 #

在R环境中读取指定目录下的多个TIF栅格文件时，先后采用两种实现方式，均触发完全相同的报错：

• 遍历文件列表逐个读取栅格，执行代码：allrasters1 <- lapply(rastlist, raster)
• 直接批量堆叠栅格文件，执行代码：allrasters2 <- stack(rastlist)

已尝试两类排查操作均未解决问题：

• 更新所有相关R依赖包
• 重新安装最新正式版本R语言

问题复现完整代码如下：

library(sp)
library(raster)
library(rgdal)
rastlist <- list.files(path = "D:\\Learn\\R\\Clip\\rasters", pattern='.tif', 
                       all.files=TRUE, full.names=FALSE)
allrasters1 <- lapply(rastlist, raster)
allrasters2 <- stack(rastlist)

报错诱发原因 #

1. 文件路径配置错误：list.files()函数中full.names参数设为FALSE时，仅返回文件名本身，不会拼接文件所在的完整目录路径。若当前R工作目录与栅格存储目录D:\\Learn\\R\\Clip\\rasters不一致，后续读取操作会在当前工作目录下查找文件，必然因找不到目标文件触发报错。
2. 依赖包兼容性问题：代码中加载的rgdal包已于2023年正式从CRAN退役，新版R环境不再支持该包，即使路径配置正确，也可能因底层栅格驱动调用异常触发读取失败。
3. 后缀匹配逻辑不严谨：pattern='.tif'为普通字符串匹配，会命中文件名中含.tif字段的非栅格文件（如.tif.aux.xml这类栅格辅助文件），导致读取列表混入无效文件触发报错。

可行解决方法 #

方案1：最小修改适配原有raster工作流 #

仅修正参数配置，保留原有代码逻辑即可解决问题，修改后代码如下：

library(sp)
library(raster)
# 移除已退役的rgdal包加载语句，新版raster无需依赖rgdal
rastlist <- list.files(
  path = "D:\\Learn\\R\\Clip\\rasters",
  pattern = '\\.tif$', # 用正则匹配后缀为.tif的文件，过滤辅助文件
  all.files = TRUE,
  full.names = TRUE # 开启完整路径返回，确保R能定位到文件
)
allrasters1 <- lapply(rastlist, raster)
allrasters2 <- stack(rastlist)

方案2：迁移到terra包（长期稳定方案） #

terra是raster包的官方继任者，栅格读写、处理效率更高，无旧版rgdal依赖的兼容问题，适配代码如下：

# 首次运行需安装包
# install.packages("terra")
library(terra)
rastlist <- list.files(
  path = "D:\\Learn\\R\\Clip\\rasters",
  pattern = '\\.tif$',
  all.files = TRUE,
  full.names = TRUE
)
# 逐个读取栅格
allrasters1 <- lapply(rastlist, rast)
# 直接批量堆叠为多波段栅格对象
allrasters2 <- rast(rastlist)

排查提示：若修正配置后仍报错，需检查目标目录下的TIF文件是否损坏、是否存在权限限制导致R无法读取文件。

内容的提问来源于stack exchange，提问作者Mohamed Osman