计算机图形学课程设计 多边形的裁剪算法

河南理工大学

万方科技学院

课程设计报告

2010 — 2011学年第二学期

课程名称计算机图形学

设计题目多边形裁剪算法

学生姓名孙晓芳

学号0816304009

专业班级计算机科学与技术10升指导教师侯守明

2011 年6 月29 日

目录

目录

目录 .............................................................................................................................................. I 第1章程序运行环境................................................................................. 错误！未定义书签。

1.1 程序运行环境的简单介绍............................................................... 错误！未定义书签。

1.2 程序运行环境的安装....................................................................... 错误！未定义书签。

1.3 多边形裁剪算法设计的内容........................................................................... 第2章直线裁剪和多边形裁剪的简单比较 (4)

2.1 直线裁剪的介绍 (4)

2.1.1 直线裁剪的基本原理………………………………………......................................

2.1.2 直线裁剪算法的分类以及和窗口交点参数值的计算……………………………..

2.2 多边形裁剪介绍 (9)

2.2.1 多边形裁剪的基本思想……………………………………………………………..

2.2.2 多边形和窗口相交的判定方法…………………………………………..

第3章多边形裁剪方法的详细介绍 (12)

3.1 Sutherland-Hodgman算法………………………………………………………………….

3.2 多边形裁剪算法的流程图 (12)

3.3多边形裁剪算法的实现 (13)

第4章代码的实现 (14)

第5章总结 (21)

参考文献 (22)

第1章程序的运行环境

1.1 程序运行环境的简单介绍

本次设计主要是运用了程序设计语言主要以C/C++语言为主，开发平台为Visual C++。现在Windows系统的主流编译环境有Visual C++，C++ Builder，Dev-C++等，它们都是支持OpenGL的。但这次设计主要选择Visual C++ 作为学习OpenGL的实验环境。

Microsoft Visual C++，（简称Visual C++、MSVC、VC++或VC）微软公司的C++开发工具，具有集成开发环境，可提供编辑C语言，C++以及C++/CLI等编程语言。Microsoft Visual C++是Microsoft公司推出的开发Win32环境程序，面向对象的可视化集成编程系统。它不但具有程序框架自动生成、灵活方便的类管理、代码编写和界面设计集成交互操作、可开发多种程序等优点，而且通过简单的设置就可使其生成的程序框架支持数据库接口、OLE2，WinSock网络、3D控制界面。

OpenGL作为盖茨设计的主要软件它与其他软件相比主要有以下几个优点

1）与C语言紧密结合：

2）强大的可移植性：

3、高性能的图形渲染：

1.2 程序安装及步骤

1．选择一个编译环境

这里我们选择Visual C++ 作为学习OpenGL的实验环境2．安装GLUT工具包

GLUT不是OpenGL所必须的，但它会给我们的学习带来一定的方便，推荐安装。Windows 环境下的GLUT下载地址：（大小约为150k）

https://www.360docs.net/doc/567829427.html,/resources/libraries/glut/glutdlls37beta.zip

Windows环境下安装GLUT的步骤：

1）将下载的压缩包解开，将得到5个文件

2）glut.h放到GL文件夹（VC6中一般是：C:\Program Files\Microsoft Visual Studio\VC98\Include\GL，VC2005中是：C:\Program Files\Microsoft Visual Studio 8\VC\Include，新建GL文件夹，再将glut.h放到GL文件夹中）。

3）glut.lib和glut32.lib放到静态函数库所在文件夹（VC6中一般是：C:\Program Files\Microsoft Visual Studio\VC98\Lib, VC2005中是：C:\Program Files\Microsoft Visual Studio 8\VC\Lib）。4）把解压得到的glut.dll和glut32.dll放到操作系统目录下面的system32文件夹内。（其路径为：C:\Windows\System32）

3．建立一个OpenGL工程

这里以VC为例：首先从开始->所有程序->Microsoft Visual C++ 6.0菜单中打开VC，也可单击文件：C:\Program Files\Microsoft Visual Studio\Visual C++6\Common\MSDev98\Bin\msdev.exe打开VC，在VC中选择File->New->Project，然后选择Win32 Console Application，输入一个工程名,设为A，然后按OK。在谈出的对话框左边点Application Settings，找到A Simple application并勾上，选择Finish。然后打开工程代码文件:A.cpp,将其内容替换为实验示范代码.

点击运行按钮就可以执行调试程序了。

1.3 多边形裁剪算法设计的内容和要求

这次设计的主要内容；1）理解多边形裁剪与直线段裁剪的区别；

2）掌握多边形的裁剪过程；

3）理解并掌握Sutherland-Hodgman算法的裁剪思想。

第2章直线裁剪和多边形裁剪的比较2.1 直线裁剪的介绍

2.1.1直线裁剪的基本原理

图1所示的为直线与窗口边界之间可能出现的几种关系。可以通过检查直线的两个端点是否在窗口之内确定如何对此直线裁剪。如果一直线的两个端点均在窗口边界之内（如图1中P5到P6的直线），则此直线应保留。如果一条直线的一个端点在窗口外（如P9）另一个点在窗口内（如P10），则应从直线与边界的交点（P9）处裁剪掉边界之外的线段。如果直线的两个端点均在边界外，则可分为两种情况：一种情况是该直线全部在窗口之外；另一种情况是直线穿过两个窗口边界。图中从P3到P4的直线属于前一种情况，应全部裁剪掉；从P7到P8的直线属于后一种情况，应保留P7到P8的线段，其余部分均裁剪掉。

图1直线相对干窗口边界的栽剪

直线裁剪算法应首先确定哪些直线全部保留或全部裁剪，剩下的即为部分裁剪的直线。对于部分裁剪的直线则首先要求出这些直线与窗口边界的交点，把从交点开始在边界外的部分裁剪掉。一个复杂的画面中可能包含有几千条直线，为了提高算法效率，加快裁剪速度，应当采用计算量较小的算法求直线与窗口边界的交点。

2.1.2直线裁剪算法的分类以及和窗口交点参数值的计算

直线段裁剪算法是复杂图形裁剪的基础。复杂的曲线可以通过折线段来近似，从而裁剪问题也可以化为直线段的裁剪问题。

主要的四种算法

直接求交算法

Cohen-Sutherland算法

中点算法

梁友栋－barskey算法

Cohen-Sutherland算法的大意是：对于每条线段P1P2，分为3种情况处理。

①若P1P2完全在窗口内，则显示该线段P1P2，简称“取”之。

②若P1P2明显在窗口外，则丢弃该线段，简称“弃”之。

③若线段既不满足“取”的条件，也不满足“弃”的条件，则把线段分为两段。其中一段完全在窗口外，可弃之。然后对另一段重复上述处理。

1．区域码及其建立

Cohen-Sutherland直线裁剪算法的核心是把所有直线的端点均分配一个表示

其相对位置的4位二进制代码。此代码称为区域码。区域码按照端点与窗口边界

的相对位置编码，即区域码的4位分别代表端点位于窗口的上、下、左、右。区

域码从右到左的各位所代表的坐标区如下所示：

位 4 3 2 1

坐标区上下右左

上述各位中某位为1，则表示点位于此坐标区。窗口周围各坐标区的区域码

如图2所示。由图2可见，位于窗中内的点，其区域码应为0000，位于窗口左

下方的点，其区域码应为0101，其余类推。

区域码各位的值可以通过对端点坐标（x，y）与窗口边界的比较求得。如果，则区域码的第一位为1，其余各位的确定与此相似。现在的

x

计算机语言都可以进行位操作，因此，可以通过以下步骤建立区域码：

①计算出端点坐标与窗口边界的差。

②按计算出的各个差的符号把区域码的相应位置为0或1，即

区域码的第一位置为（x－x wmin）的符号位；

区域码的第二位置为（x wmin-x）的符号位；

区域码的第三位置为（y-y wmin）的符号位；

图2 区域码

区域码的第四位置为（y wmin -y ）的符号位。 2．区域码裁剪算法

对所有直线的端点都建立了区域码之后，就可按区域码判断直线在窗口之内或窗口之外。这可分为如下几种情况：

①若一直线的两个端点的区域码均为0000则此直线在窗口边界之内，应子保留。

②若一直线的两个端点的区域码的同一位同时为1，则此直线全部在窗口边界之外，应子裁剪。例如，若一直线的一个端点的区域码为1001，另一个端点的区域码为0101，则此两端点的区域码的第一位均为1，说明此两端点均在窗口边界之左，因此，直线在窗口边界之外，应予裁剪。可用将直线两个端点的区域码进行与操作的方法，判断直线是否在窗口之外，若与操作的结果为0000则两端点的区域码任何位均不同时为1，此直线不一定被裁剪。

③以上两种情况之外的直线，有可能穿过窗口，也有可能不穿过窗口，如图87所示。图中所示的两条直线都不符合情况②的要求，但一条直线（P 1P 2）穿过窗口，另一直线（P 3P 4）不守过窗口。对这类直线可以进行如下处理：取窗口外的一个端点与窗口边界比较以确定可排除直线的哪一部分，然后，把直线剩下的部分与其他边界比较，这样一直到直线全部被排除或确定直线的哪一部分在窗口之内为止。可按“左、右、下、上”的次序建立检查直线端点与窗口边界关系的算法。

下面介绍对图3所示的两条直线进行处理的过程。从直线P 1P 2的下端点P 1开始，依次检查窗口的左、上、右及下边界，可发现此点在窗口之下（因为区域码的第三位为1）。然后求得直线与下边界的交点P 1排除线段P 1 P 1这样直线就缩短为P 1 P 2。因为P 2在边界之外，将此端点与各边界比较，可发现此端点在窗口上面。计算出交点P 2线段P 1P 2保留下来。这样就完成了对这条线的处理并开始处理下一直线P 3 P 4。端点P 3在窗口之左，可计算出交点3P '讲裁剪掉线段33P P '。检查43P P '的区域码，可发现直线的这一剩余部分在窗口之下故也可排除。

由于窗口边界均平行于坐标轴，所以直线与窗口边界的交点可以用直线方程的各参数很方便地求出，对于一条端点坐标为()11,y x 及()22,y x 的直线，其与一垂直边界的交点的y 坐标可以用下式计算：

()11x x m y y -+=

这里互的取值可取x 或max w x ，斜率m 可用公式()()1212x x y y m --=计算。同样地，与水平边界交点的x 坐标可用下式计算：

()m y y x x 11-+=

这里y 的值可取几min w y 或max w y 。

本算法的优点在于简单，易于实现。可以简单的描述为将直线在窗口左边的部分删去，按左，右，下，上的顺序依次进行，处理之后，剩余部分就是可见的了。在这个算法中求交点是很重要的，决定了算法的速度。另外，本算法对于其它形状的窗口未必同样有效。

特点：用编码方法可快速判断线段的完全可见和显然不可见。

Liang －Barsky 算法

Liang （梁友栋)-Barsky 算法又称为参数方程法。首先写出端点()11,y x 及

()22,y x 之间连线的参数方程如下：

()xu x u x x x x ?+=-+=1121 ()yu y u y y y y ?+=-+=1121

其中，1212,y y y x x x -=?-=?。参数u 可取0 1之间的值，坐标()y x ,表示此范围内的u 值定义的直线上的一个点。当u ＝0时，()()11,,y x y x =，直线的另一端u ＝0时，()()22,,y x y x =。

我们发现，如果直线上的某点处于()y x ,及()min min ,w w y x 及()max max ,w w y x 所定义的窗口之内，则满足以下条：

min w x ≤xu x ?+1≤max w x min w y ≤yu y ?+1≤max w y

这四个不等式可以表示为：

k up ≤k q ， k=1,2,3,4

其中，参数q p ,定义为；

,1x p ?-= m i n 11w x x q -= ,2x p ?= 1m a x 2x x q w -=

,3y p ?-= m i n 13w y y q -= ,4y p ?= 1m a x 4y y q w -=

按照直线与窗口边界的相对位置，可分为以下几种倩况。

①任何一条与窗口边界平行的直线，其0=k p ，此处k 值表示取哪一个边界（k =1,2,3及4，分别相应于左、右、下及上边界）。如果对某一k 值，还满足0

②如果0≠k q ，则情况较复杂。此时，可把直线按从()11,y x 到()22,y x 连线方向作为正向，将此直线无限延伸，同时把各窗口边界也无限延伸（见图3）；然后分以下两种情况讨论：

a 当0

内部。

b 当0>k q 时，情况相反。即直线的延伸线是由窗口边界延伸线的内部到外部。

以图3中的直线1L 为例，说明以上两种情况。表

1.1给出1L 的延伸线相对于4个边界延伸线的方向及k q 取值。

表1.1 直线方向与当0≠k q 时，可以用下式计算出一参数u ，此u 值应于直线延伸线与窗口边界k 的延伸线的交点：

k k p q u =

对于每条直线，可计算出一对u 值（1u 及2u ），由此两参数可判断直线是如何裁剪的。其中，1u 的值可由从窗口边界外发出进人窗口边界之内的直线（0k P ）所涉及的窗口边界确定，同样可计算出这些窗口边界的k r 值，取各k r 中的最小值即为2u 。如果21u u >，则此直线全部在窗口外而可排除：否则，此直线在窗口内，可由1u 及2u 计算出彼裁剪直线的端点。图3中标出直线1L 及2L 相应于1u 及2u 的与边界的交点。1L 的12u u >，则1L 彼裁剪，并可由1u 及2u 计算出裁剪点；而2L 的21u u >，则此直线全部彼排除。

此算法可用以下过程表示。此过程中开始把交点参数置为01=u 及12=u 。对每一窗口边界计算出相应的p 及q 值，然后用函数CliPtest 由参数p 及q 决定此直线能否彼排除或者决定交点参数是否要调整。当0>p 时，用参数r 修改1u 值：当0>p 时，用参数r 修改2u 值。如果最后的结果为21u u >，则排除此直线；否则，只有当新的值使直线缩短时，才修改相应的参数u 。当时且00<=q p ，可排除此直线，因为直线与边界平行且在边界之外。如果四个p 及q 均被检查而直线被排除，则被裁剪直线的端点可由1u 及2u 的值确定。

图3 直线与窗口边界的相对

2.2多边形裁剪介绍

2.2.1多边形裁剪的基本思想

这次多边形的裁剪算法主要用的是Sutherland-Hodgman算法其基本思想是：

Sutherland-Hodgman算法也叫逐边裁剪法，该算法是萨瑟兰德(I．E．Sutherland)和霍德曼(Hodgman)在1974年提出的。这种算法采用了分割处理、逐边裁剪的方法。

一，基本思想:

一次用窗口的一条边裁剪多边形。

考虑窗口的一条边以及延长线构成的裁剪线该线把平面分成两个部分:可见一侧；不可见一侧。多边形的各条边的两端点S、P。它们与裁剪线的位置关系只有四种

情况（1）仅输出1个顶点P；

情况（2）输出0个顶点；

情况（3）输出线段SP与裁剪线的1个交点I；

情况（4）输出线段SP与裁剪线的1个交点I和1个终点P

2.2.2 多边形和窗口相交的判定方法

二、算法实现：

1、已知：多边形顶点数组src，顶点个数n，

定义新多边形顶点数组dest。

2、赋初值：用变量flag来标识：

0表示在内侧，1表示在外侧。

3、对多边形的n条边进行处理，对当前点号的考虑为：0～n-1。

for(i=0；i

{

if(当前第i个顶点是否在边界内侧？)

{

if(flag!=0) /*前一个点在外侧吗？*/

{

flag=0；/*从外到内的情况，将标志置0,作为下一次循环的前一点标志*/

(dest + j) =求出交点；/*将交点dest放入新多边形*/

j++；

}

(dest + j)= (src + i)；/*将当前点srci放入新多边形*/

j++；

}

else

{

if(flag==0) /*前一个点在内侧吗？*/

{

flag=1；/*从内到外的情况，将标志置1,作为下一次循环的前一点标志*/

(dest + j) =求出交点；/*将交点dest放入新多边形*/

j++；

}

}

s= (src + i)；/*将当前点作为下次循环的前一点*/

}

P3＝C（－x, y），

P4＝C（－x,－y），

P5＝C（y，x），

P6＝C（－y，x），

P7＝C（y，－x），

P8＝C（－y，－x）。

这种性质称为八对称性。因此，只要扫描转换八分之一圆弧，就可以通过圆弧的八对称性得到整个圆。

2.2.2圆Bresenham算法基本原理

计算误差值Pi，以圆心为平面坐标轴原点，以Y轴正半轴向第一象限变化作基准x与y

的变化关系，每点的偏移量以误差值Pi为计算基量，Pi<0则yi+1＝yi，否则yi+1＝yi+1，算出y的变化量（从Y正半轴开始到1/8弧x变化快于y变化，所以x变化量为1）。Pi的递归式，当Pi小于0，Pi+1＝Pi+4xi+6，否则Pi＝Pi+1+4（xi+yi）+10。每画一点后根据中心对称原理在4个坐标轴同时向8个方向画点。直至到1/8弧相交（此时X等于Y）。

2.3用梁友栋-Barsky算法进行线段裁剪

2.3.1梁友栋-Barsky算法特点

梁友栋-Barsky算法只能应用于矩形窗口的情形。通常梁友栋-Barsky算法比Cohen－Sutherland算法效率更高，因为需要计算的交点数目减少了。更新参数u1、u2仅仅需要一次除法；线段与窗口边界的交点仅计算一次，就计算出u1、u2最后的值。相比之下，即使一条线段完全落在裁剪窗口之外，Cohen－Sutherland算法也要对它反复求交点，而且每次求交计算都需要做乘除法。

2.3.2梁友栋-Barsky算法实现

1、初始化线段交点的参数：u1=0，u2=1；

2、计算出各个裁剪边界的p、q值；

3、根据p、q来判断：是舍弃线段还是改变交点的参数。

（1）当p<0时，参数r用于更新u1；（u1=max{u1，…，rk}）

（2）当p>0时，参数r用于更新u2。（u2=min{u2，…，rk}）

（3）如果更新了u1或u2后，使u1>u2，则舍弃该线段。

（4）当p=0且q<0时，因为线段平行于边界并且位于边界之外，则舍弃该线段。

4、p、q的四个值经判断后，如果该线段未被舍弃，则裁剪线段的端点坐标由参数u1和u2的值决定。

第3章详细设计

3.1 Bresenham算法画直线

程序画出一条端点为（x1, y1）和（x2, y2）的直线。其中变量的含义是：P是误差；const1和const2，是误差的逐点变化量；inc是y的单位递变量，值为1或-1；tmp是用作象限变换时的临时变量。程序以判断|dx|>|dy|为分支，并分别将2a, 3a象限的直线和3b, 4b象限的直线变换到1a, 4a和2b, 1b方向去，以求得程序处理的简洁。

1、画点(x1, y2); dx=x2-x1; dy=y2-y1；计算误差初值P1=2dy-dx; i=1；

2、求直线的下一点位置：

xi+1=xi+1；

if Pi>0 则yi+1=yi+1；

否则yi+1=yi；

3、画点（xi+1, yi-1）；

4、求下一个误差Pi+1；

if Pi>0 则Pi+1=Pi+2dy-2dx;

否则Pi+1=Pi+2dy；

5、i=i+1; if i

3.2 Bresenham算法画圆

给出圆心的坐标(0,0)和半径Ｒ，求圆图像的最佳逼近点。

圆是中心对称的特殊图形，所以可以将圆八等分，则只须对八分之一圆孤求解，其它圆孤可以由对称变换得到，我们求的八分之一圆孤为(0,R)-(R√2,R√2)，可知最

大位移方向是x方向，x0=0,y0=R，每次对x自增，然后判断y是否减1，直到x>=y为止。误差量由Ｆ(x,y)=x^2+y^2-R^2给出。

先找递推关系，若当前d=F(x+1,y-0.5)>0，则y须减1，则下一d值为d=F(x+2,y-1.5)=(x+2)^2+(y-1.5)^2-R^2=(x+1)^2+(x-0.5)^2-R^2+2x+3-2y+2=d+2x-2y+5，若

当前d=F(x+1,y-0.5)<0，则y不变，只有x增1，则下一d值为d=F(x+2,y-0.5)=d+2x+3。

d的初值，d0=F(1,R-0.5)=1.25-R，则可以对d-0.25进行判断，因为递推关系中只有整数运算，所以d-0.25>0即d>0.25，这和d>0等价，所以d取初值1-R。

3.3用梁友栋-Barsky算法进行线段裁剪

设窗口为下列方程所确定：

x=XL x=XR (记XL≤x≤XR的区域为Δx)

y=YB y=YT (记YB≤y≤YT的区域为Δy)

则平面上任意线段P1P2在窗口内的部分为：P1P2∩Δx∩Δy

设通过P1P2的直线与窗口左、右、上、下边界的交点分别为L、R、B、T，则P1P2可见部分亦可表为：P1P2∩LR∩BT

当直线垂直时：xp1=xp2且XL≤xp1≤XR时，P1P2∩LR∩BT 等价于求P1P2∩BT；

当直线平行时：yp1=yp2且YB≤yp1≤YT时，P1P2∩LR∩BT 等价于求P1P2∩LR；

当直线既不垂直又不平行时：xp1！= xp2且yp1！= yp2

L、R、B、T都是确定的

P1P2的斜率k=( yp2-yp1)/ ( xp2-xp1)为非零或非奇异值。

xB= xp1 +(YB- yp1)/k； xT= xp1 +(YT- yp1)/k；

且 xB 0； xT

令：X1= max(xL，min(xp1，xp2)) X2= min(xR，max(xp1，xp2))

则P1P2∩LR为非空集合的充要条件是：X1≤X2（按照：max(x1 , w1) ≤ min(x2,w2)）此时线段P1P2与窗口在X轴方向的公共部分就是：［X1，X2］

而P1P2∩LR∩BT为非空集合的充要条件为：

max(X1, min(xB , xT))≤min(X2, max(xB , xT))

第4章功能实现

4.1 Bresenham算法画直线

核心代码：

void BresenhamLine(int x0,int y0,int x1,int y1) //Bresenham画线算法{

int dx,dy,x,y,d,d1,d2,inc,tmp;

dx=x1-x0;

dy=y1-y0;

if(dx*dy>=0) inc=1;

else inc=-1;

if(abs(dx)>abs(dy))

{

if(dx<0)

{

tmp=x0;x0=x1;x1=tmp;

tmp=y0;y0=y1;y1=tmp;

dx=-dx;dy=-dy;

}

d=2*dy-dx;

d1=2*dy;

d2=2*(dy-dx);

x=x0;

y=y0;

putpixel(x,y);

while(x

{

x++;

if(d<0)

d+=d1;

else

{

y+=inc;

d+=d2;

}

putpixel(x,y);

}

}

else

{

if(dy<0)

{

tmp=x0;x0=x1;x1=tmp; tmp=y0;y0=y1;y1=tmp; dx=-dx;dy=-dy;

}

d=2*dx-dy;

d1=2*dx;

d2=2*(dx-dy);

x=x0;

y=y0;

putpixel(x,y);

while(y

{

y++;

if(d<0)

d+=d1;

else

{

x+=inc;

d+=d2;

}

putpixel(x,y);

}

}

}

4.2用Bresenham算法画圆

核心代码：

void plotC(int x,int y,int xc,int yc)

{

putpixel(xc+x,yc+y);

putpixel(xc+x,yc-y);

putpixel(xc-x,yc+y);

putpixel(xc-x,yc-y);

putpixel(xc+y,yc+x);

putpixel(xc+y,yc-x);

putpixel(xc-y,yc+x);

putpixel(xc-y,yc-x);

}

void Bresenham_Circle(int xc,int yc,int r)//Bresenham函数画圆{

int x,y,d;

y=r;

d=3-2*r;

x=0;

while(x<=y)

{

plotC(x,y,xc,yc);

if(d<0)

d+=4*x+6;

else

{

d+=4*(x-y)+10;

y=y-1;

}

x=x+1;

}

}

4.3梁友栋-Barsky裁剪算法剪裁前：

剪裁后：

核心代码：

int clipTest(float p,float q,float *u1,float *u2) //裁剪算法{

int retVal;

retVal=1; //*retVal为标志变量，0：表示舍弃；1：表示可见。

float r;

if(p<0.0)

{

r=q/p;

if(r>*u2){

retVal=0;

accept=false;

}

else if(r>*u1)

*u1=r; //*u1取"进入"点的最大参数值

}

else if(p>0.0)

{

r=q/p;

if(r<*u1) {

retVal=0;

accept=false;

}

else if(r<*u2)

*u2=r; //*u2取"离开" 点的最小参数值

}

else if(q<0.0){ //*p=0，且q<0。平行于边界，而且在界外的线

{retVal=0;

accept=false;}

}

return retVal;

}

void clipLine(int xwmin,int ywmin,int xwmax,int ywmax,int *x1,int *y1,int *x2,int *y2) //裁剪

{

float u1=0.0,u2=1.0,dx=*x2-*x1,dy;

if(clipTest(-dx,*x1-xwmin,&u1,&u2))

if(clipTest(dx,xwmax-*x1,&u1,&u2))

{ dy=*y2-*y1;

if(clipTest(-dy,*y1-ywmin,&u1,&u2))

if(clipTest(dy,ywmax-*y1,&u1,&u2))

{

if(u2<1.0)

{

*x2=*x1+u2*dx; //通过u2求得裁剪后的p2端点

*y2=*y1+u2*dy;

}

if(u1>0.0)

{

*x1=*x1+u1*dx; //通过u1求得裁剪后的p1端点

*y1=*y1+u1*dy;

}

}

}

}

void keyboard(unsigned char key,int x,int y)//键盘事件，按键盘'C'键进行裁剪