Perl 绘制3D模型文件 - 斯坦福兔子 (ply格式)

523066680

模型文件下载地址
http://www.cc.gatech.edu/projects/large_models/bunny.html

以及针对这只兔子的渲染代码，没有做成通用的，在加载不同模型的时候需要自己调整大小和位置以适应窗口

使用说明：需要 OpenGL 模块的支持，成功渲染后，按 w a s d, j k 可以旋转观察角度。

1. =Copy is not Right
   
2.     Author: 523066680@163.com
   
3.       Date: 2016-12
   
4. =cut
   
5. 
   
6. use v5.16;
   
7. use IO::Handle;
   
8. use OpenGL qw/ :all /;
   
9. use OpenGL::Config;
   
10. use Time::HiRes 'sleep';
    
11. use feature 'state';
    
12. STDOUT->autoflush(1);
    
13. 
    
14. our $WinID;
    
15. our $width  = 500;
    
16. our $height = 500;
    
17. 
    
18. our $rx = 0.0;
    
19. our $ry = 0.0;
    
20. our $rz = 0.0;
    
21. 
    
22. my @vtx;
    
23. my @faces;
    
24. 
    
25. open READ, "<:raw", "bunny.ply";
    
26. my $models = 0;
    
27. 
    
28. for  my $e ( <READ> )
    
29. {
    
30.     $e =~ s/\r?\n$//;
    
31.     if ($e =~ /v ([^\s]+) ([^\s]+) ([^\s]+)/i )
    
32.     {
    
33.         #考虑 科学表示法的情况 5.5e-05
    
34.         push @vtx, [ sprintf("%f", $1) * 1500.0,
    
35.                      sprintf("%f", $2) * 1500.0,
    
36.                      sprintf("%f", $3) * 1500.0 ] ;
    
37.     }
    
38.     elsif ($e =~ /f (\d+) (\d+) (\d+)/i )
    
39.     {
    
40.         #某些模型文件的三角形索引从 1 开始
    
41.         push @faces,
    
42.             [
    
43.                 $1 - 1,
    
44.                 $2 - 1,
    
45.                 $3 - 1
    
46.             ];
    
47.     }
    
48.     else
    
49.     {
    
50.         printf("Other: $e\n");
    
51.     }
    
52. 
    
53.     if ($e =~ /feature/i )
    
54.     {
    
55.         $models++;
    
56.         if ($models > 1)
    
57.         {
    
58.             #last;
    
59.         }
    
60.     }
    
61. }
    
62. close READ;
    
63. 
    
64. printf("%d %d\n", $#vtx, $#faces);
    
65. 
    
66. printf("%f %f %f\n", $vtx[0]->[0], $vtx[0]->[1], $vtx[0]->[2] );
    
67. my @nums;
    
68. grep { push @nums, @{$_} } @faces;
    
69. @nums = (sort { $a <=> $b } @nums);
    
70. printf("\n %d\n", $nums[$#nums]);
    
71. 
    
72. &Main();
    
73. 
    
74. sub normalize
    
75. {
    
76.     my $v = shift;
    
77.     my $d = sqrt($v->[0]*$v->[0] + $v->[1]*$v->[1] + $v->[2]*$v->[2]);
    
78.     if ($d == 0.0)
    
79.     {
    
80.         printf("length zero!\n");
    
81.         return;
    
82.     }
    
83. 
    
84.     $v->[0] /= -$d;
    
85.     $v->[1] /= -$d;
    
86.     $v->[2] /= -$d;
    
87. }
    
88. 
    
89. sub normcrossprod
    
90. {
    
91.     my ( $v1, $v2, $out ) = @_;
    
92. 
    
93.     $out->[0] = $v1->[1] * $v2->[2] - $v1->[2] * $v2->[1];
    
94.     $out->[1] = $v1->[2] * $v2->[0] - $v1->[0] * $v2->[2];
    
95.     $out->[2] = $v1->[0] * $v2->[1] - $v1->[1] * $v2->[0];
    
96. 
    
97.     normalize( $out );
    
98. }
    
99. 
    
100. sub initlist
     
101. {
     
102.     my $idx;
     
103.     my $base = glGenLists(1);
     
104.     my ($a, $b, $c);
     
105. 
     
106.     my @tpa;
     
107.     my @tpb;
     
108.     my @norm = ();
     
109. 
     
110.     glNewList($base, GL_COMPILE);
     
111. 
     
112.     glBegin(GL_TRIANGLES);
     
113. 
     
114.     for my $i ( 0 ..  $#faces )
     
115.     {
     
116.         ($a, $b, $c) = @{ $faces[$i] } ;
     
117. 
     
118.         for ( 0 .. 2 )
     
119.         {
     
120.             $tpa[$_] = $vtx[$b]->[$_] - $vtx[$a]->[$_] ;
     
121.             $tpb[$_] = $vtx[$c]->[$_] - $vtx[$a]->[$_] ;
     
122.         }
     
123. 
     
124.         glColor3f( 0.9, 0.6, 0.2 );
     
125.         normcrossprod( \@tpa, \@tpb, \@norm );
     
126. 
     
127.         glNormal3f( @norm );
     
128.         glVertex3f( @{$vtx[$a]} );
     
129.         glVertex3f( @{$vtx[$b]} );
     
130.         glVertex3f( @{$vtx[$c]} );
     
131.     }
     
132.     glEnd();
     
133. 
     
134.     glEndList();
     
135. }
     
136. 
     
137. sub display
     
138. {
     
139.     my $idx;
     
140.     glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
     
141. 
     
142.     glPushMatrix();
     
143.     glRotatef($rx, 1.0, 0.0, 0.0);
     
144.     glRotatef($ry, 0.0, 1.0, 0.0);
     
145.     glRotatef($rz, 0.0, 0.0, 1.0);
     
146.     glTranslatef(0.0, -150.0, 0.0);
     
147.     glColor4f(0.8, 0.8, 0.5, 0.5);
     
148. 
     
149.     glCallList(1);
     
150. 
     
151.     glPopMatrix();
     
152.     glutSwapBuffers();
     
153. }
     
154. 
     
155. sub init
     
156. {
     
157.     glClearColor(0.0, 0.0, 0.0, 0.5);
     
158.     glPointSize(1.0);
     
159.     glLineWidth(1.0);
     
160.     glEnable(GL_DEPTH_TEST);
     
161.     glEnable(GL_POINT_SMOOTH);
     
162.     glEnable(GL_LINE_SMOOTH);
     
163.    
     
164.     my @mat_ambient = ( 0.5, 0.5, 0.5, 1.0 );
     
165.     my @mat_specular = ( 1.0, 1.0, 1.0, 1.0 );
     
166.     my @mat_shininess = ( 10.0 );
     
167.     my @mat_diffuse = (0.5, 0.5, 0.5, 1.0);
     
168. 
     
169.     my @light_specular = ( 0.5, 0.5, 0.5, 1.0 );
     
170.     my @light_position = ( 200.0, 400.0, 0.0, 0.0 );
     
171. 
     
172.     my $ambient  = OpenGL::Array->new( 4, GL_FLOAT);
     
173.     my $specular = OpenGL::Array->new( 4, GL_FLOAT);
     
174.     my $diffuse  = OpenGL::Array->new( 4, GL_FLOAT);
     
175.     my $shininess = OpenGL::Array->new( 1, GL_FLOAT);
     
176. 
     
177.     my $light_position = OpenGL::Array->new( 4, GL_FLOAT);
     
178.     my $light_specular = OpenGL::Array->new( 4, GL_FLOAT);
     
179. 
     
180.     $ambient->assign(0,  ( 0.5, 0.5, 0.5, 1.0 ) );
     
181.     $specular->assign(0, ( 1.0, 1.0, 1.0, 1.0 ) );
     
182.     $diffuse->assign(0,  ( 0.5, 0.5, 0.5, 1.0 ) );
     
183.     $shininess->assign(0,  10.0 );
     
184. 
     
185.     $light_position->assign(0, ( 200.0, 200.0, 0.0, 0.0 ) );
     
186. 
     
187.     glMaterialfv_c(GL_FRONT_AND_BACK, GL_AMBIENT, $ambient->ptr );
     
188.     glMaterialfv_c(GL_FRONT_AND_BACK, GL_SPECULAR, $specular->ptr );
     
189.     glMaterialfv_c(GL_FRONT_AND_BACK, GL_DIFFUSE, $diffuse->ptr );
     
190.     glMaterialfv_c(GL_FRONT_AND_BACK, GL_SHININESS, $shininess->ptr );
     
191. 
     
192.     glLightfv_c(GL_LIGHT0, GL_POSITION, $light_position->ptr);
     
193. 
     
194.     glEnable(GL_LIGHTING);
     
195.     glEnable(GL_LIGHT0);
     
196. 
     
197.     glColorMaterial(GL_FRONT, GL_DIFFUSE);
     
198.     glEnable(GL_COLOR_MATERIAL);
     
199. 
     
200.     initlist();
     
201. }
     
202. 
     
203. sub idle
     
204. {
     
205.     sleep 0.001;
     
206.     glutPostRedisplay();
     
207. }
     
208. 
     
209. sub Reshape
     
210. {
     
211.     my ($w, $h) = (shift, shift);
     
212.     my $half = 200.0;
     
213.     my $fa = 250.0;
     
214. 
     
215.     glViewport(0, 0, $w, $h);
     
216.     glMatrixMode(GL_PROJECTION);
     
217.     glLoadIdentity();
     
218.     #glOrtho(-$half, $half, -$half, $half, 0.0, $fa*2.0);
     
219.     gluPerspective( 90.0, 1.0, 1.0, $fa*2.0 );
     
220.     glMatrixMode(GL_MODELVIEW);
     
221.     glLoadIdentity();
     
222.     gluLookAt(0.0,0.0,$fa,0.0,0.0,0.0, 0.0,1.0, $fa);
     
223. }
     
224. 
     
225. sub hitkey
     
226. {
     
227.     our $WinID;
     
228.     my $k = lc(chr(shift));
     
229. 
     
230.     given ($k)
     
231.     {
     
232.         when ('q') { glutDestroyWindow( $WinID ); }
     
233.         when ('a') { $ry-=3.0; }
     
234.         when ('d') { $ry+=3.0; }
     
235.         when ('w') { $rx+=3.0; }
     
236.         when ('s') { $rx-=3.0; }
     
237.         when ('j') { $rz+=3.0; }
     
238.         when ('k') { $rz-=3.0; }
     
239.     }
     
240. 
     
241. }
     
242. 
     
243. sub Main
     
244. {
     
245.     glutInit();
     
246.     glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE |GLUT_DEPTH |GLUT_MULTISAMPLE );
     
247.     glutInitWindowSize(500, 500);
     
248.     glutInitWindowPosition(1,1);
     
249.     our $WinID = glutCreateWindow("bunny");
     
250.     &init();
     
251.     glutDisplayFunc(\&display);
     
252.     glutReshapeFunc(\&Reshape);
     
253.     glutKeyboardFunc(\&hitkey);
     
254.     glutIdleFunc(\&idle);
     
255.     glutMainLoop();
     
256. }
     

复制代码

523066680

强行转 Perl 代码发帖。

转换过程中感到的明显缺点是，大量的敲打 $ 符号。
爽快的地方：

• given when 语句比 switch case 整洁多了。
• 从纯文本文件中读取数据、处理数据相对方便
  

codegay

尼玛，好丑啊。哈哈。

523066680

回复 3# codegay


    另一种扫描没这么精细的三角形网格反而好看一些。
http://graphics.stanford.edu/pub/3Dscanrep/bunny.tar.gz

bun_zipper_res3.ply



代码稍作修改