如何隐藏底图上特定区域的等高线/数据

Question

我正在使用一些气象数据在底图上绘制等高线。我之前完成的完整工作示例代码在这里How to remove/omit smaller contour lines using matplotlib。一切正常，我不抱怨等高线图。但是有一种特殊情况，我必须隐藏底图上特定区域（不规则纬度和经度）上的所有等高线。

我能想到的唯一可能的解决方案是在所需区域上绘制多边形线并填充与底图相同的颜色。经过大量搜索，我找到了这个链接How to draw rectangles on a Basemap（下面的代码）

from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.patches import Polygon

def draw_screen_poly( lats, lons, m):
    x, y = m( lons, lats )
    xy = zip(x,y)
    poly = Polygon( xy, facecolor='red', alpha=0.4 )
    plt.gca().add_patch(poly)

lats = [ -30, 30, 30, -30 ]
lons = [ -50, -50, 50, 50 ]

m = Basemap(projection='sinu',lon_0=0)
m.drawcoastlines()
m.drawmapboundary()
draw_screen_poly( lats, lons, m )

plt.show()

它似乎部分工作。但是，我想绘制一个不规则的区域。

感谢任何解决方案。

编辑：1

我已经明白问题出在哪里了。似乎在多边形区域内填充的任何颜色（facecolor）都不会使它隐藏在下面的任何东西。无论alpha 使用与否，它始终只是透明的。为了说明这个问题，我裁剪了具有所有三个区域的图像，即。等高线、底图区域和多边形区域。多边形区域用红色填充，但如您所见，等高线始终可见。我在上面的代码中使用的特定行是：-

poly = Polygon(xy, facecolor='red', edgecolor='b')

因此问题不在于上面的代码。多边形填充似乎有问题。但仍然没有解决这个问题。生成的图像（裁剪后的图像）如下（在附件图像下方查看我的第二次编辑）：-

编辑 2： 从这个与我有类似要求的http://matplotlib.1069221.n5.nabble.com/Clipping-a-plot-inside-a-polygon-td41950.html 中获取线索，我可以删除一些数据。但是，删除的数据仅来自多边形区域外部而不是内部。这是我从中获得线索的代码：-

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.patches import RegularPolygon

data = np.arange(100).reshape(10, 10)
fig = plt.figure()
ax = fig.add_subplot(111)
ax.contourf(data)
poly = RegularPolygon([ 0.5,  0.5], 6, 0.4, fc='none', 
                      ec='k', transform=ax.transAxes)
for artist in ax.get_children():
    artist.set_clip_path(poly)

现在我的问题是使用什么命令来删除多边形区域内的数据？

Answer 1

没有注意到有人对此提出索赔，所以我可能只是给出已经提出的解决方案here。您可以修改 zorder 以隐藏多边形后面的东西：

import matplotlib
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt

matplotlib.rcParams['xtick.direction'] = 'out'
matplotlib.rcParams['ytick.direction'] = 'out'

delta = 0.025
x = np.arange(-3.0, 3.0, delta)
y = np.arange(-2.0, 2.0, delta)
X, Y = np.meshgrid(x, y)
Z1 = mlab.bivariate_normal(X, Y, 1.0, 1.0, 0.0, 0.0)
Z2 = mlab.bivariate_normal(X, Y, 1.5, 0.5, 1, 1)
# difference of Gaussians
Z = 10.0 * (Z2 - Z1)


# Create a simple contour plot with labels using default colors.  The
# inline argument to clabel will control whether the labels are draw
# over the line segments of the contour, removing the lines beneath
# the label
fig = plt.figure()
ax = fig.add_subplot(111)
CS = plt.contour(X, Y, Z,zorder=3)
plt.clabel(CS, inline=1, fontsize=10)
plt.title('Simplest default with labels')

rect1 = matplotlib.patches.Rectangle((0,0), 2, 1, color='white',zorder=5)

ax.add_patch(rect1)

plt.show()

，结果是：