直线Bresenham算法

生成直线的Bresenham算法

    DDA算法由于在循环中涉及实型数据的加减运算，因此直线的生成速度较慢。由Bresenham在1965年提出。设直线从起点（x1, y1）到终点（x2, y2）。直线可表示为方程y=mx+b。Bresenham算法是一种基于误差判别式来生成直线的方法。

 

一，直线Bresenham算法描述：

    设(x1，y1)和(x2，y2)分别为所求直线的起点和终点坐标，由直线的微分方程

                  y2 – y1 / x2 – x1 = m =直线的斜率                 (1)

    可通过计算由x方向的增量△x引起y的改变来生成直线：

                  xi+1=xi+△x                                                （2）

                  yi+1=yi+△y=yi+△x·m                               （3）

    也可通过计算由y方向的增量△y引起x的改变来生成直线：

                  yi+1=yi+△y                                                 (4)

                  xi+1=xi+△x=xi+△y/m                                (5)

    它也是采用递推步进的办法，令每次最大变化方向的坐标步进一个象素，同时另一个方向的坐标依据误差判别式的符号来决定是否也要步进一个象素。

我们首先讨论m=△y/△x，当0≤m≤1且x1<x2时的Bresenham算法。从DDA直线算法可知这些条件成立时，公式(2)、(3)可写成：

                   xi+1=xi+1                                                 （6）

                  yi+1=yi+m                                                （7）

     有两种Bresenham算法思想，它们各自从不同角度介绍了Bresenham算法思想，得出的误差判别式都是一样的。

 

二、直线Bresenham算法思想之一：

    由于显示直线的象素点只能取整数值坐标，可以假设直线上第i个象素点坐标为(xi，yi)，它是直线上点(xi，yi)的最佳近似，并且xi=xi(假设m<1），如下图所示。那么，直线上下一个象素点的可能位置是(xi+1，yi)或(xi+1，yi+1)。

 

     由图中可以知道，在x=xi+1处，直线上点的y值是y=m(xi+1)+b，该点离象素点(xi+1，yi)和象素点(xi+1，yi+1)的距离分别是d1和d2：

     d1=y-yi=m(xi+1)+b-yi                        (8)

     d2=(yi+1)-y=(yi+1)-m(xi+1)-b                (9)

这两个距离差是

     d1-d2=2m(xi+1)-2yi＋2b-1                    (10)
我们来分析公式(10)：
    (1)当此值为正时，d1>d2，说明直线上理论点离(xi+1，yi+1)象素较近，下一个象素点应取(xi+1，yi+1)。
    (2)当此值为负时，d1<d2，说明直线上理论点离(xi+1，yi)象素较近，则下一个象素点应取(xi+1，yi)。
    (3)当此值为零时，说明直线上理论点离上、下两个象素点的距离相等，取哪个点都行，假设算法规定这种情况下取(xi+1，yi+1)作为下一个象素点。

因此只要利用(d1-d2)的符号就可以决定下一个象素点的选择。为此，我们进一步定义一个新的判别式：

     pi=△x×(d1-d2)=2△y·xi-2△x·yi+c         (11)
式(11)中的△x=(x2-x1)>0，因此pi与(d1-d2)有相同的符号；

     这里△y=y2-y1，m=△y/△x；c=2△y+△x(2b-1)。

     下面对式(11)作进一步处理，以便得出误差判别递推公式并消除常数c。

将式(11)中的下标i改写成i+1，得到：

     pi+1=2△y·xi+1-2△x·yi+1+c                (12)

将式(2-12)减去(2-11)，并利用xi+1=xi+1，可得：

     pi+1= pi+2△y-2△x·(yi+1-yi)               (13)

再假设直线的初始端点恰好是其象素点的坐标，即满足：

     y1=mx1+b                                    (14)

由式(11)和式(14)得到p1的初始值：

     p1=2△y-△x                                 (15)

这样，我们可利用误差判别变量，得到如下算法表示： 

     初始 p1=2△y-△x

     当pi≥0时：

       yi+1=yi+1,

       xi+1=xi+1，
       pi+1=pi+2(△y-△x)                          (16)

     否则：

       yi+1=yi,
       xi+1=xi+1，
       pi+1=pi+2△y 

从式(16)可以看出，第i+1步的判别变量pi+1仅与第i步的判别变量pi、直线的两个端点坐标分量差△x和△y有关，运算中只含有整数相加和乘2运算，而乘2可利用算术左移一位来完成，因此这个算法速度快并易于硬件实现。

 

三、直线Bresenham算法思想之二：

     由于象素坐标的整数性，数学点(xi，yi)与所取象素点(xi，yir)间会引起误差(εi)，当xi列上已用象素坐标(xi，yir)表示直线上的点(xi，yi)，下一直线点B(xi+1，yi+1)，是取象素点C(xi+1，yir )，还是D(xi＋1，y(i+1)r)呢？

设A为CD边的中点，正确的选择：

若B点在A点上方，选择D点； 否则，选C点。

 

 

用误差式描述为：

       ε(xi+1)=BC-AC=(yi+1-yir)-0.5                                        (8')

求递推公式：

       ε(xi+2)=(yi+2-y(i+1)r)-0.5 = yi+1+m-y(i+1)r-0.5          (9')

当ε(xi+1)≥0时，选D点，y(i+1)r = yir+1

       ε(xi+2)= yi+1+m-yir-1-0.5=ε(xi+1)+m-1                       (10')

当ε(xi+1)﹤0时，选C点，y(i+1)r = yir

      ε(xi+2)= yi+1+m－yir-0.5=ε(xi+1)+m                           (11')

初始时：

      ε(xs+1)=BC-AC=m-0.5                                                   (12')

 

 

为了运算中不含实型数，同时不影响不等式的判断，将方程两边同乘一正整数。

令方程两边同乘2·Δx，即d=2·Δx·ε，则：

     初始时：d = 2·Δy-Δx                                                      (13')

递推式：

     当d≥0时：

       {

          d=d+2·(Δy－Δx)；
          y++；
          x++；
       }                                                                                   (14')
     否则：

       {

          d=d+2·Δy；
          x++；
        }

 

四、直线Bresenham算法完善：

    现在我们修正(16)公式，以适应对任何方向及任何斜率线段的绘制。如下图所示，线段的方向可分为八种，从原点出发射向八个区。由线段按图中所示的区域位置可决定xi+1和yi+1的变换规律。

容易证明：当线段处于①、④、⑧、⑤区时，以|△x|和|△y|代替前面公式中的△x和△y，当线段处于②、③、⑥、⑦区时，将公式中的|△x|和|△y|对换，则上述两公式仍有效。

在线段起点区分线段方向

 

五、直线Bresenham算法特点：

     由于程序中不含实型数运算，因此速度快、效率高，是一种有效的画线算法。

 

六、实现代码：

/************************************************************************* Function Name: rtLine Function Action: 直线Bresenham算法实现 Input Parameter: image--图像数据 pt1--直线段的端点之一 pt2--直线段的端点之二 thickness--线段的粗细程度 color--线段的颜色 Output Parameter: image--加了直线的图像数据 Return Value: NULL Author : WangFei 2010/5/28 12:58 Remark: *************************************************************************/ void rtLine(RtImage* image, RtPoint pt1, RtPoint pt2, int thickness, unsigned char* color) { int x = 0, y = 0, stepx = 0, stepy = 0, p = 0, i = 0, j = 0, k = 0, change = 0/*表示斜率大于45度，交换了dx, dy*/; RtPoint dp; dp.x = pt2.x - pt1.x; dp.y = pt2.y - pt1.y; if (dp.x > 0) stepx = 1; if (dp.y > 0) stepy = 1; dp.x = abs(dp.x); dp.y = abs(dp.y); if (dp.y > dp.x)//斜率大于45度 { swap(&dp.x, &dp.y);//交换x,y change = 1; } //初始 p1=2△y-△x p = (dp.y << 1) - dp.x; x = pt1.x; y = pt1.y; for(i = 1;i <= dp.x; i++)//从开始点到结束点 { for (j = 0; j < image->nChannels; j++)//图像通道，保证每个通道都可以包含数据 { for (k = 0; k < thickness; k++)//直线的绘制宽度 { *(image->pData + (y * image->width + x) * image->nChannels + j) = color[j]; if (change == 0) { *(image->pData + (max(y - k, 0) * image->width + x) * image->nChannels + j) = color[j]; *(image->pData + (min(y + k, image->height - 1) * image->width + x) * image->nChannels + j) = color[j]; } else { *(image->pData + (y * image->width + max(x - k, 0)) * image->nChannels + j) = color[j]; *(image->pData + (y * image->width + min(x + k, image->width - 1)) * image->nChannels + j) = color[j]; } } } /* 当pi≥0时： yi+1=yi+1, xi+1=xi+1， pi+1=pi+2(△y-△x) 否则： yi+1=yi, xi+1=xi+1， pi+1=pi+2△y */ if(p >= 0) { if(change == 0) y = y + stepy; else x = x + stepx; p = p - (dp.x << 1); } if(change == 0) x = x + stepx; else y = y + stepy; p = p + (dp.y << 1); } }

 

详细知识：http://jpkc.lit.edu.cn/09/jsjtxx/zxjx/index.htm