《地形模型的三维可视化》程序设计(共11页).doc

《《地形模型的三维可视化》程序设计(共11页).doc》由会员分享，可在线阅读，更多相关《《地形模型的三维可视化》程序设计(共11页).doc（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上地形模型的三维可视化程序设计一、题目请用OpenGL图形库，编制程序实现地形模型的三维可视化，并进行地形模型的纹理设置，并可以通过键盘进行交互操作（移动和旋转）。 （注：给定的地形模型数据是已经建好的不规则三角网。）二、提交资料 1程序的详细步骤和所有相关数据；步骤：录入三角形数据获得地形的范围绘制三角网设置纹理坐标计算参考点的位置移动相关数据文件：数据文件DHS.dat三角网文件DHS.tri图片文件 TERRAIN.BMP2. 程序清单(包括程序说明)；#include stdafx.h#include #include #include #include #in

2、clude #include void CALLBACK elbowAdd (void);void CALLBACK elbowSubtract (void);void CALLBACK shoulderAdd (void);void CALLBACK shoulderSubtract (void);/点结构定义typedef struct long id; double x; double y; double z;POINTXYZ;/三角形结构定义typedef structlong id;/三角形号long p3;/三角形三个顶点的序号long xl3;/拓扑关系TRIANGLE;/顶点变

3、量POINTXYZ *pPt;/三角形变量TRIANGLE *pTri;/顶点总数int iTotalNum;/三角形总数int iTotalTriNum;/外围结构typedef struct vrtagBOXdouble minx;double miny;double minz;double maxx;double maxy;double maxz;vrBOX;/外围结构vrBOX _box;unsigned int m_nID;/存储纹理的索引号int m_nWidth;/纹理图片的宽度int m_nHeight;/纹理图片的高度float (*_ptexture)2;/存储三角网的纹理

4、数据/视点参考点的偏移量double lookx, looky, lookz;/移动的速度float _speed;/垂直方向视点和所要看实体的间的距离float _height;/视点的位置double _eyeposition3;/参考点的位置double _referencepoint3;/鼠标位移与象素之比的分母，即象素数int _Pixels;/设置纹理坐标void SetTextureCoord();/读取纹理数据bool LoadTriangleBMP( LPSTR szFileNameD);/向前移动void Move(bool is_forward);/获得鼠标在x，y方向的

5、相对位移量bool Rotate(int dx, int dy);/计算视点与参考点的位置void _CalcuEyePositionAndReferencePoint(int type,float amount);void myinit(void);/读取点坐标、三角形数据void ReadData();void Render();/获得外围void SetBOX();void CALLBACK myReshape(GLsizei w, GLsizei h);void CALLBACK display(void);/读取纹理数据bool LoadTriangleBMP(char* szFil

6、eName)AUX_RGBImageRec *pBitmap = NULL;if(!szFileName) return false;/装载数据并存储pBitmap = auxDIBImageLoad(szFileName);if(pBitmap = NULL) return false;glGenTextures(1, &m_nID);/ 绑定纹理到纹理数组，并进行初始化。glBindTexture(GL_TEXTURE_2D, m_nID);/控制纹理，映射到片元(fragment)时怎样对待纹理。glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP

7、_S, GL_REPEAT);glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_REPEAT);glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);/定义二维纹理glTexImage2D(GL_TEXTURE_2D, 0, 3, pBitmap-sizeX, pBitmap-sizeY, 0, GL_RGB, GL_UNSIGNED_

8、BYTE, pBitmap-data);m_nWidth = pBitmap-sizeX;m_nHeight = pBitmap-sizeY;if (pBitmap)delete pBitmap;return true;int main(int argc, char* argv)/* 设置双缓存和RGBA颜色模式*/ auxInitDisplayMode (AUX_DOUBLE|AUX_RGBA); auxInitPosition (0, 0, 1600, 800); auxInitWindow (作业); myinit(); auxReshapeFunc (myReshape); auxKe

9、yFunc (AUX_LEFT, shoulderSubtract); auxKeyFunc (AUX_RIGHT, shoulderAdd); auxKeyFunc (AUX_UP, elbowAdd); auxKeyFunc (AUX_DOWN, elbowSubtract); auxMainLoop(display);return 0;void CALLBACK elbowAdd (void) Move(true);void CALLBACK elbowSubtract (void) Move(false);void CALLBACK shoulderAdd (void)Rotate(5

10、0,0);void CALLBACK shoulderSubtract (void)Rotate(-50,0);/读取点坐标、三角形数据void ReadData()ifstream ifile(f:qqDHS.dat);ifile iTotalNum;pPt =new POINTXYZiTotalNum;/循环读入定顶点坐标数据;for (int i=0;ipPti.id pPti.x pPti.y pPti.z;ifile.close();ifstream ifile2(f:qqDHS.tri);ifile2 iTotalTriNum;pTri =new TRIANGLEiTotalTri

11、Num;/循环读入定顶点坐标数据;for ( i=0;ipTrii.id pTrii.p0 pTrii.p1 pTrii.p2;ifile2.close();SetBOX();/功 能：获得地形的范围void SetBOX()_box.minx = ;_box.miny = ;_box.minz = ;_box.maxx = -;_box.maxy = -;_box.maxz = -;/对所有的顶点循环for (int i = 0; i pPti.x) _box.minx = pPti.x;if(_box.miny pPti.y) _box.miny = pPti.y;if(_box.minz

12、 pPti.z) _box.minz = pPti.z;if(_box.maxx pPti.x) _box.maxx = pPti.x;if(_box.maxy pPti.y) _box.maxy = pPti.y;if(_box.maxz pPti.z) _box.maxz = pPti.z;void myinit(void)/读取点坐标、三角形数据 ReadData();char *_bmpfilename; /纹理图片的文件名_bmpfilename=f:qqTERRAIN.BMP; /terrain.bmp; /suban_1.bmp;/读取纹理数据LoadTriangleBMP(_b

13、mpfilename);/计算纹理坐标SetTextureCoord();/设置视点的一些初始参数lookx = 5.0f;looky = 0.0f;lookz = 0.0f;_Pixels= 3000;_height=2.50;_speed= 5.0f; _eyeposition0 = _box.minx -300; _eyeposition1 = _box.miny -300; _eyeposition2 = _box.maxz+50; _referencepoint0 = _box.maxx;_referencepoint1 = _box.maxy;_referencepoint2 =

14、_box.minz;void CALLBACK myReshape(GLsizei w, GLsizei h) glViewport(0, 0, w, h); glMatrixMode(GL_PROJECTION); glLoadIdentity(); gluPerspective(60.0, 1.0*(GLfloat)w/(GLfloat)h, 1.0, 90000.0); glMatrixMode(GL_MODELVIEW); glLoadIdentity();void CALLBACK display(void)glClearColor(1.0,1.0,1.0,0.0); glClear

15、(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT); glPushMatrix();/投影变换gluLookAt(_eyeposition0 , _eyeposition1 , _eyeposition2,_referencepoint0, _referencepoint1 , _referencepoint2, 0,0,1); glPushMatrix(); glTranslated (_box.maxx+300 ,_box.maxy+100 , _box.maxz+300); glColor3f (1.0, 0.0, 0.0); auxSolidSphe

16、re (50.0); glPopMatrix (); GLfloat mat_ambient= 1.0, 1.0, 1.0, 1.0 ; GLfloat mat_diffuse= 1.0, 1.0, 1.0, 1.0 ; GLfloat mat_specular = 1.0, 1.0, 1.0, 1.0 ; GLfloat mat_shininess = 50.0 ; GLfloat light0_diffuse= 1.0, 1.0, 1.0, 1.0;/定义光源的位置 GLfloat light0_position = _box.maxx,_box.maxy, _box.maxz-800,

17、1.0 ; GLfloat light1_ambient= 1.0, 0.0,0.0, 1.0 ; GLfloat light1_diffuse= 1.0, 0.0, 0.0, 1.0 ; GLfloat light1_specular = 1.0, 0.0, 0.0, 1.0 ; GLfloat light1_position = (_box.minx +_box.maxx)/2.0 ,(_box.miny +_box.maxy)/2.0 , _box.maxz+500, 1.0 ; GLfloat spot_direction = 0.0, 0.0, -1.0 ; glMaterialfv

18、(GL_FRONT, GL_AMBIENT, mat_ambient); glMaterialfv(GL_FRONT, GL_DIFFUSE, mat_diffuse); glMaterialfv(GL_FRONT, GL_SPECULAR, mat_specular); glMaterialfv(GL_FRONT, GL_SHININESS,mat_shininess);glLightfv(GL_LIGHT0, GL_AMBIENT, light0_diffuse);glLightfv(GL_LIGHT0, GL_DIFFUSE, light0_diffuse);glLightfv(GL_L

19、IGHT0, GL_POSITION,light0_position);glLightfv(GL_LIGHT1, GL_AMBIENT, light1_ambient);glLightfv(GL_LIGHT1, GL_DIFFUSE, light1_diffuse);glLightfv(GL_LIGHT1, GL_SPECULAR,light1_specular);glLightfv(GL_LIGHT1, GL_POSITION,light1_position);glLightf (GL_LIGHT1, GL_SPOT_CUTOFF, 10.0);glLightfv(GL_LIGHT1, GL

20、_SPOT_DIRECTION,spot_direction); glEnable(GL_LIGHT0);/ glEnable(GL_LIGHT1); glDepthFunc(GL_LESS); glEnable(GL_DEPTH_TEST);/启动光照 glEnable(GL_LIGHTING);/绘制地面模型Render();glDisable(GL_LIGHTING); glPopMatrix(); glFlush();auxSwapBuffers(); /功 能：绘制三角网void Render()glBindTexture(GL_TEXTURE_2D, m_nID); glEnabl

21、e(GL_TEXTURE_2D); /画三角形线面glBegin(GL_POLYGON); /GL_LINE_LOOP);/for( int i = 0; i iTotalTriNum; i+)glTexCoord2f(_ptexturepTrii.p00, _ptexturepTrii.p01 );glVertex3f(pPtpTrii.p0.x, pPtpTrii.p0.y, pPtpTrii.p0.z);glTexCoord2f(_ptexturepTrii.p10, _ptexturepTrii.p11 );glVertex3f(pPtpTrii.p1.x, pPtpTrii.p1.y

22、, pPtpTrii.p1.z);glTexCoord2f(_ptexturepTrii.p20, _ptexturepTrii.p21 );glVertex3f(pPtpTrii.p2.x, pPtpTrii.p2.y, pPtpTrii.p2.z);glEnd();glDisable(GL_TEXTURE_2D); /设置纹理坐标void SetTextureCoord()_ptexture = new floatiTotalNum2;float _min_max = fabs(_box.maxx -_box.minx);float _min_maxy = fabs(_box.maxy-_

23、box.miny);for (int i = 0; i iTotalNum; i+)_ptexturei0 = (pPti.x-_box.minx)/_min_max;_ptexturei1 = (pPti.y-_box.miny)/_min_maxy;/计算视点的位置和参考点的位置void _CalcuEyePositionAndReferencePoint(int type,float amount)float a;double _h;switch(type) /向前移动就执行这步case 1:_eyeposition0 += lookx * amount;_eyeposition1 +=

24、 looky * amount;break; /水平面内旋转就执行这步case 2:lookx = lookx*cos(amount) + looky*sin(amount);looky = looky*cos(amount) - lookx*sin(amount);a = 1/sqrt(lookx*lookx + looky*looky);lookx *= a;looky *= a;break; /计算参考点的位置 _referencepoint0 = _eyeposition0 + lookx; _referencepoint1 = _eyeposition1 + looky; _refe

25、rencepoint2 = _eyeposition2 + lookz;/向前移动void Move(bool is_forward)/使视点向前移动if(is_forward) _CalcuEyePositionAndReferencePoint(1, _speed);else _CalcuEyePositionAndReferencePoint(1, - _speed);/获得鼠标在x，y方向的相对位移量bool Rotate(int dx, int dy) /调用_CalcuEyePositionAndReferencePoint()函数完成旋转操作 _CalcuEyePositionAndReferencePoint(2, (float)dx/ _Pixels); return true;3程序运行结果。专心-专注-专业