%?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