话不多说,直接上代码
def haversine(lng1, lat1, lng2, lat2):
lng1, lat1, lng2, lat2 = map(radians, [lng1, lat1, lng2, lat2])
dlon = lng2 - lng1
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 * 1000
# 经度1,纬度1,经度2,纬度2 (十进制度数)计算两经纬度坐标点的距离
文章摘自https://blog.csdn.net/vernice/article/details/46581361