已知经纬度坐标求两点间距离，用python表示

已知经纬度坐标求两点间距离，用python表示

已知地球上任意两点（lon1, lat1）,（ lon2, lat2）的经纬度坐标，求两点间的距离用haversine表示：
1.首先先将经纬度坐标的角度化成弧度（mlon1,mlat1）,(mlon2,mlat2)
2.利用公式：
其中a表示两点维度的差值，即 a = mlat1 - mlat2
其中a表示两点经度的差值，即 b = mlon1 - mlon2
其中r表示地球的半径

---

用python表示为

def geo_distance(lon1, lat1, lon2, lat2):  #地理距离
    """
    Calculate the great circle distance between two points
    on the earth (specified in decimal degrees)
    """
    lon1, lat1, lon2, lat2 = map(radians, map(float, [lon1, lat1, lon2, lat2]))  #根据提供的函数对指定序列做映，radians:将角度转换为弧度。
    dlon = lon2 - lon1
    dlat = lat2 - lat1
    a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
    c = 2 * asin(sqrt(a))
    r = 6371
    return c * r
print(geo_distance(111.5,36.08,121.47,31.23))  #我老家到上海的距离

结果

1067.9651712083848 #单位 km

其实python自带haversine库，可以直接调用，这样就可以不记公式啦！

from haversine import haversine
linfen = (111.5,36.08)
shanghai = (121.47,31.23)
dis = haversine(linfen,shanghai)
print(dis)

结果

1133.5471931316163 #单位 km

可以看出还是有差距的，那那个准确呢？ ^_

用那种方法知道了吧！嘿嘿