基于Shapely的多边形平滑：膨胀腐蚀拓扑问题解决方案咨询

核心思路 #

你的需求是通过「膨胀+腐蚀」实现保守型平滑，但Shapely的buffer函数默认会自动修复自相交几何，导致拓扑变化时失效。我们可以通过手动提取边界并生成偏移曲线的方式绕过这个限制，保留中间膨胀阶段的自相交状态，再进行腐蚀操作。

Shapely内实现方案 #

使用offset_curve函数生成偏移后的边界线（允许自相交），再构造无效多边形，最后通过负缓冲完成腐蚀并修复拓扑。

代码实现

from shapely.geometry import Polygon, LinearRing
from shapely.ops import offset_curve
from shapely.validation import make_valid

def smooth_conservative(poly, buffer_distance):
    # 提取多边形的外边界和所有内洞环
    rings = [poly.exterior] + list(poly.interiors)
    offset_boundaries = []

    for ring in rings:
        # 生成偏移曲线：外边界逆时针，正距离向外膨胀；内洞顺时针，正距离向内膨胀
        # JOIN_STYLE_ROUND（值为1）确保平滑圆角
        offset_line = offset_curve(ring, buffer_distance, join_style=1)
        # 尝试转为LinearRing，自相交时直接保留LineString（构造多边形时允许无效）
        try:
            offset_ring = LinearRing(offset_line)
        except ValueError:
            offset_ring = offset_line
        offset_boundaries.append(offset_ring)
    
    # 构造膨胀后的多边形（可能自相交、无效）
    dilated = Polygon(offset_boundaries[0], offset_boundaries[1:])
    # 执行腐蚀操作，Shapely会自动修复拓扑生成有效几何
    eroded = dilated.buffer(-buffer_distance, join_style=1)
    # 确保最终结果有效
    smoothed = make_valid(eroded)

    # 验证保守性：确保原始多边形完全被平滑后的图形覆盖
    if not poly.within(smoothed):
        #  fallback到传统buffer方法（拓扑未变化时更可靠）
        smoothed = poly.buffer(buffer_distance).buffer(-buffer_distance)
    
    return smoothed

关键细节

• offset_curve不自动修复自相交，完美保留膨胀后的边界状态；
• 最终的buffer(-distance)会从无效多边形中生成拓扑一致的有效几何；
• join_style=1（圆角连接）确保凹面被平滑，凸起部分因先向外再向内偏移，形状得以保留；
• 加入保守性验证，避免极端情况下原始图形未被覆盖。

替代框架方案（当Shapely无法满足时） #

若需要更精确控制自相交的处理逻辑，推荐使用CGAL的Python绑定，它的多边形偏移函数支持显式允许自相交的参数。

代码示例（CGAL）

from CGAL.CGAL_Polygon_offset import offset_polygon_2
from CGAL.CGAL_Kernel import Point_2
from shapely.geometry import Polygon
from shapely.validation import make_valid

def shapely_to_cgal(poly):
    """将Shapely多边形转换为CGAL格式"""
    exterior = [Point_2(p.x, p.y) for p in poly.exterior.coords]
    interiors = [[Point_2(p.x, p.y) for p in ring.coords] for ring in poly.interiors]
    return (exterior, interiors)

def smooth_with_cgal(poly, buffer_distance):
    cgal_poly = shapely_to_cgal(poly)
    
    # 膨胀：外边界向外偏移，允许自相交
    dilated_exterior = offset_polygon_2(
        cgal_poly[0], 
        buffer_distance, 
        join_style=1,  # 圆角连接
        allow_self_intersections=True
    )[0]
    # 膨胀：内洞向内偏移（负距离对应向内）
    dilated_interiors = []
    for interior in cgal_poly[1]:
        dilated_int = offset_polygon_2(
            interior, 
            -buffer_distance, 
            join_style=1, 
            allow_self_intersections=True
        )[0]
        dilated_interiors.append(dilated_int)
    
    # 转回Shapely多边形（可能无效）
    dilated_shapely = Polygon(
        [(p.x(), p.y()) for p in dilated_exterior],
        [[(p.x(), p.y()) for p in ring] for ring in dilated_interiors]
    )
    # 腐蚀并修复拓扑
    smoothed = make_valid(dilated_shapely.buffer(-buffer_distance, join_style=1))
    return smoothed

优势

CGAL的offset_polygon_2提供了对自相交的精确控制，适合复杂拓扑的多边形，平滑效果更稳定。

内容的提问来源于stack exchange，提问作者Stefan Karlsson