R4.2下raster::writeRaster写入栅格堆栈报错，求排查解决

问题本质 #

你遇到的不是代码逻辑错误，而是raster包旧版本与R4.2+的兼容性bug。R4.2开始对if条件的长度检查更严格，旧版raster的writeRaster在bylayer=TRUE的处理逻辑中，内部判断未正确处理向量长度，从而触发the condition has length > 1错误。

可行解决方法 #

方法1：指定完整文件名路径 #

不要只传图层名给filename参数，构造包含路径和后缀的完整文件名向量，可绕过内部的后缀判断bug：

dsn <- here("Clipped_ENVData/Mask2022/")
# 用file.path构造规范路径，避免分隔符问题
nameT <- file.path(dsn, paste0(names(rastList3), ".tif"))
# 传入完整文件名，保持bylayer=TRUE批量写入
raster::writeRaster(stack(rastList3), filename = nameT, bylayer = TRUE, overwrite = TRUE, format = "GTiff")

方法2：升级raster包 #

开发者已针对R4.2修复该bug，直接升级到最新版本即可解决：

install.packages("raster")

方法3：循环批量写入（兜底方案） #

如果前两种方法无法生效，用循环自动化处理所有图层，效果和手动逐个写入一致：

dsn <- here("Clipped_ENVData/Mask2022/")
nameT <- file.path(dsn, paste0(names(rastList3), ".tif"))

# 用lapply遍历所有图层完成批量写入
lapply(seq_along(rastList3), function(i) {
  raster::writeRaster(rastList3[[i]], filename = nameT[[i]], overwrite = TRUE, format = "GTiff")
})

内容的提问来源于stack exchange，提问作者Lj freer