R语言批量处理栅格时CRS参数报错问题求助

我来帮你梳理下问题所在，你遇到的Error in CRS(x) : PROJ4 argument-value pairs must begin with +错误，核心原因是你直接把文件路径字符串传给了projectRaster和resample函数，但这两个函数需要的是已经读入内存的RasterLayer对象，不是文件路径。我们一步步修正你的代码：

代码里的关键问题： #

• 没有先读取栅格文件就直接进行投影/重采样操作，导致函数误把路径字符串当成CRS参数解析，触发错误
• 投影和重采样的结果没有赋值给变量，相当于白操作了，后续无法调用结果保存
• sub函数的参数顺序搞反了，而且最终要保存的变量r根本没定义
• writeRaster的语法有多余的逗号，会导致语法错误

修正后的完整代码： #

# 设置工作目录（可选，方便后续文件管理）
setwd("E:\\Proyecto PM2.5\\2_PM_2.5_Processing\\Test")

# 获取目录下所有tif文件的完整路径（用\\.精准匹配后缀，避免误判带tif字符的其他文件）
dir <- "E:\\Proyecto PM2.5\\2_PM_2.5_Processing\\Test"
files <- list.files(path = dir, pattern = '\\.tif$', full.names = TRUE)

# 读取参考栅格（提前确认这个栅格的CRS是EPSG:32614，否则重采样结果会不符合预期）
r_ref <- raster("E:\\Proyecto PM2.5\\3_PM_2.5_Entrega\\temporal\\Raster_C.tif")

# 批量处理每个栅格
for (file in files) {
  # 1. 读取当前栅格文件到内存
  r <- raster(file)
  
  # 2. 重投影到EPSG:32614（新版本R推荐用sf::st_crs(32614)$wkt替代+init写法，需先安装sf包）
  r_proj <- projectRaster(r, crs = "+init=epsg:32614")
  
  # 3. 根据参考栅格执行重采样（最近邻法）
  r_resampled <- resample(r_proj, r_ref, method = "ngb")
  
  # 4. 生成保存文件名（这里假设你想把原文件名里的"CDTDB"替换成"ZMVM"，按需调整参数）
  savename <- sub("CDTDB", "ZMVM", basename(file))
  
  # 5. 保存处理后的栅格，指定格式并允许覆盖已存在文件
  writeRaster(r_resampled, file = savename, format = "GTiff", overwrite = TRUE)
}

额外实用提示： #

• 处理500多幅栅格时，用pbapply包的pblapply替代for循环，可以看到实时处理进度，更直观
• 若你的R版本较新，+init=epsg:xxxx的写法会弹出警告，建议改用crs = sf::st_crs(32614)$wkt（需先运行install.packages("sf")安装包）
• 处理前可以先拿1-2个测试文件跑一遍代码，确认结果符合预期后再批量处理，避免浪费时间

内容的提问来源于stack exchange，提问作者Judá García