|| || || || || || || ||

---

HeightMap.h

Go to the documentation of this file.

#ifndef Impala_Visualisation_Plot_HeightMap_h
#define Impala_Visualisation_Plot_HeightMap_h

#include "Visualization/Plot/Plottable.h"
#include "Core/Array/PixMinMax.h"

namespace Impala
{
namespace Visualization
{
namespace Plot
{


class HeightMap : public Plottable
{
public:
    HeightMap()
    {
        mData = 0;
        mScanX = 0;
        mScanY = 0;
    }

    virtual void
    Draw(Plot* plot)
    {
        if (mData == 0)
            return;

        // adjust rendering settings
        glPushAttrib(GL_ALL_ATTRIB_BITS);
        glEnable(GL_LIGHTING);
        glEnable(GL_LIGHT0);
        glEnable(GL_NORMALIZE);
        // draw surface sligtly behind lines
        glEnable(GL_POLYGON_OFFSET_FILL);
        glPolygonOffset(1,1);

        // setup light independent of model view
        glMatrixMode (GL_MODELVIEW);
        glPushMatrix();
        glLoadIdentity();
        GLfloat materialSpecular[] = { 1.0, 1.0, 1.0, 1.0 };
        GLfloat materialShininess[] = { 10.0 };
        GLfloat lightPosition[] = { 0.0, 0.0, 1000.0, 0.0 };
        glClearColor (0.0, 0.0, 0.0, 0.0);
        glShadeModel (GL_SMOOTH);
        glMaterialfv(GL_FRONT, GL_SPECULAR, materialSpecular);
        glMaterialfv(GL_FRONT, GL_SHININESS, materialShininess);
        glLightfv(GL_LIGHT0, GL_POSITION, lightPosition);
        glPopMatrix();


        int h = mData->H()-1;
        int w = mData->W()-1;
        double x,z;
        double dx = ((double)mMaxX-mMinX) / (double)w;
        double dz = ((double)mMaxZ-mMinZ) / (double)h;
        double v1,v2,v3,v4;
        z = mMinZ;
        for (int i=0 ; i<h ; i++)
        {
            x = mMinX;
            for (int j=0 ; j<w ; j++)
            {
                v1 = mData->Val(j,i);
                v2 = mData->Val(j,i+1);
                v3 = mData->Val(j+1,i+1);
                v4 = mData->Val(j+1,i);

                // outlines
                glBegin(GL_LINE_LOOP);
                glColor3f(0, 0, 0);
                glMaterialfv(GL_FRONT, GL_AMBIENT, materialSpecular);
                glVertex3d(x,   v1,z);
                glVertex3d(x,   v2,z+dz);
                glVertex3d(x+dx,v3,z+dz);
                glVertex3d(x+dx,v4,z);
                glEnd();

                // surface
                glBegin(GL_QUADS);
                glNormal3d(-dz*(v4-v1), -dx*(v2-v1), dx*dz);
                Val2Material(v1); glVertex3d(x,   v1,z);
                Val2Material(v2); glVertex3d(x,   v2,z+dz);
                Val2Material(v3); glVertex3d(x+dx,v3,z+dz);
                Val2Material(v4); glVertex3d(x+dx,v4,z);
                glEnd();

                x += dx;
            }
            z += dz;
        }

        // prevent our render settings from affecting the rest
        glPopAttrib();
    }

    void
    SetData(Core::Array::Array2dScalarReal64* data)
    {
        mData = data;
        double min, max;
        Core::Array::PixMinMax(data, &min, &max);
        int w = data->W();
        int h = data->H();
        SetDimensionsX(-w/2, w/2);
        SetDimensionsY(min, max);
        SetDimensionsZ(-h/2, h/2);
    }

    void
    Val2Color(double v)
    {
        v -= mMinY;
        v /= (mMaxY - mMinY);
        double r,g,b;
        r = v;
        g = 1.-v;
        b = -fabs(v - .5)*2.;
        glColor3f(r, g, b);
    }

    void
    Val2Material(double v)
    {
        GLfloat material[4] = { 1.0, 1.0, 1.0, 1.0 };
        v -= mMinY;
        v /= (mMaxY - mMinY);
        material[0] = v;
        material[1] = 1.-v;
        material[2] = -fabs(v - .5)*2.;
        glMaterialfv(GL_FRONT, GL_AMBIENT, material);
    }

protected:
    Core::Array::Array2dScalarReal64* mData;
    int mScanX, mScanY;
};

}//namespace Impala
}//namespace Visulization
}//namespace Plot

#endif Impala_Visualisation_Plot_HeightMap_h

---