MATLAB | 绘图复刻(九) | 泰勒图及组合泰勒图
2024-05-07 05:34:56
有粉丝问我这个图咋画：
我一看，这不就泰勒图嘛，就fileexchange上搜了一下泰勒图绘制代码，但是有的代码比较新的版本运行要改很多地方，有的代码需要包含一些压缩包没并没有的别人写的函数，于是我干脆自己写了个工具，想着如果大家用其他泰勒图绘制代码来不及修改的话，有多一次尝试机会吧。
复刻效果：
其他使用效果：
本文将分为大约四个部分：
- 1 讲解泰勒图原理
- 2 我写的工具函数基础使用教程
- 3 属性修饰教程
- 4 绘图复刻
- 5 工具函数完整代码
1 泰勒图原理
泰勒图是Karl E. Taylor于2001年首先提出，主要用来比较几个气象模式模拟的能力，因此该表示方法在气象领域使用最多，但是在其他自然科学领域也有一定的应用。
这个图巧妙地利用了相关系数(COR)、均方根误差(RMSE/RMSD)和标准差(STD)三个评价指标之间的关系，巧妙地将三个变量的数据安排在了二维图上。
讲解一下公式和关系哈，为了避免误会后文中各个评价指标均用首字母代写( C , R , S C,R,S C,R,S)，参考"标尺"使用小写字母 r r r，当前观察对象使用小写字母 f f f。
标准差
S r = ∑ n = 1 N ( r n − r ˉ ) 2 N S f = ∑ n = 1 N ( f n − f ˉ ) 2 N \begin{aligned} &S_r = \sqrt{\frac{\sum_{n = 1}^{N}(r_n-\bar{r})^2}{N}}\\ &S_f = \sqrt{\frac{\sum_{n = 1}^{N}(f_n-\bar{f})^2}{N}} \end{aligned} Sr=N∑n=1N(rn−rˉ)2 Sf=N∑n=1N(fn−fˉ)2
均方根误差
R = ∑ n = 1 N ( ( f n − f ˉ ) − ( r n − r ˉ ) ) 2 N R = \sqrt{\frac{\sum_{n = 1}^{N}((f_n-\bar{f})-(r_n-\bar{r}))^2}{N}} R=N∑n=1N((fn−fˉ)−(rn−rˉ))2
相关系数
C = ∑ n = 1 N ( f n − f ˉ ) ( r n − r ˉ ) N S r S f C = \frac{\sum_{n = 1}^{N}(f_n-\bar{f})(r_n-\bar{r})}{NS_rS_f} C=NSrSf∑n=1N(fn−fˉ)(rn−rˉ)
余弦关系
就是很简单的差的平方展开式 ( a − b ) 2 = a 2 + b 2 − 2 a b (a-b)^2=a^2+b^2-2ab (a−b)2=a2+b2−2ab，很容易能看出来：
R 2 = S r 2 + S f 2 − 2 S r S f C R^2=S_r^2+S_f^2-2S_rS_fC R2=Sr2+Sf2−2SrSfC
那我们如果把 C C C 看成 cos ( θ ) \cos(\theta) cos(θ) 这不就余弦公式嘛，因此可以将其安排在一个图中：
就可以同时看出两组数据的各种统计量的相近程度。
2 我写的工具函数使用教程
2.1 数据预处理
数据首先要与标准数据对比计算标准差、中心均方根误差、相关系数，使用SStats函数，该函数会放在文末。
clc; clear
Data = load('testData.mat');
Data = Data.Data(:,1:end-1);% Calculate STD RMSD and COR(计算标准差、中心均方根误差、相关系数)
STATS = zeros(4,size(Data,2));
for i = 1:size(Data,2)STATS(:,i) = SStats(Data(:,1),Data(:,i));
end
STATS(1,:) = [];
2.2 创建泰勒图坐标区域
% Create taylor axes(生成泰勒图坐标区域)
TD = STaylorDiag(STATS);
2.3 绘制散点
使用SPlot函数绘制，大概像下面这样(TD.SPlot)，改形状就改Marker属性即可：
% Color list(颜色列表)
colorList = [0.3569 0.0784 0.07840.6784 0.4471 0.17250.1020 0.3882 0.51760.1725 0.4196 0.43920.2824 0.2275 0.2902];
MarkerType={'o','diamond','pentagram','^','v'};
% Plot(绘制散点图)
for i = 1:size(Data,2)TD.SPlot(STATS(1,i),STATS(2,i),STATS(3,i),'Marker',MarkerType{i},'MarkerSize',15,...'Color',colorList(i,:),'MarkerFaceColor',colorList(i,:));
end
2.4 添加图例
添加图例并标注标准数据组(使用SText函数)：
% Legend
NameList = {'AAA','BBB','CCC','DDD','EEE'};
legend(NameList,'FontSize',13,'FontName','Times New Roman')% Annotation
TD.SText(STATS(1,1),STATS(2,1),STATS(3,1),{'reference';' '},'FontWeight','bold',...'FontSize',20,'FontName','Times New Roman','Color',colorList(1,:),...'VerticalAlignment','bottom','HorizontalAlignment','center')
2.5 标注每个点
循环使用SText函数即可：
for i = 1:size(Data,2)TD.SText(STATS(1,i),STATS(2,i),STATS(3,i)," "+string(NameList{i}),'FontWeight','bold',...'FontSize',14,'FontName','Times New Roman',...'VerticalAlignment','middle','HorizontalAlignment','left')
end
3 属性修饰
统一使用set函数进行修改，大概这么用：
- TD.set(对象名, 属性名, 属性值,…)
对象名就是由各个评价指标首字母( C , R , S C,R,S C,R,S)以及这个玩意的MATLAB内常用名称组合而成的，例如标准差轴就是SAxis，相关系数网格就是CGrid
3.1 基础属性
那些刻度长度啊，各个轴范围呀啥的要最先修饰。例如刻度长度：
TD.set('TickLength',[.015,.05])
各个轴范围：
TD.set('SLim',[0,300])
TD.set('RLim',[0,175])
主刻度及次刻度值：
TD.set('STickValues',0:50:300)
TD.set('SMinorTickValues',0:25:300)
TD.set('RTickValues',0:25:175)
TD.set('CTickValues',[.1,.2,.3,.4,.5,.6,.7,.8,.9,.95,.99])
3.2 网格属性
设置网格粗细及颜色：
% Set Grid(修饰各个网格)
TD.set('SGrid','Color',[.7,.7,.7],'LineWidth',1.5)
TD.set('RGrid','Color',[.77,.6,.18],'LineWidth',1.5)
TD.set('CGrid','Color',[0,0,.8],'LineStyle',':','LineWidth',.8);
3.3 标签属性
修饰轴标签及刻度标签：
% Set Tick Label(修饰刻度标签)
TD.set('STickLabelX','Color',[.8,0,0],'FontWeight','bold')
TD.set('STickLabelY','Color',[.8,0,0],'FontWeight','bold')
TD.set('RTickLabel','Color',[.77,.6,.18],'FontWeight','bold')
TD.set('CTickLabel','Color',[0,0,.8],'FontWeight','bold')% Set Label(修饰标签)
TD.set('SLabel','Color',[.8,0,0],'FontWeight','bold')
TD.set('CLabel','Color',[0,0,.8],'FontWeight','bold')
3.4 轴属性
设置轴颜色、粗细：
% Set Axis(修饰各个轴)
TD.set('SAxis','Color',[.8,0,0],'LineWidth',2)
TD.set('CAxis','Color',[0,0,.8],'LineWidth',2)
3.5 刻度属性
刻度粗细、颜色、线型啥的都能设置：
% Set Tick and MinorTick(修饰主次刻度)
TD.set('STick','Color',[.8,0,0],'LineWidth',8)
TD.set('CTick','Color',[0,0,.8],'LineWidth',8)
TD.set('SMinorTick','Color',[.8,0,0],'LineWidth',5)
TD.set('CMinorTick','Color',[0,0,.8],'LineWidth',5)
4 绘图复刻
4.1 数据说明
没去找原始论文数据，这里随便用了组数据且四个坐标区域用的数据也一样，到时候大家可自行更换：
clc; clear
Data = load('testData.mat');
Data = Data.Data;
% % 不想下数据的可以用随便捏造的数据试试
% % 随便捏造了点数据(Made up some data casually)
% Data = randn(100,6).*.2+[(linspace(-1,.5,100)').*ones(1,2),...
% (linspace(-.5,.7,100)').*ones(1,2),...
% (linspace(-.9,.2,100)').*ones(1,2)];% Calculate STD RMSD and COR(计算标准差、中心均方根误差、相关系数)
STATS = zeros(4,size(Data,2));
for i = 1:size(Data,2)STATS(:,i) = SStats(Data(:,1),Data(:,i));
end
STATS(1,:) = [];
4.2 子图创建
硬创建呗：
%% Create figure and axes(生成基础布局)
fig=figure('Units','normalized','Position',[.2,.1,.54,.74]);
bkgAx=axes(fig,'Position',[.13,.11,.67,.88],'XTick',[],'YTick',[],'Box','on',...'Color',[253,228,203]./255,'XLim',[0,100],'YLim',[0,88],'NextPlot','add');
plot(bkgAx,[50,50],[0,88],'Color','k','LineWidth',.8)
text(bkgAx,25,42,'Autumn','FontSize',16,'FontName','Times New Roman','HorizontalAlignment','center')
text(bkgAx,75,42,'Winter','FontSize',16,'FontName','Times New Roman','HorizontalAlignment','center')
text(bkgAx,25,86,'Spring','FontSize',16,'FontName','Times New Roman','HorizontalAlignment','center')
text(bkgAx,75,86,'Summer','FontSize',16,'FontName','Times New Roman','HorizontalAlignment','center')
bkgAx.XLabel.String='Standard Deviation(mm month^{-1})';
bkgAx.XLabel.FontSize=18;
bkgAx.XLabel.FontName='Times New Roman';
bkgAx.XLabel.Position=[50,-3,-1];
bkgAx.YLabel.String='Standard Deviation(mm month^{-1})';
bkgAx.YLabel.FontSize=18;
bkgAx.YLabel.FontName='Times New Roman';
bkgAx.YLabel.Position=[-7,44,-1];
% -------------------------------------------------------------------------
ax1 = axes(fig,'Position',[.13,.11,.335,.4],'Box','on');
ax2 = axes(fig,'Position',[.465,.11,.335,.4],'Box','on');
ax3 = axes(fig,'Position',[.13,.55,.335,.4],'Box','on');
ax4 = axes(fig,'Position',[.465,.55,.335,.4],'Box','on');
4.3 创建泰勒图坐标区域
可以将第一个参数换成坐标区域对象，就能将特定坐标区域修改成泰勒图坐标区域，以下是生成后做了点属性调整：
%% Create taylor axes(生成泰勒图坐标区域)
TD1 = STaylorDiag(ax1,STATS);
set(ax1,'Box','on','XTick',[],'YTick',[],'XColor','k','YColor','k',...'XLim',TD1.SLim.*1.13,'YLim',TD1.SLim.*1.15);
TD1.set('SLabel','Color','none');
TD1.set('RGrid','Color',[.77,.6,.18],'LineWidth',1.2)
TD1.set('RTickLabel','Color',[.77,.6,.18],'FontWeight','bold')
TD1.set('CLabel','FontSize',16)TD2 = STaylorDiag(ax2,STATS);
set(ax2,'Box','on','XTick',[],'YTick',[],'XColor','k','YColor','k',...'XLim',TD2.SLim.*1.13,'YLim',TD2.SLim.*1.15);
TD2.set('SLabel','Color','none');
TD2.set('STickLabelY','Color','none');
TD2.set('RGrid','Color',[.77,.6,.18],'LineWidth',1.2)
TD2.set('RTickLabel','Color',[.77,.6,.18],'FontWeight','bold')
TD2.set('CLabel','FontSize',16)TD3 = STaylorDiag(ax3,STATS);
set(ax3,'Box','on','XTick',[],'YTick',[],'XColor','k','YColor','k',...'XLim',TD2.SLim.*1.13,'YLim',TD2.SLim.*1.15);
TD3.set('SLabel','Color','none');
TD3.set('STickLabelX','Color','none');
TD3.set('RGrid','Color',[.77,.6,.18],'LineWidth',1.2)
TD3.set('RTickLabel','Color',[.77,.6,.18],'FontWeight','bold')
TD3.set('CLabel','FontSize',16)TD4 = STaylorDiag(ax4,STATS);
set(ax4,'Box','on','XTick',[],'YTick',[],'XColor','k','YColor','k',...'XLim',TD2.SLim.*1.13,'YLim',TD2.SLim.*1.15);
TD4.set('SLabel','Color','none');
TD4.set('STickLabelX','Color','none');
TD4.set('STickLabelY','Color','none');
TD4.set('RGrid','Color',[.77,.6,.18],'LineWidth',1.2)
TD4.set('RTickLabel','Color',[.77,.6,.18],'FontWeight','bold')
TD4.set('CLabel','FontSize',16)
4.4 绘制散点
最后一组散点绘制后对象被存到lgdPltHdl的矩阵，这是为了之后生成图例方便(只显示第2-6组数据的图例，不显示标准数据)：
%% 绘制散点
% Color list(颜色列表)
colorList = [145,81,155;217,34,30;68,127,183;76,181,75;145,81,155;248,130,7]./255;
% Plot(绘制散点图)
for i = 1:size(Data,2)TD1.SPlot(STATS(1,i),STATS(2,i),STATS(3,i),'Marker','o','MarkerSize',10,...'Color',colorList(i,:),'MarkerFaceColor',colorList(i,:));TD2.SPlot(STATS(1,i),STATS(2,i),STATS(3,i),'Marker','o','MarkerSize',10,...'Color',colorList(i,:),'MarkerFaceColor',colorList(i,:));TD3.SPlot(STATS(1,i),STATS(2,i),STATS(3,i),'Marker','o','MarkerSize',10,...'Color',colorList(i,:),'MarkerFaceColor',colorList(i,:));lgdPltHdl(i)=TD4.SPlot(STATS(1,i),STATS(2,i),STATS(3,i),'Marker','o','MarkerSize',10,...'Color',colorList(i,:),'MarkerFaceColor',colorList(i,:));
end
4.5 添加图例和标注
% legend
lgdHdl = legend(ax4,lgdPltHdl(2:end),{'# 1 1/x','# 2 e\^-x(x/0.1\^2)','# 3 e\^-x(x/0.5\^2)','# 4 e\^-x(x/0.2\^2)','# 5 e\^-x(x/0.25\^2)'},...'Position',[.825,.47,0.15,0.2]);
lgdHdl.Title.String='Weighting Scheme';
lgdHdl.Title.FontSize=14;
lgdHdl.Title.FontName='Times New Roman';
lgdHdl.FontSize=12;
lgdHdl.FontName='Times New Roman';
lgdHdl.Box='off';
% Annotation
TD1.SText(STATS(1,1),STATS(2,1),STATS(3,1),{'observed';''},'FontWeight','bold',...'FontSize',13,'FontName','Times New Roman','Color',colorList(1,:),...'VerticalAlignment','bottom','HorizontalAlignment','center');
TD2.SText(STATS(1,1),STATS(2,1),STATS(3,1),{'observed';''},'FontWeight','bold',...'FontSize',13,'FontName','Times New Roman','Color',colorList(1,:),...'VerticalAlignment','bottom','HorizontalAlignment','center');
TD3.SText(STATS(1,1),STATS(2,1),STATS(3,1),{'observed';''},'FontWeight','bold',...'FontSize',13,'FontName','Times New Roman','Color',colorList(1,:),...'VerticalAlignment','bottom','HorizontalAlignment','center');
TD4.SText(STATS(1,1),STATS(2,1),STATS(3,1),{'observed';''},'FontWeight','bold',...'FontSize',13,'FontName','Times New Roman','Color',colorList(1,:),...'VerticalAlignment','bottom','HorizontalAlignment','center');
text(ax1,TD1.SLim(2),TD1.SLim(2),'RMS error','Color',[.77,.6,.18],...'VerticalAlignment','bottom','HorizontalAlignment','right',...'FontSize',15,'FontName','Times New Roman','FontWeight','bold')
text(ax2,TD1.SLim(2),TD1.SLim(2),'RMS error','Color',[.77,.6,.18],...'VerticalAlignment','bottom','HorizontalAlignment','right',...'FontSize',15,'FontName','Times New Roman','FontWeight','bold')
text(ax3,TD1.SLim(2),TD1.SLim(2),'RMS error','Color',[.77,.6,.18],...'VerticalAlignment','bottom','HorizontalAlignment','right',...'FontSize',15,'FontName','Times New Roman','FontWeight','bold')
text(ax4,TD1.SLim(2),TD1.SLim(2),'RMS error','Color',[.77,.6,.18],...'VerticalAlignment','bottom','HorizontalAlignment','right',...'FontSize',15,'FontName','Times New Roman','FontWeight','bold')
5 工具函数完整代码
5.1 SStats
function Data=SStats(Cr,Cf)
% Copyright (c) 2023, Zhaoxu Liu / slandarer
% Zhaoxu Liu / slandarer (2023). taylor diagram class
% (https://www.mathworks.com/matlabcentral/fileexchange/130889-taylor-diagram-class),
% MATLAB Central File Exchange.
Cr = Cr(:);
Cf = Cf(:);
nanInd = isnan(Cr)|isnan(Cf);
Cr(nanInd) = [];
Cf(nanInd) = [];% N = length(Cf);
MEAN = mean(Cf);
STD = std(Cf,1);
RMSD = std(Cf-Cr,1);
COR = corrcoef(Cf,Cr);
Data = [MEAN, STD, RMSD, COR(1,2)].';%% calculation formula of STD RMSD and COR
% N = length(Cf);
%
% / sum[ {C-mean(C)} .^2] \
% STD = sqrt| --------------------- |
% \ N /
%
% Equivalent calculation formula:
% STD = sqrt(sum((Cf-mean(Cf)).^2)./N);
% STD = norm(Cf-mean(Cf))./sqrt(N);
% STD = rms(Cf-mean(Cf));
% STD = sqrt(var(Cf,1));
% STD = std(Cf,1);% / sum[ { [C-mean(C)] - [Cr-mean(Cr)] }.^2 ] \
% RMSD = sqrt| ------------------------------------------- |
% \ N /
%
% Equivalent calculation formula:
% RMSD = sqrt(sum((Cf-mean(Cf)-(Cr-mean(Cr))).^2)./N);
% RMSD = norm(Cf-mean(Cf)-Cr+mean(Cr))./sqrt(N);
% RMSD = rmse(Cr-mean(Cr),Cf-mean(Cf));
% RMSD = rms(Cf-mean(Cf)-Cr+mean(Cr));
% RMSD = sqrt(var(Cf-Cr,1));
% RMSD = std(Cf-Cr,1);% sum( [C-mean(C)].*[Cr-mean(Cr)] )
% COR = ---------------------------------
% N*STD(C)*STD(Cr)
%
% Equivalent calculation formula:
% COR = sum((Cf-mean(Cf)).*(Cr-mean(Cr)))./N./std(Cf,1)./std(Cr,1);
% COR = (Cf-mean(Cf)).'*(Cr-mean(Cr))./N./std(Cf,1)./std(Cr,1);
% COR = [1,0]*cov(Cf,Cr,1)./std(Cf,1)./std(Cr,1)*[0;1];
% COR = [1,0]*corrcoef(Cf,Cr)*[0;1];
end
5.2 STaylorDiag
classdef STaylorDiag < handle
% Copyright (c) 2023, Zhaoxu Liu / slandarer
% Zhaoxu Liu / slandarer (2023). taylor diagram class
% (https://www.mathworks.com/matlabcentral/fileexchange/130889-taylor-diagram-class),
% MATLAB Central File Exchange.STATS; % STD RMSD and CORax; TickLength=[.015,.008];% S-Axis means STD AxisSAxis; RAxis; CAxis;SGrid; RGrid; CGrid;SLabel; RLabel; CLabel; STickLabelX; STickLabelY; RTickLabel; CTickLabel;% STD Tick and MinorTickSTick; SMinorTick; STickValues; SMinorTickValues; SLim;%RTickValues; RLim;%CTick; CMinorTick; CLim = [0,1];CTickValues=[.1,.3,.5,.7,.9,.95,.99];CMinorTickValues=[.05,.15,.2,.25,.35,.4,.45,.55,.6,.65,.75,.8,.85,.91,.92,.93,.94,.96,.97,.98,1]; endmethodsfunction obj = STaylorDiag(varargin)if isa(varargin{1},'matlab.graphics.axis.Axes')obj.ax = varargin{1}; varargin(1) = [];elseobj.ax = gca;endobj.STATS = varargin{1};obj.ax.Parent.Color = [1,1,1];obj.ax.NextPlot = 'add';obj.ax.XGrid = 'off';obj.ax.YGrid = 'off';obj.ax.Box = 'off';obj.ax.DataAspectRatio = [1,1,1];% -------------------------------------------------------------obj.SLim = [0,max(obj.STATS(1,:)).*1.15];obj.ax.XLim = obj.SLim;obj.STickValues = obj.ax.XTick;obj.RLim = [0,max(obj.STATS(2,:)).*1.15];obj.ax.XLim = obj.RLim;obj.RTickValues = obj.ax.XTick;obj.RTickValues(1) = [];obj.ax.XLim = obj.SLim; obj.ax.YLim = obj.SLim;obj.SMinorTickValues = ...linspace(obj.STickValues(1),obj.STickValues(end),(length(obj.STickValues)-1).*5+1);obj.SMinorTickValues=setdiff(obj.SMinorTickValues,obj.STickValues);obj.CMinorTickValues=setdiff(obj.CMinorTickValues,obj.CTickValues);% -------------------------------------------------------------% STD Tick, MinorTick, Grid, TickLabel[tYTickX,tYTickY,tYTickMX,tYTickMY,tRGridX,tRGridY] = obj.getSValue();obj.SGrid = plot(obj.ax,tRGridX(:),tRGridY(:),'LineWidth',.5,'Color','k');obj.SAxis = plot(obj.ax,[0,0,obj.SLim(2)],[obj.SLim(2),0,0],'LineWidth',1.2,'Color','k');obj.STick = plot(obj.ax,[tYTickX(:);tYTickY(:)],[tYTickY(:);tYTickX(:)],'LineWidth',.8,'Color','k');obj.SMinorTick = plot(obj.ax,[tYTickMX(:);tYTickMY(:)],[tYTickMY(:);tYTickMX(:)],'LineWidth',.8,'Color','k');obj.STickLabelX = text(obj.ax,obj.STickValues,0.*obj.STickValues,string(obj.STickValues),...'FontName','Times New Roman','FontSize',13,'VerticalAlignment','top','HorizontalAlignment','center');obj.STickLabelY = text(obj.ax,0.*obj.STickValues(2:end),obj.STickValues(2:end),string(obj.STickValues(2:end))+" ",...'FontName','Times New Roman','FontSize',13,'VerticalAlignment','middle','HorizontalAlignment','right');obj.SLabel = text(obj.ax,obj.ax.YLabel.Position(1),obj.ax.YLabel.Position(2),'Standard Deviation',...'FontName','Times New Roman','FontSize',20,'VerticalAlignment','bottom','HorizontalAlignment','center','Rotation',90);obj.SGrid.Annotation.LegendInformation.IconDisplayStyle='off';obj.SAxis.Annotation.LegendInformation.IconDisplayStyle='off';obj.STick.Annotation.LegendInformation.IconDisplayStyle='off';obj.SMinorTick.Annotation.LegendInformation.IconDisplayStyle='off';% -------------------------------------------------------------% RMSD Grid TickLabel[tR2GridX,tR2GridY] = obj.getRValue();obj.RGrid = plot(obj.ax,tR2GridX(:),tR2GridY(:),'LineWidth',.5,'Color','k','LineStyle','--');obj.RTickLabel = text(obj.ax,cos(pi*5/6).*obj.RTickValues+obj.STATS(1,1),sin(pi*5/6).*obj.RTickValues,string(obj.RTickValues),...'FontName','Times New Roman','FontSize',13,'VerticalAlignment','bottom','HorizontalAlignment','center','Rotation',60);obj.RGrid.Annotation.LegendInformation.IconDisplayStyle='off';% -------------------------------------------------------------% RMSD^2 = STD_r^2 + STD_f^2 - 2*STD_r*STD_f*COR%% STD_r^2 + STD_f^2 - RMSD^2% COR = ----------------------------% 2*STD_r*STD_f%[tCAxisX,tCAxisY,tTRMSD,tCTickX,tCTickY,tCTickMX,tCTickMY,tRMSDGridX,tRMSDGridY]=obj.getCValue();obj.CGrid = plot(obj.ax,tRMSDGridX(:),tRMSDGridY(:),'LineWidth',.5,'Color',[.8,.8,.8],'LineStyle','-');obj.CAxis = plot(obj.ax,tCAxisX,tCAxisY,'LineWidth',1.2,'Color','k');obj.CTick = plot(obj.ax,tCTickX(:),tCTickY(:),'LineWidth',.8,'Color','k');obj.CMinorTick = plot(obj.ax,tCTickMX(:),tCTickMY(:),'LineWidth',.8,'Color','k');for i=1:length(tTRMSD)obj.CTickLabel{i} = text(obj.ax,cos(tTRMSD(i)).*obj.SLim(2),sin(tTRMSD(i)).*obj.SLim(2)," "+string(obj.CTickValues(i)),...'FontName','Times New Roman','FontSize',13,'VerticalAlignment','middle','HorizontalAlignment','left','Rotation',tTRMSD(i)./pi.*2.*90);endobj.CLabel = text(obj.ax,cos(pi/4).*(obj.SLim(2)-obj.ax.YLabel.Position(1)),...sin(pi/4).*(obj.SLim(2)-obj.ax.YLabel.Position(1)),'Correlation',...'FontName','Times New Roman','FontSize',20,'VerticalAlignment','bottom','HorizontalAlignment','center','Rotation',-45);obj.CGrid.Annotation.LegendInformation.IconDisplayStyle='off';obj.CAxis.Annotation.LegendInformation.IconDisplayStyle='off';obj.CTick.Annotation.LegendInformation.IconDisplayStyle='off';obj.CMinorTick.Annotation.LegendInformation.IconDisplayStyle='off';% -------------------------------------------------------------obj.ax.XColor = 'none';obj.ax.YColor = 'none';end% =================================================================function [tYTickX,tYTickY,tYTickMX,tYTickMY,tRGridX,tRGridY]=getSValue(obj)tTPi_2 = linspace(0,pi/2,200).';tYTickX = [0.*obj.STickValues;obj.STickValues.*0+obj.TickLength(1).*obj.SLim(2);nan.*obj.STickValues];tYTickY = [obj.STickValues;obj.STickValues;nan.*obj.STickValues];tYTickMX = [0.*obj.SMinorTickValues;obj.SMinorTickValues.*0+obj.TickLength(2).*obj.SLim(2);nan.*obj.SMinorTickValues];tYTickMY = [obj.SMinorTickValues;obj.SMinorTickValues;nan.*obj.SMinorTickValues];tRGridX = [obj.STickValues.*cos(repmat(tTPi_2,[1,length(obj.STickValues)]));nan.*obj.STickValues];tRGridY = [obj.STickValues.*sin(repmat(tTPi_2,[1,length(obj.STickValues)]));nan.*obj.STickValues];endfunction [tR2GridX,tR2GridY]=getRValue(obj)tTPi = linspace(0,pi,200).';tR2GridX = [obj.RTickValues.*cos(repmat(tTPi,[1,length(obj.RTickValues)]));nan.*obj.RTickValues]+obj.STATS(1,1);tR2GridY = [obj.RTickValues.*sin(repmat(tTPi,[1,length(obj.RTickValues)]));nan.*obj.RTickValues];tR2GridN = sqrt(tR2GridX.^2+tR2GridY.^2)>obj.SLim(2);tR2GridX(tR2GridN) = nan; tR2GridY(tR2GridN) = nan;endfunction [tCAxisX,tCAxisY,tTRMSD,tCTickX,tCTickY,tCTickMX,tCTickMY,tRMSDGridX,tRMSDGridY]=getCValue(obj)tTPi_2 = linspace(0,pi/2,200).';tCAxisX = cos(tTPi_2).*obj.SLim(2);tCAxisY = sin(tTPi_2).*obj.SLim(2);tRMSD = sqrt(2.*(obj.STATS(1,1).^2)-2.*(obj.STATS(1,1).^2).*obj.CTickValues);tTRMSD = asin(tRMSD./2./obj.STATS(1,1)).*2;tCTickX = [cos(tTRMSD).*obj.SLim(2);cos(tTRMSD).*obj.SLim(2).*(1-obj.TickLength(1));tTRMSD.*nan];tCTickY = [sin(tTRMSD).*obj.SLim(2);sin(tTRMSD).*obj.SLim(2).*(1-obj.TickLength(1));tTRMSD.*nan];tRMSDM = sqrt(2.*(obj.STATS(1,1).^2)-2.*(obj.STATS(1,1).^2).*obj.CMinorTickValues);tTRMSDM = asin(tRMSDM./2./obj.STATS(1,1)).*2;tCTickMX = [cos(tTRMSDM).*obj.SLim(2);cos(tTRMSDM).*obj.SLim(2).*(1-obj.TickLength(2));tTRMSDM.*nan];tCTickMY = [sin(tTRMSDM).*obj.SLim(2);sin(tTRMSDM).*obj.SLim(2).*(1-obj.TickLength(2));tTRMSDM.*nan];tRMSDGridX = [cos(tTRMSD).*obj.SLim(2);0.*tTRMSD;nan.*tTRMSD];tRMSDGridY = [sin(tTRMSD).*obj.SLim(2);0.*tTRMSD;nan.*tTRMSD];end% =================================================================function pltHdl=SPlot(obj,STD,RMSD,~,varargin)tTheta = 2.*asin(RMSD./2./STD);pltHdl = plot(obj.ax,cos(tTheta).*STD,sin(tTheta).*STD,varargin{:},'LineStyle','none');endfunction txtHdl=SText(obj,STD,RMSD,~,varargin)tTheta = 2.*asin(RMSD./2./STD);txtHdl = text(obj.ax,cos(tTheta).*STD,sin(tTheta).*STD,varargin{:});end% =================================================================function set(obj,target,varargin)if isa(varargin{1},'char')||isa(varargin{1},'string')switch targetcase {'SAxis','CAxis','SLabel','CLabel','STick','CTick','SMinorTick','CMinorTick','RTickLabel','SGrid','RGrid','CGrid'}set(obj.(target),varargin{:});case 'STickLabelX'set(obj.STickLabelX,varargin{:});case 'STickLabelY'set(obj.STickLabelY,varargin{:});case 'CTickLabel'for i=1:length(obj.CTickLabel)set(obj.CTickLabel{i},varargin{:});endendelseoriRLim = obj.RLim;obj.(target) = varargin{1};obj.ax.XColor = 'k';obj.ax.YColor = 'k';if abs(obj.SLim(2)-obj.ax.XLim(2))>eps||abs(obj.RLim(2)-oriRLim(2))>epsobj.ax.XLim = obj.SLim;obj.STickValues = obj.ax.XTick;obj.ax.XLim = obj.RLim;obj.RTickValues = obj.ax.XTick;obj.RTickValues(1) = [];obj.ax.XLim = obj.SLim; obj.ax.YLim = obj.SLim;obj.SMinorTickValues = ...linspace(obj.STickValues(1),obj.STickValues(end),(length(obj.STickValues)-1).*5+1);endobj.(target) = varargin{1};obj.SMinorTickValues=setdiff(obj.SMinorTickValues,obj.STickValues);obj.CMinorTickValues=setdiff(obj.CMinorTickValues,obj.CTickValues); %[tYTickX,tYTickY,tYTickMX,tYTickMY,tRGridX,tRGridY] = obj.getSValue();set(obj.SGrid,'XData',tRGridX(:),'YData',tRGridY(:));set(obj.SAxis,'XData',[0,0,obj.SLim(2)],'YData',[obj.SLim(2),0,0]);set(obj.STick,'XData',[tYTickX(:);tYTickY(:)],'YData',[tYTickY(:);tYTickX(:)]);set(obj.SMinorTick,'XData',[tYTickMX(:);tYTickMY(:)],'YData',[tYTickMY(:);tYTickMX(:)]);delete(obj.STickLabelX);delete(obj.STickLabelY);delete(obj.SLabel)obj.STickLabelX = text(obj.ax,obj.STickValues,0.*obj.STickValues,string(obj.STickValues),...'FontName','Times New Roman','FontSize',13,'VerticalAlignment','top','HorizontalAlignment','center');obj.STickLabelY = text(obj.ax,0.*obj.STickValues(2:end),obj.STickValues(2:end),string(obj.STickValues(2:end))+" ",...'FontName','Times New Roman','FontSize',13,'VerticalAlignment','middle','HorizontalAlignment','right');obj.SLabel = text(obj.ax,obj.ax.YLabel.Position(1),obj.ax.YLabel.Position(2),'Standard Deviation',...'FontName','Times New Roman','FontSize',20,'VerticalAlignment','bottom','HorizontalAlignment','center','Rotation',90);%[tR2GridX,tR2GridY] = obj.getRValue();set(obj.RGrid,'XData',tR2GridX(:),'YData',tR2GridY(:));delete(obj.RTickLabel)obj.RTickLabel = text(obj.ax,cos(pi*5/6).*obj.RTickValues+obj.STATS(1,1),sin(pi*5/6).*obj.RTickValues,string(obj.RTickValues),...'FontName','Times New Roman','FontSize',13,'VerticalAlignment','bottom','HorizontalAlignment','center','Rotation',60);%[tCAxisX,tCAxisY,tTRMSD,tCTickX,tCTickY,tCTickMX,tCTickMY,tRMSDGridX,tRMSDGridY]=obj.getCValue();set(obj.CGrid,'XData',tRMSDGridX(:),'YData',tRMSDGridY(:));set(obj.CAxis,'XData',tCAxisX,'YData',tCAxisY);set(obj.CTick,'XData',tCTickX(:),'YData',tCTickY(:));set(obj.CMinorTick,'XData',tCTickMX(:),'YData',tCTickMY(:));for i=length(obj.CTickLabel):-1:1delete(obj.CTickLabel{i});enddelete(obj.CLabel);for i=1:length(tTRMSD)obj.CTickLabel{i} = text(obj.ax,cos(tTRMSD(i)).*obj.SLim(2),sin(tTRMSD(i)).*obj.SLim(2)," "+string(obj.CTickValues(i)),...'FontName','Times New Roman','FontSize',13,'VerticalAlignment','middle','HorizontalAlignment','left','Rotation',tTRMSD(i)./pi.*2.*90);endobj.CLabel = text(obj.ax,cos(pi/4).*(obj.SLim(2)-obj.ax.YLabel.Position(1)),...sin(pi/4).*(obj.SLim(2)-obj.ax.YLabel.Position(1)),'Correlation',...'FontName','Times New Roman','FontSize',20,'VerticalAlignment','bottom','HorizontalAlignment','center','Rotation',-45);obj.ax.XColor = 'none';obj.ax.YColor = 'none';endendend
end
完
未经允许本代码请勿作商业用途，引用的话可以引用我file exchange上的链接，可使用如下格式：
Zhaoxu Liu / slandarer (2023). taylor diagram class (https://www.mathworks.com/matlabcentral/fileexchange/130889-taylor-diagram-class), MATLAB Central File Exchange. 检索来源 2023/6/11.
若转载请保留以上file exchange链接及本文链接！！！！!
本文全部代码已同时上传gitee仓库，可从gitee或者fileexchange下载：
https://gitee.com/slandarer/PLTreprint
MATLAB | 绘图复刻(九) | 泰勒图及组合泰勒图相关推荐
- MATLAB | 绘图复刻(二) | 折线图+误差棒+柱状图+散点抖动+灰色背景+图片叠加
看到gzh R语言ggplot2科研绘图发布了一篇绘图复刻类文章,复刻了: Nature(IF=49.962)文章(Gut microbiota modulates weight gain in mi ...
- MATLAB | 绘图复刻(三) | 分层聚类分析图：树状图+热图
好久不见啊,今天时绘图复刻第三期,这期画的比较难应该文章也不会太短... 前段时间发现公众号SCIPainter发布了一期名为<如何对基因和蛋白质的表达丰度进行相关性分析>,其中有一幅图很 ...
- MATLAB | 绘图复刻(十) | 旋转45度的相关系数热图绘制
我们在论文里经常能见到这样三角形的相关性热图,这样的热图绘制出来会更节省空间: 这期就教大家如何绘制三角形的热图,绘制效果如下: 以下是教程部分: 0 数据准备 这里随便生成了一组数据,计算了相关系数 ...
- funkyheatmap | 用这个包来完美复刻Nature Biotechnology的高颜值神图吧！~
1写在前面 天气开始暖和了☀️,发现旅游的人好多啊!~
- MATLAB | 全网唯一 MATLAB双向弦图(有向弦图)绘制
先赞后看,养成习惯~~ 先赞后看,养成习惯~~ 先赞后看,养成习惯~~ 绘制效果 下面这款弦图我已经出了很久了,也陆陆续续增添了新的功能和修了一些bug: 甚至还用它做出了一些复刻,分成两组的弦图有了 ...
- zk4元年拆解_科比zk4复刻前掌没有zoom zk4选秀日复刻中底拆解测评
科比zk4的复刻回归本来是一件很兴奋的事情,但是网上就有很多小伙伴一直在讨论zk4复刻是否和元年的配置一样,有的小伙伴希望是元年的配置,有的希望是全掌zoom,然而最近拆解图出来了,说实话小编看了还是 ...
- R语言：ggplot2画带误差棒的组合折线图教程。
1.cowplotR包的安装,输入代码install.packages("cowplot"),安装成功后,用library("cowplot")调用该R包. 2 ...
- MATLAB | 我用MATLAB复刻了这场颜色图生成大赛的作品
我用MATLAB复刻了这场颜色图绘制大赛的作品 这场比赛居然已经是七年前的事情了,大佬A Frayed Knot在StackExchange上发起了挑战,举办了名为Tweetable Mathemat ...
- 超火的举牌加油小人生成器MATLAB复刻版(app designer)
突然发现cla函数也可以应用到app designer控件上,因而对部分内容做出更改. 手痒,用matlab app designer复刻了一款举牌加油小人生成器,效果如下: 该程序可以生成正方形和长 ...
最新文章
- DVWA-暴力破解-对‘g0tmi1k’文章的学习笔记
- 电脑任务栏跑到右边去了_电脑没有声音怎么解决 电脑没有声音解决方法【详解】...
- 用Markdown写博客快速入门
- 基于迁移学习的农作物病虫害检测方法研究与应用
- asp.net web.config配置节说明
- 替换url中某个参数的值或是添加某个参数的方法(js 分页上下页可以使用)
- Coolite 换肤
- java实现杨辉三角
- 12306网上买火车票选择上中下铺的方法
- 视频转码（Java）
- Linux系统如何查看内存
- 家用简易NAS系统搭建
- 吃瓜 || 一文看懂BCH分叉事件始末
- 鬼压床、清明梦、外星劫持……睡眠麻痹症能够解释这些都市传说吗？
- 大数据实训整体解决方案
- 电信运营商移动互联网发展分析
- JAVA毕业设计共享充电宝管理系统计算机源码+lw文档+系统+调试部署+数据库
- 11.22 点餐APP第一阶段总结
- 邻接表生成模型头文件C语言
- 安装widows XP时蓝屏代码0X0000007B的问题的解决方案