第三章 matlab作图3课件.ppt

6.14 立体图形的切片处理;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 光线对立体的影响;创建光线对象命令light，light命令的语法为： light(PropertyName, PropertyValue, ...) 这里特性propertyName 为： Color — 光源的颜色 Style — 光源的距离类型，即是无限远还是某距离的光源 Position — 光源的方位，某个角度 创建材质对象的命令为material，其语法为： material shiny — 抛光材质，即具有强烈的反光效应 material dull — 具有吸光效应的材质，反光度差 material metal — 金属材质，具有反光效应 利用参数ka，kd，ks自定义反光效果的命令为： material([ka kd ks]) material([ka kd ks n]) material([ka kd ks n sc]);【例5.4.6】显示材质的不同，对光线的反映不同 cc=[0,1,0]; p=30 subplot(2,2,1) sphere(p); title(material default ) subplot(2,2,2) sphere(p); shading interp set(gcf,Renderer,zbuffer) set(findobj(gca,type,surface),... FaceLighting,phong,... AmbientStrength,.3,Diffus