TINCreateIsolines.cpp

Go to the documentation of this file.

/*!
\file terralib/mnt/core/TINCreateIsolines.cpp

\brief This file contains a class to generate isolines from TIN.

*/

#include "TINCreateIsolines.h"
#include "Utils.h"

#include "../../common/progress/TaskProgress.h"


bool te::mnt::TINCreateIsolines::run()
{
  if (!m_values.size())
    return false;

  int32_t nodesid[3];
  te::gm::Point p3da[3] = {te::gm::Point(0, te::gm::PointZType),
                           te::gm::Point(0, te::gm::PointZType),
                           te::gm::Point(0, te::gm::PointZType)};
  std::vector<te::gm::Point> ptline;
  te::gm::LineString isoline(te::gm::LineStringZType, m_srid);
  std::vector<te::gm::Point> pts;

  double cvalue;

  std::sort(m_values.begin(), m_values.end());

  double zmin = m_values[0];
  double zmax = m_values[m_values.size() - 1];

  LoadTin(m_inDsrc, m_inDsetName, zmin, zmax);

  std::size_t ntri = m_triang.size();
  te::common::TaskProgress task("Creating Isolines...");
  task.setTotalSteps((int)(m_values.size()*ntri));

  for (size_t v = 0; v < m_values.size(); v++)
  {
    cvalue = m_values[v];

    for (size_t i = 0; i < ntri; i++)
    {
      if (!task.isActive())
        return false;
      task.pulse();
      if (!NodesId((int32_t)i, nodesid))
        continue;
      for (size_t j = 0; j < 3; j++)
      {
        p3da[j].setX(m_node[(unsigned int)nodesid[j]].getNPoint().getX());
        p3da[j].setY(m_node[(unsigned int)nodesid[j]].getNPoint().getY());
        p3da[j].setZ(m_node[(unsigned int)nodesid[j]].getZ());
      }
      if (!testVertexValues(cvalue, p3da))
        continue;

      // Define intersection between contour and triangle
      defineInterTriangle(pts, p3da, cvalue);
    }

    if (pts.size() > 2)
    {
      if (!extractLines(pts, m_isolist, m_tol))
      {
        return false;
      }
       
      pts.clear();
    }
  }

  SaveIso(m_outDsetName, m_outDsrc, m_isolist, m_guidevalues, m_srid);

  return true;
}

void te::mnt::TINCreateIsolines::setInput(te::da::DataSourcePtr inDsrc,
  std::string inDsetName,
  std::unique_ptr<te::da::DataSetType> inDsetType)
{
  m_inDsrc = inDsrc;
  m_inDsetName = inDsetName;
  m_inDsetType = std::move(inDsetType);
}

void te::mnt::TINCreateIsolines::setOutput(te::da::DataSourcePtr outDsrc, std::string dsname)
{
  m_outDsrc = outDsrc;
  m_outDsetName = dsname;

}

void te::mnt::TINCreateIsolines::setParams(std::vector<double> &val, std::vector<double> &gval,  double tol)
{
  m_values = val;
  m_guidevalues = gval;
  m_tol = tol;
}