在Matplotlib中叠加轮廓图

时间:2017-03-04 23:06:43

标签: matplotlib contour

我需要比较2组的2维分布。

当我使用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)

enter image description here

1 个答案:

答案 0 :(得分:4)

您需要为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)

生成如下图像: enter image description here

有关色域图值上/下的值行为的详细信息,请参阅here