Reprint of illustrations in nature, issue 3 - area map (part2-100)

In the previous post , Reprinted 《Nature》 HD paper insertions Part2 No 456 Multiple groups of histograms at the bottom left of illustration No ：

Part1 pass the civil examinations 305 The horizontal single group multicolor histogram on the left side of illustration No ：

This issue , Let's go back and forth 《Nature》 HD paper insertions Part2 pass the civil examinations 100 Illustration No （ following abbreviation 100） Area map at the top right .

Let's take a look first Replica effect ：

hot tip ：Nature Thesis illustration reprint series , It aims to reduce the threshold of scientific research mapping , Let everyone no longer feel “ The illustrations of papers in top journals are out of reach for me ”. If it works , Share it with the people around you .

1 Factor analysis

100 in Every curve object utilize Different colors with certain transparency are filled , These colors are on the top left (a) Of Legend Can be found in .

Besides , Its coordinate axis is left 、 On 、 The right part has no scale , And the scale label of the axis below it is vertical .

These are the points that need attention .

2 Replica tool

• Matlab R2022a

• Matlab Scientific research drawing plug-in -ColorCopy Color scheme extraction tool

• Matlab Scientific research drawing plug-in -ColorMap Gradient generation / Editing tools

• Matlab Scientific research drawing plug-in -Fig2Data Image raw data extraction tool

3 Data preparation and color scheme extraction

This part is mainly extract 100 Data in and color scheme , For drawing .

First , Use Fig2Data The tool extracts curve data in the graph （ For convenience , Only 20 One of the curves 16 For subsequent processing ）.

Due to the use Fig2Data Tool manually extracted curve Insufficient sampling density , There will be some ‘ angular ’：

So use ‘interp1’ command , use ‘spline’ Interpolation method The manual sampling points are encryption ：

The encrypted data is saved in datax.mat and datay.mat In file .

%%  Data preparation %  Reading data %  The independent variables dataxload datax.mat;%  The dependent variable datayload datay.mat;

then , utilize ColorCopy Tools extract color schemes ：

%%  Color extraction C = ColorCopy;close

Because the original curve is 20 strip , And I only extracted one of them 16 strip , therefore Need to reassign colors .

However , The original author used Gradient color matching , If arbitrarily assigned , May cause uneven color changes .

therefore , Here I am Used ColorMap Gradient generation / Editor pair ColorCopy The resulting color matrix C Downsampling ：

C = ColorMap(C,N);

among ,N Represents the number of colors after downsampling , namely 16.

4 mapping

First , utilize ‘area’ command Draw system Initial area map .

for i = 1:N    area(datax(i,:),datay(i,:),'LineWidth',2,'FaceColor',C(i,:),'EdgeColor',C(i,:),...          'FaceAlpha',.3,'EdgeAlpha',1);    hold onendhXLabel = xlabel('Cortical depth (normalized)');hYLabel = ylabel('cell density');hTitle = title('MOp cells');

To avoid axis occlusion , Use ‘plot’ Command to draw another coordinate axis .

hold onplot(0:0.2:1.2,zeros(1,7),'k','LineWidth', 1)

And then , Compare with the original drawing , Further adjust the details of the coordinate area ：

%  Coordinate axis parameter adjustment set(gca, 'LineWidth', 1,...                                        %  Line width          'Box', 'off', ...                                         %  Frame          'XGrid', 'off', 'YGrid', 'off', ...                       %  grid          'TickDir', 'out', 'TickLength', [.01 .01], ...            %  scale          'XMinorTick', 'off', 'YMinorTick', 'off', ...             %  Small scale          'XColor', [.1 .1 .1],  'YColor', [.1 .1 .1],...           %  Axis color          'XTick', 0:0.2:1.2,...                                    %  Scale and range          'XLim', [-0.05 1.1],...         'XTickLabels',{num2str([0:0.2:1]','%.1f')},...         'XTickLabelRotation',90,...         'YTick', [],...         'YLim', [0 1])%  add 、 Right border hold onXL = get(gca,'xlim'); XR = XL(2);YL = get(gca,'ylim'); YT = YL(2);xc = get(gca,'XColor');yc = get(gca,'YColor');plot(XL,YT*ones(size(XL)),'color', xc,'linewidth',1)plot(XR*ones(size(YL)),YL,'color', yc,'linewidth',1)%  Font size set(gca, 'FontName', 'Arial', 'FontSize', 12)set([hXLabel,hYLabel], 'FontName',  'Arial', 'FontSize', 13)set(hTitle, 'FontSize', 14)%  The background color set(gcf,'Color',[1 1 1])

Last , Output pictures at the desired resolution ：

figW = figureWidth;figH = figureHeight;set(figureHandle,'PaperUnits',figureUnits);set(figureHandle,'PaperPosition',[0 0 figW figH]);fileout = 'test';print(figureHandle,[fileout,'.png'],'-r300','-dpng');

It can be said that the reduction degree is still relatively high .

You also want to see which picture I engrave , Can you tell me the number .

above .

Access method ：

Nature Of area map Matlab Reprint code