两多边形扣除重叠区域后出现多余点的原因排查

核心原因 #

你遇到的问题本质是浮点数精度误差的累积效应，搭配Shapely处理几何运算的内部逻辑差异导致的：

• 当你先计算overlap = geo_1.intersection(geo_2)时，由于地理坐标以浮点数存储，运算后交集的边界会产生极其细微的精度偏差（比如某个顶点坐标存在十亿分之一级的误差）。
• 接着用geo_1.difference(overlap)切割时，这个微小偏差会被放大，导致Shapely在处理几何边界切割时，误将原多边形上本属于geo_2区域的一个点保留下来，形成所谓的“多余点”。
• 而直接调用geo_1.difference(geo_2)时，Shapely内部会采用更高效的一次性差集算法，直接对比两个原始几何的完整边界，避免了两次独立运算（交集+差集）带来的精度误差叠加，因此结果是正确的。

操作建议 #

• 不需要拆分“先求交集再减交集”的步骤，直接使用单次difference运算即可得到可靠结果，这也是几何差集运算的标准用法。
• 如果必须拆分运算（比如需要复用交集结果），可以对交集结果做一次几何修复，比如调用overlap = geo_1.intersection(geo_2).buffer(0)，buffer(0)可以自动修复因精度问题产生的几何拓扑错误，再执行差集运算就能避免多余点。

验证示例 #

你可以添加以下代码验证差异：

# 查看两个结果的WKT
print("分步运算结果WKT:\n", geo_1_overlap_removed.wkt)
print("\n直接差集结果WKT:\n", geo_1_minus_geo_2.wkt)

# 检查多余点是否属于geo_2
extra_point = list(geo_1_overlap_removed.exterior.coords)[-2]  # 假设多余点是倒数第二个坐标
print("\n多余点是否在geo_2内:", geo_2.contains(shapely.Point(extra_point)))

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