梁友栋-Barsky裁剪算法

    Cyrus和Beck用参数化方法提出了比Cohen-Sutherland更有效的算法。后来梁友栋和Barsky独立地提出了更快的参数化线段裁剪算法，也称为Liany-Barsky（LB）算法。

     一、梁友栋-Barsky裁剪算法思想：

　　我们知道，一条两端点为P1（x1，y1）、P2（x2，y2）的线段可以用参数方程形式表示：

x= x1+ u·（x2-x1）= x1+ u·Δx y= y1+ u·（y2-y1）= y1+ u·Δy

0≤u≤1

（3-9）

　　式中，Δx=x2-x1，Δy=y2-y1，参数u在0～1之间取值，P（x，y）代表了该线段上的一个点，其值由参数u确定，由公式可知，当u=0时，该点为P1（x1，y1），当u=1时，该点为P2（x2，y2）。如果点P（x，y）位于由坐标（xwmin，ywmin）和（xwmax，ywmax）所确定的窗口内，那么下式成立：

xwmin≤x1+ u·Δx≤xwmax ywmin≤y1+ u·Δy≤ywmax

（3-10）

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

u·pk ≤qk ， k=1，2，3，4

（3-11）

　　其中，p、q定义为：

p1=-Δx， q1=x1-xwmin p2=Δx， q2=xwmax-x1 p3=-Δy， q3=y1-ywmin p4=Δy， q4=ywmax-y1

（3-12）

　　从（3-12）式可以知道：任何平行于窗口某边界的直线，其pk=0（但并不是所有的Pk均为0，是存在pk=0的意思。平行于窗口某边界的图片，会出现 （p1&&p2）||（p3&&p4）=0的情况），k值对应于相应的边界（k=1，2，3，4对应于左、右、下、上边界）。如果还满足qk<0（默认x1为最左点?默认斜率大于0小于1？），则线段完全在边界外，应舍弃该线段。如果pk=0并且qk≥0，则线段平行于窗口某边界并在窗口内，见图中所示。公式（3-12）式还告诉我们：

　　1、当pk<0时，线段从裁剪边界延长线的外部延伸到内部；

　　2、当pk>0时，线段从裁剪边界延长线的内部延伸到外部；



　　例如，当Δx≥0时，对于左边界p1<0（p1=-Δx），线段从左边界的外部到内部；

　　　　 　　　　　　对于右边界p2>0（p2=Δx），线段从右边界的内部到外部。

　　　　　当Δy<0时，对于下边界p3>0（p3=-Δy），线段从下边界的内部到外部；

　　　　　　　　　　对于上边界p4<0（p4=Δy），线段从上边界的外部到内部。

　　　　　当pK≠0时，可以计算出参数u的值，它对应于无限延伸的直线与延伸的窗口边界k的交点，即：

（3-13）

　　对于每条直线，可以计算出参数u1和u2，该值定义了位于窗口内的线段部分：

　　1、u1的值由线段从外到内遇到的矩形边界所决定（pk<0），对这些边界计算rk=qk/pk，u1取0和各个r值之中的最大值。

　　2、u2的值由线段从内到外遇到的矩形边界所决定（pk>0），对这些边界计算rk=qk/pk，u2取0和各个r值之中的最小值。

　　3、如果u1>u2，则线段完全落在裁剪窗口之外，应当被舍弃；否则，被裁剪线段的端点可以由u1和u2计算出来。

　二、梁友栋-Barsky裁剪算法实现：

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

　　2、计算出各个裁剪边界的pk、qk值(k=1,2,3,4)；

　　3、[调用函数clipTest()]，在函数中根据p、q来判断：是舍弃线段还是改变交点的参数。（计算rk=qk/pk，根据(1)、(2)··· 执行）

　　　（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的值决定。

　　三、梁友栋-Barsky裁剪算法特点：

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

　　梁友栋-Barsky和Cohen－Sutherland算法还可以扩展为三维裁剪算法。

　　四、梁友栋-Barsky裁剪算法程序：

        MFC编写的应用程序（包含源代码）：http://download.csdn.net/detail/marcnuth/4938143

        梁友栋算法核心源代码：

#include "graphics.h"
int clipTest(float p,float q,float *u1,float *u2)
{
　int flag=1; /*flag为标志变量，0：表示舍弃；1：表示可见。*/
　float r;
　if(p<0.0)
　{
　　r=q/p;
　　if(r>*u2)
　　　flag=0;
　　else if(r>*u1)
　　　*u1=r; /*u1取"进入"点的最大参数值*/
　}
　else if(p>0.0)
　{
　　r=q/p;
　　if(r<*u1)
　　　flag=0;
　　else if(r<*u2)
　　　*u2=r; /*u2取"离开" 点的最小参数值*/
　}
　else if(q<0.0) /*p=0，且q<0。平行于边界，而且在界外的线*/
　　flag=0;
　return(flag);
}
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;
　　　　　}
　　　　　line(x1,y1,x2,y2); /*裁剪后的线段输出*/
　　　　}
　　}
}

本文来源： http://course.cug.edu.cn/cugThird/CGOL_NET/CLASS/course/3-2-3-a.htm