Python气象绘图——一页多图以及不同投影坐标系

一、作业要求

一页多图，利用所给数据绘制：1. Robinson投影的全球气温图；2.极地投影的50°N以北区域的位势高度图；3.兰伯特投影的东亚地区(70°E-140°E, 15°N-55°N)气温图；4.PlateCarree投影的欧亚地区(50°W-140°E, 0°-90°N)风场图。需要设置colorbar，经纬度坐标，参考矢量等参数。每个子图用a，b，c，d标注

二、数据来源

数据来源于欧洲中期天气预报中心（ECMWF）的再分析数据（ERA5）。数据包括位势高度，气温以及纬向风和经向风。

数据读取以及数据信息。

# 数据读取
filename = r'ERA5_GP_T_UV.nc'
data = xr.open_dataset(filename)
hgt, temp = data.z, data.t
u, v = data.u, data.v
lon = data.longitude.values
lat = data.latitude.values
print(data)

ECMWF再分析数据的下载网址如下：

https://cds.climate.copernicus.eu/cdsapp#!/search?type=dataset

 三、代码

import numpy as np
import matplotlib.pyplot as plt
import xarray as xr
from matplotlib import gridspec
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import cartopy.mpl.ticker as cticker
import matplotlib.path as mpath
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER

plt.rcParams['font.sans-serif']= ['Microsoft YaHei'] # 设置“微软雅黑”，图上显示出中文
plt.rcParams['