tiqu.m

数据变换

figure;
colordef black
fid=fopen('man.txt','r');
e=fread(fid,[192,inf],'int8')
h=(e==12)*12
i=find(h~=0); 
n=length(i);
f=reshape(h,192,96,498)
[x y z f] = subvolume(f, [nan nan nan nan nan nan]);
p = patch(isosurface(x,y,z,f, 5), 'FaceColor', 'red', 'EdgeColor', 'none');
p2 = patch(isocaps(x,y,z,f,5), 'FaceColor', 'blue', 'EdgeColor', 'none');
isonormals(x,y,z,f,p);
view(3);
daspect('auto')
colormap(gray(100))
camva(9);
box on
camlight(40, 40);
camlight(-20,-10);
lighting gouraud