矩阵中的颜色编码值

Question

我有一个 500x500 的单元矩阵，每次迭代都会更新。我正在模拟随机过程，结果新值出现或消失。我想要做的是给值-1 和0 一个固定的颜色代码（比如说绿色和蓝色）。所有其他值（从1 到矩阵中的最大值）可以是与包含-1 或0 的单元格不同的颜色。可以为所有大于0 的值插入颜色。我知道caxis 函数，但这让我只排除-1 和0 的值或从-1 开始进行插值。这个问题有什么解决办法吗？它也需要是一个快速的解决方案，因为每次迭代都会打印矩阵......

[解决方案]

tic
a = randi([-1,10],100,100);
cint = [-1,0,linspace(1,10,10)];
cmap = [0,0,1;0,1,0;autumn(10)];
[~,c] = histc(a,cint);
d = cmap(reshape(c,10000,1),:);
for k=1:3
    im(:,:,k) = reshape(d(:,k),100,100);
end
image(im)    
toc

Answer 1

好吧，既然你特别要求一个单元格。

a = randi([-1,10],10,10); % Generate a 10x10 matrix of random integers
a = num2cell(a);
b = [-1,0,linspace(1,10,5)]; % color interval
cmap = [0,1,0;0,0,1;hot(5)]; % colormap with 7 colors
[~,c] = cellfun(@(x) histc(x,b), a, 'un', false); % Find color index
imdataCells = cellfun(@(x) cmap(x,:), c,'un', false); % Get colors

我还没有测试创建数据的图像，但据我所知（通过查看数据并将 a 与区间进行比较）它应该可以工作。

编辑

要获得正确格式的图像数据，可以使用cellfun 转换图像数据

foo = zeros(size(a,1),size(a,2),3);
for k = 1:3
foo(:,:,k) = cellfun(@(x) x(k),imdataCells);
end

Answer 2

例如surf() 的第四个参数是值颜色的矩阵，您可以更改它或更改使用的颜色图。这种方法 simpy 为您的需求定义了一个颜色图：

% Your Data of -1's, 0's and positive values
cdata = num2cell(peaks(500));
cdata = cell2mat(cdata);
cdata(cdata<0) = -1;
cdata(cdata>0 & cdata <1) = 0;

% the plot
figure()
h = imagesc(cdata);

% the colors
N = 20; % number of colors for interpolation
blue2red = zeros(N,3);
for k = 1:N
    blue2red(k,:) =[k/N, 0, 1-k/N];
end

% adjust the colormap
cmap = zeros(N+2,3);
cmap(1,:) = [0,1,0]; % first is green (for -1)
cmap(2,:) = [0,0,1]; % second is blue (for 0)
cmap(3:end,:) = blue2red; % all the rest is an interpolation from blue to red
colormap(cmap)
colorbar

% Test the picture whilst values are updating
for k = 1:1000;
% Test the picture whilst values are updating
for k = 1:1000;
    old_cdata = get(h, 'CData');        
    cdata = old_cdata + (rand(size(cdata))-0.5)*0.1; % assume values get updated
    cdata(cdata<0) = 0; % assume positive values are at least 0
    cdata(old_cdata == -1) = -1; % assume -1 stays -1
    cdata(old_cdata == 0) = 0; % assume 0 stays 0
    set(h, 'CData', cdata);
    drawnow
end
end