radarsat-1 RD成像算法

%% RADARSAT_RDA_Imaging.m
% 用RD算法对RADARSAT-1的回波数据成像
% author: 四院五队 王建（07042185）
% 为了节省程序的内存开销，代码中包含了较多的clear语句
% 为了减少程序运行时的内存开销峰值，长表达式分多步实现

%% 抽取数据

% 在文件specify_parameters.m中指定相关参数
specify_parameters;
% 从文件中抽取数据
extract_data;
% 对数据块1进行成像
AzimuthBlock=1;
% 加载数据
load（strcat（output_path‘CDdata‘num2str（AzimuthBlock）））;
file_pre = strcat（ output_path output_prefix ‘_‘ num2str（AzimuthBlock） ）;
% 得到自动增益控制（AGC）值
AGC_values = load_AGC_block（ file_pre first_rg_line ...
    Nrg_lines_blk AzimuthBlock  UseMATfiles ）;
% AGC还原，得到真实的回波数据
org_data = load_DATA_block（ file_pre output_path Nrg_lines_blk ...
    Nrg_cells AGC_values AzimuthBlock UseMATfiles ）;
% 以下所有的clear语句都是为了节省内存开销
clear AGC_values

%% 进一步明确RADARSAT-1的相关参数

length_replica  =  2880;         % Total length （I&Q） of replica record
tot_Nrg_cells   =  9288;         % Total number of range cells per line
tot_Nrg_lines   = 19432;         % Total number of range lines （records）
first_replica   =     7;         % First record that contains the replica
PRF             = 1256.98;       % Pulse Reputation Frequency （Hz）
Fa              = PRF;           % 方位向采样频率
Fr              = 32.317e+6;     % Radar sampling rate （Hz）
f0              = 5.300e+9;      % Radar center frequency （Hz）
c               = 2.9979e+8;     % Speed of light （m/s）
t0              = 0.0065956;     % data window start time
R0              = 0.0065956*c/2; % Slant range of first radar sample （m）
Nrepl           = 1349;          % No. of valid samples in the replica
Kr              = 0.72135e+12;   % FM rate of radar pulse （Hz/s）
Tr              = 41.75e-6;      % Chirp duration （s）
Vr              = 7062;          % Effective radar velocity （m/s）
WaveLength      = c/f0;          % Radar wavelength （m）
Ka              = 1733;          % Azimuth FM rate （Hz/s）
f_dc            = -6900;         % Doppler centroid （Hz）
Lr              = 15;            % Range Aperture
La              = 1.5;           % Azimuth aperture
Na_mf           = 705;           % 方位向上匹配滤波器长度
Ta              = Na_mf/Fa;      % the duration of the azimuth match filter

%% 由RADARSAT-1的参数进一步得到所需参数

% 得到待成像的回波数据