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

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

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

lon lat data1 data2 data3
117 38 0.1 1 100
python 复制代码
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);
相关推荐
乌鸦94437 分钟前
《STL--stack 和 queue 的使用及其底层实现》
开发语言·c++·priority_queue·适配器stack、queue
Chef_Chen1 小时前
从0开始学习R语言--Day12--泊松分布
开发语言·学习·r语言
廿二松柏木1 小时前
MATLAB实现井字棋
开发语言·matlab
小智学长 | 嵌入式1 小时前
SOC-ESP32S3部分:22-分区表
开发语言·单片机·esp32
gou123412341 小时前
【Golang进阶】第八章:并发编程基础——从Goroutine调度到Channel通信实战
开发语言·后端·golang
一只自律的鸡1 小时前
STL之vector
开发语言·c++·算法
此刻我在家里喂猪呢2 小时前
qt之开发大恒usb3.0相机三
开发语言·qt
哈哈哈哈哈哈哈哈哈...........2 小时前
【Java】ForkJoin 框架
java·开发语言
zh_xuan2 小时前
c++ typeid运算符
开发语言·c++
✎ ﹏梦醒͜ღ҉繁华落℘3 小时前
WPF学习
c语言·开发语言·笔记