Python 根据经纬度坐标计算位置

高德API 经纬度地址查询：https://lbs.amap.com/tools/picker

示例代码 

from math import radians, cos, sin, asin, sqrt


# 公式计算两点间距离（m）
def geodistance(lng1, lat1, lng2, lat2):
    # lng1,lat1,lng2,lat2 = (116.37, 40.04, 116.43, 40.08)
    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


if __name__ == "__main__":
    t = geodistance(116.37, 40.04, 116.43, 40.08)
    print t # 6.772KM

 

#调用geopy包中的方法

from geopy.distance import geodesic

print(geodesic((30.28708,120.12802999999997), (28.7427,115.86572000000001)).m) #计算两个坐标直线距离
print(geodesic((30.28708,120.12802999999997), (28.7427,115.86572000000001)).km) #计算两个坐标直线距离
# 返回 447.2497993542003 千米


# 南昌：华东交通大学（120.12802999999997,30.28708）
# 杭州：浙江工商大学（115.86572000000001,28.7427）