MostCenteredMaximalSegmentEstimator.h

#pragma once

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

#define MostCenteredMaximalSegmentEstimator_RECURSES

#if !defined MostCenteredMaximalSegmentEstimator_h

#define MostCenteredMaximalSegmentEstimator_h

// Inclusions
#include <iostream>
#include <list>

#include "DGtal/base/Common.h"
#include "DGtal/base/Exceptions.h"
#include "DGtal/base/Circulator.h"

#include "DGtal/geometry/curves/estimation/CSegmentComputerEstimator.h"
#include "DGtal/geometry/curves/CForwardSegmentComputer.h"
#include "DGtal/geometry/curves/SaturatedSegmentation.h"


namespace DGtal
{

  // template class MostCenteredMaximalSegmentEstimator
  template <typename SegmentComputer, typename SCEstimator>
  class MostCenteredMaximalSegmentEstimator
  {

    BOOST_CONCEPT_ASSERT(( CForwardSegmentComputer<SegmentComputer> )); 
    BOOST_CONCEPT_ASSERT(( CSegmentComputerEstimator<SCEstimator> )); 
    BOOST_STATIC_ASSERT(( boost::is_same< SegmentComputer, 
                          typename SCEstimator::SegmentComputer >::value ));

    // ----------------------- Types ------------------------------
  public:

    typedef typename SegmentComputer::ConstIterator ConstIterator;
    typedef typename SCEstimator::Quantity Quantity;

    typedef SaturatedSegmentation<SegmentComputer> Segmentation; 
    typedef typename Segmentation::SegmentComputerIterator SegmentIterator; 

    // ----------------------- Standard services ------------------------------
  public:

    MostCenteredMaximalSegmentEstimator();

    MostCenteredMaximalSegmentEstimator(const SegmentComputer& aSegmentComputer, 
                                        const SCEstimator& aSCEstimator);

    ~MostCenteredMaximalSegmentEstimator() {};

    // ----------------------- Interface --------------------------------------
  public:

    void init(const double h, const ConstIterator& itb, const ConstIterator& ite);

    Quantity eval(const ConstIterator& it);

    template <typename OutputIterator>
    OutputIterator eval(const ConstIterator& itb, const ConstIterator& ite, 
                        OutputIterator result); 


    bool isValid() const;

    // ------------------------- Protected Datas ------------------------------
  protected:

    // ------------------------- Private Datas --------------------------------
  private:

    double myH; 

    ConstIterator myBegin,myEnd;

    SegmentComputer mySC; 

    SCEstimator mySCEstimator;

    // ------------------------- Internal services ------------------------------

  private:

    template <typename OutputIterator>
    OutputIterator endEval(const ConstIterator& itb, const ConstIterator& ite, ConstIterator& itCurrent, 
                           SegmentIterator& first, SegmentIterator& last, 
                           OutputIterator result); 

    template <typename OutputIterator>
    OutputIterator endEval(const ConstIterator& /*itb*/, const ConstIterator& ite, ConstIterator& itCurrent, 
                           SegmentIterator& /*first*/, SegmentIterator& last, 
                           OutputIterator result, IteratorType); 

    template <typename OutputIterator>
    OutputIterator endEval(const ConstIterator& itb, const ConstIterator& ite, ConstIterator& itCurrent, 
                           SegmentIterator& first, SegmentIterator& last, 
                           OutputIterator result, CirculatorType); 

    // ------------------------- Hidden services ------------------------------

  private:

    MostCenteredMaximalSegmentEstimator ( const MostCenteredMaximalSegmentEstimator & other );

    MostCenteredMaximalSegmentEstimator & operator= ( const MostCenteredMaximalSegmentEstimator & other );


  }; // end of class MostCenteredMaximalSegmentEstimator

} // namespace DGtal


// Includes inline functions.
#include "DGtal/geometry/curves/estimation/MostCenteredMaximalSegmentEstimator.ih"
#include "DGtal/geometry/curves/estimation/SegmentComputerEstimators.h"

//                                                                           //

#endif // !defined MostCenteredMaximalSegmentEstimator_h

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