BasicConvolutionWeights.h

#pragma once

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

#define BasicConvolutionWeights_RECURSES

#if !defined BasicConvolutionWeights_h

#define BasicConvolutionWeights_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
 #include "DGtal/kernel/NumberTraits.h"

namespace DGtal
{

  // template class ConstantConvolutionWeight
  template <typename TDistance>
  class ConstantConvolutionWeights
  {
    // ----------------------- Standard services ------------------------------
  public:
    
    typedef TDistance Distance;


    inline
    double operator()(const Distance &/*aDisplacment*/) const
    {
      return 1.0;
    }
  };

  // template class ConstantConvolutionWeight
  template <typename TDistance>
  class GaussianConvolutionWeights
  {
    // ----------------------- Standard services ------------------------------
  public:

    typedef TDistance Distance;

    GaussianConvolutionWeights(const double sigma): mySigma(sigma)
    {
      myCoef = 1.0/(mySigma * sqrt(2.0*M_PI));
      myCoef2 = 1.0/(2.0*M_PI);
    }

    inline
    double operator()(const Distance &aDisplacment) const
    {
      return myCoef*exp(-NumberTraits<Distance>::castToDouble(aDisplacment)*
        NumberTraits<Distance>::castToDouble(aDisplacment)*myCoef2);
    }

    double mySigma;

    double myCoef;
    double myCoef2;
  };




} // namespace DGtal


#endif // !defined BasicConvolutionWeights_h

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