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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 传真(Fax)编程 > 查看源码
ogl3d.cpp
资源名称：GTIFFSample.zip　[点击查看]
上传用户：lbr_007
上传日期：2019-05-31
资源大小：282k
文件大小：19k
源码类别：

传真(Fax)编程

开发平台：

Visual C++

ogl3d.cpp：源码内容
  1.  // ogl3d.cpp
  2. //
  4. #include "stdafx.h"
  5. #include "ogl3d.h"
  6. #include "dibbmp.h"
  8. #include <math.h>
  9. #include <vector>
  10. #include <arraycontainer.h>
  12. const int OGL3D_EARTH_LIST = 111;
  13. const int OGL3D_EARTH_LINES = 222;
  14. const int SOLAR_SYSTEM_LIST = 333;
  16. OGL3D::OGL3D(void)
  17. {
  18. InitOGL3D();
  19. }
  21. OGL3D::~OGL3D(void)
  22. {
  23. }
  25. void OGL3D::SetTexture(DIBSection& dib)
  26. {
  27. if (dib.IsCreated())
  28. {
  29. GenTexture(dib);
  30. }
  31. }
  33. void OGL3D::GenTexture(DIBSection& dib)
  34. {
  35. UINT32 width = dib.Width();
  36. UINT32 height = dib.Height();
  38. int w, h;
  40. UINT32 max_dimension = width > height ? width : height;
  42. if (width >= 512) w = 512;
  43. else if (width >= 256) w = 256;
  44. else if (width >= 128) w = 128;
  45. else if (width >= 64) w = 64;
  46. else if (width >= 32) w = 32;
  47. else w = 16;
  49. if (height >= 512) h = 512;
  50. else if (height >= 256) h = 256;
  51. else if (height >= 128) h = 128;
  52. else if (height >= 64) h = 64;
  53. else if (height >= 32) h = 32;
  54. else h = 16;
  56. dib.ResizeImage(m_texture_dib,w,h);
  58. if (m_texture_dib.IsCreated())
  59. {
  60. m_have_texture = 1;
  61. }
  62. }
  64. void OGL3D::RenderSpecialEffects(void)
  65. {
  66. double origin[] = { 0.0, 0.0, 0.0};
  67. double v1[] = { -0.25, 0.65, -0.95 };
  68. double v2[] = { 0.70, -0.60, 0.80 };
  69. double v3[3];
  71. v3[0] = (v1[1]*v2[2] - v1[2]*v2[1]);
  72. v3[1] = (v1[2]*v2[0] - v1[0]*v2[1]);
  73. v3[2] = (v1[0]*v2[1] - v1[1]*v2[0]);
  75. glBegin( GL_LINES );
  77. glColor3ub(64,64,64);
  79. glVertex3dv((GLdouble *)origin);
  80. glVertex3dv((GLdouble *)v1);
  82. glVertex3dv((GLdouble *)origin);
  83. glVertex3dv((GLdouble *)v2);
  85. glVertex3dv((GLdouble *)origin);
  86. glVertex3dv((GLdouble *)v3);
  88. glEnd();
  90. DrawCone(v1);
  91. DrawCone(v2);
  93. glColor3ub(192,192,192);
  94. DrawCone(v3);
  95. }
  97. void OGL3D::RenderPath(void)
  98. {
  99. GLdouble u[3], v[3], w[3];
  100. GLdouble R = 1.005;
  101. double lat1 = 39.0;
  102. double lon1 = 77.0;
  103. double lat2 = 0.0;
  104. double lon2 = 110.0;
  106. const double delta = 500.0;
  108. double lat, lon, phi1, phi2, theta1, theta2, lon_incr, lat_incr, d;
  110. glColor3ub(128,128,128);
  112. u[0] = 0.0; u[1] = 1.25; u[2] = 0.0;
  113. v[0] = 0.0; v[1] = -1.25; v[2] = 0.0;
  114. DrawCone(u);
  115. DrawCone(v);
  117. glBegin( GL_LINES );
  118. {
  121. phi1 = phi2 = lon1 * GEO::DE2RA;
  123. glVertex3dv(u);
  124. glVertex3dv(v);
  126. glColor3ub(32,32,32);
  128. lat_incr = (90.0 - lat1) / delta;
  130. for (lat = 90.0; lat >= (lat1+lat_incr); lat -= lat_incr){
  131. theta1 = lat * GEO::DE2RA;
  132. theta2 = (lat - lat_incr) * GEO::DE2RA;
  134. u[0] = R * cos(phi1) * cos(theta1);//x
  135. u[1] = R * sin(theta1);//y
  136. u[2] = R * sin(phi1) * cos(theta1);//z
  138. v[0] = R * cos(phi2) * cos(theta2);//x
  139. v[1] = R * sin(theta2);//y
  140. v[2] = R * sin(phi2) * cos(theta2);//z
  142. glVertex3dv(u);
  143. glVertex3dv(v);
  144. }
  146. v[0] = v[1] = v[2] = 0.0;
  147. glVertex3dv(u);
  148. glVertex3dv(v);
  150. phi1 = phi2 = lon2 * GEO::DE2RA;
  151. lat_incr = (90.0 - lat2) / delta;
  153. for (lat = 90.0; lat >= (lat2+lat_incr); lat -= lat_incr){
  154. theta1 = lat * GEO::DE2RA;
  155. theta2 = (lat - lat_incr) * GEO::DE2RA;
  157. u[0] = R * cos(phi1) * cos(theta1);//x
  158. u[1] = R * sin(theta1);//y
  159. u[2] = R * sin(phi1) * cos(theta1);//z
  161. v[0] = R * cos(phi2) * cos(theta2);//x
  162. v[1] = R * sin(theta2);//y
  163. v[2] = R * sin(phi2) * cos(theta2);//z
  165. glVertex3dv(u);
  166. glVertex3dv(v);
  167. }
  169. lat_incr = (lat2 - lat1) / delta;
  170. lon_incr = (lon2 - lon1) / delta;
  172. for (d = 0.0; d < (delta-1.0); d += 1.0000){
  174. theta1 = (lat1 + d * lat_incr) * GEO::DE2RA;
  175. theta2 = (lat1 + (d+1.00) * lat_incr) * GEO::DE2RA;
  176. phi1 = (lon1 + d * lon_incr) * GEO::DE2RA;
  177. phi2 = (lon1 + (d+1.00) * lon_incr) * GEO::DE2RA;
  179. u[0] = R * cos(phi1) * cos(theta1);//x
  180. u[1] = R * sin(theta1);//y
  181. u[2] = R * sin(phi1) * cos(theta1);//z
  183. v[0] = R * cos(phi2) * cos(theta2);//x
  184. v[1] = R * sin(theta2);//y
  185. v[2] = R * sin(phi2) * cos(theta2);//z
  187. glVertex3dv(u);
  188. glVertex3dv(v);
  189. }
  191. R = 2.0;
  192. u[0] = R * cos(phi2) * cos(theta2);//x
  193. u[1] = R * sin(theta2);//y
  194. u[2] = R * sin(phi2) * cos(theta2);//z
  195. glVertex3dv(u);
  196. glVertex3dv(v);
  198. R = 1.005;
  199. theta1 = lat1 * GEO::DE2RA;
  200. phi1 = lon1 * GEO::DE2RA;
  201. v[0] = 0.0;//R * cos(phi1) * cos(theta1);//x
  202. v[1] = 0.0;//R * sin(theta1);//y
  203. v[2] = 0.0;//R * sin(phi1) * cos(theta1);//z
  204. glVertex3dv(u);
  205. glVertex3dv(v);
  206. }
  207. glEnd();
  209. glColor3ub(200,200,255);
  211. glBegin( GL_TRIANGLES );
  212. {
  213. R = 1.005;
  215. phi1 = lon1 * GEO::DE2RA;
  216. theta1 = lat1 * GEO::DE2RA;
  217. phi2 = lon2 * GEO::DE2RA;
  218. theta2 = lat2 * GEO::DE2RA;
  220. u[0] = R * cos(phi1) * cos(theta1);//x
  221. u[1] = R * sin(theta1);//y
  222. u[2] = R * sin(phi1) * cos(theta1);//z
  224. R = 2.0;
  225. v[0] = R * cos(phi2) * cos(theta2);//x
  226. v[1] = R * sin(theta2);//y
  227. v[2] = R * sin(phi2) * cos(theta2);//z
  229. w[0] = w[1] = w[2] = 0.0;
  231. glVertex3dv(u);
  232. glVertex3dv(v);
  233. glVertex3dv(w);
  234. }
  235. glEnd();
  237. }
  239. void OGL3D::RenderPlane(void)
  240. {
  241. GLdouble u[3], v[3], w[3], n[3], uu[3], vv[3];
  242. double delta = 1.0;
  243. double theta;
  244. double R = 1.010;
  245. UINT32 count;
  247. w[0] = w[1] = w[2] = 0.0;
  249. glColor3ub(200,200,255);
  251. /*
  252. glPushMatrix();
  253. {
  254. glRotated(20.0, 0.0, 0.0, 1.0);
  256. glBegin( GL_TRIANGLES );
  257. {
  258. count = 0;
  259. for (theta=0.0;theta < 360.0; theta += delta){
  260. u[0] = R * cos(theta*GEO::DE2RA);
  261. u[1] = 0.0;
  262. u[2] = R * sin(theta*GEO::DE2RA);
  264. v[0] = R * cos((theta+delta)*GEO::DE2RA);
  265. v[1] = 0.0;
  266. v[2] = R * sin((theta+delta)*GEO::DE2RA);
  268. glVertex3dv(w);
  269. glVertex3dv(v);
  270. glVertex3dv(u);
  272. if (count == 0)
  273. {
  274. memcpy(uu,u,(size_t)(sizeof(double)*3));
  275. }
  276. if (count == 35)
  277. {
  278. memcpy(vv,u,(size_t)(sizeof(double)*3));
  279. }
  280. if (count == 12)
  281. {
  282. n[0] = 1.2 * cos(theta*GEO::DE2RA);
  283. n[1] = 0.0;
  284. n[2] = 1.2 * sin(theta*GEO::DE2RA);
  285. }
  286. count++;
  287. }
  288. }
  289. glEnd();
  291. glColor3ub(32,32,32);
  293. DrawCone(uu,false);
  294. DrawCone(vv,false);
  296. glBegin( GL_LINES );
  297. {
  298. u[0] = 0.0; u[1] = 0.75; u[2] = 0.0;
  299. glVertex3dv(w);
  300. glVertex3dv(u);
  301. //glVertex3dv(w);
  302. //glVertex3dv(n);
  303. }
  304. glEnd();
  306. DrawCone(u);
  307. //DrawCone(n);
  308. }
  309. glPopMatrix();
  310. */
  312. glColor3ub(200,200,255);
  314. glPushMatrix();
  315. {
  316. glRotated(60.0, 0.0, 1.0, 0.0);
  317. glRotated(50.0, 0.0, 0.0, 1.0);
  319. glBegin( GL_TRIANGLES );
  320. {
  321. count = 0;
  323. for (theta=0.0;theta < 360.0; theta += delta){
  324. u[0] = R * cos(theta*GEO::DE2RA);
  325. u[1] = 0.0;
  326. u[2] = R * sin(theta*GEO::DE2RA);
  327. v[0] = R * cos((theta+delta)*GEO::DE2RA);
  328. v[1] = 0.0;
  329. v[2] = R * sin((theta+delta)*GEO::DE2RA);
  331. glVertex3dv(w);
  332. glVertex3dv(v);
  333. glVertex3dv(u);
  335. if (count == 25)
  336. {
  337. memcpy(uu,u,(size_t)(sizeof(double)*3));
  338. }
  339. if (count == 75)
  340. {
  341. memcpy(vv,u,(size_t)(sizeof(double)*3));
  342. }
  343. count++;
  344. }
  345. }
  346. glEnd();
  348. glColor3ub(32,32,32);
  350. DrawCone(uu,false);
  351. DrawCone(vv,false);
  353. glBegin( GL_LINES );
  354. {
  355. u[0] = 0.0; u[1] = 0.75; u[2] = 0.0;
  356. glVertex3dv(w);
  357. glVertex3dv(u);
  358. }
  359. glEnd();
  361. DrawCone(u);
  362. }
  363. glPopMatrix();
  364. }
  366. void OGL3D::RenderPrimaryImage(void)
  367. {
  368. if ( m_draw_mode == DrawLines )
  369. {
  370. glPolygonMode(GL_FRONT_AND_BACK,GL_LINE);
  371. }
  372. else
  373. {
  374. glPolygonMode(GL_FRONT_AND_BACK,GL_FILL);
  375. }
  377. #ifndef RENDER_FUNCTION
  378. glColor3ub(96,96,96);
  379. RenderGlobe();
  380. //RenderCylinder();
  381. RenderRectangle();
  383. //RenderCone();
  385. //GLView3D::RenderPrimaryImage();
  386. //RenderSpecialEffects();
  387. //RenderEarthImage();
  388. //RenderPath();
  389. //RenderPlane();
  390. #else
  391. //RenderFunction();
  392. #endif
  393. }
  395. namespace FXYDomain {
  396. const double xmin = 0.0;
  397. const double xmax = 8.0;
  398. const double ymin = 0.0;
  399. const double ymax = 8.0;
  400. const double xincr = 0.50;
  401. const double yincr = 0.50;
  402. };
  404. #define FXY(x,y) (2.0 * exp(-x) * exp(-2.0 * y))
  406. void OGL3D::RenderFunction(void)
  407. {
  408. //GLView3D::RenderPrimaryImage();
  410. UINT32 ndx1, ndx2;
  411. double x, y, z, zmin, zmax;
  412. double scale;
  413. GLdouble u[3], v[3], w[3], n[3];
  415. UINT32 nrows = (UINT32)(0.50 + (FXYDomain::ymax - FXYDomain::ymin)/FXYDomain::yincr);
  416. UINT32 ncols = (UINT32)(0.50 + (FXYDomain::xmax - FXYDomain::xmin)/FXYDomain::xincr);
  418. zmin = 1e10;
  419. zmax = -1e10;
  421. ArrayContainer<double> grid;
  422. grid.SetBufferSize(nrows*ncols);
  424. for (UINT32 row = 0; row < nrows; row++){
  425. y = FXYDomain::ymin + row * FXYDomain::yincr;
  426. for (UINT32 col = 0; col < ncols; col++){
  427. x = FXYDomain::xmin + col * FXYDomain::xincr;
  428. z = FXY(x,y);
  429. if (z < zmin) zmin = z;
  430. if (z > zmax) zmax = z;
  431. grid.AddElement(z);
  432. }
  433. }
  435. double xscale = 1.0 / (FXYDomain::xmax - FXYDomain::xmin);
  436. double yscale = 1.0 / (FXYDomain::ymax - FXYDomain::ymin);
  437. double zscale = 1.0 / (zmax - zmin);
  439. scale = (xscale < yscale) ? xscale : yscale;
  440. scale = (scale < zscale) ? scale : zscale;
  442. double * ptr = grid.GetData();
  444. glColor3ub(0,0,224);
  446. if ( m_draw_mode == DrawLines )
  447. {
  448. glPolygonMode(GL_FRONT_AND_BACK,GL_LINE);
  449. }
  450. else
  451. {
  452. glPolygonMode(GL_FRONT_AND_BACK,GL_FILL);
  453. }
  455. /*
  456. glBegin( GL_TRIANGLES );
  457. {
  458. for (UINT32 row = 0; row < (nrows-1); row++){
  459. y = FXYDomain::ymin + row * FXYDomain::yincr;
  461. for (UINT32 col = 0; col < (ncols-1); col++){
  462. x = FXYDomain::xmin + col * FXYDomain::xincr;
  464. ndx1 = row * ncols + col;
  465. ndx2 = (row + 1) * ncols + col;
  467. u[0] = x*scale; u[1] = ptr[ndx1]*scale; u[2] = y*scale;
  468. v[0] = (x + FXYDomain::xincr)*scale; v[1] = ptr[ndx2+1]*scale; v[2] = (y + FXYDomain::yincr)*scale;
  469. w[0] = x*scale; w[1] = ptr[ndx2]*scale; w[2] = (y + FXYDomain::yincr)*scale;
  471. NormalVector(u,v,w,n);
  472. glNormal3dv(n);
  473. glVertex3dv(u);
  474. glVertex3dv(v);
  475. glVertex3dv(w);
  477. u[0] = x*scale; u[1] = ptr[ndx1]*scale; u[2] = y*scale;
  478. v[0] = (x + FXYDomain::xincr)*scale; v[1] = ptr[ndx1+1]*scale; v[2] = y*scale;
  479. w[0] = (x + FXYDomain::xincr)*scale; w[1] = ptr[ndx2+1]*scale; w[2] = (y + FXYDomain::yincr)*scale;
  481. NormalVector(u,w,v,n);
  482. glNormal3dv(n);
  483. glVertex3dv(u);
  484. glVertex3dv(w);
  485. glVertex3dv(v);
  486. }
  487. }
  488. }
  489. glEnd();
  490. */
  492. glColor3ub(224,0,0);
  494. glBegin(GL_LINES);
  495. {
  496. for (int i=1;i<19;i++){
  497. y = (double)i/2.0;
  498. u[0] = sqrt(y)*scale; u[1] = 0.0; u[2] = y*scale;
  499. v[0] = 8.0*scale;; v[1] = 0.0; v[2] = y*scale;
  501. glVertex3dv(u);
  502. glVertex3dv(v);
  503. }
  504. }
  505. glEnd();
  506. }
  508. void OGL3D::RenderRectangle()
  509. {
  510. glColor3ub(224,0,0);
  511. glPushMatrix();
  512. {
  513. GLdouble D = 0.50;
  514. glRotated(-45.0, 0.0, 1.0, 0.0);
  515. glRotated(-60.0, 1.0, 0.0, 0.0);
  516. glTranslated(0.0, 0.0, 1.0);
  518. glBegin(GL_LINES);
  519. {
  520. GLdouble u[3], v[3], w[3];
  522. u[0] = -D; u[1] = -D; u[2] = 0.0;
  523. v[0] = D; v[1] = -D; v[2] = 0.0;
  524. glVertex3dv(u);
  525. glVertex3dv(v);
  526. u[0] = D; u[1] = D; u[2] = 0.0;
  527. glVertex3dv(v);
  528. glVertex3dv(u);
  529. v[0] = -D; v[1] = D; v[2] = 0.0;
  530. glVertex3dv(u);
  531. glVertex3dv(v);
  532. u[0] = -D; u[1] = -D; u[2] = 0.0;
  533. glVertex3dv(v);
  534. glVertex3dv(u);
  535. }
  536. glEnd();
  537. }
  538. glPopMatrix();
  539. }
  541. void OGL3D::RenderCylinder()
  542. {
  543. glPushMatrix();
  544. {
  546. glRotated(90.0, 0.0, 1.0, 0.0);
  547. glTranslated(0.0, 0.0, -1.0);
  549. GLUquadricObj * q = gluNewQuadric();
  550. gluQuadricDrawStyle(q, GLU_LINE);
  551. gluQuadricNormals(q, GLU_SMOOTH);
  552. gluCylinder(q, 1.0, 1.0, 2.0, 36, 36);
  553. gluDeleteQuadric(q);
  554. }
  555. glPopMatrix();
  556. }
  558. void OGL3D::RenderCone()
  559. {
  560. glPushMatrix();
  561. {
  562. glTranslated(0.0, 1.414, 0.0);
  563. glRotated(90.0, 1.0, 0.0, 0.0);
  565. GLUquadricObj * q = gluNewQuadric();
  566. //gluQuadricDrawStyle(q, m_draw_mode==DrawLines ? GLU_LINE : GLU_FILL );
  567. gluQuadricDrawStyle(q, GLU_LINE);
  568. gluQuadricNormals(q, GLU_SMOOTH);
  569. gluCylinder(q, 0.00, 1.414, 1.414, 36, 36);
  570. gluDeleteQuadric(q);
  571. }
  572. glPopMatrix();
  573. }
  575. void OGL3D::RenderGlobe()
  576. {
  577. /*
  578. glPushMatrix();
  579. {
  580. glRotated(-90.0, 1.0, 0.0, 0.0);
  582. GLUquadricObj * q = gluNewQuadric();
  583. //gluQuadricDrawStyle(q, m_draw_mode==DrawLines ? GLU_LINE : GLU_FILL);
  584. gluQuadricDrawStyle(q, GLU_FILL);
  585. gluQuadricNormals(q, GLU_SMOOTH);
  586. gluSphere(q,1.0,16,16);
  587. gluDeleteQuadric(q);
  588. }
  589. glPopMatrix();
  590. */
  592. const int NumLatitudes = 18;
  593. const int NumLongitudes = 36;
  595. double start_lat, start_lon;
  596. double theta1, phi1, theta2, phi2, lat_incr, lon_incr;
  598. GLdouble u[3], v[3], w[3], n[3];
  599. GLdouble R;
  601. if (!m_model_built)
  602. {
  603. glNewList( OGL3D_EARTH_LIST, GL_COMPILE);
  605. glColor3ub(192,192,192);
  607. if (m_have_texture)
  608. {
  609. glTexImage2D(GL_TEXTURE_2D,0,3,m_texture_dib.Width(),m_texture_dib.Height(),
  610. 0,GL_BGR_EXT,GL_UNSIGNED_BYTE,(GLvoid *)m_texture_dib.GetBits());
  612. glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_REPEAT);
  613. glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_REPEAT);
  614. glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
  615. glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
  616. glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_DECAL);
  618. //glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_BLEND);
  619. //glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_MODULATE);
  621. glEnable( GL_TEXTURE_2D );
  622. }
  623. else
  624. {
  625. glDisable( GL_TEXTURE_2D );
  626. }
  628. m_model_built = 1;
  630. glBegin( GL_TRIANGLES );
  632. start_lat = -90;
  633. start_lon = 0.0;
  634. R = 1.0;
  636. lat_incr = 180.0 / NumLatitudes;
  637. lon_incr = 360.0 / NumLongitudes;
  639. int row, col;
  641. for (col = 0; col < NumLongitudes; col++){
  642. phi1 = (start_lon + col * lon_incr) * GEO::DE2RA;
  643. phi2 = (start_lon + (col + 1) * lon_incr) * GEO::DE2RA;
  645. for (row = 0; row < NumLatitudes; row++){
  646. theta1 = (start_lat + row * lat_incr) * GEO::DE2RA;
  647. theta2 = (start_lat + (row + 1) * lat_incr) * GEO::DE2RA;
  649. u[0] = R * cos(phi1) * cos(theta1);//x
  650. u[1] = R * sin(theta1);//y
  651. u[2] = R * sin(phi1) * cos(theta1);//z
  653. v[0] = R * cos(phi1) * cos(theta2);//x
  654. v[1] = R * sin(theta2);//y
  655. v[2] = R * sin(phi1) * cos(theta2);//z
  657. w[0] = R * cos(phi2) * cos(theta2);//x
  658. w[1] = R * sin(theta2);//y
  659. w[2] = R * sin(phi2) * cos(theta2);//z
  661. NormalVector(u,v,w,n);
  662. glNormal3dv(n);
  663. glTexCoord2d((180.0 - phi1*GEO::RA2DE)/360.0,(theta1 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  664. glVertex3dv(u);
  665. glTexCoord2d((180.0 - phi1*GEO::RA2DE)/360.0,(theta2 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  666. glVertex3dv(v);
  667. glTexCoord2d((180.0 - phi2*GEO::RA2DE)/360.0,(theta2 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  668. glVertex3dv(w);
  670. v[0] = R * cos(phi2) * cos(theta1);//x
  671. v[1] = R * sin(theta1);//y
  672. v[2] = R * sin(phi2) * cos(theta1);//z
  674. NormalVector(u,w,v,n);
  675. glNormal3dv(n);
  676. glTexCoord2d((180.0 - phi1*GEO::RA2DE)/360.0,(theta1 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  677. glVertex3dv(u);
  678. glTexCoord2d((180.0 - phi2*GEO::RA2DE)/360.0,(theta2 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  679. glVertex3dv(w);
  680. glTexCoord2d((180.0 - phi2*GEO::RA2DE)/360.0,(theta1 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  681. glVertex3dv(v);
  682. }
  683. }
  685. glEnd();
  687. if (m_have_texture)
  688. {
  689. glDisable( GL_TEXTURE_2D );
  690. }
  693. glEndList();
  694. }
  696. glCallList( OGL3D_EARTH_LIST );
  697. }
  699. void OGL3D::RenderEarthImage(void)
  700. {
  701. const int NumLatitudes = 36;
  702. const int NumLongitudes = 72;
  704. double start_lat, start_lon;
  705. double theta1, phi1, theta2, phi2, lat_incr, lon_incr;
  707. GLdouble u[3], v[3], w[3], n[3];
  708. GLdouble R;
  710. if (!m_model_built)
  711. {
  712. glNewList( OGL3D_EARTH_LIST, GL_COMPILE);
  714. glColor3ub(192,192,192);
  716. if (m_have_texture)
  717. {
  718. glTexImage2D(GL_TEXTURE_2D,0,3,m_texture_dib.Width(),m_texture_dib.Height(),
  719. 0,GL_BGR_EXT,GL_UNSIGNED_BYTE,(GLvoid *)m_texture_dib.GetBits());
  721. glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_REPEAT);
  722. glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_REPEAT);
  723. glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
  724. glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
  725. glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_DECAL);
  727. //glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_BLEND);
  728. //glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_MODULATE);
  730. glEnable( GL_TEXTURE_2D );
  731. }
  732. else
  733. {
  734. glDisable( GL_TEXTURE_2D );
  735. }
  737. m_model_built = 1;
  739. glBegin( GL_TRIANGLES );
  741. start_lat = -90;
  742. start_lon = 0.0;
  743. R = 1.0;
  745. lat_incr = 180.0 / NumLatitudes;
  746. lon_incr = 360.0 / NumLongitudes;
  748. int row, col;
  750. for (col = 0; col < NumLongitudes; col++){
  751. phi1 = (start_lon + col * lon_incr) * GEO::DE2RA;
  752. phi2 = (start_lon + (col + 1) * lon_incr) * GEO::DE2RA;
  754. for (row = 0; row < NumLatitudes; row++){
  755. theta1 = (start_lat + row * lat_incr) * GEO::DE2RA;
  756. theta2 = (start_lat + (row + 1) * lat_incr) * GEO::DE2RA;
  758. u[0] = R * cos(phi1) * cos(theta1);//x
  759. u[1] = R * sin(theta1);//y
  760. u[2] = R * sin(phi1) * cos(theta1);//z
  762. v[0] = R * cos(phi1) * cos(theta2);//x
  763. v[1] = R * sin(theta2);//y
  764. v[2] = R * sin(phi1) * cos(theta2);//z
  766. w[0] = R * cos(phi2) * cos(theta2);//x
  767. w[1] = R * sin(theta2);//y
  768. w[2] = R * sin(phi2) * cos(theta2);//z
  770. NormalVector(u,v,w,n);
  771. glNormal3dv(n);
  772. glTexCoord2d((180.0 - phi1*GEO::RA2DE)/360.0,(theta1 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  773. glVertex3dv(u);
  774. glTexCoord2d((180.0 - phi1*GEO::RA2DE)/360.0,(theta2 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  775. glVertex3dv(v);
  776. glTexCoord2d((180.0 - phi2*GEO::RA2DE)/360.0,(theta2 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  777. glVertex3dv(w);
  779. v[0] = R * cos(phi2) * cos(theta1);//x
  780. v[1] = R * sin(theta1);//y
  781. v[2] = R * sin(phi2) * cos(theta1);//z
  783. NormalVector(u,w,v,n);
  784. glNormal3dv(n);
  785. glTexCoord2d((180.0 - phi1*GEO::RA2DE)/360.0,(theta1 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  786. glVertex3dv(u);
  787. glTexCoord2d((180.0 - phi2*GEO::RA2DE)/360.0,(theta2 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  788. glVertex3dv(w);
  789. glTexCoord2d((180.0 - phi2*GEO::RA2DE)/360.0,(theta1 + GEO::PIOVER2)*GEO::RA2DE/180.0);
  790. glVertex3dv(v);
  791. }
  792. }
  794. glEnd();
  796. if (m_have_texture)
  797. {
  798. glDisable( GL_TEXTURE_2D );
  799. }
  801. //GLView3D::RenderPrimaryImage();
  803. glEndList();
  805. /*
  806. glNewList( OGL3D_EARTH_LINES, GL_COMPILE );
  808. start_lat = -85.0;
  809. lat_incr = 5.0000;
  810. lon_incr = 5.0 / 60.0; // 5 minutes of arc
  811. R = 1.00100;
  813. glBegin( GL_LINES );
  815. for (col = 0; col < NumLongitudes; col++){
  816. phi1 = (start_lon + col * lon_incr) * GEO::DE2RA;
  817. phi2 = (start_lon + (col + 1) * lon_incr) * GEO::DE2RA;
  819. for (row = 0; row < NumLatitudes; row++){
  820. theta1 = (start_lat + row * lat_incr) * GEO::DE2RA;
  821. theta2 = (start_lat + (row + 1) * lat_incr) * GEO::DE2RA;
  823. u[0] = R * cos(phi1) * cos(theta1);//x
  824. u[1] = R * sin(theta1);//y
  825. u[2] = R * sin(phi1) * cos(theta1);//z
  827. v[0] = R * cos(phi1) * cos(theta2);//x
  828. v[1] = R * sin(theta2);//y
  829. v[2] = R * sin(phi1) * cos(theta2);//z
  831. w[0] = R * cos(phi2) * cos(theta2);//x
  832. w[1] = R * sin(theta2);//y
  833. w[2] = R * sin(phi2) * cos(theta2);//z
  835. glVertex3dv(u);
  836. glVertex3dv(v);
  838. glVertex3dv(v);
  839. glVertex3dv(w);
  841. v[0] = R * cos(phi2) * cos(theta1);//x
  842. v[1] = R * sin(theta1);//y
  843. v[2] = R * sin(phi2) * cos(theta1);//z
  845. glVertex3dv(u);
  846. glVertex3dv(v);
  848. glVertex3dv(v);
  849. glVertex3dv(w);
  850. }
  851. }
  853. glEnd();
  855. glEndList();
  856. */
  857. }
  859. glCallList( OGL3D_EARTH_LIST );
  860. //glCallList( OGL3D_EARTH_LINES );
  861. }