rdisplay.m

rang doppler imaging and motion compensation中的源代码

 function imaged = rdisplay(disp,data,sflag,dsize,window,resol,beta)
 %
 % Displays the image in various ways.
 % disp is either 'display', 'graydisp', 'logdisp' or 'none'.
 % sflag is either 'imagery', 'rangeprofile' or 'crossrange'.
 % window is either 'rect', 'hamming', 'hanning', 'SVA' or 'kaiser'.
 % dsize is size of matrix to display.  The matrix is zero-padded to meet the 
 % required dimension.  dsize(1)=number of rows.  dsize(2)=number of columns.
 %
 % display: Linear amplitude display.
 % graydisp: Display in gray.
 % logdisp: Log scale display.
 % imagery: Displays the image.
 % rangeprofile: Displays the magnitude of a range profile.
 % crossrange: Displays the magnitude of a cross-range profile.
 % 'rect': Rectangular window.
 % 'hamming' : Hamming window.
 % 'hanning': Han window.
 % 'SVA': Spatially variant apodization.
 %
 
   [m n] = size(data);
   x=-n/2*resol(1):resol(1):(n/2-1)*resol(1);
   y=-m/2*resol(2):resol(2):(m/2-1)*resol(2);


   if (strcmp(window,'rect')  == 1 | strcmp(window,'SVA')  == 1)
   win = ones(m,n);
   end

   if strcmp(window,'hamming')  == 1
   winm= hamming(m);
   winn= hamming(n)';
   win = winm*winn;
   end

   if strcmp(window,'hanning')  == 1;
   winm= hanning(m);
   winn= hanning(n)';
   win = winm*winn;
   end

   if strcmp(window,'kaiser')  == 1
   win = kaiser2(m,n,beta);
   end

   windata=data.*win;

   if strcmp(sflag,'none') ~= 1
   figure
   end

   if strcmp(sflag,'imagery') == 1
   imagep=fft(ifft(windata',dsize(2)).',dsize(1));
   imaged=fftshift(imagep);
   if strcmp(window,'SVA')  == 1
   imaged = SVA2(imaged);
   end
   end

   if strcmp(sflag,'rangeprofile') == 1
   rangep=fftshift2(ifft(windata',dsize(2)));
   imaged=abs(rangep.');
   end

   if strcmp(sflag,'crossrange') == 1
   crossrangep=fftshift2(fft(windata,dsize(1)));
   imaged=abs(crossrangep);
   end

   if strcmp(disp,'display') == 1
   imagesc(abs(imaged))
   end

   if strcmp(disp,'graydisp') == 1
   numimag=abs(imaged);
   maxnum=max(max(numimag));
   normnum=numimag/maxnum;
   invimag=-normnum+1;
   imagesc(x,y,invimag); colormap(gray)
   end

   if strcmp(disp,'contdisp') == 1
   numimag=contour(abs(imaged));
   end
   
   if strcmp(disp,'logdisp') == 1
   numimag=abs(imaged);
   maxnum=max(max(numimag));
   normnum=numimag/maxnum;
   invimag=-normnum+1;
   logimag=10.^(invimag);
 %  invimag=-normnum+2;
 %  logimag=10*log10(invimag);

   imagesc(logimag); colormap(gray)
   end