ReverseDistanceTransformation.h

 
#pragma once
 
#if defined(ReverseDistanceTransformation_RECURSES)
#error Recursive header files inclusion detected in ReverseDistanceTransformation.h
#else // defined(ReverseDistanceTransformation_RECURSES)
#define ReverseDistanceTransformation_RECURSES
 
#if !defined ReverseDistanceTransformation_h
#define ReverseDistanceTransformation_h
 
// Inclusions
#include <iostream>
#include <vector>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/NumberTraits.h"
#include "DGtal/geometry/volumes/distance/CPowerSeparableMetric.h"
#include "DGtal/geometry/volumes/distance/PowerMap.h"
#include "DGtal/images/DefaultConstImageRange.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"
#include "DGtal/images/ImageContainerBySTLVector.h"
#include "DGtal/images/CImage.h"
#include "DGtal/base/ConstAlias.h"
 
namespace DGtal
{
 
  // template class ReverseDistanceTransformation
  template < typename TWeightImage,
             typename TPSeparableMetric,
             typename TImageContainer =
               ImageContainerBySTLVector<HyperRectDomain<typename TWeightImage::Domain::Space>,
             typename TWeightImage::Domain::Space::Vector>
           >
  class ReverseDistanceTransformation
    : public PowerMap<TWeightImage, TPSeparableMetric, TImageContainer>
  {
 
  public:
 
    typedef TWeightImage WeightImage;
 
    typedef TPSeparableMetric PowerSeparableMetric;
 
    typedef TImageContainer ImageContainer;
    BOOST_CONCEPT_ASSERT(( concepts::CImage<ImageContainer> ));
 
    typedef typename PowerSeparableMetric::Value Value;
 
    typedef typename TWeightImage::Domain::Space::Point Point;
 
    typedef typename TWeightImage::Domain::Space::Vector Vector;
 
    typedef typename PowerSeparableMetric::Weight Weight;
 
    typedef  ReverseDistanceTransformation<TWeightImage,
                                           TPSeparableMetric,
                                           TImageContainer> Self;
 
    typedef PowerMap<TWeightImage,TPSeparableMetric> Parent;
 
    typedef  DefaultConstImageRange<Self> ConstRange;
 
 
    typedef typename PowerMap< TWeightImage,TPSeparableMetric,
                               TImageContainer >::Domain  Domain;
 
 
    ReverseDistanceTransformation(ConstAlias<Domain> aDomain,
                                  ConstAlias<WeightImage> aWeightImage,
                                  ConstAlias<PowerSeparableMetric> aMetric):
      PowerMap<TWeightImage,TPSeparableMetric,TImageContainer>(aDomain,
                                                               aWeightImage,
                                                               aMetric)
    {}
 
    ReverseDistanceTransformation(ConstAlias<Domain> aDomain,
                                  ConstAlias<WeightImage> aWeightImage,
                                  ConstAlias<PowerSeparableMetric> aMetric,
                                  typename Parent::PeriodicitySpec const & aPeriodicitySpec)
      : PowerMap<TWeightImage,TPSeparableMetric,TImageContainer>(aDomain,
                                                                 aWeightImage,
                                                                 aMetric,
                                                                 aPeriodicitySpec)
    {}
 
    ~ReverseDistanceTransformation() {};
 
    // ------------------- Private functions ------------------------
  public:
 
    const Domain & domain() const
    {
      return Parent::domain();
    }
 
    ConstRange constRange() const
    {
      return ConstRange(*this);
    }
 
    Value operator()(const Point &aPoint) const
    {
      return this->myMetricPtr->powerDistance(aPoint,
                                              this->myImagePtr->operator()(aPoint),
                                              this->myWeightImagePtr->operator()( this->myImagePtr->operator()(aPoint)));
    }
 
    Vector getPowerVector(const Point &aPoint) const
    {
      return this->myImagePtr->operator()(aPoint);
    }
 
    const PowerSeparableMetric* metricPtr() const
    {
      return Parent::metricPtr();
    }
 
    void selfDisplay ( std::ostream & out ) const
    {
      out << "[ReverseDistanceTransformation] underlying PowerMap={";
      Parent::selfDisplay(out);
      out << "}";
    }
 
    // ------------------- protected methods ------------------------
  protected:
 
    ReverseDistanceTransformation();
 
 
    // ------------------- Private members ------------------------
  private:
 
  }; // end of class ReverseDistanceTransformation
 
 
// //                                                                           //
// ///////////////////////////////////////////////////////////////////////////////
 
  template <typename W,typename TSep>
  inline
  std::ostream&
  operator<< ( std::ostream & out,
               const ReverseDistanceTransformation<W,TSep> & object )
  {
    object.selfDisplay( out );
    return out;
  }
 
 
 
} // namespace DGtal
 
//                                                                           //
 
#endif // !defined ReverseDistanceTransformation_h
 
#undef ReverseDistanceTransformation_RECURSES
#endif // else defined(ReverseDistanceTransformation_RECURSES)