Matlab进阶绘图第1期—山脊图/嵴线图/峰峦图
2023/2/24 15:28:31 阅读:199 发布者:
山脊图又称峰峦图、嵴线图,英文常称为Ridgeline Plot、Joy Plot。
山脊图一般由共享X轴的多组彼此重叠的填充折线图组成,其常用于对不同组数据间的分布情况进行比较,在一些顶级期刊中也能看见它的身影,比如《Science》:
《Nature》:
对于Matlab用户,由于Matlab中未收录山脊图的绘制函数,因此需要大家自行想办法(例如使用fill函数)。
本文使用Santiago Benito制作的joyPlot工具(Matlab Central, 2020)进行山脊图的绘制,先来看一下成品效果:
1. 数据准备
此部分主要是读取原始数据,定义自变量和因变量。
% 读取数据
load data.mat
% 自变量
X = x;
% 因变量
Y = data;
2. 颜色定义
作图不配色就好比做菜不放盐,总让人感觉少些味道。
但颜色搭配比较考验个人审美,需要多加尝试。
这里直接使用TheColor配色工具中的SCI权威配色库中的渐变色:
%% 颜色定义
map = TheColor('sci',2064);
map = flipud(map);
(点击上图查看TheColor具体功能)
获取方式:公众号(阿昆的科研日常)后台回复 TC
3. 山脊图绘制
调用‘joyPlot’函数,绘制初始山脊图。
joyPlot(data',x,0.7, ...
'constant', ...
false,...
'LineWidth',1.2, ...
'FaceColor','position', ...
'FaceAlpha',0.8)
hTitle = title('Ridgeline Plot');
hXLabel = xlabel('K (w)');
hYLabel = ylabel('Samples');
4. 细节优化
为了插图的美观,将初始山脊图赋上之前选择的颜色:
% 赋色
colormap(map)
然后,对坐标轴细节等进行美化,设置完毕后,以期刊所需分辨率、格式输出图片。
%%
输出
figW = figureWidth;
figH = figureHeight;
set(figureHandle,'PaperUnits',figureUnits);
set(figureHandle,'PaperPosition',[0 0 figW figH]);
fileout = 'test';
print(figureHandle,[fileout,'.png'],'-r300','-dpng');
转自:“阿昆的科研日常”微信公众号
如有侵权,请联系本站删除!