[Perl+OpenGL]模仿贝塞尔曲线屏保效果

523066680

首先当然是要安装OpenGL模块。
旧版本的OpenGL自带Bezier求值器，但是并不返回坐标值，而是求值并帮你进行了顶点绘制操作。
为了能够更自由地操纵顶点数据，还是到wikipedia挖了一段函数：
zh.wikipedia.org/wiki/貝茲曲線 ，将 PointOnCubicBezier函数 转成Perl语言的格式

然后用50x50x3的结构数组对顶点进行存储，迭代替换，效果就做出来啦，应该还有很大优化空间


图片都超过1MB，其他的发到网络相册里了 paktc.lofter.com/post/45586_4f19e97

1. =info
   
2.     Code By: vicyang/5230
   
3.     E-mail : paktcmail@gmail.com
   
4.     DATE   : 2015-01-02
   
5. =cut
   
6. 
   
7. use v5.16;
   
8. use IO::Handle;
   
9. use OpenGL qw/ :all /;
   
10. use OpenGL::Config;
    
11. use Time::HiRes 'sleep';
    
12. 
    
13. STDOUT->autoflush(1);
    
14. 
    
15. our $WinID;
    
16. &Main();
    
17. 
    
18. sub delta {
    
19.     my ($pta, $ptb) = @_;
    
20.     my $dtx = ($ptb->[0] - $pta->[0]);
    
21.     my $dty = ($ptb->[1] - $pta->[1]);
    
22.     my $delta = sqrt ($dtx ** 2 + $dty ** 2);
    
23.     return ($delta, $dtx, $dty);
    
24. }
    
25. 
    
26. sub PointOnCubicBezier {
    
27.     my ($cp0, $cp1, $cp2, $cp3, $t) = @_;
    
28. 
    
29.     my ($cx, $bx, $ax, $cy, $by, $ay);
    
30.     my ($tSquared, $tCubed, $result);
    
31. 
    
32.     $cx = 3.0 * ($cp1->[0] - $cp0->[0]);
    
33.     $bx = 3.0 * ($cp2->[0] - $cp1->[0]) - $cx;
    
34.     $ax = $cp3->[0] - $cp0->[0] - $cx - $bx;
    
35. 
    
36.     $cy = 3.0 * ($cp1->[1] - $cp0->[1]);
    
37.     $by = 3.0 * ($cp2->[1] - $cp1->[1]) - $cy;
    
38.     $ay = $cp3->[1] - $cp0->[1] - $cy - $by;
    
39. 
    
40.     $tSquared = $t **2;
    
41.     $tCubed = $t **3;
    
42. 
    
43.     $result->[0] = (
    
44.             ($ax * $tCubed) +
    
45.             ($bx * $tSquared) +
    
46.             ($cx * $t) +
    
47.             $cp0->[0]
    
48.         );
    
49.     $result->[1] = (
    
50.             ($ay * $tCubed) +
    
51.             ($by * $tSquared) +
    
52.             ($cy * $t) +
    
53.             $cp0->[1]
    
54.         );
    
55. 
    
56.     return $result;
    
57. }
    
58. 
    
59. sub getRevPoint {
    
60.     my ($A, $O, $B) = (shift, shift, undef);
    
61. 
    
62.     $B->[0] = $O->[0] + ($O->[0] - $A->[0]);
    
63.     $B->[1] = $O->[1] + ($O->[1] - $A->[1]);
    
64. 
    
65.     return $B;
    
66. }
    
67. 
    
68. sub display {
    
69.     state $count = 0;
    
70.     state $loops = 0;
    
71. 
    
72.     my ($i, $j);
    
73.     my $coord;
    
74.    
    
75.     my ($dt, $dtx, $dty);
    
76. 
    
77.     state $n_cpts = 4;            #四个控制点
    
78.     state $cpts_path_parts = 50;  
    
79.     state $bezier_parts = 50;     #每条贝塞尔曲线分50部分
    
80.     state $n_bezier = 50;         #贝塞尔曲线的堆叠数量
    
81.     state $curve = [];            #每个cpt的曲线路径存储
    
82.     state $pti = 0;               #pti数组记录4个cpt的其自身的cpt'坐标
    
83.     state $pt = [
    
84.         [0.0, 0.0],
    
85.         [0.0, 0.0],
    
86.         [0.0, 0.0],
    
87.         [0.0, 0.0],
    
88.     ];
    
89. 
    
90.     state $ptarr = [];
    
91. 
    
92.     if (($loops == 0) and ($count == 0)) {
    
93.         for $i (0 .. $n_bezier-1) {
    
94.             for $j (0 .. $bezier_parts) {
    
95.                 $ptarr->[$i][$j] = [1000.0, -1000.0];
    
96.             }
    
97. 
    
98.             $curve->[$i] = [
    
99.                 [rand(200)-100.0, rand(200)-100.0],
    
100.                 [rand(200)-100.0, rand(200)-100.0],
     
101.                 [rand(200)-100.0, rand(200)-100.0],
     
102.                 [rand(200)-100.0, rand(200)-100.0]
     
103.             ]
     
104.         }
     
105.     }
     
106.    
     
107.     #处理四个控制点的演变
     
108. RECOUNT:
     
109.     if ( $pti <= $cpts_path_parts ) {
     
110.         for $i (0 .. $n_cpts-1) {
     
111.             $coord = &PointOnCubicBezier(
     
112.                     @{$curve->[$i]}[0..3],
     
113.                     $pti/$cpts_path_parts
     
114.             );
     
115.             $pt->[$i] = $coord;  #第 $i+1 个 cpt 的坐标
     
116.         }
     
117.         $pti++;
     
118.     } else {
     
119.         #控制点超过50次的时候，创建新的随机控制点，从头开始
     
120.         #起点是上一次的第50个点，同时也是新曲线的第0点
     
121.         #为了避免点重合，新的曲线点下标从1开始，而不是0开始
     
122.         for $i (0 .. $n_cpts-1) {
     
123.             $curve->[$i] = [
     
124.                 $pt->[$i],
     
125.                 getRevPoint($curve->[$i][2], $curve->[$i][3]),
     
126.                 #Control Point 3 is A, Point 4 is O, Find Point B
     
127.                 [rand(200)-100.0, rand(200)-100.0],
     
128.                 [rand(200)-100.0, rand(200)-100.0]
     
129.             ];
     
130.         }
     
131.         $pti = 1;
     
132.         goto RECOUNT;
     
133.     }
     
134. 
     
135. 
     
136.     glClear(GL_COLOR_BUFFER_BIT);
     
137.     glBlendFunc(GL_SRC_ALPHA, GL_DST_ALPHA);
     
138. 
     
139.     glPushMatrix();
     
140.         for $i (0 .. $bezier_parts ) {
     
141.             $coord = &PointOnCubicBezier(
     
142.                 @{$pt}[0, 1, 2, 3],
     
143.                 $i/$bezier_parts
     
144.             );
     
145.             $ptarr->[$count][$i] = $coord;
     
146.         }
     
147. 
     
148.         glBegin(GL_POINTS);
     
149.         
     
150. 
     
151.         for $i (0 .. $n_bezier-1) {
     
152.             glColor4f(0.4, 0.6, 0.8, 0.5);
     
153.             for $j (0.0 .. $bezier_parts) {
     
154.                 glVertex3f($ptarr->[$i][$j][0], $ptarr->[$i][$j][1], 0.0);
     
155.             }
     
156.         }
     
157.         glEnd();
     
158. 
     
159.     glPopMatrix();
     
160.     glutSwapBuffers();
     
161. 
     
162.     if ($count < ($n_bezier-1)) {
     
163.         $count++;
     
164.     } else {
     
165.         $count = 0;
     
166.         $loops++;
     
167.     }
     
168. 
     
169. }
     
170. 
     
171. sub init {
     
172.     glClearColor(0.0, 0.0, 0.0, 1.0);
     
173.     glPointSize(5.0);
     
174.     glLineWidth(10.0);
     
175.     glEnable(GL_BLEND);
     
176.     glEnable(GL_POINT_SMOOTH);
     
177.     glEnable(GL_LINE_SMOOTH);
     
178. }
     
179. 
     
180. sub idle {
     
181.     sleep 0.05;
     
182.     glutPostRedisplay();
     
183. }
     
184. 
     
185. sub Reshape {
     
186.     glViewport(0.0,0.0,500.0,500.0);
     
187.     glMatrixMode(GL_PROJECTION);
     
188.     glLoadIdentity();
     
189.     glOrtho(-100.0,100.0,-100.0,100.0,0.0,200.0);
     
190.     glMatrixMode(GL_MODELVIEW);
     
191.     glLoadIdentity();
     
192.     gluLookAt(0.0,0.0,100.0,0.0,0.0,0.0, 0.0,1.0,100.0);
     
193. }
     
194. 
     
195. sub hitkey {
     
196.     my $key = shift;
     
197.     if (lc(chr($key)) eq 'q') {
     
198.         glutDestroyWindow($WinID);
     
199.     } elsif ($key == 27) {
     
200.         glutDestroyWindow($WinID);
     
201.     }
     
202. }
     
203. 
     
204. sub Main {
     
205.     glutInit();
     
206.     glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE );#| GLUT_MULTISAMPLE);
     
207.     glutInitWindowSize(500, 500);
     
208.     glutInitWindowPosition(1,1);
     
209.     our $WinID = glutCreateWindow("Curve_Vic");
     
210.     &init();
     
211.     glutDisplayFunc(\&display);
     
212.     glutReshapeFunc(\&Reshape);
     
213.     glutKeyboardFunc(\&hitkey);
     
214.     glutIdleFunc(\&idle);
     
215.     glutMainLoop();
     
216. }

复制代码

huang6894

(๑ŐдŐb哇~点赞

死水站

好炫酷，牛！