計算機圖形學(xué)作業(yè)-給定圓心,半徑,弧度和起始點繪制圓弧

計算機圖形學(xué)第六組繪制圓弧中點畫圓算法假設(shè)x坐標(biāo)為的各像素點中，與該圓弧最近者已確定，為,那么，下一個與圓弧最近的像素只能是正右方的，或右下方的兩者之一。令M為和的中點，易知M的坐標(biāo)為。顯然，若M在圓內(nèi)，則離圓弧近，應(yīng)取為下一個像素；否則應(yīng)取。

對于圓上的點，F(xiàn)=0；對于圓內(nèi)的點，F(xiàn)<0；對于圓外的點，F(xiàn)>0；令函數(shù)判別式d：d的初始值為：在d≥0的情況下，取右下方像素，在d<0的情況下，取正右方像素，

為了消除在計算判別式初始值產(chǎn)生的浮點數(shù)運算，將用4d來代替d。下述程序使用中點畫圓算法繪制一個1/8圓弧。voidMidPointCircle(intr,intcolor){intx,y,dx=0;y=r;d=5-4r;CirclePoint(x,y,color);while(x<=y){if(d<=0)d+=8x+12;else{d+=8(x-y)+20;y--;}x++;CirclePoint(x,y,color);}}考慮到中點畫圓法畫的是單調(diào)區(qū)間，

我們做以下四種分析：

（我們只做逆時針繪制）

令函數(shù)領(lǐng)判別式為

如果di<0,說明M在圓內(nèi)，下一點?。▁i+1，yi）即xi+1=xi+1，

yi+1=yi

新的判別式為如果di>=0,說明M在圓外或圓上，下一點?。▁i+1，yi-1）即xi+1=xi+1，

yi+1=yi-1

新的判別式為初始去浮點令hi=di-0.25則如果hi<0.25，即hi<0，下一點?。▁i+1，yi）即xi+1=xi+1，

yi+1=yi

新的判別式為如果hi>=0.25,即hi>=0,說明M在圓外或圓上，下一點?。▁i+1，yi-1）即xi+1=xi+1，

yi+1=yi-1

新的判別式為voidMidPointCircle(intx0,inty0,intxc,intyc,intR,intcolor){intx,y,h;x=x0;y=y0;d=x*x+xc*xc-2x*xc+2x-2xc+y*y+yc*yc-y+yc-2*y*yc-R*R+1;CirclePoint(x,y,color);while(L(x,y)>1&&y>yc)//L(x,y)當(dāng)前點到終點的距離{if(h<0)h+=2*x-2xc+3;else

{h+=2*x-2xc-2*y+2*yc+5;

y--;

}x++;CirclePoint(x,y,color);}}令函數(shù)領(lǐng)判別式為

如果di<0,說明M在圓內(nèi)，下一點取（xi，yi-1）即xi+1=xi，

yi+1=yi-1,

新的判別式為如果di>=0,說明M在圓外或圓上，下一點取（xi-1，yi-1）即xi+1=xi-1，

yi+1=yi-1

新的判別式為初始去浮點令hi=di-0.25則如果hi<0.25，即hi<0，下一點?。▁i，yi-1）即

xi+1=xi，

yi+1=yi-1,

新的判別式為如果hi>=0.25,即hi>=0,說明M在圓外或圓上，下一點?。▁i-1，yi-1）即xi+1=xi-1，

yi+1=yi-1

新的判別式為voidMidPointCircle(intx0,inty0,intxc,intyc,intR,intcolor){intx,y,h;x=x0;y=y0;d=x*x+xc*xc-2x*xc-x+xc+y*y+yc*yc-2y+2yc-2*y*yc-R*R+1;CirclePoint(x,y,color);while(L(x,y)>1&&x>xc){if(h<0)h+=2*yc-2y+3;else

{h+=2*xc-2x-2*y+2*yc+5;

y--;

}x++;CirclePoint(x,y,color);}}令函數(shù)領(lǐng)判別式為

如果di<0,說明M在圓內(nèi)，下一點?。▁i-1，yi）即xi+1=xi-1，

yi+1=yi

新的判別式為如果di>=0,說明M在圓外或圓上，下一點?。▁i-1，yi+1）即xi+1=xi-1，

yi+1=yi+1

新的判別式為初始去浮點令hi=di-0.25則如果hi<0.25，即hi<0，下一點取（xi-1，yi）即

xi+1=xi-1，

yi+1=yi

新的判別式為如果hi>=0.25,即hi>=0,說明M在圓外或圓上，下一點取（xi-1，yi+1）即xi+1=xi-1，

yi+1=yi+1

新的判別式為voidMidPointCircle(intx0,inty0,intxc,intyc,intR,intcolor){intx,y,h;x=x0;y=y0;d=x*x+xc*xc-2x*xc+2x-2xc+y*y+yc*yc-y+yc-2*y*yc-R*R+1;CirclePoint(x,y,color);while(L(x,y)>1&&x>xc){if(h<0)h+=2*xc+2x+3;else

{h+=2*xc-2x-2*yc+2*y+5;

y--;

}x++;CirclePoint(x,y,color);}}令函數(shù)令判別式為

如果di<0,說明M在圓內(nèi)，下一點?。▁i，yi+1）即xi+1=xi，

yi+1=yi+1，

新的判別式為如果di>=0,說明M在圓外或圓上，下一點?。▁i+1，yi,+1）即xi+1=xi+1，

yi+1=yi+1

新的判別式為初始去浮點令hi=di-0.25則如果hi<0.25，即hi<0，下一點取（xi，yi+1）即

xi+1=xi，

yi+1=yi+1

新的判別式為如果hi>=0.25,即hi>=0,說明M在圓外或圓上，下一點?。▁i+1，yi+1）即xi+1=xi+1，

yi+1=yi+1

新的判別式為voidMidPointCircle(intx0,inty0,intxc,intyc,intR,intcolor){intx,y,h;x=x0;y=y0;d=x*x+xc*xc-2x*xc+x-xc+y*y+yc*yc+2y-2yc-2*y*yc-R*R+1;CirclePoint(x,y,color);while(L(x,y)>1&&y>yc){if(h<0)h+=2*y-2yc+3;else