看过很多大佬绘制过中国地图，有R-ggplot的，也有python-geopandas等，还有其他cartopy等的作图包，我就不一一举例了，今天用proplot这个包，这个包的简介可以参考链接proplot，主要的知识点如下：
-
proplot绘制中国地图 -
ax.inset添加南海子图
导入需要的包
import os
import numpy as np
import xarray as xr
import proplot as plot
import cartopy.crs as ccrs
from cartopy.io.shapereader import Reader
import maskout
shp_path = 'D:/China_shp'
导入数据
- 以历史气温数据为例（historical）
##### 数据读取
ds = xr.open_dataset('*historical*.nc')
tas = ds['tas'][75] # 随意选择的年份
#截取一段时间内的降水作为绘图的数据源
tas -= 273.15
lon = ds.lon
lat = ds.lat
lon_range = lon[(lon > 70) & (lon < 140)]
lat_range = lat[(lat > 0) & (lat < 60)]
tas_region = tas.sel(lon=lon_range, lat=lat_range)
levels = plot.arange(-5,30,5)
主图&白化[1]
fig, axs = plot.subplots(width=6,proj=ccrs.PlateCarree())
# First axis
ax = axs[0]
m = ax.contourf(tas_region,
cmap='CoolWarm',
levels=levels,
extend='both')
# 搜白化教程
# http://bbs.06climate.com/forum.php?mod=viewthread&tid=42437&highlight=%CD%EA%C3%C0%B0%D7%BB%AF
con_mask = maskout.shp2clip(m, ax, r'E:\maskout\country1.shp', ['China'])
ax.add_geometries(Reader(os.path.join(shp_path,
'bou2_4l.shp')).geometries(),
ccrs.PlateCarree(),
facecolor='none',
edgecolor='black',
linewidth=0.5)
ax.format(
coast=False, labels=True, borders=False,
latlim=(10, 60),
lonlim=(70, 140),
lonlines=10, latlines=10,
)
[1]地图白化
添加南海子图
# Create inset axes representing a "zoom-in"
# https://proplot.readthedocs.io/en/latest/insets_panels.html
iax = ax.inset([0.83,0,0.175,0.3]) # 不同的主图参数，这个参数设置也不一样，自己摸索
n = iax.contourf(tas_region,
cmap='CoolWarm',
levels=levels,
extend='both')
con_mask = maskout.shp2clip(n, iax, r'E:\maskout\country1.shp', ['China'])
iax.add_geometries(Reader(os.path.join(shp_path,
'bou2_4l.shp')).geometries(),
ccrs.PlateCarree(),
facecolor='none',
edgecolor='black',
linewidth=0.5)
iax.format(
coast=False, labels=False, borders=False,
latlim=(0, 26),
lonlim=(105, 125),
lonlines=10, latlines=10,
)
fig.colorbar(m, loc='r')
出图
China
总结
Proplot这个包太强大了，在这里就不细说了！共同学习，那就这样吧！(时间紧，任务重/奋斗脸)，写这个就当是抛砖引玉吧！希望大家多多指教！
