前言
在数学上,理想的直线是没有宽度的,它是由无数个点构成的集合。对直线进行光栅化时,只能在显示器说给定的有限个像素组成的矩阵中,确定最佳逼近于该直线的一组像素,并且按扫描线顺序。
本节介绍绘制线宽为一个像素的直线的三个常用算法:数值微分,中点画线和Bresenham算法。
数值微分法
已知过端点 P0(x0, y0),P1(x1, y1) 的直线段 L(P0, P1);直线斜率为 k = (y1 - y0) / (x1 - x0)。
于是 yi+1 = kxi+1 + b。
于是,x每增加1,y就增加k。画点的时候还需要判断 int(y+0.5) 向下取整。
1 // 数值微分法,伪代码 2 void DDAline(int x0, int y0, int x1, int y1) { 3 int dx, dy, x=x0, y=y0; 4 double k; 5 dx = x0 - x1; dy = y0 - y1; 6 k = dy / dx; 7 draw(x, y) 8 while (x <= x1) { 9 x += 1; 10 y += k; 11 draw(x, int(y+0.5)); 12 } 13 }
效果图如下(0,1)(5, 4):
中点画线法
假设线段 F(M) = ax + by + c = 0。(a=y0-y1; b=x1-x0; c=x0y1-x1y0)。
中点画线法的思想就是对于点(xp,yp)的下一个点 M(xp+1,y+0.5),拿这个中点和实际点比较,如果实际点在中点上方(F(M) < 0),则取(xp+1,yp+1)为下一个点。如果实际点在中点下方(即中点代入直线方程的值大于0,F(M) >= 0),则取(xp+1, yp)。
为了加速计算,我们通常采用增量的方法。
假设从(xp,yp)开始画线,d的初值d0 = F(x0+1, y0+0.5)= F(x0,y0)+a+0.5。(d = a+0.5b)
- 若 d>=0 , 即中点代入原直线方程中的值大于0,即中点在目标直线上方,我们应该取下面的点(xp+1, yp)。判断下一个像素的位置时,应计算 d = F(xp+2,yp+0.5) = d + a,增量为 a。
- 若 d < 0,即中点在目标直线的下方,即取中点上面的点(xp+1,yp+1)。判断下一个像素的位置时,应计算 d = F(xp+2, yp+1.5) = d + a + b,增量为 a+b。
因为我们只需要知道d的正负,所以可以调整 d‘ = 2a + b。
1 void Midpoint(int x0, int y0, int x1, int y1) { 2 int a, b, d1, d2, d, x, y; 3 a = y0 - y1; 4 b = x1 - x0; 5 d = 2*a + b; 6 d1 = 2*a, d2 = 2 * (a+b); 7 x = x0, y = y0; 8 draw(x, y); 9 while (x < x1) { 10 if (d < 0) 11 {x++, y++, d+=d2;} 12 else 13 {x++; d+=d1;} 14 draw(x, y); 15 } 16 }
效果如下图(0,1)(5, 4):
Bresenham算法
Bresenham算法也是采用增量的方法,y每次累加k,超过0.5的中点位置就取上面的点,然后增量减一。如果我们增量的起始值设置为 d-0.5。就可以直接和0比较。
每次和0.5比较都比较麻烦,所以假设一个格子的宽度为2dx,高为2dx,那么一个 k 就是2dy。
所以如果我们设置起始增量为 -dx(-0.5*2dx)。每次移动单位为1的长度,增量就增加2dy,然后拿着这个增量和 0 比较,如果大于0,取上面的点,然后增量要减少 2dx。
1 void IntegerBresenham(int x0, int y0, int x1, int y1) { 2 int x, y, dx, dy; 3 dx = x0-x1, dy = y0 - y1, e = -dx; 4 x = x0, y = y0; 5 for(i=0; i<=dx; ++i) { 6 draw(x, y); 7 x++, e=e+2*dy; 8 if (e >= 0) {y++; e = e-2*dx;} 9 } 10 }
下面是效果图:
完整代码
1 import matplotlib.pyplot as plt 2 from matplotlib.ticker import MultipleLocator 3 4 # 数值微分法 5 def DDA(x0, y0, x1, y1): 6 dy = y0 - y1 7 dx = x0 - x1 8 k = dy / dx 9 x = x0 10 y = y0 11 12 a = [] 13 b = [] 14 while (x <= x1): 15 a.append(x) 16 x += 1 17 b.append(int(y+0.5)) 18 y = y+k 19 plt.scatter(a, b, color=‘r‘) 20 21 # 中点画线法 22 def midpoint(x0, y0, x1, y1): 23 a = y0 - y1 24 b = x1 - x0 25 d = 2 * a + b 26 d1 = 2 * a 27 d2 = 2 * (a + b) 28 print(a, b, d1, d2, d) 29 x = x0 30 y = y0 31 32 a = [x0] 33 b = [y0] 34 while x < x1: 35 if d < 0: 36 x += 1 37 y += 1 38 d += d2 39 else: 40 x += 1 41 d += d1 42 a.append(x) 43 b.append(y) 44 plt.scatter(a, b, color=‘r‘) 45 46 def Bresenham(x0, y0, x1, y1): 47 dx = x1 - x0 48 dy = y1 - y0 49 e = -dx 50 x = x0 51 y = y0 52 53 a = [] 54 b = [] 55 while (x <= x1): 56 a.append(x) 57 b.append(y) 58 x += 1 59 e = e + 2 * dy 60 if e >= 0: 61 y += 1 62 e = e - 2 * dx 63 plt.scatter(a, b, color=‘r‘) 64 65 x_target = [0, 5] 66 y_target = [1, 4] 67 68 ax = plt.subplot(111); 69 plt.plot(x_target, y_target) 70 ax.xaxis.grid(True, which=‘major‘) 71 ax.yaxis.grid(True, which=‘major‘) 72 ax.xaxis.set_major_locator(MultipleLocator(1)) 73 ax.yaxis.set_major_locator(MultipleLocator(1)) 74 75 # DDA(0, 1, 5, 4) 76 # midpoint(0, 1, 5, 4) 77 Bresenham(0, 1, 5, 4) 78 plt.show()
时间: 2024-10-12 12:36:29