clear
clc
ds = tabularTextDatastore('D:\milk_example','FileExtensions','.txt');
[size_ds,~] = size(ds.Files);
for i = 1: size_ds
%% Import data from text file
%% Set up the Import Options and import the data
opts = delimitedTextImportOptions("NumVariables", 4);
% Specify range and delimiter
opts.DataLines = [2, Inf];
opts.Delimiter = ",";
% Specify column names and types
opts.VariableNames = ["Var1", "VarName2", "Var3", "VarName4"];
opts.SelectedVariableNames = ["VarName2", "VarName4"];
opts.VariableTypes = ["string", "double", "string", "double"];
% Specify file level properties
opts.ExtraColumnsRule = "ignore";
opts.EmptyLineRule = "read";
% Specify variable properties
opts = setvaropts(opts, ["Var1", "Var3"], "WhitespaceRule", "preserve");
opts = setvaropts(opts, ["Var1", "Var3"], "EmptyFieldRule", "auto");
% Import the data
asciispec = readtable(ds.Files{i,1}, opts);
%% Convert to output type
asciispec = table2array(asciispec);
tempt_data = asciispec;
tempt_data(find(tempt_data(:,2)<-0.5001),:) =[];
tempt_data(find(tempt_data(:,2)>9.5001),:) =[];
ppm_1 = tempt_data(1:2:end,2);
table_1(:,i) = tempt_data(1:2:end,1);
%% Clear temporary variables
clear opts
end
final_table = [ppm_1, table_1];%with both ppm and intensity
save('final_table','final_table')
save('table_1','table_1')
save('ppm_1','ppm_1')
names = char(ds.Files);
dlmwrite('sample_names.txt',names,'delimiter',''); %samples_names
figure (1)
plot (ppm_1,table_1)
set(gca,'XDir','reverse')