这种密度散点图可以用Matlab画?

文摘   教育   2024-11-08 10:55   山东  

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

除了作图,密度散点图绘制的关键还在于密度的计算,这个我在之前的文章中有分享过(Matlab基于KD树的离散点密度特征提取)。

本期就来分享一下密度散点图的绘制方法,先来看一下成品效果:

特别提示:本期内容『数据+代码』已上传资源群中,加群的朋友请自行下载。有需要的朋友可以关注同名公号【阿昆的科研日常】,后台回复关键词【全家桶】查看加入方式


 

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);

(点击上图查看TheColor具体功能)

获取方式:公众号(阿昆的科研日常)后台回复 TC


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');

也可以尝试其它配色:

以上。


5. 全家桶

全家桶计划是指将过去以及未来发布的各种相关资源统统上传到一个群中,大家加入全家桶计划后,可以随时随地地下载群中的所有资源

全家桶中收录有……

Matlab论文插图绘制模板系列(已完结,包含3种版本*139种经典数据图,后台回复关键词【139查看单独获取方式):

Matlab进阶绘图系列除Matlab默认绘图外的复杂绘图,现包72种经典数据图):

Nature/Science论文插图复刻系列(让Matlab论文插图达到顶刊效果,现包括5种数据图):

顶刊论文高清图集系列,包括Nature高清图集(学习Nature绘图思路,包括2400张左右):

Science高清图集(学习Science绘图思路,超1000张):

Matlab配色补充包系列(包括配色强化addcolorplus、450种常用色colorplus、SCI专用配色工具Rggsci以及中国、日本、美国、英国、法国传统颜色等):

Matlab科研绘图精品工具包(包括渐变色生成/编辑器ColorMap、渐变色生成器专业版ColorMapPro、光影渲染器Shadow、顶级期刊配色拷贝工具ColorCopy和图片数据提取工具Fig2Data):

以及Matlab终极配色工具TheColor(2套):

公众号后台回复关键词【全家桶】查看加入方式。

以上。


如果你觉得我的分享对你有帮助的话,欢迎大家在这里点赞、在看、分享。当然,也欢迎大家在这里打赏。互动越多,更新越快哦~

声明:本公众号的所有原创内容,在未经允许的情况下,不得用于任何商业用途,违者必究。

阿昆的科研日常
测绘科普,经验分享,科研日常,蜜汁脑洞
 最新文章