Shapely交集函数未达预期：多边形与线交集返回结果异常

核心原因拆解 #

Shapely的几何运算基于GEOS库，结果完全由拓扑关系和浮点数精度决定，你遇到的右侧线返回Point的情况，基本逃不出这几个场景：

1. 右侧红线仅与交集多边形单点接触 #

如果右侧红线刚好只落在交集多边形的顶点上，而非与多边形的某条边重叠，Shapely会严格返回Point类型——因为拓扑上这属于"触碰"（touch）而非"相交重叠"。
可以用right_line.touches(intersection_polygon)验证，返回True就实锤了这个情况。

2. 浮点数精度引发的拓扑异常 #

GEOS的默认精度阈值约为1e-8，当两条几何的坐标差小于这个值时，会被判定为重合或单点相交。比如本该和多边形边重叠的红线，因为坐标计算时的微小误差，被识别成仅交于一点。
解决办法是用buffer(0)修复多边形拓扑：

# 修复交集多边形的拓扑问题（自动处理自相交、顶点重合等）
cleaned_polygon = intersection_polygon.buffer(0)
# 重新计算交集
right_result = cleaned_polygon.intersection(right_line)

3. 右侧红线的几何定义错误 #

检查右侧红线的坐标是否符合预期：比如红线是不是只画到了多边形的顶点，没有覆盖整条边；或者端点坐标写错了。直接打印红线坐标对比多边形边界：

print("右侧红线坐标：", list(right_line.coords))
print("交集多边形边界：", list(intersection_polygon.boundary.coords))

快速验证代码 #

用这段代码定位具体原因：

# 检查拓扑关系
print("右侧线与多边形触碰：", right_line.touches(intersection_polygon))
print("右侧线与多边形相交：", right_line.intersects(intersection_polygon))

# 输出关键坐标
print("\n交集多边形顶点：")
for coord in intersection_polygon.exterior.coords:
    print(f"  {coord}")

如果touches为True，说明是单点接触；如果intersects为True但结果是Point，那就是精度问题，用buffer(0)修复即可。

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