给定两点经纬度坐标求距离的python代码

以下是一个示例代码，可以计算给定两点经纬度坐标之间的距离（单位为千米）：

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

def distance(lat1, lon1, lat2, lon2):
    """
    Calculate the great circle distance between two points 
    on the earth (specified in decimal degrees)
    """
    # 将十进制度数转化为弧度
    lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])

    # haversine公式
    dlon = lon2 - lon1 
    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

# 示例使用
print(distance(39.9042, 116.4074, 31.2304, 121.4737)) # 北京(39.9042, 116.4074)到上海(31.2304, 121.4737)的距离

该代码使用了Haversine公式来计算两点之间的距离。需要注意的是，这个公式假设地球是一个完美的球体，因此在实际使用中可能存在一定的误差。