datareading.m

POLSAR data reader in matlab

%THIS PROGRAM READS DATA OF POLSAR APPLICATION
function I=DataReading()
% clear all;close all;clc
fid=fopen('T11.bin');T11 = fread(fid,'float');
fid=fopen('T12_imag.bin');T12_imag = fread(fid,'float');
fid=fopen('T12_real.bin');T12_real = fread(fid,'float');
fid=fopen('T13_real.bin');T13_real = fread(fid,'float');
fid=fopen('T13_imag.bin');T13_imag = fread(fid,'float');
fid=fopen('T22.bin');T22 = fread(fid,'float');
fid=fopen('T23_imag.bin');T23_imag = fread(fid,'float');
fid=fopen('T23_real.bin');T23_real = fread(fid,'float');
fid=fopen('T33.bin');T33 = fread(fid,'float');
fclose('all');
clear fid;
%==============================================================
%reshape the data in matrix form
temp=size(T11);
Rows=sqrt(temp(1));Cols=Rows;
T11=reshape(T11,Rows,Cols);
T12_imag=reshape(T12_imag,Rows,Cols);
T12_real=reshape(T12_real,Rows,Cols);
T13_real=reshape(T13_real,Rows,Cols);
T13_imag=reshape(T13_imag,Rows,Cols);
T22=reshape(T22,Rows,Cols);
T23_real=reshape(T23_real,Rows,Cols);
T23_imag=reshape(T23_imag,Rows,Cols);
T33=reshape(T33,Rows,Cols);
%==============================================================
%CONVERTING COMPLEX IMAGES TO REAL
T12=complex(T12_real,T12_imag);
T13=complex(T13_real,T13_imag);
T23=complex(T23_real,T23_imag);
%==============================================================
%FORMING COMPOSITE IMAGE
I=T11';
I(:,:,2)=abs(T12');
I(:,:,3)=abs(T13');
I(:,:,4)=T22';
I(:,:,5)=abs(T23');
I(:,:,6)=T33';
%==============================================================
%DISPLAYING IMAGE
% for i=1:6
%     figure,imshow(I(:,:,i))
% end
%plotting the data as RGB image
M(:,:,1)=I(:,:,1);
M(:,:,2)=I(:,:,4);
M(:,:,3)=I(:,:,6);
M=10*log(M+1);
figure,imshow(M);title('Approximated PauliRGB image reconstructed from T11,T22,T33')
end