从该代码处绘制wrf中的土地利用报错内容及其解决方法

1.报错内容：

微信公众平台 (qq.com)https://mp.weixin.qq.com/s/Cn0vhvfroVADPnT237LXNw

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
Cell In[9], line 7
      4     start(file_in,shp_path)
      6 if __name__ == "__main__":
----> 7     main()

Cell In[9], line 4, in main()
      2 file_in = "F:/pythonplot/geo_em.d01.nc"
      3 shp_path = "F:/毕业数据处理/中国shp/3省市县行政区划界线/chinaProvince.shp"
----> 4 start(file_in,shp_path)

Cell In[5], line 43, in start(file_in, shp_path)
     41 ax.xaxis.set_major_formatter(LONGITUDE_FORMATTER)
     42 ax.yaxis.set_major_formatter(LATITUDE_FORMATTER)
---> 43 lambert_xticks(ax, xticks)
     44 lambert_yticks(ax, yticks)
     46 # 叠加shp

Cell In[8], line 19, in lambert_xticks(ax, ticks)
     17 te = lambda xy: xy[0]
     18 lc = lambda t, n, b: np.vstack((np.zeros(n) + t, np.linspace(b[2], b[3], n))).T
---> 19 xticks, xticklabels = _lambert_ticks(ax, ticks, 'bottom', lc, te)
     20 ax.xaxis.tick_bottom()
     21 ax.set_xticks(xticks)

Cell In[8], line 35, in _lambert_ticks(ax, ticks, tick_location, line_constructor, tick_extractor)
     33 def _lambert_ticks(ax, ticks, tick_location, line_constructor, tick_extractor):
     34     """Get the tick locations and labels for an axis of a Lambert Conformal projection."""
---> 35     outline_patch = sgeom.LineString(ax.outline_patchs.get_path().vertices.tolist())
     36     axis = find_side(outline_patch, tick_location)
     37     n_steps = 30

AttributeError: 'GeoAxes' object has no attribute 'outline_patch'

2.解决方法： 

源代码：

def _lambert_ticks(ax, ticks, tick_location, line_constructor, tick_extractor):
    """Get the tick locations and labels for an axis of a Lambert Conformal projection."""
    outline_patch = sgeom.LineString(ax.outline_patch.get_path().vertices.tolist())
    axis = find_side(outline_patch, tick_location)
    n_steps = 30
    extent = ax.get_extent(ccrs.PlateCarree())
    _ticks = []
    for t in ticks:
        xy = line_constructor(t, n_steps, extent)
        proj_xyz = ax.projection.transform_points(ccrs.Geodetic(), xy[:, 0], xy[:, 1])
        xyt = proj_xyz[..., :2]
        ls = sgeom.LineString(xyt.tolist())
        locs = axis.intersection(ls)
        if not locs:
            tick = [None]
        else:
            tick = tick_extractor(locs.xy)
        _ticks.append(tick[0])
    # Remove ticks that aren't visible:
    ticklabels = copy(ticks)
    while True:
        try:
            index = _ticks.index(None)
        except ValueError:
            break
        _ticks.pop(index)
        ticklabels.pop(index)
    return _ticks, ticklabels

改变代码

def _lambert_ticks(ax, ticks, tick_location, line_constructor, tick_extractor):
    """Get the tick locations and labels for an axis of a Lambert Conformal projection."""
    outline_patch = sgeom.LineString(ax.patch.get_path().vertices.tolist())
    axis = find_side(outline_patch, tick_location)
    n_steps = 30
    extent = ax.get_extent(ccrs.PlateCarree())
    _ticks = []
    for t in ticks:
        xy = line_constructor(t, n_steps, extent)
        proj_xyz = ax.projection.transform_points(ccrs.Geodetic(), xy[:, 0], xy[:, 1])
        xyt = proj_xyz[..., :2]
        ls = sgeom.LineString(xyt.tolist())
        locs = axis.intersection(ls)
        if not locs:
            tick = [None]
        else:
            tick = tick_extractor(locs.xy)
        _ticks.append(tick[0])
    # Remove ticks that aren't visible:
    ticklabels = copy(ticks)
    while True:
        try:
            index = _ticks.index(None)
        except ValueError:
            break
        _ticks.pop(index)
        ticklabels.pop(index)
    return _ticks, ticklabels