我正在制作2D numpy meshgrid的颜色图:
X, Y = np.meshgrid(fields, frequencies)
cs = ax.contourf(X, Y, fields_freqs_abs_grid, cmap="viridis", N=256)
fields_freqs_abs_grid中按颜色绘制的值已经过对数缩放。
python的matplotlib生成的色彩图很粗糙 - 即使我使用" N = 256"它也会缩放8种颜色。对于RGB像素的数量。将N增加到2048并没有改变任何东西。在相同数据上使用MatLab语言的图可生成具有明显更高颜色分辨率的色彩图。如何增加Python中映射的颜色数量?
谢谢!
答案 0 :(得分:2)
Warren Weckesser's comments绝对有效,可以为您提供高分辨率图像。我在下面的例子中实现了他的想法。
关于使用contourf()
,我不确定这是否是版本相关问题,但在最新版本中,
N
contourf()
没有kwarg。
正如您在文档中看到的那样,您希望使用N
作为arg(语法:contourf(X,Y,Z,N)
)来指定要绘制的级别数而不是RGB像素数。 contourf()
绘制填充的轮廓,分辨率取决于要绘制的级别数。您的N=256
将不会执行任何操作,contourf()
会自动选择7 levels。
以下代码是从官方example修改的,用不同的N
比较分辨率。如果存在版本问题,此代码会使用python 3.5.2; matplotlib 1.5.3
:
import numpy as np
import matplotlib.pyplot as plt
delta = 0.025
x = y = np.arange(-3.0, 3.01, delta)
X, Y = np.meshgrid(x, y)
Z1 = plt.mlab.bivariate_normal(X, Y, 1.0, 1.0, 0.0, 0.0)
Z2 = plt.mlab.bivariate_normal(X, Y, 1.5, 0.5, 1, 1)
Z = 10 * (Z1 - Z2)
fig, ((ax1, ax2), (ax3, ax4)) = plt.subplots(2, 2)
fig.set_size_inches(8, 6)
# Your code sample
CS1 = ax1.contourf(X, Y, Z, cmap="viridis", N=256)
ax1.set_title('Your code sample')
ax1.set_xlabel('word length anomaly')
ax1.set_ylabel('sentence length anomaly')
cbar1 = fig.colorbar(CS1, ax=ax1)
# Contour up to N=7 automatically-chosen levels,
# which should give the same as your code.
N = 7
CS2 = ax2.contourf(X, Y, Z, N, cmap="viridis")
ax2.set_title('N=7')
ax2.set_xlabel('word length anomaly')
ax2.set_ylabel('sentence length anomaly')
cbar2 = fig.colorbar(CS2, ax=ax2)
# Contour up to N=100 automatically-chosen levels.
# The resolution is still not as high as using imshow().
N = 100
CS3 = ax3.contourf(X, Y, Z, N, cmap="viridis")
ax3.set_title('N=100')
ax3.set_xlabel('word length anomaly')
ax3.set_ylabel('sentence length anomaly')
cbar3 = fig.colorbar(CS3, ax=ax3)
IM = ax4.imshow(Z, cmap="viridis", origin='lower', extent=(-3, 3, -3, 3))
ax4.set_title("Warren Weckesser's idea")
ax4.set_xlabel('word length anomaly')
ax4.set_ylabel('sentence length anomaly')
cbar4 = fig.colorbar(IM, ax=ax4)
fig.tight_layout()
plt.show()