DigitalMetricAdapter.h

#pragma once

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

#define DigitalMetricAdapter_RECURSES

#if !defined DigitalMetricAdapter_h

#define DigitalMetricAdapter_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/geometry/volumes/distance/CMetricSpace.h"
#include "DGtal/kernel/CInteger.h"
#include "DGtal/base/ConstAlias.h"
#include "DGtal/kernel/NumberTraits.h"

namespace DGtal
{

  // template class DigitalMetricAdapter
  template <typename TMetric, typename TInteger = typename TMetric::Space::Integer>
  class DigitalMetricAdapter
  {
    // ----------------------- Standard services ------------------------------
  public:

    typedef TMetric Metric;
    BOOST_CONCEPT_ASSERT(( concepts::CMetricSpace<Metric> ));

    typedef TInteger Value;
    BOOST_CONCEPT_ASSERT(( concepts::CInteger<Value> ));

    // ----------------------- CDigitalMetric------------------------------------
    typedef typename Metric::Point Point;

    typedef typename Metric::Space Space;

    typedef Value RawValue;
    
       
    DigitalMetricAdapter(ConstAlias<Metric> aMetric):
      myMetric(&aMetric) {}

    
    DigitalMetricAdapter(const DigitalMetricAdapter &anotherMetric)
    {
       if (this != &anotherMetric)
         this->myMetric = anotherMetric.myMetric;
    }
      

    ~DigitalMetricAdapter() {}

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

    // ----------------------- CDigitalMetricSpace --------------------------------------

    DigitalMetricAdapter & operator= ( const DigitalMetricAdapter & other )
    {
      myMetric = other.myMetric;
      return *this;
    }

    Value operator() (const Point &p, const Point &q) const
    {
      return static_cast<Value>( std::ceil ( NumberTraits<typename Metric::Value>::castToDouble( myMetric->operator()(p,q))));
    }

    Value rawDistance(const Point &p, const Point &q) const
    {
      return this->operator()(p,q);
    }

    DGtal::Closest closest(const Point &a,
                           const Point &p,
                           const Point &q) const
    {
      if (this->operator()(a,p) < this->operator()(a,q))
        return DGtal::ClosestFIRST;
      if (this->operator()(a,p) > this->operator()(a,q))
        return DGtal::ClosestSECOND;
      return DGtal::ClosestBOTH;
    }

  

    void selfDisplay ( std::ostream & out ) const
    {
      out << "[DigitalMetricAdapter] adapted from "<< *myMetric;
    }

    bool isValid() const { return (myMetric != NULL); }

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

    const Metric *myMetric;
    
    // ------------------------- Hidden services ------------------------------
  protected:

    DigitalMetricAdapter();

    // ------------------------- Internals ------------------------------------
  private:

  }; // end of class DigitalMetricAdapter


  template <typename T>
  std::ostream&
  operator<< ( std::ostream & out, const DigitalMetricAdapter<T> & object );

} // namespace DGtal

//                                                                           //

#endif // !defined DigitalMetricAdapter_h

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