脚本之家 服务器常用软件
快捷导航
您的位置：首页 → 软件编程 → C 语言 → Matlab山脊图

基于Matlab实现山脊图的绘制

 更新时间：2022年05月10日 09:22:27   作者：slandarer  
这篇文章主要介绍了如何利用Matlab实现山脊图的绘制，文中的示例代码讲解详细，对我们学习Matlab有一定的帮助，需要的可以参考一下

写了一个用于绘制山脊图的模板，仅需要往需要填写数据的区域填入数据点击运行即可，以下提供两款模板，第一款为纯色模板而第二款为渐变色模板：

纯色山脊图模板

模板代码：

function ridgeMapTMPL1
% @author: slandarer

% 在这里放入你的数据=======================================================
X1=normrnd(2,2,1,1000);
X2=[normrnd(4,4,1,1000),normrnd(5,2,1,200)];
X3=[normrnd(2.5,2,1,1000),normrnd(6,4,1,200)];
X4=[normrnd(1,1,1,300),normrnd(2,4,1,200)];
X5=[normrnd(4,2,1,300),normrnd(2,4,1,600)];

% 把数据放到元胞数组，要是数据太多可写循环放入
dataCell={X1,X2,X3,X4,X5};    

% 各个数据类的名称，可空着
dataName={'A','B','C','D','E'};

% 各个山脊的颜色，可空着也可只放一个颜色
colorList=[];                  

% 非必要属性
% -------------------------------------------------------------------------
legendOn=false;  % 是否绘制图例 true/false
xTickOn=true;    % 开启X轴

sep=1/6;         % 设置山脊距离，可空着
xLim=[];         % 设置X轴范围距离，可空着
axColor=[];      % 设置背景颜色，可空着
fontName='';     % 设置X，Y轴标签字体，可空着
fontSize=14;     % 设置X，Y轴标签字号，可空着
faceAlpha=1;           % 面透明度
edgeColor=[1,1,1].*.3; % 边缘颜色
% =========================================================================

% 设置颜色
classNum=length(dataCell);
if size(colorList,1)==0
    colorList=repmat([130,170,172]./255,[classNum,1]);
else
    colorList=repmat(colorList,[ceil(classNum/size(colorList,1)),1]);
end

% 设置间隙距离
if isempty(sep)
sep=1/6;
end

hold on
ax=gca;
ax.YGrid='on';
ax.YLim=[0,sep*classNum+sep/2];
ax.YTick=0:sep:sep*(classNum-1);
ax.YColor='none';
ax.LineWidth=1.2;
% ax.YTickLabel=dataName;
if isempty(fontName),fontName='Helvetica';end
if isempty(fontSize),fontSize=14;end
ax.FontName=fontName;
ax.FontSize=fontSize;
if ~xTickOn,ax.XColor='none';end
if ~isempty(xLim),ax.XLim=xLim;end
if ~isempty(axColor),ax.Color=axColor;end
if isempty(dataName)
    for i=1:classNum
        dataName{i}=num2str(i);
    end
end

% 绘制山脊
for i=1:classNum
    tX=dataCell{i};tX=tX(:);
    [F,Xi]=ksdensity(tX);
    patchCell(i)=fill([Xi(1),Xi,Xi(end)],[0,F,0]+(sep).*(classNum-i).*ones(1,length(F)+2),...
        colorList(i,:),'EdgeColor','none','FaceAlpha',faceAlpha);
    plot([Xi(1),Xi,Xi(end)],[0,F,0]+(sep).*(classNum-i).*ones(1,length(F)+2),...
        'Color',edgeColor,'LineWidth',1.2)
end

% 绘制图例
if legendOn
legend(patchCell,dataName)
end

% 绘制字符
ax.UserData.classNum=classNum;
ax.UserData.sep=sep;
for k=1:classNum
    ax.UserData.(['t',num2str(k)])=text(ax.XLim(1),(sep).*(classNum-k),[dataName{k},' '],...
        'FontSize',fontSize,'FontName',fontName,'HorizontalAlignment','right','VerticalAlignment','bottom');
end
    function reTXT(~,~)
        for kk=1:ax.UserData.classNum
            ax.UserData.(['t',num2str(kk)]).Position=...
                [ax.XLim(1),(ax.UserData.sep).*(ax.UserData.classNum-kk),0];
        end
    end

set(ax.Parent,'WindowButtonMotionFcn',@reTXT); 
% 额外的属性设置===========================================================
% ax.(...)=...

% =========================================================================
end

绘制效果：

更该部分配置后的绘制效果：

将每个山脊赋予不同颜色，删除X轴并添加图例

% 在这里放入你的数据=======================================================
X1=normrnd(2,2,1,1000);
X2=[normrnd(4,4,1,1000),normrnd(5,2,1,200)];
X3=[normrnd(2.5,2,1,1000),normrnd(6,4,1,200)];
X4=[normrnd(1,1,1,300),normrnd(2,4,1,200)];
X5=[normrnd(4,2,1,300),normrnd(2,4,1,600)];

% 把数据放到元胞数组，要是数据太多可写循环放入
dataCell={X1,X2,X3,X4,X5};    

% 各个数据类的名称，可空着
dataName={'A','B','C','D','E'};

% 各个山脊的颜色，可空着也可只放一个颜色
colorList=[0.98,0.70,0.95
0.69,0.81,1.00
0.50,0.88,0.89
0.50,0.87,0.61
0.86,0.82,0.50];                  

% 非必要属性
% -------------------------------------------------------------------------
legendOn=true;  % 是否绘制图例 true/false
xTickOn=false;   % 开启X轴

sep=1/6;         % 设置山脊距离，可空着
xLim=[];         % 设置X轴范围距离，可空着
axColor=[];      % 设置背景颜色，可空着
fontName=[];     % 设置X，Y轴标签字体，可空着
fontSize=14;     % 设置X，Y轴标签字号，可空着
faceAlpha=1;           % 面透明度
edgeColor=[1,1,1].*.3; % 边缘颜色
% ===============================================================

渐变色山脊图模板

模板代码：

function ridgeMapTMPL2
% @author: slandarer

% 在这里放入你的数据=======================================================
X1=normrnd(2,2,1,1000);
X2=[normrnd(4,4,1,1000),normrnd(5,2,1,200)];
X3=[normrnd(2.5,2,1,1000),normrnd(6,4,1,200)];
X4=[normrnd(1,1,1,300),normrnd(2,4,1,200)];
X5=[normrnd(4,2,1,300),normrnd(2,4,1,600)];

% 把数据放到元胞数组，要是数据太多可写循环放入
dataCell={X1,X2,X3,X4,X5};    

% 各个数据类的名称，可空着
dataName={'A','B','C','D','E'};

% 山脊的渐变颜色，可空着也放数组也可放颜色名称
cmap=[];
% cmap=PYCM().plasma();   
% cmap='colorcube'

% 非必要属性
% -------------------------------------------------------------------------
xTickOn=true;    % 开启X轴

sep=1/6;         % 设置山脊距离，可空着
xLim=[];         % 设置X轴范围距离，可空着
fontName='';     % 设置X，Y轴标签字体，可空着
fontSize=14;     % 设置X，Y轴标签字号，可空着
% =========================================================================


classNum=length(dataCell);

% 设置间隙距离
if isempty(sep)
sep=1/6;
end

hold on
ax=gca;
ax.YGrid='on';
ax.YLim=[0,sep*classNum+sep/2];
ax.YTick=0:sep:sep*(classNum-1);
ax.YColor='none';
ax.LineWidth=1.2;
% ax.YTickLabel=dataName;
if isempty(fontName),fontName='Helvetica';end
if isempty(fontSize),fontSize=14;end

ax.FontName=fontName;
ax.FontSize=fontSize;
if ~xTickOn,ax.XColor='none';end
if ~isempty(xLim),ax.XLim=xLim;end
if isempty(dataName)
    for i=1:classNum
        dataName{i}=num2str(i);
    end
end

% 绘制山脊
for i=1:classNum
    tX=dataCell{i};tX=tX(:);
    [F,Xi]=ksdensity(tX);
    patchCell(i)=patch([Xi(1),Xi,Xi(end)],[0,F,0]+(sep).*(classNum-i).*ones(1,length(F)+2),...
        [Xi(1),Xi,Xi(end)],'FaceColor','interp','EdgeColor','none');
    plot([Xi(1),Xi,Xi(end)],[0,F,0]+(sep).*(classNum-i).*ones(1,length(F)+2),...
        'Color',[.3,.3,.3],'LineWidth',1.2)
end
if isempty(cmap)
colormap()
else
colormap(cmap)
end

% 绘制字符
ax.UserData.classNum=classNum;
ax.UserData.sep=sep;
for k=1:classNum
    ax.UserData.(['t',num2str(k)])=text(ax.XLim(1),(sep).*(classNum-k),[dataName{k},' '],...
        'FontSize',fontSize,'FontName',fontName,'HorizontalAlignment','right','VerticalAlignment','bottom');
end
    function reTXT(~,~)
        for kk=1:ax.UserData.classNum
            ax.UserData.(['t',num2str(kk)]).Position=...
                [ax.XLim(1),(ax.UserData.sep).*(ax.UserData.classNum-kk),0];
        end
    end

set(ax.Parent,'WindowButtonMotionFcn',@reTXT); 
% 额外的属性设置===========================================================
% ax.(...)=...

% =========================================================================
end

绘制效果：

使用MATLAB自带颜色：

例：将cmap=[]改为：

cmap='summer'

cmap='colorcube'

使用PYCM函数：

例：将cmap=[]改为：

cmap=PYCM().plasma()

到此这篇关于基于Matlab实现山脊图的绘制的文章就介绍到这了,更多相关Matlab山脊图内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

您可能感兴趣的文章:

相关文章

  • C语言用栈和队列实现的回文检测功能示例

    C语言用栈和队列实现的回文检测功能示例

    这篇文章主要介绍了C语言用栈和队列实现的回文检测功能,结合具体实例形式分析了C语言栈和队列的定义及使用栈和队列进行回文检测的操作技巧,需要的朋友可以参考下
    2017-06-06
  • C语言位段(位域)机制结构体的特殊实现及解析

    C语言位段(位域)机制结构体的特殊实现及解析

    这篇文章主要为大家介绍了C语言位段位域机制结构体的特殊实现讲解有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步早日升职加薪
    2022-02-02
  • C++超详细分析函数重载的使用

    C++超详细分析函数重载的使用

    C++ 允许多个函数拥有相同的名字，只要它们的参数列表不同就可以，这就是函数的重载（Function Overloading），借助重载，一个函数名可以有多种用途
    2022-04-04
  • C++使用string的大数加法运算（1）

    C++使用string的大数加法运算（1）

    这篇文章主要为大家详细介绍了C++使用string的大数加法运算，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
    2019-09-09
  • 算法详解之分治法具体实现

    算法详解之分治法具体实现

    这篇文章主要介绍了算法详解之分治法具体实现,需要的朋友可以参考下
    2014-02-02
  • 详解C++句柄类

    详解C++句柄类

    本篇文章给大家详细分析了C++句柄类的相关知识点，对此有需要的朋友跟着学习参考下吧。
    2018-06-06
  • C语言构建连连看游戏(矩阵方式)

    C语言构建连连看游戏(矩阵方式)

    这篇文章主要为大家详细介绍了C语言构建连连看游戏，采用矩阵方式，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
    2021-09-09
  • C语言简易扫雷游戏

    C语言简易扫雷游戏

    这篇文章主要为大家详细介绍了C语言简易扫雷游戏，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
    2019-10-10
  • 函数指针与指针函数的学习总结

    函数指针与指针函数的学习总结

    函数指针是指向函数的指针，指针函数是指一个函数的返回值是一个指针。以下就是对函数指针与指针函数的应用进行了详细的分析介绍，需要的朋友可以参考下
    2013-07-07
  • C/C++编译报错printf was not declared in this scope问题及解决

    C/C++编译报错printf was not declared in 

    这篇文章主要介绍了C/C++编译报错printf was not declared in this scope问题及解决方案,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教
    2023-08-08

最新评论