如何列出同模式栅格对象组？多城市栅格文件批量处理问询

嘿，我懂你现在的痛点——手动给150个城市一个个加载不同模式的栅格文件太麻烦了对吧？之前我处理过类似的批量栅格任务，给你两个实用的解决方案，既能节省时间又不容易出错：

方法1：基础循环（新手友好，逻辑清晰） #

这种方法用最基础的R循环实现，容易理解和调试，适合刚接触批量处理的朋友：

library(raster)

# 设置你的文件夹路径
folder <- "D:/Folder_example"

# 筛选出所有符合你指定模式的tif文件，获取完整路径
all_tifs <- list.files(path = folder, 
                       pattern = "(_E|_ED|_SC_t|_SH_t|_WH_t)\\.tif$", 
                       full.names = TRUE)

# 提取所有唯一的城市名称前缀（去掉文件名里的模式后缀）
city_names <- unique(sub("(_E|_ED|_SC_t|_SH_t|_WH_t)\\.tif$", "", basename(all_tifs)))

# 创建一个空列表，用来存储每个城市的所有栅格对象
city_rasters <- list()

# 循环处理每个城市
for(city in city_names) {
  # 筛选出当前城市的所有模式文件
  city_files <- grep(city, all_tifs, value = TRUE)
  
  # 分别加载每种模式的栅格
  E <- raster(grep("_E.tif", city_files, value = TRUE))
  ED <- raster(grep("_ED.tif", city_files, value = TRUE))
  SC_t <- raster(grep("_SC_t.tif", city_files, value = TRUE))
  SH_t <- raster(grep("_SH_t.tif", city_files, value = TRUE))
  WH_t <- raster(grep("_WH_t.tif", city_files, value = TRUE))
  
  # 把当前城市的所有栅格存入列表，用城市名作为键
  city_rasters[[city]] <- list(E = E, ED = ED, SC_t = SC_t, SH_t = SH_t, WH_t = WH_t)
  
  # 可选：打印进度，方便跟踪处理情况
  cat("完成城市", city, "的栅格加载\n")
}

处理完成后，你可以通过city_rasters[["city1"]]$E这样的方式，快速调用某个城市的特定模式栅格。

方法2：用purrr实现简洁的批量处理（适合熟悉tidyverse的用户） #

如果你习惯用tidyverse风格的代码，purrr包的map函数能让代码更简洁优雅：

library(raster)
library(purrr)

folder <- "D:/Folder_example"
all_tifs <- list.files(path = folder, 
                       pattern = "(_E|_ED|_SC_t|_SH_t|_WH_t)\\.tif$", 
                       full.names = TRUE)
city_names <- unique(sub("(_E|_ED|_SC_t|_SH_t|_WH_t)\\.tif$", "", basename(all_tifs)))

# 定义一个加载单个城市栅格的函数
load_city_rasters <- function(city_name) {
  city_files <- grep(city_name, all_tifs, value = TRUE)
  list(
    E = raster(grep("_E.tif", city_files, value = TRUE)),
    ED = raster(grep("_ED.tif", city_files, value = TRUE)),
    SC_t = raster(grep("_SC_t.tif", city_files, value = TRUE)),
    SH_t = raster(grep("_SH_t.tif", city_files, value = TRUE)),
    WH_t = raster(grep("_WH_t.tif", city_files, value = TRUE))
  )
}

# 批量处理所有城市
city_rasters <- map(city_names, load_city_rasters)
# 给列表元素命名为城市名，方便调用
names(city_rasters) <- city_names

几个实用的注意事项 #

• 缺失文件处理：如果某个城市缺少某类模式的文件，代码会报错。你可以给每个栅格加载步骤加个判断，比如if(length(grep("_E.tif", city_files)) > 0) {E <- raster(...)} else {E <- NA}，避免中断整个批量任务；
• 大文件优化：如果你的栅格文件很大，建议考虑用terra包（raster的替代包，性能更优），或者把每个城市的栅格合并为一个stack，节省内存；
• 进度跟踪：如果处理150个城市耗时较长，方法1里的cat语句能帮你实时看到处理进度，避免焦虑~

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