ciaphase.m

来自「rang doppler imaging and motion compensa」· M 代码 · 共 41 行

     function [startpp,endpp] = CIAphase(data,startp,endp,maxdev,stdev)
     % 
     % Finds starting and ending points of a segment using
     % the linearity of the phase for complex image analysis.
     %
     
     [m,n]=size(endp);
     if n == 1
     n = m;
     m = 1;
     endp = endp(:);
     end
     sizend=n;

     [m,n]=size(data);
     if n == 1
     n = m;
     m = 1;
     data = data(:);
     end

     for an=1:sizend
     startend=startp(an):endp(an);
     sizestend=size(startend,2);
     auxfd=abs(data(startend)); 
     ampv=mean(auxfd)*sizestend;
     pwr=unwrap(angle(data(startend)));
     dpwr=diff(pwr);

     auxpwr=mean(dpwr);
     devpwr=std(dpwr);
     maxdp=max(abs(dpwr-auxpwr));
     if (maxdp > maxdev | devpwr > stdev)
     startp(an)=0;
     endp(an)=0;
     end
     end
   
     startpp = nonzeros(startp);
     endpp = nonzeros(endp);