计算两个纬度 - 经度点之间的距离？(Haversine配方)

a = sin²((纬度2 - 纬度1) / 2) + cos(纬度1) * cos(纬度2) * sin²((经度2 - 经度1) / 2)
c = 2 * atan2(√a, √(1-a))
d = R * c

import math

def haversine(lat1, lon1, lat2, lon2):
    R = 6371  # 地球半径，单位为千米
    # 将经纬度的度数转换为弧度
    lat1, lon1, lat2, lon2 = map(math.radians, [lat1, lon1, lat2, lon2])
    # 计算Haversine公式中的部分变量
    a = math.sin((lat2 - lat1) / 2) ** 2 + math.cos(lat1) * math.cos(lat2) * math.sin((lon2 - lon1) / 2) ** 2
    c = 2 * math.asin(math.sqrt(a))
    d = R * c
    return d