MATLAB | 绘图复刻 (二十七) | 多子图堆叠面积图

今天复刻的是：

Barton, A., Afrad, M.H., Taylor-Brown, A. et al. Evolution of pandemic cholera at its global source. Nature 653, 491--498 (2026). https://doi.org/10.1038/s41586-026-10340-x

的 Fig.2a，原图长这样：

复刻效果(这里用的是我随机生成的数据)

正文部分

数据读取

matlab 复制代码

% 数据读取
T1 = readtable('Group_AAAA.csv');
T2 = readtable('Group_BBBB.csv');
T3 = readtable('Group_CCCC.csv');
Data1 = T1{:, 4:end};
Data2 = T2{:, 4:end};
Data3 = T3{:, 4:end};

% 获取变量名
varNames = T1.Properties.VariableNames(4:end);
varNames = strrep(varNames, '_', '.');
gNames = {'Group-AAAA', 'Group-BBBB', 'Group-CCCC'};

% 获取日期
Year = T1.Year; Month = T1.Month; Day = T1.Day;
DT = datetime(Year, Month, Day);

这三组数据是我随机生成的，三组数据都长这样：

前三列分别是日期的年月日，后面的列都是用来绘图的变量。

配色获取

这里使用我编写的常用colormap合集工具slanCM:

https://www.mathworks.com/matlabcentral/fileexchange/120088-200-colormaps

matlab 复制代码

% 获取配色
cmap = [flipud(slanCM(17, 18)); flipud(slanCM(19, 18))];

使用其中的 17 号配色 RdPu 及 19 号配色 GnBu：

绘图及相关技巧

由于 area 函数支持使用时间直接作为 x 轴，所以绘图非常简单，就创建三个子图然后依次：

matlab 复制代码

aHdl1 = area(ax1, DT, Data1);
aHdl2 = area(ax2, DT, Data2);
aHdl3 = area(ax3, DT, Data3);

这里讲一个修改颜色的技巧，可以将颜色使用 num2cell 函数分割成元胞数组列向量，然后使用 set 函数批量设置属性：

复制代码

set(aHdl1, {'FaceColor'}, num2cell(cmap, 2));
set(aHdl2, {'FaceColor'}, num2cell(cmap, 2));
set(aHdl3, {'FaceColor'}, num2cell(cmap, 2));

这样的写法对比使用 for 循环一个一个设置颜色快十几倍，这里只是设置 36 个对象可能看起来没那么明显，但是数量再多的话，比如成百上千个对象会快非常多。

数值堆叠面积图-绘图部分完整代码

matlab 复制代码

% 图窗及坐标区域创建
fig = figure('Units','normalized', 'Position',[.02,.05,.9,.88], 'Color','w');
ax1 = axes('Parent',fig, 'Position',[.1, 2/3 + .065, .7, 1/3 - .1], 'NextPlot','add', ...
    'TickLength',[.005,.001], 'LineWidth',1.5, 'TickDir','out', 'FontSize',13, 'FontName','Arial');
ax2 = axes('Parent',fig, 'Position',[.1, 1/3 + .065, .7, 1/3 - .1], 'NextPlot','add', ...
    'TickLength',[.005,.001], 'LineWidth',1.5, 'TickDir','out', 'FontSize',13, 'FontName','Arial');
ax3 = axes('Parent',fig, 'Position',[.1, 0   + .065, .7, 1/3 - .1], 'NextPlot','add', ...
    'TickLength',[.005,.001], 'LineWidth',1.5, 'TickDir','out', 'FontSize',13, 'FontName','Arial');
% 绘图
aHdl1 = area(ax1, DT, Data1, 'LineWidth',1.5); set(aHdl1, {'FaceColor'}, num2cell(cmap, 2));
aHdl2 = area(ax2, DT, Data2, 'LineWidth',1.5); set(aHdl2, {'FaceColor'}, num2cell(cmap, 2));
aHdl3 = area(ax3, DT, Data3, 'LineWidth',1.5); set(aHdl3, {'FaceColor'}, num2cell(cmap, 2));

% 添加标签及图例
set(ax1.YLabel, 'String',gNames{1}, 'FontSize',16)
set(ax2.YLabel, 'String',gNames{2}, 'FontSize',16)
set(ax3.YLabel, 'String',gNames{3}, 'FontSize',16)
set(ax3.XLabel, 'String','Year', 'FontSize',16)
ax1.XTickLabel = regexprep(ax1.XTickLabel, '\D', '');
ax2.XTickLabel = regexprep(ax2.XTickLabel, '\D', '');
ax3.XTickLabel = regexprep(ax3.XTickLabel, '\D', '');
annotation('textbox', [.03, .4995, .001, .001], 'String','Value', 'FontSize',21, 'FontName','Arial', ...
    'HorizontalAlignment','center', 'VerticalAlignment','middle', 'Rotation',90, 'EdgeColor','none');
lgdHdl = legend(ax1, varNames, 'AutoUpdate','off', 'Position', [.825, .05, .15, .9], 'NumColumns',2, 'Box','off');
lgdHdl.ItemTokenSize = [22, 22];
set(lgdHdl.Title, 'String','Sublineage', 'FontSize',21, 'FontName','Arial')

比例堆叠面积图-绘图部分完整代码

matlab 复制代码

% 图窗及坐标区域创建
fig = figure('Units','normalized', 'Position',[.02,.05,.9,.88], 'Color','w');
ax1 = axes('Parent',fig, 'Position',[.1, 2/3 + .065, .7, 1/3 - .1], 'NextPlot','add', 'YLim',[0,1], ...
    'TickLength',[.005,.001], 'LineWidth',1.5, 'TickDir','out', 'FontSize',13, 'FontName','Arial');
ax2 = axes('Parent',fig, 'Position',[.1, 1/3 + .065, .7, 1/3 - .1], 'NextPlot','add', 'YLim',[0,1], ...
    'TickLength',[.005,.001], 'LineWidth',1.5, 'TickDir','out', 'FontSize',13, 'FontName','Arial');
ax3 = axes('Parent',fig, 'Position',[.1, 0   + .065, .7, 1/3 - .1], 'NextPlot','add', 'YLim',[0,1], ...
    'TickLength',[.005,.001], 'LineWidth',1.5, 'TickDir','out', 'FontSize',13, 'FontName','Arial');
% 绘图
N = length(varNames);
aHdl1 = area(ax1, DT, Data1./repmat(sum(Data1, 2), [1, N]), 'LineWidth',1.5); set(aHdl1, {'FaceColor'}, num2cell(cmap, 2));
aHdl2 = area(ax2, DT, Data2./repmat(sum(Data2, 2), [1, N]), 'LineWidth',1.5); set(aHdl2, {'FaceColor'}, num2cell(cmap, 2));
aHdl3 = area(ax3, DT, Data3./repmat(sum(Data3, 2), [1, N]), 'LineWidth',1.5); set(aHdl3, {'FaceColor'}, num2cell(cmap, 2));

% 添加标签及图例
set(ax1.YLabel, 'String',gNames{1}, 'FontSize',16)
set(ax2.YLabel, 'String',gNames{2}, 'FontSize',16)
set(ax3.YLabel, 'String',gNames{3}, 'FontSize',16)
set(ax3.XLabel, 'String','Year', 'FontSize',16)
ax1.XTickLabel = regexprep(ax1.XTickLabel, '\D', '');
ax2.XTickLabel = regexprep(ax2.XTickLabel, '\D', '');
ax3.XTickLabel = regexprep(ax3.XTickLabel, '\D', '');
annotation('textbox', [.03, .4995, .001, .001], 'String','Proportion', 'FontSize',21, 'FontName','Arial', ...
    'HorizontalAlignment','center', 'VerticalAlignment','middle', 'Rotation',90, 'EdgeColor','none');
lgdHdl = legend(ax1, varNames, 'AutoUpdate','off', 'Position', [.825, .05, .15, .9], 'NumColumns',2, 'Box','off');
lgdHdl.ItemTokenSize = [22, 22];
set(lgdHdl.Title, 'String','Sublineage', 'FontSize',21, 'FontName','Arial')

% 尝试绘制边框
try
    xline(ax1, ax1.XLim, 'Color','k', 'Alpha',1, 'LineWidth',1.5)
    xline(ax2, ax2.XLim, 'Color','k', 'Alpha',1, 'LineWidth',1.5)
    xline(ax3, ax3.XLim, 'Color','k', 'Alpha',1, 'LineWidth',1.5)
catch
end

结

完整代码、相关工具函数、及本文所使用数据已经放在如下gitee仓库：

https://gitee.com/slandarer/PLTreprint/

网盘链接，不做后续更新：

https://pan.baidu.com/s/1ZW8jeE87n8dTGilDHlANHA?pwd=slan