python利用地图两个点的经纬度计算两点间距离

这种文章本来网上一搜一大把，但有的计算结果是错的，遂找到个计算正确的，记录一下。

参考原文：LBS 球面距离公式

#!/usr/bin/python3from math import sin, asin, cos, radians, fabs, sqrtEARTH_RADIUS=6371           # 地球平均半径，6371kmdef hav(theta):s = sin(theta / 2)return s * sdef get_distance_hav(lat0, lng0, lat1, lng1):"用haversine公式计算球面两点间的距离。"# 经纬度转换成弧度lat0 = radians(lat0)lat1 = radians(lat1)lng0 = radians(lng0)lng1 = radians(lng1)dlng = fabs(lng0 - lng1)dlat = fabs(lat0 - lat1)h = hav(dlat) + cos(lat0) * cos(lat1) * hav(dlng)distance = 2 * EARTH_RADIUS * asin(sqrt(h))return distancelon1,lat1 = (22.599578, 113.973129) #深圳野生动物园(起点）lon2,lat2 = (22.6986848, 114.3311032) #深圳坪山站 (百度地图测距：38.3km)d2 = get_distance_hav(lon1,lat1,lon2,lat2)print(d2)lon2,lat2 = (39.9087202, 116.3974799) #北京天安门(1938.4KM)d2 = get_distance_hav(lon1,lat1,lon2,lat2)print(d2)lon2,lat2 =(34.0522342, -118.2436849) #洛杉矶(11625.7KM)d2 = get_distance_hav(lon1,lat1,lon2,lat2)print(d2)

输出结果很接近：

38.35252313201571938.215463747781511640.132941713944

附自己做查距离的站点距离查询