th(size(M))
disp(info.dime.dim(1))
disp(length(size(M)))
%%
%   info.dime.dim(2:5)
%   size(M)
disp(info.dime.dim(2:5))
disp(size(M))


%% Visualization in 2D
% View a 2D slice of the image using the Image Toolbox function _imagesc_
%
%   im = flipud(squeeze(M(65,:,:,1))');
%   imagesc(im)
%   colormap gray, axis tight

  im = flipud(squeeze(M(65,:,:,1))');
  imagesc(im)
  colormap gray, axis tight
  
%%
% Note the use of _flipud_ and _transpose (')_ to orient the
% image and _squeeze_ to reduce the image t