如何在R中堆叠各子文件夹中的独立栅格图层？

嗨，别担心，作为R新手遇到这种批量处理的问题太正常了，我来帮你调整代码，轻松拿到每个子文件夹对应的栅格堆叠！

先说说你原来代码里的小问题 #

• 路径没加引号（path/parentfolder/需要用双引号括起来才合法）
• 循环里用setwd()切换工作目录容易出混乱，其实完全没必要
• 没有保留每个子文件夹的堆叠结果，最后只会剩下最后一次循环的输出

修改后的完整代码 #

# 1. 获取父文件夹下的所有子文件夹（注意路径要加引号！）
list_dirs <- list.dirs(path = "path/parentfolder/", recursive = FALSE)

# 2. 初始化一个空列表，用来存储每个子文件夹的栅格堆叠
stack_list <- list()

# 3. 遍历每个子文件夹
for (i in seq_along(list_dirs)) {
  # 获取当前子文件夹的完整路径
  current_dir <- list_dirs[i]
  
  # 列出当前文件夹下符合"cool"模式的栅格文件，full.names=TRUE获取完整路径
  raster_files <- list.files(path = current_dir, pattern = "cool", full.names = TRUE, recursive = FALSE)
  
  # 检查是否找到文件，避免空文件导致报错
  if (length(raster_files) > 0) {
    # 生成当前文件夹的栅格堆叠
    current_stack <- stack(raster_files)
    
    # 将堆叠存入列表，同时用子文件夹名称给列表元素命名，方便后续识别
    stack_list[[basename(current_dir)]] <- current_stack
  } else {
    # 如果没找到文件，给出警告提示
    warning(paste("⚠️ 没有找到符合模式的文件：", current_dir))
  }
}

代码细节解释 #

• 不用setwd()：通过full.names=TRUE直接获取栅格文件的完整路径，彻底避免切换工作目录带来的路径混乱
• 用列表存结果：stack_list会按子文件夹名称保存10个栅格堆叠，你可以通过stack_list[[1]]或者stack_list[["子文件夹名称"]]快速访问对应堆叠
• 容错处理：加了文件数量判断，防止某个子文件夹没有符合条件的文件时脚本直接崩溃

后续使用小技巧 #

运行完代码后，stack_list就是你要的10个栅格堆叠的集合啦！比如想查看第一个子文件夹的堆叠信息，直接输入：

stack_list[[1]]

或者用子文件夹名称精准查找：

stack_list[["subfolder_name"]]

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