Skip to content
Snippets Groups Projects
Commit 54dc0b03 authored by Julien Malik's avatar Julien Malik
Browse files

ENH: remove Matlab dir

parent d63ff568
No related branches found
No related tags found
No related merge requests found
function otb_save_hsi_tiff(hsi, path)
otb_save_cube_tiff(hsi.getDataCube(), [path '/hsi_cube.tif']);
otb_save_matrix_Spectras_tiff(hsi.getEndmembers(), [path '/endmembers.tif']);
otb_save_cube_tiff(hsi.getAbundanceCube(), [path '/abundance.tif']);
otb_save_cube_tiff(hsi.getNoiseCube(), [path '/noise.tif']);
return;
function otb_save_matrix_2D_tiff(matrix, filename)
t = Tiff(filename,'w');
tagstruct.ImageLength = size(matrix,1);
tagstruct.ImageWidth = size(matrix,2);
tagstruct.Photometric = Tiff.Photometric.MinIsBlack;
tagstruct.BitsPerSample = 64;
tagstruct.SampleFormat = 3;
tagstruct.SamplesPerPixel = 1;
tagstruct.RowsPerStrip = 1;
tagstruct.PlanarConfiguration = Tiff.PlanarConfiguration.Chunky;
tagstruct.Software = 'MATLAB';
t.setTag(tagstruct);
t.write(matrix);
return;
% each column of matrix is saved as a multiband pixel
function otb_save_matrix_Spectras_tiff(matrix, filename)
t = Tiff(filename,'w');
tagstruct.ImageLength = 1;
tagstruct.ImageWidth = size(matrix,2);
tagstruct.Photometric = Tiff.Photometric.MinIsBlack;
tagstruct.BitsPerSample = 64;
tagstruct.SampleFormat = 3;
tagstruct.SamplesPerPixel = size(matrix,1);
tagstruct.RowsPerStrip = 1;
tagstruct.PlanarConfiguration = Tiff.PlanarConfiguration.Chunky;
tagstruct.Software = 'MATLAB';
t.setTag(tagstruct);
t.write(reshape(matrix',1,size(matrix,2), size(matrix,1)));
return;
% save a cube as a multiband image
function otb_save_cube_tiff(cube, filename)
t = Tiff(filename,'w');
tagstruct.ImageLength = size(cube,1);
tagstruct.ImageWidth = size(cube,2);
tagstruct.Photometric = Tiff.Photometric.MinIsBlack;
tagstruct.BitsPerSample = 64;
tagstruct.SampleFormat = 3;
tagstruct.SamplesPerPixel = size(cube,3);
tagstruct.RowsPerStrip = 1;
tagstruct.PlanarConfiguration = Tiff.PlanarConfiguration.Chunky;
tagstruct.Software = 'MATLAB';
t.setTag(tagstruct);
t.write(cube);
return;
function hsi = otb_synthetic_hsi()
J = 5; % nb endmembers
iota = 0.9; % purity index - 0< iota < 1 with iota~<1
snr = 40; % noise
eta = 0; %parameter to control noise : white noise (eta = 0), color noise (eta > 0)
noise_type = 'additive'; % 'additive' or 'poisson'
% verbose = 0;
nbRows = 100;
nbColumns = 100;
libName = 'usgs';
specLib = class_specLib(libName);
hsi = class_lsmmSynthHsi(J, specLib, iota, snr, eta, noise_type, nbRows, nbColumns);
%disp(hsi.getEndmembers());
%disp(hsi.getSpecLib());
%disp(hsi.getNoiseMatrix());
%disp(hsi.getNbEndmembers());
%disp(hsi.getLibraryIdces());
%disp(hsi.nbColumns());
return;
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment