山脊图又称峰峦图、嵴线图,英文常称为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);
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');
也可以尝试其它配色: