如何用Shapely优化大多边形去除边缘小多边形后的闭合形状

当使用Shapely的difference操作处理边缘包含小多边形的大多边形时，结果会出现带凹口的多边形（而非独立闭合形状）。以下是适配简单/复杂多边形的通用解决方案，包含拆分凹口多边形和面积过滤两步：

核心思路 #

1. 拆分凹口多边形：提取小多边形与大多边形的非共享边缘作为切割线，将带凹口的差集结果拆分为多个独立的闭合多边形。
2. 过滤小面积多边形：根据面积阈值自动保留有效多边形，剔除过小的碎片。

完整实现代码 #

from shapely.geometry import Polygon
from shapely.ops import split, unary_union

# 示例多边形（可替换为任意复杂多边形）
large_polygon = Polygon([(0, 0), (0, 4), (6, 4), (6, 0)])
small_polygon = Polygon([(4, 1), (4, 3), (6, 3), (6, 1)])

# 计算初始差集
result = large_polygon.difference(small_polygon)

# 步骤1：提取切割线并拆分多边形
# 获取小多边形边界中不与大多边形重合的部分
small_boundary = small_polygon.boundary
large_boundary = large_polygon.boundary
cut_lines = small_boundary.difference(large_boundary)

# 处理多段切割线的情况
if hasattr(cut_lines, 'geoms'):
    cut_lines = unary_union(cut_lines)

# 拆分凹口多边形为独立闭合形状
split_polygons = split(result, cut_lines)

# 步骤2：过滤面积过小的多边形（阈值可自定义）
area_threshold = 0.5
filtered = [poly for poly in split_polygons.geoms if poly.area > area_threshold]

# 生成最终结果
final_result = unary_union(filtered)

# 输出验证
print(f"拆分后有效多边形数量：{len(filtered)}")
for idx, poly in enumerate(filtered):
    print(f"多边形{idx+1}面积：{round(poly.area, 2)}")

适配复杂多边形的说明 #

• 该方案不依赖直线/垂直等规则判断，完全基于几何边界运算，自动识别非共享边缘作为切割线，适配任意曲线、非规则多边形场景。
• 如果小多边形完全位于大多边形内部（无共享边缘），差集结果会是带孔洞的多边形，此时无需拆分，直接提取外环和内环生成独立多边形后再做面积过滤即可。

注意事项 #

• 面积阈值需根据业务需求调整，避免误删有效形状。
• 若切割线为空（小多边形完全与大多边形重合），差集结果为空，需提前做判断处理。

内容的提问来源于stack exchange，提问作者Vincent Rougeau-Moss