在《Matlab论文插图绘制模板第53期—网格曲面图（Mesh）》中，我分享过网格曲面图的绘制模板。

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

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

本文利用freezeColors工具（John Iversen, MathWork, 2023），以及我自己制作的colorar_k2工具，进行双网格曲面的绘制，先来看一下成品效果：

1. 数据准备

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

% 读取数据

load data.mat

% 预处理/初始化绘图参数

% 网格1抽稀

x1 = X(1:5:end,1:5:end);

y1 = Y(1:5:end,1:5:end);

z1 = Z1(1:5:end,1:5:end);

% 网格2抽稀

x2 = X(1:5:end,1:5:end);

y2 = Y(1:5:end,1:5:end);

z2 = Z2(1:5:end,1:5:end);

2. 颜色定义

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

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

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

%% 颜色定义

map1 = TheColor('sci',2068);

% map1 = flipud(map1);

map2 = TheColor('sci',2073);

3. 双网格曲面图绘制

调用‘mesh’和‘freezColors’命令，绘制初始双网格曲面图。

ax = gca;

% 网格曲面1绘制

s = mesh(x1,y1,z1,'LineWidth',1);

caxis([min(z1(:)) max(z1(:))]);

colormap(map1)

freezeColors;

hold on

% 网格曲面2绘制

s2 = mesh(x2,y2,z2,'LineWidth',1);

caxis([min(z2(:)) max(z2(:))]);

colormap(map2)

freezeColors;

% 标题、标签、视角

hTitle = title('DoubleMesh 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 tight

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],... % 坐标轴颜色

'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');

以上。

转自：“阿昆的科研日常”微信公众号

如有侵权，请联系本站删除！