ElementaryConvolutionNormalVectorEstimator.h

 
#pragma once

 
#if defined(ElementaryConvolutionNormalVectorEstimator_RECURSES)
#error Recursive header files inclusion detected in ElementaryConvolutionNormalVectorEstimator.h
#else // defined(ElementaryConvolutionNormalVectorEstimator_RECURSES)
#define ElementaryConvolutionNormalVectorEstimator_RECURSES
 
#if !defined ElementaryConvolutionNormalVectorEstimator_h
#define ElementaryConvolutionNormalVectorEstimator_h

// Inclusions
#include <iostream>
#include <DGtal/base/Common.h>
#include <DGtal/topology/SCellsFunctors.h>
#include <vector>

 
namespace DGtal
{
  namespace functors
  {
  // template class ElementaryConvolutionNormalVectorEstimator
  template <typename TSurfel, typename TEmbedder>
  class ElementaryConvolutionNormalVectorEstimator
  {
  public:
 
    typedef TSurfel Surfel;
    typedef TEmbedder SCellEmbedder;
    typedef typename SCellEmbedder::RealPoint RealPoint;
    typedef RealPoint Quantity;

    ElementaryConvolutionNormalVectorEstimator(ConstAlias<SCellEmbedder> anEmbedder,
                                               const double h):
      myEmbedder(&anEmbedder), myH(h)
    {
    }

    void pushSurfel(const Surfel & aSurf,
                    const double aDistance)
    {
      RealPoint elementary;
      Dimension i = myEmbedder->space().sOrthDir ( aSurf );
      elementary[ i ] = myEmbedder->space().sDirect ( aSurf, i ) ? 1 : -1;
 
      myWeightedVector += aDistance *elementary;
    }

    Quantity eval( )
    {
      return myWeightedVector.getNormalized();
    }

    void reset()
    {
      myWeightedVector = RealPoint().diagonal(0.0);
    }
 
 
  private:

    const SCellEmbedder * myEmbedder;

    RealPoint myWeightedVector;
 
    //Grid step
    double myH;
 
  }; // end of class ElementaryConvolutionNormalVectorEstimator
  } //namespace functors
} // namespace DGtal
 
 
//                                                                           //
 
#endif // !defined ElementaryConvolutionNormalVectorEstimator_h
 
#undef ElementaryConvolutionNormalVectorEstimator_RECURSES
#endif // else defined(ElementaryConvolutionNormalVectorEstimator_RECURSES)