Shapely缓冲操作异常：几何对象出现意外孔洞问题咨询

我在使用Shapely处理闭合LineString时，执行buffer(80)操作后生成的几何对象出现意外孔洞。该LineString已通过is_valid和is_simple验证，移除一对坐标后孔洞消失；使用Sedona的ST_Buffer时该问题出现频次更高。

相关代码如下：

from shapely import LineString
from geopandas import GeoDataFrame

l = LineString([
  (250,447),
  (319,446),
  (325,387),
  (290,374),
  (259,378),
  (254,385),
  (240,409),
  (244,440),
  (250,447),
])

assert l.is_valid
assert l.is_simple

GeoDataFrame({'geometry': [
  l,
  l.buffer(80),
]}).plot(column='geometry')

生成的带孔洞的缓冲几何效果如下：

问题原因 #

这种情况本质是闭合线的局部线段间距过小，缓冲操作后线段的缓冲区域发生自相交，被几何引擎误判为内部孔洞。虽然原LineString通过了is_valid和is_simple验证（说明线本身无自相交），但当缓冲距离足够大时，原本独立的线段缓冲范围会重叠交叉，触发拓扑计算的异常。Sedona基于JTS引擎，和Shapely依赖的GEOS同源，这类拓扑敏感操作在复杂线结构下更容易出现此类问题。

解决办法 #

• 简化线几何：移除冗余或过于接近的坐标点，可使用shapely.simplify()设置合理容差，在保留整体形状的前提下减少顶点数量：

  simplified_l = l.simplify(tolerance=1.0)
  buffered_simplified = simplified_l.buffer(80)
  

• 调整缓冲参数：修改buffer方法的join_style参数，改用斜接或斜切风格替代默认圆角，避免局部缓冲区域重叠：

  from shapely.geometry.polygon import JOIN_BEVEL
  buffered_bevel = l.buffer(80, join_style=JOIN_BEVEL)
  

• 修复缓冲后几何：对带孔洞的多边形执行buffer(0)操作，利用几何引擎的拓扑修复能力自动消除孔洞：

  fixed_buffer = l.buffer(80).buffer(0)
  

内容的提问来源于stack exchange，提问作者A. West