MATLAB使用多个扇形颜色变化表示空间一个点的多种数值

MATLAB使用多个扇形颜色变化表示空间一个点的多种数值

excel中表格中数据格式，多行

lon	lat	data1	data2	data3
117	38	0.1	1	100

clear;close all;
figure('Position',[100 100 800 800]);

num_points = 14; 
[num,txt,raw] = xlsread('test.xlsx');
x = num(:,1);
y = num(:,2);


data1 = num(:,3);  % 维度1 
data2 = num(:,4);   % 维度2 
data3 = num(:,5);  % 维度3
% 归一化数据到 [0,1]
data1_norm = (data1 - min(data1)) / (max(data1) - min(data1));
data2_norm = (data2 - min(data2)) / (max(data2) - min(data2));
data3_norm = (data3 - min(data3)) / (max(data3) - min(data3));


cmap_red = hot(256);   
cmap_green = cool(256);
cmap_blue = copper(256); 


% 扇形的半径和角度
radius = 0.06;  % 扇形半径（根据点间距调整）
theta = linspace(0, 2*pi, 100);  % 完整圆周
sector_angle = 120;  % 每个扇形占120度

for i = 1:num_points
    % 当前点的中心坐标
    center_x = x(i);
    center_y = y(i);
    
    % 绘制三个扇形（每个扇形120度）
    for sector = 1:3
        % 计算扇形起始和结束角度
        start_angle = (sector-1) * sector_angle;
        end_angle = sector * sector_angle;
        
        % 生成扇形边缘坐标（极坐标转笛卡尔坐标）
        angles = linspace(start_angle, end_angle, 50);
        angles_rad = deg2rad(angles);
        x_sector = center_x + radius * [0, cos(angles_rad), 0];
        y_sector = center_y + radius * [0, sin(angles_rad), 0];
        
        % 根据数据维度选择颜色
        switch sector
            case 1  % 第一个扇形对应data1
                color_value = cmap_red(ceil(data1_norm(i)*255 + 1), :);
            case 2  % 第二个扇形对应data2
                color_value = cmap_green(ceil(data2_norm(i)*255 + 1), :);
                
            case 3  % 第三个扇形对应data3
                color_value = cmap_blue(ceil(data3_norm(i)*255 + 1), :);
               
        end
        
        % 绘制扇形
        patch(x_sector, y_sector, color_value, 'EdgeColor', 'none','Markersize',12,'FaceAlpha',1);
    end
end

hold on;
plot(x,y,'k.','MarkerSize',10)
axis equal;
xlim([117, 121]);
ylim([34, 37]);
box on;
set(gca,'FontSize',16);
set(gca,'Layer','top','linewidth',1.2,'xminortick','on','yminortick','on');
mt = 0.4;

ax_red = axes('Position', [0.1+mt 0.25 0.2 0.2]);
colormap(ax_red, cmap_red);
caxis(ax_red, [min(data1), max(data1)]);
c = colorbar(ax_red, 'Location', 'east');
set(get(c,'title'),'string','气温','fontsize',10);
ax_red.Visible = 'off';



ax_green = axes('Position', [0.2+mt 0.25 0.2 0.2]);
colormap(ax_green, cmap_green);
caxis(ax_green, [min(data2), max(data2)]);
c = colorbar(ax_green, 'Location', 'east');
ax_green.Visible = 'off';
set(get(c,'title'),'string','气压','fontsize',10);


ax_blue = axes('Position', [0.3+mt 0.25 0.2 0.2]);
colormap(ax_blue, cmap_blue);
caxis(ax_blue, [min(data3), max(data3)]);
c = colorbar(ax_blue, 'Location', 'east');
ax_blue.Visible = 'off';
set(get(c,'title'),'string','降水','fontsize',10);