计算机图形学学习(二) Bresenham画圆算法讲解及matlab实现

Bresenham画圆算法介绍

• 先讲解1/8圆：设圆的半径为r,考虑圆心在(0,0),从x=0、y=r开始的顺时针方向的1/8圆周的生成过程。

• x每次增加1个单位,从x=0开始,到x=y结束,即有x_i+1=x_i+1；另外，y_i+1则在y_i或者y_i-1中选择(如图)
  

• y_i+1的取值取决于d₁和d₂的大小
  y²= r²-(x_i+1)²
  d₁= y_i²-y²= y_i²-r²+(x_i+ 1)²
  d₂ = y²-(y_i-1)² = r²-(x_i+1)²- (y_i一1)²
  令p_i=d₁- d₂,代人d₁、d₂得
  p_i = 2(x_i+1)²+y_i²+(y_i- 1)² - 2r²

• p_i称为误差项。如果p_i<0,则下一个像素点y_i+1=y_i,
  否则y_i+1=y_i-1(其中p_i=0时选择y_i+1=y_i-1)

• p_i的递归式为:
  p_i+1 = p_i+4x_i+6+2(y_i² +1-y_i²)- 2(y_i+1-y_i)
  p_i的初值代x_i=0,y_i=r可得:
  p₀= 2(0+1)²+r² +(r- 1)2² = 3一2r

---

Bresenham画圆算法实现

根据上面的推导,圆生成算法思路如下:
(1)求误差初值,p₀=3- 2r,i=1,画点(0,r)
(2)求下一个点的y坐标,其中x_i+1=x_i+1,如果p_i<0则y_i+1=y_i,否则y_i+1=y_i-1
(3)画点(x_i+1,y_i+1)
(4)计算下一个误差,如果p_i<0则p