如何调整Basemap使完整大圆航线显示在地图中央？

为工作文档开发了一款可实时展示两座城市间范围的程序，使用Basemap绘制新加坡至旧金山的大圆航线时，当前地图无法完整显示整条航线，需要调整地图参数，让完整的大圆航线显示在地图中央。

当前地图仅能显示部分航线：新加坡位于地图右侧，旧金山位于左侧，中间的大圆航线被地图边界截断。

问题根源是当前Basemap的投影中心和边界范围未匹配大圆航线的覆盖区域，可通过以下调整解决：

• 调整地图中心：以新加坡和旧金山的经纬度中点作为地图中心（处理跨180度经线的情况），确保航线居中显示。
• 修改边界范围：根据航线的经纬度覆盖范围，设置合适的地图边界参数，确保整条航线被包含在地图内。
• 修复标记坐标：更新滑块回调函数，将原始经纬度转换为Basemap的投影坐标后再设置标记位置，避免标记偏移。

from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.widgets import Slider

# 创建画布和地图轴
fig = plt.figure()
ax = fig.add_axes([0.1, 0.3, 0.8, 0.6])

# 定义城市坐标
sglat = 1.35
sglon = 103.82
sflat = 37.77
sflon = -122.42

# 计算地图中心经纬度（处理跨180度经线的情况）
center_lon = (sglon + sflon) / 2
if abs(sglon - sflon) > 180:
    center_lon += 180
center_lat = (sglat + sflat) / 2

# 初始化Basemap，调整参数适配航线
m = Basemap(llcrnrlon=center_lon - 120, llcrnrlat=-30, 
            urcrnrlon=center_lon + 120, urcrnrlat=50,
            rsphere=(6378137.00, 6356752.3142),
            resolution='l', projection='merc',
            lat_0=center_lat, lon_0=center_lon, lat_ts=20.)

# 计算大圆航线的点集
num_points = 1000
lons, lats = m.gcpoints(sglon, sglat, sflon, sflat, num_points)

# 绘制大圆航线、海岸线和大陆填充
m.drawgreatcircle(sglon, sglat, sflon, sflat, linewidth=2, color='b')
m.drawcoastlines()
m.fillcontinents()

# 绘制经纬线
m.drawparallels(np.arange(-90, 90, 20), labels=[1, 1, 0, 1])
m.drawmeridians(np.arange(-180, 180, 30), labels=[1, 1, 0, 1])

ax.set_title('新加坡至旧金山大圆航线')

# 创建滑块控件
slider_ax = fig.add_axes([0.25, 0.1, 0.5, 0.03])
slider = Slider(slider_ax, '位置', 0, 1, valinit=0)

# 初始化标记点
marker, = ax.plot([], [], marker='o', markersize=8, color='r')

# 更新标记位置的回调函数
def update_marker(val):
    position = int(val * (num_points - 1))
    # 将经纬度转换为地图投影坐标
    x, y = m(lons[position], lats[position])
    marker.set_data([x], [y])
    fig.canvas.draw_idle()

# 绑定滑块事件
slider.on_changed(update_marker)

plt.show()

1. 地图中心优化：计算两个城市的经纬度中点，当经度差超过180度时调整中心经度，确保地图聚焦在太平洋方向的航线上。
2. 边界范围调整：设置经度范围为中心经度±120度，纬度范围从-30到50度，完整覆盖新加坡到旧金山的大圆航线区域。
3. 标记坐标修复：原代码直接使用原始经纬度更新标记，这会导致标记位置错误，现在通过m(lon, lat)转换为投影坐标后再设置，确保标记准确跟随航线。

内容的提问来源于stack exchange，提问作者M K