huitumri.m

来自「包括MATLAB7.0实用指南（1）书中所有的例子的源代码。」· M 代码 · 共 13 行

load mri
D=squeeze(D);
[x y z D]=reducevolume(D,[1,1,1]);
D=smooth3(D);
p1=patch(isosurface(x,y,z,D,5,'verbose'),...
   '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 .4])
colormap(gray(100))