Read, process and display ESA proposal sample data

read_and_display.m

clear spectral_collage 

cd /media/martin/data/ESA_GANGA/TASI-examples/

%read and display
file = fopen('TASI_2018_07_18_124001.cc.dc.rc2bb.sc.rY.sum12.Norb_g.dat.atmY.24-bands.BBT.dat');
im = fread(file,'float');

im = reshape(im,[1754,1713])';
im(im==0)=315;

subplot(2,3,1)
imagesc(im)
colorbar
title('BBT')

fclose(file);

%read and display
file = fopen('TASI_2018_07_18_124001.cc.dc.rc2bb.sc.rY.sum12.Norb_g.dat.atmY.24-bands.LST.dat');
im = fread(file,'float');

im = reshape(im,[1754,1713])';
im(im==0)=315;

subplot(2,3,2)
imagesc(im)
colorbar
title('LST')

fclose(file);

%read and display
file = fopen('TASI_2018_07_18_124001.cc.dc.rc2bb.sc.rY.sum12.Norb_g.dat.atmY.24-bands.LSE.dat');
im = fread(file,'float');

wavelengths = [8602.250,  8711.750,  8821.250,  8930.750,  9040.250,  9149.750,  9259.250,  9368.750,  9806.750,  9916.250, 10025.750, 10135.250, 10244.750, 10354.250, 10463.750, 10573.250, 10682.750, 10792.250, 10901.750, 11011.250, 11120.750, 11230.250, 11339.750, 11449.250];

im = reshape(im,[1754,1713,24]);
im=permute(im,[2 1 3]);
%im(im==0)=min(im(im>0));

subplot(2,3,3)

divisions=linspace(1,size(im,2),length(wavelengths));
for i = 1:23
    spectral_collage(:,int32(divisions(i):divisions(i+1)))=im(:,int32(divisions(i):divisions(i+1)),i);
end

imagesc(spectral_collage)
colorbar
title(['LSE spectral collage (' int2str(size(im,3)) ' channels)'])

fclose(file);

%read and display
file = fopen('TASI_2018_07_18_124001.cc.dc.rc2bb.sc.rY.sum12.Norb_g.dat.atmY.24-bands.LLL.dat');
im = fread(file,'float');

im = reshape(im,[1754,1713,24]);
im=permute(im,[2 1 3]);
im(im==0)=min(im(im>0));

subplot(2,3,4)

divisions=linspace(1,size(im,2),length(wavelengths));
for i = 1:23
    spectral_collage(:,int32(divisions(i):divisions(i+1)))=im(:,int32(divisions(i):divisions(i+1)),i);
end

imagesc(spectral_collage)
colorbar
title(['LLL spectral collage (' int2str(size(im,3)) ' channels)'])
fclose(file);

im = load('TASI_2018_07_18_124001.cc.dc.rc2bb.sc.rY.sum12.Norb_g.pix.hdr.mat');
im = im.arr;
im(im==0)=min(im(im>0));

subplot(2,3,5)

divisions=floor(linspace(1,size(im,2),size(im,3)));
for i = 1:size(im,3)-1
    spectral_collage(:,int32(divisions(i):divisions(i+1)))=im(:,int32(divisions(i):divisions(i+1)),i);
end

imagesc(spectral_collage)
colorbar
title(['Landsat(?) spectral collage (' int2str(size(im,3)) ' channels)'])

cd /media/martin/data/ESA_GANGA/HyPlant-examples/
file = fopen('20180718-S3IT-1240-3050-L3-S-DUAL_radiance_img_atm_pol_smcorr-rect.bsq');
im=fread(file,'int16');
im=reshape(im,[1566,1421,626]);
im=permute(im,[2 1 3]);
divisions=floor(linspace(1,size(im,2),size(im,3)));
clear spectral_collage
for i = 1:size(im,3)-1
    spectral_collage(:,int32(divisions(i):divisions(i+1)))=im(:,int32(divisions(i):divisions(i+1)),i);
end
subplot(2,3,6)
imagesc(spectral_collage)
title(['HyPlant spectral collage (' int2str(size(im,3)) ' channels)'])
colorbar

read_pix.ipynb