如何在R中自动化对文件夹内所有栅格数据执行focal模糊处理流程

嗨，我完全懂你现在的困扰——手动逐个读取、处理栅格实在太耗时了，尤其是当你有大量研究站点的栅格要处理时。咱们可以把整个流程自动化，只需要设置好工作目录，剩下的交给代码自动完成！

下面是我整理的完整自动化方案，完全基于你现有的处理逻辑，只是把单个栅格的操作封装成函数，再批量应用到文件夹里的所有栅格：

步骤1：加载包并设置基础参数 #

首先还是加载terra包，然后设置你的工作目录（存放所有待处理栅格的文件夹）和输出目录（如果不想覆盖原文件，建议单独建一个输出文件夹）：

library(terra)

# 设置工作目录（存放待处理栅格的文件夹）
wd <- "C:/Users/nikos/OneDrive/Desktop/psf_paris2/"
# 设置输出目录（可与工作目录相同，也可单独指定）
output_dir <- wd

步骤2：获取文件夹内所有栅格文件 #

用list.files筛选出文件夹里所有的.tif文件，记得开启full.names = TRUE来获取完整文件路径：

# 获取所有tif格式的栅格文件路径
raster_files <- list.files(path = wd, pattern = "\\.tif$", full.names = TRUE)

步骤3：编写单个栅格的处理函数 #

把你原来对单个栅格的处理逻辑封装成一个函数，这样我们就可以反复调用它处理每个栅格：

process_raster <- function(raster_path, output_dir) {
  # 读取栅格
  r <- rast(raster_path)
  # 提取栅格基础名称（比如从"pop.tif"提取"pop"）
  raster_name <- tools::file_path_sans_ext(basename(raster_path))
  
  # 执行你原来的focal循环逻辑
  for (i in seq(from = 0.2, to = 0.8, by = 0.2)) {
    print(paste0("处理中: ", raster_name, " | 参数i = ", i))
    
    # 创建高斯核
    gf <- focalMat(r, i * 400, "Gauss")
    # 执行focal模糊
    r_gf <- focal(r, w = gf, na.rm = TRUE)
    # 聚合栅格
    r_gf <- aggregate(r_gf, fact = 4, fun = "mean", cores = 8)
    
    # 处理i值的命名格式（去掉小数点，比如0.2变成"02"）
    stringedi <- gsub("\\.", "", toString(format(i, nsmall = 2)))
    # 拼接输出文件路径
    output_path <- file.path(output_dir, paste0(raster_name, stringedi, ".tif"))
    
    # 保存处理后的栅格
    writeRaster(r_gf, filename = output_path, overwrite = TRUE)
    
    # 释放内存（处理大量栅格时避免内存溢出）
    rm(r_gf, gf)
    gc()
  }
  
  # 释放当前栅格的内存
  rm(r)
  gc()
}

步骤4：批量处理所有栅格 #

最后，用lapply或者for循环把处理函数应用到所有栅格文件上，两种方式任选：

# 方式1：用lapply批量处理（简洁高效）
lapply(raster_files, process_raster, output_dir = output_dir)

# 方式2：用for循环处理（适合需要逐行查看进度时）
for (file in raster_files) {
  process_raster(file, output_dir = output_dir)
}

额外小提示： #

• 我添加了tools::file_path_sans_ext来自动提取栅格名称，不管你的栅格是pop.tif还是landuse.tif，都能生成对应的输出文件名（比如landuse02.tif）
• 加入了rm()和gc()来及时释放内存，避免处理大量栅格时出现内存不足的问题
• 你只需要修改wd和output_dir这两个参数，剩下的流程完全自动执行

这样一来，你只需要把所有待处理的栅格放到指定文件夹里，运行这段代码就能自动完成所有处理啦！

备注：内容来源于stack exchange，提问作者Nikos