terra在Rscript模式下计算数千栅格图层sum时文件不存在报错求助

问题现象 #

• 使用terra::rast加载2000个小体积GeoTIFF（单个<100kb），R交互模式下可正常加载并通过sum计算总值
• Rscript脚本模式下，rast能正确加载图层、打印对象，但调用sum时抛出错误：Error: [sum] file does not exist: ...
• 已验证所有文件路径真实存在，且限制加载1000个图层时sum可正常运行
• 环境信息：terra版本1.8.60，Ubuntu 22.04.5 LTS服务器，内存1TB，临时目录所在SSD剩余空间>5TB

原因推测 #

核心差异来自terra在交互/非交互模式下的临时文件管理机制：

• 交互模式中，R会话的临时目录上下文更稳定，临时文件会被保留至会话结束
• 非交互模式（Rscript）下，当栅格数量超过阈值（如1000）时，terra内部生成的临时文件可能被提前清理，导致sum调用时无法找到依赖文件
• 排除临时空间不足、文件路径错误的可能，更倾向于临时文件生命周期管理的模式差异

解决方案 #

已验证可行方案：用terra::vrt替代rast加载大量栅格 #

VRT（虚拟栅格）是GDAL的虚拟格式，不会实际合并文件，仅创建源文件的引用数据集，从根源避免了临时文件管理问题：

library(terra)

# 根据运行模式选择加载逻辑
if (interactive()) {
    # 交互模式下用rast直接加载，操作更直观
    layers <- rast(file.path(preproc_folder, scen_list$file))
} else {
    # 非交互模式下用vrt，-separate参数确保每个源文件作为独立波段加载
    layers <- vrt(file.path(preproc_folder, scen_list$file), options = c("-separate"))
}

# 计算所有图层的总和
sum_layers <- sum(layers)

其他可选方案 #

若需继续使用rast，可尝试以下调整：

1. 指定稳定的自定义临时目录

# 设置不会被自动清理的临时目录
terraOptions(tempdir = "/path/to/stable/temp/folder")
layers <- rast(files)
sum_layers <- sum(layers)

2. 分批次计算后合并

# 按500个为一批拆分文件列表
batch_size <- 500
batches <- split(files, ceiling(seq_along(files)/batch_size))

# 计算每个批次的局部和
batch_sums <- lapply(batches, function(batch) {
    sum(rast(batch))
})

# 合并所有批次的结果得到总和
sum_layers <- sum(rast(batch_sums))

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