# 引用部分
import numpy as np
import pandas as pd
from scipy.interpolate import Rbf # 径向基函数 : 将站点信息插到格点上 用于绘制等值线
import matplotlib.pyplot as plt
import matplotlib.colors as colors
import matplotlib as mpl
import cartopy.crs as ccrs
import cartopy.io.shapereader as shpreader
from cartopy.mpl.ticker import LongitudeFormatter,LatitudeFormatter
import maskout #只显示某个地区
plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文
plt.rcParams['axes.unicode_minus']=False #用来正常显示负号
# 数据准备
df = pd.read_excel('工作簿1.xlsx') #读取Excel
df.columns = ['stid','lon','lat','height','rain']
# df = df[(df['rain']>0)] #筛选坐标在某个范围内,由于数据量少,直接在EXCEL中做处理 &(df['lat']>35)&(df['lon']>115)&(df[