如何在Shapely中生成不自相交的自交叉LineString缓冲区

问题描述 #

我有一条自交叉的LineString，使用LineString().buffer(2, cap_style=BufferCapStyle.flat)生成缓冲区后，结果是带内部洞的自相交多边形，需要生成两个不相交的缓冲区效果。

解决方案 #

核心思路是先拆分自交叉线为无交叉的子线段，再逐个生成缓冲区，最后合并结果。以Shapely库为例，步骤如下：

1. 拆分自交叉线：用unary_union将自交叉的LineString转换为MultiLineString，自动拆分出所有无交叉的子线段。
2. 生成独立缓冲区：遍历每个子线段，调用buffer生成各自的缓冲区。
3. 合并缓冲区：再次用unary_union将所有独立缓冲区合并为MultiPolygon。

代码示例 #

from shapely.geometry import LineString
from shapely.ops import unary_union
from shapely.geometry.polygon import BufferCapStyle

# 定义自交叉的LineString示例
cross_line = LineString([(0, 0), (2, 2), (0, 2), (2, 0)])

# 拆分自交叉线为无交叉子线段
split_lines = unary_union(cross_line)

# 为每个子线段生成缓冲区并合并
non_overlapping_buffer = unary_union([
    line.buffer(2, cap_style=BufferCapStyle.flat) 
    for line in split_lines.geoms
])

原理说明 #

自交叉线直接缓冲时，GIS算法会把交叉重叠区域处理为多边形的内部洞，因此得到带洞的单个多边形。拆分后逐个缓冲，每个子线段的缓冲区相互独立，不会产生交叉重叠，最终合并后就是两个不相交的多边形，符合目标效果。

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