xianshi.m

数据变换

fid=fopen('man.txt','r')
c=fread(fid,[192,inf],'int8')
d=reshape(c,192,96,498)
[x y z d] = subvolume(d, [nan nan nan nan nan nan]);
p = patch(isosurface(x,y,z,d, 5), 'FaceColor', 'red', 'EdgeColor', 'none');
p2 = patch(isocaps(x,y,z,d, 5), 'FaceColor', 'blue', 'EdgeColor', 'none');
isonormals(x,y,z,d,p);
figure;
colordef white
view(3)
daspect('auto')
colormap(hsv(128))
camva(9);
box on
camlight(40, 40);
camlight(-20,-10);
lighting gouraud