CreateIsolinesCore.h

Go to the documentation of this file.

/*  Copyright (C) 2001-2009 National Institute For Space Research (INPE) - Brazil.
 
    This file is part of the TerraLib - a Framework for building GIS enabled applications.
 
    TerraLib is free software: you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as published by
    the Free Software Foundation, either version 3 of the License,
    or (at your option) any later version.
 
    TerraLib is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
    GNU Lesser General Public License for more details.
 
    You should have received a copy of the GNU Lesser General Public License
    along with TerraLib. See COPYING. If not, write to
    TerraLib Team at <terralib-team@terralib.org>.
 */
 
/*!
  \file terralib/mnt/core/CreateIsolines.h
 
  \brief This file contains a class that represents the create isolines.
 
  \reference Adapted from TerraLib4.
*/
 
#ifndef __TERRALIB_MNT_INTERNAL_CREATEISOLINES_H
#define __TERRALIB_MNT_INTERNAL_CREATEISOLINES_H
 
// Terralib Includes
#include "Config.h"
#include "CreateIsolinesParams.h"
 
#include "../../../../src/terralib/dataaccess.h"
#include "../../../../src/terralib/dataaccess/datasource/DataSourceFactory.h"
#include "../../../../src/terralib/memory/DataSet.h"
#include "../../../../src/terralib/common/Holder.h"
#include "../../../../src/terralib/statistics/core/Enums.h"
 
// STL Includes
#include <map>
#include <memory>
 
 
namespace te
{
  //forward declarations
  namespace da  { class DataSetType; }
  namespace gm  { class Geometry; }
  namespace mem { class DataSet; }
}
 
struct  RasterBlockSize
{
  int m_idThread;
  int m_numRows;
  int m_initalRow;
  int m_finalRow;
  RasterBlockSize() : m_idThread(0),
   m_numRows(0),
   m_initalRow(0),
   m_finalRow(0) {}
};
 
class ConnectLinesParams
{
  public:
 
  /*! \brief Default constructor. */
    ConnectLinesParams();
 
  /*! \brief Virtual destructor. */
    ~ConnectLinesParams();
 
  public:
 
    double m_quota;
    int m_srid;
    std::vector< te::gm::LineString* >  m_vecSegments;
    std::vector< te::gm::LineString* > m_lsOut;
};
 
class GenerateSegmentsParams
{
  public:
 
  /*! \brief Default constructor. */
  GenerateSegmentsParams();
 
  /*! \brief Virtual destructor. */
  ~GenerateSegmentsParams();
 
  public:
 
  std::unique_ptr<te::rst::Raster> m_rasterPtr;
  std::vector< std::vector< te::gm::LineString* > > m_vecSegments;
  std::vector<double> m_nvals;
};
 
namespace te
{
  namespace mnt
  {
    class TEMNTEXPORT CreateIsolines
    {
      /*!
      \class CreateIsolines
 
      \brief Class used to execute the create isolines
 
      */
 
    public:
 
      /*! \brief Default constructor. */
      CreateIsolines();
 
      /*! \brief Virtual destructor. */
      ~CreateIsolines();
 
    public:
 
      void setInput(te::da::DataSourcePtr inRasterDsrc, 
        std::string inRasterName, 
        std::unique_ptr<te::da::DataSetType> inDsetType);
 
      std::unique_ptr<te::rst::Raster> getPrepareRaster();
 
      void setParams(std::vector<double> &nval, std::vector<double> &gval, double vmax, double vmin, double dummy, bool hasDummy);
  
      void setOutput(te::da::DataSourcePtr outDsrc, std::string dsname);
 
      bool run(std::unique_ptr<te::rst::Raster> raster);
      
      static bool generateSegments(std::unique_ptr<te::rst::Raster> raster, std::vector<double> nvals, std::vector< std::vector< te::gm::LineString* > >& vecSegments);
 
      static bool connectLines(std::vector< te::gm::LineString* >  vec, int srid, std::vector< te::gm::LineString* >& lsOut);
 
      static void interpolacao(int direction, te::gm::LineString* line, double quota, double coord, double c_inf, double c_sup, double z_inf, double z_sup);
 
      static bool connectLinesThreaded(ConnectLinesParams* params);
 
      static void firstCaseSpecial(double quota, double& lineSupLeft, double& lineSupRigth, double& lineInfLeft, double& lineInfRigth);
 
      static void segments(int idQuota,   std::unique_ptr<te::gm::LineString>  line, double quota, double ylg_sup, double xlg_ant, double xlg_pos, double ylg_inf,
      double& lineSupLeft, double& lineSupRigth, double& lineInfLeft, double& lineInfRigth, std::vector< std::vector<te::gm::LineString*> >& vecSegments);
 
      static bool generateSegmentsThreaded(GenerateSegmentsParams* params);
 
      std::vector<RasterBlockSize> calculateBlocks(unsigned int numRows, unsigned int numThreads);
 
      void rstMemoryBlock(std::unique_ptr<te::rst::Raster> raster, std::vector<RasterBlockSize> vecBlocks, std::vector<GenerateSegmentsParams*>& vecGenerateParams);
 
      bool connectLines1(std::vector<te::gm::LineString*>  vecSegments, double scale, std::vector<te::gm::LineString*>& lsOut);
      bool connectLines(std::vector<te::gm::LineString*>  vecSegments, std::vector<te::gm::LineString*>& lsOut);
     
    protected:
 
       te::da::DataSourcePtr m_inRasterDsrc;
      std::string m_inRasterName;
      std::unique_ptr<te::da::DataSetType> m_inRasterDsType;
 
      te::da::DataSourcePtr m_outDsrc;
      std::string m_outDsetName;
   
      std::vector<double> m_values;
      std::vector<double> m_guidevalues;
      int m_srid;
     
      static double m_vmax;
      static double m_vmin;
      static double m_dummy;
      static bool m_hasDummy;
    };
  }
}
 
#endif //__TERRALIB_MNT_INTERNAL_CREATEISOLINES_H