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