上一次最后给出了一个成像算法的示例,下面是它对应的代码,里面已经有了初步的注释。
当然这个是个简单情况,但已经大概说明了流程。代码具体的解释,我打算先放一放,毕竟it is a long story,我们需要从长计议。后续章节将结合代码逐渐讲解一下成像的步骤。在合适的时间,我会具体讲解这些代码的意义,以及他忽略的部分。
代码来源已经注明,在此表示感谢!
数据下载地址
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(raw,fn{
1});
size_azimuth = size(Raw_data,1);
size_range = size(Raw_data,2);
%--------------------------------------------------------------------------
% B.) Sensor parameters (ERS satellite)
fs=18.962468*10^6; % Range Sampling Frequency [Hz]
K_r=4.18989015*10</