圆的扫描转换-中点Bresenham算法

圆的扫描转换是在屏幕像素点阵中确定最佳逼近于理想圆的像素点集的过程。圆的绘制可以使用简单方程画圆算法或极坐标画圆算法，但这些算法涉及开方运算或三角运算，效率很低。
仅包含加减运算的顺时针绘制1/8圆的中点Bresenham算法，根据对称性可以绘制整圆 。

默认的圆是圆心位于坐标系原点，半径为R的圆。
屏幕设备坐标系的原点位于左上角，绘制结果为1/4圆，需要进行圆心平移或使用自定义坐系可以绘制整圆。圆是椭圆的特例，使用椭圆中点Bresenham算法也可绘制。

算法原理

圆心在原点、半径为R的圆方程的隐函数表达式为：F(x,y)=x2+y2−R2=0

圆将平面划分成三个区域：对于圆上的点，F(x，y)＝0；对于圆外的点，F（x，y）＞0；对于圆内的点，F（x，y）＜0。
根据圆的对称性，可以用四条对称轴x＝0，y＝0，x＝y，x＝－y将圆分成8等份。只要绘制出第一象限内的1/8圆弧，根据对称性就可绘制出整圆，这称为八分法画圆算法。
假定第一象限内的任意点为P(x,y)，可以顺时针确定另外7个点：P(y,x)，P(y,-x)，P(x,-y)，P(-x,-y)，P(-y,-x)，P(-y, x)，P(-x,y)。

中点Bresenham算法要从(0,R) 顺时针确定最佳逼近于该段圆弧的像素点集。此段圆弧的斜率k处处满足|k|＜1，即|Δx|＞|Δy|，所以x方向为主位移方向，因此中点Bresenham算法的原理简化如下：x方向上每次加1，y方向上减不减1取决于中点误差项的值。
假定圆上当前点是Pi(xi,yi)
，下一像素只能在Pu(xi+1,yi)和Pd(xi+1,yi−1)中选取。Pu和Pd的中点为M(xi+1,yi−0.5)显然，若M点在理想圆弧的下方，则Pu点离圆弧近，选取Pu；否则应选取Pd。

构造中点误差项

中点误差项的递推公式

中点误差项的初始值

代码

double x,y,d;    
d=1.25-R;x=0;y=R;
for(x=0;x<=y;x++)
{
    CirclePoint(x,y,pDC);//调用八分法画圆子函数
    if (d<0)
        d+=2*x+3;
    else
    {
        d+=2*(x-y)+5;
        y--;
    } 
 }