MATLAB | 这些花里胡哨的热图怎么画
2024-05-16 02:56:04
好早之前写过一个绘制相关系数矩阵的代码，但是会自动求相关系数，而且画出来的热图只能是方形，这里写一款允许nan值出现，任意形状的热图绘制代码，绘制效果如下：
如遇到bug请后台提出，并去gitee下载更新版本
教程部分
1 基础绘制
假设有个随机数矩阵，绘图：
Data=rand(15,15);
SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();
假设存在负数：
Data=rand(15,15)-.5;
SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();
假设存在nan:
Data=rand(12,12)-.5;
Data([4,5,13])=nan;
SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();
假设非方阵：
Data=rand(25,30);
SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();
颜色怎么换，在后面会讲到。
2 显示文字
使用setText方法显示文字：
Data=rand(12,12)-.5;
Data([4,5,13])=nan;SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();SHM.setText();
如果要修饰或显示文字请在修改颜色后进行。
3 修改形状
形状有以下几个可选项：
- sq : 方形(默认)
- pie : 饼图
- circ : 圆形
- oval : 椭圆形
- hex : 六边形
- asq : 自带调整大小的方形
- acirc : 自带调整大小的圆形
figure()
Data=rand(12,12)-.5;
SHM=SHeatmap(Data,'Format','oval');
SHM=SHM.draw(); figure()
Data=rand(12,12)-.5;
SHM=SHeatmap(Data,'Format','pie');
SHM=SHM.draw();
展示一下全部类型：
sq
pie
circ
oval
hex
asq
acirc
4 颜色范围
新版本使用clim老版本使用caxis调整映射范围：
Data=rand(12,12)-.5;
Data([4,5,13])=nan;SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();clim([-.8,.8])
SHM.setText();
调整前：
调整后：
5 colormap
matlab自带的colormap都能用：
Data=rand(14,14);
SHM=SHeatmap(Data,'Format','sq');
SHM.draw();colormap(bone)
可以联合我写的slanCM工具函数：
https://slandarer.blog.csdn.net/article/details/127719784
举点例子：
colormap(flipud(slanCM(20)))
colormap(flipud(slanCM(21)))
colormap(slanCM(61))
colormap(slanCM(177))
再比如：
Data=rand(12,12)-.5;
Data([4,5,13])=nan;SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();clim([-.7,.7])
colormap(slanCM(141))
SHM.setText();
clim([-.7,.7])
colormap(slanCM(136))
clim([-.7,.7])
colormap(slanCM(134))
clim([-.8,.8])
colormap(slanCM(110))
clim([-1,1])
colormap(slanCM(97))
6 调整colormap位置
假设编写如下代码：
Data=rand(3,12);
SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();
改成这样：
Data=rand(3,12);
SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();CB=colorbar;
CB.Location='southoutside';
7 修改标签
改编x轴y轴标签就还是设置XTickLabel
和YTickLabel:
Data=rand(10,10);
SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw(); ax=gca;
ax.XTickLabel={'X-1','X-2','X-3','X-4','X-5','X-6','X-7','X-8','X-9','X-10'};
ax.YTickLabel={'Y-1','Y-2','Y-3','Y-4','Y-5','Y-6','Y-7','Y-8','Y-9','Y-10'};
ax.FontSize=14;
8 整体修饰
使用：
- setBox ：修饰边框
- setPatch ：修饰图形
- setText ：修饰文本
可对各个部分进行整体修改。
修改图形边缘及框颜色：
Data=rand(10,10);
Data([4,5,13])=nan;SHM=SHeatmap(Data,'Format','pie');
SHM=SHM.draw(); SHM.setBox('Color',[0,0,.8])
SHM.setPatch('EdgeColor',[.8,0,0])
修改全部文本颜色：
Data=rand(10,10);
Data([4,5,13])=nan;SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw(); SHM.setText('Color',[0,0,.8],'FontSize',14)
9 局部修改
使用：
- setTextMN
- setPatchMN
设置第m行n列对象的属性：
举个例子，把数值大于0.9部分文字变为星号，并用红色框框起来：
Data=rand(9,9);
Data([4,5,13])=nan;SHM=SHeatmap(Data,'Format','sq');
SHM=SHM.draw();SHM.setText(); for i=1:size(Data,1)for j=1:size(Data,2)if Data(i,j)>.9SHM.setTextMN(i,j,'String','**','FontSize',20)SHM.setPatchMN(i,j,'EdgeColor',[1,0,0],'LineWidth',2)endend
end
局部设置可设置nan部分颜色：
SHM.setPatchMN(4,1,'FaceColor',[.8,.6,.6])
工具函数完整代码
classdef SHeatmap
% =========================================================================
% Format
% -------------------------------------------------------------------------
% sq : 方形(默认)
% pie : 饼图
% circ : 圆形
% oval : 椭圆形
% hex : 六边形
% asq : 自带调整大小的方形
% acirc : 自带调整大小的圆形
% =========================================================================
% @author : slandarer
% gzh : slandarer随笔
% -------------------------------------------------------------------------
% Zhaoxu Liu / slandarer (2023). special heatmap
% (https://www.mathworks.com/matlabcentral/fileexchange/125520-special-heatmap),
% MATLAB Central File Exchange. 检索来源 2023/3/1.
% -------------------------------------------------------------------------propertiesax,arginList={'Format'}Format='sq' % sq : 方形(默认)% pie : 饼图% circ : 圆形% oval : 椭圆形% hex : 六边形% asq : 自带调整大小的方形% acirc : 自带调整大小的圆形DatadfColor1=[0.9686 0.9882 0.9412; 0.9454 0.9791 0.9199; 0.9221 0.9700 0.8987; 0.8988 0.9609 0.8774;0.8759 0.9519 0.8560; 0.8557 0.9438 0.8338; 0.8354 0.9357 0.8115; 0.8152 0.9276 0.7892;0.7909 0.9180 0.7685; 0.7545 0.9039 0.7523; 0.7180 0.8897 0.7361; 0.6816 0.8755 0.7199;0.6417 0.8602 0.7155; 0.5962 0.8430 0.7307; 0.5507 0.8258 0.7459; 0.5051 0.8086 0.7610;0.4596 0.7873 0.7762; 0.4140 0.7620 0.7914; 0.3685 0.7367 0.8066; 0.3230 0.7114 0.8218;0.2837 0.6773 0.8142; 0.2483 0.6378 0.7929; 0.2129 0.5984 0.7717; 0.1775 0.5589 0.7504;0.1421 0.5217 0.7314; 0.1066 0.4853 0.7132; 0.0712 0.4488 0.6950; 0.0358 0.4124 0.6768;0.0314 0.3724 0.6364; 0.0314 0.3319 0.5929; 0.0314 0.2915 0.5494; 0.0314 0.2510 0.5059]dfColor2=[0.6196 0.0039 0.2588; 0.6892 0.0811 0.2753; 0.7588 0.1583 0.2917; 0.8283 0.2354 0.3082;0.8706 0.2966 0.2961; 0.9098 0.3561 0.2810; 0.9490 0.4156 0.2658; 0.9660 0.4932 0.2931;0.9774 0.5755 0.3311; 0.9887 0.6577 0.3690; 0.9930 0.7266 0.4176; 0.9943 0.7899 0.4707;0.9956 0.8531 0.5238; 0.9968 0.9020 0.5846; 0.9981 0.9412 0.6503; 0.9994 0.9804 0.7161;0.9842 0.9937 0.7244; 0.9526 0.9810 0.6750; 0.9209 0.9684 0.6257; 0.8721 0.9486 0.6022;0.7975 0.9183 0.6173; 0.7228 0.8879 0.6325; 0.6444 0.8564 0.6435; 0.5571 0.8223 0.6448;0.4698 0.7881 0.6460; 0.3868 0.7461 0.6531; 0.3211 0.6727 0.6835; 0.2553 0.5994 0.7139;0.2016 0.5261 0.7378; 0.2573 0.4540 0.7036; 0.3130 0.3819 0.6694; 0.3686 0.3098 0.6353]Colormap;maxVpatchHdl;boxHdl;pieHdl;textHdlendmethodsfunction obj=SHeatmap(Data,varargin)obj.Data=Data;obj.maxV=max(max(abs(obj.Data)));% 获取其他数据disp(char([64 97 117 116 104 111 114 32 58 32,...115 108 97 110 100 97 114 101 114]))for i=1:2:(length(varargin)-1)tid=ismember(obj.arginList,varargin{i});if any(tid)obj.(obj.arginList{tid})=varargin{i+1};endend% 设置配色if any(any(obj.Data<0))obj.Colormap=obj.dfColor2;% tX=linspace(0,1,size(obj.Colormap,1));% tXi=linspace(0,1,256);% tR=interp1(tX,obj.Colormap(:,1),tXi);% tG=interp1(tX,obj.Colormap(:,2),tXi);% tB=interp1(tX,obj.Colormap(:,3),tXi);% obj.Colormap=[tR(:),tG(:),tB(:)];elseobj.Colormap=obj.dfColor1(end:-1:1,:);endendfunction obj=draw(obj)obj.ax=gca;hold on;box onobj.ax.FontName='Times New Roman';obj.ax.FontSize=12;obj.ax.LineWidth=.8;obj.ax.XLim=[.5,size(obj.Data,2)+.5];obj.ax.YLim=[.5,size(obj.Data,1)+.5];obj.ax.YDir='reverse';obj.ax.TickDir='out';obj.ax.TickLength=[0.002,0.002];obj.ax.DataAspectRatio=[1,1,1];obj.ax.YTick=1:size(obj.Data,1);obj.ax.XTick=1:size(obj.Data,2);colormap(obj.Colormap)colorbarif any(any(obj.Data<0))try caxis(obj.maxV.*[-1,1]),catch,endtry clim(obj.maxV.*[-1,1]),catch,endelsetry caxis(obj.maxV.*[0,1]),catch,endtry clim(obj.maxV.*[0,1]),catch,endend% 调整初始界面大小fig=obj.ax.Parent;fig.Color=[1,1,1];if max(fig.Position(3:4))<600fig.Position(3:4)=[1.6,1.8].*fig.Position(3:4);fig.Position(1:2)=fig.Position(1:2)./4;endbX1=repmat([.5,size(obj.Data,2)+.5,nan],[size(obj.Data,1)+1,1])';bY1=repmat((.5:1:(size(obj.Data,1)+.5))',[1,3])';bX2=repmat((.5:1:(size(obj.Data,2)+.5))',[1,3])';bY2=repmat([.5,size(obj.Data,1)+.5,nan],[size(obj.Data,2)+1,1])';obj.boxHdl=plot([bX1(:);bX2(:)],[bY1(:);bY2(:)],'LineWidth',.8,'Color',[1,1,1].*.85);if isequal(obj.Format,'sq')set(obj.boxHdl,'Color',[1,1,1,0])endbaseT=linspace(0,2*pi,200);hexT=linspace(0,2*pi,7);thetaMat=[1,-1;1,1].*sqrt(2)./2;for row=1:size(obj.Data,1)for col=1:size(obj.Data,2) if isnan(obj.Data(row,col))obj.patchHdl(row,col)=fill([-.5,.5,.5,-.5].*.98+col,[-.5,-.5,.5,.5].*.98+row,[.8,.8,.8],'EdgeColor','none');obj.pieHdl(row,col)=fill([0,0,0,0],[0,0,0,0],[0,0,0]);obj.textHdl(row,col)=text(col,row,'×','FontName','Times New Roman','HorizontalAlignment','center','FontSize',20);elsetRatio=abs(obj.Data(row,col))./obj.maxV;switch obj.Formatcase 'sq'obj.patchHdl(row,col)=fill([-.5,.5,.5,-.5].*.98+col,[-.5,-.5,.5,.5].*.98+row,...obj.Data(row,col),'EdgeColor','none');case 'asq'obj.patchHdl(row,col)=fill([-.5,.5,.5,-.5].*.98.*tRatio+col,[-.5,-.5,.5,.5].*.98.*tRatio+row,...obj.Data(row,col),'EdgeColor','none');case 'pie'baseCircX=cos(baseT).*.92.*.5;baseCircY=sin(baseT).*.92.*.5;obj.pieHdl(row,col)=fill(baseCircX+col,baseCircY+row,...[1,1,1],'EdgeColor',[1,1,1].*.3,'LineWidth',.8);baseTheta=linspace(pi/2,pi/2+obj.Data(row,col)./obj.maxV.*2.*pi,200);basePieX=[0,cos(baseTheta).*.92.*.5];basePieY=[0,sin(baseTheta).*.92.*.5];obj.patchHdl(row,col)=fill(basePieX+col,-basePieY+row,...obj.Data(row,col),'EdgeColor',[1,1,1].*.3,'lineWidth',.8);case 'circ'baseCircX=cos(baseT).*.92.*.5;baseCircY=sin(baseT).*.92.*.5;obj.patchHdl(row,col)=fill(baseCircX+col,baseCircY+row,...obj.Data(row,col),'EdgeColor','none','lineWidth',.8);case 'acirc'baseCircX=cos(baseT).*.92.*.5;baseCircY=sin(baseT).*.92.*.5;obj.patchHdl(row,col)=fill(baseCircX.*tRatio+col,baseCircY.*tRatio+row,...obj.Data(row,col),'EdgeColor','none','lineWidth',.8);case 'oval'tValue=obj.Data(row,col)./obj.maxV;baseA=1+(tValue<=0).*tValue;baseB=1-(tValue>=0).*tValue;baseOvalX=cos(baseT).*.98.*.5.*baseA;baseOvalY=sin(baseT).*.98.*.5.*baseB;baseOvalXY=thetaMat*[baseOvalX;baseOvalY];obj.patchHdl(row,col)=fill(baseOvalXY(1,:)+col,-baseOvalXY(2,:)+row,...obj.Data(row,col),'EdgeColor',[1,1,1].*.3,'lineWidth',.8);case 'hex'obj.patchHdl(row,col)=fill(cos(hexT).*.5.*.98.*tRatio+col,sin(hexT).*.5.*.98.*tRatio+row,...obj.Data(row,col),'EdgeColor',[1,1,1].*.3,'lineWidth',.8);endobj.textHdl(row,col)=text(col,row,sprintf('%.2f',obj.Data(row,col)),'FontName','Times New Roman','HorizontalAlignment','center','Visible','off');endendendend% 修饰文本function setText(obj,varargin)graymap=mean(get(obj.ax,'Colormap'),2);climit=get(obj.ax,'CLim');for row=1:size(obj.Data,1)for col=1:size(obj.Data,2) set(obj.textHdl(row,col),'Visible','on','Color',...[1,1,1].*(interp1(linspace(climit(1),climit(2),size(graymap,1)),graymap,obj.Data(row,col))<.5),varargin{:})endendendfunction setTextMN(obj,m,n,varargin)set(obj.textHdl(m,n),varargin{:})end% 设置图形样式function setPatch(obj,varargin)for row=1:size(obj.Data,1)for col=1:size(obj.Data,2) if ~isnan(obj.Data(row,col))set(obj.patchHdl(row,col),varargin{:})if isequal(obj.Format,'pie')set(obj.pieHdl(row,col),varargin{:}) endendendendendfunction setPatchMN(obj,m,n,varargin)set(obj.patchHdl(m,n),varargin{:})if isequal(obj.Format,'pie')set(obj.pieHdl(m,n),varargin{:})endend% 设置框样式function setBox(obj,varargin)set(obj.boxHdl,varargin{:})endend
% =========================================================================
% @author : slandarer
% gzh : slandarer随笔
% -------------------------------------------------------------------------
% Zhaoxu Liu / slandarer (2023). special heatmap
% (https://www.mathworks.com/matlabcentral/fileexchange/125520-special-heatmap),
% MATLAB Central File Exchange. 检索来源 2023/3/1.
% -------------------------------------------------------------------------
end
完
编程不易，希望大家该点赞的点赞，该在看的在看！！
未经允许本代码请勿作商业用途，引用的话可以引用我file exchange上的链接，可使用如下格式：
Zhaoxu Liu / slandarer (2023). special heatmap (https://www.mathworks.com/matlabcentral/fileexchange/125520-special-heatmap), MATLAB Central File Exchange. 检索来源 2023/3/1.
若转载请保留以上file exchange链接及本文链接！！！！!
MATLAB | 这些花里胡哨的热图怎么画相关推荐
- MATLAB | 如何用MATLAB绘制花里胡哨的山脊图
本期推送教大家如何绘制各种样式的山脊图,这里做了一个工具函数用来实现好看的山脊图的绘制,编写不易请多多点赞,大体绘制效果如下: 依旧工具函数放在文末. 教程部分 0 数据准备 数据为多个一维向量放在元 ...
- matlab 修改heatmap colorbar 热图
坑 1.准好数据 2.点开color map 3.打开属性检查器 4.右击color bar 5. 手工交互式调整heatmap ================================ 还有 ...
- r语言聚类分析_技术贴 | R语言pheatmap聚类分析和热图
点击蓝字↑↑↑"微生态",轻松关注不迷路 本文由阿童木根据实践经验而整理,希望对大家有帮助. 原创微文,欢迎转发转载. 导读 pheatmap默认会对输入矩阵数据的行和列同时进行聚 ...
- origin做相关性分析图_高分文章的相关性热图，十分钟教你学会无代码复现
相关性热图可以便捷的显示多个变量之间的相关性,因此在组学数据分析中得到广泛应用,例如基因的共表达分析.样本重复检验.微生物群落的共发生网络分析等. 在昨天推文下就有人问热图怎么画?画热图最常规的办法自 ...
- 差异表达基因热图怎么看_多变的热图1（新手专用）
热图(heatmap)用不同的颜色和颜色的深浅来直观的展示数据之间的差异.在测序类的文章里,几乎必有一幅热图用来展示差异表达基因.很多工具都可以完成热图的制作,今天这篇文章主要介绍利用R语言的 phe ...
- 使用Matlab画出复杂网络博弈中热图(三维图)(学术论文)
如图所示,有的时候叫热图(可能根据不同的色彩分别吧),下面简单介绍下,有点累人的图.(Matlab所作,后面会给出代码) 具体需要看这幅图的横坐标与纵坐标的粒度是多少,比如这幅图,每幅子图选取的粒度为 ...
- matlab画动态聚类图,聚类热图，你用什么软件画？
原标题:聚类热图,你用什么软件画? 聚类热图是我们展示芯片.或者测序结果比较常用的方式,在文章中出现的频率非常高,一般来说,这张图的每行对应一个基因,每列对应一个样本,行和列形成的是基因数*样本数的表 ...
- MATLAB画聚类热图,使用ComplexHeatmap简单绘制热图
ComplexHeatmap其实是一个很全面的R包,它除了可以绘制简单热图还有其他复杂实用的热图,这里主要简单的介绍一下如何用这个R包来绘制简单热图 1. 首先按照http://bioconducto ...
- MATLAB画聚类热图,使用ComplexHeatmap绘制热图
初学R,慢慢摸索. 安装: if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages(" ...
最新文章
- 线段树 ---- 牛客多校4 ETree Xor 区间异或分段
- list存入mysql乱序_MySQL案例-并行复制乱序提交引起的同步异常
- 你不知道的JavaScript运算符
- sql数据迁移到oracle数据库,从Oracle到SQL Server数据库主键的迁移
- 《Java程序性能优化》之设计优化
- php中冒号是什么,在PHP中：(双冒号)和-(箭头)有什么区别？
- 最长公共前缀_最长的公共前缀
- c语言用switch字母判断星期几,c语言程序，输入年月日，判断这一天是这一年的第几天，同时判断这一天是星期几。（用switch语句）...
- D3 Geographies
- Feature flag，一个让软件发布轻松不掉发的神技
- Raid、lvm知识
- Camera Shutter Sound can't restore to default
- 部署Java Web项目到Heroku
- spark 程序配置日志级别 通过log4j
- kettle 报错【Maximum wait time of 10 seconds exceed while acquiring lock】
- 基于 SaaS 的 RADIUS 认证在无线网络准入中的应用和优化
- 如何使用逆分布函数模拟MEDIAN（）聚合函数
- 【微机原理与接口技术】学习笔记5 I/O接口和并行接口芯片8255A
- leanote 蚂蚁云笔记 源码编译
- 标准引领市场：四方光电参编《便携式甲醛检测仪》团体标准
热门文章
- 告别低效，戴尔官网高效编程电脑OptiPlex低至5折购！换新不等待！
- 教你关闭iTunes的自动更新功能
- Win7无法打开本地连接属性如何解决？
- Hololens连接mysql_Hololens开发笔记：UDP接收数据
- wfb-ng Release 23.01镜像无头烧录配置(1)
- 名编辑电子杂志大师教程 | 电子杂志打印功能
- 海马汽车经销商管理系统技术解析（六）顾客生日祝贺
- 聊聊最近使用Electron开发互动直播的一点心得
- 计算机 毕设课设 之含文档+PPT+源码等]精品微信小程序springboot服装企业人事管理系统+后台管理系
- 迭代硬阈值类算法总结||IHT/NIHT/CGIHT/HTP