from 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 distance lon1, 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) print(get_distance_hav(113.6917100000,34.8279700000,113.674250000, 34.8087700000))