OpenGL基础图形编程（八）变换

八、OpenGL变换

　　OpenGL变换是本篇的重点内容，它包含计算机图形学中最主要的三维变换，即几何变换、投影变换、裁剪变换、视口变换，以及针对OpenGL的特殊变换概念理解和使用方法，如相机模拟、矩阵堆栈等。

学好了这章，才開始真正走进三维世界。
8.1、从三维空间到二维平面

　　8.1.1 相机模拟
　　在真实世界里，全部的物体都是三维的。

可是。这些三维物体在计算机世界中却必须以二维平面物体的形式表现出来。

那么，这些物体是如何从三维变换到二维的呢？以下我们採用相机（Camera）模拟的方式来讲述这个概念，如图8-1所看到的。

图8-1 相机模拟

　　实际上，从三维空间到二维平面，就如同用相机拍照一样，通常都要经历下面几个步骤 （括号内表示的是对应的图形学概念）：
　　第一步，将相机置于三角架上，让它对准三维景物（视点变换，Viewing Transformation）。
　　第二步。将三维物体放在适当的位置（模型变换。Modeling Transformation）。
　　第三步，选择相机镜头并调焦，使三维物体投影在二维胶片上（投影变换，Projection Transformation）。

　　第四步，决定二维像片的大小（视口变换，Viewport Transformation）。

　　这样。一个三维空间里的物体就能够用对应的二维平面物体表示了。也就能在二维的电脑屏幕上正确显示了。

　　8.1.2 三维图形显示流程
　　运用相机模拟的方式比較通俗地解说了三维图形显示的基本过程，但在详细应用OpenGL函数库编程时。还必须了解三维图形世界中的几个特殊坐标系的概念。以及用这些概念表达的三维图形显示流程。

　　计算机本身仅仅能处理数字。图形在计算机内也是以数字的形式进行加工和处理的。

大家都知道。坐标建立了图形和数字之间的联系。为了使被显示的物体数字化，要在被显示的物体所在的空间中定义一个坐标系。这个坐标系的长度单位和坐标轴的方向要适合对被显示物体的描写叙述，这个坐标系称为世界坐标系。

　　计算机对数字化的显示物体作了加工处理后，要在图形显示器上显示，这就要在图形显示器屏幕上定义一个二维直角坐标系，这个坐标系称为屏幕坐标系。

这个坐标系坐标轴的方向通常取成平行于屏幕的边缘。坐标原点取在左下角，长度单位常取成一个象素的长度，大小能够是整型数。
　　为了使显示的物体能以合适的位置、大小和方向显示出来，必需要通过投影。投影的方法有两种，即正射投影和透视投影。

　　有时为了突出图形的一部分，仅仅把图形的某一部分显示出来，这时能够定义一个三维视景体（Viewing Volume）。正射投影时通常是一个长方体的视景体。透视投影时通常是一个棱台似的视景体。仅仅有视景体内的物体能被投影在显示平面上。其它部分则不能。在屏幕窗体内能够定义一个矩形，称为视口（Viewport），视景体投影后的图形就在视口内显示。

　　为了适应物理设备坐标和视口所在坐标的区别，还要作一适应物理坐标的变换。这个坐标系称为物理设备坐标系。依据上面所述。三维图形的显示流程应如图8-2所看到的。

图8-2 三维图形的显示流程

　　8.1.3 基本变换简单分析
　　以下举一个简单的变换样例。cube.c：

　　例8-4 简单变换例程（cube.c）　　

#include "glos.h"
　　#include <GL/gl.h>
　　#include <GL/glu.h>
　　#include <GL/glaux.h>

　　void myinit(void);
　　void CALLBACK myReshape(GLsizei w, GLsizei h);
　　void CALLBACK display(void);

　　void CALLBACK display (void)
　　{
　　　　glClear(GL_COLOR_BUFFER_BIT);
　　　　glColor3f (1.0, 1.0, 1.0);
　　　　glLoadIdentity (); /* clear the matrix */
　　　　glTranslatef (0.0, 0.0, -5.0); /* viewing transformation */
　　　　glScalef (1.0, 2.0, 1.0); /* modeling transformation */
　　　　auxWireCube(1.0); /* draw the cube */
　　　　glFlush();
　　}

　　void myinit (void)
　　{
　　　　glShadeModel (GL_FLAT);
　　}

　　void CALLBACK myReshape(GLsizei w, GLsizei h)
　　{
　　　　glMatrixMode (GL_PROJECTION); /* prepare for and then */
　　　　glLoadIdentity (); /* define the projection */
　　　　glFrustum (-1.0, 1.0, -1.0, 1.0, 1.5, 20.0); /* transformation */
　　　　glMatrixMode (GL_MODELVIEW); /* back to modelview matrix */
　　　　glViewport (0, 0, w, h); /* define the viewport */
　　}

　　void main(void)
　　{
　　　　auxInitDisplayMode (AUX_SINGLE | AUX_RGBA);
　　　　auxInitPosition (0, 0, 500, 500);
　　　　auxInitWindow ("Perspective 3-D Cube");
　　　　myinit ();
　　　　auxReshapeFunc (myReshape);
　　　　auxMainLoop(display);
　　}

　　以上程序执行结果就是绘制一个三维的正面透视立方体。当中已经用到了相机模拟中提到的四种基本变换，即视点变换、模型变换、投影变换和视口变换。

图8-3 三维的正面透视立方体

　　以下简单分析一下整个程序过程：
　　1）视点变换。视点变换是在视点坐标系中进行的。视点坐标系于一般的物体所在的世界坐标系不同。它遵循左手法则，即左手大拇指指向Z正轴，与之垂直的四个手指指向X正轴，四指弯曲90度的方向是Y正轴。

而世界坐标系遵循右手法则的。如图8-4所看到的。当矩阵初始化glLoadIdentity()后。调用glTranslatef()作视点变换。函数參数(x, y, z)表示视点或相机在视点坐标系中移动的位置。这里z=-5.0。意思是将相机沿Z负轴移动5个单位。

　　通常相机位置缺省值同场景中的物体一样，都在原点处，并且相机初始方向都指向Z负轴。
　　这里相机移走后。仍然对准立方体。假设相机须要指向还有一方向。则调用glRotatef()能够改变。

图8-4 视点坐标系与世界坐标系

　　2）模型变换。模型变换是在世界坐标系中进行的。在这个坐标系中，能够对物体实施平移 glTranslatef()、旋转glRotatef()和放大缩小glScalef()。样例里仅仅对物体进行比例变换。glScalef(sx, sy, sz)的三个參数各自是X、Y、Z轴向的比例变换因子。缺省时都为1.0。即物体没变化。程序中物体Y轴比例为2.0，其余都为1.0，就是说将立方体变成长方体。
　　3）投影变换。

投影变换类似于选择相机的镜头。

本例中调用了一个透视投影函数 glFrustum()，在调用它之前先要用glMatrixMode()说明当前矩阵方式是投影GL_PROJECTION。这个投影函数一共同拥有六个參数。由它们能够定义一个棱台似的视景体。

即视景体内的部分可见，视景体外的部分不可见，这也就包括了三维裁剪变换。
　　4）视口变换。视口变换就是将视景体内投影的物体显示在二维的视口平面上。通常。都调用函数glViewport()来定义一个视口，这个过程类似于将照片放大或缩小。

　　总而言之，一旦全部必要的变换矩阵被指定后，场景中物体的每个顶点都要依照被指定的变换矩阵序列逐一进行变换。

注意：OpenGL 中的物体坐标一律採用齐次坐标，即(x, y, z, w)。故全部变换矩阵都採用4X4矩阵。一般说来，每一个顶点先要经过视点变换和模型变换，然后进行指定的投影。假设它位于视景体外。则被裁剪掉。最后。余下的已经变换过的顶点x、y、z坐标值都用比例因子w除。即x/w、y/w、z/w。再映射到视口区域内，这样才干显示在屏幕上。

8.2、几何变换
　　实际上。上述所说的视点变换和模型变换本质上都是一回事，即图形学中的几何变换。
　　仅仅是视点变换一般仅仅有平移和旋转，没有比例变换。当视点进行平移或旋转时。视点坐标系中的物体就相当于在世界坐标系中作反方向的平移或旋转。

因此，从某种意义上讲，二者能够统一。仅仅是各自出发点不一样而已。

读者能够依据详细情况。选择当中一个角度去考虑，这样便于理解。

　　8.2.1 两个矩阵函数解释
　　这里先解释两个基本OpenGL矩阵操作函数，便于以后章节的讲述。函数解释例如以下：

void glLoadMatrix{fd}(const TYPE *m)

　　设置当前矩阵中的元素值。函数參数*m是一个指向16个元素(m0, m1, ..., m15)的指针。这16个元素就是当前矩阵M中的元素。其排列方式例如以下：

M =	| m0 m4 m8　m12 |
        | m1 m5 m9　m13 |
        | m2 m6 m10 m14 |
        | m3 m7 m11 M15 |

void glMultMatrix{fd}(const TYPE *m)

　　用当前矩阵去乘*m所指定的矩阵。并将结果存放于*m中。

当前矩阵能够是用glLoadMatrix() 指定的矩阵，也能够是其他矩阵变换函数的综合结果。
　　当几何变换时，调用OpenGL的三个变换函数glTranslate*()、glRotate*()和glScale*()，实质上相当于产生了一个近似的平移、旋转和比例矩阵，然后调用glMultMatrix()与当前矩阵相乘。可是直接调用这三个函数程序执行得快一些，因OpenGL自己主动能计算矩阵。

　　8.2.2 平移
　　平移变换函数例如以下：　

void glTranslate{fd}(TYPE x,TYPE y,TYPE z)

　　三个函数參数就是目标分别沿三个轴向平移的偏移量。这个函数表示用这三个偏移量生成的矩阵乘以当前矩阵。

当參数是(0.0,0.0,0.0)时。表示对函数glTranslate*()的操作是单位矩阵，也就是对物体没有影响。

平移示意如图8-5所看到的。

图8-5 平移示意图

　　8.2.3 旋转
　　旋转变换函数例如以下：　　

void glRotate{fd}(TYPE angle,TYPE x,TYPE y,TYPE z)

　　函数中第一个參数是表示目标沿从点(x, y, z)到原点的方向逆时针旋转的角度，后三个參数是旋转的方向点坐标。

这个函数表示用这四个參数生成的矩阵乘以当前矩阵。当角度參数是0.0时，表示对物体没有影响。旋转示意如图8-6所看到的。

图8-6 旋转示意图

　　8.2.3 缩放和反射
　　缩放和反射变换函数例如以下：　

void glScale{fd}(TYPE x,TYPE y,TYPE z)

　　三个函数參数值就是目标分别沿三个轴向缩放的比例因子。这个函数表示用这三个比例因子生成的矩阵乘以当前矩阵。这个函数能完毕沿对应的轴对目标进行拉伸、压缩和反射三项功能。当參数是(1.0, 1.0, 1.0)时。表示对函数glScale*()操作是单位矩阵。也就是对物体没有影响。当当中某个參数为负值时，表示将对目标进行对应轴的反射变换，且这个參数不为1.0。则还要进行对应轴的缩放变换。最好不要令三个參数值都为零。这将导致目标沿三轴都缩为零。缩放和反射示意如图8-7所看到的。

图8-7 缩放和反射示意图

　　8.2.5 几何变换举例
　　以上介绍了三个基本几何变换函数，以下举一个简单的样例进一步说明它们的使用方法。

程序例如以下：

　　例 8-5 几何变换例程（geomtrsf.c）　　

#include "glos.h"
　　#include <GL/gl.h>
　　#include <GL/glu.h>
　　#include <GL/glaux.h>

　　void myinit(void);
　　void draw_triangle(void);
　　void CALLBACK display(void);
　　void CALLBACK myReshape(GLsizei w, GLsizei h);

　　void draw_triangle(void)
　　{
　　　　glBegin(GL_LINE_LOOP);
　　　　　　glVertex2f(0.0, 25.0);
　　　　　　glVertex2f(25.0, -25.0);
　　　　　　glVertex2f(-25.0, -25.0);
　　　　glEnd();
　　}

　　void CALLBACK display(void)
　　{
　　　　glClearColor (0.0, 0.0, 0.0, 1.0);
　　　　glClear (GL_COLOR_BUFFER_BIT); 

　　　　/* draw an original triangle */
　　　　glLoadIdentity ();
　　　　glColor3f (1.0, 1.0, 1.0);　　/* white */
　　　　draw_triangle ();

　　　　/* translating a triangle along X_axis */
　　　　glLoadIdentity ();
　　　　glTranslatef (-20.0, 0.0, 0.0);
　　　　glColor3f(1.0,0.0,0.0);　　　/* red */
　　　　draw_triangle ();

　　　　/* scaling a triangle along X_axis by 1.5 and along Y_axis by 0.5 */
　　　　glLoadIdentity();
　　　　glScalef (1.5, 0.5, 1.0);
　　　　glColor3f(0.0,1.0,0.0);　　　/* green */
　　　　draw_triangle ();

　　　　/* rotating a triangle in a counterclockwise direction about Z_axis */
　　　　glLoadIdentity ();
　　　　glRotatef (90.0, 0.0, 0.0, 1.0);
　　　　glColor3f(0.0,0.0,1.0);　　　/* blue */
　　　　draw_triangle ();

　　　　/* scaling a triangle along Y_axis and reflecting it about Y_axis */
　　　　glLoadIdentity();
　　　　glScalef (1.0, -0.5, 1.0);
　　　　glColor3f(1.0,1.0,0.0); 　　/* yellow */
　　　　draw_triangle ();

　　　　glFlush();
　　}

　　void myinit (void)
　　{
　　　　glShadeModel (GL_FLAT);
　　}

　　void CALLBACK myReshape(GLsizei w, GLsizei h)
　　{
　　　　glViewport(0, 0, w, h);
　　　　glMatrixMode(GL_PROJECTION);
　　　　glLoadIdentity();
　　　　if (w <= h)
　　　　　　glOrtho(-50.0, 50.0, -50.0*(GLfloat)h/(GLfloat)w, 50.0*(GLfloat)h/(GLfloat)w,-1.0,1.0); 
　　　　else
　　　　　　glOrtho(-50.0*(GLfloat)w/(GLfloat)h, 50.0*(GLfloat)w/(GLfloat)h, -50.0, 50.0,-1.0,1.0); 
　　　　glMatrixMode(GL_MODELVIEW);
　　}

　　void main(void)
　　{
　　　　auxInitDisplayMode (AUX_SINGLE | AUX_RGBA);
　　　　auxInitPosition (0, 0, 500, 500);
　　　　auxInitWindow ("Geometric Transformations");
　　　　myinit ();
　　　　auxReshapeFunc (myReshape);
　　　　auxMainLoop(display);
　　}

　　以上程序执行结果：第一个白色三角形是原始三角形，第二个红色三角形是白三角沿X 负轴平移后的三角形。第三个绿色三角形是白三角分别沿X轴和Y轴比例变换后的三角形，第四个蓝色三角形是白三角绕Z正轴逆时针转90度后的三角形。第五个黄色三角形是白三角沿Y轴方向缩小一倍且相对于X轴作反射后形成的三角形。

图8-8 三角形的几何变换

8.3、投影变换

　　投影变换是一种非常关键的图形变换。OpenGL中仅仅提供了两种投影方式，一种是正射投影，还有一种是透视投影。无论是调用哪种投影函数，为了避免不必要的变换，其前面必须加上下面两句：

glMAtrixMode(GL_PROJECTION);
glLoadIdentity();

　　其实。投影变换的目的就是定义一个视景体，使得视景体外多余的部分裁剪掉，终于图像仅仅是视景体内的有关部分。

本指南将具体讲述投影变换的概念以及使用方法。

　　8.3.1 正射投影（Orthographic Projection）
　　正射投影，又叫平行投影。这样的投影的视景体是一个矩形的平行管道，也就是一个长方体，如图8-9所看到的。

正射投影的最大一个特点是不管物体距离相机多远，投影后的物体大小尺寸不变。这样的投影通经常使用在建筑蓝图绘制和计算机辅助设计等方面，这些行业要求投影后的物体尺寸及相互间的角度不变。以便施工或制造时物体比例大小正确。

图8-9 正射投影视景体

　　OpenGL正射投影函数共同拥有两个。这在前面几个样例中已用过。一个函数是：　

void glOrtho(GLdouble left,GLdouble right,GLdouble bottom,GLdouble top,
　GLdouble near,GLdouble far)

　　它创建一个平行视景体。实际上这个函数的操作是创建一个正射投影矩阵，而且用这个矩阵乘以当前矩阵。当中近裁剪平面是一个矩形。矩形左下角点三维空间坐标是（left，bottom，-near）。右上角点是（right，top。-near）；远裁剪平面也是一个矩形，左下角点空间坐标是（left，bottom，-far）。右上角点是（right。top，-far）。全部的near和far值同一时候为正或同一时候为负。假设没有其它变换。正射投影的方向平行于Z轴，且视点朝向Z负轴。
　　这意味着物体在视点前面时far和near都为负值，物体在视点后面时far和near都为正值。

还有一个函数是：　　

void gluOrtho2D(GLdouble left,GLdouble right,GLdouble bottom,GLdouble top)

　　它是一个特殊的正射投影函数，主要用于二维图像到二维屏幕上的投影。

它的near和far缺省值分别为-1.0和1.0，全部二维物体的Z坐标都为0.0。因此它的裁剪面是一个左下角点为（left，bottom）、右上角点为（right，top）的矩形。

　　8.3.2 透视投影（Perspective Projection）
　　透视投影符合人们心理习惯，即离视点近的物体大，离视点远的物体小，远到极点即为消失，成为灭点。它的视景体类似于一个顶部和底部都被切除掉的棱椎。也就是棱台。这个投影通经常使用于动画、视觉仿真以及其他很多具有真实性反映的方面。
　　OpenGL透视投影函数也有两个，当中函数glFrustum()在8.1.3节中提到过。它所形成的视景体如图8-10所看到的。

 

图8-10 函数glFrustum()透视投影视景体

　　这个函数原型为：　　

void glFrustum(GLdouble left,GLdouble Right,GLdouble bottom,GLdouble top,GLdouble near,GLdouble far);

　　它创建一个透视视景体。其操作是创建一个透视投影矩阵。而且用这个矩阵乘以当前矩阵。这个函数的參数仅仅定义近裁剪平面的左下角点和右上角点的三维空间坐标，即（left。bottom，-near）和（right，top，-near）。最后一个參数far是远裁剪平面的Z负值。其左下角点和右上角点空间坐标由函数依据透视投影原理自己主动生成。near和far表示离视点的远近，它们总为正值。

　　还有一个函数是：　

void gluPerspective(GLdouble fovy,GLdouble aspect,GLdouble zNear, GLdouble zFar);

　　它也创建一个对称透视视景体。但它的參数定义于前面的不同，如图8-11所看到的。其操作是创建一个对称的透视投影矩阵，而且用这个矩阵乘以当前矩阵。

參数 fovy定义视野在X-Z平面的角度，范围是[0.0, 180.0]；參数aspect是投影平面宽度与高度的比率。參数zNear和Far各自是远近裁剪面沿Z负轴到视点的距离，它们总为正值。

图8-11 函数gluPerspective()透视投影视景体

　　以上两个函数缺省时。视点都在原点，视线沿Z轴指向负方向。二者的应用实例将在兴许章节中介绍。
8.4、裁剪变换
　　在OpenGL中，空间物体的三维裁剪变换包括两个部分：视景体裁剪和附加平面裁剪。视景体裁剪已经包括在投影变换里，前面已述。这里不再反复。

以下简单讲一下平面裁剪函数的使用方法。
　　除了视景体定义的六个裁剪平面（上、下、左、右、前、后）外，用户还可自己再定义一个或多个附加裁剪平面。以去掉场景中无关的目标，如图8-12所看到的。

图8-12 附加裁剪平面和视景体

　　附加平面裁剪函数为：

　void glClipPlane(GLenum plane,Const GLdouble *equation);

　　函数定义一个附加的裁剪平面。当中參数equation指向一个拥有四个系数值的数组。这四个系数各自是裁剪平面Ax+By+Cz+D=0的A、B、 C、D值。因此，由这四个系数就能确定一个裁剪平面。參数plane是GL_CLIP_PLANEi(i=0,1,...)。指定裁剪面号。

　　在调用附加裁剪函数之前。必须先启动glEnable(GL_CLIP_PLANEi)。使得当前所定义的裁剪平面有效；当不再调用某个附加裁剪平面时，可用glDisable(GL_CLIP_PLANEi)关闭对应的附加裁剪功能。

　　以下这个样例不仅说明了附加裁剪函数的使用方法，并且调用了gluPerspective()透视投影函数。读者能够细细体会当中的使用方法。例程例如以下：

　　例8-6 裁剪变换例程（clipball.c）

　　

#include "glos.h"
　　#include <GL/gl.h>
　　#include <GL/glu.h>
　　#include <GL/glaux.h>

　　void myinit(void);
　　void CALLBACK myReshape(GLsizei w, GLsizei h);
　　void CALLBACK display(void);

　　void CALLBACK display(void)
　　{
　　　　GLdouble eqn[4] = {1.0, 0.0, 0.0, 0.0};

　　　　glClear(GL_COLOR_BUFFER_BIT);

　　　　glColor3f (1.0, 0.0, 1.0);
　　　　glPushMatrix();
　　　　glTranslatef (0.0, 0.0, -5.0);

　　　　/* clip the left part of wire_sphere : x<0 */
　　　　glClipPlane (GL_CLIP_PLANE0, eqn);
　　　　glEnable (GL_CLIP_PLANE0);
　　　　glRotatef (-90.0, 1.0, 0.0, 0.0);
　　　　auxWireSphere(1.0);
　　　　glPopMatrix();
　　　　glFlush();
　　}

　　void myinit (void)
　　{
　　　　glShadeModel (GL_FLAT);
　　}

　　void CALLBACK myReshape(GLsizei w, GLsizei h)
　　{
　　　　glViewport(0, 0, w, h);
　　　　glMatrixMode(GL_PROJECTION);
　　　　glLoadIdentity();
　　　　gluPerspective(60.0, (GLfloat) w/(GLfloat) h, 1.0, 20.0);
　　　　glMatrixMode(GL_MODELVIEW);
　　}

　　void main(void)
　　{
　　　　auxInitDisplayMode (AUX_SINGLE | AUX_RGB);
　　　　auxInitPosition (0, 0, 500, 500);
　　　　auxInitWindow ("Arbitrary Clipping Planes");
　　　　myinit ();
　　　　auxReshapeFunc (myReshape);
　　　　auxMainLoop(display);
　　}

图8-13 剪取后的网状半球体

8.5、视口变换
　　在前面几节内容中已相继提到过视口变换。这一节将针对OpenGL来讲述视口变换的原理及其相关函数的使用方法。运用相机模拟方式，我们非常easy理解视口变换就是类似于照片的放大与缩小。在计算机图形学中，它的定义是将经过几何变换、投影变换和裁剪变换后的物体显示于屏幕窗体内指定的区域内。这个区域通常为矩形，称为视口。OpenGL中相关函数是：　　

glViewport(GLint x,GLint y,GLsizei width, GLsizei height);

　　这个函数定义一个视口。函数參数(x, y)是视口在屏幕窗体坐标系中的左下角点坐标。參数width和height各自是视口的宽度和高度。

缺省时，參数值即(0, 0, winWidth, winHeight) 指的是屏幕窗体的实际尺寸大小。

全部这些值都是以象素为单位，全为整型数。

　　注意：在实际应用中。视口的长宽比率总是等于视景体裁剪面的长宽比率。假设两个比率不相等，那么投影后的图像显示于视口内时会发生变形。如图8-14所看到的。另外，屏幕窗体的改变一般不明显影响视口的大小。因此，在调用这个函数时，最好实时检測窗体尺寸，及时修正视口的大小，保证视口内的图像能随窗体的变化而变化。且不变形。

图8-14 视景体到视口的映射

8.6 矩阵堆栈
　　学过计算机的人或许都知道这个使用频率极高的名词 — “堆栈”。顾名思义，堆栈指的是一个顶部打开底部封闭的柱状物体，通经常使用来存放经常使用的东西。这些东西从顶部依次放入，但取出时也仅仅能从顶部取出。即“先进后出，后进先出”。

在计算机中，它常指在内存中开辟的一块存放某些变量的连续区域。

因此，OpenGL的矩阵堆栈指的就是内存中专门用来存放矩阵数据的某块特殊区域。
　　实际上，在创建、装入、相乘模型变换和投影变换矩阵时。都已用到堆栈操作。一般说来。矩阵堆栈经常使用于构造具有继承性的模型，即由一些简单目标构成的复杂模型。

比如，一辆自行车就是由两个轮子、一个三角架及其他一些零部件构成的。它的继承性表如今当自行车往前走时。首先是前轮旋转，然后整个车身向前平移，接着是后轮旋转，然后整个车身向前平移，如此进行下去，这样自行车就往前走了。
　　矩阵堆栈对复杂模型运动过程中的多个变换操作之间的联系与独立十分有利。

由于全部矩阵操作函数如glLoadMatrix()、glMultMatrix()、 glLoadIdentity()等仅仅处理当前矩阵或堆栈顶部矩阵，这样堆栈中以下的其他矩阵就不受影响。堆栈操作函数有以下两个：

void glPushMatrix(void);
void glPopMatrix(void);

　　第一个函数表示将全部矩阵依次压入堆栈中，顶部矩阵是第二个矩阵的备份；压入的矩阵数不能太多，否则出错。

第二个函数表示弹出堆栈顶部的矩阵，令原第二个矩阵成为顶部矩阵。接受当前操作，故原顶部矩阵被破坏。当堆栈中仅存一个矩阵时。不能进行弹出操作，否则出错。由此看出，矩阵堆栈操作与压入矩阵的顺序刚好相反，编程时要特别注意矩阵操作的顺序。
　　为了更好地理解这两个函数。我们能够形象地觉得glPushMatrix()就是“记住自己在哪”，glPopMatrix()就是“返回自己原来所在地”。

请看以下一例：

　　例8-7 堆栈操作例程（arm.c）　

#include "glos.h"
　　#include <GL/gl.h>
　　#include <GL/glu.h>
　　#include <GL/glaux.h>

　　void myinit(void);
　　void drawPlane(void);
　　void CALLBACK elbowAdd (void);
　　void CALLBACK elbowSubtract (void);
　　void CALLBACK shoulderAdd (void);
　　void CALLBACK shoulderSubtract (void);
　　void CALLBACK display(void);
　　void CALLBACK myReshape(GLsizei w, GLsizei h);

　　static int shoulder = 0, elbow = 0;

　　void CALLBACK elbowAdd (void)
　　{
　　　　elbow = (elbow + 5) % 360;
　　}

　　void CALLBACK elbowSubtract (void)
　　{
　　　　elbow = (elbow - 5) % 360;
　　}

　　void CALLBACK shoulderAdd (void)
　　{
　　　　shoulder = (shoulder + 5) % 360;
　　}

　　void CALLBACK shoulderSubtract (void)
　　{
　　　　shoulder = (shoulder - 5) % 360;
　　}

　　void CALLBACK display(void)
　　{
　　　　glClear(GL_COLOR_BUFFER_BIT);
　　　　glColor3f(0.0, 1.0, 1.0);

　　　　glPushMatrix();
　　　　glTranslatef (-0.5, 0.0, 0.0);
　　　　glRotatef ((GLfloat)
　　　　shoulder, 0.0, 0.0, 1.0);
　　　　glTranslatef (1.0, 0.0, 0.0);
　　　　auxWireBox(2.0, 0.2, 0.5);

　　　　glTranslatef (1.0, 0.0, 0.0);
　　　　glRotatef ((GLfloat) elbow, 0.0, 0.0, 1.0);
　　　　glTranslatef (0.8, 0.0, 0.0);
　　　　auxWireBox(1.6, 0.2, 0.5);

　　　　glPopMatrix();
　　　　glFlush();
　　}

　　void myinit (void)
　　{
　　　　glShadeModel (GL_FLAT);
　　}

　　void CALLBACK myReshape(GLsizei w, GLsizei h)
　　{
　　　　glViewport(0, 0, w, h);
　　　　glMatrixMode(GL_PROJECTION);
　　　　glLoadIdentity();
　　　　gluPerspective(65.0, (GLfloat) w/(GLfloat) h, 1.0, 20.0);
　　　　glMatrixMode(GL_MODELVIEW);
　　　　glLoadIdentity(); glTranslatef (0.0, 0.0, -5.0); /* viewing transform */
　　}

　　void main(void)
　　{
　　　　auxInitDisplayMode (AUX_SINGLE | AUX_RGBA);
　　　　auxInitPosition (0, 0, 400, 400);
　　　　auxInitWindow ("Composite Modeling Transformations");

　　　　myinit ();

　　　　auxKeyFunc (AUX_LEFT, shoulderSubtract);
　　　　auxKeyFunc (AUX_RIGHT, shoulderAdd);
　　　　auxKeyFunc (AUX_UP, elbowAdd);
　　　　auxKeyFunc (AUX_DOWN, elbowSubtract);
　　　　auxReshapeFunc (myReshape);
　　　　auxMainLoop(display);
　　}

　　从以上例程能够看出，复杂的机械手臂是由两个简单的长方体根据一定的继承关系构成的，而这个继承关系是由矩阵堆栈的顺序决定的。

图8-15 简单机械手臂的符合运动