Matplotlib中的叠加轮廓图

Question

我需要比较两组二维分布。

当我使用matplotlib.pyplot.contourf并叠加图形时，每个等高线图的背景颜色填充整个图形空间。有没有办法使每个等高线图的最低等高线透明，以便更容易看到每个等高线的中心?

import numpy as np
import matplotlib.pyplot as plt
from matplotlib import cm 
import scipy.stats as st

def make_cloud(x, y, std, n=100):
    x = np.random.normal(x, std, n)
    y = np.random.normal(y, std, n)
    return np.array(zip(x, y))

def contour_cloud(x, y, cmap):
    xmin, xmax = -4, 4
    ymin, ymax = -4, 4

    xx, yy = np.mgrid[xmin:xmax:100j, ymin:ymax:100j]
    positions = np.vstack([xx.ravel(), yy.ravel()])
    values = np.vstack([x, y])
    kernel = st.gaussian_kde(values)
    f = np.reshape(kernel(positions).T, xx.shape)

    plt.contourf(xx, yy, f, cmap=cmap, alpha=0.5)


cloud1 = make_cloud(-1, 1, 1)
cloud2 = make_cloud(1, -1, 1)

plt.scatter(x=cloud1[:,0], y=cloud1[:,1])
plt.scatter(x=cloud2[:,0], y=cloud2[:,1], color='red')

fig = plt.gcf()
ax = plt.gca()

contour_cloud(x=cloud1[:, 0], y=cloud1[:, 1], cmap=cm.Blues)
contour_cloud(x=cloud2[:, 0], y=cloud2[:, 1], cmap=cm.Reds)

- Chris

Answer 1

有几个控件是您需要查看的，以便使用 contourf。您可以手动更改不同的级别，并更改颜色映射的上/下限规格。默认情况下，最低级别（或最高级别以上）面积的填充似乎是透明的。

因此，实现您想要的最简单的方法是手动指定级别，并指定这些级别使得存在低于最低级别的点，但不存在高于最高级别的点。

如果您将以下内容替换：

plt.contourf(xx, yy, f, cmap=cmap, alpha=0.5)

使用：

step = 0.02
m = np.amax(f)
levels = np.arange(0.0, m, step) + step
plt.contourf(xx, yy, f, levels, cmap=cmap, alpha=0.5)

生成的图片如下所示：

有关超出/低于颜色映射值时的行为的更多细节，请参见此处。

- Douglas Dawson