matlab绘制雷达图

Posted 2022-04-17 studyer_domi

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')