Bresenham画线算法

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/sinat_41104353/article/details/82858375

今天复习这个算法的时候，有一点没有立马反应过来，故特此记此篇博客作为备忘。

描点原理

原理就是描实距离D点最近的那个点，距离的判断通过比较d1和d2的大小来确定。d1大，说明距离 $(x_i + 1, y_i)$更远，要描实 $(x_i + 1, y_i + 1)$。相反同理。

需要说明的是此图将横线和竖线的交点看作像素点。

不妨假设斜率 k 在 0 到 1 之间，并且 $x_2>x_1$。
于是我们可以有上图这样一条直线。设在第 i 步已经确定第 i 个像素点是 $(x_i, y_i)$，现在看第 i+1 步如何确定第 i+1 个像素点的位置：
$d_1 = y - y_i = k(x_i + 1) + b -y_i$
$d_2 = y_i + 1 - y = (y_i + 1) - k(x_i + 1) - b$
$d_1 > d_2$，下一个像素点取 $(x_i + 1, y_i + 1)$
$d_1 < d_2$，下一个像素点取 $(x_i + 1, y_i)$
$d_1 = d_2$，取两个像素点中的任意一个

我们可以知道 $d_1 和d_2$ 大小的比较可以通过作差来得出。
$d_1 - d_2 = 2k(x_i + 1) - 2y_i + 2b - 1$，其中的 k 值为 $\dfrac{\Delta{y}}{\Delta{x}}$，其中 $\Delta{x}$ 为 $x_2 - x_1$， $\Delta{y}$ 为 $y_2 - y_1$。

用 $p_i$ 代替 $d_1 - d_2$

这里有个问题是，k 值是通过除法得出来的，我们知道计算机计算除法都是不精确的，所以这里，我们要做一点处理，来避免除法运算。这里，我们做的是两边同乘 $\Delta{x}$
$\Delta{x}(d_1 - d_2) = 2\Delta{y}(x_i + 1) - \Delta{x}·2y_i + \Delta{x}·(2b - 1)$

由于 $\Delta{x}$ 在我们假设的时候就是大于0的，于是 $d_1 - d_2$ 的符号并不会因此而改变，于是我们不妨假设 $p_i = \Delta{x}(d_1 - d_2)$。
$p_i = 2\Delta{y}(x_i + 1) - \Delta{x}·2y_i + \Delta{x}·(2b - 1)$
我们将其中不变的一部分提取出来，并用常数 c 来代替：
$p_i = 2\Delta{y}·x_i - 2\Delta{x}·y_i + c$
这样，我们就可以通过判断 $p_i$ 的正负来判断是描哪个点。

$p_i$ 递推

当然，我们仅仅知道一个 $p_i$ 是不够的，我们得进一步的递推下去：
$p_{i+1} = 2\Delta{y}·x_{i+1} - 2\Delta{x}·y_{i+1} + c$
这里的 $p_{i+1}$ 代表下下一个描点的判别式。

通过做一个减法，我们可以得出 $p_{i+1} 和 p_i$ 之间的数值关系：
$p_{i+1} - p_i = 2\Delta{y} - 2\Delta{x}(y_{i+1} - y_i)$

于是我们可以得出下面结论：

1. $p_i \geqslant 0$，应取 $y_{i+1} = y_i + 1, p_{i+1} = p_i + 2(\Delta{y} - \Delta{x})$
2. $p_i < 0$，应取 $y_{i+1} = y_i , p_{i+1} = p_i + 2\Delta{y}$

如何确定 $p_1$ 呢？

$x_1 = x_1$， $y_1 = \dfrac{\Delta{y}}{\Delta{x}}x_1 + b$ => $\Delta{x}·y_1 = \Delta{y}·x_1 + \Delta{x}·b$

在 $p_i$ 的公式中，不妨令 i = 1，有：
$p_i = 2\Delta{y}(x_i + 1) - \Delta{x}·2y_i + \Delta{x}·(2b - 1)$
$p_1 = 2\Delta{y}(x_1 + 1) - \Delta{x}·2y_1 + \Delta{x}·(2b - 1)$
将刚才推导出来的式子带入化简得：
$p_1 = 2\Delta{y} - \Delta{x}$

程序代码

void BresenhamLine(int x1, int y1, int x2, int y2)
{
	int x, y, dx, dy, p;  //dx和dy是Δx和Δy
	x = x1;
	y = y1;
	
	dx = x2 - x1;
	dy = y2 - y1;
	
	p = 2*dy - dx;
	for (; x <= x2; x++) {
		SetPixel(x, y);
		if (p > 0) {
			y++;
			p += 2*(dy-dx);
		} else {
			p += 2*dy;
		}
	}
}