SAR成像算法程序

SAR RD成像算法 二维脉压 附实验数据

clear all
%% 频域脉压&方位匹配
C=3.0e8;                    %光速
c=C;
PRF=1395.0;                 %脉冲重复频率
Fs=1.12492800e+7;           %采样频率
Tp=4.00025602e-5;           %发射脉冲宽度
fc=9.600000000000000e+009;  %发射信号中心频率
gama=-2.3724394e+11;        % Chirp rate
B=-gama*Tp                  %发射带宽
lambda=C/fc;                %中心频率波长
f0=4.7346377e+6;            %采样中频 

u=3.986008*1e5*1e9;
Re=6378.135e3;
h=215e3;
Vs=sqrt（u/（h+Re））;
Vg=Re/（h+Re）*Vs;         %载机速度
v=sqrt（Vs*Vg）;

fdc=1397.0               %提供参考
kr=-9357.0               %提供参考

Rs=-2*v.^2/（kr*lambda）;
fai=acos（（（Re+h）.^2-Rs.^2-Re.^2）/（2*Re*Rs））;        %入射角
beta=acos（（（Re+h）.^2+Rs.^2-Re.^2）/（2*（Re+h）*Rs））;   %视角

Da=12;
azibeamwidth=lambda/Da; %方位波束宽度 
Ta=Rs*azibeamwidth/v;
DeltaF=2*v*sin（azibeamwidth）/lambda;
Band_Ratio=DeltaF/PRF;


tp_num=fix（Tp*Fs）;
tp_nan=fix（Ta*PRF）;


fd1=fopen（‘raws.dat‘‘rb‘）;
s=fread（fd1‘uchar‘）;
fclose（fd1）
s=s（1:2:2097152）+j*s（2:2:2097152）;
s=reshape（s10241024）;
s=s-mean（mean（s））;
[ran_numazi_num]=size（s）;
%% 成像
x=s;
[nrnnan]=size（x）;
tnrn=[-nrn/2:nrn/2-1].‘/Fs;
fr=[-nrn/2:nrn/2-1].‘/nrn*Fs;%距离频率
H_match=exp（-1j*pi/gama*fr.^2）;
% H=exp（-1j*pi/gama*fr.^2）;
H1=exp（1j*4*pi*Rs/c*fr）;
x1=fftshift（fft（fftshift（x1）[]1）1）;
x2=zeros（nrnnan）;
%距离脉压
for n=1:nan
   x2（:n）=x1（:n）.*H_match.*H1;
end
x3=ifftshift（ifft（ifftshift（x21）[]1）1）;
figure;imagesc（abs（x3））
%方位脉压
deltaR=c/2/Fs;
R=Rs+[-nrn/2:nrn/2-1]*deltaR;
faM=2*v/lambda;
fa=[-nan/2:nan/2-1]/nan*PRF;
x4=fftshift（fft（fftshift（x32）[]2）2）;
for nn=1:nrn
%    Ha=exp（1j*2*pi*Rs/v*sqrt（faM.^2-fa.^2）） ;
   Ha=exp（1j*2*pi*R（nn）/v*sqrt（faM.^2-fa.^2）） ;%场景很小，直接用场景中心匹配都不会散焦
   x5（nn:）=x4（nn:）.*Ha;
end
x6=ifftshift（ifft（ifftshift（x52）[]2）2）;
figure;imagesc（abs（x6）[0 600]）

imwrite（abs（x6）/150‘YD_result1.jpg‘‘jpg‘）
imshow（‘YD_result.jpg‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    2097152  1999-05-06 22:03  raws.dat

     文件       1975  2017-06-17 22:01  SAR_image_simple.m

----------- ---------  ---------- -----  ----

              2099127                    2