原理：哈弗賽恩公式

import math

def calculate_distance(lat1, lon1, lat2, lon2):
    R = 6371  # 地球半徑，單位為千米
    rad = lambda angle: angle * math.pi / 180  # 將角度轉(zhuǎn)換為弧度

    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  # 北京的經(jīng)緯度
lon1 = 112.257778
lat2 = 34.419  # 上海的經(jīng)緯度
lon2 = 111.876

distance = calculate_distance(lat1, lon1, lat2, lon2)
print(distance)  # 輸出直線距離（單位：千米）