huitumri2.m

包括MATLAB7.0实用指南（1）书中所有的例子的源代码。

load mri
D=squeeze(D);
[x y z D]=subvolume(D,[60,80,nan,80,nan,nan]);
%D=smooth3(D);
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(60,30);
axis tight;
daspect([1 1 .4])
colormap(gray(100))