cia.m

rang doppler imaging and motion compensation中的源代码

function [startp,endp] = CIA(amp,amplen,thre)
%
% Finds a starting and ending points of a segment using
% complex image analysis method.
%

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

     zerovec=zeros(1,amplen);
     onevec=ones(1,amplen);
     meanamp=mean(amp);
     mtt4=thre*meanamp;
     threth=amp > mtt4;
     thretha=double(threth);

     jflag=0;
     startp=[];
     flen=[];
     totlen=n-amplen+1;          
     for al=1:totlen;
     alm=al+amplen-1;
     thrm=thretha(al:alm);
     sthrm=sum(thrm);
     if jflag > 0 & sthrm < amplen
     if al < totlen
     flen=[flen,1];
     jflag=jflag-1;
     end
     if al == totlen
     flen=[flen,ones(1,jflag)];
     sflen=size(flen,2);
     flen=[flen,zeros(1,n-sflen)];
     end
     else
       if (sthrm == amplen & al == 1) | (sthrm == amplen & flen(al-1) == 0)
         if al == totlen
         flen=[flen,onevec];
         startp=[startp,al];
         else
         flen=[flen,1];
         startp=[startp,al];
         jflag=amplen;
         end
           elseif sthrm == amplen & al == totlen & flen(al-1) == 1
           flen(totlen:n)=1;   
             elseif sthrm == amplen & flen(al-1) == 1
             flen=[flen,1];
                elseif al == totlen
                flen=[flen,zerovec];
       else
       flen=[flen,0];
       end
     end
     end
       
     endp=[];
     for am=amplen:n-1
     if flen(am) == 1 & flen(am+1) == 0
     endp=[endp,am];
     end
     end
     if flen(n) == 1
     endp=[endp,n];
     end