椭圆算法
#include <GL/glut.h>
#include <stdlib.h>
#include <math.h>
#include <windows.h>
void init(void)
{//初始化函数,也可以不要
glClearColor (1.0, 1.0, 1.0, 0.0); // 指定清空颜色(背景色)为白色
glMatrixMode (GL_PROJECTION); //指定投影矩阵
gluOrtho2D (0.0, 600.0, 0.0, 600.0); //指定坐标系上显示的区域
}
class scrPt {
public:
GLint x, y;
};
void setPixel (GLint x, GLint y)
{
glBegin (GL_POINTS);
glVertex2i (x, y);
glEnd ( );
Sleep(10);
glFlush();
}
inline int round (const float a) { return int (a + 0.5); }
/* The following procedure accepts values for an ellipse
* center position and its semimajor and semiminor axes, then
* calculates ellipse positions using the midpoint algorithm.
*/
void ellipseMidpoint (int xCenter, int yCenter, int Rx, int Ry)
{
int Rx2 = Rx * Rx;
int Ry2 = Ry * Ry;
int twoRx2 = 2 * Rx2;
int twoRy2 = 2 * Ry2;
int p;
int x = 0;
int y = Ry;
int px = 0;
int py = twoRx2 * y;
void ellipsePlotPoints (int, int, int, int);
/* Plot the initial point in each quadrant. */
ellipsePlotPoints (xCenter, yCenter, x, y);
/* Region 1 */
p = round (Ry2 - (Rx2 * Ry) + (0.25 * Rx2));
while (px < py) {
x++;
px += twoRy2;
if (p < 0)
p += Ry2 + px;
else {
y--;
py -= twoRx2;
p += Ry2 + px - py;
}
ellipsePlotPoints (xCenter, yCenter, x, y);
}
/* Region 2 */
p = round (Ry2 * (x+0.5) * (x+0.5) + Rx2 * (y-1) * (y-1) - Rx2 * Ry2);
while (y > 0) {
y--;
py -= twoRx2;
if (p > 0)
p += Rx2 - py;
else {
x++;
px += twoRy2;
p += Rx2 - py + px;
}
ellipsePlotPoints (xCenter, yCenter, x, y);
}
}
void ellipsePlotPoints (int xCenter, int yCenter, int x, int y)
{
setPixel (xCenter + x, yCenter + y);
setPixel (xCenter - x, yCenter + y);
setPixel (xCenter + x, yCenter - y);
setPixel (xCenter - x, yCenter - y);
}
void myDisplay(void)
{// 图形绘制回调函数,通常把图形绘制代码写到这个函数里面
glClear(GL_COLOR_BUFFER_BIT);
glColor3f (0.0, 0.0, 1.0); // 指定前景色(当前绘制颜色)为蓝色
ellipseMidpoint (300, 300, 200, 100);
glEnd ( );
glFlush ( ); // 使绘制立即反映到屏幕上
}
void main(int argc, char *argv[])
{
glutInit(&argc, argv);
glutInitDisplayMode(GLUT_RGB | GLUT_SINGLE);
glutInitWindowPosition(0, 0);
glutInitWindowSize(600, 600);
glutCreateWindow("ellipseMidpoint");
init();
glutDisplayFunc(myDisplay);
glutMainLoop();
}