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)