ossimDemGrid.h

Go to the documentation of this file.

//*******************************************************************
//
// License: MIT
//
// Author: Ken Melero
//         Orginally written by Jamie Moyers (jmoyers@geeks.com)
//         Adapted from the package KDEM.
// Description: This class defines a DEM grid.
//
//********************************************************************
// $Id$

#ifndef ossimDemGrid_HEADER
#define ossimDemGrid_HEADER 1

#include <ossim/base/ossimConstants.h>
#include <ossim/base/ossimIosFwd.h>
#include <ossim/support_data/ossimDemHeader.h>
#include <ossim/support_data/ossimDemProfile.h>
#include <vector>

class ossimDemGrid
{
public:

   ossimDemGrid(ossim_float32 missingDataValue = -32768.0);

   ~ossimDemGrid();

   long read(ossim::istream& dem, bool incrementalRead = false);

   // Accessors
   ossimDemHeader const& getHeader() const;

   long getWidth() const;

   long getHeight() const;

   ossim_float32 getElevation(long x, long y) const;

   ossim_float32 getMissingDataValue() const;

   void getGroundCoords(long x, long y, double& ground_x, double& ground_y);

private:

   ossim_float32   _missDataVal;
   ossimDemHeader  _header;
   long            _width;
   long            _height;
   ossim_float32*  _grid;
   bool            _firstTime; // Whether or not we've started reading the DEM.
   long            _curProfile;
   std::vector<ossimDemProfile> _profiles;  // Used by fillUTM()
   double _northwest_x, _northwest_y;

   void setElevation(long x, long y, ossim_float32 val);
   long fillGeographic(ossim::istream& dem, bool incrementalRead);
   long fillUTM(ossim::istream& dem, bool incrementalRead);
};

inline void ossimDemGrid::setElevation(long x, long y, ossim_float32 val)
{
   if (_grid == 0)
      return;

   _grid[(_width * y) + x] = val;
}

#endif  // #ifndef ossimDemGrid_HEADER