如何批量导出R全局环境中的所有sf数据框为Shapefile？

步骤1：加载必要包 #

确保已安装并加载sf包：

library(sf)

步骤2：筛选全局环境中的所有sf对象 #

用以下代码提取全局环境内所有类型为sf的对象：

# 获取全局环境中所有对象名称
all_objects <- ls(envir = .GlobalEnv)

# 筛选出sf类型的对象
sf_objects <- all_objects[sapply(all_objects, function(x) inherits(get(x), "sf"))]

步骤3：批量导出为Shapefile #

指定导出文件夹（例如此处设为当前目录下的shapefiles文件夹），遍历所有sf对象完成批量导出：

# 指定导出目录，不存在则自动创建
output_dir <- "shapefiles"
if (!dir.exists(output_dir)) dir.create(output_dir)

# 循环导出每个sf对象
for (obj_name in sf_objects) {
  # 获取当前sf对象
  sf_obj <- get(obj_name)
  # 拼接导出路径：文件夹+对象名（自动生成.shp等配套文件）
  output_path <- file.path(output_dir, obj_name)
  # 导出Shapefile，已存在则覆盖
  st_write(sf_obj, dsn = output_path, driver = "ESRI Shapefile", delete_layer = TRUE)
  # 可选：打印导出进度
  cat(sprintf("已导出：%s\n", obj_name))
}

关键说明 #

• delete_layer = TRUE：若目标路径已有同名Shapefile，会直接覆盖旧文件，避免因文件存在报错。
• 导出后会自动生成.shp、.shx、.dbf等Shapefile配套文件，无需手动处理。
• 若需自定义导出文件名，可修改output_path中的文件名部分，替换obj_name为自定义字符串即可。

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