密度散点图本质上是一种特征渲染的散点图，其颜色表示某一点所在区域的密度信息。

除了作图，密度散点图绘制的关键还在于密度的计算。

当然，不管是作图还是密度的计算，这些在《Matlab论文插图绘制模板》和《Matlab点云处理及可视化》往期文章中都有介绍。

为了方便起见，本期就直接分享一下密度散点图的绘制方法，先来看一下成品效果：

1. 数据准备

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

% 读取数据

load data.mat

% 初始化绘图参数

data = [x,y];

% 密度计算

radius = 1.5; % 定义半径

density_2D = density2D_KD(data(:,1:2),radius); % 2D平面密度

2. 颜色定义

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

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

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

%% 颜色定义

map = TheColor('sci',2064);

map = flipud(map);

3. 密度散点图绘制

调用‘scatter’命令，绘制初始密度散点图。

scatter(data(:,1), data(:,2), 5, density_2D, 'filled')

hTitle = title('Satellite-derived bathymetry');

hXLabel = xlabel('ICESat-2 bathymetric points in depth (m)');

hYLabel = ylabel('Estimated depth (m)');

4. 细节优化

为了插图的美观，将初始密度散点图赋上之前选择的颜色：

% 赋色

colormap(map)

colorbar

进一步，对坐标轴细节等进行美化：

% 坐标轴美化

set(gca, 'Box', 'off', ... % 边框

'LineWidth',1,... % 线宽

'XGrid', 'on', 'YGrid', 'on', ... % 网格

'TickDir', 'out', 'TickLength', [.005 .005], ... % 刻度

'XMinorTick', 'off', 'YMinorTick', 'off', ... % 小刻度

'XColor', [.1 .1 .1], 'YColor', [.1 .1 .1],... % 坐标轴颜色

'XTick', 0:40:160,... % 坐标区刻度、范围

'XLim', [0 160],...

'YTick', 0:40:160,...

'YLim', [0 160])

% 字体和字号

set(gca, 'FontName', 'Arial', 'FontSize', 10)

set([hXLabel, hYLabel], 'FontSize', 11, 'FontName', 'Arial')

set(hTitle, 'FontSize', 12, 'FontWeight' , 'bold')

% 背景颜色

set(gcf,'Color',[1 1 1])

% 添加上、右框线

xc = get(gca,'XColor');

yc = get(gca,'YColor');

unit = get(gca,'units');

ax = axes( 'Units', unit,...

'Position',get(gca,'Position'),...

'XAxisLocation','top',...

'YAxisLocation','right',...

'Color','none',...

'XColor',xc,...

'YColor',yc);

set(ax, 'linewidth',1,...

'XTick', [],...

'YTick', []);

设置完毕后，以期刊所需分辨率、格式输出图片。

%% 图片输出

figW = figureWidth;

figH = figureHeight;

set(figureHandle,'PaperUnits',figureUnits);

set(figureHandle,'PaperPosition',[0 0 figW figH]);

fileout = 'test';

print(figureHandle,[fileout,'.png'],'-r300','-dpng');

也可以尝试其它配色：

以上。

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

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