绘图函数
- 绘图效果
- 复制粘贴存好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
- 函数调用
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分别是手动分级和自动分级的效果,都还可以,相对来说自动分级更方便些,手动方便自己调整,按需取用即可