%% Extract Subset of Volume Data Set
% Load the mri data set, which is a collection of MRI slices of a human skull. Then:
% Squeeze the four-dimensional array into three dimensions and extract a subset using subvolume.
% Use a patch to display an isosurface of the skull outline.
% Add a second patch with interpolated face colors.
% Set the view of the object.
% Change the colormap to a grayscale colormap.
% Add lights to the right and left of the camera to illuminate the object.

load mri
D = squeeze(D);
[x,y,z,D] = subvolume(D,[60,80,nan,80,nan,nan]);
p1 = patch(isosurface(x,y,z,D, 5),...
     'FaceColor','red','EdgeColor','none');
isonormals(x,y,z,D,p1);
p2 = patch(isocaps(x,y,z,D, 5),...
     'FaceColor','interp','EdgeColor','none');
view(3); 
axis tight; 
daspect([1 1 0.4])
colormap(gray(100))
camlight right; 
camlight left; 
lighting gouraud

% Copyright 2012 The MathWorks, Inc.