BTDCalculator.asv

%% This script calculates the brightness temperature difference of the 
% bands of Himawari-8 Data. 
% Band 13 minus Band 15 (10.4 um - 12.4um)

clear

addpath '/Users/denny/OneDrive - Nanyang Technological University/Y4/FYP/H8_Codes'

%------------------------Config to change-------------

Volcano = 'Sinabung';
YYYYMM = '201905';
DD = {'12','13','14','15','16',...
    '17','18','19','20','21','22','23','24','25','26','27','28','29','30'};

DayNight = 'Night';

%------------------------------------------------------

Matfile_Data_Folder = '/Users/denny/OneDrive - Nanyang Technological University/Y4/FYP/H8_Processed_Data/';

for i = 1:length(DD)

foldername = [Volcano,'_',YYYYMM,'/',Volcano,'_',YYYYMM,DD{i},'_',DayNight];
matfilename = [Volcano,'_',YYYYMM,DD{i},'_',DayNight,'.mat'];
%matfilename = [Volcano,'_',YYYYMM,DD,'_',DayNight,'_Median.mat'];
%stackfilename = [Volcano,'_',YYYYMM,DD,'_',DayNight,'_Stacked.mat'];

matfiletoread = ([Matfile_Data_Folder,'/',foldername,'/',matfilename]);
%stackfiletoread = ([Matfile_Data_Folder,'/',foldername,'/',stackfilename]);

load(matfiletoread)
%load(stackfiletoread)

cd([Matfile_Data_Folder,foldername])

%%

BTD_13_15 = BTD_calculator_function(tbb_13,tbb_15);

BTD_7_13 = BTD_calculator_function(tbb_07,tbb_13);

BTD_13_14 = BTD_calculator_function(tbb_13,tbb_14);

BTD_14_15 = BTD_calculator_function(tbb_14,tbb_15);

ratio_14_15 = BTD_ratio_function(tbb_14,tbb_15);

BTD_7_13 = BTD_calculator_function(tbb_07,tbb_14);

% tbb_14_15_ratio = BTD_ratio_function


%% creating median file name. 
BTDfilename = [Volcano,'_',YYYYMM,DD{i},'_',DayNight,'_BTD.mat'];


% specify which variables to be saved depending on what is to be read.
save([Matfile_Data_Folder,'/',foldername,'/',BTDfilename],'BTD_13_15',...
    'BTD_7_13','BTD_13_14','BTD_14_15','ratio_14_15')

end