用Python根据经纬度计算距离

做房产数据分析的时候在数据集遇到有经纬度，好奇怎么用Python计算两地的距离。
先查查算距离的公式（地球是立体的是圆哒，算欧几里得距离是要做地鼠的哦）

在这个计算距离网站有，输入两点的经纬度可以给你算出来距离，下面也有给出计算公式：
Haversine公式：计算球面距离

$$\begin{array}{rl}a& ={\sin }^2(\Delta \phi /2)+\cos {\phi }_1\cdot \cos {\phi }_2\cdot {\sin }^2(\Delta \lambda /2)\\ c& =2\cdot \text{atan}2(\sqrt{a},\sqrt{1-a})\\ d& =R\cdot c\end{array}$$

$\phi$是纬度 (latitude)，$\lambda$是经度 (longitude)，$R$在这里就是地球的半径6,371km。atan2是一个三角函数，详情可以看这里，Python的 math.atan2 直接用，不怕不怕。
需要注意：note that angles need to be in radians to pass to trig functions! 经纬度需要转换成弧度，用 math.radians 转换就行。

def haversine(lat1, long1, lat2, long2):
    from math import radians, sin, cos, atan2, sqrt
    lat1, long1, lat2, long2 = map(radians, [lat1, long1, lat2, long2])
    a = sin((lat1-lat2)/2)**2 + cos(lat1)*cos(lat2)*(sin((long1-long2)/2)**2)
    c = 2 * atan2(sqrt(a), sqrt(1-a))
    return 6371 * c

算出来的结果是以km为单位的~