我正在努力通过不同的情节保持相同的颜色条范围。
例如,我有这些可视化:
使用此代码生成的内容:
def plot_contour(x_dim, y_dim, x_steps, y_steps, scalar_field, file_path):
plt.figure()
x, y = numpy.mgrid[-x_dim:x_dim/:x_steps*1j, -y_dim:y_dim:y_steps*1j]
cs = plt.contourf(x, y, scalar_field, zorder=1, extent=[-x_dim, x_dim, -y_dim, y_dim])
plt.colorbar(cs)
plt.savefig(file_path + '.png', dpi=Vc.dpi)
plt.close()
我希望能够比较两个字段,因此,我想对它们使用相同的颜色映射。
我的第一种方法是使用参数v_min
和v_max
,使用数据的最小值/最大值。
cs = plt.contourf(x, y, scalar_field, zorder=1, extent=[-x_dim, x_dim, -y_dim, y_dim], vmin=-1.00, vmax=1.05) # Manual setting to test
然后我得到了相同的颜色映射:
但我也想在图中显示相同的颜色条范围。我试着用
cb = plt.colorbar(cs)
cb.set_clim(vmin=-1.00, vmax=1.05)
没有成功。
这个完整的例子产生了相同的行为:
import matplotlib
import numpy as numpy
import matplotlib.cm as cm
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt
matplotlib.rcParams['xtick.direction'] = 'out'
matplotlib.rcParams['ytick.direction'] = 'out'
delta = 0.025
x = numpy.arange(-3.0, 3.0, delta)
y = numpy.arange(-2.0, 2.0, delta)
X, Y = numpy.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
Za = 10.0 * (Z2 - Z1)
Zb = 5.0 * (Z2 - Z1)
def bounds(scalar_fields):
"""
Get the bounds of a set of scalar_fields
:param scalar_fields : the scalar field set
:return: a set of normalized vector field components
"""
max_bound = -numpy.inf
min_bound = numpy.inf
for scalar_field in scalar_fields:
max_lim = numpy.max(scalar_field)
min_lim = numpy.min(scalar_field)
if max_lim > max_bound:
max_bound = max_lim
if min_lim < min_bound:
min_bound = min_lim
return min_bound, max_bound
def plot_contour(x_dim, y_dim, x_steps, y_steps, scalar_field, v_min, v_max, file_path):
plt.figure()
x, y = numpy.mgrid[-x_dim/2:x_dim/2:x_steps*1j, -y_dim/2:y_dim/2:y_steps*1j]
cs = plt.contourf(x, y, scalar_field, zorder=1, extent=[-x_dim/2.0, x_dim/2.0, -y_dim/2.0, y_dim/2.0],
vmin=v_min, vmax=v_max)
cb = plt.colorbar(cs)
plt.savefig(file_path + '.png')
plt.close()
v_min, v_max = bounds([Za, Zb])
x_dim = y_dim = 6
y_steps = x.shape[0]
x_steps = y.shape[0]
plot_contour(x_dim, y_dim, x_steps, y_steps, Za, v_min, v_max, 'Za')
plot_contour(x_dim, y_dim, x_steps, y_steps, Zb, v_min, v_max, 'Zb')
我怎么能这样做?
提前谢谢。
答案 0 :(得分:1)
如果您希望颜色条中的颜色与两个等高线图中的相同值对应,则您不仅需要控制颜色条,还需要控制等高线图中的级别。也就是说,为了比较图之间的相同水平,图应该具有相同的轮廓水平。这很容易做到。这是该情节的一个例子:
有两种方法:1)提前计算水平; 2)使用一个图中的水平来设置另一个图中的水平。我会做第二个,因为从中应该清楚如何做第一个(例如,使用levels = numpy.linspace(v_min, vmax, 10)
,但要明确,我不在这里使用它,但是我让mpl计算水平。)
首先,我在这里也使用:
Za = 10.0 * (Z2 - Z1)
Zb = 6.0 * (Z2 - Z1) # 6, rather than 5
然后,绘制:
def plot_contour(x_dim, y_dim, x_steps, y_steps, scalar_field, file_path, v_min, v_max, levels=None):
x, y = numpy.mgrid[-x_dim/2:x_dim/2:x_steps*1j, -y_dim/2:y_dim/2:y_steps*1j]
cs = plt.contourf(x, y, scalar_field, zorder=1, cmap=cm.jet, extent=[-x_dim/2.0, x_dim/2.0, -y_dim/2.0, y_dim/2.0], vmin=v_min, vmax=v_max, levels=levels)
plt.colorbar(cs)
return cs.levels
v_min, v_max = bounds([Za, Zb])
plt.figure()
plt.subplot(121)
levels = plot_contour(x_dim, y_dim, x_steps, y_steps, Za, 'Za', v_min, v_max)
plt.subplot(122)
plot_contour(x_dim, y_dim, x_steps, y_steps, Zb, 'Zb', v_min, v_max, levels=levels)
plt.show()