RADARSAT_RDA.rar

RADARSAT_RDA.rar

% RADARSAT-1的数据处理之RDA成像SRC2
% 精确实现二次距离压缩的RDA
% 2016.10.08

clc;clear all;close all;

%% 提取原始数据
run specify_parameters.m
run extract_data.m
clear all;

%% 参数设置
% 雷达参数
C = 2.9979e8;                                                               % 光速
f0 = 5.3e9;                                                                 % 中心频率
lambda = C/f0;                                                              % 波长
% 飞行平台
V = 7062;                                                                   % 平台速度
% 距离向
Kr = -0.72135e12;                                                           % 调频率
Tr = 41.75e-6;                                                              % 脉冲持续时间
Fr = 32.317e6;                                                              % 距离向采样率
dtau = 1/Fr;                                                                % 距离向采样时间间隔
% 方位向
Ka = 1733;                                                                  % 方位向调频率
Fa = 1256.98;                                                               % 方位向采样率
fc = -6900;                                                                 % 多普勒中心频率
deta = 1/Fa;                                                                % 距离向采样时间间隔
% 其他参数
t0 = 6.5956e-3;                                                             % 获取数据时间
R0 = t0*C/2;                                                                % 最短斜距

%% 获取原始数据
oldFolder = cd（‘.\scene01‘）;
load CDdata1.mat
cd（oldFolder）;
data=double（data）;
[NaNr] = size（data）;

%% 确定参考距离和参考时间
Rref = （t0+Nr/2*dtau）*C/2;                                                  % 参考距离选在成像区域中心

%% 数据末尾补零
Za = 800;                                                                   % 方位向补零数
Zr = ceil（Tr/dtau）;                                                         % 距离向补零数
data = cat（2datazeros（NaZr））;                                            % 距离向补零
data = cat（1zeros（ZaNr+Zr）data）;                                         % 方位向补零
Na = Na+Za;
Nr = Nr+Zr;

figureimagesc（abs（data））;axis image;set（gcf‘Color‘‘w‘）;
title（‘时域：补零后的原始信号幅度‘）;
xlabel（‘距离向（采样点）‘）;ylabel（‘方位向（采样点）‘）;

%% 时间轴、频率轴设置
tau = t0+（0:Nr-1）*dtau;                                                     % 距离向时间轴
ftau = （（0:Nr-1）-Nr/2）/Nr*Fr;                                               % 距离向频率轴
eta = （（0:Na-1）-Na/2）*deta;                                                 % 方位向时间轴
feta = fc+（（0:Na-1）-Na/2）/Na*Fa;                                            % 方位向频率轴

%% 对齐到方位向多普勒中心
S0 = data.*exp（-1i*2*pi*fc*（eta‘*ones（1Nr）））;                              % 搬移至多普勒频率中心
clear data

%% RC
Tref = t0+Nr/2*dtau;                                                        % 距离向参考时间
ht_rc = （abs（（ones（Na1）*tau）-Tref）Hf_rc = conj（fftshift（fft（fftshift（ht_rc2）[]2）2））;                      % 滤波器的距离向FFT变换
clear ht_rc
S0_rc = fftshift（fft（fftshift（S02）[]2）2）;                               % 信号

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

     文件       5940  2016-10-12 18:41  RADARSAT_RDA_3.m

     文件       6692  2016-10-17 19:56  RADARSAT_RDA_2.m

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

                12632                    2