Matlab：contourf或colormap根据半径绘制填充的椭圆

Question

我正在尝试使用contourf或colormap函数根据其arcsin（b / a）值（a =长轴，b =短轴）绘制带有颜色的填充椭圆。

clearvars -except data colheaders
close all
clc
data(:,9)=data(:,9)*pi/180; % Convers Column 9 (angle of rotation) in rad
data(:,6)=1196-data(:,6); % Reset the Y coordinate axis to bottom left

theta = 0 : 0.01 : 2*pi; % Converts phi in rad 

imax=29;
% Define colors
cvalues=asin(data(1:imax,8)./data(1:imax,7))./asin(1);
cm = colormap; % returns the current color map

% Sort and get their index to access the color array 
[~,idx] = sort(cvalues);

% Create colormap
%ColorMap=jet;

for i=1:imax

x = data(i,7)/2 * cos(theta) * cos(data(i,9)) - data(i,8)/2 * sin(theta) * sin(data(i,9)) + data(i,5);
y = data(i,8)/2 * sin(theta) * cos(data(i,9)) + data(i,7)/2 * cos(theta) * sin(data(i,9)) + data(i,6);

colorID = max(1, sum(cvalues(i) > [0:1/length(cm(:,1)):1])); 
ColorMap(i,:) = cm(colorID, :); % returns your color
hold on
% Columns (5,6) are the centre (x,y) of the ellipse
% Columns (7,8) are the major and minor axes (a,b)
% Column 9 is the rotation angle with the x axis

%% TRYING A FASTER WAY OF PLOTTING
%A(:,i)=x'; 
%B(:,i)=y';
%%
fill(x,y,ColorMap(i,:),'EdgeColor', 'None')
text(data(i,5),data(i,6),[num2str(asin(1)*180*cvalues(i)/pi)]) % Assigns number to each ellipse
end

%%
%fill(A,B,ColorMap(1:200,3)','EdgeColor', 'None')
%%

% Adds colorbar to plot
colorbar('SouthOutside')
caxis([0 90])
axis equal;
%xlim([0 7649]);
%ylim([0 15927]);
grid on;

我得到的照片看起来像我认为效果很好：

我没有在椭圆中添加数字，而是添加了我获得的角度（圆圈为90，非常长的椭圆为0）。现在是我真正的实验，我将需要绘制几千个椭圆，我们发现绘制它们需要很多时间，你会看到我们尝试了另一种方法来基本记录数据并一次性绘制。但是如果你有任何建议我们到目前为止还没有成功：）

Answer 1

以下是在Matlab中使用内置色彩映射的方法（查看here以获取完整列表）。

诀窍是创建一个n x 3数组，其中n是您希望表示的颜色数。这是省略号的数量。

您可以像这样创建一个色彩映射：

MyColorMap = jet(n)); %// jet or anything listed in the link above. You can also create your own colormap.

这就是我们将在以下代码中做的事情。为了使订单正确，我们需要对asin(b/a)的值进行排序并获取右椭圆的每个索引。代码被注释，因此很容易遵循。我在图中添加了4个椭圆，以更好地查看颜色差异：

clear
clc
close all

%// Define dummy data
data = zeros(8,9);

data(:,5) = [3;5;12;8;2;7;4;6]; % Centre location X
data(:,6) = [1; -5 ;-2; 4;2;-3;9;5]; % Centre location Y
data(:,7) = [6 ;7;8;6;1;4;2;5]; % Major Axis a
data(:,8) = [2;5;4;2;2;5;3;7]; % Minor axis b

data(:,9) = [10;40;45;90;35;70;85;110]; % Angle of rotation phi

data(:,9)=data(:,9)*pi/180; % Converts phi in rads

theta = 0 : 0.01 : 2*pi;


%// Define colors here
cvalues = asin(data(:,8)./data(:,7));

%// Sort and get their index to access the color array
[~,idx] = sort(cvalues);

%// Create colormap with the "jet" colors
ColorMap = jet(numel(cvalues));

这是包含“颜色”的数组。它看起来像这样：

ColorMap =
                         0                         0                         1
                         0                       0.5                         1
                         0                         1                         1
                       0.5                         1                       0.5
                         1                         1                         0
                         1                       0.5                         0
                         1                         0                         0
                       0.5                         0                         0

因此每一行代表红色，蓝色和绿色（RGB）的组合，并且有多个行可以绘制椭圆。现在填写省略号：

hold all
for i=1:numel(idx)

    k = idx(i);
    x = data(k,8)/2 * cos(theta) * cos(data(k,9)) - data(k,7)/2 * sin(theta) * sin(data(k,9)) + data(k,5);
    y = data(k,7)/2 * sin(theta) * cos(data(k,9)) + data(k,8)/2 * cos(theta) * sin(data(k,9)) + data(k,6);

    fill(x,y,ColorMap(i,:))

    plot(x, y, 'LineWidth', 1);

    % Label each ellipse with a number
    text(data(i,5),data(i,6),num2str(i),'Color','k','FontSize',12)


end
axis equal;
grid on;

使用jet色彩图输出：蓝色是第一个绘制的椭圆，红色是最后一个。您可以根据需要添加颜色栏（添加colorbar）

使用另一个色彩图，bone色彩图，提供以下内容：

希望有所帮助！