RADARSAT_wKA.rar

RADARSAT_wKA.rar

% RADARSAT-1的数据处理之wKA成像1
% wKA的精确实现
% 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

%% 参考函数相乘
Sff = fftshift（fft2（fftshift（S0）））;                                         % 二维FFT
clear S0
WindowR = ones（Na1）*kaiser（Nr2.5）‘;                                       % 距离窗
WindowA = kaiser（Na2.5）*ones（1Nr）;                                        % 方位窗
Sff = Sff.*WindowR.*WindowA;                                                % 加窗
clear WindowR WindowA
Hrfm = exp（1i*

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

     文件       5618  2016-10-13 22:32  RADARSAT_wKA_1.m

     文件       4885  2016-10-12 18:42  RADARSAT_wKA_2.m

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

                10503                    2