在使用cartopy绘制极地投影(北极或者南极地区时),会出现不少问题,最早发现的极地投影无法添加坐标系问题已在新版cartopy20.0+里解决,具体可参见我的博客:artopy20.0+解决python极地投影问题
但是,新版的cartopy并没有解决绘制等值线与等值线填充投影重叠导致绘图错误的问题,我花了一天时间,终于在GITHUB上找到了解决办法: Plotting curvilinear with wrapping issue: suggestion on how to do it,这个方法可以很好地解决cartopy绘制极地投影的问题。下面两张图就是修改前后的对比:绘制图形为2018年春季北极地区500hPa风场与位势高度场,可以看到未修改前的古怪多边形和错误的风场,以及修改后明显该方法对于基地投影的改善相当有效。
具体代码除了参考帖子,我也写了博客:python cartopy极地投影重叠问题解决,有兴趣可以观看,也给自己博客打个小广告:我会不定期分享一些个人经验(包括代码、模式debug和一些文献方法阅读在博客里,欢迎大家围观指正)。
希望能帮助到其他被cartopy极地投影折磨的人。