Bresenham画圆算法 与中点画圆法

Bresenham画圆算法

不失一般性,假设圆的圆心位于坐标原点(如果圆心不在原点,可以通过坐标平移使其与原点重合),半径为R。以原点为圆心的圆C有四条对称轴:x = 0, y = 0, x = yx = -y。若已知圆弧上一点P1Cx, y),利用其对称性便可以得到关于四条对称轴的其它7个点,即:

  P2Cx,y),

  P3C(-x, y),

  P4C(-x,y),

  P5Cyx),

  P6C(-yx),

  P7Cy,-x),

  P8C(-y,-x)。

这种性质称为八对称性。

因此,只要扫描转换八分之一圆弧,就可以通过圆弧的八对称性得到整个圆。

 

 

 

Bresenham算法】

    简单图形的扫描转换常用算法是Bresenham算法。它的思想在于用误差量来衡量点选取的逼近程度。其过程如下:

以平面二维图形的扫描转换为例,设要画的图形方程为F(x, y)=0,要画的区域为[x0, x](不妨设x方向是最大位移方向,即△x > y),则F(x,y) 也是一个误差度量函数,我们拿离散的点值代入如果大于0则正向偏离,否则负向偏离,等于0的情况比较少,它表示的是不偏离即恰好与真实点重合。既然x是最大位移方向,那每次对x自增1,相应的y可以选择不增或增1(-1,具体问题具体分析),选择的方法就是d = F(x + 1, y ± 0.5)的正负情况进行判断从而选择y的值。

实际情况中还要考虑到浮点数的计算问题,因为基本的图形扫描转换算法最好能够硬件实现,所以摆脱浮点数是最好的,常用的方法是对d进行递推,而不是直接由F(x,y)给出(直接给出速度会慢)

 

【圆的扫描转换算法】

    以画圆为例,给出圆心的坐标(0, 0)和半径R,求圆图像的最佳逼近点。

圆是中心对称的特殊图形,所以可以将圆八等分,则只须对八分之一圆孤求解,其它圆孤可以由对称变换得到,我们求的八分之一圆孤为(0, R) -(R2,R2),可知最大位移方向是x方向,x0 = 0, y0 = R每次对x自增,然后判断y是否减1,直到x >= y为止(从点(0, R)到圆的八分之一处就有这种情况)。误差量由F(x, y) = x^2 + y^2 - R^2给出

先找递推关系,若当前d = F(x + 1, y - 0.5) > 0,则y须减1,则下一d值为

d =  F(x + 2, y - 1.5) = (x + 2)^2 + (y - 1.5)^2 - R^2 = (x + 1)^2 + (x - 0.5)^2 - R^2 + 2x + 3 - 2y + 2 = d + 2x - 2y + 5,若当前d = F(x + 1, y - 0.5) < 0,则y不变,只有x1,则下一d值为d = F(x + 2, y - 0.5) = d + 2x + 3

    d的初值,d0 = F(1, R - 0.5) = 1.25 - R,则可以对d - 0.25进行判断,因为递推关系中只有整数运算,所以d - 0.25 > 0d > 0.25,这和d > 0等价,所以d取初值1 - R

【C代码】

[仅以八分圆孤为例]

void roundMid(int x1,int y1,int R,CDC* pDC)

{

int x = 0,y = R;

int d = 1 - R;                   // 起点(0,R),下一点中点(1,R-0.5),d=1*1+(R-0.5)*(R-0.5)-R*R=1.25-R,

// d只参与整数运算,所以小数部分可省略

while(y > x)                     // y>x即第一象限的第区八分圆

{

pDC->SetPixel(x+x1,y+y1,RGB(255,0,0));         // 圆心(x1,y1),画点时直接相加平移,画区

pDC->SetPixel(y+x1,x+y1,RGB(255,0,0));         // 画区

pDC->SetPixel(-x+x1,y+y1,RGB(255,0,0));        // 画区

pDC->SetPixel(-y+x1,x+y1,RGB(255,0,0));        // 画区

pDC->SetPixel(-x+x1,-y+y1,RGB(255,0,0));       // 画区

pDC->SetPixel(-y+x1,-x+y1,RGB(255,0,0));       // 画区

pDC->SetPixel(x+x1,-y+y1,RGB(255,0,0));        // 画区

pDC->SetPixel(y+x1,-x+y1,RGB(255,0,0));        // 画区

if(d < 0)

d = d + 2 * x + 3;                        // d的变化

else

{

d = d + 2 * (x - y) + 5;                  // d <= 0,d的变化

y--;                                      // y坐标减

}

x++;                                           // x坐标加

}

}

 

http://home.srvee.com/space-1021-do-blog-id-946.html

http://hi.baidu.com/sink_cup/blog/item/5acac7166be8cc1d962b4319.html

http://www.cnblogs.com/phinecos/archive/2007/07/28/834407.html (比较详细)

 

 

  • 2
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值