python求球面距离

python求球面距离

今天我们用python代码写一个求球面距离的公式，可以通过经纬度，计算球面间两点的距离。

代码如下：

def geodistance(lng1,lat1,lng2,lat2):
    lng1, lat1, lng2, lat2 = map(radians, [float(lng1), float(lat1), float(lng2), float(lat2)]) # 经纬度转换成弧度
    dlon=lng2-lng1
    dlat=lat2-lat1
    a=sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
    distance=2*asin(sqrt(a))*6371*1000 # 地球平均半径，6371km
    distance=round(distance/1000,3)
    return distance

lng1,lat1，lng2，lat2分别是两对经纬度。通过上述代码可以计算这两对经纬度地点之间的距离。