otbOBIAMuParserFunctor.h

Go to the documentation of this file.

/*
 * Copyright (C) 1999-2011 Insight Software Consortium
 * Copyright (C) 2005-2022 Centre National d'Etudes Spatiales (CNES)
 *
 * This file is part of Orfeo Toolbox
 *
 *     https://www.orfeo-toolbox.org/
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
 
 
#ifndef otbOBIAMuParserFunctor_h
#define otbOBIAMuParserFunctor_h
 
#include "otbParser.h"
#include "otbMacro.h"
#include <string>
 
 
namespace otb
{
namespace Functor
{
 
template <class TLabelObject>
// class ITK_EXPORT OBIAMuParserFunctor : public itk::LightObject
class ITK_EXPORT OBIAMuParserFunctor
{
public:
  typedef Parser              ParserType;
  typedef OBIAMuParserFunctor Self;
 
  std::string GetNameOfClass()
  {
    return "OBIAMuParserFunctor";
  }
 
  inline bool operator()(const TLabelObject& a)
  {
 
    double value;
 
    if (a.GetNumberOfAttributes() != m_AAttributes.size())
    {
      this->SetAttributes(a);
    }
 
    for (unsigned int i = 0; i < m_AAttributes.size(); ++i)
    {
 
      std::string name = (m_AttributesName[i]);
      m_AAttributes[i] = a.GetAttribute(name.c_str());
    }
 
    try
    {
      value = m_Parser->Eval();
    }
    catch (itk::ExceptionObject& err)
    {
      itkExceptionMacro(<< err);
    }
    return static_cast<bool>(value);
  }
 
  void SetExpression(const std::string expression)
  {
    m_Expression = expression;
    m_Parser->SetExpr(m_Expression);
  }
 
  std::string GetExpression() const
  {
    return m_Expression;
  }
 
  void ParseAttributeName(std::string& attributeName)
  {
 
    for (unsigned int i = 0; i < attributeName.size(); ++i)
    {
      if (attributeName[i] == ':')
      {
        attributeName.erase(i, 1);
        attributeName[i] = '_';
      }
    }
    // TODO JGU
    // replace "Band" by "b" for homogeneity with other functors
  }
 
  void SetAttributes(const TLabelObject& a)
  {
 
    unsigned int nbOfAttributes = a.GetNumberOfAttributes();
 
    m_AAttributes.resize(nbOfAttributes, 0.0);
    m_AttributesName.resize(nbOfAttributes, "");
    m_AttributesName = a.GetAvailableAttributes();
    for (unsigned int i = 0; i < nbOfAttributes; ++i)
    {
      std::string attributeName = m_AttributesName.at(i);
      ParseAttributeName(attributeName); // eliminate '::' from string name
      m_Parser->DefineVar(attributeName, &(m_AAttributes[i]));
    }
  }
 
  void SetAttributes(std::vector<std::string> shapeAttributes, std::vector<std::string> statAttributes, unsigned int nbOfBands)
  {
    int          index          = 0;
    unsigned int nbOfAttributes = shapeAttributes.size() + statAttributes.size() * nbOfBands;
 
    m_AAttributes.resize(nbOfAttributes, 0.0);
    m_AttributesName.resize(nbOfAttributes, "");
    std::ostringstream varName;
    for (unsigned int i = 0; i < shapeAttributes.size(); ++i)
    {
 
      varName << "SHAPE::" << shapeAttributes.at(i);
      m_AttributesName.at(index) = varName.str();
      varName.str("");
      varName << "SHAPE_" << shapeAttributes.at(i);
 
      m_Parser->DefineVar(varName.str(), &(m_AAttributes[index]));
      varName.str("");
      index++;
    }
    for (unsigned int i = 0; i < statAttributes.size(); ++i)
    {
      for (unsigned int bandIndex = 1; bandIndex <= nbOfBands; bandIndex++)
      {
        varName << "STATS::Band" << bandIndex << "::" << statAttributes.at(i);
        m_AttributesName.at(index) = varName.str();
        varName.str("");
        varName << "STATS_Band" << bandIndex << "_" << statAttributes.at(i);
        m_Parser->DefineVar(varName.str(), &(m_AAttributes[index]));
        varName.str("");
        index++;
      }
    }
  }
 
  bool CheckExpression()
  {
    return m_Parser->CheckExpr();
  }
 
  const std::map<std::string, Parser::ValueType*>& GetVar() const
  {
    return this->m_Parser->GetVar();
  }
 
  Parser::FunctionMapType GetFunList() const
  {
    return this->m_Parser->GetFunList();
  }
 
  OBIAMuParserFunctor()
  {
    m_Parser = ParserType::New();
    m_AAttributes.resize(0);
  };
 
  ~OBIAMuParserFunctor(){};
 
protected:
private:
  OBIAMuParserFunctor(const Self&) = delete;
  void operator=(const Self&) = delete;
 
  std::string              m_Expression;
  ParserType::Pointer      m_Parser;
  std::vector<double>      m_AAttributes;
  std::vector<std::string> m_AttributesName;
  double                   m_ParserResult;
};
} // end of Functor namespace
 
 
} // end namespace otb
 
 
#endif