Cartopy绘图：绘制盒状等值线的pcolor替代方案及pcolormesh维度不兼容问题排查

一、报错原因分析 #

你遇到的TypeError: Dimensions of C (15, 21) are incompatible with X (15) and/or Y (15)，本质是**pcolor/pcolormesh的坐标维度规则和你的数据不匹配**。

这类块图的绘制逻辑是：每个数据点C[i,j]对应一个由坐标点(X[i,j], Y[i,j])、(X[i+1,j], Y[i,j])、(X[i,j+1], Y[i,j+1])、(X[i+1,j+1], Y[i+1,j+1])围成的矩形方块。所以它要求：

• 如果C是(M, N)的二维数组，那么X和Y的维度必须是(M+1, N+1)（二维网格），或者一维的X长度为N+1、Y长度为M+1。

看你的代码：

• 数据df是(15,21)，对应15行纬度、21列经度的网格中心值
• 你用meshgrid生成的lon2d和lat2d是(15,21)，这是网格中心的坐标，不是网格边界的坐标，维度比要求少了1，自然报错。

二、解决方法 #

最直接的方式是生成网格边界坐标，替换原来的中心坐标：

# 假设你原来的lon、lat是一维数组（对应网格中心）
# 计算网格边界：每个中心向两侧扩展0.125（步长0.25的一半）
lon_edges = np.linspace(lon[0]-0.125, lon[-1]+0.125, len(lon)+1)
lat_edges = np.linspace(lat[0]-0.125, lat[-1]+0.125, len(lat)+1)
# 生成边界的二维网格
lon2d_edges, lat2d_edges = np.meshgrid(lon_edges, lat_edges)

# 用边界网格绘图
plot = plt.pcolor(lon2d_edges, lat2d_edges, df, cmap='jet', transform=ccrs.PlateCarree())

另外，如果你用xarray处理数据，可以直接调用xarray的plot.pcolormesh，它会自动处理坐标维度，省去手动生成meshgrid的麻烦：

# 把你的data_mean转换成xarray的DataArray，带上lon、lat坐标
da = xr.DataArray(data_mean['mean'].values.reshape(15,21), 
                  dims=['latitude', 'longitude'],
                  coords={'latitude': lat, 'longitude': lon})
# 直接绘图
da.plot.pcolormesh(ax=ax, cmap='jet', transform=ccrs.PlateCarree())

三、绘制盒状块图的替代方案 #

除了pcolor/pcolormesh，还有这些常用选项：

• plt.imshow：适用于规则网格，通过extent参数绑定地理坐标，设置origin='lower'匹配地理坐标系的方向：

  plt.imshow(df, cmap='jet', origin='lower', 
             extent=[lon[0]-0.125, lon[-1]+0.125, lat[0]-0.125, lat[-1]+0.125],
             transform=ccrs.PlateCarree())
  

• xarray内置绘图：如上面提到的，xarray会自动关联坐标和投影，尤其适合处理NetCDF数据的场景，代码更简洁。
• plt.matshow：和imshow类似，专注于矩阵可视化，同样通过extent调整地理范围，适合小范围规则网格的快速绘制。

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