这期做了个可能有用的小工具，一般论文中热图很少给出数据，于是就想写个小工具通过热图上的颜色估计出数据值来，目前写了个初版的工具分享给大家！

工具函数

由于只是初版，要手动改的地方还是不少的，要设置好要提取数据的图片地址，结果存储地址，数据热图行数、列数以及colorbar的范围。

function getHeatmapData
% 要提取的图像地址
picPath='1.png';
% 提取结果存储路径
HMPath='data1.mat';
% 热图行数、列数
heatMapSize=[8,8];
% colorbar所表示范围
climColorbar=[-1,1];
% 提取colorbar上颜色点数量
pntNumColorBar=20;
% 数据计算所用点数(数值越大越精确)
pntNumInterp=1000;% figure窗口初始化
pic=imread(picPath);
figure()
ax=gca;hold on
ax.XLim=[0,size(pic,2)];
ax.YLim=[0,size(pic,1)];
imshow(pic)
% 绘制红蓝色点的句柄
baHdl=plot(ax,0,0,'b+','MarkerSize',12,'LineWidth',1.5);
rxHdl=plot(ax,0,0,'rx','MarkerSize',12,'LineWidth',1.5);
pntSet=zeros(0,2);
colorSet=zeros(4,3);
set(gcf,'WindowButtonDownFcn',@buttondown)function buttondown(~,~)xy=get(gca,'CurrentPoint');xp=xy(1,2);yp=xy(1,1);pos=[yp,xp];if strcmp(get(gcf,'SelectionType'),'normal')pntSet=[pntSet;pos];elseif size(pntSet,1)>0pntSet(end,:)=[];endpntSet=round(pntSet);pntSet(pntSet<1)=1;if size(pntSet,1)>4set(gcf,'WindowButtonDownFcn',[])saveColor();close allelse% 标记采样点位置和数值baHdl.XData=pntSet(:,1);baHdl.YData=pntSet(:,2);rxHdl.XData=pntSet(:,1);rxHdl.YData=pntSet(:,2);delete(findobj('type','text'))for i=1:size(pntSet,1)colorSet(i,:)=[pic(pntSet(i,2),pntSet(i,1),1),...pic(pntSet(i,2),pntSet(i,1),2),...pic(pntSet(i,2),pntSet(i,1),3)];text(ax,pntSet(i,1),pntSet(i,2),['    ',num2str(i),':[',num2str(colorSet(i,1)),',',...num2str(colorSet(i,2)),',',num2str(colorSet(i,3)),']'])endendendfunction saveColor(~,~)% 获取colorbar上颜色cbPntSet=pntSet(1:2,:);cbPntSetdiff=abs(cbPntSet(1,:)-cbPntSet(2,:));if cbPntSetdiff(1)>cbPntSetdiff(2)tPntSet=cbPntSet(1,1)+linspace(0,1,pntNumColorBar)'.*(cbPntSet(2,1)-cbPntSet(1,1));tPntSet=round(tPntSet);CMap=[pic(round((cbPntSet(1,2)+cbPntSet(2,2))/2),tPntSet,1)',...pic(round((cbPntSet(1,2)+cbPntSet(2,2))/2),tPntSet,2)',...pic(round((cbPntSet(1,2)+cbPntSet(2,2))/2),tPntSet,3)'];elsetPntSet=cbPntSet(1,2)+linspace(0,1,pntNumColorBar)'.*(cbPntSet(2,2)-cbPntSet(1,2));tPntSet=round(tPntSet);CMap=[pic(tPntSet,round((cbPntSet(1,1)+cbPntSet(2,1))/2),1),...pic(tPntSet,round((cbPntSet(1,1)+cbPntSet(2,1))/2),2),...pic(tPntSet,round((cbPntSet(1,1)+cbPntSet(2,1))/2),3)];endCMap=double(CMap);% -----------------------------------------------------------------% 提取heatmap上颜色hmPntSet=pntSet(3:4,:);hmXLim=[min(hmPntSet(:,1)),max(hmPntSet(:,1))];hmYLim=[min(hmPntSet(:,2)),max(hmPntSet(:,2))];YList=linspace(hmYLim(1),hmYLim(2),heatMapSize(1));XList=linspace(hmXLim(1),hmXLim(2),heatMapSize(2));[XMesh,YMesh]=meshgrid(round(XList),round(YList));CMesh=zeros(heatMapSize);for i=1:heatMapSize(1)for j=1:heatMapSize(2)CMesh(i,j,1)=pic(YMesh(i,j),XMesh(i,j),1);CMesh(i,j,2)=pic(YMesh(i,j),XMesh(i,j),2);CMesh(i,j,3)=pic(YMesh(i,j),XMesh(i,j),3);endend% 计算heatmap数值范围CMesh=double(CMesh);t1=linspace(climColorbar(1),climColorbar(2),pntNumColorBar)';t2=linspace(climColorbar(1),climColorbar(2),pntNumInterp)';CMapInterp=[interp1(t1,CMap(:,1),t2,'linear'),interp1(t1,CMap(:,2),t2,'linear'),interp1(t1,CMap(:,3),t2,'linear')];Data=zeros(heatMapSize);for i=1:heatMapSize(1)for j=1:heatMapSize(2)tRGB=[CMesh(i,j,1),CMesh(i,j,2),CMesh(i,j,3)];tnorm2=sum((CMapInterp-tRGB).^2,2);[ind,~]=find(tnorm2==min(tnorm2),1);Data(i,j)=t2(ind);endend% 存储数值
%         timestr=char(datetime('now'));
%         timestr(timestr==' ')='_';
%         timestr(timestr==':')='_';
%         nowStr=[timestr,'.mat'];CMapInterp=CMapInterp./255;CMap=CMap./255;save(HMPath,'Data','CMap','CMapInterp','climColorbar')end
end

基本使用

假设有这么一张图(压缩包内提供的示例图7)

要提取的数据8行8列，colorbar范围-1到1，因此将工具函数最前面的代码做出如此改动：

% 要提取的图像地址
picPath='7.png';
% 提取结果存储路径
HMPath='data7.mat';
% 热图行数、列数
heatMapSize=[8,8];
% colorbar所表示范围
climColorbar=[-1,1];
% 提取colorbar上颜色点数量
pntNumColorBar=20;
% 数据计算所用点数(数值越大越精确)
pntNumInterp=1000;

依次点击colorbar最小值处，最大值处，以及热图的对角线两个点处，之后随便点击图片任意处(共点击五次)，就能将数据存储在mat文件中：

大概这样取四个点，并随便点击图像某处(若是点错了可鼠标右键撤回)：

之后查看mat文件：

其中：

Data 提取的热图数据
CMap colormap
CMapInterp 插值后colormap
climColor colorbar范围

获取了这些数据后，可以很轻松的靠以下代码复刻出来图像：

HM=load('data7.mat');
HMHdl=heatmap(HM.Data);
colormap(HM.CMapInterp)
% 新版本可以将caxis换为clim
caxis(HM.climColorbar)% 不显示数值
HMHdl.CellLabelColor='none';

看起来提取的很完美，虽然肯定有误差，但是对于一般的论文复现已经足够了：

行列数不相同

几乎完全相同的处理手段,比如此图是15行20列，colorbar范围是-0.6到0.8，将最前面的代码修改为：

% 要提取的图像地址
picPath='6.png';
% 提取结果存储路径
HMPath='data6.mat';
% 热图行数、列数
heatMapSize=[15,20];
% colorbar所表示范围
climColorbar=[-.6,.8];
% 提取colorbar上颜色点数量
pntNumColorBar=20;
% 数据计算所用点数(数值越大越精确)
pntNumInterp=1000;

运行后进行取点：

之后还是运行绘图代码：

HM=load('data6.mat');
HMHdl=heatmap(HM.Data);
colormap(HM.CMapInterp)
% 新版本可以将caxis换为clim
caxis(HM.climColorbar)% 不显示数值
HMHdl.CellLabelColor='none';

数据提取结果：

非方形色块热图

依旧依据热图性质修改一下参数：

% 要提取的图像地址
picPath='4.png';
% 提取结果存储路径
HMPath='data4.mat';
% 热图行数、列数
heatMapSize=[11,11];
% colorbar所表示范围
climColorbar=[-1,1];
% 提取colorbar上颜色点数量
pntNumColorBar=20;
% 数据计算所用点数(数值越大越精确)
pntNumInterp=1000;

这种圆点形状的尽量要取点时取在正中央：

数据提取结果：

色块上有文本

依旧依据热图性质修改一下参数：

% 要提取的图像地址
picPath='8.png';
% 提取结果存储路径
HMPath='data8.mat';
% 热图行数、列数
heatMapSize=[8,8];
% colorbar所表示范围
climColorbar=[-1,1];
% 提取colorbar上颜色点数量
pntNumColorBar=20;
% 数据计算所用点数(数值越大越精确)
pntNumInterp=1000;

有文本的热图，在热图色块上点击时，尽量俩点点击在同一位置且不要碰到文字，比如俩色块都点击同一边角（下图所示即为点击两个色块左下角）：

提取效果：

完

以上已给出完整代码，想要更多例子练练手：

【热图数据提取】

更新时会跟进更新以下连接：
【链接】：
https://pan.baidu.com/s/1zoStpQ54LqmR7fSXMWWyOA?pwd=slan

【提取码】：slan

MATLAB | 如何从热图中提取数据相关推荐

MATLAB之从fig文件中提取数据
1 如果你的fig文件中的图像每个像素都有数据,可以通过以下方式获得图像每个点的值,输出data是矩阵,大小是图像像素的行列数 open('figname.fig'); h=get(gca,'Chil ...
matlab从fig图文件中提取数据
matlab从fig图文件中提取数据 lh=findall(gca,'type','line'); cell_l = cellfun('length',get(lh,'xdata')); i_lh=f ...
MATLAB figure中提取数据
MATLAB figure中提取数据 (2011-10-26 14:26:21) 转载▼ 标签: 杂谈分类: matlab figure画出来,提取数据有很多好处,方便保存,计算,加工,还可以导入到 ...
matlab程序模拟汽车理论中的数据特性图动力特性图,matlab程序模拟汽车理论中的数据特性图1.doc...
matlab程序模拟汽车理论中的数据特性图1.doc 车辆理论作业第三组第一大题确定一辆轻型货车的动力性能1.绘制发动机的使用外特性曲线N60014000NN/1000T1931329527N/100 ...
如何从卫星图中提取水系数据
01 概述相信大家之前有看过很多通过电子地图提取水系的教程,该方法虽然可以提取水系,但是有一定的限制,如果在一些比较偏远的地方,电子地图上没有水系数据就无法提取,存在一定的局限性,之前已经为大家介绍 ...
如何用matlab将交通实时路况图中路从图中提取出来，保留颜色
如何用matlab将交通实时路况图中路从图中提取出来,保留颜色欢迎使用Markdown编辑器你好! 这是你第一次使用 Markdown编辑器所展示的欢迎页.如果你想学习如何使用Markdown编 ...
从.fig文件中提取数据
matlab的绘图窗口可以保存为.fig文件,方便随后的修改和再加工.那么,当需要提取图中曲线的数据时,该怎么办?(2018b版matlab) 平移.缩放,使图中所有数据处于图窗内: 工具->刷 ...
【Python】使用日历热图进行时序数据可视化
相信很多人都会在 Github 中看到这么一个热图,该热图记录的是 Github 平台使用的日常贡献.在每个日历年的热图中以天为单位采样的时间序列数据.GitHub 的贡献图表示用户在过去几年中所做的 ...
matlab读Excel表格数据画图,matlab读Excel表格数据画图-matlab如何从excel表格中读取数据？...
如何将excel表格中大量数据导入matlab中并作图 1.打开matlab,点击主页下面的数据导入,你可以导入excel数据,这里导入你自己命名的huitushuju文件. 2.单击"打开 ...

MATLAB | 如何从热图中提取数据