2

我有一个代表路线的坐标列表（纬度/经度）。给定一定的半径和另一个坐标，我需要检查坐标是否在路线中（在任何点的给定半径内）及其与路线起点的距离。

我查看了 Shapely，它看起来是一个很好的解决方案。

我开始创建一个StringLine

from shapely.geometry import LineString
route = LineString[(x, y), (x1, y1), ...]

然后检查该点是否在路线附近，我添加了一个缓冲区并检查了交叉点

from shapely.geometry import Point
p = Point(x, y)
r = 0.5
intersection = p.buffer(r).intersection(route)
if intersection.is_empty: 
    print "Point not on route"
else: 
    # Calculate P distance from the begning of route

我一直在计算距离。我想过分割路线p并测量前半部分的长度，但我得到的交叉点结果是一个HeterogeneousGeometrySequence我不确定我能做什么。

我相信我找到了解决方案：

if p.buffer(r).intersects(route):
 return route.project(p)
4

1 回答 1

3

与其缓冲一个昂贵且不完美的几何图形（因为缓冲需要许多段和许多其他选项），只需查看该点是否在距离阈值内：

if route.distance(p) <= r:
    return route.project(p)

此外，您现在可能已经意识到您的距离单位是度数。如果您想要线性距离，例如米，则需要使用不同的库使其变得更加复杂。

于 2015-03-08T20:44:00.250 回答