SAR成像算法小例子matlab
代码片段和文件信息
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 67109056 2014-09-25 19:41 ExecutablesERS.mat
文件 7086 2016-03-29 10:26 ExecutablesSAR_focusing_tutorial_solution.m
文件 6386 2014-09-25 19:41 ExecutablesSAR_focusing_tutorial_temp.m
文件 28745364 2014-09-25 19:41 ExecutablesTest_Image.mat
目录 0 2016-03-29 15:19 Executables
----------- --------- ---------- ----- ----
95867892 5
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 ExecutablesERS.mat
文件 7086 2016-03-29 10:26 ExecutablesSAR_focusing_tutorial_solution.m
文件 6386 2014-09-25 19:41 ExecutablesSAR_focusing_tutorial_temp.m
文件 28745364 2014-09-25 19:41 ExecutablesTest_Image.mat
目录 0 2016-03-29 15:19 Executables
----------- --------- ---------- ----- ----
95867892 5
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件举报,一经查实,本站将立刻删除。
评论列表(条)