noisesect.m

rang doppler imaging and motion compensation中的源代码

function noisemat = noisesect(noisesig,sect,desiren)
% 
% Select the region composed by noise in spatial domain.
%

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


    noisemat = [];
    numrow=m/sect;
    for it = 1:n
    reshapesig=reshape(noisesig(:,it),numrow,sect);
    meanres=mean(reshapesig);
    [smean,indmean]=sort(meanres);
    findmean=find(indmean<=desiren);

    noisevect=[];
    for itt=1:desiren
    noisevect=[noisevect;reshapesig(:,indmean(findmean(itt)))];
    end
    noisemat = [noisemat,noisevect];
    end