%%Range Migration Algorithm,RMA
%NUDT,Hezhihua
%E-mail:skynismile@yahoo.com.cn
%last update:24/6/2005

%%圖形顯示
clear all
load sarecho.mat  kx ky x y Xc

figure('name','Range Migration Algorithm','units','normal',...
          'NumberTitle','off','position',[0.2,0.2,0.7,0.6])
load sarecho.mat sxc_kxky;
subplot(221)
plot_img(ky,kx,sxc_kxky,40),grid on,
axis([-1,1,-1,1])
xlabel('ky domain(cross-range),m^-^1')
ylabel('kx domain(range),m^-^1')
title('|s_x_c(k_x,k_y)|')
clear sxc_kxky
drawnow 

subplot(222)
load sarecho.mat Stolt_kxky;
plot_img(ky,kx,Stolt_kxky,40),grid on
axis([-1,1,-1,1])
xlabel('ky domain(cross-range),m^-^1')
ylabel('kx domain(range),m^-^1')
title('|Stolt(k_x,k_y)|')
clear Stolt_kxky
drawnow 

subplot(223)
load sarecho.mat Stolt_xky;
plot_img(ky,x-Xc,Stolt_xky,35),grid on
axis([-1,1,-20,70])
set(gca,'Ytick',[0,50])
xlabel('ky domain(cross-range),m^-^1')
ylabel('x domain(range),m')
title('|Stolt(x,k_y)|')
clear Stolt_xky
drawnow 

subplot(224)
load sarecho.mat Stolt_xy;
plot_img(y,x-Xc,Stolt_xy,50),grid on
axis([-100,100,-20,70])
set(gca,'Xtick',[-100,-50,0,50,100])
set(gca,'Ytick',[0,50])
xlabel('y domain(cross-range),m')
ylabel('x domain(range),m')
title('|Stolt(x,y)|')
clear Stolt_xy
drawnow