如何在R中填充带孔洞的time trade area shapefile

环境准备 #

首先加载处理矢量所需的依赖包：

# 加载空间处理核心包
library(sf)
library(lwgeom)
library(dplyr)
# 若包未安装，可先执行：install.packages(c("sf", "lwgeom", "dplyr"))

读取与拓扑校验 #

首先读取生成的初始可达区域文件，先修复基础拓扑错误（多数不规则穿孔是拓扑错误导致的）：

# 读取初始shapefile，替换为本地实际存储路径
trade_area <- st_read("初始time_trade_area.shp")
# 自动修复无效拓扑、重叠边界等基础问题
trade_area_valid <- st_make_valid(trade_area)

一键填充内部孔洞 #

lwgeom包提供了专门的多边形孔洞填充接口，可直接批量处理所有要素的内环孔洞：

# 默认填充所有多边形内部的孔洞
trade_area_filled <- st_fill_holes(trade_area_valid)

注：如果需要保留指定面积以上的孔洞（比如大型湖泊等不需要填充的区域），可传入threshold参数指定孔洞面积阈值，仅填充小于阈值的孔洞。示例：st_fill_holes(trade_area_valid, threshold = 10000) 仅填充面积小于10000平方米的孔洞。

修复碎片化不完整问题 #

如果填充后存在离散的小碎片、同属一个时间范围的要素被拆分的情况，可按属性做融合处理：

# 按时间范围标识字段融合要素，需替换为你的shapefile中实际的属性字段名
trade_area_final <- trade_area_filled %>% 
  group_by(time_range) %>% 
  summarise(geometry = st_union(geometry)) %>% 
  ungroup()

兼容提示：如果使用传统sp格式的空间数据，可先通过st_as_sf()转换为sf格式处理，完成后再通过as(处理后的数据, "Spatial")转回sp对象。

导出最终结果 #

# 导出处理完成的无孔洞shapefile，设置delete_layer=TRUE可覆盖同名旧文件
st_write(trade_area_final, "time_trade_area_fixed.shp", delete_layer = TRUE)

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