SAR成像算法小matlab

SAR成像算法小例子matlab

clear all; close all; clc;

% Tutorial: Focusing of SAR raw data
% Programmed by Diana Weihing & Stefan Auer
% Remote Sensing Technology 
% Technische Universit鋞 M黱chen
% 2011

%--------------------------------------------------------------------------

% A.） Load raw data

% Dummy image （‘Test_Image‘）
%raw = load（‘Test_Image‘）; % type: field
%max_power = 2.1807e+008;

% Real SAR data （‘ERS‘）
raw = load（‘ERS‘）; % type: field
max_power = 10000;

% Multilooking on/off?
multilook = 1; % 0: no; 1: yes
ml_factor = 5; % multilooking factor （spatial）

% Get image size
fn = fieldnames（raw）;
Raw_data = getfield（rawfn{1}）;
size_azimuth = size（Raw_data1）;
size_range = size（Raw_data2）;
%--------------------------------------------------------------------------

% B.） Sensor parameters （ERS satellite）
fs=18.962468*10^6;          % Range Sampling Frequency [Hz]      
K_r=4.18989015*10^（11）;     % FM Rate Range Chirp [1/s^2] --> up-chirp  
tau_p=37.12*10^（-6）;        % Chirp duration [s]
V=7098.0194;                % Effective satellite velocity [m/s]
lambda=0.05656;             % Length of carrier wave [m]
R_0=852358.15;              % Range to center of antenna footprint [m]
ta=0.6;                     % Aperture time [s]  
prf=1679.902;               % Pulse Repitition Frequency [Hz]

%-------------------------------------------------------------------------- 

% C.） Display raw data
figure;
imagesc（abs（Raw_data））;colormap（‘gray‘）;
title（‘Raw data （unfocused）‘）; xlabel（‘Range pixels‘）; ylabel（‘Azimuth pixels‘）;
axis equal; axis off;

%--------------------------------------------------------------------------

% D.） Define correlation chirp in range

% Basic definitions
range_chirp=zeros（1size_range）; % empty vector to be filled with chirp values
tau=-tau_p/2:1/fs:tau_p/2; % time axis in range
omega = -fs/2:1/tau_p:fs/2; % frequency axis in range ?

% Define chirp in range    
ra_chirp_temp = exp（1i.*pi.*K_r.*tau.^2）; 

% Get size of chirp
size_chirp_r = length（tau）;

%--------------------------------------------------------------------------

% Display range chirp （time + frequency domain）
figure;
subplot（211）; plot（taureal（ra_chirp_temp）‘b‘）; hold on;
plot（tauimag（ra_chirp_temp）‘r‘）;
title（‘Range chirp‘）; xlabel（‘Time axis [sec]‘）; ylabel（‘Amplitude‘）;
xlim（[min（tau） max（tau）]）;
subplot（212）; plot（omega abs（fftshift（fft（ra_chirp_temp））））;
title（‘Spectrum of range chirp‘）; xlabel（‘Frequency [Hz]‘）; ylabel（‘Absolute‘）;
xlim（[min（omega） max（omega）]）;

%--------------------------------------------------------------------------

% Position chirp in range vector （centered）
% --> used for correlation procedure
range_chirp（ceil（（size_range-size_chirp_r）/2）:size_chirp_r+ceil（（size_range-size_chirp_r）/2）-1）=ra_chirp_temp;

% Transform vector in frequency domain （fourier transform）
RANGE_CHIRP=fft（range_chirp）;

% Define chirp for correlation 
% --> conjugate complex of 

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

     文件   67109056  2014-09-25 19:41  Executables\ERS.mat

     文件       7086  2016-03-29 10:26  Executables\SAR_focusing_tutorial_solution.m

     文件       6386  2014-09-25 19:41  Executables\SAR_focusing_tutorial_temp.m

     文件   28745364  2014-09-25 19:41  Executables\Test_Image.mat

     目录          0  2016-03-29 15:19  Executables

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

             95867892                    5