遍历文件夹CSV执行空间相交操作时遇列未定义错误求助

核心错误原因 #

• 管道逻辑混乱+变量混淆：原代码在st_transform后直接赋值裁剪范围e，打断了管道链；后续st_intersection(x, e)中的x是文件名（字符串），并非处理后的空间对象，导致sf包试图从字符串中读取LON/LAT列，触发"undefined columns selected"错误。
• 文件读取适配问题：使用read.table读取CSV时，默认分隔符为空白符，无法正确解析逗号分隔的CSV结构，可能导致LON/LAT列无法被识别。
• 代码结构错误：匿名函数中管道与赋值操作混用，破坏了数据处理的连续逻辑。

修正后的代码 #

library(sf)
library(dplyr)

setwd("C:\\mydirectory\\")

# 精准匹配以.csv结尾的文件，避免误匹配
filelist <- list.files(pattern = "\\.csv$")

lapply(filelist, function(filename) {
  # 读取CSV文件，read.csv更适配逗号分隔格式
  df <- read.csv(filename, header = TRUE)
  
  # 提前定义裁剪范围，独立于管道链
  clip_extent <- st_as_sfc(st_bbox(c(xmin = 331380, xmax = 338080, 
                                     ymax = 4698700, ymin = 4691000), 
                                   crs = st_crs(32619)))
  
  # 连续的空间处理流程
  df %>%
    st_as_sf(coords = c("LON", "LAT"), crs = 4326, remove = FALSE) %>%
    st_transform(crs = 32619) %>%
    st_intersection(clip_extent) %>%
    st_drop_geometry() %>% # 导出CSV时移除空间几何列（若需保留可删除此行）
    write.csv(paste0("new_", filename), row.names = FALSE)
})

关键修改说明 #

1. 文件读取优化：替换read.table为read.csv，确保正确解析逗号分隔的CSV文件，避免列识别错误。
2. 变量命名与逻辑分离：将函数参数名改为filename，读取后的数据存为df，避免与空间对象混淆；裁剪范围提前定义，不干扰管道流程。
3. 管道逻辑修复：保持管道链的连续性，st_intersection使用管道传递的空间对象与裁剪范围相交，不再误用文件名。
4. 几何列处理：添加st_drop_geometry()，导出CSV时自动移除空间几何列（若需要保留几何信息，可删除该行）。
5. 文件匹配精准化：修改list.files的匹配规则为\\.csv$，仅匹配以.csv结尾的文件，避免误匹配其他含.csv的文件名。

内容的提问来源于stack exchange，提问作者Tiffany Rae Duhl