MATLAB | 全网唯一,双变量及三变量映射图表的MATLAB绘制

最新推荐文章于 2024-07-29 14:25:51 发布
最新推荐文章于 2024-07-29 14:25:51 发布
阅读量8.4k 收藏 59
点赞数 11
分类专栏： # MATLAB 科学绘图 文章标签： matlab 开发语言 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/slandarer/article/details/126702967
版权

应粉丝要求写了一个多变量映射图绘制的工具函数，依旧全网唯一hiahiahia，一般图像只有一个变量需要映射到颜色，其他变量映射到点的大小啥的，但这两个变量之间的联系就没有那么直观，而且如果变量更多例如三个变量就比较难找到表现形式，但是，多变量映射图可以解决这些问题。

本篇先讲解双变量映射图表绘制再讲解三变量映射图表绘制。两个工具函数由于过长均放在文末。

示例图

双变量映射图表

输入输出说明

需要A,B两个完全相同大小的数值矩阵，通过

[CMapData,CMapHdl]=multiVarMap2(A,B)

的格式调用双变量映射图绘制工具函数，返回值为一个CMapDataRGB颜色数组，用于给图形对象赋予颜色，以及一个CMapHdl图例。

基本使用

这里第一个例子是一个网格数据,multi2_3.mat中存储着变量X1,X2,A,B,C这里只用AB数据，分别是当前格点属于类A和B的概率，简单编写函数如下：

multiData=load('multi2_3.mat');
X1=multiData.X1;
X2=multiData.X2;
A=multiData.A;
B=multiData.B;

% 先绘制一个没有赋予颜色的图形
surfHdl=surf(X1,X2,zeros(size(X1)),'EdgeColor','none');
hold on;view(2);axis tight

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap2(A,B); 
set(surfHdl,'CData',CMapData);

颜色划分

上面每个变量颜色划分了四段，我们可以调整pieceNum属性划分更多段颜色，例如分成8段：

multiData=load('multi2_3.mat');
X1=multiData.X1;
X2=multiData.X2;
A=multiData.A;
B=multiData.B;

% 先绘制一个没有赋予颜色的图形
surfHdl=surf(X1,X2,zeros(size(X1)),'EdgeColor','none');
hold on;view(2);axis tight

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap2(A,B,'pieceNum',8); 
set(surfHdl,'CData',CMapData);

其他配色

们可以调整colorList属性调整配色，需要4x3大小的RGB数组，分别代表左下，右下，左上，右上的颜色，当然函数内置了12种默认配色，通过1-12的编号调用：

multiData=load('multi2_3.mat');
X1=multiData.X1;
X2=multiData.X2;
A=multiData.A;
B=multiData.B;

% 先绘制一个没有赋予颜色的图形
surfHdl=surf(X1,X2,zeros(size(X1)),'EdgeColor','none');
hold on;view(2);axis tight

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap2(A,B,'pieceNum',8,'colorList',10); 
set(surfHdl,'CData',CMapData);

全部配色展示：

图例移动

可以随便拖拽：

修改变量名和标题

multiData=load('multi2_3.mat');
X1=multiData.X1;
X2=multiData.X2;
A=multiData.A;
B=multiData.B;

% 先绘制一个没有赋予颜色的图形
surfHdl=surf(X1,X2,zeros(size(X1)),'EdgeColor','none');
hold on;view(2);axis tight

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap2(A,B,'pieceNum',8,'colorList',10,'title','XXXXX','varName',{'AAAAA','BBBBB'}); 
set(surfHdl,'CData',CMapData);

三维绘图

不管二维还是三维图例都可以添加

multiData=load('multi2_3.mat');
X1=multiData.X1;
X2=multiData.X2;
A=multiData.A;
B=multiData.B;

% 先绘制一个没有赋予颜色的图形
surfHdl=surf(X1,X2,max(A,B),'EdgeColor','none');
hold on;view(3);axis tight

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap2(A,B,'colorList',10,'pieceNum',8);
set(surfHdl,'CData',CMapData);

示例2地图绘制

本文展示mapping toolbox绘制的图双变量映射，因此需要提前安装mapping toolbox工具箱，当然其他地图绘制工具箱绘制的图用本文的双变量映射也差不多，以下内容懒得收集数据，于是随便生成了一些数据：

axm=usamap("conus");

states=shaperead("usastatelo.shp",'UseGeoCoords',true);
% 俩州离太远画不开，不要
for i=length(states):-1:1
    if states(i).Name=="Alaska"||states(i).Name=="Hawaii"
        states(i)=[];
    end
end

faceColors=makesymbolspec('Polygon',{'INDEX',[1 numel(states)],'FaceColor',...
  polcmap(numel(states))});
usaMap=geoshow(states, 'DisplayType','polygon','SymbolSpec', faceColors);

% 不想去收集数据了，随便生成点随机数
A=rand(numel(states),1);
B=rand(numel(states),1);

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap2(A,B,'colorList',11,'pieceNum',8);
% 循环修改每个州颜色
for i=1:length(usaMap.Children)
    usaMap.Children(i).FaceColor=[CMapData(i,1,1),CMapData(i,1,2),CMapData(i,1,3)];
end

示例3气泡图绘制

x=1:20;
y=rand(1,20);
sz=rand(1,20);
BC=bubblechart(x,y,sz); 

% 不想去收集数据了，随便生成点随机数
A=rand(1,20);
B=rand(1,20);


% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap2(A,B,'colorList',10,'pieceNum',8);
BC.CData=[CMapData(:,:,1)',CMapData(:,:,2)',CMapData(:,:,3)'];

三变量映射图表

使用方式和双变量映射区别不大：

multiData=load('multi2_3.mat');
X1=multiData.X1;
X2=multiData.X2;
A=multiData.A;
B=multiData.B;
C=multiData.C;

% 先绘制一个没有赋予颜色的图形
surfHdl=surf(X1,X2,zeros(size(X1)),'EdgeColor','none');
hold on;view(2);axis tight

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap3(A,B,C,'colorList',1,'pieceNum',8);
set(surfHdl,'CData',CMapData);

其他配色

只准备了俩配色，以下内容因为画完后背景色和图例文字色相同因此修改了图例文本颜色：

multiData=load('multi2_3.mat');
X1=multiData.X1;
X2=multiData.X2;
A=multiData.A;
B=multiData.B;
C=multiData.C;

% 先绘制一个没有赋予颜色的图形
surfHdl=surf(X1,X2,zeros(size(X1)),'EdgeColor','none');
hold on;view(2);axis tight

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap3(A,B,C,'colorList',2,'pieceNum',8);
set(surfHdl,'CData',CMapData);
CMapHdl.XColor=[1,1,1];
CMapHdl.YColor=[1,1,1];
CMapHdl.ZColor=[1,1,1];
CMapHdl.Title.Color=[1,1,1];

图例移动

可以随便拖拽、旋转：

示例2地图绘制

随机生成颜色太乱，因此第二个图把所有数据都从大到小排列，生成可一组不太乱的随机数据。

axm=usamap("conus");

states=shaperead("usastatelo.shp",'UseGeoCoords',true);
% 俩州离太远画不开，不要
for i=length(states):-1:1
    if states(i).Name=="Alaska"||states(i).Name=="Hawaii"
        states(i)=[];
    end
end

faceColors=makesymbolspec('Polygon',{'INDEX',[1 numel(states)],'FaceColor',...
  polcmap(numel(states))});
usaMap=geoshow(states, 'DisplayType','polygon','SymbolSpec', faceColors);

% 不想去收集数据了，随便生成点随机数
A=sort(rand(numel(states),1));
B=sort(rand(numel(states),1));
C=sort(rand(numel(states),1));

% A=rand(numel(states),1);
% B=rand(numel(states),1);
% C=rand(numel(states),1);

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap3(A,B,C,'colorList',1,'pieceNum',8);
% 循环修改每个州颜色
for i=1:length(usaMap.Children)
    usaMap.Children(i).FaceColor=[CMapData(i,1,1),CMapData(i,1,2),CMapData(i,1,3)];
end

示例3气泡图绘制

x=1:20;
y=rand(1,20);
sz=rand(1,20);
BC=bubblechart(x,y,sz); 

% 不想去收集数据了，随便生成点随机数
A=rand(1,20);
B=rand(1,20);
C=rand(1,20);

% 调用工具函数生成图例和映射表
[CMapData,CMapHdl]=multiVarMap3(A,B,C,'colorList',1,'pieceNum',8);
BC.CData=[CMapData(:,:,1)',CMapData(:,:,2)',CMapData(:,:,3)'];

工具函数完整代码

双变量

function [CMapData,CMapHdl]=multiVarMap2(A,B,varargin) 
% @author : slandarer
% gzh  : slandarer随笔
help multiVarMap2
% colorList应为4x3大小RGB数组
% [R1,G1,B1;R2,G2,B2;R3,B3,G3;R4,G4,B4]
% 分别对应图例左下，右下，左上，右上的颜色
% 或数值编号1-12

% =========================================================================
% 基本属性提取
obj.arginList={'colorList','pieceNum','varName','title'};
obj.colorList=1;
obj.pieceNum=4;
obj.varName={'var1','var2'};
obj.title='ColorMap';
for i=1:2:(length(varargin)-1)
    tid=ismember(obj.arginList,varargin{i});
    if any(tid)
        obj.(obj.arginList{tid})=varargin{i+1};
    end
end
% =========================================================================
% 基础配色表
baseAtla{1}=[233,233,233;201, 90, 90;100,173,191; 86, 64, 71]./255;
baseAtla{2}=[249,244,248;139,182,224;235,130,144;122,123,171]./255;
baseAtla{3}=[243,243,243;140,227,176;231,164,209;124,143,176]./255;
baseAtla{4}=[233,233,233;201,180, 90;154,115,176;129, 76, 51]./255;
baseAtla{5}=[222,222,222;  0,110,175;205,  0, 31; 74, 33, 76]./255;
baseAtla{6}=[233,233,233;108,132,182;116,175,129; 37, 90, 91]./255;
baseAtla{7}=[233,233,233; 89,201,201;191,100,173; 56, 71,149]./255;
baseAtla{8}=[243,243,243; 79,158,195;243,180,  0; 40, 40, 40]./255;
baseAtla{9}=[233,231,242; 78,174,209;223, 78,167; 37,19,139]./255;
baseAtla{10}=[227,228,223;127,135,124; 57,120,164; 20, 72, 83]./255;
baseAtla{11}=[243,232,156;248,161,127;207,102,148; 92, 82,166]./255;
baseAtla{12}=[248,225,152;181,167,124;108,119,149; 65, 92,170]./255;

if numel(obj.colorList)==1,atla=baseAtla{obj.colorList};end
% =========================================================================
% 数据处理及图例色卡数据构建
% nX1=X1(:);nX2=X2(:);
CMapData=getCMap(A,B,atla,obj.pieceNum);
    function CMap=getCMap(A,B,atla,n)
        nA=(A(:)-min(A(:)))./(max(A(:))-min(A(:)));
        nB=(B(:)-min(B(:)))./(max(B(:))-min(B(:)));
        nA=(ceil(nA.*n)-1)./(n-1);nB=(ceil(nB.*n)-1)./(n-1);
        nA(nA<0)=0;nB(nB<0)=0;

        colorAB=nA.*nB.*atla(4,:)+...
                nA.*(1-nB).*atla(2,:)+...
                (1-nA).*nB.*atla(3,:)+...
                (1-nA).*(1-nB).*atla(1,:);

        CMap=zeros([size(A),3]);
        CMap(:,:,1)=reshape(colorAB(:,1),size(A));
        CMap(:,:,2)=reshape(colorAB(:,2),size(A));
        CMap(:,:,3)=reshape(colorAB(:,3),size(A));
    end
% =========================================================================
% 坐标区域初定位及修饰
ax1=gca;fig=ax1.Parent;
ax1Pos=ax1.Position;
ax2=axes('Parent',fig,'Position',[ax1Pos(1)+3/4*ax1Pos(3),...
    ax1Pos(2)+6/9*ax1Pos(4),ax1Pos(3)/4,ax1Pos(4)*2.5/9]);
ax1.UserData=1;
ax2.UserData=2;
ax2.XLim=[min(A(:)),max(A(:))];
ax2.YLim=[min(B(:)),max(B(:))];
ax2.TickDir='out';ax2.TickLength=[0.02,0.025];
ax2.PlotBoxAspectRatio=[1,1,1];
ax2.Layer='top';
ax2.LineWidth=1;
ax2.FontName='Arial';
ax2.FontSize=9;
ax2.Title.String=obj.title;
ax2.Title.FontSize=12;
ax2.XLabel.String=obj.varName{1};
ax2.XLabel.FontSize=11;
ax2.YLabel.String=obj.varName{2};
ax2.YLabel.FontSize=11;
ax2.Toolbar.Visible='off';
hold(ax2,'on');view(2);
% =========================================================================
% 绘制右上角图例
[AA,BB]=meshgrid(linspace(0,1,obj.pieceNum),...
                 linspace(0,1,obj.pieceNum));
tCMap=getCMap(AA,BB,atla,obj.pieceNum);
legendCMap=zeros([size(AA)+1,3]);
legendCMap(1:end-1,1:end-1,1)=tCMap(:,:,1);
legendCMap(1:end-1,1:end-1,2)=tCMap(:,:,2);
legendCMap(1:end-1,1:end-1,3)=tCMap(:,:,3);
[XX,YY]=meshgrid(linspace(ax2.XLim(1),ax2.XLim(2),obj.pieceNum+1),...
                 linspace(ax2.YLim(1),ax2.YLim(2),obj.pieceNum+1));
surf(ax2,XX,YY,zeros(size(XX)),'CData',legendCMap,'EdgeColor','none','FaceColor','flat');
% =========================================================================
% 图例拖动回调
set(fig,'WindowButtonDownFcn',@BtnDownFunc, ...
        'WindowButtonUpFcn',@BtnUpFunc, ...
        'WindowButtonMotionFcn',@BtnMotionFunc)
selectedAx=false;
    % 选择图例时的回调
    function BtnDownFunc(~,~)
        taxes=get(fig,'currentAxes');
        if taxes.UserData==2
            selectedAx=true;
        end
    end
    % 松开图例时的回调
    function BtnUpFunc(~,~)
        selectedAx=false;
        set(fig,'currentAxes',ax1)
    end
    % 选择图例并拖动时的回调
    function BtnMotionFunc(~,~)
        if selectedAx
            tpos=get(fig,'CurrentPoint');
            figPos=fig.Position;
            normPos=tpos./figPos(3:4);
            ax2.Position(1:2)=-ax2.Position(3:4)./2+normPos;
        end
    end
set(fig,'currentAxes',ax1)
CMapHdl=ax2;
end

三变量

function [CMapData,CMapHdl]=multiVarMap3(A,B,C,varargin) 
% @author : slandarer
% gzh  : slandarer随笔

help multiVarMap3
% colorList应为8x3大小RGB数组
% [R1,G1,B1;R2,G2,B2;R3,B3,G3;R4,G4,B4,...]

% =========================================================================
% 基本属性提取
obj.arginList={'colorList','pieceNum','varName','title'};
obj.colorList=1;
obj.pieceNum=6;
obj.varName={'var1','var2','var3'};
obj.title='ColorMap';
for i=1:2:(length(varargin)-1)
    tid=ismember(obj.arginList,varargin{i});
    if any(tid)
        obj.(obj.arginList{tid})=varargin{i+1};
    end
end
% =========================================================================
% 基础配色表
% [0;X;Y;Z;XY;YZ;XZ;XYZ]
baseAtla{1}=[222,222,222;201,90,90;116,175,129;108,132,182;
             243,180,0;37,90,91;74,33,76;40 40 40]./255; 
baseAtla{2}=[0 0 0;1 0 0;0 1 0;0 0 1;1 1 0;0 1 1;1 0 1;1 1 1]; 

if numel(obj.colorList)==1,atla=baseAtla{obj.colorList};end
% =========================================================================
% 数据处理及图例色卡数据构建
% nX1=X1(:);nX2=X2(:);
CMapData=getCMap3(A,B,C,atla,obj.pieceNum);
    function CMap=getCMap2(A,B,atla,n)
        nA=(A(:)-min(A(:)))./(max(A(:))-min(A(:)));
        nB=(B(:)-min(B(:)))./(max(B(:))-min(B(:)));
        nA=(ceil(nA.*n)-1)./(n-1);
        nB=(ceil(nB.*n)-1)./(n-1);
        nA(nA<0)=0;nB(nB<0)=0;

        colorAB=nA.*nB.*atla(4,:)+...
                nA.*(1-nB).*atla(2,:)+...
                (1-nA).*nB.*atla(3,:)+...
                (1-nA).*(1-nB).*atla(1,:);

        CMap=zeros([size(A),3]);
        CMap(:,:,1)=reshape(colorAB(:,1),size(A));
        CMap(:,:,2)=reshape(colorAB(:,2),size(A));
        CMap(:,:,3)=reshape(colorAB(:,3),size(A));
    end
    function CMap=getCMap3(A,B,C,atla,n)
        nA=(A(:)-min(A(:)))./(max(A(:))-min(A(:)));
        nB=(B(:)-min(B(:)))./(max(B(:))-min(B(:)));
        nC=(C(:)-min(C(:)))./(max(C(:))-min(C(:)));
        nA=(ceil(nA.*n)-1)./(n-1);
        nB=(ceil(nB.*n)-1)./(n-1);
        nC=(ceil(nC.*n)-1)./(n-1);
        nA(nA<0)=0;nB(nB<0)=0;nC(nC<0)=0;

        colorAB=atla(1,:).*(1-nA).*(1-nB).*(1-nC)+...
                atla(2,:).*(nA).*(1-nB).*(1-nC)+...
                atla(3,:).*(1-nA).*(nB).*(1-nC)+...
                atla(4,:).*(1-nA).*(1-nB).*(nC)+...
                atla(5,:).*(nA).*(nB).*(1-nC)+...
                atla(6,:).*(1-nA).*(nB).*(nC)+...
                atla(7,:).*(nA).*(1-nB).*(nC)+...
                atla(7,:).*(nA).*(nB).*(nC);

        CMap=zeros([size(A),3]);
        CMap(:,:,1)=reshape(colorAB(:,1),size(A));
        CMap(:,:,2)=reshape(colorAB(:,2),size(A));
        CMap(:,:,3)=reshape(colorAB(:,3),size(A));
    end
% =========================================================================
% 坐标区域初定位及修饰
ax1=gca;fig=ax1.Parent;
ax1Pos=ax1.Position;
ax2=axes('Parent',fig,'Position',[ax1Pos(1)+6/9*ax1Pos(3),...
    ax1Pos(2)+6/9*ax1Pos(4),ax1Pos(3)*2.5/9,ax1Pos(4)*2.5/9]);
ax1.UserData=1;
ax2.UserData=2;
ax2.XLim=[min(A(:)),max(A(:))];
ax2.YLim=[min(B(:)),max(B(:))];
ax2.ZLim=[min(C(:)),max(C(:))];
ax2.TickDir='out';ax2.TickLength=[0.02,0.025];
ax2.PlotBoxAspectRatio=[1,1,1];
ax2.Layer='top';
ax2.LineWidth=1;
ax2.FontName='Arial';
ax2.FontSize=9;
ax2.Title.String=obj.title;
ax2.Title.FontSize=12;
ax2.XLabel.String=obj.varName{1};
ax2.XLabel.FontSize=11;
ax2.YLabel.String=obj.varName{2};
ax2.YLabel.FontSize=11;
ax2.ZLabel.String=obj.varName{3};
ax2.ZLabel.FontSize=11;
ax2.Toolbar.Visible='off';
hold(ax2,'on');view(3);
% =========================================================================
% 绘制右上角图例
[AA,BB]=meshgrid(linspace(0,1,obj.pieceNum),...
                 linspace(0,1,obj.pieceNum));

    function newMap=expMap(oriMap)
        newMap=zeros([size(oriMap,1)+1,size(oriMap,2)+1,3]);
        newMap(1:end-1,1:end-1,1)=oriMap(:,:,1);
        newMap(1:end-1,1:end-1,2)=oriMap(:,:,2);
        newMap(1:end-1,1:end-1,3)=oriMap(:,:,3);
    end

[XX,YY]=meshgrid(linspace(ax2.XLim(1),ax2.XLim(2),obj.pieceNum+1),linspace(ax2.YLim(1),ax2.YLim(2),obj.pieceNum+1));
tCMap=getCMap2(AA,BB,atla([1,2,3,5],:),obj.pieceNum);tCMap=expMap(tCMap);
surf(ax2,XX,YY,ones(size(XX)).*ax2.ZLim(1),'CData',tCMap,'EdgeColor','none','FaceColor','flat');
tCMap=getCMap2(AA,BB,atla([4,7,6,8],:),obj.pieceNum);tCMap=expMap(tCMap);
surf(ax2,XX,YY,ones(size(XX)).*ax2.ZLim(2),'CData',tCMap,'EdgeColor','none','FaceColor','flat');
% -------------------------------------------------------------------------
[XX,ZZ]=meshgrid(linspace(ax2.XLim(1),ax2.XLim(2),obj.pieceNum+1),linspace(ax2.ZLim(1),ax2.ZLim(2),obj.pieceNum+1));
tCMap=getCMap2(AA,BB,atla([1,2,4,7],:),obj.pieceNum);tCMap=expMap(tCMap);
surf(ax2,XX,ones(size(XX)).*ax2.YLim(1),ZZ,'CData',tCMap,'EdgeColor','none','FaceColor','flat');
tCMap=getCMap2(AA,BB,atla([3,5,6,8],:),obj.pieceNum);tCMap=expMap(tCMap);
surf(ax2,XX,ones(size(XX)).*ax2.YLim(2),ZZ,'CData',tCMap,'EdgeColor','none','FaceColor','flat');
% -------------------------------------------------------------------------
[YY,ZZ]=meshgrid(linspace(ax2.YLim(1),ax2.YLim(2),obj.pieceNum+1),linspace(ax2.ZLim(1),ax2.ZLim(2),obj.pieceNum+1));
tCMap=getCMap2(AA,BB,atla([1,3,4,6],:),obj.pieceNum);tCMap=expMap(tCMap);
surf(ax2,ones(size(YY)).*ax2.XLim(1),YY,ZZ,'CData',tCMap,'EdgeColor','none','FaceColor','flat');
tCMap=getCMap2(AA,BB,atla([2,5,7,8],:),obj.pieceNum);tCMap=expMap(tCMap);
surf(ax2,ones(size(YY)).*ax2.XLim(2),YY,ZZ,'CData',tCMap,'EdgeColor','none','FaceColor','flat');

% 图例拖动回调
set(fig,'WindowButtonDownFcn',@BtnDownFunc, ...
        'WindowButtonUpFcn',@BtnUpFunc, ...
        'WindowButtonMotionFcn',@BtnMotionFunc, ...
        'KeyPressFcn',@keyFunc)
selectedAx=false;
    % 选择图例时的回调
    function BtnDownFunc(~,~)
        taxes=get(fig,'currentAxes');
        if taxes.UserData==2
            selectedAx=true;
        end
    end
    % 松开图例时的回调
    function BtnUpFunc(~,~)
        selectedAx=false;
        set(fig,'currentAxes',ax1)
    end
    % 选择图例并拖动时的回调
    function BtnMotionFunc(~,~)
        if selectedAx
            tpos=get(fig,'CurrentPoint');
            figPos=fig.Position;
            normPos=tpos./figPos(3:4);
            ax2.Position(1:2)=-ax2.Position(3:4)./2+normPos;
        end
    end
set(fig,'currentAxes',ax1)
CMapHdl=ax2;
end

完

如果有用请留个赞叭~

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

Zhaoxu Liu (2022). bivariate/ternary choropleth map,多变量映射图表 (https://www.mathworks.com/matlabcentral/fileexchange/117165-bivariate-ternary-choropleth-map), MATLAB Central File Exchange. 检索来源 2022/9/5.

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

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

slandarer
关注
  • 11
    点赞
  • 踩
  • 59
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
专栏目录
MATLAB｜科研绘图｜山脊图
猪猪侠
11-14 617
也被称为Joy Plot，是一种用于可视化数据分布的图表，特别是用于显示多个组的分布情况。在这种图表中，每个组的数据分布都通过平滑的密度曲线来表示，这些曲线沿着垂直轴堆叠，形成类似山脊的视觉效果。山脊图是核密度估计（KDE）的一个应用，它提供了比传统的条形图或直方图更平滑的数据分布视图。：山脊图主要用于展示和比较不同组或类别内数据的分布情况。
MATLAB画四维图，四维曲面图
06-05
很多问题会设计到三个自变量，一个因变量，一般做图方法只有三个坐标轴，MATLAB可以用颜色作为第四维进行画图，画四维图，X轴，Y轴， Z轴分别代表三个自变量，颜色的变化代表因变量，可以很直观的体现三个自变量和一个因变量的对应关系！
R语言绘制三变量分区地图
qq_32649321的博客
06-13 514
其中，shape = 21表示将点的形状设置为圆环，fill = "transparent"表示圆环内部透明填充，color = "red"表示圆环的颜色为红色，stroke = 2表示圆环的线宽为2。aes是"aesthetic"的缩写，用于ggplot2包中的函数，用来映射数据到图形属性，例如颜色、形状、大小等。学习笔记：Ternary choropleth maps。还可以设置不同行之间的比例：（如第三行是第一行的3倍）（2）加载点矢量数据，并显示在地图上。（3）根据某变量大小显示点的大小。
matlab基础--二维图和三维图
weixin_44605190的博客
01-17 1678
##线图 使用plot函数创建二位线图 例如： 绘制从 0 到 2π 的正弦函数值，命令： x=0:pi/100:2*pi; y=sin(x); plot(x,y) 运行结果： 同时，可以对图像做标记，例如，标记轴并添加标题，命令： xlabel('x') ylabel('sin(x)') title('plot of the function') 运行结果： 此外，向plot函数添加第三个输入参数，可用红色虚线绘制相同的变量（此时，定义的标题和标签不再被用于当前的图窗窗...
MATLAB科研绘图与学术图表绘制从入门到精通
热门推荐
努力是为了站在万人之中，成为别人的光
04-02 7万+
通过以上介绍，相信您已经对MATLAB科研绘图与学术图表绘制有了一定的了解。实践是最好的老师，希望您能通过不断尝试和实践，掌握MATLAB绘图技巧，为您的科研工作增添光彩。
MATLAB科技绘图与数据分析
weixin_43178406的博客
05-14 2万+
本文主要介绍了MATLAB科技绘图与数据分析，希望能对学习和使用MATLAB科技绘图和数据分析的同学们有所帮助。 文章目录 1. 前言 2. 书籍推荐 2.1 内容简介 2.2 本书作者 2.3 本书目录 2.4 适合读者 3. 购买链接
【深度学习-第4篇】使用MATLAB快速实现CNN多变量回归预测
fengzhuqiaoqiu的博客
10-20 2134
上一篇我们讲了使用CNN进行分类的MATLAB代码。 这一篇我们讲CNN的多变量回归预测。 是的，同样是傻瓜式的快速实现。
OFDM调制及OTFS调制（matlab仿真）
weixin_53750855的博客
11-23 5857
正交频分复用技术（OFDM）是多载波调制技术之一，基本原理就是将传输信道分成若干个正交的子信道，将高速数据流转换成并行的低速子数据流，调制到每个子信道上进行传输。它可以利用快速傅里叶反变换（IFFT）和快速傅里叶变换（FFT）来完成数据的调制和解调。由于数据在子信道上传输的速率较低，在经过无线频率选择性衰落信道时，每个子信道都可以看成是平坦性衰落，只要进行适当的频率均衡就能够克服由衰落产生的影响。
matlab四维图形,MATLAB二维、三维、四维绘图
weixin_35275073的博客
03-16 4291
二维图形绘制直角坐标极坐标统计图面域图填色图羽毛图、射线图Voronoi图彩带图二维半图三维图形绘制线状图数组图网格图面状图统计图瀑布图散点图函数图四维图形绘制切片图颜色权重图Reference二维作图一般使用plot, polar函数，辅以axis, xlabel, ylabel, legend, legend, title等添加坐标轴标题、图表标题等信息。也可以打开数据表格，选中需要绘图的数据...
Matlab GUI编程技巧（九）：详解 uitable 函数显示表格及相关操作（创建表用户界面）
左手の明天的博客
05-25 1万+
🔗运行环境：Matlab 🚩撰写作者：左手の明天 🥇精选专栏：《python》 🔥推荐专栏：《算法研究》 💗 大家好🤗🤗🤗，我是左手の明天！💗 📆 最近更新：2022 年 5 月 25 日，左手の明天的第246篇原创博客 📚更新于专栏：Matlab GUI编程技巧 目录 🚩创建表用户界面组件：uitable ⭐️语法 ⭐️ColumnWidth-表列的宽度 ⭐️Position-表的位置和大小 ⭐️Units-测量单位 ⭐️显示表数组数据 ...
bubbleplot:绘制 3 或 4 个变量 - 第三个变量是符号大小，第四个（通过保持）是符号颜色。-matlab开发
06-01
在 MATLAB 中，`bubbleplot` 是一种非常有用的可视化工具，尤其适用于展示三维或四维数据。这个函数允许你同时展示两个定量变量（X 和 Y）的关系，并通过第三个变量（符号大小）和第四个变量（符号颜色）来进一步...
Matlab各种图表绘图演示1
08-04
自 Matlab 2014a 和 2014b 版本开始，Matlab 引入了全新的绘图系统，显著提升了图表的颜色配置和输出图片的质量，使得数据可视化更加美观且专业。下面将详细介绍几种常见的 Matlab 图表类型及其使用方法。 1. **...
MATLAB图形绘制-三维密度图
12-12
这种图表在科学研究、工程分析和数据分析等领域都有广泛的应用，特别是在处理多变量数据时。 三维密度图的核心在于将四维数据（x、y、z和强度/颜色）转换成一个三维空间中的点云，每个点的颜色代表其对应的第四维度...
plotgantt:从 Matlab 结构绘制甘特图。-matlab开发
05-31
在 MATLAB 开发中，创建可视化图表是数据理解和分析的关键步骤之一。甘特图是一种流行的时间线图表，用于展示项目进度、任务安排或者时间序列数据。本文将深入探讨如何使用自定义函数 `plotgantt` 从 MATLAB 结构...
成对相关可视化：该函数可视化多个变量之间的成对相关。-matlab开发
06-01
4. 绘制单变量直方图：为每个变量绘制单独的直方图，用以展示其分布特性。 5. 布局与同步：调整子图布局，确保轴限制和色标一致，消除重复的标签。 6. 添加标题和图例：为了增加可读性，不要忘记添加合适的图例和...
matlab 三维核密度图_python数据分布型图表直方图系列核密度估计图
weixin_39855186的博客
11-22 3842
核密度估计图核密度估计图(kernel density plot)用于显示数据在X轴连续数据段内的分布状况；直方图的变种，使用平滑曲线来绘制数值水平，从而得出更平滑的分布；较直方图，核密度估计图，不受所使用分组数量的影响，更好地界定分布形状；核密度估计(kernel density estimation)在概率论中用来估计未知的密度函数，属于非参数检验方法之一；由Rosenblatt和E...
MATLAB R2017a 安装与破解 ​编辑
最新发布
Hello-FPGA
07-29 345
转载自。
MATLAB基础：函数与函数控制语句
Zjs15630471263的博客
07-26 1124
今天我们继续学习Matlab中函数相关知识。help 命令是最基本的查询方法，可查询所有目录、指定目录、命令、函数。我们直接点击帮助菜单即可查询所需的API函数。lookfor 关键字用于搜索相关的命令和函数。如，我们输入lookfor input，MATLAB就会显示所有关于input的API函数type M 文件名用于显示文件代码。如果想显示以.m扩展名结尾的就要加通配符 which M 文件名用于显示指定文件的路径。返回一个路径demo 用于演示matlab功能。如图，我们可以获取各种演示代码来解
matlab用slice绘制
07-28
slice函数是MATLAB中用于绘制三维数据切片的函数。它可以显示在空间体横截面处所测得的变量的值。使用slice函数，你可以定义切片的位置和方向，并将其与颜色映射关联以表示测量的幅值。以下是使用slice函数绘制切片的示例代码： ```matlab load fluidtemp.mat % 加载数据 xslice = \[5 9.9\]; % 定义切片的x位置范围 yslice = 3; % 定义切片的y位置 zslice = \[-3 0\]; % 定义切片的z位置范围 pl = slice(x, y, z, temp, xslice, yslice, zslice); % 绘制切片 for i = 1:length(pl) pl(i).FaceColor = 'interp'; % 设置切片的颜色 pl(i).EdgeColor = 'none'; % 去除切片的边界线 end xlabel('x'); ylabel('y'); zlabel('z'); % 设置坐标轴标签 cb = colorbar; % 创建并标记颜色条 cb.Label.String = 'Temperature (℃)'; % 设置颜色条标签 title('Temperature Distribution'); % 设置图表标题 view(-34, 24); % 调整视角 ``` 这段代码将根据给定的数据和切片位置绘制出温度分布的切片图。切片图中的颜色表示温度的幅值，切片位置可以根据需要进行调整。你可以根据自己的数据和需求修改代码中的参数来绘制不同的切片图。\[1\] 另外，如果你想绘制单独的切片或全空间的立体切片，可以参考以下示例代码： ```matlab % 单独切片 v = csvread('vdata.txt'); % 加载数据 v = reshape(v, 5, 27, 6); % 重新调整数据的维度 \[x, y, z\] = meshgrid(1:27, 1:5, 1:6); % 创建网格坐标 figure; subplot(2, 2, 1); slice(x, y, z, v, \[\], \[1\], \[\]); shading interp; set(gca, 'zdir', 'reverse'); axis equal; grid on; subplot(2, 2, 2); slice(x, y, z, v, \[\], \[2\], \[\]); shading interp; colormap('jet'); set(gca, 'zdir', 'reverse'); axis equal; grid on; subplot(2, 2, 3); slice(x, y, z, v, \[\], \[3\], \[\]); shading interp; set(gca, 'zdir', 'reverse'); axis equal; grid on; subplot(2, 2, 4); slice(x, y, z, v, \[\], \[4\], \[\]); shading interp; set(gca, 'zdir', 'reverse'); axis equal; grid on; % 全空间立体切片 figure; slice(x, y, z, v, \[1:2:27\], \[2 3 4\], \[2 3 4 5\]); shading interp; colorbar; colormap('jet'); set(gca, 'zdir', 'reverse'); axis equal; grid on; box on; ``` 这段代码将根据给定的数据绘制出单独切片和全空间立体切片。你可以根据自己的数据和需求修改代码中的参数来绘制不同的切片图。\[2\]\[3\] #### 引用[.reference_title] - *1* [MATLAB slice函数绘制三维n×m×k矩阵的立体图或剖面图 —可视化四维数据](https://blog.csdn.net/u010510134/article/details/128036889)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insert_down28v1,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* *3* [matlab slice](https://blog.csdn.net/weixin_30814329/article/details/98533320)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insert_down28v1,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

slandarer

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值