資源簡介
SAR成像算法小例子matlab
代碼片段和文件信息
clear?all;?close?all;?clc;
%?Tutorial:?Focusing?of?SAR?raw?data
%?Programmed?by?Diana?Weihing?&?Stefan?Auer
%?Remote?Sensing?Technology?
%?Technische?Universit鋞?M黱chen
%?2011
%--------------------------------------------------------------------------
%?A.)?Load?raw?data
%?Dummy?image?(‘Test_Image‘)
%raw?=?load(‘Test_Image‘);?%?type:?field
%max_power?=?2.1807e+008;
%?Real?SAR?data?(‘ERS‘)
raw?=?load(‘ERS‘);?%?type:?field
max_power?=?10000;
%?Multilooking?on/off?
multilook?=?1;?%?0:?no;?1:?yes
ml_factor?=?5;?%?multilooking?factor?(spatial)
%?Get?image?size
fn?=?fieldnames(raw);
Raw_data?=?getfield(rawfn{1});
size_azimuth?=?size(Raw_data1);
size_range?=?size(Raw_data2);
%--------------------------------------------------------------------------
%?B.)?Sensor?parameters?(ERS?satellite)
fs=18.962468*10^6;??????????%?Range?Sampling?Frequency?[Hz]??????
K_r=4.18989015*10^(11);?????%?FM?Rate?Range?Chirp?[1/s^2]?-->?up-chirp??
tau_p=37.12*10^(-6);????????%?Chirp?duration?[s]
V=7098.0194;????????????????%?Effective?satellite?velocity?[m/s]
lambda=0.05656;?????????????%?Length?of?carrier?wave?[m]
R_0=852358.15;??????????????%?Range?to?center?of?antenna?footprint?[m]
ta=0.6;?????????????????????%?Aperture?time?[s]??
prf=1679.902;???????????????%?Pulse?Repitition?Frequency?[Hz]
%--------------------------------------------------------------------------?
%?C.)?Display?raw?data
figure;
imagesc(abs(Raw_data));colormap(‘gray‘);
title(‘Raw?data?(unfocused)‘);?xlabel(‘Range?pixels‘);?ylabel(‘Azimuth?pixels‘);
axis?equal;?axis?off;
%--------------------------------------------------------------------------
%?D.)?Define?correlation?chirp?in?range
%?Basic?definitions
range_chirp=zeros(1size_range);?%?empty?vector?to?be?filled?with?chirp?values
tau=-tau_p/2:1/fs:tau_p/2;?%?time?axis?in?range
omega?=?-fs/2:1/tau_p:fs/2;?%?frequency?axis?in?range??
%?Define?chirp?in?range????
ra_chirp_temp?=?exp(1i.*pi.*K_r.*tau.^2);?
%?Get?size?of?chirp
size_chirp_r?=?length(tau);
%--------------------------------------------------------------------------
%?Display?range?chirp?(time?+?frequency?domain)
figure;
subplot(211);?plot(taureal(ra_chirp_temp)‘b‘);?hold?on;
plot(tauimag(ra_chirp_temp)‘r‘);
title(‘Range?chirp‘);?xlabel(‘Time?axis?[sec]‘);?ylabel(‘Amplitude‘);
xlim([min(tau)?max(tau)]);
subplot(212);?plot(omega?abs(fftshift(fft(ra_chirp_temp))));
title(‘Spectrum?of?range?chirp‘);?xlabel(‘Frequency?[Hz]‘);?ylabel(‘Absolute‘);
xlim([min(omega)?max(omega)]);
%--------------------------------------------------------------------------
%?Position?chirp?in?range?vector?(centered)
%?-->?used?for?correlation?procedure
range_chirp(ceil((size_range-size_chirp_r)/2):size_chirp_r+ceil((size_range-size_chirp_r)/2)-1)=ra_chirp_temp;
%?Transform?vector?in?frequency?domain?(fourier?transform)
RANGE_CHIRP=fft(range_chirp);
%?Define?chirp?for?correlation?
%?-->?conjugate?complex?of?
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???67109056??2014-09-25?19:41??Executables\ERS.mat
?????文件???????7086??2016-03-29?10:26??Executables\SAR_focusing_tutorial_solution.m
?????文件???????6386??2014-09-25?19:41??Executables\SAR_focusing_tutorial_temp.m
?????文件???28745364??2014-09-25?19:41??Executables\Test_Image.mat
?????目錄??????????0??2016-03-29?15:19??Executables
-----------?---------??----------?-----??----
?????????????95867892????????????????????5
- 上一篇:matlab深度學習算法合集
- 下一篇:小波分析理論與matlabR2007實現
評論
共有 條評論