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 Data ------------------------------
  private:
    // ------------------------- Private Data --------------------------------
  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)