Shapely.intersection方法返回错误结果问题咨询

问题原因 #

你遇到的是Shapely v2.0.1的浮点数精度问题——虽然线段b完全在a上，但浮点数计算的微小误差让Shapely误判了两者的重合关系，导致返回中点而非整条线段b。

验证共线关系：

• 线段a的参数方程是 x=30+6t，y=30t（t∈[0,1]）
• 代入b的起点(32.8,14)：t=(32.8-30)/6≈0.4667，y=30*0.4667≈14，完全匹配
• 代入b的终点(35.2,26)：t=(35.2-30)/6≈0.8667，y=30*0.8667≈26，也完全匹配

可行解决方案 #

1. 先判断包含关系再处理 #

用relate方法获取DE-9IM矩阵，确认b是否完全在a内，再手动返回正确结果：

from shapely.geometry import LineString

a = LineString([[30.0,0.0],[36.0,30.0]])
b = LineString([[32.8,14.0],[35.2,26.0]])

rel_matrix = a.relate(b)
# 完全包含的DE-9IM矩阵为"1F2FFF102"或"F12FFF102"
if rel_matrix in ("1F2FFF102", "F12FFF102"):
    intersection = b
else:
    intersection = a.intersection(b)

print(intersection)  # 输出LineString((32.8, 14), (35.2, 26))

2. 用缓冲规避精度问题 #

给线段a加个极小缓冲，再计算交集后合并结果：

from shapely.geometry import LineString
from shapely.ops import linemerge

a = LineString([[30.0,0.0],[36.0,30.0]])
b = LineString([[32.8,14.0],[35.2,26.0]])

# 添加1e-9的缓冲，不影响实际结果但能解决精度误判
buffered_a = a.buffer(1e-9)
temp_intersection = buffered_a.intersection(b)
# 把结果合并成LineString
intersection = linemerge(temp_intersection)

print(intersection)

3. 升级Shapely版本 #

Shapely v2.1及以上版本修复了不少共线线段的精度问题，直接升级大概率能解决：

pip install --upgrade shapely>=2.1.0

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