python气象数据处理与绘图_Python气象数据处理与绘图(6)：极地投影绘图

当我们研究高纬(极地)区域时，往往需要极地投影。在前文中通过EOF介绍了圆柱投影，实际上极地投影的绘制仅需改变其中几个参数便可以实现。

还是先上图：

图1 极地投影填色图

填色图的意义并不重要(我之前研究求过的两个时间序列与海冰密集度的相关系数分布)，接下来分享一下绘图的经验和代码。

proj =ccrs.NorthPolarStereo(central_longitude=90)

#在圆柱投影中proj = ccrs.PlateCarree(central_longitude=xx)

leftlon, rightlon, lowerlat, upperlat = (-180,180,60,90)

#仅画60°E-90°E部分

img_extent = [leftlon, rightlon, lowerlat, upperlat]

fig1 = plt.figure(figsize=(12,8))

#以下我仅展示了左半部分，右半部分基本一致，在此省略

f1_ax1 = fig3.add_axes([0.1, 0.1, 0.5, 0.5],projection = ccrs.NorthPolarStereo())

#注意此处添加了projection = ccrs.NorthPolarStereo()，指明该axes为北半球极地投影

f1_ax1.gridlines()

f1_ax1.set_extent(img_extent, ccrs.PlateCarree())

#通过圆柱投影的范围限制地图范围，这样设置地图参数较为方便

f