详述 Bresenham 线生成算法

给定直线的起始点和结束点，在显示设备上绘制该线段是计算机图形学中最基本的操作之一。本文将介绍的 Bresenham Line Drawing Algorithm 即为执行该操作的其中一个较为优化的算法。

计算机显示设备是由一个个像素组成的，该算法的目的即为“点亮”由起始像素点 $(x_0,x_0)$ 至结束像素点 $(x_1,y_1)$ 的线段上的像素点，当分辨率足够大，像素点足够小，看起来就会像一条直线。

举例：

为简化后续讨论，设置如下限定条件：

1. 起始点在结束点左下方；
2. $x_0<x_1$ 且 $y_0<=y_1$ ；
3. 线段的斜率 $0\le k\le 1$ ；
4. 对于区间 $[x_0,x_1]$ 内的所有 $x$ ，该位置仅存在一个像素点。

设置以上限定条件之后，我们可以得到一个十分粗糙的算法：

def easyway(x0,y0,x1,y1):
  m = (y2-y1)/(x2-x1)
  for x in range(x0, x1+1, 1):
    # find y (integer)
    y = round(mx+b)
    colorify(x,y)

该算法是可行的，但在循环的每一步都需要计算 $mx+b$，执行速度较慢，因此，我们需要寻找一个快速得到 $y$ 值的方法。

推导

由于线段的斜率在0到1之间，因此对于点 $(x_k,y_k)$ ，下个点只可能是 $(x_k+1,y_k)$ 或 $(x_k+1,y_k+1)$ ，其中 $k$ 表示横坐标 $x$ 移动的次数，如图所示。

因此，我们需要从 $y_k$ 和 $y_k+1$ 中获取离线段较近的值，假设直线方程为 $y=mx+b$ ，设 $y$ 为直线上实际的值，$d_1(k)$ 和 $d_2(k)$ 分别为距离 $y_k$ 和 $y_k+1$ 的距离，则有：

$$y=m(x_k+1)+b$$

$$d_1(k)=y-y_k$$

$$d_2(k)=(y_k+1)-y$$

比较 $d_1(k)$ 和 $d_2(k)$ 的值，即可判断下个点的具体位置：

1. 若 $d_1(k)\le d_2(k)$ ，认为下个点应为 ( x k + 1 , y k