根据两点经纬度计算球面上大圆距离

原理：哈弗赛恩公式

import math

def calculate_distance(lat1, lon1, lat2, lon2):
    R = 6371  # 地球半径，单位为千米
    rad = lambda angle: angle * math.pi / 180  # 将角度转换为弧度

    lat1_rad = rad(lat1)
    lon1_rad = rad(lon1)
    lat2_rad = rad(lat2)
    lon2_rad = rad(lon2)

    dLat = lat2_rad - lat1_rad
    dLon = lon2_rad - lon1_rad

    a = math.sin(dLat / 2) ** 2 + math.cos(lat1_rad) * math.cos(lat2_rad) * math.sin(dLon / 2) ** 2
    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1 - a))

    distance = R * c  # 返回单位为千米的距离
    return distance

# 示例用法
lat1 = 34.444444  # 北京的经纬度
lon1 = 112.257778
lat2 = 34.419  # 上海的经纬度
lon2 = 111.876

distance = calculate_distance(lat1, lon1, lat2, lon2)
print(distance)  # 输出直线距离（单位：千米）