光栅图形学(一)——直线段的扫描转换算法

1. 数值微分法

  已知过端点$P_0(x_0,y_0)$，$P_1(x_1,y_1)$的直线段$L(P_0,P_1)$；直线斜率$k=\frac{y_1-y_0}{x_1-x_0}$。画线的过程为：从$x$的左端点$x_0$开始，向右端点步进，步长为1(像素)，按$y=kx+b$计算相应的$y$坐标，并取像素点$(x,round(y))$作为当前坐标。但这样做，计算每一个点需要做一个乘法、一个加法。设步长为$\Delta x$，有$x_{i+1}=x_i+\Delta x$，于是
$$y_{i+1}=k{x}_{i+1}+b=k{x}_i+k\Delta x+b=y_i+k\Delta x\text{\hspace{1em}(1.1)}$$
  当$\Delta x=1$时，有$y_{i+1}=y_i+k$。即$x$每递增1，$y$递增$k$(即直线的斜率)。这样，计算就由一个乘法和加法减少为一个加法。

• 算法程序：DDA画线算法程序

void DDALine(int x0, int y0, int x1, int y1, int color)
{
   
	int x;
	float dx, dy, y, k;
	dx = x1 - x0, dy = y1 - y0;
	k = dy/dx, y = y0;
	for(x = x0; x < x1; x++,y+=k)
	{
   
		drawpiexl(x, int(y + 0.5), color);
	}
}

注意： 用int(y + 0.5)取整的目的是为取离真正交单近的像素网格点作为光栅化的点。
  应当注意的是，上述算法仅适用于$k\le 1$的情形。在这种情况下，$x$每增加1，$y$最多增加1。当$|k|\ge 1$时，必须把$x,y$地位互换，$y$每增加1，$x$相应增加$1/k$。本节的其它算法也是如此。在这个算法中，$y$与$k$必须用浮点数表示，而且每一步都要对$y$进行四舍五入后取整，这使得该算法不利于硬件的实现。

2.中点画线法

  通过观察可以发现，画直线段的过程中，当前像素点$(x_p,y_p)$，下一个像素点有两种选择点$P_1(x_p+$