是否有任何颜色映射,或者是否有一种简单的方法来转换matplotlib颜色映射,从而在0.5附近提供更大的颜色范围,在极端情况下提供更小的颜色范围?我正在创建一组子块,其中一个颜色值是其他颜色值的10倍左右,所以它的值占主导地位,其余的图看起来都一样。举个简单的例子,我们有:
import numpy as np
import matplotlib.pyplot as plt
x = np.linspace(1,10,10)
y = np.linspace(1,10,10)
t1 = np.random.normal(2,0.3,10)
t2 = np.random.normal(9,0.01,10)
t2_max = max(t2)
plt.figure(figsize=(22.0, 15.50))
p = plt.subplot(1,2,1)
colors = plt.cm.Accent(t1/t2_max)
p.scatter(x, y, edgecolors=colors, s=15, linewidths=4)
p = plt.subplot(1,2,2)
colors = plt.cm.Accent(t2/t2_max)
p.scatter(x, y, edgecolors=colors, s=15, linewidths=4)
plt.subplots_adjust(left=0.2)
cbar_ax = plt.axes([0.10, 0.15, 0.05, 0.7])
sm = plt.cm.ScalarMappable(cmap=plt.cm.Accent, norm=plt.Normalize(vmin=0, vmax=t2_max))
sm._A = []
cbar = plt.colorbar(sm,cax=cbar_ax)
plt.show()
t1的变异比t2的变异大得多,但由于t2值高,这种变异不明显。我想要的是一张地图,它将在t1的平均值附近提供更大的颜色梯度,而不需要转换数据本身。我在这里找到了一个解决方案http://protracted-matter.blogspot.co.nz/2012/08/nonlinear-colormap-in-matplotlib.html,但无法使它适用于我的散点图。
编辑: 从下面的答案可以修改为带负数的类,并固定边界。
import numpy as np
import matplotlib.pyplot as plt
x = y = np.linspace(1, 10, 10)
t1mean, t2mean = -6, 9
sigma1, sigma2 = .3, .01
t1 = np.random.normal(t1mean, sigma1, 10)
t2 = np.random.normal(t2mean, sigma2, 10)
class nlcmap(object):
def __init__(self, cmap, levels):
self.cmap = cmap
self.N = cmap.N
self.monochrome = self.cmap.monochrome
self.levels = np.asarray(levels, dtype='float64')
self._x = self.levels
self.levmax = self.levels.max()
self.levmin = self.levels.min()
self.transformed_levels = np.linspace(self.levmin, self.levmax,
len(self.levels))
def __call__(self, xi, alpha=1.0, **kw):
yi = np.interp(xi, self._x, self.transformed_levels)
return self.cmap(yi / (self.levmax-self.levmin)+0.5, alpha)
tmax = 10
tmin = -10
#the choice of the levels depends on the data:
levels = np.concatenate((
[tmin, tmax],
np.linspace(t1mean - 2 * sigma1, t1mean + 2 * sigma1, 5),
np.linspace(t2mean - 2 * sigma2, t2mean + 2 * sigma2, 5),
))
levels = levels[levels <= tmax]
levels.sort()
print levels
cmap_nonlin = nlcmap(plt.cm.jet, levels)
fig, (ax1, ax2) = plt.subplots(1, 2)
ax1.scatter(x, y, edgecolors=cmap_nonlin(t1), s=15, linewidths=4)
ax2.scatter(x, y, edgecolors=cmap_nonlin(t2), s=15, linewidths=4)
fig.subplots_adjust(left=.25)
cbar_ax = fig.add_axes([0.10, 0.15, 0.05, 0.7])
#for the colorbar we map the original colormap, not the nonlinear one:
sm = plt.cm.ScalarMappable(cmap=plt.cm.jet,
norm=plt.Normalize(vmin=tmin, vmax=tmax))
sm._A = []
cbar = fig.colorbar(sm, cax=cbar_ax)
#here we are relabel the linear colorbar ticks to match the nonlinear ticks
cbar.set_ticks(cmap_nonlin.transformed_levels)
cbar.set_ticklabels(["%.2f" % lev for lev in levels])
plt.show()
你的link为colormap提供了很好的解决方案。我编辑了一点,但里面包含了所有必要的东西。你需要为你的非线性颜色映射选择一些合理的级别。我使用了两个以平均值为中心的范围,介于样本的标准偏差之间。通过将其更改为另一个数字,可以获得两个平均值周围颜色的不同局部渐变。
对于颜色条,你
第二种方法可以在查看数据时获得更高的分辨率,看起来更好,实现方式如下:
在结果中,请注意,颜色栏的刻度不是等距的:
您可以使用LinearSegmentedColormap:
这样,您需要在字典中设置颜色查找表,例如下面的“cdict”。
这显示了值之间的转换。我将红色设置为在
t1/t2_max
(0.15到0.35)的值周围变化很大,而蓝色设置为在t2/t2_max
(0.9到1.0)的值周围变化很大。格林什么也不做。我建议阅读docs来了解它是如何工作的。(请注意,这可以自动在您的值周围自动变化)。然后我调整了代码以显示图表:相关问题 更多 >
编程相关推荐