绘制三维人脸的代码分享。

Colbie

这几天整理之前写过的matlab程序，偶然发现绘制三维人脸的代码。
通过参数方程分别绘制了头部、眼睛、鼻子、嘴巴四个部分，概略绘制效果还好。
6 `/ c7 \/ K5 g0 ?) B在此向坛友们分享

( b& j8 U- Y7 N0 _2 A代码粘贴如下：
% R# H! h) a0 Z: H- P6 u! L0 k

• <P>
• t=-pi:0.1:pi;
• [t,a]=meshgrid(t);
• y=5.*sin(a).*cos(t);
• z=10.*sin(a).*sin(t);
• x=cos(a);
• mesh(x,y,z)
• xlabel('X');
• ylabel('Y');
• zlabel('Z');
• hold on
• y1=sin(a).*cos(t)+2;
• z1=sin(a).*sin(t)+3;
• x1=cos(a)+0.5;
• mesh(x1,y1,z1)
• view(90,0)
• y2=sin(a).*cos(t)-2;
• z2=sin(a).*sin(t)+3;
• x2=cos(a)+0.5;
• hold on
• mesh(x2,y2,z2)
• x3=1/6*a.*cos(t)+0.85;
• y3=1/3*a.*sin(t);
• z33=-2*a+4.5;
• i=find(z33>4.5);
• z33(i)=NaN;
• z3=z33;
• hold on
• mesh(x3,y3,z3)</P>
• <P>x4=2/3*a.*cos(t)+0.6;
• y4=2*a.*sin(t);
• z44=a.*a-5;
• j=find(z44>-4);
• z44(j)=NaN;
• z4=z44;
• hold on
• mesh(x4,y4,z4)</P>
• <P>x5=2/3*a.*cos(t)+0.6;
• y5=2*a.*sin(t);
• z5=-a.*a-2.999;
• k=find(z5<-3.999);
• z5(k)=NaN;
• hold on
• mesh(x5,y5,z5)
• </P>
  0 g  h0 c6 \& X5 B1 L8 w1 j

4 E1 w! o( n( F
; w  L" |1 L1 ~. @7 V* S
, U9 C) ]! F) r  e0 V6 r           
2 _/ `  ~: H# e9 F8 _% U1 C
9 |) i# }- c3 r运行结果：


. }+ z! w, d4 j; g. T
4 p7 O6 d0 c4 D. x8 S' g6 Y4 j3 ^, r" ~9 Z



+ h2 d8 p/ z5 ~$ c, V1 d7 u8 P

2 n* Y7 C2 E4 l4 i, v7 _7 L9 a  O

somethingabc

去试试，谢谢楼主分享