• MATLAB | 全网唯一,使用MATLAB绘制好看的韦恩图(venn)


    韦恩图绘制的代码MATHWORKS上也有,不过最多只支持四组绘制,于是我参考R语言ggVennDiagram包制作了一款支持7组样本的韦恩图绘制MATLAB函数

    注意:仅代码无法运行,需要导入.mat格式的素材包,素材包及m文件获取放在文末,

    另:由于代码开发基于polyshape对象,本韦恩图绘制函数理论上需要至少R2017b及之后版本才能使用。

    由于工具函数过长,将被放在最后展示,以下将先展示函数用法。

    示例图


    使用教程

    1 数据格式

    数值格式需要如下形式的数值向量,目前还不支持字符串元胞数组格式,请自行将其改编为数值向量,如下例子所示若数值400同时出现在X1及X2中,即说明该元素属于X1,X2的交集。

    X1=randi([1,500],[100,1]);
    X2=randi([1,500],[100,1]);
    X3=randi([1,500],[100,1]);
    X4=randi([1,500],[100,1]);
    X5=randi([1,500],[100,1]);
    X6=randi([1,500],[100,1]);
    X7=randi([1,500],[100,1]);
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7

    2 基本使用

    以下两种形式均可:

    XX={X1,X2,X3,X4,X5,X6,X7};
    
    VN=venn(XX{1:7});
    VN=VN.draw();
    
    • 1
    • 2
    • 3
    • 4
    VN=venn(X1,X2,X3,X4,X5,X6,X7);
    VN=VN.draw();
    
    • 1
    • 2

    3 添加标签

    draw函数绘图前使用labels函数添加标签:

    VN=venn(X1,X2,X3,X4,X5,X6,X7);
    VN=VN.labels('AAA','BBB','CCC','DDD','EEE','FFF','GGG');
    VN=VN.draw();
    
    • 1
    • 2
    • 3

    4 样本数量

    支持样本数量为2-7的任意值,例如绘制两样本则使用如下格式调用:

    VN=venn(X1,X2);
    VN=VN.draw();
    
    • 1
    • 2

    以下展示2-7组样本绘制效果:

    在这里插入图片描述

    5 修改多边形格式

    使用setPatch批量修改多边形格式,Patch对象具有的格式均可以修改。

    使用setPatchN修改第N个多边形格式,Patch对象具有的格式均可以修改。

    批量修改格式

    全部改为灰黑色:

    VN.setPatch('FaceColor',[0,0,0],'EdgeColor',[0,0,0])
    
    • 1

    修改第N个格式

    单独将第1个对象改为灰黑色:

    VN.setPatchN(1,'FaceColor',[0,0,0],'EdgeColor',[0,0,0])
    
    • 1

    循环修改格式

    循环修改各个多边形颜色:

    colorList=[78 101 155;
              138 140 191;
              184 168 207;
              231 188 198;
              253 207 158;
              239 164 132;
              182 118 108]./255;
    for i=1:7
        VN.setPatchN(i,'FaceColor',colorList(i,:),'EdgeColor',colorList(i,:))
    end
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9
    • 10

    6 字体修改

    使用:

    • setFont
    • setLabel

    分别修改数字和标签的格式,举个例子:

    VN.setFont('Color',[.9,0,0],'FontSize',14)
    VN.setLabel('Color',[0,0,.9],'FontSize',25,'FontName','Cambria')
    
    • 1
    • 2


    完整代码

    注意:仅代码无法运行,需要导入.mat格式的素材包,素材包及m文件获取放在文末,

    classdef venn
    % @author : slandarer
    % 公众号  : slandarer随笔
    
    % 使用示例:
    % =========================================================================
    % X1=randi([1,500],[100,1]);
    % X2=randi([1,500],[100,1]);
    % X3=randi([1,500],[100,1]);
    % X4=randi([1,500],[100,1]);
    % X5=randi([1,500],[100,1]);
    % X6=randi([1,500],[100,1]);
    % X7=randi([1,500],[100,1]);
    % XX={X1,X2,X3,X4,X5,X6,X7};
    % 
    % 
    % VN=venn(XX{1:7});
    % VN.draw();
    
        properties
            ax % 绘图坐标区域
            % -----------------------------------------------------------------
            linePnts
            labelSet={' ',' ',' ',' ',' ',' ',' '};
    %         labels={'AAA','BBB','CCC','DDD','EEE','FFF','GGG'};
            labelPos
            % -----------------------------------------------------------------
            classNum    % 多边形数量
            dataList    % 数据列表
            pshapeHdl   % polyshape对象
            fillHdl     % fill绘制的半透明多边形
            textHdl     % 绘制文本的句柄
            labelHdl    % 绘制标签的句柄
        end
    
        methods
            function obj = venn(varargin)
                if isa(varargin{1},'matlab.graphics.axis.Axes')
                    obj.ax=varargin{1};varargin(1)=[];
                else
                    obj.ax=gca;
                end
                hold on
                obj.classNum=length(varargin);
                obj.dataList=varargin;
    
                % 获取mat文件中的线条数据
                obj.linePnts=load('LD.mat');
                obj.linePnts=obj.linePnts.lineData;
    
                % 初始定义标签位置
                obj.labelPos{2}=[-.38,.3;.38,.3];
                obj.labelPos{3}=[-.38,.3;-.38,-.4;.38,-.4];
                obj.labelPos{4}=[-.38,.2;.38,.2;-.15,.3;.15,.3];
                obj.labelPos{5}=[cos(linspace(2*pi/5,2*pi,5)+2*pi/5-pi/7).*.47;
                                 sin(linspace(2*pi/5,2*pi,5)+2*pi/5-pi/7).*.47]';
                obj.labelPos{6}=[cos(linspace(2*pi/6,2*pi,6)+2*pi/3-pi/6).*.49;
                                 sin(linspace(2*pi/6,2*pi,6)+2*pi/3-pi/6).*.49]';
                obj.labelPos{6}=obj.labelPos{6}+[0,+.09;-.01,-.04;0,+.015;0,-.1;0,0;0,-.015];
                obj.labelPos{7}=[cos(linspace(2*pi/7,2*pi,7)+2*pi/5-pi/7).*.47;
                                 sin(linspace(2*pi/7,2*pi,7)+2*pi/5-pi/7).*.47]';
                help venn
            end
    
            function obj=draw(obj)
                warning off
                % 坐标区域修饰
                obj.ax.XLim=[-.5,.5];
                obj.ax.YLim=[-.5,.5];
                obj.ax.XTick=[];
                obj.ax.YTick=[];
                obj.ax.XColor='none';
                obj.ax.YColor='none';
                obj.ax.PlotBoxAspectRatio=[1,1,1];
                % 循环绘制半透明多边形
                tcolorList=lines(7);
                for i=1:obj.classNum
                    tPData=obj.linePnts(obj.classNum).pnts{i};
                    obj.pshapeHdl{i}=polyshape(tPData(:,1),tPData(:,2));
                    obj.fillHdl(i)=fill(tPData(:,1),tPData(:,2),tcolorList(i,:),...
                        'FaceAlpha',.2,'LineWidth',1.5,'EdgeColor',tcolorList(i,:));
                end
                % 构造初始bool集合
                baseData=[];
                for i=1:obj.classNum
                    baseData=[baseData;obj.dataList{i}(:)];
                end
                baseShpae=polyshape([-.5,-.5,.5,.5],[.5,-.5,-.5,.5]);
                pBool=abs(dec2bin((1:(2^obj.classNum-1))'))-48;
                % 循环绘制标签
                for i=1:obj.classNum
                    tPos=obj.labelPos{obj.classNum};
                    obj.labelHdl(i)=text(tPos(i,1),tPos(i,2),obj.labelSet{i},...
                        'HorizontalAlignment','center','FontName','Arial','FontSize',16);
                end
                % 循环计算数字位置
                for i=1:size(pBool,1)
                    tShpae=baseShpae;
                    tData=baseData;
                    for j=1:size(pBool,2)
                        switch pBool(i,j)
                            case 1
                                tShpae=intersect(tShpae,obj.pshapeHdl{j});
                                tData=intersect(tData,obj.dataList{j});
                            case 0
                                tShpae=subtract(tShpae,obj.pshapeHdl{j});
                                tData=setdiff(tData,obj.dataList{j});
                        end                 
                    end
                    [cx,cy]=centroid(tShpae);
                    obj.textHdl(i)=text(cx,cy,num2str(length(tData)),...
                        'HorizontalAlignment','center','FontName','Arial');
                end  
            end
            % =================================================================
            % 设置标签文本内容
            function obj=labels(obj,varargin)
                tlabel{length(varargin)}=' ';            
                for i=1:length(varargin)
                    tlabel{i}=varargin{i};
                end
                obj.labelSet=tlabel;
            end
            % 批量设置多边形格式
            function setPatch(obj,varargin)
                for i=1:obj.classNum
                    set(obj.fillHdl(i),varargin{:})
                end
            end
            % 单独设置多边形格式
            function setPatchN(obj,N,varargin)
                for i=1:obj.classNum
                    set(obj.fillHdl(N),varargin{:})
                end
            end
            % 设置数值字体
            function setFont(obj,varargin)
                for i=1:length(obj.textHdl)
                    set(obj.textHdl(i),varargin{:})
                end
            end
            % 设置标签字体
            function setLabel(obj,varargin)
                for i=1:length(obj.labelHdl)
                    set(obj.labelHdl(i),varargin{:})
                end
            end
        end
    % @author : slandarer
    % 公众号  : slandarer随笔
    end
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9
    • 10
    • 11
    • 12
    • 13
    • 14
    • 15
    • 16
    • 17
    • 18
    • 19
    • 20
    • 21
    • 22
    • 23
    • 24
    • 25
    • 26
    • 27
    • 28
    • 29
    • 30
    • 31
    • 32
    • 33
    • 34
    • 35
    • 36
    • 37
    • 38
    • 39
    • 40
    • 41
    • 42
    • 43
    • 44
    • 45
    • 46
    • 47
    • 48
    • 49
    • 50
    • 51
    • 52
    • 53
    • 54
    • 55
    • 56
    • 57
    • 58
    • 59
    • 60
    • 61
    • 62
    • 63
    • 64
    • 65
    • 66
    • 67
    • 68
    • 69
    • 70
    • 71
    • 72
    • 73
    • 74
    • 75
    • 76
    • 77
    • 78
    • 79
    • 80
    • 81
    • 82
    • 83
    • 84
    • 85
    • 86
    • 87
    • 88
    • 89
    • 90
    • 91
    • 92
    • 93
    • 94
    • 95
    • 96
    • 97
    • 98
    • 99
    • 100
    • 101
    • 102
    • 103
    • 104
    • 105
    • 106
    • 107
    • 108
    • 109
    • 110
    • 111
    • 112
    • 113
    • 114
    • 115
    • 116
    • 117
    • 118
    • 119
    • 120
    • 121
    • 122
    • 123
    • 124
    • 125
    • 126
    • 127
    • 128
    • 129
    • 130
    • 131
    • 132
    • 133
    • 134
    • 135
    • 136
    • 137
    • 138
    • 139
    • 140
    • 141
    • 142
    • 143
    • 144
    • 145
    • 146
    • 147
    • 148
    • 149
    • 150
    • 151

    MATLAB韦恩图绘制能画成这样属实不易,如果有用请留个叭~

    未经允许本代码请勿作商业用途,引用的话可以引用我file exchange上的链接,可使用如下格式:

    Zhaoxu Liu (2022). venn diagram 韦恩图 (https://www.mathworks.com/matlabcentral/fileexchange/116760-venn-diagram), MATLAB Central File Exchange. 检索来源 2022/8/26.

    若转载请保留以上file exchange链接及本文链接!!!

    全部m文件及素材包获取:
    链接:https://pan.baidu.com/s/1w1YfaXN6d5j0Ras60MVQ2Q?pwd=slan
    提取码:slan

  • 相关阅读:
    实用调试小技巧
    2022牛客暑期多校训练营1 个人题解
    Redis系列之常见数据类型应用场景
    带你深入理解3.4.2的版本更新,对用户带来了什么?
    python计算当天零点时间
    Nuxt3框架局部文件引用外部JS/CSS文件的相关配置方法
    【已验证-直接用】微信小程序wx.request请求服务器json数据并渲染到页面
    【Nuxt3】Vue3 + Element-plus 打包后报错 @popperjs/core
    【javaSE】 Lambda表达式与Lambda表达式的使用
    手把手教你搭建农产品商城小程序:详细步骤解析
  • 原文地址:https://blog.csdn.net/slandarer/article/details/126544953