polygon - Pyproj 点之间以及点与多边形之间的距离

Question

我想计算两点之间以及点与多边形之间的投影距离。所有坐标都在相同的投影纬度、经度 (WGS84) 下指定。

我使用 pyproj 计算了一个点和一个多边形之间的距离，如下所示：

from pyproj import Proj, transform, Geod
geod = Geod(ellps='WGS84')
angle1,angle2,dist1 = geod.inv(wLong1, sLat1, wLong2, sLat2)
#this returns distance in m

我想使用相同的函数来计算点和边界框之间的距离。

bbox = box(wLong1, sLat1, eLong1, nLat1)
point = Point(wLong2,sLat2)
dist2 = (point.distance(bbox))

与第一个示例（以米为单位的 dist1）不同，我认为第二个示例（dist2）以度为单位返回距离。如何将此值转换为示例 1 中的仪表？

Answer 1

您需要地球曲率的平均半径 (rm) 进行计算。

from pyproj import Geod
from math import radians

# ... some code

a = Geod.a
b = Geod.b
rm = (2.0*a + b)/3.0   # simple mean radius, as defined by IUGG
rm * radians(dist2)    # your dist in meters

存在更准确的 rm 公式，但以上是很好的近似值。