开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > OpenGL > 查看源码
LiFeng_Terrain.cpp
资源名称：虚拟地形建模.rar　[点击查看]
上传用户：dfjhuyju
上传日期：2013-03-13
资源大小：11035k
文件大小：6k
源码类别：

OpenGL

开发平台：

Visual C++

LiFeng_Terrain.cpp：源码内容
  1. #include "stdafx.h"
  2. #include "LiFeng.h"
  3. #include "LiFeng_Terrain.h"
  4. #include <fstream.h>
  6. LiFeng_TIN::LiFeng_TIN()
  7. {
  8. m_Point=NULL;
  9. m_Triangle=NULL;
  10. Point_Count=0;
  11. Triangle_Count=0;
  12. Max_X=0;
  13. Min_X=0;
  14. Max_Y=0;
  15. Min_Y=0;
  16. Max_Z=0;
  17. Min_Z=0;
  19. }
  21. LiFeng_TIN::~LiFeng_TIN()
  22. {
  23. delete [] m_Point;
  24. delete [] m_Triangle;
  25. }
  27. BOOL LiFeng_TIN::LiFeng_ReadTin(CString filename)
  28. {
  29. ifstream fin(filename);
  30. fin >>Point_Count>>Triangle_Count;
  31. m_Point=new Cvpoint [Point_Count];
  32. m_Triangle=new Triangle [Triangle_Count];
  35. for(int i=0;i<Point_Count;i++)
  36. {
  37. fin>>m_Point[i].x>>m_Point[i].y>>m_Point[i].z;
  39. }
  42. for(int j=0;j<Triangle_Count;j++)
  43. {
  44. fin>>m_Triangle[j].p1>>m_Triangle[j].p2>>m_Triangle[j].p3;
  46. }
  48. return true;
  50. }
  53. BOOL LiFeng_TIN::LiFeng_Limit()
  54. {
  55. Max_X=m_Point[0].x;
  56. Min_X=m_Point[0].x;
  57. Max_Y=m_Point[0].y;
  58. Min_Y=m_Point[0].y;
  59. Max_Z=m_Point[0].z;
  60. Min_Z=m_Point[0].z;
  62. for(int i=0;i<Point_Count;i++)
  63. {
  64. if(m_Point[i].x>Max_X)
  65. Max_X=m_Point[i].x;
  67. if(m_Point[i].x<Min_X)
  68. Min_X=m_Point[i].x;
  70. if(m_Point[i].y>Max_Y)
  71. Max_Y=m_Point[i].y;
  73. if(m_Point[i].y<Min_Y)
  74. Min_Y=m_Point[i].y;
  76. if(m_Point[i].z>Max_Z)
  77. Max_Z=m_Point[i].z;
  79. if(m_Point[i].z<Min_Z)
  80. Min_Z=m_Point[i].z;
  83. }
  84. return true;
  85. }
  87. BOOL LiFeng_TIN::LiFeng_Compile()
  88. {
  89. glNewList(1,GL_COMPILE);
  90. glBegin(GL_TRIANGLES);
  93. Cvpoint p[3];
  94. Cvpoint vect[3];
  95. float r[3];
  96. float s[3];
  97. float dx=Max_X-Min_X;
  98. float dy=Max_Y-Min_Y;
  100. glColor3f(1.0,0,0);
  101. for(int i=0;i<Triangle_Count;i++)
  102. {
  103. p[0].x=m_Point[m_Triangle[i].p1].x;
  104. p[0].y=m_Point[m_Triangle[i].p1].y;
  105. p[0].z=m_Point[m_Triangle[i].p1].z;
  107. p[1].x=m_Point[m_Triangle[i].p2].x;
  108. p[1].y=m_Point[m_Triangle[i].p2].y;
  109. p[1].z=m_Point[m_Triangle[i].p2].z;
  111. p[2].x=m_Point[m_Triangle[i].p3].x;
  112. p[2].y=m_Point[m_Triangle[i].p3].y;
  113. p[2].z=m_Point[m_Triangle[i].p3].z;
  115. vect[0].x=p[2].x-p[0].x;
  116. vect[0].y=p[2].y-p[0].y;
  117. vect[0].z=p[2].z-p[0].z;
  119. vect[1].x=p[2].x-p[1].x;
  120. vect[1].y=p[2].y-p[1].y;
  121. vect[1].z=p[2].z-p[1].z;
  123. vect[2].x=vect[0].y*vect[1].z-vect[0].z*vect[1].y;
  124. vect[2].y=vect[0].z*vect[1].x-vect[0].x*vect[1].z;
  125. vect[2].z=vect[0].x*vect[1].y-vect[0].y*vect[1].x;
  127. r[0]=p[0].x/dx;
  128. r[1]=p[1].x/dx;
  129. r[2]=p[2].x/dx;
  131. s[0]=p[0].y/dy;
  132. s[1]=p[1].y/dy;
  133. s[2]=p[2].y/dy;
  136. glNormal3f(vect[2].x,vect[2].z,vect[2].y);
  138. glTexCoord2f(r[0],s[0]);
  139. glVertex3f(p[0].x,p[0].z,p[0].y);
  141. glTexCoord2f(r[1],s[1]);
  142. glVertex3f(p[1].x,p[1].z,p[1].y);
  144. glTexCoord2f(r[2],s[2]);
  145. glVertex3f(p[2].x,p[2].z,p[2].y);
  146.   
  149. }
  150. glEnd();
  151. glEndList();
  152. return true;
  153. }
  155. BOOL LiFeng_TIN::LiFeng_DisplayTin()
  156. {
  157. glCallList(1);
  160. return true;
  161. }
  165. LiFeng_DEM::LiFeng_DEM()
  166. {
  167. DEMHEADER header;
  168. header.dem_Height=0;
  169. header.LeftDown_X=0;
  170. header.LeftDown_Y=0;
  171. header.Max_Height=0;
  172. header.Min_Height=0;
  173. header.Real_Distance=0;
  174. header.RightUp_X=0;
  175. header.RightUp_Y=0;
  176. header.X_Count=0;
  177. header.X_Interval=0;
  178. header.Y_Count=0;
  179. header.Y_Interval=0;
  181. Dem_Header=&header;
  182. }
  185. LiFeng_DEM::~LiFeng_DEM()
  186. {
  187. delete  Dem_Height;
  188. }
  191. BOOL LiFeng_DEM::LiFeng_ReadDem(CString filename)
  192. {
  195. CFile findem;
  196. findem.Open(filename,CFile::modeRead);
  198. findem.Read(&Dem_Header->LeftDown_X,8);
  199. findem.Read(&Dem_Header->LeftDown_Y,8);
  201. findem.Read(&Dem_Header->RightUp_X,8);
  202. findem.Read(&Dem_Header->RightUp_Y,8);
  204. findem.Read(&Dem_Header->X_Count,4);
  205. findem.Read(&Dem_Header->Y_Count,4);
  207. findem.Read(&Dem_Header->X_Interval,8);
  208. findem.Read(&Dem_Header->Y_Interval,8);
  210. findem.Read(&Dem_Header->Real_Distance,8);
  212. findem.Read(&Dem_Header->Max_Height,4);
  213. findem.Read(&Dem_Header->Min_Height,4);
  215. Dem_Height=new float [Dem_Header->X_Count*Dem_Header->Y_Count];
  217. findem.ReadHuge(Dem_Height,Dem_Header->X_Count*Dem_Header->Y_Count*4);
  219. findem.Close();
  221. return true;
  222. }
  225. BOOL LiFeng_DEM::LiFeng_Compile()
  226. {
  227. Cvpoint p[4];
  228. Cvpoint vect[5];
  229. float r[4];
  230. float s[4];
  232. float dx=Dem_Header->X_Interval*Dem_Header->Real_Distance*Dem_Header->X_Count;
  233. float dy=Dem_Header->Y_Interval*Dem_Header->Real_Distance*Dem_Header->Y_Count;
  235. glNewList(1,GL_COMPILE);
  236. glBegin(GL_TRIANGLES);
  239. glColor3f(1.0,0,0);
  240. for(int j=0;j<Dem_Header->Y_Count-1;j++)
  241. {
  242. for(int i=0;i<Dem_Header->X_Count-1;i++)
  243. {
  244. p[0].x=i*Dem_Header->X_Interval*Dem_Header->Real_Distance;
  245. p[0].y=Dem_Height[j*Dem_Header->X_Count+i];
  246. p[0].z=j*Dem_Header->Y_Interval*Dem_Header->Real_Distance;
  248. p[1].x=(i+1)*Dem_Header->X_Interval*Dem_Header->Real_Distance;
  249. p[1].y=Dem_Height[j*Dem_Header->X_Count+i+1];
  250. p[1].z=j*Dem_Header->Y_Interval*Dem_Header->Real_Distance;
  252. p[2].x=(i+1)*Dem_Header->X_Interval*Dem_Header->Real_Distance;
  253. p[2].y=Dem_Height[(j+1)*Dem_Header->X_Count+i+1];
  254. p[2].z=(j+1)*Dem_Header->Y_Interval*Dem_Header->Real_Distance;
  256. p[3].x=i*Dem_Header->X_Interval*Dem_Header->Real_Distance;
  257. p[3].y=Dem_Height[(j+1)*Dem_Header->X_Count+i];
  258. p[3].z=(j+1)*Dem_Header->Y_Interval*Dem_Header->Real_Distance;
  262. vect[0].x=p[2].x-p[0].x;
  263. vect[0].y=p[2].y-p[0].y;
  264. vect[0].z=p[2].z-p[0].z;
  266. vect[1].x=p[1].x-p[0].x;
  267. vect[1].y=p[1].y-p[0].y;
  268. vect[1].z=p[1].z-p[0].z;
  270. vect[3].x=p[3].x-p[0].x;
  271. vect[3].y=p[3].y-p[0].y;
  272. vect[3].z=p[3].z-p[0].z;
  275. vect[2].x=vect[0].y*vect[1].z-vect[0].z*vect[1].y;
  276. vect[2].y=vect[0].z*vect[1].x-vect[0].x*vect[1].z;
  277. vect[2].z=vect[0].x*vect[1].y-vect[0].y*vect[1].x;
  280. vect[4].x=vect[0].y*vect[3].z-vect[0].z*vect[1].y;
  281. vect[4].y=vect[0].z*vect[3].x-vect[0].x*vect[1].z;
  282. vect[4].z=vect[0].x*vect[3].y-vect[0].y*vect[1].x;
  285. r[0]=p[0].x/dx;
  286. r[1]=p[1].x/dx;
  287. r[2]=p[2].x/dx;
  288. r[3]=p[3].x/dx;
  290. s[0]=p[0].y/dy;
  291. s[1]=p[1].y/dy;
  292. s[2]=p[2].y/dy;
  293. s[3]=p[3].y/dy;
  296. glNormal3f(vect[2].x,vect[2].z,vect[2].y);
  298. glTexCoord2f(r[0],s[0]);
  299. glVertex3f(p[0].x,p[0].z,p[0].y);
  301. glTexCoord2f(r[1],s[1]);
  302. glVertex3f(p[1].x,p[1].z,p[1].y);
  304. glTexCoord2f(r[2],s[2]);
  305. glVertex3f(p[2].x,p[2].z,p[2].y);
  308. glNormal3f(vect[4].x,vect[4].z,vect[4].y);
  310. glTexCoord2f(r[0],s[0]);
  311. glVertex3f(p[0].x,p[0].z,p[0].y);
  313. glTexCoord2f(r[2],s[2]);
  314. glVertex3f(p[2].x,p[2].z,p[2].y);
  316. glTexCoord2f(r[3],s[3]);
  317. glVertex3f(p[3].x,p[3].z,p[3].y);
  319. }
  320. }
  321. glEnd();
  322. glEndList();
  323. return true;
  325. }
  327. BOOL LiFeng_DEM::LiFeng_DisplayDem()
  328. {
  329. glCallList(1);
  332. return true;
  334. }