Nature风格的北半球投影图，几行代码就能绘制...

Python优雅地绘制北半球投影地图

经常在Nature正刊和子刊中看到北半球或南半球投影的球形地图，非常美观~

这种图ArcGIS往往不方便绘制

在子刊***Nat. Clim. Chang.***中，通过北极俯视投影表示NAO和AO的EOF分解。

Hamouda, M.E., Pasquero, C. & Tziperman, E. Decoupling of the Arctic Oscillation and North Atlantic Oscillation in a warmer climate. Nat. Clim. Chang. 11, 137–142 (2021). https://doi.org/10.1038/s41558-020-00966-8

如子刊Nature Geosci中表示了海冰的漂移，还添加了DEM底图，非常美观。

Darby, D., Ortiz, J., Grosch, C. et al. 1,500-year cycle in the Arctic Oscillation identified in Holocene Arctic sea-ice drift. Nature Geosci 5, 897–900 (2012). https://doi.org/10.1038/ngeo1629

接下来让我们复现一下试试

数据

以海温数据为例，ERA5月温度平均，是一张721*1440的0.25degTIF

import rioxarray as rxr
file_name='/home/mw/project/SST.tif' 
ds=rxr.open_rasterio(file_name)
data = ds[0]

如果你有其它nc数据，也可以通过xarray读取

投影设置

其实cartopy提供了一些投影使用