ArithmeticOperations.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/rp/ArithmeticOperations.h
  \brief Performs arithmetic operation over raster data.
 */

#ifndef __TERRALIB_RP_INTERNAL_ARITHMETICOPERATIONS_H
#define __TERRALIB_RP_INTERNAL_ARITHMETICOPERATIONS_H

#include "Algorithm.h"
#include "../raster/Grid.h"
#include "../raster/Interpolator.h"
#include "../common/progress/TaskProgress.h"

//Boost
#include <boost/math/constants/constants.hpp>

#include <map>
#include <memory>
#include <stack>
#include <string>
#include <vector>
#include <cmath>

namespace te
{
  namespace rst
  {
    class Raster;
  }

  namespace rp
  {
    /*!
      \class ArithmeticOperations
      
      \brief Performs arithmetic operation over raster data.
      
      \details The expression (arithmetic stringIt) must be formed by a combination of the following elements, 
               separated by spaces:

        Operators:
        +(Sum of Rasters or a raster and a real number),
        -(Subtraction of Rasters or a raster and a real number),
        *(Multiplication of Rasters or a raster and a real number),
        /(Division of Rasters or a raster and a real number)
        ^(Exponential of a raster and a real number)
        sqrt() (Square root of a raster)
        sin() (Sine of a raster)
        asin() (Arc Sine of a raster)
        cos() (Cosine of a raster)
        acos() (Arc Cosine of a raster)
        log() (Common (base-10) logarithm of a raster)
        tan() (Tangent of a raster)
        atan() (Arc tangent of a raster)
        ln() (Natural (base e) logarithm of a raster)

        Real Numbers (negative numbers must follow the form "-1.0")

        Raster bands: R0:1, R0:1, R1:0, .... (R0:1 is a reference to the first raster - with index 0 - from feeder, second band - with index 1).
      
      \ingroup rp_gen
     */
    class TERPEXPORT ArithmeticOperations : public Algorithm
    {
      public:

        /*!
          \class InputParameters
          \brief ArithmeticOperations input parameters
         */        
        class TERPEXPORT InputParameters : public AlgorithmInputParameters
        {
          public:

            std::vector< te::rst::Raster* > m_inputRasters; //!< Input rasters vector.

            std::string m_arithmeticString;  //!< Arithmetic string.

            bool m_normalize; //!< Output values normalization will be performed to fit the original input raster values range (default:true).

            bool m_enableProgress; //!< Enable/Disable the progress interface (default:false).
            
            te::rst::Interpolator::Method m_interpMethod; //!< The raster interpolator method (default:NearestNeighbor).
            
            InputParameters();
            
            InputParameters( const InputParameters& );
            
            ~InputParameters();
            
            //overload
            void reset() _NOEXCEPT_OP(false);
            
            //overload
            const  InputParameters& operator=( const InputParameters& params );
            
            //overload
            AbstractParameters* clone() const;
        };

        /*!
          \class OutputParameters
          \brief ArithmeticOperations output parameters
         */        
        class TERPEXPORT OutputParameters : public AlgorithmOutputParameters
        {
          public:
            
            std::string m_rType; //!< Output raster data source type (as described in te::raster::RasterFactory ).
            
            std::map< std::string, std::string > m_rInfo; //!< The necessary information to create the output rasters (as described in te::raster::RasterFactory). 
            
            std::unique_ptr< te::rst::Raster > m_outputRasterPtr; //!< The generated output registered raster.
            
            OutputParameters();
            
            OutputParameters( const OutputParameters& );
            
            ~OutputParameters();
            
            //overload
            void reset() _NOEXCEPT_OP(false);
            
            //overload
            const  OutputParameters& operator=( const OutputParameters& params );
            
            //overload
            AbstractParameters* clone() const;
        };

        ArithmeticOperations();
        
        ~ArithmeticOperations();
       
        //overload
        bool execute( AlgorithmOutputParameters& outputParams ) _NOEXCEPT_OP(false);
        
        //overload
        void reset() _NOEXCEPT_OP(false);
        
        //overload
        bool initialize( const AlgorithmInputParameters& inputParams ) _NOEXCEPT_OP(false);
        
        bool isInitialized() const;

      protected:

        /*!
          \class ExecStackElement
          \brief Execution stack element.
         */ 
        class TERPEXPORT ExecStackElement
        {
          public :
        
            bool m_isRaster; //!< true if this is a raster pointer element.

            bool m_isRealNumber; //!< true if this is a real number element.

            bool m_isBinaryOp; //!< true if this is a binary operator pointer element.

            bool m_isUnaryOp; //!< true if this is a unary operator pointer element.

            double m_realNumberValue; //!< Real number value.

            unsigned int m_rasterBand; //!< Raster band index.

            te::rst::Raster* m_rasterNPtr; //!< Raster pointer.

            std::string m_binaryOp; //!< Binary operator.

            std::string m_unaryOp; //!< Unary operator.
            
            mutable std::unique_ptr< te::rst::Raster > m_rasterHandler; //!< Raster handler.

            ExecStackElement( const ExecStackElement& rhs )
              : m_isRaster(false), m_isRealNumber(false), m_isBinaryOp(false), m_isUnaryOp(false),
              m_realNumberValue(0), m_rasterBand(0), m_rasterNPtr(0), m_binaryOp(""), m_unaryOp("")
            {
              operator=( rhs );
            }
            
            ExecStackElement()
              : m_isRaster(false), m_isRealNumber(false), m_isBinaryOp(false), m_isUnaryOp(false),
              m_realNumberValue(0), m_rasterBand(0), m_rasterNPtr(0), m_binaryOp(""), m_unaryOp("")
            {}
            
            ExecStackElement& operator=( const ExecStackElement& rhs )
            {
              m_isRaster = rhs.m_isRaster;
              m_isRealNumber = rhs.m_isRealNumber;
              m_isBinaryOp = rhs.m_isBinaryOp;
              m_isUnaryOp = rhs.m_isUnaryOp;
              m_realNumberValue = rhs.m_realNumberValue;
              m_rasterBand = rhs.m_rasterBand;
              m_rasterNPtr = rhs.m_rasterNPtr;
              m_binaryOp = rhs.m_binaryOp;
              m_unaryOp = rhs.m_unaryOp;
              m_rasterHandler.reset( rhs.m_rasterHandler.release() );
              
              return *this;
            }
          
            ~ExecStackElement() {}
        };

        typedef std::stack< ExecStackElement > ExecStackT; //!< Execution stack type definition.
        
        /*!
          \brief Type definition for a operation function pointer.
         */         
        typedef void (ArithmeticOperations::*BinOpFuncPtrT)( const double& inputValue1, const double& inputValue2,
          double& outputValue ) const;

        typedef void (ArithmeticOperations::*UnaryOpFuncPtrT)(const double& inputValue1, double& outputValue) const;

        ArithmeticOperations::InputParameters m_inputParameters; //!< Input execution parameters.

        bool m_isInitialized; //!< Tells if this instance is initialized.
        
        /*!
          \brief Addition binary operator function.
         */             
        inline void additionBinOp( const double& inputValue1, const double& inputValue2,
          double& outputValue ) const { outputValue = inputValue1 + inputValue2; }
      
        /*!
          \brief Subtraction binary operator function.
         */             
        inline void subtractionBinOp( const double& inputValue1, const double& inputValue2,
          double& outputValue ) const { outputValue = inputValue1 - inputValue2; }

        /*!
          \brief Multiplication binary operator function.
         */             
        inline void multiplicationBinOp( const double& inputValue1, const double& inputValue2,
          double& outputValue ) const { outputValue = inputValue1 * inputValue2; }
          
        /*!
          \brief Division binary operator function.
         */             
        inline void divisionBinOp( const double& inputValue1, const double& inputValue2,
          double& outputValue ) const {
          ( inputValue2 == 0.0 ) ? ( outputValue = 0 ) 
          : ( outputValue = inputValue1 / inputValue2 ); 
        }

        /*!
          \brief Exponencial binary operator function.
         */             
        inline void exponencialBinOp( const double& inputValue1, const double& inputValue2,
          double& outputValue ) const { outputValue = pow(inputValue1, inputValue2); }

        /*!
        \brief Square root unary operator function.
        */
        inline void sqrtUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = sqrt(inputValue); }

        /*!
        \brief Sine unary operator function.
        */
        inline void sinUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = sin(inputValue); }

        /*!
        \brief Arc sine unary operator function.
        */
        inline void asinUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = asin(inputValue); }

        /*!
        \brief Cosine unary operator function.
        */
        inline void cosUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = cos(inputValue); }

        /*!
        \brief Arc cosine unary operator function.
        */
        inline void acosUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = acos(inputValue); }

        /*!
        \brief Common logarithm unary operator function.
        */
        inline void logUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = log10(inputValue); }

        /*!
        \brief Tangent unary operator function.
        */
        inline void tanUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = tan(inputValue); }

        /*!
        \brief Arc Tangent unary operator function.
        */
        inline void atanUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = atan(inputValue); }

        /*!
        \brief  unary operator function.
        */
        inline void lnUnaryOp(const double& inputValue,
          double& outputValue) const { outputValue = log(inputValue); }
          
        /*!
          \brief Execute the automata parsing the given input string.
          \param aStr The input arithmetic expression string.
          \param inRasters Input rasters pointers.
          \param outRaster Output raster pointer (pre-initiated).
          \param generateOutput If true, the output raster data will be generated, if false only the automata execution will be performed..
          \param progressPtr A pointer to a progress interface to be pulsed on each operation or a null pointer.
        */
        bool executeString( const std::string& aStr, 
          const std::vector< te::rst::Raster* >& inRasters,
          std::unique_ptr<te::rst::Raster>& outRaster,
          bool generateOutput,
          te::common::TaskProgress* const progressPtr ) const;
        
        /*!
          \brief Convert the input tokens vector from the infix notation to postfix notation.
          \param input The input tokens vector.
          \param output The output tokens vector.
        */
        void inFix2PostFix( const std::vector< std::string >& input, 
          std::vector< std::string >& output ) const;
        
        /*!
          \brief Print tokens to stout.
          \param input The input tokens vector.
        */
        void printTokens( const std::vector< std::string >& input ) const;        
        
        /*!
          \brief Returns true if the given token is an operator.
          \param inputToken Input token.
          \return Returns true if the given token is an operator.
        */         
        bool isOperator( const std::string& inputToken ) const;
      
        /*!
          \brief Returns true if the given token is a binary operator.
          \param inputToken Input token.
          \return Returns true if the given token is a binary operator.
        */               
        bool isBinaryOperator( const std::string& inputToken ) const;
      
        /*!
          \brief Returns true if the given token is a unary operator.
          \param inputToken Input token.
          \return Returns true if the given token is a unary operator.
        */            
        bool isUnaryOperator( const std::string& inputToken ) const;
      
        /*!
          \brief Returns true if operator1 has greater of equal precedence over operator2.
          \param operator1 Operator1 input token.
          \param operator2 Operator2 input token.
          \return Returns true if operator1 has greater of equal precedence over operator2.
        */      
        bool op1HasGreaterOrEqualPrecedence( const std::string& operator1, 
          const std::string& operator2 ) const;
        
        /*!
          \brief Returns true if the given token is a raster data token.
          \param token Input token.
          \param rasterIdx The output converted raster index value.
          \param bandIdx The output converted band index value.
          \return Returns true if the given token is a raster data token.
        */           
        bool isRasterBandToken( const std::string& token, unsigned int& rasterIdx,
          unsigned int& bandIdx ) const;
      
        /*!
          \brief Execute the given binary operator using the current given execution stack.
          \param token Operator token.
          \param execStack Execution stack.
          \param generateOutput if true the execution will generate valid output data, if false only dummy stack elements will be generated.
          \return true if OK, false on errors..
        */       
        bool execBinaryOperator( const std::string& token, ExecStackT& 
          execStack, bool generateOutput ) const;
          
        /*!
          \brief Execute the given binary operator using the given input rasters
          \param inRaster1 Input raster 1.
          \param band1 Input raster 1 band.
          \param inRaster2 Input raster 2.
          \param band2 Input raster 2 band.
          \param binOptFunctPtr The binary operation function pointer.
          \param outRasterPtr The generated output raster.
          \return true if OK, false on errors..
        */       
        bool execBinaryOperatorRasterXRaster( const te::rst::Raster& inRaster1, 
          const unsigned int band1, const te::rst::Raster& inRaster2, 
          const unsigned int band2,
          const BinOpFuncPtrT binOptFunctPtr,
          std::unique_ptr<te::rst::Raster>& outRasterPtr ) const; 
          
        /*!
          \brief Execute the given binary operator using the given input raster and a real number
          \param inRaster Input raster.
          \param bandIdx Input raster band.
          \param binOptFunctPtr The binary operation function pointer.
          \param outRasterPtr The generated output raster.
          \param realNumberIsRigthtTerm true if the real number is the right term.
          \return true if OK, false on errors..
        */       
        bool execBinaryOperatorRasterXReal( const te::rst::Raster& inRaster, 
          const unsigned int bandIdx, const double value, 
          const BinOpFuncPtrT binOptFunctPtr,
          std::unique_ptr<te::rst::Raster>& outRasterPtr,
          const bool realNumberIsRigthtTerm ) const;  
      
        /*!
          \brief Execute the given unary operator using the current given execution stack.
          \param token Operator token.
          \param execStack Execution stack.
          \param generateOutput if true the execution will generate valid output data, if false only dummy stack elements will be generated.
          \return true if OK, false on errors..
        */         
        bool execUnaryOperator( const std::string& token, ExecStackT& 
          execStack, bool generateOutput ) const;

        /*!
        \brief Execute the given unary operator using the given input raster
        \param inRaster Input raster.
        \param band1 Input raster band.
        \param unaryOptFunctPtr The unary operation function pointer.
        \param outRasterPtr The generated output raster.
        \return true if OK, false on errors..
        */
        bool execUnaryOperatorRaster(const te::rst::Raster& inRaster,
          const unsigned int band, const UnaryOpFuncPtrT unaryOptFunctPtr,
          std::unique_ptr<te::rst::Raster>& outRasterPtr) const;

        /*!
        \brief Execute the given unary operator using the given a real number
        \param unaryOptFunctPtr The unary operation function pointer.
        \param outRasterPtr The generated output raster.
        \param realNumberIsRigthtTerm true if the real number is the right term.
        \return true if OK, false on errors..
        */
        bool execUnaryOperatorReal(const double value,
          const UnaryOpFuncPtrT unaryOptFunctPtr,
          std::unique_ptr<te::rst::Raster>& outRasterPtr,
          const bool realNumberIsRigthtTerm) const;
        
        /*!
          \brief Returns true if the given token is a real number.
          \param token Input token.
          \param realValue The output converted value.
          \return Returns true if the given token is a real number.
        */         
        bool isRealNumberToken( const std::string& token, double& realValue ) const;
      
        /*!
          \brief Allocate a new RAM memory raster.
          \param grid The output raster grid.
          \param rasterPtr The output raster pointer.
          \return Returns true if OK, false on errors.
        */        
        bool allocResultRaster( const te::rst::Grid& grid, 
          std::unique_ptr<te::rst::Raster>& rasterPtr ) const;
        
        /*!
          \brief Split the input string into a vector of token strings
          \param inputStr The input string.
          \param outTokens The generated output tokens vector.
        */         
        void getTokensStrs( const std::string& inputStr,
          std::vector< std::string >& outTokens ) const;
    };
  } // end namespace rp
}   // end namespace te

#endif // __TERRALIB_RP_INTERNAL_ARITHMETICOPERATIONS_H