线段裁剪二：梁友栋Barsky算法

梁友栋-Barsky裁剪算法

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

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

　　我们知道，一条两端点为P₁（x₁，y₁）、P₂（x₂，y₂）的线段可以用参数方程形式表示：

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

0≤u≤1

（1）

　　式中，Δx=x₂-x₁，Δy=y₂-y₁，参数u在0～1之间取值，P（x，y）代表了该线段上的一个点，其值由参数u确定，由公式可知，当u=0时，该点为P₁（x₁，y₁），当u=1时，该点为P₂（x₂，y₂）。如果点P（x，y）位于由坐标（x_min，y_min）和（x_max，y_max）所确定的窗口内，那么下式成立：

xmin≤x1+ u·Δx≤xmax ymin≤y1+ u·Δy≤ymax

（2）

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

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

（3）

　　其中，p、q定义为：

p1=-Δx， q1=x1-xmin p2= Δx， q2=xmax-x1 p3=-Δy， q3=y1-ymin p4= Δy， q4=ymax-y1

（4）

　　从（4）式可以知道：任何平行于窗口某边界的直线，其p_k=0，k值对应于相应的边界（k=1，2，3，4对应于左、右、下、上边界）。如果还满足q_k<0，则线段完全在边界外，应舍弃该线段。如果p_k=0并且q_k≥0，则线段平行于窗口某边界并在窗口内，见图中所示。公式（4）式还告诉我们：

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

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

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

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

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

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

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

u = qk/pk

（5）

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

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

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

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

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

　　1、初始化线段交点的参数：u₁=0，u₂=1；

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

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

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

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

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

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

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

三、算法代码

bool _rtPruneLB(RtVector& vector, RtRect rect) { RtVector dest; bool flag = false; float u1 = 0, u2 =1; int p[4], q[4]; float r; p[0] = vector.sp.x - vector.ep.x; p[1] = vector.ep.x - vector.sp.x; p[2] = vector.sp.y - vector.ep.y; p[3] = -vector.sp.y + vector.ep.y; q[0] = vector.sp.x - rect.x; q[1] = rect.x + rect.w - vector.sp.x; q[2] = vector.sp.y - rect.y; q[3] = rect.y + rect.h - vector.sp.y; for(int i = 0; i < 4; i++) { r = (float)q[i] / (float)p[i]; if(p[i] < 0) { u1 = max(u1,r); if(u1 > u2) { flag = true; } } if(p[i] > 0) { u2 = min(u2, r); if(u1 > u2) { flag = true; } } if(p[i] == 0 && p[i] < 0) { flag = true; } } if ( flag ) { return; } dest.sp.x = vector.sp.x - u1 *(vector.sp.x - vector.ep.x); dest.sp.y = vector.sp.y - u1 *(vector.sp.y - vector.ep.y); dest.ep.x = vector.sp.x - u2 *(vector.sp.x - vector.ep.x); dest.ep.y = vector.sp.y - u2 *(vector.sp.y - vector.ep.y); vector =dest; }