LambdaMST3DBy2D.h

 
#pragma once
 
#if defined(LAMBDAMST3DBy2D_RECURSES)
#error Recursive header files inclusion detected in LambdaMST3DBy2D.h
#else // defined(LAMBDAMST3DBy2D_RECURSES)
#define LAMBDAMST3DBy2D_RECURSES
 
#if !defined LAMBDAMST3DBy2D_h
#define LAMBDAMST3DBy2D_h
 
#include <algorithm>
#include <cmath>
#include <list>
#include <iterator>
#include <DGtal/base/Common.h>
#include <DGtal/helpers/StdDefs.h>
#include "DGtal/kernel/CSpace.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/kernel/BasicPointFunctors.h"
#include "DGtal/geometry/curves/estimation/LambdaMST2D.h"
#include "DGtal/geometry/curves/ArithmeticalDSSComputer.h"
#include "DGtal/geometry/curves/SaturatedSegmentation.h"
 
namespace DGtal {
  template < typename Iterator3D, typename Functor, typename LambdaFunctor, int CONNECTIVITY = 8 >
  class LambdaMST3DBy2DEstimator
  {
  public:
    // ----------------------- Types ------------------------------
  public:
    typedef PointVector < 3, double > RealVector3D;
    typedef PointVector < 3, int > Point3D;
    typedef PointVector < 2, int > Point2D;
    typedef PointVector < 2, double > RealVector2D;
    typedef std::vector < Point2D > TCurve2D;
    typedef ArithmeticalDSSComputer < typename TCurve2D::const_iterator, int, CONNECTIVITY > SegmentComputer2D;
    typedef SaturatedSegmentation < SegmentComputer2D > Segmentation2D;
    typedef typename Functor::MAIN_AXIS MAIN_AXIS;
 
    // ----------------------- Private types ------------------------------
  private:
    typedef LambdaMST2D < Segmentation2D, LambdaFunctor > TEstimator;
    typedef functors::Projector < SpaceND < 2, int > > Projector2d;
    
    // ----------------------- Standard services ------------------------------
  public:
    LambdaMST3DBy2DEstimator();
    
    void init ( Iterator3D itB, Iterator3D itE, MAIN_AXIS axis );
    
    bool isValid ( ) const;
    
    RealVector3D eval ( const Point3D & point );
 
    template <typename OutputIterator>
    OutputIterator eval ( Iterator3D itb, Iterator3D ite, OutputIterator result );
 
    // ------------------------- Internals ------------------------------------
  protected:
    
    RealVector2D Estimate2DTangent ( TCurve2D::const_iterator itb, TCurve2D::const_iterator ite, const Point2D & point );
    template < typename OutputIterator >
    OutputIterator Estimate2DTangent ( TCurve2D::const_iterator itb, TCurve2D::const_iterator ite, OutputIterator result );
 
    // ------------------------- Private Datas --------------------------------
  private:
    Iterator3D myBegin;
    Iterator3D myEnd;
    Functor myFunctor;
    MAIN_AXIS myAxis;
    TCurve2D tXY, tXZ, tYZ;
    Projector2d myProjXY, myProjXZ, myProjYZ;
  }; // end of class LambdaMST3DBy2DEstimator
  
  
  class TangentFromDSS3DBy2DFunctor
  {
  public:
    // ----------------------- Types ------------------------------
    typedef PointVector<3, double> Vector3D;
    typedef PointVector<2, double> Vector2D;
    enum MAIN_AXIS {X = 0, Y = 1, Z = 2};
    
    // ----------------------- Interface --------------------------------------
    Vector3D operator() ( MAIN_AXIS mainAxis, const Vector2D & v0, const Vector2D & v1 ) const;
  };
  
  //-------------------------------------------------------------------------------------------
  
  // Template class LambdaMST3D
  template < 
  typename Iterator3D, typename LambdaFunctor = functors::Lambda64Function, int CONNECTIVITY = 8 >
  class LambdaMST3DBy2D:
  public LambdaMST3DBy2DEstimator < Iterator3D, TangentFromDSS3DBy2DFunctor, LambdaFunctor, CONNECTIVITY >
  {
    typedef LambdaMST3DBy2DEstimator < Iterator3D, TangentFromDSS3DBy2DFunctor, LambdaFunctor, CONNECTIVITY > Super;
    
  public:
    LambdaMST3DBy2D() : Super() {}
  };
}// namespace DGtal
 
// Includes inline functions
#include "DGtal/geometry/curves/estimation/LambdaMST3DBy2D.ih"
 
//                                                                           //
 
#endif // !defined LAMBDAMST3DBy2D_h
 
#undef LAMBDAMST3DBy2D_RECURSES
#endif // else defined(LAMBDAMST3DBy2D_RECURSES)