o 2D


%% Visualization in 3D
% View a 3D volume
%
%   S = smooth3(M(:,:,:,1));
%   isosurface(S,0)
%
% <<isosurface.png>>

%% Section plane through a 3D volume
% Create a section plane through a 3D volume, using information from the 
% NIFTI header to scale the axes
%
%   colormap gray
%   D = M(round(1:info.dime.dim(2)/2),:,:,1); % split along sagittal plane
%   Ds = smooth3(D);
%   hiso = patch(isosurface(Ds,0),...
%     'FaceColor',[.5,.5,.8],...
%     'EdgeColor','none');
%   isonormals(Ds,hiso)
%   hc