%?RADARSAT-1的數(shù)據(jù)處理之wKA成像1
%?wKA的精確實現(xiàn)
%?2016.10.08

clc;clear?all;close?all;

%%?提取原始數(shù)據(jù)
run?specify_parameters.m
run?extract_data.m
clear?all;

%%?參數(shù)設(shè)置
%?雷達參數(shù)
C?=?2.9979e8;???????????????????????????????????????????????????????????????%?光速
f0?=?5.3e9;?????????????????????????????????????????????????????????????????%?中心頻率
lambda?=?C/f0;??????????????????????????????????????????????????????????????%?波長
%?飛行平臺
V?=?7062;???????????????????????????????????????????????????????????????????%?平臺速度
%?距離向
Kr?=?-0.72135e12;???????????????????????????????????????????????????????????%?調(diào)頻率
Tr?=?41.75e-6;??????????????????????????????????????????????????????????????%?脈沖持續(xù)時間
Fr?=?32.317e6;??????????????????????????????????????????????????????????????%?距離向采樣率
dtau?=?1/Fr;????????????????????????????????????????????????????????????????%?距離向采樣時間間隔
%?方位向
Ka?=?1733;??????????????????????????????????????????????????????????????????%?方位向調(diào)頻率
Fa?=?1256.98;???????????????????????????????????????????????????????????????%?方位向采樣率
fc?=?-6900;?????????????????????????????????????????????????????????????????%?多普勒中心頻率
deta?=?1/Fa;????????????????????????????????????????????????????????????????%?距離向采樣時間間隔
%?其他參數(shù)
t0?=?6.5956e-3;?????????????????????????????????????????????????????????????%?獲取數(shù)據(jù)時間
R0?=?t0*C/2;????????????????????????????????????????????????????????????????%?最短斜距

%%?獲取原始數(shù)據(jù)
oldFolder?=?cd（‘.\scene01‘）;
load?CDdata1.mat
cd（oldFolder）;
data=double（data）;
[NaNr]?=?size（data）;

%%?確定參考距離
Rref?=?（t0+Nr/2*dtau）*C/2;??????????????????????????????????????????????????%?參考距離選在成像區(qū)域中心

%%?數(shù)據(jù)末尾補零
Za?=?800;???????????????????????????????????????????????????????????????????%?方位向補零數(shù)
Zr?=?ceil（Tr/dtau）;?????????????????????????????????????????????????????????%?距離向補零數(shù)
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（‘方位向（采樣點）‘）;

%%?時間軸、頻率軸設(shè)置
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

%%?參考函數(shù)相乘
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