python - 确定 Shapely 点是否在 LineString/MultiLineString 内

Question

我正在尝试使用 Shapely 的within函数对 LineString 和 Point 文件进行“空间连接”（仅供参考，点文件是使用 上的interpolate函数生成的LineString）。问题是 - 没有返回任何内容。

# this condition is never satisfied
if point.within(line):
    # here I write stuff to a file

在哪里：

point = POINT (-9763788.9782693591000000 5488878.3678984242000000)
line = LINESTRING (-9765787.998118492 5488940.974948905, -9748582.801636808 5488402.127570709)

我错过了什么？

Answer 1

在直线上查找点时存在浮点精度错误。改为使用具有适当阈值的距离。

from shapely.geometry import Point, LineString

line = LineString([(-9765787.9981184918, 5488940.9749489054), (-9748582.8016368076, 5488402.1275707092)])
point = Point(-9763788.9782693591, 5488878.3678984242)

line.within(point)  # False
line.distance(point)  # 7.765244949417793e-11
line.distance(point) < 1e-8  # True