How to map a vector to a colourmap ?

77 次查看（过去 30 天）
Harish Rajan
Harish Rajan 2020-3-24
评论： Harish Rajan 2020-3-30
Hello,
I am trying to map a vector containing energy of modes of a dynamic system to a vector containing colour values so that when i plot this I want the markers to have the colour corresponding to their energy levels. For example the highest energy value will correspond to a Solid black colour and the next one a little less black then the others have grey and the least energy ones have white colour (thus making them invisible in the plot).
Best Regards
Harish
  2 个评论
Adam
Adam 2020-3-24
Does regular plotting and applying a colourmap not give what you want?
e.g.
A = magic(25);
figure; imagesc( A );
colorbar
You can certainly do mapping onto a colourmap manually and produce a true RGB image if you wish, but it is obviously more effort and data cursor information would reflect that it is now RGB rather than just give you the raw value.
Alternatively you can manipulate your data in some way if a straight linear mapping onto the colourmap is not what you want. Or create a custom colourmap.
Hard to say without seeing an example though.
Harish Rajan
Harish Rajan 2020-3-24
But the point is I dont want to have an RGB image. The objective is to create a vector with colour values corresponding to the energy values, so that i can plot them normally using the plot command and when a issue a fill command for the markers, instead of all them being same colour they follow the one given by the colour map.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

J. Alex Lee
J. Alex Lee 2020-3-24
A brute force way would be to interpolate...not sure if there's a better way
x = % ... data vector
indices = 0:255;
xref = linspace(min(x),max(x),length(indices));
indexed_x = interp1(xref,indices,'nearest');
% then you can map to whatever colormap you want, e.g.
graymap = gray(length(indices));
xcolors = graymap(indexed_x,:);
or something like that?
  5 个评论
显示 3更早的评论
J. Alex Lee
J. Alex Lee 2020-3-24
编辑：J. Alex Lee 2020-3-24
upon looking at docs of scatter(), i see that it is doing more or less exactly what you want, AND it can use the axes underlying colormap, to Adam's point ( you don't have to manually index). Take a look at the docs https://www.mathworks.com/help/matlab/ref/scatter.html#d120e1106626
scatter(x,y,sz,energy)
colormap gray
J. Alex Lee
J. Alex Lee 2020-3-24
Just for completeness, Image Analyst's full example but with using scatter:
theta = linspace(0,(2*pi),30);
x=cos(theta);
y=sin(theta);
energy = linspace(0,1000,30);
figure(1);
scatter(x,y,120,energy,'filled')
grid on
axis square
colormap jet
% to do the black and white example:
% colormap gray

请先登录，再进行评论。

更多回答（1 个）

Image Analyst
Image Analyst 2020-3-24
Do you mean where the marker color is the color taken from the "energy" array, like this:
% Define our data.
theta = linspace(0,(2*pi),30);
x=cos(theta);
y=sin(theta);
energy = linspace(0,1000,30);
numPoints = length(x);
% Get a colormap, a unique color for every energy level
cmap = jet(numPoints); % Initialize jet colormap.
% Get energy in the range 1 to numPoints so we can use that to get a row from the colormap.
qEnergy = imquantize(energy, numPoints);
for k = 1 : numPoints
% Get the color for this energy level:
thisEnergy = qEnergy(k);
thisColor = cmap(thisEnergy);
fprintf('Plotting point #%d at (%.3f, %.3f) with color (%.3f, %.3f, %.3f)\n',...
k, x(k), y(k), cmap(k, 1), cmap(k, 2), cmap(k, 3));
plot(x(k), y(k), '.', 'Color', cmap(k, :), 'MarkerSize', 40);
hold on;
end
grid on;
axis square;
fprintf('Done running %s.m ...\n', mfilename);
  4 个评论
显示 2更早的评论
Image Analyst
Image Analyst 2020-3-30
Harish, you can see how I defined the colormap before the loop. You can make it be anything you want. You can also see how I said the 'MarkerSize' was 40 in the call to plot(). You can change that 40 to be whatever you want. It can even be a function of k, the loop iterator if you want.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Orange 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by