geo/tcElevationData.cpp

   1 /***************************************************************************
   2  *   This file is part of Tecorrec.                                        *
   3  *   Copyright 2008 James Hogan <james@albanarts.com>                      *
   4  *                                                                         *
   5  *   Tecorrec is free software: you can redistribute it and/or modify      *
   6  *   it under the terms of the GNU General Public License as published by  *
   7  *   the Free Software Foundation, either version 2 of the License, or     *
   8  *   (at your option) any later version.                                   *
   9  *                                                                         *
  10  *   Tecorrec is distributed in the hope that it will be useful,           *
  11  *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
  12  *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
  13  *   GNU General Public License for more details.                          *
  14  *                                                                         *
  15  *   You should have received a copy of the GNU General Public License     *
  16  *   along with Tecorrec.  If not, write to the Free Software Foundation,  *
  17  *   Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.   *
  18  ***************************************************************************/
  19
  20 /**
  21  * @file tcElevationData.cpp
  22  * @brief Terrain elevation data.
  23  */
  24
  25 #include "tcElevationData.h"
  26 #include "tcSrtmData.h"
  27 #include <glVector.h>
  28
  29 #include <GL/gl.h>
  30
  31 #include <cmath>
  32
  33 /*
  34  * Constructors + destructor
  35  */
  36
  37 /// Initialise from SRTM data.
  38 tcElevationData::tcElevationData(tcSrtmData* srtm)
  39 : tcGeoData(tcGeoData::ElevationData)
  40 , m_srtm(srtm)
  41 {
  42 }
  43
  44 /// Destructor.
  45 tcElevationData::~tcElevationData()
  46 {
  47 }
  48
  49 /*
  50  * Rendering
  51  */
  52
  53 void tcElevationData::renderSchematic(double meanRadius)
  54 {
  55   tcGeoData::renderSchematic(meanRadius);
  56
  57   for (int i = 0; i < m_srtm->lines(); i += 2)
  58   {
  59     glBegin(GL_LINE_STRIP);
  60     {
  61       for (int j = 0; j < m_srtm->length(); j += 2)
  62       {
  63         uint32_t alt = m_srtm->pixelAltitude(i,j);
  64         if (alt == 32768)
  65         {
  66           glEnd();
  67           glBegin(GL_LINE_STRIP);
  68           continue;
  69         }
  70         tcGeo coord = m_srtm->pixelCoordinate(i,j);
  71         GLvec3d dir = coord;
  72         double rad = meanRadius + alt;
  73         glColor3f(1.0f, 1.0f-(float)alt/3278.0f, 0.0f);
  74         glVertex3(dir * rad);
  75       }
  76     }
  77     glEnd();
  78   }
  79 }