创建R循环以从目录中读取shapefile并对每个shapefile执行分区统计

# 安装并加载sf包
install.packages("sf")
library(sf)

# 设置目录路径
shapefile_directory <- "/path/to/shapefile_directory"

# 获取目录中所有的shapefile文件
shapefile_files <- list.files(shapefile_directory, pattern = "\\.shp$", full.names = TRUE)

# 循环处理每个shapefile文件
for (shapefile in shapefile_files) {
  # 读取shapefile数据
  shapefile_data <- st_read(shapefile)
  
  # 执行分区统计
  # 这里仅作为示例，你可以根据具体需求进行相关统计分析
  partition_statistics <- table(shapefile_data$partition_column)
  
  # 输出分区统计结果
  print(partition_statistics)
}