PlaneProbingLNeighborhood.h

 
#pragma once
 
#if defined(PlaneProbingLNeighborhood_RECURSES)
#error Recursive header files inclusion detected in PlaneProbingLNeighborhood.h
#else // defined(PlaneProbingLNeighborhood_RECURSES)
#define PlaneProbingLNeighborhood_RECURSES
 
#if !defined PlaneProbingLNeighborhood_h
#define PlaneProbingLNeighborhood_h
 
// Inclusions
#include <cassert>
#include <iostream>
#include <vector>
#include "DGtal/base/Common.h"
#include "DGtal/geometry/helpers/PlaneProbingEstimatorHelper.h"
#include "DGtal/geometry/surfaces/estimation/PlaneProbingRNeighborhood.h"
#include "DGtal/kernel/CPointPredicate.h"
 
namespace DGtal
{
 
  // template class PlaneProbingLNeighborhood
  template <typename TPredicate>
  class PlaneProbingLNeighborhood: public DGtal::PlaneProbingRNeighborhood<TPredicate>
  {
    BOOST_CONCEPT_ASSERT((DGtal::concepts::CPointPredicate<TPredicate>));
 
    // ----------------------- Public types ------------------------------
  public:
    using Predicate         = TPredicate;
    using Point             = typename Predicate::Point;
    using Vector            = Point;
    using Integer           = typename Point::Coordinate;
    using Triangle          = std::array<Vector, 3>;
 
    using HexagonState        = typename PlaneProbingNeighborhood<TPredicate>::HexagonState;
    using UpdateOperation     = typename PlaneProbingNeighborhood<TPredicate>::UpdateOperation;
 
    using Index                 = typename PlaneProbingNeighborhood<TPredicate>::Index;  
    using PointOnProbingRay     = typename PlaneProbingNeighborhood<TPredicate>::PointOnProbingRay;
    using GridPoint             = typename detail::GridPoint<Integer, Index>;
    using GridPointOnProbingRay = typename detail::GridPointOnProbingRay<Integer, Index>;
    
    // ----------------------- Internal type -------------------------------
  private: 
    struct ClosestGridPoint
    {
      ClosestGridPoint () = default;
 
      ClosestGridPoint (const GridPoint& aGridPoint,
                        const bool& aFirst, const bool& aSecond )
        : myGridPoint(aGridPoint), myPair(std::make_pair(aFirst,aSecond)) {}
 
      GridPoint myGridPoint; 
      std::pair<bool,bool> myPair; 
    };
    
    // ----------------------- Standard services ------------------------------
  public:
    PlaneProbingLNeighborhood() = delete;
 
    PlaneProbingLNeighborhood(Predicate const& aPredicate, Point const& aQ, Triangle const& aM);
 
    ~PlaneProbingLNeighborhood();
 
    PlaneProbingLNeighborhood ( const PlaneProbingLNeighborhood & other ) = delete;
 
    PlaneProbingLNeighborhood ( PlaneProbingLNeighborhood && other ) = delete;
 
    PlaneProbingLNeighborhood & operator= ( const PlaneProbingLNeighborhood & other ) = delete;
 
    PlaneProbingLNeighborhood & operator= ( PlaneProbingLNeighborhood && other ) = delete;
 
    // ----------------------- Plane-Probing services ------------------------------
  public:
    
    HexagonState hexagonState () override;
 
    UpdateOperation closestCandidate () override;
    
    // ------------------------- Protected Datas ------------------------------
  protected:
 
    std::vector<ClosestGridPoint> myGrids; 
    // ------------------------- Helpers to find a closest point --------------
  protected:
    
    ClosestGridPoint closestInGrid (const Index& aIdx) const;
 
    void candidatesInGrid (const GridPoint& y1, const GridPoint& y2,
                           std::back_insert_iterator<std::vector<GridPoint> > out) const;
 
    GridPointOnProbingRay closestOnBoundedRayLinearWithPredicate (GridPointOnProbingRay const& aRay, Integer const& aBound) const;
    
    GridPointOnProbingRay closestOnBoundedRayLogWithPredicate (GridPointOnProbingRay const& aRay, Integer const& aBound) const; 
   
    UpdateOperation getOperationFromGridPoint (GridPoint const& aClosest) const;
   
    void updateGrid (const Index& aIdx);
 
    Point direction (GridPoint const& aP) const;
 
    // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
  }; // end of class PlaneProbingLNeighborhood
 
 
  template <typename TPredicate>
  std::ostream&
  operator<< ( std::ostream & out, const PlaneProbingLNeighborhood<TPredicate> & object );
 
} // namespace DGtal
 
 
  // Includes inline functions.
#include "DGtal/geometry/surfaces/estimation/PlaneProbingLNeighborhood.ih"
 
//                                                                           //
 
#endif // !defined PlaneProbingLNeighborhood_h
 
#undef PlaneProbingLNeighborhood_RECURSES
#endif // else defined(PlaneProbingLNeighborhood_RECURSES)