Example Script for COmputing on Hazelnut Data

example.m

filename = 'macadamia.mat'
if ~exist( filename, 'file' );
urlwrite( ... Write a url to a file
    'https://dl.dropboxusercontent.com/u/22455492/natural-3d/macadamia.mat', ... The requested url
    filename )
end % if:~exist

%%

load( filename ); % Variable is input as structure

%% Calculate Volume Fraction

volumefraction = mean( structure(:) );

return

%%

shift = [eye(3); -1*eye(3)];

edge = zeros( size( structure ) );
for ii = 1 : size( shift,1);
    temp = structure - circshift( structure, shift(ii,:));
    temp(:) = abs(temp);
    edge = edge | temp;
end



%% Visualize the dataset

if ~exist( 'slices', 'dir' )
    mkdir( 'slices' );
end
clc;
for ii = 1 : 80 : 500
    pcolor( double(structure(:,:,ii)) );
    figure(gcf)
    axis off
    colormap gray
    title(sprintf('Black::Pore and White::Solid @ Slice::%i',ii),...
        'Fontsize',16 );
    saveas( gcf,fullfile('.','slices',sprintf('Slice-%i.png',ii) ) )
end