三维密度散点图本质上是一种特征渲染的三维散点图，其颜色表示某一点所在区域的密度信息。
除了作图，三维密度散点图绘制的关键还在于密度的计算。
当然，不管是作图还是密度的计算，这些在《Matlab论文插图绘制模板》和《Matlab点云处理及可视化》往期文章中都有介绍。
为了方便起见，本期就直接分享一下三维密度散点图的绘制方法，先来看一下成品效果：
1. 数据准备
此部分主要是读取原始数据并初始化绘图参数。
% 读取数据
load data.mat
% 初始化绘图参数
data = [x,y,z];
% 密度计算
radius = 1; % 定义半径
density_3D = density3D_KD(data(:,1:3),radius); % 3D密度
2. 颜色定义
作图不配色就好比做菜不放盐，总让人感觉少些味道。
但颜色搭配比较考验个人审美，需要多加尝试。
这里直接使用TheColor配色工具中的SCI权威配色库：
%% 颜色定义
map = TheColor('sci',2097);
% map = flipud(map);
3. 密度散点图绘制
调用‘scatter3’命令，绘制初始三维密度散点图。
scatter3(data(:,1), data(:,2), data(:,3), 5, density_3D, 'filled')
hTitle = title('Density Scatter3');
hXLabel = xlabel('XAxis');
hYLabel = ylabel('YAxis');
hZLabel = zlabel('ZAxis');
4. 细节优化
为了插图的美观，将初始三维密度散点图赋上之前选择的颜色：
% 赋色
colormap(map)
colorbar
进一步，对坐标轴细节等进行美化：
% 坐标轴美化
view(45,27)
% 俯视
% view(0,90)
% set(gca,'xlim',[-8.5 8.5],...
% 'ylim',[-8.5 8.5])
set(gca, 'Box', 'on', ... % 边框
'XGrid', 'on', 'YGrid', 'on', 'ZGrid', 'on',... % 网格
'TickDir', 'out', 'TickLength', [.005 .005], ... % 刻度
'XColor', [.1 .1 .1], 'YColor', [.1 .1 .1], 'ZColor', [.1 .1 .1]) % 坐标轴颜色
% 字体和字号
set(gca, 'FontName', 'Arial', 'FontSize', 10)
set([hXLabel, hYLabel, hZLabel], 'FontSize', 11, 'FontName', 'Arial')
set(hTitle, 'FontSize', 12, 'FontWeight' , 'bold')
% 背景颜色
set(gcf,'Color',[1 1 1])
设置完毕后，以期刊所需分辨率、格式输出图片。
%% 图片输出
figW = figureWidth;
figH = figureHeight;
set(figureHandle,'PaperUnits',figureUnits);
set(figureHandle,'PaperPosition',[0 0 figW figH]);
fileout = 'test';
print(figureHandle,[fileout,'.png'],'-r300','-dpng');
也可以尝试其它配色：
以上。
转自：“阿昆的科研日常”微信公众号
如有侵权，请联系本站删除！
