matlab绘图（5）：分级散点色度图

绘图函数

1. 绘图效果
   

2. 复制粘贴存好plot_color_classified_scatter.m函数
   注意这里包含一个通用matlab绘图保存函数plot_save.m，详见帖子：matlab绘图（6）:模块化图像保存，https://blog.csdn.net/CatfishYalu/article/details/141387516?spm=1001.2014.3001.5502

function plot_color_classified_scatter(A, B, C,plot_info,class_info)
    % 函数 plot_Classified_Scatter
    % 绘制一个以A为横轴，B为纵轴的散点图，颜色按照C的值变化，并允许自定义坐标轴和色度条的标签
    % A, B, C: 都应该是相同大小的矩阵
    % class_info: 包含范围和分级信息的向量 [range_min, range_max, range_d]
    % plot_info: 包含绘图标签和保存信息的结构体
    % 例如：
    % class_num: 色块个数
    % decimal: 控制色条刻度标签的小数位数
    % xLabel, yLabel, cLabel: 分别为x轴、y轴和色度条的标签
    % filePath: (可选) 保存图像的文件夹路径
    % imageName: (可选) 保存图像的文件名

    % 将A, B, C矩阵转换为列向量
    x = A(:);
    y = B(:);
    z = C(:); % C的值将用于色彩映射
    
    if nargin<5 % 如果没有class_info矩阵，即手动调节
        range_max = max(z, [], 'omitnan');
        range_min = min(z, [], 'omitnan');
        class_num=plot_info.class_num;
    else
        [range_min,range_max,range_d]=deal(class_info(1),class_info(2),class_info(3));
        class_num=round((range_max-range_min)/range_d);
    end 

    % 定义颜色梯度
    colors = jet(class_num); % 使用10个颜色的'jet'色图，可以根据需要更改颜色图
    edges = linspace(range_min, range_max, class_num+1); % 定义范围0到1，分成10段
      
    % 生成散点图
    [width, height] = deal(plot_info.size(1),plot_info.size(2));
    figure('Units','centimeter','Position',[10 5 width height]);  
    hold on;
    
    % 根据c的值设置每个点的颜色
    for i = 1:length(edges)-1
        % 找到z值在当前范围内的索引
        idx = z >= edges(i) & z < edges(i + 1);
        
        % 绘制相应颜色的散点
        scatter(x(idx), y(idx), 6, colors(i, :), 'filled'); % 6为点的大小，可调整
    end
    
    % 在图形上显示自定义色条
    colormap(colors);
    cb = colorbar('Ticks', edges);
    tickFormat = sprintf('%%.%df', plot_info.decimal); % 根据输入的小数位数设置格式
    cb.TickLabels = arrayfun(@(x) num2str(x, tickFormat), edges, 'UniformOutput', false);
    caxis([min(edges), max(edges)]);
    
    % 设置色条的刻度及对应数字
    ylabel(cb, plot_info.cLabel); % 设置色度条的标题

    % 设置轴标签字体和字号
    set(gca, 'FontName', 'Times New Roman', 'FontSize', 8);
    %% 绘制辅助信息及图像保存
    plot_save(plot_info)
end

3. 函数调用

plot_info=table;
plot_info.size=[8,6];
plot_info.xLabel = 'Fc';
plot_info.yLabel = 'LST(K)';
plot_info.cLabel = 'albedo';
plot_info.filePath = fullfile(wdir,'\results\6_scatter_AB_C\');
plot_info.imageName = strcat(cal_date,'_AB_C'); 
plot_info.title=strcat(cal_date,'-AB_C');
plot_info.save=1;
plot_info.decimal=2;%控制色条刻度标签的小数位数
%% 若自动分级用这几句代码
% plot_info.class_num=10;
% plot_color_classified_scatter(Fc, LST, albedo, plot_info);
%% 若手动分级用这几句代码
class_info=[0,0.3,0.05];
plot_color_classified_scatter(Fc, LST, albedo,plot_info,class_info);  

附图1和2分别是手动分级和自动分级的效果，都还可以，相对来说自动分级更方便些，手动方便自己调整，按需取用即可