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

Question

有一条折线，其中折线的顶点列表为[[x1，y1），（x2，y2），（x3，y3），…]和一个点（x，y）。在Shapely中，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()计算距离时，该方法应具有坐标。有什么方法可以从此方法返回它吗？

Answer 1

您正在描述的GIS术语是线性参考，Shapely拥有这些方法。

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

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

另一种方法是使用nearest_points：

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

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