ossimLasReader.h

Go to the documentation of this file.

//----------------------------------------------------------------------------
//
// File: ossimLasReader.h
//
// License: MIT
// 
// See LICENSE.txt file in the top level directory for more details.
//
// Author:  David Burken
//
//----------------------------------------------------------------------------
// $Id$
#ifndef ossimLasReader_HEADER
#define ossimLasReader_HEADER 1

#include <ossim/imaging/ossimImageHandler.h>
#include <ossim/base/ossimDpt.h>
#include <ossim/base/ossimRefPtr.h>
#include <ossim/base/ossimUnitConversionTool.h>
#include <ossim/imaging/ossimImageData.h>
#include <ossim/imaging/ossimImageGeometry.h>
#include <ossim/projection/ossimProjection.h>
#include <mutex>

class ossimLasHdr;
class ossimLasPointRecordInterface;

class ossimLasReader : public ossimImageHandler
{
public:

   ossimLasReader();

   virtual ~ossimLasReader();

   virtual bool open();

   virtual bool isOpen()const;

   virtual void close();

   virtual ossimRefPtr<ossimImageData> getTile(const  ossimIrect& rect,
                                               ossim_uint32 resLevel=0);
   
   virtual bool getTile(ossimImageData* result, ossim_uint32 resLevel=0);
   
   virtual ossim_uint32 getNumberOfInputBands() const;

   virtual ossim_uint32 getNumberOfLines(ossim_uint32 resLevel = 0) const;

   virtual ossim_uint32 getNumberOfSamples(ossim_uint32 resLevel = 0) const;
    
   virtual ossim_uint32 getImageTileWidth() const;

   virtual ossim_uint32 getImageTileHeight() const;

   virtual ossim_uint32 getTileWidth() const;
   
   virtual ossim_uint32 getTileHeight() const;

   ossimScalarType getOutputScalarType() const;

   virtual void getEntryList(std::vector<ossim_uint32>& entryList) const;

   virtual ossim_uint32 getCurrentEntry() const;

   virtual bool setCurrentEntry(ossim_uint32 entryIdx);

   virtual ossimString getShortName() const;
   
   virtual ossimString getLongName()  const;

   virtual ossimRefPtr<ossimImageGeometry> getImageGeometry();

   virtual double getMinPixelValue(ossim_uint32 band=0) const;

   virtual double getMaxPixelValue(ossim_uint32 band=0) const;

   virtual double getNullPixelValue(ossim_uint32 band=0) const;

   virtual bool saveState(ossimKeywordlist& kwl, const char* prefix=0)const;
  
   virtual bool loadState(const ossimKeywordlist& kwl, const char* prefix=0);

   virtual void setProperty(ossimRefPtr<ossimProperty> property);

   virtual ossimRefPtr<ossimProperty> getProperty(const ossimString& name) const;

   virtual void getPropertyNames(std::vector<ossimString>& propertyNames) const;

   virtual ossim_uint32 getNumberOfDecimationLevels() const;

   virtual ossim_uint32 getNumberOfOutputBands() const;

protected:
   virtual void completeOpen();
   
private:

   class Bucket
   {
   public:
      Bucket(): a(0.0), c(0), red(0), green(0), blue(0) {}
         
      void add(const ossim_float64& point) { a += point; ++c; }
      ossim_float64 getValue() const { return ( c ? a/c : -99999.0 ); }
      void setRed(const ossim_uint16& value) { red = value; }
      void setGreen(const ossim_uint16& value) { green = value; }
      void setBlue(const ossim_uint16& value) { blue = value; }
      void setIntensity(const ossim_uint16& value) { intensity = value; }
      ossim_uint16 getRed() const { return red; }
      ossim_uint16 getGreen() const { return green; }
      ossim_uint16 getBlue() const { return blue; }
      ossim_uint16 getIntensity() const { return intensity; }

      ossim_float64 a; // accumulation
      ossim_uint32  c; // count
      ossim_uint16 red;
      ossim_uint16 green;
      ossim_uint16 blue;
      ossim_uint16 intensity;
   };

   bool init();
   
   void initValues(); // m_ul, m_lr, m_minZ
   bool initProjection();
   void initTile();
   void initUnits(const ossimKeywordlist& geomKwl);
   bool parseVarRecords();

   bool initFromExternalMetadata();
   
   void getScale(ossimDpt& scale, ossim_uint32 resLevel) const;

   void setGsd( const ossim_float64& gsd );

   void convertToMeters(ossim_float64& value) const;

   ossimLasPointRecordInterface* getNewPointRecord() const;

   std::ifstream                m_str;
   ossimLasHdr*                 m_hdr;
   ossimRefPtr<ossimProjection> m_proj;
   ossimDpt                     m_ul;
   ossimDpt                     m_lr;
   ossim_float64                m_maxZ;
   ossim_float64                m_minZ;
   ossimDpt                     m_gsd;
   ossimRefPtr<ossimImageData>  m_tile;
   ossim_uint8                  m_entry;
   std::mutex                   m_mutex;
   bool                         m_scan;  // Scan for bounds at open.
   ossimUnitType                m_units;
   ossimUnitConversionTool*     m_unitConverter;
TYPE_DATA
};

inline bool ossimLasReader::isOpen()const
{
   return ( m_str.is_open() && m_hdr );
}

inline void ossimLasReader::convertToMeters(ossim_float64& value) const
{
   if ( value )
   {
      m_unitConverter->setValue(value, m_units);
      value = m_unitConverter->getMeters();
   }
}

#endif /* #ifndef ossimLasReader_HEADER */