import math
<p>def calculate_distance(lat1, lon1, lat2, lon2):</p><h1>地球半径（千米）</h1><pre class="brush:python;toolbar:false;">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
<p>from geopy.distance import geodesic</p><h1>坐标格式：(纬度, 经度)</h1><p>point1 = (39.9042, 116.4074)  # 北京
point2 = (31.2304, 121.4737)  # 上海</p><h1>计算距离</h1><p>distance = geodesic(point1, point2).kilometers
print(f"两地距离为: {distance:.2f} 千米")</p>