Import_milk_example

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')