用Python根据经纬度计算距离

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

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

a = sin ⁡ 2 ( Δ φ / 2 ) + cos ⁡ φ 1 ⋅ cos ⁡ φ 2 ⋅ sin ⁡ 2 ( Δ λ / 2 ) c = 2 ⋅ atan 2 ( a , 1 − a ) d = R ⋅ c \begin{aligned}a &= \sin^2(\Delta\varphi/2)+\cos\varphi_1\cdot\cos\varphi_2\cdot\sin^2(\Delta\lambda/2)\\ c&=2\cdot \text{atan}2(\sqrt{a},\sqrt{1-a})\\ d&=R\cdot c \end{aligned} acd=sin2(Δφ/2)+cosφ1cosφ2sin2(Δλ/2)=2atan2(a ,1a )=Rc

φ \varphi φ是纬度 (latitude), λ \lambda λ是经度 (longitude), R R 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为单位的~

你可能感兴趣的:(Python)