LocalConvolutionNormalVectorEstimator.h

#pragma once

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

#define LocalConvolutionNormalVectorEstimator_RECURSES

#if !defined LocalConvolutionNormalVectorEstimator_h

#define LocalConvolutionNormalVectorEstimator_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/base/CountedPtr.h"
#include "DGtal/topology/BreadthFirstVisitor.h"
#include "DGtal/topology/DigitalSurface.h"
#include "DGtal/geometry/surfaces/estimation/CConvolutionWeights.h"

namespace DGtal
{
  // template class LocalConvolutionNormalVectorEstimator
  template <typename TDigitalSurface, typename TKernelFunctor>
  class LocalConvolutionNormalVectorEstimator
  {

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

    typedef TDigitalSurface DigitalSurface;
    typedef TKernelFunctor KernelFunctor;
    typedef DigitalSurface Surface;
    typedef typename Surface::ConstIterator ConstIterator;
    typedef typename Surface::KSpace::Space::RealVector Quantity;
    typedef typename Surface::SCell SCell;

    BOOST_CONCEPT_ASSERT(( CConvolutionWeights<TKernelFunctor>));

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

    LocalConvolutionNormalVectorEstimator(const DigitalSurface & aSurface,
                                          const KernelFunctor & aFunctor);

    ~LocalConvolutionNormalVectorEstimator() {};

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

    const Surface & surface() const;

    void init(const double h,
              const unsigned int radius);

    Quantity eval(const SCell & scell) const;

    Quantity eval(const ConstIterator& it) const;

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

    template <typename OutputIterator>
    OutputIterator evalAll( OutputIterator result ) const;


    bool isValid() const;

    // ------------------------- Protected Datas ------------------------------
  protected:
    LocalConvolutionNormalVectorEstimator()
    {
      myFlagIsInit = false;
    }


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

    double myH;

    bool myFlagIsInit;

    unsigned int myRadius;

    const DigitalSurface & mySurface;

    const KernelFunctor & myKernelFunctor;

    // ------------------------- Hidden services ------------------------------
  private:

    LocalConvolutionNormalVectorEstimator ( const LocalConvolutionNormalVectorEstimator & other );

    LocalConvolutionNormalVectorEstimator & operator= ( const LocalConvolutionNormalVectorEstimator & other );


  }; // end of class LocalConvolutionNormalVectorEstimator

} // namespace DGtal


// Includes inline functions.
#include "DGtal/geometry/surfaces/estimation/LocalConvolutionNormalVectorEstimator.ih"

//                                                                           //

#endif // !defined LocalConvolutionNormalVectorEstimator_h

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