如何通过数据环绕处理消除Basemap地图投影的等值线不连续

问题根因 #

• 原始网格经度从-180°排到+180°，两个端点在地理上是同一个位置，但网格没有做实首尾衔接关系。Mollweide投影本身是全局椭圆投影，会自动识别这种经向边界，所以渲染不会断；但npstere、spstere这类圆形极投影不会自动处理跨±180°边界的网格连线，就会在圆形轮廓边（也就是你看到的0°纬度位置）把等值线、填充色拆成两段，出现明显断裂。
• 原代码里手动切lon[:,:-1]的操作反而丢了边界上的整列数据，让断裂问题更严重。

具体处理步骤 #

1. 对原始网格做经向环绕补全 #

沿经度方向拼接边界列数据，让网格形成经向闭环：

# 假设lon、lat、dat的形状为(纬度采样数, 经度采样数)
# 将经度=-180°的第一列数据，拼接到数组末尾，形成闭合环
lon_wrapped = np.concatenate([lon, lon[:, 0:1] + 360], axis=1)
lat_wrapped = np.concatenate([lat, lat[:, 0:1]], axis=1)
dat_wrapped = np.concatenate([dat, dat[:, 0:1]], axis=1)

注意：拼接时给末尾的经度值加360°，是为了和原始网格的经度值形成连续递增的序列，避免Basemap转换坐标时出现重复坐标点。

2. 用补全后的网格替换原绘制代码中的数组 #

• 北半球投影适配代码：

ax1 = plt.subplot(321)
m = Basemap(projection='npstere', lon_0=0, resolution='c', boundinglat=25, round=True)
# 全部使用环绕补全后的网格
m.contour(lon_wrapped, lat_wrapped, dat_wrapped, 8, colors='k', latlon=True, linewidths=0.5)
m.contourf(lon_wrapped, lat_wrapped, dat_wrapped,
            levels=np.linspace(0,10,100), vmin=0, vmax=10, extend='both',
            cmap=reversed_color_map, latlon=True)
# 极区投影的纬线刻度范围要和boundinglat匹配，避免绘制无效刻度
m.drawparallels(np.arange(25, 90, 30), labels=[0,0,0,0])
m.drawmeridians(np.arange(-180, 180, 30), labels=[0,0,0,0], zorder=100)

• 南半球投影适配代码：

plt.subplot(322)
m = Basemap(projection='spstere', lon_0=0, resolution='c', boundinglat=-25, round=True)
m.contour(lon_wrapped, lat_wrapped, dat_wrapped, 8, colors='k', latlon=True, linewidths=0.5)
m.contourf(lon_wrapped, lat_wrapped, dat_wrapped,
            levels=np.linspace(0,10,100), vmin=0, vmax=10, extend='both',
            cmap=reversed_color_map, latlon=True)
m.drawparallels(np.arange(-90, -25, 30), labels=[0,0,0,0])
m.drawmeridians(np.arange(0, 360, 30), labels=[0,0,0,0], zorder=100)

3. 缝隙兜底处理 #

如果补全后仍有细微渲染缝隙，可在contour/contourf方法中增加tri=True参数，强制将网格拆分为三角形单元计算等值线，消除边界处的渲染错位。

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