第三讲 matlab作图3.pptVIP

* * 6.14 立体图形的切片处理 在科学研究中我们常常将物体进行切片处理，观察物体内部的结构。典型的有对人体的磁共震获得人体某部位的立体数据，然后方便地截取某截面进行深入研究。或将人体组织进行冷冻，然后进行病理切片进行病理研究。2001年全国大学生数学建模竞赛就是有关切片问题的例子。MATLAB提供了各种命令可以对立体进行切片处理。 【例6.14.1】 对射流数据进行处理，现在我们对它进行切片处理，并从不同的角度来看切片的图像。结果见图5.4.10 [x,y,z,v] = flow; ax=[-38.5,16;0,45] % 建立视角数据 for m=1:2 % 建立切片大小 xmin = min(x(:)); ymin = min(y(:)); zmin = min(z(:)); xmax = max(x(:)); ymax = max(y(:)); zmax = max(z(:)); subplot(1,2,m) % 同时作两幅图 hslice = surf(linspace(xmin,xmax,100),... linspace(ymin,ymax,100),... zeros(100)); % 建立切片 rotate(hslice,[-2,0,0],-45) %旋转切片45度，形成新的切片。 xd = get(hslice,XData); yd = get(hslice,YData); zd = get(hslice,ZData); delete(hslice) % 将原切片删除 % 对切片上色、插值加密等附加修饰 h = slice(x,y,z,v,xd,yd,zd); set(h,FaceColor,interp,... EdgeColor,none,... DiffuseStrength,.8) daspect([1,1,1]) axis tight box off view(ax(m,:)) camzoom(1.4) camproj perspective lightangle(-45,45) colormap (jet(24)) set(gcf,Renderer,zbuffer) end 见下图 【例5.4.9】本例为2001全国大学生数学建模竞赛A题，题例见第十章。再现血管的原状并切片进行比较来验证恢复的正确性。图5.4.11为血管的再现图，图5.4.12为切片比较图。 flops(0);tic load (f:\cmcmtu\Result) %读取数据，中轴线上100点的坐标，半径 imgnumber=100 R=sum(Result(:,1))/imgnumber; %Result(:,2:3)=Result(:,2:3) z=1:imgnumber; Poly_zx=polyfit(z,Result(z,2),6); Poly_zy=polyfit(z,Result(z,3),6); x=polyval(Poly_zx,z); y=polyval(Poly_zy,z);% 用多项式拟合中轴方程并在上面取100个点 % 作三维图 for k=1:imgnumber [X,Y,Z]=ellipsoid(x(k),y(k),z(k),R,R,R,10); h=surf(X,Y,Z); set(h,FaceColor,interp,... EdgeColor,none,... DiffuseStrength,.8) hold on end hold on alpha(0.2) hslice = surf(linspace(0,512,100),... linspace(0,512,100),... zeros(100)+30) % 建立切片 colormap (jet(24)) delete(h) 三维血管复原图 对血管进行切片 切片比较图 6.15 光线对立体的影响 我们都知道光线对三维图形有影响，如一个球正对光线的部分发亮，而背光处则发暗。对光的反映还根据物体的材料不同而有不同，例如强光照在一个网球上没有非常强烈的反光，而如果强光照在一个抛光的金属求上时则有非常强烈的反光。所有这些现象在MATLAB中都可以模拟，既有各种光线强度、位置的