Bresenham算法详解

Bresenham算法

一、思想

• 通过各行、各列像素中心构造一组虚拟网格线，按照直线起点到终点的顺序，计算直线与各垂直网格线的交点，然后根据误差项的符号确定该列象素中与此交点最近的象素。

二、推导

1. 假设每次$x$增加$1$，$y$的递增（减）量为$0$或$1$，它取决于实际直线与最近光栅网格点的距离，这个距离的最大误差为$0.5$。
2. 误差项$d$的初值$d_0＝0$，且
   $$d＝d＋k(0<k<1)$$
   一旦$d\ge 1$，就把它减去$1$，保证$d$的相对性，使$d$在$0$、$1$之间。
   
3. 推导结果
   $d_{i+1}=d_i+k$ $(d_0=0,k=\frac{\Delta y}{\Delta x},i=0,1,2,3\cdots ,0<k<1)$
   当$d>1$时，$d=d-1$；
   $$\{\begin{array}{ll}{x}_{i+1}& =x_i+1\\ y_{i+1}& =\{\begin{array}{ll}{y}_i+1& (d>0.5)\\ y_i& (d\le 0.5)\end{array}\end{array}$$

三、改进

把个算法的效率提高到整数加法上。

1. 改进1：令 $e=d-0.5$，即$e_{i+1}=e_i+k$，$e_0=-0.5$。如果$e>0$，$e=e-1$
   $$\{\begin{array}{l}{x}_{i+1}=x_i+1\\ y_{i+1}=\{\begin{array}{ll}{y}_i+1& (e>0)\\ y_i& (e\le 0)\end{array}\end{array}$$

2. 改进2：用$e\ast 2\ast △x$来替换$e$，即$e_{i+1}=e_i+2△y$，$e_0=-△x$。如果$e>0$，$e=e-2△x$
   $$\{\begin{array}{l}{x}_{i+1}=x_i+1\\ y_{i+1}=\{\begin{array}{ll}{y}_i+1& (e>0)\\ y_i& (e\le 0)\end{array}\end{array}$$

四、总结

1. 拓展
   （1）以上讨论的是$0\le k\le 1$的情况，即$0＜△y＜△x$的情况；
   （2）若是$1<k$，即$0＜△x＜△y$的情况，则需将$x$和$y$的位置交换；
   （3）若$△y＜0$或$△x＜0$时，要将算法中的$y＝y＋1$换成$y＝y－1$，$x＝x＋1$换成$x＝x－1$。

2. 算法步骤
   （1）输入直线的两端点$P_0(x_0,y_0)$和$P_1(x_1,y_1)$；
   （2）计算初始值$△x=x_1-x_0$、$△y=y_1-y_0$、$e=-△x$、$k=\frac{△y}{△x}$；
   （3）通过$k$值，判断$k$的情况；（此处默认$k$为$0\le k\le 1$）
   （4）绘制点$(x,y)$;
   （5）$e$更新为$e+2△y$，判断$e$的符号。若$e>0$，则$(x,y)$更新为$(x+1,y+1)$，同时将$e$更新为$e-2△x$，否则 $(x,y)$更新为$(x+1,y)$；
   （6）当直线没有画完时，重复步骤5和6。否则结束。