用cartopy绘制大圆路径

import cartopy.crs as ccrs
import matplotlib.pyplot as plt

start_lon, start_lat = 116.3974, 39.9092  # 起始点经纬度
end_lon, end_lat = -73.9787, 40.7488    # 终点经纬度

# 创建大圆路径样式
gc = ccrs.Geodetic()
path = ccrs.Geodesic().great_circle(start_lon, start_lat, end_lon, end_lat, 100)

# 创建地图画布
fig = plt.figure(figsize=(10, 6))
ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())

# 添加地图背景
ax.stock_img()

# 绘制大圆路径
ax.plot([start_lon, end_lon], [start_lat, end_lat], transform=gc, linewidth=2, color='red')

# 设置地图显示范围
ax.set_extent([-180, 180, -90, 90], crs=ccrs.PlateCarree())

# 添加经纬度网格线和标签
ax.gridlines(draw_labels=True, color='gray', alpha=0.5, linestyle='--')

# 添加起始点和终点标记
ax.plot(start_lon, start_lat, 'bo', transform=ccrs.PlateCarree(), label='起始点')
ax.plot(end_lon, end_lat, 'bo', transform=ccrs.PlateCarree(), label='终点')

# 添加标题和图例
ax.set_title('大圆路径示例')
ax.legend()

# 显示地图
plt.show()

pip install cartopy