VoronoiMap.h

#pragma once

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

#define VoronoiMap_RECURSES

#if !defined VoronoiMap_h

#define VoronoiMap_h

// Inclusions
#include <iostream>
#include <vector>
#include "DGtal/base/Common.h"
#include "DGtal/images/ImageContainerBySTLVector.h"
#include "DGtal/kernel/CPointPredicate.h"

#include "DGtal/geometry/volumes/distance/SeparableMetricHelper.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"

namespace DGtal
{

  // template class VoronoiMap
  template < typename TSpace,
             typename TPointPredicate,
             DGtal::uint32_t p>
  class VoronoiMap
  {

  public:
    BOOST_CONCEPT_ASSERT(( CSpace< TSpace > ));
    BOOST_CONCEPT_ASSERT(( CPointPredicate<TPointPredicate> ));

    BOOST_STATIC_ASSERT ((boost::is_same< typename TSpace::Point,
                                          typename TPointPredicate::Point >::value )); 
    
    typedef TSpace Space;

    typedef TPointPredicate PointPredicate;

    typedef HyperRectDomain<Space> Domain;
   
    typedef DGtal::int64_t IntegerLong;

    typedef typename Space::Vector Vector;
    typedef typename Space::Point Point;
    typedef typename Space::Dimension Dimension;
    typedef typename Space::Size Size;
    typedef typename Space::Point::Coordinate Abscissa;
 
    typedef SeparableMetricHelper<  Point ,  IntegerLong , p > SeparableMetric;
  
    typedef ImageContainerBySTLVector<  Domain,
                                        Point > OutputImage;

    VoronoiMap(const Domain & aDomain,
               const PointPredicate & predicate);

    ~VoronoiMap();

  public:

    OutputImage compute ( ) ;

    
    // ------------------- Private functions ------------------------
  private:    
    
    void computeOtherSteps(OutputImage & output,
                           const Dimension dim) const;
    void computeOtherStep1D (OutputImage & output, 
                             const Point &row, 
                             const Size dim,
                             std::vector<Point> &Sites) const;
    
    // ------------------- protected methods ------------------------
  protected:

    VoronoiMap();
   
    
    // ------------------- Private members ------------------------
  private:

    SeparableMetric myMetric;

    const Domain & myDomain;
    
    const PointPredicate  & myPointPredicate;
    
    Point myLowerBoundCopy;
    
    Point myUpperBoundCopy;
    
    Point myInfinity;


  }; // end of class VoronoiMap

} // namespace DGtal


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

//                                                                           //

#endif // !defined VoronoiMap_h

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