ap = patch(isocaps(D,0),...
%     'FaceColor','interp',...
%     'EdgeColor','none');
%   view(35,30) 
%   axis tight 
%   pixdim = info.dime.pixdim(2:4); % get voxel dimensions
%   daspect(1./[pixdim(1),pixdim(2),pixdim(3)]) % scale proportional to voxel shape
%   lightangle(90,0);
%   set(gcf,'Renderer','opengl'); lighting phong
%   set(hcap,'AmbientStrength',1)
%   set(hiso,'SpecularColorReflectance',0,'SpecularExponent',50)
%
% <<brainsection.png>>

%% MREView
% Use _MREView_ to visualize volumes of up 