wtv02img.m

rang doppler imaging and motion compensation中的源代码

%		WTV02IMG.M
% this program uses the adaptive wavelet time-freq decomposition to generate 
% the ISAR image. it generates the wavelet coefficients B0 thru B5
%
figure
subplot(1,1,1)
% +++++++
% +++++  Find b0, the inner product between s0 and g0 		++++++
% +++++  INNER PRODUCT is considered as the orthogonal projection ++++
% +++++++
xpos1=3;
ypos1=5;

b0 = g0_256*th_256;

f=(0:255)'*(256*ts);		
b0_img256=abs(b0)/max(max(abs(b0)));   	% normalize
b0_img256=fftshift(b0_img256);         	% center image

% 			   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%++++-----------------------------
%+++ Find  INNER RPODUCT between g1 and  +++
%++++

b1 = g1_256*th_256;
b1_img256=abs(b1)/max(max(abs(b1)));   	% normalize
b1_img256=fftshift(b1_img256);         	% center image

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%++++-----------------------------
%+++ Find  INNER RPODUCT between g1 and  +++
%++++

b1 = g2_256*th_256;
b2_img256=abs(b1)/max(max(abs(b1)));   	% normalize
b2_img256=fftshift(b2_img256);         	% center image

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%++++-----------------------------
%+++ Find  INNER RPODUCT between g1 and  +++
%++++

b1 = g3_256*th_256;
b3_img256=abs(b1)/max(max(abs(b1)));   	% normalize
b3_img256=fftshift(b3_img256);         	% center image

subplot(2,2,4)
X=0:42/63:42;
Y=0:48/63:48;
clevel=[-3,3,0];
% contour(X,Y,10*log10(b0_img256),clevel);
colormap(gray(128))
b0_img256 = (1-b0_img256);
image(X,Y,(128*b0_img256));
xlabel('ws (m)')
ylabel('wc (m)')
grid
tps = sprintf('t = %5.3f sec',tp4)
text(xpos1,ypos1,tps)

subplot(2,2,3)
clevel=[-3,3,0];
% contour(X,Y,10*log10(b1_img256),clevel);
colormap(gray(128))
b1_img256 = (1-b1_img256);
image(X,Y,(128*b1_img256));
xlabel('ws (m)')
ylabel('wc (m)')
tps = sprintf('t = %5.3f sec',tp3)
text(xpos1,ypos1,tps)
grid

subplot(2,2,2)
clevel=[-3,3,0];
% contour(X,Y,10*log10(b2_img256),clevel);
colormap(gray(128))
b2_img256 = (1-b2_img256);
image(X,Y,(128*b2_img256));
xlabel('ws (m)')
ylabel('wc (m)')
tps = sprintf('t = %5.3f sec',tp2)
text(xpos1,ypos1,tps)
grid
title('GWT approach')

subplot(2,2,1)
clevel=[-3,3,0];
% contour(X,Y,10*log10(b3_img256),clevel);
colormap(gray(128))
b3_img256 = (1-b3_img256);
image(X,Y,(128*b3_img256));
xlabel('ws (m)')
ylabel('wc (m)')
tps = sprintf('t = %5.3f sec',tp1)
text(xpos1,ypos1,tps)
grid
etime(clock,tiempo0)
%keyboard