在 Shapely 或 GeoPandas 中查找多边形上最近点的坐标

Question

我想找到一个点到一个多边形的距离，但真正的近距离

我正在执行以下操作：

from shapely.geometry import Point, Polygon
import shapely
import shapely.ops
from shapely.ops import nearest_points

coords2 = [(50.848677, 4.338074), (50.833264, 4.344961), (50.840809, 4.366227), (50.852455, 
4.367945 ), (50.858306, 4.346693)]
poly_2 = Polygon(coords2)
p2 = Point(50.811948, 4.382617)
print(p2.within(poly_2))
p2__ = nearest_points(poly_2, p2)
print("****:",p2__[0])

结果是： ****: POINT (50.840809 4.366227) 这是距离几何最近的点。 我想得到真正的最近点

任何想法？

Answer 1

现在这很容易，您知道结果点将位于连接多边形上 2 个点的线段上，并且这些点之一是您正在寻找的最近点。 这应该很容易，因为您只需计算两点之间的中点，如果您的原始结果点位于中点的任一侧，则该点更接近该特定点。 例如：如果该点位于线段 AB 的前半部分，那么最近的点将是 A，但是，如果它位于第二部分，则它更靠近 B 点。

现在剩下的唯一任务是找到穿过计算结果点的线段 AB。 为此，只需遍历多边形的每条边，根据 y=mx+c 获取线的 eq 并计算您的答案是否满足线的 eq。

希望这可以帮助！

Answer 2

nearest_points算法可以按照您的意愿工作，它返回形状上最接近另一个形状的任何位置的点。 只是在你的情况下，多边形的角点恰好是最近的点。 这是一个插图：

您给出的点p2是左上角的蓝色点。 如蓝色圆圈所示，红色多边形的角恰好是离p2最近的点。 如果我们假设一个新的假设点p3的坐标为(50.82, 4.36) （图中橙色），则nearest_points点算法返回多边形上的一个点，而不是角点（在这种情况下为 ~ (50.837, 4.354) ）。