%spectrum Estimation --ladmm
spectEstimation('data/DRCSI_data_format_v1.mat', 'data/DRCSI_whole_spatmask.mat', ...
'data/DRCSI_dict.mat', 'demos/DRCSI_ladmm.ini', 'Result/DRCSI_inj_mouse_data_ladmm_spect.mat');
%spectrum Estimation --admm
% spectEstimation('Phantom2D/Phantom_data.mat', 'Phantom2D/Phantom_mask.mat', ...
% 'Phantom2D/Phantom_spectrm_info.mat', 'demos/Phantom_admm.ini', 'Result/Phantom_admm_spect.mat')
%
%spectrum Estimation --nnls
% spectEstimation('Phantom2D/Phantom_data.mat', 'Phantom2D/Phantom_mask.mat', ...
% 'Phantom2D/Phantom_spectrm_info.mat', 'demos/Phantom_nnls.ini', 'Result/Phantom_nnls_spect.mat')
%plot Average spectra
plotAvgSpectra('Result/DRCSI_inj_mouse_data_ladmm_spect.mat', ...
'data/DRCSI_whole_spatmask.mat', 'Result/DRCSI_inj_mouse_spectra_ladmm_spectra', 'png','pdf');
% Plot spectroscopic image
idx=1;%encoding sample chosen to show back ground MR data intensity
plotspectIm('Result/DRCSI_inj_mouse_data_ladmm_spect.mat', 'data/DRCSI_data_format_v1.mat',idx, ...
'data/DRCSI_whole_spatmask.mat','Result/DRCSI_inj_mouse_spectroscopic_Im', 'png')
% Plot component Maps
plotCompMaps('Result/DRCSI_inj_mouse_data_ladmm_spect.mat', 'data/spectrum_mask_inj_mouse.mat', ...
'data/four_color.mat', 'Result/DRCSI_inj_mouse_component_maps', 'png', 'epsc');