第05部分_计算机图形学_圆、椭圆的扫描转换_060323.pptVIP

圆、椭圆的扫描转换 Ray ray@mail.buct.edu.cn 内 容 圆的扫描转换 中点画圆法 Bresenham画圆法 椭圆的扫描转换 圆的性质与中点画圆法 圆的特征 八对称性 中点画圆法 考虑中心在原点，半径为R的第二个8分圆，构造判别式（圆方程） 若 d0, 则取P1为下一象素，而且再下一象素的判别式为 若d=0, 则应取P2为下一象素，而且下一象素的判别式为 第 一个象素是（0,R），判别式d的初始值为 中点画圆法程序代码 MidPointCircle(int r int color) { int x,y; float d; x=0; y=r; d=1.25-r; circlepoints (x,y,color); //显示圆弧上的八个对称点 while(x=y) { if(d0) d+=2*x+3; else { d+=2*(x-y)+5; y--;} x++; circlepoints (x,y,color); } } 改进(1) 为了进一步提高算法的效率，可以将上面的算法中的浮点数改写成整数，仅用整数实现中点画圆法。 e=d-0.25代替d e0=1-R 改进(1) 的程序代码 MidPointCircle(int r int color) { int x,y,d; x=0; y=r; d=1-r; circlepoints (x,y,color); //显示圆弧上的八个对称点 while(x=y) { if(d0) d+=2*x+3; else { d+=2*(x-y)+5; y--;} x++; circlepoints (x,y,color); } } 改进(2) 判别式d的增量是x、y的线性函数，因此可引进增量的变量deltax和deltay，使算法中不出现乘法。 改进(2)后的程序代码 MidPointCircle(int r int color) { int x,y,d; x=0; y=r; deltax=3;deltay=2-r-r;d=1-r; circlepoints (x,y,color); //显示圆弧上的八个对称点 while(x=y) { if(d0) {d+=deltax; deltax+=2;} else { d+=deltax+deltay; deltax+=2; deltay+=2; y--;} x++; circlepoints (x,y,color); } } Bresenham画圆算法 考虑4分圆的情形。 △H=(x+1)2+y2-R2 △D=(x+1)2+(y-1)2-R2 △V=x2+(y-1)2-R2 分3种情况： H在圆外，D、V在圆内； D在圆上，H在圆外，V在圆内； H、D在圆外，V在圆内。 当△D0时，δHD≤0，则取H，否则取D； 当△D=0时，取D； 当△D0时，δDV≤0，则取D，否则取V； δHD=| △H |-| △D | δDV=| △D |-| △V | δHD和 δDV均可以通过△D 推算出来 可采用增量算法计算△’D ： 如果下一个象素是H： （x’,y’)=(x+1,y) △’D=((x+1)+1)2+ (y-1)2-R2=△D +2x+1 如果下一个象素是D： △’D= △D +2x-2y+2 如果下一个象素是V： △’D=△D –2y+1 Bresenham画圆法程序代码 Bresenham_Circle(r,color) int r, color; { int x,y,delta,delta1,delta2,direction; x=0;y=r; delta=2*(1-r); while(y=0){ circlepoints (x,y,color); if (delta0) {delta1=2*(delta+y)-1; if (delta1=0) direction=1; else direction=2;} else if (delta0) {delta2=2*(delta-x)-1; if(delta2=0) direction=2; else direction =3;} else direction=2; switch (direction) { case 1: x++;delta+=2*x+1;break; case 2:x++;y--;delta+=2*(x-y+1);break; case 3:y--;delta+=(-2*y+1);break; }/*switch*/ }/*while*/ }/*Bresenham_Ci