这种文章本来网上一搜一大把，但有的计算结果是错的，遂找到个计算正确的，记录一下。
参考原文：LBS 球面距离公式
#!/usr/bin/python3from math import sin, asin, cos, radians, fabs, sqrt EARTH_RADIUS=6371 # 地球平均半径，6371km def hav(theta): s = sin(theta / 2) return s * s def 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.3525231320157
1938.2154637477815
11640.132941713944
附自己做查距离的站点距离查询
