【光栅图形学】直线和圆的扫描线算法（代码及推导过程）

直线和圆的扫描转化算法

给定两点$A(x_a,y_a),B(x_b,y_b),(x_a\le x_b)$,在像素整点上画出该线段。

我们不妨先假设$0\le k\le 1$，我们探讨如下三种算法：

数值微分算法DDA（Digital Differential Analyzer）

显然直线方程$AB:\frac{y-y_a}{x-x_a}=\frac{y_b-y_a}{x_b-x_a}$

转化为直线$y=kx+b(k=\frac{y_b-y_a}{x_b-x_a},b=y_a-k{x}_a=\frac{x_b{y}_a-x_a{y}_b}{x_b-x_a})$，

我们可以枚举$x=x_a,x_a+1,x_a+2,\dots ,x_b$得要相应的$y=kx+b$，计算量为一次加法一次乘法和一次取整运算。如何优化？

用增量思想代替原式中的乘法，递推计算$y_i=k{x}_i+b$，有

$y_{i+1}=k{x}_{i+1}+b=k(x_i+1)+b=(k{x}_i+b)+k=y_i+k$

我们得到算法的递推式：
$$P_i(x_i,y_i)=\{\begin{array}{rlr}& (x_a,y_a)& i=1\\ & (x_{i-1}+1,y_{i-1}+k)& i\ge 2\end{array}$$

完整的代码如下（包括$|k|\ge 1$的部分）

void DDA(int x1,int y1,int x2,int y2)
{
	int dx=std::abs(x1-x2),dy=std::abs(y1-y2);
	// if |k| greater than 1, enumerate x and calculate y=kx+b
	if(dx>dy)
	{
		// we assume that x1<=x2
		if(x1>x2) std::swap(x1,x2),std::swap(y1,y2);
		double k=double(y2-y1)/(x2-x1),y=y1;
		// x_{i+1}=x_i+1,y_{i+1}=y_i+k for x \in [x1+1,x2]
		for(int x=x1;x<=x2;x++,y+=k)
		{
			// draw current pixel
			draw_pixel(x,int(y));
		}
	}
	// if |k| less than 1, enumerate y
	else
	{
		// same as above
		if(y1>y2) std::swap(y1,y2),std::swap(x1,x2);
		double k=double(x2-x1)/(y2-y1),x=x1;
		for(int y=y1;y<=y2;y++,x+=k)
		{
			draw_pixel(int(x),y);
		}
	}
}

中点画线法

将直线转化为一般式$f(x,y)=Ax+By+C=0$，即 ( y a − y b