R语言中如何循环网格单元经纬度并匹配读取对应文件？

原代码无法得到结果的核心问题：

• 正则表达式错误：list.files的pattern写法不正确，无法匹配所有目标文件
• 条件判断逻辑混乱：直接用数值和文件名字符串做比较，且逻辑运算符使用错误
• 未处理部分文件无.txt后缀的情况

步骤1：正确获取所有目标文件 #

先修正文件列表的匹配规则，覆盖带/不带.txt后缀的文件：

setwd("/Users/salah95/FD")
# 匹配所有以flux_开头的文件，兼容有无.txt后缀的情况
flux_data <- list.files(pattern = "^flux_.*")

步骤2：提取文件名经纬度并匹配读取 #

放弃低效的嵌套循环，先从文件名中解析出经纬度数值，再和gc中的数据精准匹配，同时完成数据读取与计算：

# 读取网格单元数据
gc <- read.table(file='/Users/salah95/streamflow/bcccsm11/grid_cell_ids.txt', header=TRUE, sep = "")

# 初始化列表存储结果
result_list <- list()

# 遍历每个flux文件
for (file in flux_data) {
  # 从文件名中提取经纬度：去除前缀flux_和后缀.txt，拆分出数值
  clean_file <- sub("^flux_|\\.txt$", "", file)
  file_parts <- strsplit(clean_file, "_")[[1]]
  file_lat <- as.numeric(file_parts[1])
  file_long <- as.numeric(file_parts[2])
  
  # 找到gc中经纬度完全匹配的行
  match_idx <- which(gc$Lat == file_lat & gc$Long == file_long)
  
  if (length(match_idx) > 0) {
    # 读取文件数据
    raw_data <- read.delim2(file, header=TRUE, sep = "")
    
    # 计算气温和降水变化（示例：均值/总和，可根据需求修改）
    temp_mean <- mean(raw_data$temperature, na.rm = TRUE)
    precip_total <- sum(raw_data$precipitation, na.rm = TRUE)
    
    # 关联流域信息并存储结果
    for (idx in match_idx) {
      result_key <- paste(gc$Drainage[idx], file_lat, file_long, sep = "_")
      result_list[[result_key]] <- list(
        流域名称 = gc$Drainage[idx],
        纬度 = file_lat,
        经度 = file_long,
        气温均值 = temp_mean,
        降水总量 = precip_total,
        原始数据 = raw_data # 不需要可删除
      )
    }
  }
}

# 将结果转换为数据框，方便后续分析
result_df <- do.call(rbind, lapply(result_list, as.data.frame))

关键修正说明 #

• 正则与文件名处理：用^flux_.*匹配所有目标文件，通过sub清理文件名冗余部分，确保能正确提取经纬度数值
• 匹配逻辑优化：用which定位匹配行，避免嵌套循环的冗余计算
• 结果结构化存储：用列表存储多维度结果，最后转成数据框，便于后续统计或可视化

额外优化建议 #

• 若存在浮点数精度误差，可改用dplyr::near(gc$Lat, file_lat, tol = 1e-6)替代==做匹配
• 文件数量较多时，可使用purrr包的函数简化循环逻辑
• 可先预先生成gc的经纬度索引，进一步提升匹配效率

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