Отрезать pcolormesh за пределами круглого set_boundary в Cartopy

Я использую Cartopy для своих полярных исследований и хотел бы обрезать круговую границу вокруг моих данных, которые я рисую в проекции NorthPolarStereo(). Я использую set_extent, чтобы указать, с какой широты я хотел бы нанести свои данные, и использую set_boundary для создания круговой границы, как объяснено в галерее. Затем я использую matplotlib.pyplot.pcolormesh для построения фактических данных. Однако, скажем, я использую set_extent для определения минимальной широты 55 градусов, некоторые из моих данных ниже 55 градусов все еще отображаются за пределами моего set_boundary. Как мне обрезать эти данные?

map_crs = ccrs.NorthPolarStereo(central_longitude=0.0, globe=None)

# Build axes
fig     = plt.figure()
ax      = plt.axes(projection=map_crs)

plotfield = ax.pcolormesh(lons, lats, data, transform=ccrs.PlateCarree())
ax.set_extent((-180, 180, 55, 90), crs=ccrs.PlateCarree())
gl = ax.gridlines()

# Circular clipping
theta = np.linspace(0, 2*np.pi, 400)
center, radius = [0.5, 0.5], 0.5
verts = np.vstack([np.sin(theta), np.cos(theta)]).T
circle = mpath.Path(verts * radius + center)
ax.set_boundary(circle, transform=ax.transAxes)