matlab绘制雷达图MO组及NORMAL组 RNFL层、GCL层、IPL层、GCC层、ORL层做雷达图(共10张),比较NMO及normal组各层厚度差别,每层NMO及normal组颜色条的范围及颜色对应厚度值一样;
clearvars;close all;clc;
%% RNFL
nmo = [7.833333333 18.38888889 19.11111111 17.38888889 16.77777778 23.22222222 24.16666667 25.38888889 17.55555556];
normal = [11.18421053 23.86842105 24.57894737 20.13157895 17.10526316 39.39473684 40.89473684 47.10526316 18.76315789];
% M S1 I1 N1 T1 S2 I2 N2 T2
nmo = roundn(nmo,-4);
normal = roundn(normal,-4);
R = [0 0.5 1.5 3]; % (distance in km)
Az = [90 180 270 360 450]; % in degrees
diametre = 0:0.02:3;
degree = 90:1:450;
[xi,yi]=meshgrid(diametre,degree);
[xi1,yi1]=meshgrid(R,Az);
i=4;
valuenmo = [nmo(1) nmo(1) nmo(4) nmo(4+i);
nmo(1) nmo(1) nmo(2) nmo(2+i);
nmo(1) nmo(1) nmo(5) nmo(5+i);
nmo(1) nmo(1) nmo(3) nmo(3+i);
nmo(1) nmo(1) nmo(4) nmo(4+i);
];
valuenormal = [normal(1) normal(1) normal(4) normal(4+i);
normal(1) normal(1) normal(2) normal(2+i);
normal(1) normal(1) normal(5) normal(5+i);
normal(1) normal(1) normal(3) normal(3+i);
normal(1) normal(1) normal(4) normal(4+i);
];
% value_new=griddata(R,Az,value,xi,yi); %全场沿坐标面插值
valuenmo_new=griddata(xi1,yi1,valuenmo,xi,yi); %全场沿坐标面插值
valuenormal_new=griddata(xi1,yi1,valuenormal,xi,yi); %全场沿坐标面插值
figure(1)
[~,c1]= polarPcolor(diametre,degree,valuenmo_new,min(min(nmo),min(normal))-2,max(max(nmo),max(normal))+2);
% ylabel(c1,'nmo数值','fontSize',14);
% title('RNFL','fontSize',14);
set(gcf,'color','w')
figure(2)
[~,c2]= polarPcolor(diametre,degree,valuenormal_new,min(min(nmo),min(normal))-2,max(max(nmo),max(normal))+2);
% ylabel(c2,'normal数值','fontSize',14);
% title('RNFL','fontSize',14);
set(gcf,'color','w')
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删