形状优美的 LineString 相交是否有问题？

Question

在与 Shapely 合作时，我遇到了一个奇怪的问题。有两个点 p1 和 p2，其中第一个点属于多边形，第二个点不属于多边形。当我尝试查找包含这两个点作为端点的 LineString 与多边形边界线之间的交点时，我收到消息称未找到交点。我想知道，怎么可能？

from shapely.geometry import Polygon as SPolygon, Point, LineString

p1 = Point(5.414213562373095, 2.585786437626905)
p2 = Point(15.17279752753168, -7.172797527531679)

l = LineString([p1, p2])

l1 = LineString([(2, 2), (2, 6)])
l2 = LineString([(2, 6), (6, 6)])
l3 = LineString([(6, 6), (6, 2)])
l4 = LineString([(6, 2), (2, 2)])


sp = SPolygon([(2, 2), (2, 6), (6, 6), (6, 2)])

print "Polygon contains p1:", sp.contains(p1)
print "Polygon contains p2:", sp.contains(p2)

for i, line in enumerate((l1, l2, l3, l4)):
    res = l.intersects(line)
    print "Line {0} intersects l1: {1}".format(i, res) 

这是输出：

Polygon contains p1: True
Polygon contains p2: False
Line 0 intersects l1: False
Line 1 intersects l1: False
Line 2 intersects l1: False
Line 3 intersects l1: False

Answer 1

我换到l.intersects(line)了l.intersection(line)，并且确实在以下位置找到了一个交叉路口：

LINESTRING (6 6, 6 2) at Point (6, 2) 

不知道为什么.instersects()表现不同。

然后我舍入了 p1 和 p2

p1 = Point(round(5.414213562373095, 2), round(2.585786437626905, 2))
p2 = Point(round(15.17279752753168, 2), round(-7.172797527531679, 2))

我有 2 个十字路口

LINESTRING (6 6, 6 2) at POINT (6 2)
LINESTRING (6 2, 2 2) at POINT (6 2)

此修复程序也适用于.intersects()(2 True's)

Shapely 对浮点精度可能有点敏感，我通常可以通过舍入来解决问题。尽管这对你来说可能无法接受。