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
{
  namespace deprecated
  {
    // template class ConstantConvolutionWeight
    template <typename TDistance>
    class ConstantConvolutionWeights
    {
      // ----------------------- Standard services ------------------------------
    public:

      typedef TDistance Distance;


      inline
      double operator()(const Distance &/*aDisplacement*/) 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)