直线裁剪算法研究
摘要：直线裁剪是计算机图形学中的一个重要技术，在对常见的直经线裁剪的算法分析的基础上，针对Cohen-Sutherland算法和Liang-Barsky算法进行了分析研究。并对两种算法了计算直线与窗口边界的交点时，进行了有效有比较。
关键词：裁剪；算法；Cohen-Sutherland；Liang－Barsky；
1 引言
直线是图形系统中使用最多的一个基本元素。所以对于直线段的裁剪算法是被研究最深入的一类算法，目前在矩形窗口的直线裁剪算法中，出现了许多有效的算法。其中比较著名的有：Cohen-Sutherland算法、中点分割算法、Liang-Barsky算法、Sobkow-Pospisil-Yang算法，及Nicholl-Lee-Ncholl算法等。
2 直线裁剪的基本原理
图1所示的为直线与窗口边界之间可能出现的几种关系。可以通过检查直线的两个端点是否在窗口之确定如何对此直线裁剪。如果一直线的两个端点均在窗口边界之（如图1中P5到P6的直线），则此直线应保留。如果一条直线的一个端点在窗口外（如P9）另一个点在窗口（如P10），则应从直线与边界的交点（P9）处裁剪掉边界之外的线段。如果直线的两个端点均在边界外，则可分为两种情况：一种情况是该直线全部在窗口之外；另一种情况是直线穿过两个窗口边界。图中从P3到P4的直线属于前一种情况，应全部裁剪掉；从P7到P8的直线属于后一种情况，应保留P7到P8的线段，其余部分均裁剪掉。
图1直线相对干窗口边界的栽剪
直线裁剪算法应首先确定哪些直线全部保留或全部裁剪，剩下的即为部分裁剪的直
线。对于部分裁剪的直线则首先要求出这些直线与窗口边界的交点，把从交点开始在边界外的部分裁剪掉。一个复杂的画面中可能包含有几千条直线，为了提高算法效率，加快裁剪速度，应当采用计算量较小的算法求直线与窗口边界的交点。
3 cohen－sutherland直线裁剪算法
Cohen-Sutherland算法的大意是：对于每条线段P1P2，分为3种情况处理。 ①若P1P2完全在窗口，则显示该线段P1P2，简称“取”之。 ②若P1P2明显在窗口外，则丢弃该线段，简称“弃”之。
③若线段既不满足“取”的条件，也不满足“弃”的条件，则把线段分为两段。其中一段完全在窗口外，可弃之。然后对另一段重复上述处理。
1．区域码及其建立
Cohen-Sutherland直线裁剪算法的核心是把所有直线的端点均分配一个表示其相对位置的4位二进制代码。此代码称为区域码。区域码按照端点与窗口边界的相对位置编码，即区域码的4位分别代表端点位于窗口的上、下、左、右。区域码从右到左的各位所代表的坐标区如下所示：
位 4 3 2 1 坐标区 上 下 右 左
上述各位中某位为1，则表示点位于此坐标区。窗口周围各坐标区的区域码如图2所示。由图2可见，位于窗中的点，其区域码应为0000，位于窗口左下方的点，其区域码应为0101，其余类推。
区域码各位的值可以通过对端点坐标（x，y）与窗口边界的比较求得。如果x ①计算出端点坐标与窗口边界的差。 ②按计算出的各个差的符号把区域码的相应位置为0或1，即 区域码的第一位置为（x－xwmin）的符号位； 区域码的第二位置为（xwmin-x）的符号位； 区域码的第三位置为（y-ywmin）的符号位； 区域码的第四位置为（ywmin-y）的符号位。 2．区域码裁剪算法 对所有直线的端点都建立了区域码之后，就可按区域码判断直线在窗口之或窗口之外。这可分为如下几种情况： ①若一直线的两个端点的区域码均为0000则此直线在窗口边界之，应子保留。 ②若一直线的两个端点的区域码的同一位同时为1，则此直线全部在窗口边界之外，应子裁剪。例如，若一直线的一个端点的区域码为1001，另一个端点的区域码为0101，则此两端点的区域码的第一位均为1，说明此两端点均在窗口边界之左，因此，直线在窗口边界之外，应予裁剪。可用将直线两个端点的区域码进行与操作的方法，判断直线是否 图2 区域码 在窗口之外，若与操作的结果为0000则两端点的区域码任何位均不同时为1，此直线不一定被裁剪。 ③以上两种情况之外的直线，有可能穿过窗口，也有可能不穿过窗口，如图87所示。图中所示的两条直线都不符合情况②的要求，但一条直线（P1P2）穿过窗口，另一直线（P3P4）不守过窗口。对这类直线可以进行如下处理：取窗口外的一个端点与窗口边界比较以确定可排除直线的哪一部分，然后，把直线剩下的部分与其他边界比较，这样一直到直线全部被排除或确定直线的哪一部分在窗口之为止。可按“左、右、下、上”的次序建立检查直线端点与窗口边界关系的算法。 下面介绍对图3所示的两条直线进行处理的过程。从直线P1P2的下端点P1开始，依次检查窗口的左、上、右及下边界，可发现此点在窗口之下（因为区域码的第三位为1）。然后求得直线与下边界的交点P1排除线段P1 P1这样直线就缩短为P1 P2。因为P2在边界之外，将此端点与各边界比较，可发现此端点在窗口上面。计算出交点P2线段P1P2保留下来。这样就完成了对这条线的处理并开始处理下一直线P3 P4。端点P3在窗口之左，可计算出交点P3?讲裁剪掉线段P3P3?。检查P3?P4的区域码，可发现直线的这一剩余部分在窗口之下故也可排除。 由于窗口边界均平行于坐标轴，所以直线与窗口边界的交点可以用直线方程的各参数很方便地求出，对于一条端点坐标为?x1,y1?及?x2,y2?的直线，其与一垂直边界的交点的y坐标可以用下式计算： y?y1?m?x?x1? 这里互的取值可取x或xwmax，斜率m可用公式m??y2?y1??x2?x1?计算。同样地，与水平边界交点的x坐标可用下式计算： x?x1??y?y1?m 这里y的值可取几ywmin或ywmax。 下面是区域码直线裁剪算法的C语言描述，其中每个端点的区域码用4元素布尔数组表示。 区域码直线裁剪算法代码 Void Line_Clipping(x1,y1,x2,y2,xw_xmin,yw_ymin,xw_max,yw_max) float x1,y1,x2,y2,xw_xmin,yw_ymin,xw_max,xw_max,yw_max; /*?x1,y1?和?x2,y2?是线段端点坐标*/ { Int draw,done; Int code1,code2,code; float x,y; Draw=FALSE; done=FALSE; code1=get_code(x1,y1,xw_xmin,yw_ymin,xw_max,yw_max); code2=get_code(x2,y2,xw_xmin,yw_ymin,xw_max,yw_max); while(! Done) { if (code1= =0 && code= =2) { draw=TRUE;done=TRUE; } else if (code1 & code2 !=0) done=TRUE; else { If (code1 !=0) code=code1; else code=code2; if (code&TOP!=0) { y=yw_ymax; x=x1+(y-y1)*(x2-x1)/(y2-y1); } else if (code&BOTTOM!=0) { y=yw_ymin; x=x1+(y-y1)*(x2-x1)/(y2-y1); } else if (code&RTGHT!=0) { x=xw_xmax; y=y1+(x-x1)*(y2-y1)/(x2-x1); } else if (code&LIFT!=0) { x=xw_xmin; y=y1+(x-x1)*(y2-y1)/(x2-x1); } if (code= =code1) { x1=x; y1=y; code1=get_code{x1,y1,xw_xmin,yw_ymin,xw_max,yw_max}; } else { x2=x;y2=y; code2=get_code{x2,y2,xw_xmin,yw_ymin,xw_max,yw_max}; } If (draw) Draw_Line(x1,y1,x2,y2); } int get_code(x,y,xw_xmin,yw_ymin,xw_max,yw_max) float x,y,xw_xmin,yw_ymin,xw_max,yw_max; { int code; code=0; if (y>yw_ymax) code| =TOP; else if (y 4 Liang－Barsky算法 Liang（梁友栋)-Barsky算法又称为参数方程法。首先写出端点?x1,y1?及?x2,y2?之间连线的参数方程如下： x?x1??x2?x1?u?x1??xu y?y1??y2?y1?u?y1??yu 其中，?x?x2?x1,?y?y2?y1。参数u可取0 1之间的值，坐标?x,y?表示此围的u值定义的直线上的一个点。当u＝0时，?x,y???x1,y1?，直线的另一端u＝0时， ?x,y???x2,y2?。 我们发现，如果直线上的某点处于?x,y?及?xwmin,ywmin?及?xwmax,ywmax?所定义的窗口之，则满足以下条： xwmin≤x1??xu≤xwmax ywmin≤y1??yu≤ywmax
直线裁剪算法研究（Cohen - Sutherland算法和Liang - Barsky算法）
