rma.m

条带式合成孔径雷达的cs。RDA。RMA算法 作者何志华

%%Range Migration Algorithm,RMA
%NUDT,Hezhihua
%E-mail:[e-mail]
%last update:24/6/2005

load sarecho.mat sxc_kxky kx ky kxc Nx Ny

%%Stolt插值
P=8/2;%8点sinc插值
Stolt=zeros(Nx,Ny);%初始化
dkx=kx(2)-kx(1);%频域间隔
KX_kxky=(sqrt((kx+kxc).^2'*ones(1,Ny)+ones(Nx,1)*ky.^2)-kxc)-kx'*ones(1,Ny);
KX=floor(KX_kxky/dkx);

for j=1:Ny
    for i=P+max(max(KX)):Nx-max(max(KX))-P
       T=KX(i,j)+(-P+1:P);
       f=sxc_kxky(i+T,j).';
       Stolt(i,j)=sum(f.*sinc(KX_kxky(i,j)/dkx-T));
    end
    j
end

subplot(221)
plot_img(ky,kx,sxc_kxky,35),grid on
xlabel('ky domain(cross-range),m^-^1')
ylabel('kx domain(range),m^-^1')
title('|s_x_c(k_x,k_y)|')

subplot(222)
plot_img(ky,kx,Stolt,35),grid on
xlabel('ky domain(cross-range),m^-^1')
ylabel('kx domain(range),m^-^1')
title('|s_x_c(k_x,k_y)|')