如何减少raster包中11000个栅格镶嵌操作的内存占用？

解决大数量栅格合并内存溢出的方案 #

我之前处理过类似的上万栅格合并的需求，那种内存直接拉满的痛苦太懂了！给你几个亲测有效的解决方案，应该能帮你搞定：

1. 换用terra包（最推荐） #

raster包的内存管理确实在处理超大量栅格时力不从心，它的升级版terra在内存优化上做了极大改进，尤其是支持虚拟栅格集合，不需要把所有栅格加载到内存就能完成合并。

示例代码：

library(terra)
# 获取所有栅格文件路径（根据你的格式调整pattern）
raster_paths <- list.files("你的栅格存储目录", pattern = "\\.tif$", full.names = TRUE)
# 创建虚拟栅格集合（不会一次性加载所有栅格到内存）
raster_stack <- sprc(raster_paths)
# 执行合并，直接写入磁盘，全程内存占用极低
merged_raster <- merge(raster_stack, filename = "最终合并结果.tif", overwrite = TRUE, fun = mean)

terra::merge()会自动分块处理栅格数据，只在需要时读取部分数据到内存，完全避免一次性加载11000个栅格的噩梦。

2. 优化raster包的分批合并逻辑 #

如果你坚持要用raster包，核心是不要一次性加载所有栅格到内存，并且每次合并后都将结果写入磁盘，释放内存。你之前分批出问题可能是没做好“落地磁盘+内存回收”的步骤：

示例代码：

library(raster)
# 配置raster的内存参数（配合磁盘写入生效）
rasterOptions(chunksize = 1e6, maxmemory = 1e8, overwrite = TRUE)

# 获取栅格路径
raster_paths <- list.files("你的栅格存储目录", pattern = "\\.tif$", full.names = TRUE)
# 设定每批处理的数量（根据你的机器内存调整，比如100个一批）
batch_size <- 100
total_batches <- ceiling(length(raster_paths)/batch_size)

# 处理第一批并写入磁盘
first_batch <- lapply(raster_paths[1:batch_size], raster)
current_merge <- do.call(mosaic, c(first_batch, fun = mean, filename = "temp_batch_1.tif"))
# 强制释放内存
rm(first_batch); gc()

# 循环处理剩余批次
for (i in 2:total_batches) {
  start <- (i-1)*batch_size + 1
  end <- min(i*batch_size, length(raster_paths))
  # 加载当前批次栅格
  current_batch <- lapply(raster_paths[start:end], raster)
  # 合并当前批次并写入临时文件
  batch_merge <- do.call(mosaic, c(current_batch, fun = mean, filename = paste0("temp_batch_", i, ".tif")))
  # 合并临时批次结果与总结果，写入新的临时文件
  current_merge <- mosaic(current_merge, batch_merge, fun = mean, filename = "temp_total_merge.tif")
  
  # 清理临时文件和变量，释放内存
  file.remove(paste0("temp_batch_", i, ".tif"))
  rm(current_batch, batch_merge); gc()
  # 替换总结果文件
  file.remove("temp_batch_1.tif")
  file.rename("temp_total_merge.tif", "temp_batch_1.tif")
}

# 重命名最终结果
file.rename("temp_batch_1.tif", "最终合并结果.tif")

关键注意点：

• 每次调用mosaic都必须指定filename，让结果直接写入磁盘，而不是留在内存
• 每次处理完一批就用rm()删除变量，gc()强制回收内存
• 临时文件要及时清理，避免磁盘空间浪费

3. 检查栅格的一致性 #

确保所有待合并的栅格投影、分辨率、范围完全一致，否则mosaic会自动做重投影/重采样，这会额外占用大量内存。如果不一致，先批量统一栅格的空间参数再合并。

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