MATLAB | 两种上色方式的旭日图绘制

2023-10-27 00:01
文章标签 matlab 方式 绘制 两种 旭日

本文主要是介绍MATLAB | 两种上色方式的旭日图绘制，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

嘿，这次真的是好久不见了，好不容易才有点空写点文章，这段时间忙到后台回复都有点来不及看，很抱歉有一部分后台留言刚看到就已经超过时限没法回复了，不过根据大家的留言，需求主要集中在希望出一期旭日图的教程，今天它来啦~

旭日图要说简单也不简单，要说难也不算难，唯一难点估计就是将数据分类排好序，然后每一圈能对应起来呗，上色啥的也不算啥难点，因此这期就简单教大家如何把大框架画出来以及提供两种上色方式。

还是挺酷的，近期没啥空单独开发出一个支持更多上色模式的类，大家可以先根据本文提供的代码先自行DIY着。

教程部分

0 数据

这里懒得搜集数据，自己随机生成了一个table类型数据，前三列都是名称，最后一列是数值：

rng(1)
ULList = 'ABCD';
LLList = 'abcd';rowNum = 50;
idx1   = randi([1,length(ULList)],[1,rowNum]);
idx2   = randi([1,3],[rowNum,1]);
Name1  = ULList(idx1'*ones(1,4));
Name2  = [LLList(idx1'*ones(1,4)),num2str(idx2)];
Name3  = [LLList(idx1)',num2str(idx2),...char(45.*ones(rowNum,1)),char(randi([97,122],[rowNum,5]))];
Value  = rand([rowNum,1]);
Table  = table(Name1,Name2,Name3,Value);

大概长这样：

1 数据处理

利用 grp2idx 进行分类，利用 sortrows 将相同类归在一起

% 利用 grp2idx 进行分类，利用 sortrows 将相同类归在一起
NameList  = Table.Properties.VariableNames;
NameNum   = length(NameList)-1;
NameCell{NameNum} = ' ';
valueList = zeros(length(Table.(NameList{1})),NameNum);
for i = 1:NameNum-1tName = Table.(NameList{i});tUniq = unique(tName,'rows');NameCell{i} = tUniq;ind = grp2idx([tUniq;tName]);ind(1:length(tUniq)) = [];valueList(:,i) = ind;
end
valueList(:,end) = -Table.(NameList{end});
[VAL,IDX] = sortrows(valueList,1:size(valueList,2));
VAL(:,end) = -VAL(:,end);

2 属性预设

设置配色，字体，文本显示的阈值等信息：

% 此处可以设置配色
CList=[0.3882    0.5333    0.70591.0000    0.6824    0.20390.9373    0.4353    0.41570.5490    0.7608    0.79220.3333    0.6784    0.53730.7647    0.7373    0.24710.7333    0.4627    0.57650.7294    0.6275    0.58040.6627    0.7098    0.68240.4627    0.4627    0.4627];
% 在这可修改字体
FontProp = {'FontSize',12,'Color',[0,0,0]};
% 在这可设置比例低于多少的部分不显示文字
TextThreshold = 0.012;

3 绘图

径向渐变

% 开始绘图
figure('Units','normalized','Position',[.2,.1,.52,.72]);
ax = gca; hold on
ax.DataAspectRatio = [1,1,1];
ax.XColor = 'none';
ax.YColor = 'none';
tT = linspace(0,1,100);
for i = 1:size(VAL,2)-1tRateSum = 0;tNum = length(NameCell{i});for j = 1:tNumtRate = sum(VAL(VAL(:,i) == j,end))./sum(VAL(:,end));  tTheta = [tRateSum+tT.*tRate,tRateSum+tRate-tT.*tRate].*pi.*2;tR = [tT.*0+i,tT.*0+i+1];if i == 1fill(cos(tTheta).*tR,sin(tTheta).*tR,CList(j,:),'EdgeColor',[1,1,1],'LineWidth',1)elsetCN = VAL(find(VAL(:,i) == j,1),1);fill(cos(tTheta).*tR,sin(tTheta).*tR,...CList(tCN,:).*0.8^(i-1)+[1,1,1].*(1-0.8^(i-1)),'EdgeColor',[1,1,1],'LineWidth',1)endrotation = (tRateSum+tRate/2)*360;if tRate > TextThresholdif rotation>90&&rotation<270rotation=rotation+180;text(cos((tRateSum+tRate/2).*pi.*2)*i,sin((tRateSum+tRate/2).*pi.*2)*i,NameCell{i}(j,:)+" ",FontProp{:},...'Rotation',rotation,'HorizontalAlignment','right')elsetext(cos((tRateSum+tRate/2).*pi.*2)*i,sin((tRateSum+tRate/2).*pi.*2)*i," "+NameCell{i}(j,:),FontProp{:},...'Rotation',rotation)endendtRateSum = tRateSum+tRate;end
end
% 绘制最外圈饼状图
tRateSum = 0;
tNameCell = Table.(NameList{end-1});
for j = 1:size(VAL,1)tRate = VAL(j,end)./sum(VAL(:,end)); tTheta = [tRateSum+tT.*tRate,tRateSum+tRate-tT.*tRate].*pi.*2;tR = [tT.*0+size(VAL,2),tT.*0+size(VAL,2)+1];tCN = VAL(j,1);fill(cos(tTheta).*tR,sin(tTheta).*tR,CList(tCN,:).*0.8^(size(VAL,2)-1)+...[1,1,1].*(1-0.8^(size(VAL,2)-1)),'EdgeColor',[1,1,1],'LineWidth',1)rotation = (tRateSum+tRate/2)*360;if tRate > TextThresholdif rotation>90&&rotation<270rotation=rotation+180;text(cos((tRateSum+tRate/2).*pi.*2)*size(VAL,2),sin((tRateSum+tRate/2).*pi.*2)*size(VAL,2),tNameCell(IDX(j),:)+" ",FontProp{:},...'Rotation',rotation,'HorizontalAlignment','right')elsetext(cos((tRateSum+tRate/2).*pi.*2)*size(VAL,2),sin((tRateSum+tRate/2).*pi.*2)*size(VAL,2)," "+tNameCell(IDX(j),:),FontProp{:},...'Rotation',rotation)endendtRateSum = tRateSum+tRate;
end

切向渐变

% 开始绘图
figure('Units','normalized','Position',[.2,.1,.52,.72]);
ax = gca; hold on
ax.DataAspectRatio = [1,1,1];
ax.XColor = 'none';
ax.YColor = 'none';
tT = linspace(0,1,100);
LCList = CList(1:length(NameCell{1}),:);
for i = 1:size(VAL,2)-1tRateSum = 0;tNum = length(NameCell{i});NCList = zeros(tNum,3);for j = 1:tNumtRate = sum(VAL(VAL(:,i) == j,end))./sum(VAL(:,end));  tTheta = [tRateSum+tT.*tRate,tRateSum+tRate-tT.*tRate].*pi.*2;tR = [tT.*0+i,tT.*0+i+1];if i == 1fill(cos(tTheta).*tR,sin(tTheta).*tR,CList(j,:),'EdgeColor',[1,1,1],'LineWidth',1)elsetCN = VAL(find(VAL(:,i) == j,1),i-1);tNN = j-VAL(find(VAL(:,i-1) == tCN,1),i)+1;tPN = length(unique(VAL(VAL(:,i-1) == tCN,i)));if mod(i,2)~=0tRN = tNN;elsetRN = tPN+1-tNN;end    NCList(j,:)=LCList(tCN,:).*0.8^(tRN-1)+[1,1,1].*(1-0.8^(tRN-1));fill(cos(tTheta).*tR,sin(tTheta).*tR,NCList(j,:),'EdgeColor',[1,1,1],'LineWidth',1)endrotation = (tRateSum+tRate/2)*360;if tRate > TextThresholdif rotation>90&&rotation<270rotation=rotation+180;text(cos((tRateSum+tRate/2).*pi.*2)*i,sin((tRateSum+tRate/2).*pi.*2)*i,NameCell{i}(j,:)+" ",FontProp{:},...'Rotation',rotation,'HorizontalAlignment','right')elsetext(cos((tRateSum+tRate/2).*pi.*2)*i,sin((tRateSum+tRate/2).*pi.*2)*i," "+NameCell{i}(j,:),FontProp{:},...'Rotation',rotation)endendtRateSum = tRateSum+tRate;endif i ~=1LCList=NCList;end
end
% 绘制最外圈饼状图
tRateSum = 0;
tNameCell = Table.(NameList{end-1});
NCList = zeros(size(VAL,1),3);
for j = 1:size(VAL,1)tRate = VAL(j,end)./sum(VAL(:,end)); tTheta = [tRateSum+tT.*tRate,tRateSum+tRate-tT.*tRate].*pi.*2;tR = [tT.*0+size(VAL,2),tT.*0+size(VAL,2)+1];tCN = VAL(j,size(VAL,2)-1);tNN = j-find(VAL(:,size(VAL,2)-1) == tCN,1)+1;tPN = sum(VAL(:,size(VAL,2)-1) == tCN);if mod(size(VAL,2),2)~=0tRN = tNN;elsetRN = tPN+1-tNN;endNCList(j,:)=LCList(tCN,:).*0.8^(tRN-1)+[1,1,1].*(1-0.8^(tRN-1));fill(cos(tTheta).*tR,sin(tTheta).*tR,NCList(j,:),'EdgeColor',[1,1,1],'LineWidth',1)rotation = (tRateSum+tRate/2)*360;if tRate > TextThresholdif rotation>90&&rotation<270rotation=rotation+180;text(cos((tRateSum+tRate/2).*pi.*2)*size(VAL,2),sin((tRateSum+tRate/2).*pi.*2)*size(VAL,2),tNameCell(IDX(j),:)+" ",FontProp{:},...'Rotation',rotation,'HorizontalAlignment','right')elsetext(cos((tRateSum+tRate/2).*pi.*2)*size(VAL,2),sin((tRateSum+tRate/2).*pi.*2)*size(VAL,2)," "+tNameCell(IDX(j),:),FontProp{:},...'Rotation',rotation)endendtRateSum = tRateSum+tRate;
end

完整代码(径向)

就是前面的代码顺着复制下来，可以换成自己的数据和颜色：
篇幅有限就只先给出径向渐变完整代码，切向渐变代码就光把最后一段稍微一更换即可：

% @author:slandarer%% ========================================================================
% 随机生成一组数据，可将 Table 自行更换
rng(1)
ULList = 'ABCD';
LLList = 'abcd';rowNum = 50;
idx1   = randi([1,length(ULList)],[1,rowNum]);
idx2   = randi([1,3],[rowNum,1]);
Name1  = ULList(idx1'*ones(1,4));
Name2  = [LLList(idx1'*ones(1,4)),num2str(idx2)];
Name3  = [LLList(idx1)',num2str(idx2),...char(45.*ones(rowNum,1)),char(randi([97,122],[rowNum,5]))];
Value  = rand([rowNum,1]);
Table  = table(Name1,Name2,Name3,Value);%% ========================================================================
% 利用 grp2idx 进行分类，利用 sortrows 将相同类归在一起
NameList  = Table.Properties.VariableNames;
NameNum   = length(NameList)-1;
NameCell{NameNum} = ' ';
valueList = zeros(length(Table.(NameList{1})),NameNum);
for i = 1:NameNum-1tName = Table.(NameList{i});tUniq = unique(tName,'rows');NameCell{i} = tUniq;ind = grp2idx([tUniq;tName]);ind(1:length(tUniq)) = [];valueList(:,i) = ind;
end
valueList(:,end) = -Table.(NameList{end});
[VAL,IDX] = sortrows(valueList,1:size(valueList,2));
VAL(:,end) = -VAL(:,end);
%% ========================================================================
% 此处可以设置配色
CList=[0.3882    0.5333    0.70591.0000    0.6824    0.20390.9373    0.4353    0.41570.5490    0.7608    0.79220.3333    0.6784    0.53730.7647    0.7373    0.24710.7333    0.4627    0.57650.7294    0.6275    0.58040.6627    0.7098    0.68240.4627    0.4627    0.4627];
% 在这可修改字体
FontProp = {'FontSize',12,'Color',[0,0,0]};
% 在这可设置比例低于多少的部分不显示文字
TextThreshold = 0.012;
%% ========================================================================
% 开始绘图
figure('Units','normalized','Position',[.2,.1,.52,.72]);
ax = gca; hold on
ax.DataAspectRatio = [1,1,1];
ax.XColor = 'none';
ax.YColor = 'none';
tT = linspace(0,1,100);
for i = 1:size(VAL,2)-1tRateSum = 0;tNum = length(NameCell{i});for j = 1:tNumtRate = sum(VAL(VAL(:,i) == j,end))./sum(VAL(:,end));  tTheta = [tRateSum+tT.*tRate,tRateSum+tRate-tT.*tRate].*pi.*2;tR = [tT.*0+i,tT.*0+i+1];if i == 1fill(cos(tTheta).*tR,sin(tTheta).*tR,CList(j,:),'EdgeColor',[1,1,1],'LineWidth',1)elsetCN = VAL(find(VAL(:,i) == j,1),1);fill(cos(tTheta).*tR,sin(tTheta).*tR,...CList(tCN,:).*0.8^(i-1)+[1,1,1].*(1-0.8^(i-1)),'EdgeColor',[1,1,1],'LineWidth',1)endrotation = (tRateSum+tRate/2)*360;if tRate > TextThresholdif rotation>90&&rotation<270rotation=rotation+180;text(cos((tRateSum+tRate/2).*pi.*2)*i,sin((tRateSum+tRate/2).*pi.*2)*i,NameCell{i}(j,:)+" ",FontProp{:},...'Rotation',rotation,'HorizontalAlignment','right')elsetext(cos((tRateSum+tRate/2).*pi.*2)*i,sin((tRateSum+tRate/2).*pi.*2)*i," "+NameCell{i}(j,:),FontProp{:},...'Rotation',rotation)endendtRateSum = tRateSum+tRate;end
end
% 绘制最外圈饼状图
tRateSum = 0;
tNameCell = Table.(NameList{end-1});
for j = 1:size(VAL,1)tRate = VAL(j,end)./sum(VAL(:,end)); tTheta = [tRateSum+tT.*tRate,tRateSum+tRate-tT.*tRate].*pi.*2;tR = [tT.*0+size(VAL,2),tT.*0+size(VAL,2)+1];tCN = VAL(j,1);fill(cos(tTheta).*tR,sin(tTheta).*tR,CList(tCN,:).*0.8^(size(VAL,2)-1)+...[1,1,1].*(1-0.8^(size(VAL,2)-1)),'EdgeColor',[1,1,1],'LineWidth',1)rotation = (tRateSum+tRate/2)*360;if tRate > TextThresholdif rotation>90&&rotation<270rotation=rotation+180;text(cos((tRateSum+tRate/2).*pi.*2)*size(VAL,2),sin((tRateSum+tRate/2).*pi.*2)*size(VAL,2),tNameCell(IDX(j),:)+" ",FontProp{:},...'Rotation',rotation,'HorizontalAlignment','right')elsetext(cos((tRateSum+tRate/2).*pi.*2)*size(VAL,2),sin((tRateSum+tRate/2).*pi.*2)*size(VAL,2)," "+tNameCell(IDX(j),:),FontProp{:},...'Rotation',rotation)endendtRateSum = tRateSum+tRate;
end

完

以上以及是完整代码，若日后代码有更新可在以下gitee仓库查看：

https://gitee.com/slandarer/spdraw/

这篇关于MATLAB | 两种上色方式的旭日图绘制的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/285844。

相关文章

Debezium 与 Apache Kafka 的集成方式步骤详解

Debezium 与 Apache Kafka 的集成方式步骤详解

《Debezium与ApacheKafka的集成方式步骤详解》本文详细介绍了如何将Debezium与ApacheKafka集成,包括集成概述、步骤、注意事项等,通过KafkaConnect,D... 目录一、集成概述二、集成步骤1. 准备 Kafka 环境2. 配置 Kafka Connect3. 安装 D
阅读更多...
Springboot中分析SQL性能的两种方式详解

Springboot中分析SQL性能的两种方式详解

《Springboot中分析SQL性能的两种方式详解》文章介绍了SQL性能分析的两种方式：MyBatis-Plus性能分析插件和p6spy框架,MyBatis-Plus插件配置简单,适用于开发和测试环... 目录SQL性能分析的两种方式:功能介绍实现方式:实现步骤:SQL性能分析的两种方式:功能介绍记录
阅读更多...
SQL 中多表查询的常见连接方式详解

SQL 中多表查询的常见连接方式详解

《SQL中多表查询的常见连接方式详解》本文介绍SQL中多表查询的常见连接方式,包括内连接（INNERJOIN）、左连接（LEFTJOIN）、右连接（RIGHTJOIN）、全外连接（FULLOUTER... 目录一、连接类型图表（ASCII 形式）二、前置代码（创建示例表）三、连接方式代码示例1. 内连接（I
阅读更多...
Android里面的Service种类以及启动方式

Android里面的Service种类以及启动方式

《Android里面的Service种类以及启动方式》Android中的Service分为前台服务和后台服务,前台服务需要亮身份牌并显示通知,后台服务则有启动方式选择,包括startService和b... 目录一句话总结：一、Service 的两种类型：1. 前台服务（必须亮身份牌）2. 后台服务（偷偷干
阅读更多...
Java深度学习库DJL实现Python的NumPy方式

Java深度学习库DJL实现Python的NumPy方式

《Java深度学习库DJL实现Python的NumPy方式》本文介绍了DJL库的背景和基本功能,包括NDArray的创建、数学运算、数据获取和设置等,同时,还展示了如何使用NDArray进行数据预处理... 目录1 NDArray 的背景介绍1.1 架构2 JavaDJL使用2.1 安装DJL2.2 基本操
阅读更多...
最长公共子序列问题的深度分析与Java实现方式

最长公共子序列问题的深度分析与Java实现方式

《最长公共子序列问题的深度分析与Java实现方式》本文详细介绍了最长公共子序列（LCS）问题,包括其概念、暴力解法、动态规划解法,并提供了Java代码实现,暴力解法虽然简单,但在大数据处理中效率较低,... 目录最长公共子序列问题概述问题理解与示例分析暴力解法思路与示例代码动态规划解法DP 表的构建与意义动
阅读更多...
JS 实现复制到剪贴板的几种方式小结

JS 实现复制到剪贴板的几种方式小结

《JS实现复制到剪贴板的几种方式小结》本文主要介绍了JS实现复制到剪贴板的几种方式小结,包括ClipboardAPI和document.execCommand这两种方法,具有一定的参考价值,感兴趣的... 目录一、Clipboard API相关属性方法二、document.execCommand优点：缺点：
阅读更多...
Python创建Excel的4种方式小结

Python创建Excel的4种方式小结

《Python创建Excel的4种方式小结》这篇文章主要为大家详细介绍了Python中创建Excel的4种常见方式,文中的示例代码简洁易懂,具有一定的参考价值,感兴趣的小伙伴可以学习一下... 目录库的安装代码1——pandas代码2——openpyxl代码3——xlsxwriterwww.cppcns.c
阅读更多...
Deepseek使用指南与提问优化策略方式

Deepseek使用指南与提问优化策略方式

《Deepseek使用指南与提问优化策略方式》本文介绍了DeepSeek语义搜索引擎的核心功能、集成方法及优化提问策略,通过自然语言处理和机器学习提供精准搜索结果,适用于智能客服、知识库检索等领域... 目录序言1. DeepSeek 概述2. DeepSeek 的集成与使用2.1 DeepSeek API
阅读更多...
CSS弹性布局常用设置方式

CSS弹性布局常用设置方式

《CSS弹性布局常用设置方式》文章总结了CSS布局与样式的常用属性和技巧,包括视口单位、弹性盒子布局、浮动元素、背景和边框样式、文本和阴影效果、溢出隐藏、定位以及背景渐变等,通过这些技巧,可以实现复杂... 一、单位元素vm 1vm 为视口的1%vh 视口高的1%vmin 参照长边vmax 参照长边re
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2