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/images/CImage.h"
#include "DGtal/geometry/volumes/distance/SeparableMetricHelper.h"
#include "DGtal/kernel/NumberTraits.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"
#include "DGtal/kernel/sets/CDigitalSet.h"

namespace DGtal
{

  // template class ReverseDistanceTransformation
  template <typename Image, DGtal::uint32_t p, typename IntegerShort = DGtal::int8_t >
  class ReverseDistanceTransformation
  {

  public:
    
    //BOOST_CONCEPT_ASSERT(( CImage<Image> ));
    BOOST_CONCEPT_ASSERT(( CBoundedInteger<IntegerShort> ));
    

    typedef ImageContainerBySTLVector<  HyperRectDomain<typename Image::Domain::Space> , IntegerShort > OutputImage;
  
    typedef typename Image::Value Value;
    typedef typename Image::Point Point;
    typedef typename Image::Vector Vector;
    typedef typename Image::Dimension Dimension;
    typedef typename Image::Size Size;
    typedef typename Image::Integer Integer;
    typedef typename Image::Domain Domain;
    typedef typename Image::Domain::Space::Point::Coordinate Coordinate;
  
    typedef SeparableMetricHelper<  Point ,  Integer , p > SeparableMetric;
  

    ReverseDistanceTransformation(const IntegerShort defaultForeground = NumberTraits< IntegerShort >::ONE, 
          const IntegerShort defaultBackground = NumberTraits< IntegerShort >::ZERO);

    ~ReverseDistanceTransformation();

  public:

    OutputImage reconstruction(const Image & inputImage);

    template<typename DigitalSet>
    void reconstructionAsSet(DigitalSet &aSet, const Image &inputImage);
    
    
    

    // ------------------- Private functions ------------------------
  private:

    
    bool reconstructionInternal(const Image &aImage, Image &output, Image &swap);

    OutputImage castValues(const Image &input) const;
  
    void computeSteps(const Image & inputImage, 
          Image & output, 
          const Dimension dim)const;

    void computeSteps1D (const Image & input, 
       Image & output, 
       const Point &startingPoint, 
       const Size dim, 
       Integer s[], Integer t[]) const;


    // ------------------- Private members ------------------------
  private:
    SeparableMetric myMetric;
  
    Point myLowerBoundCopy;
    
    Point myUpperBoundCopy;

    Point myExtent;

    Vector myDisplacementVector;

    IntegerShort myForegroundValue;
  
    IntegerShort myBackgroundValue;

  }; // end of class ReverseDistanceTransformation

} // namespace DGtal


// Includes inline functions.
#include "DGtal/geometry/volumes/distance/ReverseDistanceTransformation.ih"

//                                                                           //

#endif // !defined ReverseDistanceTransformation_h

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