import math
def calculate_distance(lat1, lon1, lat2, lon2):
地球半径（千米）
R = 6371.0

# 将经纬度转换为弧度
lat1_rad = math.radians(lat1)
lon1_rad = math.radians(lon1)
lat2_rad = math.radians(lat2)
lon2_rad = math.radians(lon2)

# 计算差值
dlat = lat2_rad - lat1_rad
dlon = lon2_rad - lon1_rad

# haversine 公式
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
示例：北京和上海的经纬度
beijing = (39.9042, 116.4074)
shanghai = (31.2304, 121.4737)
dist = calculate_distance(beijing[0], beijing[1], shanghai[0], shanghai[1])
print(f"两地距离约为: {dist:.2f} 千米")

# 安装命令：pip install geopy
from geopy.distance import geodesic
坐标格式：(纬度, 经度)
point1 = (39.9042, 116.4074)  # 北京
point2 = (31.2304, 121.4737)  # 上海
计算距离
distance = geodesic(point1, point2).kilometers
print(f"两地距离为: {distance:.2f} 千米")