Circle Drawing Algorithm

Circle Drawing Algorithm
Where do we draw a circle???

Properties of a circle: A circle is defined as a set of points that are all the given distance (xc,yc). This distance relationship is expressed by the pythagorean theorem in Cartesian coordinates as (x xc)2 + (y yc) 2 = r2 We could use this equation to calculate the points on the circle circumference by stepping along x-axis in unit steps from xc-r to xc+r and calculate the corresponding y values at each position as y = yc +(- ) (r2 (xc x )2)1/2 This is not the best method: Considerable amount of computation Spacing between plotted pixels is not uniform
A Simple Circle Drawing Algorithm (cont)

y0 202 02 20
y1 202 12 20 y2 202 22 20
y19 202 192 6
y20 202 202 0
Polar co-ordinates for a circle

We could use polar coordinates r and , x = xc + r cos y = yc + r sin A fixed angular step size can be used to plot equally spaced points along the circumference A step size of 1/r can be used to set pixel positions to approximately 1 unit apart for a continuous boundary But, note that circle sections in adjacent octants within one quadrant are symmetric with respect to the 45 deg line dividing the two octants Thus we can generate all pixel positions around a circle by calculating just the points within the sector from x=0 to x=y This method is still computationally expensive
Bresenham to Midpoint
Bresenham's line algorithm for raster displays is adapted to circle generation by setting up decision parameters for finding the closest pixel to the circumference at each sampling step. Bresenham's circle algorithm avoids square-root calculations by comparing the squares of the pixel separation distances.
Midpoint Circle Algorithm

We will first calculate pixel positions for a circle centered around the origin (0,0). Then, each calculated position (x,y) is moved to its proper screen position by adding xc to x and yc to y Note that along the circle section from x=0 to x=y in the first octant, the slope of the curve varies from 0 to -1 Therefore, we can take unit steps in the positive x direction over this octant and use a decision parameter to determine which of the two possible y positions is closer to the circle path at each step. Positions in the other seven octants are then obtained by symmetry.
Circle function around the origin is given by fcircle(x,y) = x2 + y2 r2

Any point (x,y) on the boundary of the circle satisfies the equation fcircle(x,y) =0
Mid-Point Circle Algorithm (cont)

Assume that we have just plotted point (xk, yk) The next point is a choice between (xk+1, yk) and (xk+1, yk-1) We would like to choose the point that is nearest to the actual circle So how do we make this choice?
(xk, yk) (xk+1, yk)
(xk+1, yk-1)
yk
midpoint
yk
midpoint
yk-1
yk-1

6
3 1 2 3 4

6 M 5
3 1 2 3 4

6
M
5
3 1 2 3 4

For a point in the interior of the circle, the circle function is negative and for a point outside the circle, the function is positive Thus, fcircle(x,y) < 0 if (x,y) is inside the circle boundary fcircle(x,y) = 0 if (x,y) is on the circle boundary fcircle(x,y) > 0 if (x,y) is outside the circle boundary
yk
yk-1
x2+y2-r2=0 Midpoint between candidate pixels at sampling position xk+1 along a circular path
Midpoint
xk xk+1 xk+3

Assuming we have just plotted the pixel at (xk,yk) , we next need to determine whether the pixel at position (xk + 1, yk-1) is closer to the circle Our decision parameter is the circle function evaluated at the midpoint between these two pixels pk = fcircle (xk +1, yk-1/2) = (xk +1)2 + (yk -1/2)2 r2 If pk < 0 , this midpoint is inside the circle and the pixel on the scan line yk is closer to the circle boundary. Otherwise, the mid position is outside or on the circle boundary, and we select the pixel on the scan line yk-1

Successive decision parameters are obtained using incremental calculations pk+1 = fcircle(xk+1+1, yk+1-1/2) = [(xk+1)+1]2 + (yk+1 -1/2)2 r2 OR pk+1 = pk+2(xK+1) + (yK+12 yk2) (yk+1- yk)+1 Where yk+1 is either yk or yk-1, depending on the sign of pk Increments for obtaining pk+1: 2xk+1+1 if pk is negative 2xk+1+1-2yk+1 otherwise
Midpoint circle algorithm

Note that following can also be done incrementally: 2xk+1 = 2xk +2 2 yk+1 = 2yk 2 At the start position (0,r) , these two terms have the values 2 and 2r-2 respectively Initial decision parameter is obtained by evaluating the circle function at the start position (x0,y0) = (0,r) p0 = fcircle(1, r-1/2) = 1+ (r-1/2)2-r2 OR P0 = 5/4 -r If radius r is specified as an integer, we can round p0 to p0 = 1-r
The actual algorithm

1: Input radius r and circle center (xc,yc) and obtain the first point on the circumference of the circle centered on the origin as (x0,y0) = (0,r) 2: Calculate the initial value of the decision parameter as
P0 = 5/4 - r
3: At each xk position starting at k = 0 , perform the following test: If pk < 0 , the next point along the circle centered on (0,0) is (xk+1, yk) and pk+1 = pk + 2xk+1 + 1
The algorithm
Otherwise the next point along the circle is (xk+1, yk-1) and pk+1 = pk + 2xk+1 +1 -2yk+1 Where 2xk+1 = 2xk+2 and 2yk+1 = 2yk-2 4: Determine symmetry points in the other seven octants 5: Move each calculated pixel position (x,y) onto the circular path centered on (xc ,yc) and plot the coordinate values
x = x+ xc , y= y+ yc
6: Repeat steps 3 through 5 until x >= y
Mid-Point Circle Algorithm Example

To see the mid-point circle algorithm in action lets use it to draw a circle centred at (0,0) with radius 10
Mid-Point Circle Algorithm Summary

The key insights in the mid-point circle algorithm are: Eight-way symmetry can hugely reduce the work in drawing a circle Moving in unit steps along the x axis at each point along the circles edge we need to choose between two possible y coordinates
Blank Grid
Blank Grid
10 9 8 7 6 5 4
3
2 1 0 0 1 2 3 4 5 6 7 8 9 10
Blank Grid

Circle Drawing Algorithm

Uploaded by

Circle Drawing Algorithm

Uploaded by

Circle Drawing Algorithm

Where do we draw a circle???

A Simple Circle Drawing Algorithm (cont)

y19 202 192 6

y20 202 202 0

Polar co-ordinates for a circle

Midpoint Circle Algorithm

Circle function around the origin is given by fcircle(x,y) = x2 + y2 r2

Mid-Point Circle Algorithm (cont)

Mid-Point Circle Algorithm (cont)

Mid-Point Circle Algorithm (cont)

Mid-Point Circle Algorithm (cont)

Midpoint Circle Algorithm

Midpoint Circle Algorithm

Midpoint Circle Algorithm

Midpoint circle algorithm

The actual algorithm

Mid-Point Circle Algorithm Example

Mid-Point Circle Algorithm Summary

You might also like