Add in the terrain from osghanggliding to make it pretty. For the ter…

…rain to appear properly textured make sure to call

export OSG_FILE_PATH=/path/to/osg-data

compile

@@ -1 +1 @@
-g++ visualization.cpp -o visualization  -L/usr/local/lib -losg -lOpenThreads -losgDB -losgGA -losgUtil -losgViewer -losgEarth -losgEarthUtil
+g++ visualization.cpp hat.cpp terrain.cpp trees.cpp sky.cpp -o visualization  -L/usr/local/lib -losg -lOpenThreads -losgDB -losgGA -losgUtil -losgViewer -losgEarth -losgEarthUtil

hat.cpp

@@ -0,0 +1,168 @@
+/* OpenSceneGraph example, osghangglide.
+*
+*  Permission is hereby granted, free of charge, to any person obtaining a copy
+*  of this software and associated documentation files (the "Software"), to deal
+*  in the Software without restriction, including without limitation the rights
+*  to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+*  copies of the Software, and to permit persons to whom the Software is
+*  furnished to do so, subject to the following conditions:
+*
+*  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+*  IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+*  FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+*  AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+*  LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+*  OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+*  THE SOFTWARE.
+*/
+
+#ifdef _MSC_VER
+#include <Windows.h>
+#pragma warning( disable : 4244 )
+#endif
+
+#include <osg/GL>
+#include <osg/Math>
+#include <stdio.h>
+
+#include "terrain_coords.h"
+#include "hat.h"
+
+static int inited = 0;
+
+static float dbcenter[3];
+static float dbradius;
+
+static void getDatabaseCenterRadius( float dbcenter[3], float *dbradius )
+{
+    int i;
+    double n=0.0;
+    double center[3] = { 0.0f, 0.0f, 0.0f };
+    float  cnt;
+
+    cnt = 39 * 38;
+    for( i = 0; i < cnt; i++ )
+    {
+        center[0] += (double)vertex[i][0];
+        center[1] += (double)vertex[i][1];
+        center[2] += (double)vertex[i][2];
+
+        n = n + 1.0;
+    }
+
+    center[0] /= n;
+    center[1] /= n;
+    center[2] /= n;
+
+    float r = 0.0;
+
+    //    for( i = 0; i < sizeof( vertex ) / (sizeof( float[3] )); i++ )
+    for( i = 0; i < cnt; i++ )
+    {
+        double d = sqrt(
+            (((double)vertex[i][0] - center[0]) * ((double)vertex[i][0] - center[0])) +
+            (((double)vertex[i][1] - center[1]) * ((double)vertex[i][1] - center[1])) +
+            (((double)vertex[i][2] - center[2]) * ((double)vertex[i][2] - center[2]))  );
+
+        if( d > (double)r ) r = (float)d;
+
+    }
+
+    *dbradius = r;
+    dbcenter[0] = (float)center[0];
+    dbcenter[1] = (float)center[1];
+    dbcenter[2] = (float)center[2];
+
+    int index = 19 * 39 + 19;
+    dbcenter[0] = vertex[index][0] - 0.15;
+    dbcenter[1] = vertex[index][1];
+    dbcenter[2] = vertex[index][2] + 0.35;
+}
+
+
+static void init( void )
+{
+    getDatabaseCenterRadius( dbcenter, &dbradius );
+    inited = 1;
+}
+
+
+static void getNormal( float *v1, float *v2, float *v3, float *n )
+{
+    float V1[4], V2[4];
+    float   f;
+    int i;
+
+    /* Two vectors v2->v1 and v2->v3 */
+
+    for( i = 0; i < 3; i++ )
+    {
+        V1[i] = v1[i] - v2[i];
+        V2[i] = v3[i] - v2[i];
+    }
+
+    /* Cross product between V1 and V2 */
+
+    n[0] = (V1[1] * V2[2]) - (V1[2] * V2[1]);
+    n[1] = -((V1[0] * V2[2]) - ( V1[2] * V2[0] ));
+    n[2] = (V1[0] * V2[1] ) - (V1[1] * V2[0] );
+
+    /* Normalize */
+
+    f = sqrtf( ( n[0] * n[0] ) + ( n[1] * n[1] ) + ( n[2] * n[2] ) );
+    n[0] /= f;
+    n[1] /= f;
+    n[2] /= f;
+}
+
+
+float Hat( float x, float y, float z )
+{
+    int m, n;
+    int i, j;
+    float tri[3][3];
+    float norm[3];
+    float d, pz;
+
+    if( inited == 0 ) init();
+
+    // m = columns
+    // n = rows
+    m = (sizeof( vertex ) /(sizeof( float[3])))/39;
+    n = 39;
+
+    i = 0;
+    while( i < ((m-1)*39) && x > (vertex[i+n][0] - dbcenter[0]) )
+        i += n;
+
+    j = 0;
+
+    while( j < n-1 && y > (vertex[i+j+1][1] - dbcenter[1]) )
+        j++;
+
+    tri[0][0] = vertex[i+0+j+0][0] - dbcenter[0];
+    tri[0][1] = vertex[i+0+j+0][1] - dbcenter[1];
+    //tri[0][2] = vertex[i+0+j+0][2] - dbcenter[2];
+    tri[0][2] = vertex[i+0+j+0][2];
+
+    tri[1][0] = vertex[i+n+j+0][0] - dbcenter[0];
+    tri[1][1] = vertex[i+n+j+0][1] - dbcenter[1];
+    //tri[1][2] = vertex[i+n+j+0][2] - dbcenter[2];
+    tri[1][2] = vertex[i+n+j+0][2];
+
+    tri[2][0] = vertex[i+0+j+1][0] - dbcenter[0];
+    tri[2][1] = vertex[i+0+j+1][1] - dbcenter[1];
+    //tri[2][2] = vertex[i+0+j+1][2] - dbcenter[2];
+    tri[2][2] = vertex[i+0+j+1][2];
+
+    getNormal( tri[0], tri[1], tri[2], norm );
+
+    d = (tri[0][0] * norm[0]) +
+        (tri[0][1] * norm[1]) +
+        (tri[0][2] * norm[2]);
+
+    d *= -1;
+    pz = (-(norm[0] * x) - (norm[1] * y) - d)/norm[2];
+
+    return z - pz;
+}

hat.h

@@ -0,0 +1,24 @@
+/* -*-c++-*- 
+*
+*  OpenSceneGraph example, osghangglide.
+*
+*  Permission is hereby granted, free of charge, to any person obtaining a copy
+*  of this software and associated documentation files (the "Software"), to deal
+*  in the Software without restriction, including without limitation the rights
+*  to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+*  copies of the Software, and to permit persons to whom the Software is
+*  furnished to do so, subject to the following conditions:
+*
+*  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+*  IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+*  FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+*  AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+*  LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+*  OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+*  THE SOFTWARE.
+*/
+
+#ifndef __HAT_H
+#define __HAT_H
+extern float Hat( float x, float y, float z );
+#endif

sky.cpp

@@ -0,0 +1,142 @@
+/* OpenSceneGraph example, osghangglide.
+*
+*  Permission is hereby granted, free of charge, to any person obtaining a copy
+*  of this software and associated documentation files (the "Software"), to deal
+*  in the Software without restriction, including without limitation the rights
+*  to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+*  copies of the Software, and to permit persons to whom the Software is
+*  furnished to do so, subject to the following conditions:
+*
+*  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+*  IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+*  FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+*  AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+*  LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+*  OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+*  THE SOFTWARE.
+*/
+
+#include <math.h>
+
+#include <osg/Geode>
+#include <osg/Geometry>
+#include <osg/Texture2D>
+#include <osg/TexEnv>
+#include <osg/Depth>
+#include <osg/StateSet>
+
+#include <osgDB/ReadFile>
+
+#ifdef _MSC_VER
+#pragma warning( disable : 4244 )
+#pragma warning( disable : 4305 )
+#endif
+
+using namespace osg;
+
+Node *makeSky( void )
+{
+    int i, j;
+    float lev[] = { -5, -1.0, 1.0, 15.0, 30.0, 60.0, 90.0  };
+    float cc[][4] =
+    {
+        { 0.0, 0.0, 0.15 },
+        { 0.0, 0.0, 0.15 },
+        { 0.4, 0.4, 0.7 },
+        { 0.2, 0.2, 0.6 },
+        { 0.1, 0.1, 0.6 },
+        { 0.1, 0.1, 0.6 },
+        { 0.1, 0.1, 0.6 },
+    };
+    float x, y, z;
+    float alpha, theta;
+    float radius = 20.0f;
+    int nlev = sizeof( lev )/sizeof(float);
+
+    Geometry *geom = new Geometry;
+
+    Vec3Array& coords = *(new Vec3Array(19*nlev));
+    Vec4Array& colors = *(new Vec4Array(19*nlev));
+    Vec2Array& tcoords = *(new Vec2Array(19*nlev));
+
+
+    int ci = 0;
+
+    for( i = 0; i < nlev; i++ )
+    {
+        for( j = 0; j <= 18; j++ )
+        {
+            alpha = osg::DegreesToRadians(lev[i]);
+            theta = osg::DegreesToRadians((float)(j*20));
+
+            x = radius * cosf( alpha ) * cosf( theta );
+            y = radius * cosf( alpha ) * -sinf( theta );
+            z = radius * sinf( alpha );
+
+            coords[ci][0] = x;
+            coords[ci][1] = y;
+            coords[ci][2] = z;
+
+            colors[ci][0] = cc[i][0];
+            colors[ci][1] = cc[i][1];
+            colors[ci][2] = cc[i][2];
+            colors[ci][3] = 1.0;
+
+            tcoords[ci][0] = (float)j/18.0;
+            tcoords[ci][1] = (float)i/(float)(nlev-1);
+
+            ci++;
+        }
+
+
+    }
+
+    for( i = 0; i < nlev-1; i++ )
+    {
+        DrawElementsUShort* drawElements = new DrawElementsUShort(PrimitiveSet::TRIANGLE_STRIP);
+        drawElements->reserve(38);
+
+        for( j = 0; j <= 18; j++ )
+        {
+            drawElements->push_back((i+1)*19+j);
+            drawElements->push_back((i+0)*19+j);
+        }
+
+        geom->addPrimitiveSet(drawElements);
+    }
+
+    geom->setVertexArray( &coords );
+    geom->setTexCoordArray( 0, &tcoords );
+
+    geom->setColorArray( &colors );
+    geom->setColorBinding( Geometry::BIND_PER_VERTEX );
+
+
+    Texture2D *tex = new Texture2D;
+    tex->setImage(osgDB::readImageFile("Images/white.rgb"));
+
+    StateSet *dstate = new StateSet;
+
+    dstate->setTextureAttributeAndModes(0, tex, StateAttribute::OFF );
+    dstate->setTextureAttribute(0, new TexEnv );
+    dstate->setMode( GL_LIGHTING, StateAttribute::OFF );
+    dstate->setMode( GL_CULL_FACE, StateAttribute::ON );
+
+
+    // clear the depth to the far plane.
+    osg::Depth* depth = new osg::Depth;
+    depth->setFunction(osg::Depth::ALWAYS);
+    depth->setRange(1.0,1.0);   
+    dstate->setAttributeAndModes(depth,StateAttribute::ON );
+
+    dstate->setRenderBinDetails(-2,"RenderBin");
+
+    geom->setStateSet( dstate );
+
+    Geode *geode = new Geode;
+    geode->addDrawable( geom );
+
+    geode->setName( "Sky" );
+
+    return geode;
+}