Matlab进阶绘图第41期—双三角网格曲面图

在《Matlab论文插图绘制模板第67期---三角网格图(Trimesh)》中,我分享过三角 网格 曲面图的绘制模板。

然而,有的时候,需要在一张图上绘制两个及以上的三角网格曲面图,且每个三角网格曲面图使用不同的配色方案

在Matlab中,一张图上只支持一种colormap/colorbar,所以想要绘制两个及以上的三角网格曲面图,需要大家自行设法解决

本文利用freezeColors工具 (John Iversen, MathWork, 2023),以及我自己制作的colorbar_k2工具,进行双三角网格曲面的绘制,先来看一下成品效果:

特别提示:本期内容『数据+代码』已上传资源群中,加群的朋友请自行下载。有需要的朋友可以关注同名公号【阿昆的科研日常】,后台回复关键词【绘图桶】查看加入方式


1. 数据准备

此部分主要是读取原始数据并初始化绘图参数

复制代码
% 读取数据load data.mat% 初始化绘图参数inter = 8;% 网格1抽稀x1 = X(1:inter:end,1:inter:end);y1 = Y(1:inter:end,1:inter:end);z1 = Z1(1:inter:end,1:inter:end);% 网格2抽稀x2 = X(1:inter:end,1:inter:end);y2 = Y(1:inter:end,1:inter:end);z2 = Z2(1:inter:end,1:inter:end);

2. 颜色定义

作图不配色就好比做菜不放盐,总让人感觉少些味道。

但颜色搭配比较考验个人审美,需要多加尝试。

这里直接使用TheColor配色工具中的SCI权威配色库

复制代码
%% 颜色定义map1 = TheColor('sci',2068);% map1 = flipud(map1);map2 = TheColor('sci',2073);

3. 双三角网格曲面图绘制

调用'trimesh'和'freezColors'命令,绘制初始双三角网格曲面图

复制代码
ax = gca;% 三角网格曲面1绘制T1 = delaunay(x1,y1);% 三角剖分trimesh(T1,x1,y1,z1,'linewidth',0.8)caxis([min(z1(:)) max(z1(:))]);colormap(map1)freezeColors; hold on% 三角网格曲面2绘制T2 = delaunay(x2,y2);% 三角剖分trimesh(T2,x2,y2,z2,'linewidth',0.8)caxis([min(z2(:)) max(z2(:))]);colormap(map2)freezeColors; % 标题、标签、视角hTitle = title('DoubleTrimesh Plot');hXLabel = xlabel('x');hYLabel = ylabel('y');hZLabel = zlabel('z');view(-35,30)

4. 细节优化

为了插图的信息完整性,利用colorbar_k2工具添加颜色条,并 对图形细节等进行美化

复制代码
% 添加颜色条colorbar_k2('right',Z1,map1,Z2,map2)% 坐标区调整axes(ax)axis tightset(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],...       % 坐标轴颜色         'zlim',[0 700])% 字体和字号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');

以上。

相关推荐
郝学胜-神的一滴37 分钟前
Effective Python 第44条:用纯属性与修饰器取代旧式的 setter 与 getter 方法
开发语言·python·程序人生·软件工程
程序员莫小特2 小时前
老题新解|计算2的N次方
开发语言·数据结构·算法·青少年编程·信息学奥赛一本通
white-persist3 小时前
XXE 注入漏洞全解析:从原理到实战
开发语言·前端·网络·安全·web安全·网络安全·信息可视化
wearegogog1233 小时前
基于块匹配的MATLAB视频去抖动算法
算法·matlab·音视频
人生导师yxc4 小时前
Java中Mock的写法
java·开发语言
半路程序员4 小时前
Go语言学习(四)
开发语言·学习·golang
沐知全栈开发4 小时前
C# 枚举(Enum)
开发语言
秦禹辰4 小时前
轻量级开源文件共享系统PicoShare本地部署并实现公网环境文件共享
开发语言·后端·golang
脑子慢且灵4 小时前
C语言与Java语言编译过程及文件类型
java·c语言·开发语言·汇编·编辑器
蒙奇D索大4 小时前
【C语言加油站】C语言文件操作详解:从“流”的概念到文件的打开与关闭
c语言·开发语言·笔记·学习·改行学it