ProjectiveGT.h

Go to the documentation of this file.

/*  Copyright (C) 2008 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/geometry/ProjectiveGT.h

  \brief 2D Projective Geometric transformation.
*/

#ifndef __TERRALIB_GEOMETRY_INTERNAL_PROJECTIVEGT_H
#define __TERRALIB_GEOMETRY_INTERNAL_PROJECTIVEGT_H

// TerraLib
#include "Config.h"
#include "GeometricTransformation.h"

// Boost
#include <boost/concept_check.hpp>
#include <boost/numeric/ublas/matrix.hpp>

namespace te
{
  namespace gm
  {
    /*!
      \class ProjectiveGT

      \brief 2D Projective Geometric transformation.
      
      \note The transformation parameters (inside GTParameters::m_directParameters) are disposed on the following form:

      \code
      | u |   | a b c | * | x |
      | v | = | d e f |   | y |
      | 1 |   | g h 1 |   | 1 |

      u = a.x + b.y + c
          -------------------
          g.x + h.y + 1

      v = d.x + e.y + f
          -------------------
          g.x + h.y + 1

      |u0|   |x0   y0  1  0   0  0  -x0*u0 -y0*v0|.|a|
      |v0|   |0    0   0 x0  y0  1  -x0*v0 -y0*v0| |b|
      |u1|   |...................................| |c|
      |v1| = |...................................| |d|
      |..|   |...................................| |e|
      |..|   |...................................| |f|
      |..|   |...................................| |g|
      |..|   |...................................| |h|

      GTParameters::m_directParameters = [ a b c d e f g h ]
      \endcode
    */
    class TEGEOMEXPORT ProjectiveGT : public GeometricTransformation
    {
      public:

        /*! \brief Default constructor. */
        ProjectiveGT();

        /*! \brief Destructor. */
        ~ProjectiveGT();

        const std::string& getName() const;

        bool isValid( const GTParameters& params ) const;

        void directMap( const GTParameters& params, const double& pt1X, 
                        const double& pt1Y, double& pt2X, double& pt2Y ) const;
          
        void inverseMap( const GTParameters& params, const double& pt2X, 
                         const double& pt2Y, double& pt1X, double& pt1Y ) const;

        unsigned int getMinRequiredTiePoints() const;

        GeometricTransformation* clone() const;

      protected:
        
        // Variables used by method computeParameters
        mutable unsigned int m_computeParameters_tiepointsSize;
        mutable unsigned int m_computeParameters_index1;
        mutable unsigned int m_computeParameters_index2;
        mutable unsigned int m_computeParameters_blockOffset;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_L_DM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_A_DM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_L_IM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_A_IM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_A_DM_t;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_A_IM_t;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_N_DM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_N_IM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_U_DM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_U_IM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_N_DM_inv;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_N_IM_inv;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_X_DM;
        mutable boost::numeric::ublas::matrix< double > m_computeParameters_X_IM;

        bool computeParameters( GTParameters& params ) const;
    };

  } // end namespace gm
}   // end namespace te

#endif  // __TERRALIB_GEOMETRY_INTERNAL_PROJECTIVEGT_H