python等值面追踪_Python兰伯特投影中国区域等值线图（含南海小地图）

import numpy as npimport xarray as xrimport cartopy.crs as ccrsimport matplotlib.pyplot as pltimport maskout2import maskoutimport lambert_ticksfrom scipy.interpolate import Rbffrom cartopy.io.shapereader import Readerfrom cartopy.mpl.gridliner import LATITUDE_FORMATTER, LONGITUDE_FORMATTERfrom matplotlib import rcParamsconfig={"font.family":'Times New Roman',"font.size":16,"mathtext.fontset":'stix'}rcParams.update(config)plt.rcParams['figure.figsize']=(15,15)nc_data = xr.open_dataset('F:/Rpython/lp28/data/r160.nc')olon = nc_data.coords["station_lon"]olat = nc_data.coords["station_lat"]stid = nc_data.coords["station_id"]stid_data = np.random.rand(len(olon)) * 500nlon = np.linspace(70, 140, 140)nlat = np.linspace(15, 55, 60)nlon, nlat = np.meshgrid(nlon, nlat)func = Rbf(olon, olat, stid_data, function='linear')stid_pre = fu