Shapely Difference运算后结果仍存在重叠的技术咨询

阿华AIGC实验室

2026-4-28

问题解析与解决方案

1. 为什么`difference`操作后仍显示`overlaps=True`？

核心原因是浮点计算精度误差。你的poly1包含一条斜线边界（右侧边从(2008048.8398, 2487460.0675)到(2008046.0958, 2487551.9931)），而Shapely依赖的GEOS库使用浮点数进行几何计算，这类非轴对齐的边界在求差时容易产生微小的精度偏差——差集poly3的边界可能和poly1的边界存在极微小的内部重叠，刚好触发了overlaps()的判断条件（两个几何体的内部有相交区域）。

对比你的测试用例：测试用例中的多边形都是轴对齐的矩形，边界都是水平/垂直的，浮点计算误差可以忽略不计，因此差集和原多边形的边界完美贴合，内部无重叠，overlaps()自然返回False。

2. 为什么交集是`MultiLineString`但`touches=False`？

touches()的判断逻辑是：两个几何体的边界相交，但内部完全不相交。如果由于浮点误差，poly3和poly1存在哪怕极微小的内部重叠，touches()就会返回False——哪怕它们的主要交集是边界线（也就是你看到的MultiLineString）。

而你的测试用例中，差集和原多边形的边界完全贴合，内部没有任何重叠，因此touches()返回True，符合预期。

解决办法

针对这类浮点精度问题，你可以尝试以下几种方案：

修复几何有效性：对差集结果执行buffer(0)操作，它可以自动修复因浮点误差导致的自相交、微小重叠等几何问题：
```
poly3 = poly2.difference(poly1).buffer(0)
```
之后再检查poly3.intersects(poly1)、poly3.overlaps(poly1)和poly3.touches(poly1)，结果应该会符合你的预期。
基于面积判断重叠：如果buffer(0)无法完全解决问题，可以计算交集的面积，当面积小于某个极小阈值（比如1e-6）时，视为精度误差导致的“伪重叠”：
```
poly_int = poly3.intersection(poly1)
if poly_int.area < 1e-6:
    # 视为无内部重叠，按touches处理
    print("实际为边界接触，无内部重叠")
```
坐标整数化（适用场景）：如果你的坐标允许，可以将所有坐标值放大为整数（比如乘以1e6后取整），减少浮点运算的精度损失，计算完成后再缩小回去。