TrueDigitalSurfaceLocalEstimator.h

 
#pragma once

 
#if defined(TrueDigitalSurfaceLocalEstimator_RECURSES)
#error Recursive header files inclusion detected in TrueDigitalSurfaceLocalEstimator.h
#else // defined(TrueDigitalSurfaceLocalEstimator_RECURSES)
#define TrueDigitalSurfaceLocalEstimator_RECURSES
 
#if !defined TrueDigitalSurfaceLocalEstimator_h
#define TrueDigitalSurfaceLocalEstimator_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/base/ConstAlias.h"
#include "DGtal/base/CountedConstPtrOrConstPtr.h"
#include "DGtal/topology/CanonicSCellEmbedder.h"
 
namespace DGtal
{

  // template class TrueDigitalSurfaceLocalEstimator
  template <typename TKSpace, typename TShape, typename TGeometricFunctor>
  class TrueDigitalSurfaceLocalEstimator
  {
    BOOST_CONCEPT_ASSERT(( concepts::CCellularGridSpaceND< TKSpace > ));
 
  public:
    typedef TrueDigitalSurfaceLocalEstimator<TKSpace,TShape,TGeometricFunctor> Self;
    typedef TKSpace KSpace;
    typedef TShape Shape;
    typedef TGeometricFunctor GeometricFunctor;
    typedef typename KSpace::Space Space;
    typedef typename KSpace::SCell SCell;
    typedef typename KSpace::Surfel Surfel;
    typedef typename Space::RealPoint RealPoint;
    typedef typename RealPoint::Coordinate Scalar;
    typedef typename GeometricFunctor::Quantity Quantity;
    typedef CanonicSCellEmbedder<KSpace> SCellEmbedder;
    BOOST_CONCEPT_ASSERT(( concepts::CUnaryFunctor< GeometricFunctor, RealPoint, Quantity > ));
 
    // ----------------------- Standard services ------------------------------
  public:

    ~TrueDigitalSurfaceLocalEstimator();

    TrueDigitalSurfaceLocalEstimator();

    TrueDigitalSurfaceLocalEstimator( const Self & other );

    Self& operator=( const Self & other );

    Scalar h() const;

    void attach( ConstAlias<Shape> aShape );

    void setParams( ConstAlias<KSpace> ks, 
                    Clone<GeometricFunctor> fct,
                    const int maxIter = 20, 
                    const Scalar accuracy = 0.0001, 
                    const Scalar gamma = 0.5 );

    template <typename SurfelConstIterator>
    void init( const Scalar _h, 
               SurfelConstIterator itb,
               SurfelConstIterator ite );
    
    template <typename SurfelConstIterator>
    Quantity eval( SurfelConstIterator it ) const;

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

    RealPoint embed( Surfel surfel ) const;
 
    // ----------------------- Interface --------------------------------------
  public:

    void selfDisplay ( std::ostream & out ) const;

    bool isValid() const;
 
    // ------------------------- Protected Datas ------------------------------
  protected:
    CountedConstPtrOrConstPtr<KSpace> myKSpace;
    CowPtr<GeometricFunctor> myFct;
    SCellEmbedder myEmbedder;
    CountedConstPtrOrConstPtr<Shape> myShape;
    Scalar myH;
    bool myNearest;

    int myMaxIter;
    Scalar myAccuracy;
    Scalar myGamma;
 
    // ------------------------- Private Datas --------------------------------
  private:
 
    // ------------------------- Hidden services ------------------------------
  protected:
 
  private:
 
 
    // ------------------------- Internals ------------------------------------
  private:
 
  }; // end of class TrueDigitalSurfaceLocalEstimator
 

  template <typename TKSpace, typename TShape, typename TGeometricFunctor>
  std::ostream&
  operator<< ( std::ostream & out, const TrueDigitalSurfaceLocalEstimator<TKSpace, TShape, TGeometricFunctor> & object );
 
} // namespace DGtal
 

// Includes inline functions.
#include "DGtal/geometry/surfaces/estimation/TrueDigitalSurfaceLocalEstimator.ih"
 
//                                                                           //
 
#endif // !defined TrueDigitalSurfaceLocalEstimator_h
 
#undef TrueDigitalSurfaceLocalEstimator_RECURSES
#endif // else defined(TrueDigitalSurfaceLocalEstimator_RECURSES)