LambdaMST2D.h

#pragma once

#if defined(LAMBDAMST2D_RECURSES)
#error Recursive header files inclusion detected in LambdaMST2D.h
#else // defined(LAMBDAMST2D_RECURSES)

#define LAMBDAMST2D_RECURSES

#if !defined LAMBDAMST2D_h

#define LAMBDAMST2D_h

#include <algorithm>
#include <iterator>
#include <cmath>
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/kernel/CSpace.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/geometry/curves/estimation/FunctorsLambdaMST.h"
#include "DGtal/geometry/curves/CForwardSegmentComputer.h"
#include "DGtal/geometry/curves/estimation/CLMSTTangentFromDSS.h"

namespace DGtal {
  template < typename TSpace, typename TSegmentation, typename Functor >
  class LambdaMST2DEstimator
  {
    //Checking concepts
    BOOST_STATIC_ASSERT(( TSpace::dimension == 2 ));
    BOOST_CONCEPT_ASSERT(( concepts::CSpace<TSpace> ));
    BOOST_CONCEPT_ASSERT(( concepts::CLMSTTangentFromDSS<Functor> ));
    BOOST_CONCEPT_ASSERT(( concepts::CForwardSegmentComputer<typename TSegmentation::SegmentComputer> ));
    // ----------------------- Types ------------------------------
  public:
    typedef TSegmentation Segmentation;
    typedef typename TSegmentation::SegmentComputer SegmentComputer;
    typedef typename SegmentComputer::ConstIterator ConstIterator;
    typedef typename Functor::Value Value;
    typedef typename TSpace::RealVector RealVector;
    typedef typename TSpace::Point Point;
    
    // ----------------------- Interface --------------------------------------
  public:
    LambdaMST2DEstimator();
    
    void init ( const ConstIterator & itb, const ConstIterator & ite );
    
    void attach ( Alias<TSegmentation> SegmentComputer );
    
    RealVector eval ( const ConstIterator & it );
    
    template <typename OutputIterator>
    OutputIterator eval ( const ConstIterator & itb, const ConstIterator & ite, 
                         OutputIterator result );
    
    // ----------------------- Standard services ------------------------------
  public:
    
    bool isValid() const;
    // ------------------------- Internals ------------------------------------
  protected:
    
    template <typename OutputIterator>
    void accumulate ( std::vector < Value > & outValues, OutputIterator & result );
    
    // ------------------------- Private Datas --------------------------------
  private:
    
    ConstIterator myBegin;
    ConstIterator myEnd;
    Functor myFunctor;
    TSegmentation * dssSegments;
    
  }; // end of class LambdaMST2DEstimator
  
  //-------------------------------------------------------------------------------------------
  
  // Template class LambdaMST2D
  template < typename DSSSegmentationComputer, typename LambdaFunction = functors::Lambda64Function >
  class LambdaMST2D: 
  public LambdaMST2DEstimator < Z2i::Space, DSSSegmentationComputer,
    TangentFromDSS2DFunctor < typename DSSSegmentationComputer::SegmentComputer, LambdaFunction > >
    {
      typedef LambdaMST2DEstimator < Z2i::Space, DSSSegmentationComputer,
      TangentFromDSS2DFunctor < typename DSSSegmentationComputer::SegmentComputer, LambdaFunction> > Super;
      
    public: 
      LambdaMST2D (): Super() {}
    };
}// namespace DGtal

// Includes inline functions
#include "DGtal/geometry/curves/estimation/LambdaMST2D.ih"

//                                                                           //

#endif // !defined LAMBDAMST2D_h

#undef LAMBDAMST2D_RECURSES
#endif // else defined(LAMBDAMST2D_RECURSES)