Comp.Graphic 高手，帮帮忙。

ehyeo

谁会 Cubic Hermite Curve 在 C++ 里的 basic code? 帮帮忙。。

[ 本帖最后由 ehyeo 于 6-3-2007 09:55 PM 编辑 ]

ehyeo

没人救我吗？？   还是我来错论坛了？？

这个codes是两个endpoints的，如果我想要有6个点的话，要怎么modify呢？？

#include <iostream.h>
#include <stdlib.h>
#include <GL/glut.h>

float Geometry[4][3] = {
        { 10,10,0 },        //        Point1
        {-10,5,-2},//        Point2
        { 5,-5,0  },        //        Tangent1
        { 5,10,0  }        //        Tangent2
       
};
unsigned int LOD=20;

//------------------------------------------------------------        OnReshape()
//
void OnReshape(int w, int h)
{
        // prevent a division by zero when minimised
        if (h==0)
                h=1;

        // set the drawable region of the window
        glViewport(0,0,w,h);

        // set up the projection matrix
        glMatrixMode(GL_PROJECTION);
        glLoadIdentity();

        // just use a perspective projection
        gluPerspective(45,(float)w/h,0.1,100);

        // go back to modelview matrix so we can move the objects about
        glMatrixMode(GL_MODELVIEW);
        glLoadIdentity();
}

//------------------------------------------------------------        Draw()
//
void OnDraw() {

        // clear the screen & depth buffer
        glClear(GL_DEPTH_BUFFER_BIT|GL_COLOR_BUFFER_BIT);

        // clear the previous transform
        glLoadIdentity();

        // set the camera position
        gluLookAt(        1,10,30,        //        eye pos
                                0,5,0,        //        aim point
                                0,1,0);        //        up direction

        glColor3f(0.5,0.5,1);
        glBegin(GL_LINE_STRIP);

        // use the parametric time value 0 to 1
        for(int i=0;i!=LOD;++i) {

                float t = (float)i/(LOD-1);

                // calculate blending functions
                float b0 =  2*t*t*t - 3*t*t + 1;  //=(1-t)^2(1+2t) H0 in notes
                float b1 = -2*t*t*t + 3*t*t;      //=t^2(3-2t)  H3 in notes
                float b2 = t*t*t - 2*t*t + t;     //=t(1-t)^2   H1 in notes
                float b3 = t*t*t - t*t;           //=-t^2(1-t)  H2 in notes

                // calculate the x,y and z of the curve point
                float x = b0*Geometry[0][0] +
                                  b1*Geometry[1][0] +
                                  b2*Geometry[2][0] +
                                  b3*Geometry[3][0] ;

                float y = b0*Geometry[0][1] +
                                  b1*Geometry[1][1] +
                                  b2*Geometry[2][1] +
                                  b3*Geometry[3][1] ;

                float z = b0*Geometry[0][2] +
                                  b1*Geometry[1][2] +
                                  b2*Geometry[2][2] +
                                  b3*Geometry[3][2] ;

                // specify the point
                glVertex3f( x,y,z );
        }
        glEnd();


        glColor3f(0,1,0.5);
        glPointSize(3);
        glBegin(GL_POINTS);
                glVertex3fv( Geometry[0] );
                glVertex3fv( Geometry[1] );
        glEnd();

        glColor3f(0.5,1,0.5);
        glPushMatrix();
                glTranslatef(Geometry[0][0],Geometry[0][1],Geometry[0][2]);
                glBegin(GL_LINES);
                        glVertex3f( 0,0,0 );
                        glVertex3fv( Geometry[2] );
                glEnd();
        glPopMatrix();

        glPushMatrix();
                glTranslatef(Geometry[1][0],Geometry[1][1],Geometry[1][2]);
                glBegin(GL_LINES);
                        glVertex3f( 0,0,0 );
                        glVertex3fv( Geometry[3] );
                glEnd();
        glPopMatrix();

        // currently we've been drawing to the back buffer, we need
        // to swap the back buffer with the front one to make the image visible
        glutSwapBuffers();
}

//------------------------------------------------------------        OnInit()
//
void OnInit() {
        // enable depth testing
        glEnable(GL_DEPTH_TEST);
}

//------------------------------------------------------------        OnExit()
//
void OnExit() {
}

//------------------------------------------------------------        OnKeyPress()
//
void OnKeyPress(unsigned char key,int,int) {
        switch(key) {

        // increase the LOD
        case '+':
                ++LOD;
                break;

        // decrease the LOD
        case '-':
                --LOD;

                // have a minimum LOD value
                if (LOD<3)
                        LOD=3;
                break;
        default:
                break;
        }

        // ask glut to redraw the screen for us...
        glutPostRedisplay();
}

//------------------------------------------------------------        main()
//
int main(int argc,char** argv) {

        // initialise glut
        glutInit(&argc,argv);

        // request a depth buffer, RGBA display mode, and we want double buffering
        glutInitDisplayMode(GLUT_DEPTH|GLUT_RGBA|GLUT_DOUBLE);

        // set the initial window size
        glutInitWindowSize(640,480);

        // create the window
        glutCreateWindow("Hermite Curve: +/- to Change Level of Detail";

        // set the function to use to draw our scene
        glutDisplayFunc(OnDraw);

        // set the function to handle changes in screen size
        glutReshapeFunc(OnReshape);

        // set the function for the key presses
        glutKeyboardFunc(OnKeyPress);
       
        // run our custom initialisation
        OnInit();

        // set the function to be called when we exit
        atexit(OnExit);

        // this function runs a while loop to keep the program running.
        glutMainLoop();
        return 0;
}