RADARSAT_CSA.rar

RADARSAT_CSA.rar

% RADARSAT-1的数据处理之CSA成像
% 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;                                                  % 参考距离
Vref = V;                                                                   % 参考速度
fref = fc;                                                                  % 参考频率

%% 数据末尾补零
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;                                            % 方位向频率轴

%% 中间变量设置
D_feta_V = sqrt（1-C^2*feta.^2/（4*V^2*f0^2））;                                % D（fetaV）
D_feta_Vref = sqrt（1-C^2*feta.^2/（4*Vref^2*f0^2））;                          % D（fetaVref）
D_fref_V = sqrt（1-C^2*fref^2/（4*V^2*f0^2））;                                 % D（frefV）
D_fref_Vref = sqrt（1-C^2*fref^2/（4*Vref^2*f0^2））;                           % D（fre

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

     文件       6284  2016-10-13 16:08  RADARSAT_CSA.m

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

                 6284                    1