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() {};
        
    // ------------------- 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)