mradial.m

mradial对帕克变换有着重要的作用

function [th,mr,thm,mrm,thl,mlm,a,b]=mradial(id,iq)
%求平均曲线
long=length(id);
long=floor(long/16);
[theta,radial]=cart2pol(id(1:long),iq(1:long));
mr=zeros(1,41);th=zeros(1,41);thm=zeros(1,5);mrm=zeros(1,5);
thl=zeros(1,5);mlm=zeros(1,5);
for i=1:40;
    sum1=0;n1=0;
    for i1=1:length(theta)
        if theta(i1)>((i-21)*pi/20)&theta(i1)<=((i-20)*pi/20)
            sum1=sum1+radial(i1);n1=n1+1;
        end
    end
    th(i)=pi*(2*i-41)/40;
    if n1>0
        mr(i)=sum1/n1;
    end
end
for i=1:40
    if mr(i)==0
            j=i;k=i;
             if i<=20
                 while mr(j)==0
                        j=j+1;
                 end
                 while mr(k)==0
                     k=k-1;
                    if k<=0
                         k=40;
                    end
                end
            else
                while mr(j)==0
                    j=j+1;
                    if j>40
                        j=1;
                    end
                end
                while mr(k)==0
                    k=k-1;
                end
            end
        mr(i)=(mr(j)+mr(k))/2;
    end
end
mr(41)=mr(1);th(41)=th(1);
mmm=0;
for i=1:20
    if (mr(i)+mr(20+i))>mmm
        mmm=mr(i)+mr(20+i);
        mmmth=th(i);mi=i;
    end
    a=mmm;
end
for i=1:3
    thm(i)=mmmth;thm(6-i)=mmmth+pi;
end

mrm(1)=mr(mi);
mrm(2)=mr(mi)/2;
mrm(3)=0;
mrm(5)=mr(mi+20);
mrm(4)=mrm(5)/2;    
%++++++++++++++++++++++++++++++++++++++++++++++++++++短轴
mmin=mmm;
for i=1:20
      if (mr(i)+mr(20+i))<mmin
          mmin=mr(i)+mr(20+i);
          mmmth=th(i);mi=i;
      end
      b=mmin;
  end
  for i=1:3
      thl(i)=mmmth;thl(6-i)=mmmth+pi;
  end
mlm(1)=mr(mi);
mlm(2)=mr(mi)/2;
mlm(3)=0;
mlm(5)=mr(mi+20);
mlm(4)=mlm(5)/2;