根据经度纬度计算距离 帮我写一个程序

要编写一个根据经纬度计算两点间距离的程序，通常我们会使用球面三角学原理，比如地球上两点之间的大圆距离（也称为“Haversine”公式）。这里是一个简单的Python示例，使用了haversine库来简化计算：

from haversine import haversine, Unit

# 定义一个计算经纬度之间距离的函数
def calculate_distance(lat1, lon1, lat2, lon2):
    # 将角度转换为弧度
    lat1, lon1, lat2, lon2 = map(radians, [lat1, lon1, lat2, lon2])
    
    # 使用haversine函数计算距离
    distance = haversine((lat1, lon1), (lat2, lon2), unit=Unit.KILOMETERS)
    
    return distance

# 示例使用
lat1, lon1 = 39.9042, 116.4074  # 北京坐标
lat2, lon2 = 51.5074, -0.1278   # 伦敦坐标
distance = calculate_distance(lat1, lon1, lat2, lon2)
print(f"北京到伦敦的距离大约是 {distance:.2f} km.")