radarsat.rar_RADARSAT_Radarsat-2_Radarsat1_radarsat-1_radarsat成像

clear; length_replica = 2880; % Total length (I&Q) of replica record tot_Nrg_cells = 9288; % Total number of range cells per line tot_Nrg_lines = 19432; % Total number of range lines (records) first_replica = 7; % First record that contains the replica PRF = 1256.98; % Pulse Reputation Frequency (Hz) Fr = 32.317e+6; % Radar sampling rate (Hz) f0 = 5.300e+9; % Radar center frequency (Hz) c = 2.9979e+8; % Speed of light (m/s) R0 = 0.0065956*c/2; % Slant range of first radar sample (m) Nrepl = 1349; % No. of valid samples in the replica Kr = 0.72135e+12; % FM rate of radar pulse (Hz/s) Tr = 41.75e-6; % Chirp duration (s) B=30.11e6; %数据格式 num_aux = 50; % RSI - D4 Pg.32 num_header = 192; % RSI - D4 Pg.32 file_header_length = 16252; % RSI - D4 Pg.103 h = num_header; % 192 ha = num_header + num_aux; % 242 num_pixel_data = 2*tot_Nrg_cells; hpa = num_header + num_aux + num_pixel_data; % 18818 numlength=8*hpa ; length_replica = 2880; % Total length (I&Q) of replica record rep_block_length = length_replica + 8*hpa ; % 18576 fid=fopen('E:\合成孔怪雷达数据\DAT_01.001','r'); fp_w=fopen('E:\合成孔怪雷达数据\radasat1data.raw','a+'); fseek(fid,16252,'cof'); data1=zeros(hpa,8); data2=zeros(18576,8); for k=1:2429 data1=fread(fid,[hpa,8],'uint8'); fseek(fid,length_replica,'cof'); data2=data1(243:hpa,1:8); fwrite(fp_w,data2,'uint8'); end fclose('all') nr=tot_Nrg_cells ; % nr_b=tot_Nrg_cells-Nrepl ; % nr_s=(nr-nr_b)/2+1; %%%%%%%%%%%距离门起始%%%%%%%%%%%%%%%%%%% % nr_e=nr_s+nr_b-1; %%%%%%%%%%%%距离门终止%%%%%%%%%%%%%%%%%%% dt=1/Fr; nr_b=8192; nr_s=(nr-nr_b)/2+1; %%%%%%%%%%%距离门起始%%%%%%%%%%%%%%%%%%% nr_e=nr_s+nr_b-1; %%%%%%%%%%%距离门终止%%%%%%%%%%%%%%%%%%% fd=1733; fdc=-6900; lambda=c/f0; vel=7062; % vel = 7.8509e+003; r=R0+((1:nr))*c/2.0/Fr; fr=-2*vel*vel/lambda./r; Ntr=floor(Tr*Fr); transig(1:Ntr) = exp(sqrt(-1)*pi*Kr*(dt*(1:Ntr)-Tr/2).^2); transig=conj(fft(transig,nr)); f=[0:nr/2-1 -nr/2+1:0]*Fr/nr; %%%%%%%%%%%设置频率轴%%%%%%%%%%%%%%%%%%% hr_f=exp(j*pi*f.^2/Kr); fp=fopen('E:\合成孔怪雷达数据\radasat1data.raw','r'); na=tot_Nrg_lines; na=16384; % na=8192; Sr=zeros(nr,na,'single'); fa=fdc+fd*[-na/2+1:0 0:na/2-1 ]/PRF; lambda=c/f0; fseek(fp,9288*2*8,'cof'); delt_beta = sqrt(1-(lambda/2/vel*fa).^2)-1; Sr=zeros(nr_b,na,'single'); for h=1:na data=fread(fp,[9288*2 1],'uint8'); % fseek(fp,9288*2*7,'cof'); % Gr= 1/Kr-r*2*lambda.*(fa(h)*r.^2-1)/delt_beta(h)/c^2; % Hr=exp(-j*pi*Gr.*f.^2+exp(r(h)*2*lambda/delt_beta(h)/c.*f)); % fseek(fp,(subsamp-1)*pulse_leng,'cof'); %%%%%%%%%%降采样处理%%%%%%%%%%%%%%%%%% % Sr_temp(1:1:nr)=data(243:2:end,1)+i*data(244:2:end,1); Sr_temp(1:nr)=(data(1:2:end))+j*(data(2:2:end)); Sr_temp=ifft(fft(Sr_temp-mean(Sr_temp)).*hr_f) ; %%%%%%%%%%%距离向压缩%%%%%%%%%%%%%%%%%% % Sr_temp=ifft(fft(Sr_temp).*transig) ; Sr(:,h)=single(Sr_temp(nr_s:nr_e).'); %%%%%%%%%%%%%%%截取距离门点数%%%%%%%%%% end % fa=fdc+[0:na/2-1 -na/2+1:0]*PRF/na; fa=fdc+fd*[0:na/2-1 -na/2+1:0]/PRF; % f=[0:Nfast/2-1 -Nfast/2+1:0]*Fsr/Nfast; % fa=[0:na/2-1 -na2+1:0]*prf/na; % R=4; % kkk=1:na; % Up00=zeros(size(1:na*R+35000)); % for n=1:nr % Up00=0*Up00; % dR=round(pi*lambda*f(1:nr).^2/c/Kr*Fr); % Up00=0*Up00; % Up00(1:na*R)=single(interp(Sr(n,:),R)); % Sr(n,:)=single(Up00((kkk-1)*R+dR(n)+1)); % end % phase2 = exp((j00*4*pi/lamda)*(delt_beta'*rr)); % t=[0:na/2-1 -na/2+1:0]/PRF % f=fd*t.^2+fdc*t; % f=exp(-j*2*pi*f); % f=fft(f); max=0; for h=1:nr_b ha_f=exp(j*pi*fa.^2/fr(h)); phase2 = exp((j*4*pi/lambda)*(delt_beta'*r(h))); % Sr(h,:)=single(abs(ifft(fft(Sr(h,:)).*ha_f))); %%%%%%%%%%%%方位压缩%%%%%%%%%%%%%%%% % Sr(h,:)=single(abs(ifft(fft(Sr(h,:)).*f))); Sr(h,:)=single(abs(ifft(fft(Sr(h,:)).*phase2.'))); % if max((Sr(h,:)))>=Max % max=max(Sr(h,:)); % end end fid=fopen('e:\radasat1.raw','w'); fwrite(fid,65535*(Sr)/max,'int16') fclose('all')