这可以使用shapely形状和union操作来完成。这个想法是为掩码为真的每个单元格创建一个正方形（或矩形等），然后将相邻的正方形组合在一起。然后可以绘制这些几何图形，或使用这些几何图形创建matplotlib多边形。下图以红色显示轮廓（根据需要带有楼梯边缘）：

import shapely.geometry
import shapely.ops

geoms = []
for yidx, xidx in zip(*np.where(mask)):
    geoms.append(shapely.geometry.box(xbound[xidx], ybound[yidx], xbound[xidx+1], ybound[yidx+1]))
full_geom = shapely.ops.unary_union(geoms)

for geom in full_geom.geoms:
    plt.plot(*geom.exterior.xy, linewidth=4, color='r')
plt.imshow(mask, origin='lower', extent=(0, 10, 0, 10))

很难用颜色/图案填充每个形状，因此可以使用matplotlib多边形：

from matplotlib.patches import Polygon

polygons = []
for geom in full_geom.geoms:
    p = Polygon(list(zip(*geom.exterior.xy)), closed=True, facecolor='none', edgecolor='r', linestyle='-', hatch='//', linewidth=4)
    polygons.append(p)

plt.figure()
axes = plt.gca()
for p in polygons:
    axes.add_patch(p)
plt.imshow(mask, origin='lower', extent=(0, 10, 0, 10))

最后，可以将原始shapely多边形直接添加到制图图中：

import cartopy.crs as ccrs

fig, ax = plt.subplots(subplot_kw={'projection': ccrs.PlateCarree()})
ax.add_geometries(
    full_geom.geoms,
    crs=ccrs.PlateCarree(),
    facecolor='white',
    edgecolor='red',
    linewidth=4,
    hatch='//',
)
ax.coastlines()
ax.gridlines(draw_labels=True, dms=True, x_inline=False, y_inline=False)
ax.set_xlim((0, 10))
ax.set_ylim((0, 10))