SAR面目标成像

SAR的面目标成像 利用matlab实现了 平顶楼 四角锥等三种立体目标的slc仿真

clear all;

Tr=4e-6;                                     
Br=50e6;                                     
Kr=Br/Tr;                                       
Nr=512;
fs=100e6;
dt=1/fs;  
c=3e8;                                    
fc=75e9;                                  
lamda=c/fc; 
Rc=865000;
 

v=7098.2;                                            
Lsar=lamda*Rc/6;                                      
Ta=Lsar/v;
Na=512;
% theta=28/180*pi;
fa=2*v/6*2;
du=1/fa;                                        
 
%%%%%%%%%%%%%%%%%%%%%%%%%载入回波数据1%%%%%%%%%%%%%%%%%%%%%%%%%%%
load planechoh1.mat;
% load zhui1.mat;
% load ping1.mat;
 

tr=-Tr/2:dt:Tr/2;
pr_t=exp（i*2*pi*1/2*Kr*tr.^2）;
pr_t=[zeros（156）pr_t];                  
pr_f=fft（pr_tNr）;
wr=fftshift（hamming（Nr））.‘;
for flaga=1:Na
    echo（flaga:）=fft（fftshift（echo（flaga:）））;
    echo（flaga:）=echo（flaga:）.*conj（pr_f）.*wr;
    echo（flaga:）=ifft（echo（flaga:））;
end
 

Ka=2*v^2/lamda/Rc;
ta=-Ta/2:du:Ta/2;
ta=ta+0.8/6*du;
pa_t=exp（i*2*pi*1/2*Ka*ta.^2+i*0*2*pi*ta）;
pa_t=[zeros（175）pa_t];                
pa_f=fft（pa_tNa）;
wa=fftshift（hamming（Na））;
for flagr=1:Nr
    echo（:flagr）=fft（fftshift（echo（:flagr）））;
    echo（:flagr）=echo（:flagr）.*（pa_f）.‘.*wa;
    echo（:flagr）=ifft（echo（:flagr））;
end
% fig1=echo;
 
% a1=echo（55:105180:280）;
% figure
% imagesc（abs（a1））;

I=mat2gray（abs（echo））;
S=medfilt2（I）;
figureimshow（S）

save slc1.mat echo 


%%%%%%%%%%%%%%%%%%%%%%%%% 第二组成像  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%载入回波数据2%%%%%%%%%%%%%%%%%%%%%%%%%%%
load planechoh2.mat;
% load zhui2.mat;
% load ping2.mat;
 

tr=-Tr/2:dt:Tr/2;
pr_t=exp（i*2*pi*1/2*Kr*tr.^2）;
pr_t=[zeros（156）pr_t];                  
pr_f=fft（pr_tNr）;
wr=fftshift（hamming（Nr））.‘;
for flaga=1:Na
    echo（flaga:）=fft（fftshift（echo（flaga:）））;
    echo（flaga:）=echo（flaga:）.*conj（pr_f）.*wr;
    echo（flaga:）=ifft（echo（flaga:））;
end
 

Ka=2*v^2/lamda/Rc;
ta=-Ta/2:du:Ta/2;
ta=ta+0.8/6*du;
pa_t=exp（i*2*pi*1/2*Ka*ta.^2+i*0*2*pi*ta）;
pa_t=[zeros（175）pa_t];               
pa_f=fft（pa_tNa）;
wa=fftshift（hanning（Na））;
for flagr=1:Nr
    echo（:flagr）=fft（fftshift（echo（:flagr）））;
    echo（:flagr）=echo（:flagr）.*（pa_f）.‘.*wa;
    echo（:flagr）=ifft（echo（:flagr））;
end


I=mat2gray（abs（echo））;
S=medfilt2（I）;
figureimshow（S）
save slc2.mat echo

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

     文件    2106836  2012-06-19 00:33  成像\ping1.mat

     文件    1908754  2012-06-19 00:34  成像\ping2.mat

     文件    3286893  2012-06-15 20:34  成像\planechoh1.mat

     文件    2989756  2012-06-15 20:38  成像\planechoh2.mat

     文件       2430  2012-06-30 16:17  成像\slc.m

     文件    2106923  2012-06-19 01:13  成像\zhui1.mat

     文件    1923488  2012-06-19 01:14  成像\zhui2.mat

     目录          0  2013-03-11 11:03  成像

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

             14325080                    8