通过两点经纬度求在地球上的距离

Haversine公式

$$\frac{\Pi }{180}\ast deg=rad$$

$$haversin(\theta )=si{n}^2(\frac{\theta }{2})=\frac{1-cos(\theta )}{2}$$

$$hav(\frac{d}{r})=hav({\phi }_1-{\phi }_2)+cos({\phi }_1)cos({\phi }_2)hav({\lambda }_1-{\lambda }_2)$$

$$d=2\ast r\ast arcsin(\sqrt{si{n}^2(\frac{{\phi }_1-{\phi }_2}{2})+cos({\phi }_1)cos({\phi }_2)si{n}^2(\frac{{\lambda }_1-{\lambda }_2}{2})})$$

d is the distance between two points
r is the radius of the earth
${\phi }_1{\phi }_2$ is the latitude at point 1 and the latitude at point 2
${\lambda }_1{\lambda }_2$ is the Longitude at point 1 and longitude at point 2

python 代码实现

ng1,lat1 为开始点的经纬度
ng2,lat2 为结束点的经纬度

def geodistance(lng1,lat1,lng2,lat2):
    #lng1,lat1,lng2,lat2 = (120.12802999999997,30.28708,115.86572000000001,28.7427)
    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

Reference

公式参考
代码参考