Shapely 中两个几何体最近点的坐标

有一条折线，其顶点坐标列表 = [(x1,y1), (x2,y2), (x3,y3),...] 和一个点 (x,y)。在身材匀称中，geometry1.distance(geometry2)返回两个几何图形之间的最短距离。

>>> from shapely.geometry import LineString, Point
>>> line = LineString([(0, 0), (5, 7), (12, 6)])  # geometry2
>>> list(line.coords)
[(0.0, 0.0), (5.0, 7.0), (12.0, 6.0)]
>>> p = Point(4,8)  # geometry1
>>> list(p.coords)
[(4.0, 8.0)]
>>> p.distance(line)
1.4142135623730951

但我还需要找到最接近点（x，y）的线上点的坐标。在上面的例子中，这是该点的坐标LineString距离 1.4142135623730951 单位的对象Point(4,8)。方法distance()计算距离时应该有坐标。有什么办法让它从这个方法返回吗？

您所描述的 GIS 术语是线性参考 https://en.wikipedia.org/wiki/Linear_referencing, and Shapely有这些方法 https://shapely.readthedocs.io/en/latest/manual.html#linear-referencing-methods.

# Length along line that is closest to the point
print(line.project(p))

# Now combine with interpolated point on line
p2 = line.interpolate(line.project(p))
print(p2)  # POINT (5 7)

另一种方法是使用nearest_points https://shapely.readthedocs.io/en/latest/manual.html#nearest-points:

from shapely.ops import nearest_points
p2 = nearest_points(line, p)[0]
print(p2)  # POINT (5 7)

它提供与线性参考技术相同的答案，但可以从更复杂的几何输入（例如两个多边形）确定最近的点对。