OpenGL三点确定一个圆

OpenGL三点确定一个圆

 

 

//rorger,2011 //contact me:jhcyd112@163.com #include "windows.h" #include "glut.h" #include "math.h" class Point2{ public: double x, y ; void set(double x, double y){this->x=x ;this->y=y;} void set(Point2&p){x=p.x;y=p.y;} double operator *(Point2& rhn) { return this->x*rhn.x + this->y*rhn.y; } Point2 operator*(double rhn) { Point2 result ; result.x=this->x*rhn; result.y=this->y*rhn; return result; } double dotPorduct(Point2& rhn) { return this->x*rhn.x + this->y*rhn.y; } Point2 operator -(Point2& rhn) { Point2 result ; result.x=this->x-rhn.x; result.y=this->y-rhn.y; return result; } Point2 operator +(Point2& rhn) { Point2 result ; result.x=this->x+rhn.x; result.y=this->y+rhn.y; return result; } Point2 operator /(double rhn) { Point2 result ; result.x=this->x/rhn; result.y=this->y/rhn; return result; } double GetLength() { return sqrt(this->x*this->x+this->y*this->y); } }; typedef Point2 Vector2; int GetCircumCenter (Point2 A,Point2 B,Point2 C, Point2& S,double& radius) { Vector2 a,b,c,na; a=B-A; b=C-B; c=A-C; na.set(-a.y,a.x); S=A+ (a+na*((b*c)/(na*c))) / ((double)2.0); double temp = (b*c)/(na*c); radius=(a.GetLength())/2.0 * (sqrt( temp * temp + 1)); return 0; } void drawArc(double x,double y,double start_angle,double end_angle,double delta_angle,double radius,bool fill) { if (fill) { glBegin(GL_TRIANGLE_FAN); } else { glBegin(GL_LINE_STRIP); } for (double i=start_angle;i<=end_angle;i+=delta_angle) { double vx=x+radius * cos(i); double vy=y+radius*sin(i); glVertex2d(vx,vy); } glEnd(); } void drawCircle(double x, double y, double radius,bool fill) { double pi=acos(-1.0); drawArc(x,y,0,2*pi,pi/360,radius,fill); } void drawPie(double x,double y,double start_angle,double end_angle,double delta_angle,double radius,bool fill) { drawArc(x,y,start_angle,end_angle,delta_angle,radius,fill); if (fill) { glBegin(GL_TRIANGLE_FAN); } else { glBegin(GL_LINES); } glVertex2d(x,y); glVertex2d(x+radius*cos(start_angle),y+radius*sin(start_angle)); if (!fill) { glVertex2d(x,y); } glVertex2d(x+radius*cos(end_angle),y+radius*sin(end_angle)); glEnd(); } void testCircle(Point2 A,Point2 B,Point2 C, Point2& S,double& radius) { GetCircumCenter(A,B,C,S,radius); glColor3f(1.0,1.0,0.0); drawCircle(S.x,S.y,radius,false); glFlush(); glColor3f(0.8,1.0,0.5); glBegin(GL_LINE_STRIP); glVertex2d(A.x,A.y); glVertex2d(B.x,B.y); glVertex2d(C.x,C.y); glVertex2d(A.x,A.y); glEnd(); glEnable(GL_LINE_STIPPLE); glColor3f(0.6,0.5,1.0); glLineStipple(1,0x00e1); glBegin(GL_LINES); glVertex2d(S.x,S.y); glVertex2d(A.x,A.y); glVertex2d(S.x,S.y); glVertex2d(B.x,B.y); glVertex2d(S.x,S.y); glVertex2d(C.x,C.y); glEnd(); glDisable(GL_LINE_STIPPLE); glPointSize(6.0); glColor3f(0.0,1.0,0.0); glBegin(GL_POINTS); glVertex2d(A.x,A.y); glVertex2d(B.x,B.y); glVertex2d(C.x,C.y); glColor3f(1.0,0.5,0.5); glVertex2d(S.x,S.y); glEnd(); glFlush(); } void myDisplay(void) // plot the sinc function, using world coordinates { glColor3f(1.0,0.0,0.0); glBegin(GL_LINES); glVertex2d(-15.0,0); glVertex2d(15.0,0); glEnd(); glBegin(GL_LINES); glVertex2d(0,-15.0); glVertex2d(0,15.0); glEnd(); Point2 A,B,C; Point2 S ; double radius ; //test1,线段有交点 pass A.set(-5,0); B.set(5,0); C.set(0,9); radius=0; testCircle(A,B,C,S,radius); //test2,线段无交点，但所在直线有交点 pass // A.set(0,2); // B.set(5,0); // C.set(0,5); // testCircle(A,B,C,S,radius); } void myInit(void) { glClearColor(0.0, 0.0, 0.0, 0.0); glColor3f(0.0f, 0.0f, 0.0f); glPointSize(4.0); glMatrixMode(GL_PROJECTION); glLoadIdentity(); gluOrtho2D(-15, 15, -15, 15); glLineWidth(2.0); } void main(int argc, char** argv) { glutInit(&argc, argv); glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB); glutInitWindowSize(480,480); glutInitWindowPosition(100,150); glutCreateWindow("LineInterSectoin"); glutDisplayFunc(myDisplay); myInit( ); glutMainLoop( ); }

 

 

 

 

 

 

这是两张不同情况下的图片；

 

写好这个程序，突然意识到一个很重要的问题；

 

我把Point2 和Vector2混用了，真有点不伦不类，即使它们是多么的相似，我还是应把它们区分开来。

 

2011-04-08 22:32

I'm tired.