clear all; close all; clc
% orig_gen = 'gen';
% type = 'Berea';
types = ["sr"];
% types = ["166", "158"];
orig_gen_all = ["gen", "orig"];
for i = 1:length(types)
for j = 1:length(orig_gen_all)
type = types(i);
orig_gen = orig_gen_all(j);
% ours, berea, ketton
myFolder = sprintf('../Documents/128/%s/', type);
myFiles = dir(sprintf('../Documents/128/%s/%s*.hdf5', type, orig_gen)); %gets all wav files in struct
disp(myFiles);
meanbreadth = zeros(size(myFiles));
volume = zeros(size(myFiles));
surface = zeros(size(myFiles));
euler = zeros(size(myFiles));
for k = 1:length(myFiles)
baseFileName = myFiles(k).name;
fullFileName = fullfile(myFolder, baseFileName);
fprintf(1, 'Now reading %s\n', fullFileName);
obj = h5read(fullFileName, '/data');
meanbreadth(k) = mb(obj);
volume(k) = vol(obj);
surface(k) = surf(obj);
euler(k) = eu(obj);
% mb(obj)
% all of your actions for filtering and plotting go here
end
meanbreadth
dlmwrite(sprintf('../Documents/mink128/%s_mean_breadth_%s.txt',orig_gen, type), meanbreadth)
volume
dlmwrite(sprintf('../Documents/mink128/%s_volume_%s.txt',orig_gen, type),volume)
surface
dlmwrite(sprintf('../Documents/mink128/%s_surface_%s.txt',orig_gen, type),surface)
euler
dlmwrite(sprintf('../Documents/mink128/%s_euler_%s.txt',orig_gen, type),euler)
end
end
function f = eu(obj)
% i = 1;
% while 1
% subimage = obj.read();
% % size(subimage)
% subimages(i,:,:) = subimage;
% if (obj.lastDirectory())
% break;
% end
% obj.nextDirectory()
% i = i+1;
% end
% size(subimages)
f = imEuler3d(obj);
end
function f = vol(obj)
% i = 1;
% while 1
% subimage = obj.read();
% % size(subimage)
% subimages(i,:,:) = subimage;
% if (obj.lastDirectory())
% break;
% end
% obj.nextDirectory()
% i = i+1;
% end
% size(subimages)
f = imVolume(obj);
end
function f = surf(obj)
% i = 1;
% while 1
% subimage = obj.read();
% % size(subimage)
% subimages(i,:,:) = subimage;
% if (obj.lastDirectory())
% break;
% end
% obj.nextDirectory()
% i = i+1;
% end
% size(subimages)
f = imSurface(obj);
end
function f = mb(obj)
% i = 1;
% while 1
% subimage = obj.read();
% % size(subimage)
% subimages(i,:,:) = subimage;
% if (obj.lastDirectory())
% break;
% end
% obj.nextDirectory()
% i = i+1;
% end
% size(subimages)
f = imMeanBreadth(obj, 13);
end