Matlab进阶绘图第35期—双特征渲染三维散点图
2023-12-13 22:35:50
在《Matlab论文插图绘制模板第21期—三维散点图(特征渲染赋色)》中,我分享过特征渲染三维散点图的绘制模板。
然而,有的时候,需要在一张图上绘制两个及以上的特征渲染三维散点图,且每个特征渲染三维散点图使用不同的配色方案。
但在Matlab中,一张图上只支持一种colormap/colorbar,所以想要绘制两个及以上的特征渲染三维散点图,需要大家自行设法解决。
比如我在写博士论文的时候,由于时间比较赶,没有时间研究Matlab的实现方案,所以临时用CloudCompare替代实现,至于效果嘛,虽说一般人也挑不出啥毛病,但不够标准,也不够灵活,没办法批量化出图,感觉还是差点意思。
本文利用freezeColors工具(John Iversen, MathWork, 2023),以及我自己制作的colorar_k2工具,进行特征渲染三维散点图的绘制,先来看一下成品效果:
特别提示:本期内容『数据+代码』已上传资源群中,加群的朋友请自行下载。有需要的朋友可以关注同名公号【阿昆的科研日常】,后台回复关键词【绘图桶】查看加入方式。
1. 数据准备
此部分主要是读取原始数据并初始化绘图参数。
% 读取数据
data = load('data.txt');
% 初始化绘图参数
idx1 = find(data(:,4)==1);
idx2 = find(data(:,4)==0);
% 散点数据1
x1 = data(idx1,1);
y1 = data(idx1,2);
z1 = data(idx1,3);
f1 = data(idx1,3);
% 散点数据2
x2 = data(idx2,1);
y2 = data(idx2,2);
z2 = data(idx2,3);
f2 = data(idx2,3);
2. 颜色定义
作图不配色就好比做菜不放盐,总让人感觉少些味道。
但颜色搭配比较考验个人审美,需要多加尝试。
这里直接使用TheColor配色工具中的SCI权威配色库:
%% 颜色定义
map1 = TheColor('sci',2068);
% map1 = flipud(map1);
map2 = TheColor('sci',2044);
map2 = flipud(map2);
3. 双特征渲染三维散点图绘制
调用‘scatter3’和‘freezColors’命令,绘制初始双特征渲染三维散点图。
ax = gca;
% 特征渲染三维散点图1绘制
scatter3(x1, y1, z1, 5, f1, 'filled')
caxis([min(f1(:)) max(f1(:))]);
colormap(map1)
freezeColors;
hold on
% 特征渲染三维散点图2绘制
scatter3(x2, y2, z2, 5, f2, 'filled')
caxis([min(f2(:)) max(f2(:))]);
colormap(map2)
freezeColors;
% 标题、标签、视角
hTitle = title('DoubleFeatureScatter Plot');
hXLabel = xlabel('x');
hYLabel = ylabel('y');
hZLabel = zlabel('z');
view(-37.5,30)
4. 细节优化
为了插图的美观,利用colorbar_k2工具添加颜色条,并对图形细节等进行美化:
% 添加颜色条
colorbar_k2('right',f1,map1,f2,map2)
% 坐标区调整
axes(ax)
axis equal
set(gca, 'Box', 'off', ... % 边框
'LineWidth', 1, 'GridLineStyle', '-',... % 坐标轴线宽
'XGrid', 'on', 'YGrid', 'on', 'ZGrid', 'on',... % 网格
'TickDir', 'out', 'TickLength', [.01 .01], ... % 刻度
'XColor', [.1 .1 .1], 'YColor', [.1 .1 .1],'ZColor', [.1 .1 .1]) % 坐标轴颜色
% 字体和字号
set(gca, 'FontName', 'Arial', 'FontSize', 11)
set([hXLabel,hYLabel,hZLabel], 'FontName', 'Arial', 'FontSize', 11)
set(hTitle, 'FontSize', 12, 'FontWeight' , 'bold')
% 背景颜色
set(gcf,'Color',[1 1 1])
设置完毕后,以期刊所需分辨率、格式输出图片。
%% 图片输出
figW = figureWidth;
figH = figureHeight;
set(figureHandle,'PaperUnits',figureUnits);
set(figureHandle,'PaperPosition',[0?0?figW figH]);
fileout = 'test';
print(figureHandle,[fileout,'.png'],'-r300','-dpng');
以上。
文章来源:https://blog.csdn.net/qq_26447137/article/details/134831796
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。 如若内容造成侵权/违法违规/事实不符,请联系我的编程经验分享网邮箱:veading@qq.com进行投诉反馈,一经查实,立即删除!
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。 如若内容造成侵权/违法违规/事实不符,请联系我的编程经验分享网邮箱:veading@qq.com进行投诉反馈,一经查实,立即删除!