用scipy测量地理距离

Question

我没有使用scipy的pdist函数的结果。我感兴趣的是真实的地理距离(首选单位:公里)。取以下坐标:

from scipy.spatial.distance import pdist
coordinates = [ (42.057, -71.08), (39.132, -84.5155) ]
distance = pdist(coordinates)
print distance
# [ 13.75021037]

但单位是什么?谷歌表示，这两个点之间的距离为1179公里。从13.75021037到那里怎么走?

Answer 1

来自scipy的pdist方法不支持lon, lat坐标的距离，如注释所述。

但是，如果您想获得pdist返回的那种距离矩阵，您可以使用pdist方法和geopy包中提供的距离方法。为此，pdist允许使用带有两个参数的自定义函数(lambda函数)计算距离。

下面是一个例子:

from scipy.spatial.distance import pdist
from geopy.distance import vincenty
import numpy as np
coordinates = np.array([[19.41133431, -99.17822823],
                        [19.434514  , -99.180934],
                        [19.380412  , -99.178789])
# Using the vincenty distance function.
m_dist = pdist(coordinates, # Coordinates matrix or tuples list
               # Vicenty distance in lambda function
               lambda u, v: vincenty(u, v).kilometers)

Answer 2

使用最新的Python 3，现在会给出弃用警告。我发现@cffk的答案更容易理解:

(为了方便粘贴在这里)

>>> from geopy.distance import geodesic, great_circle
>>> p1 = (31.8300167,35.0662833) # (lat, lon) - https://goo.gl/maps/TQwDd
>>> p2 = (31.8300000,35.0708167) # (lat, lon) - https://goo.gl/maps/lHrrg
>>> geodesic(p1, p2).meters
429.1676644986777
>>> great_circle(p1, p2).meters
428.28877358686776