python 根据grb2文件出风速色斑图和风向标

因工作原因需要出气象图，下面直接上代码

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import shapefile
from matplotlib.path import Path
from matplotlib.patches import PathPatch
import xarray as xr
import numpy as np
import scipy.ndimage as ndimage
from matplotlib.font_manager import FontProperties
import matplotlib.colors as colors
from collections.abc import Iterable
from shapely.geometry import Point as ShapelyPoint
from shapely.geometry import Polygon as ShapelyPolygon

#处理地图白化函数
def shp2clip(originfig, ax, shpfile, regionlist):
    sf = shapefile.Reader(shpfile)
    vertices = []
    codes = []
    for shape_rec in sf.shapeRecords():
        pts = shape_rec.shape.points
        prt = list(shape_rec.shape.parts) + [len(pts)]
        for i in range(len(prt) - 1):
            for j in range(prt[i], prt[i + 1]):
                vertices.append((pts[j][0], pts[j][1]))
            codes += [Path.MOVETO]
            codes += [Path.L