ChordNaivePlaneComputer.h

 
#pragma once
 
#if defined(ChordNaivePlaneComputer_RECURSES)
#error Recursive header files inclusion detected in ChordNaivePlaneComputer.h
#else // defined(ChordNaivePlaneComputer_RECURSES)
#define ChordNaivePlaneComputer_RECURSES
 
#if !defined ChordNaivePlaneComputer_h
#define ChordNaivePlaneComputer_h
 
// Inclusions
#include <iostream>
#include <set>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/CSignedNumber.h"
#include "DGtal/kernel/CSpace.h"
#include "DGtal/geometry/surfaces/ParallelStrip.h"
 
namespace DGtal
{
 
  // template class ChordNaivePlaneComputer
  template < typename TSpace,
             typename TInputPoint,
             typename TInternalScalar >
  class ChordNaivePlaneComputer
  {
 
    BOOST_CONCEPT_ASSERT(( concepts::CSpace< TSpace > ));
    BOOST_CONCEPT_ASSERT(( concepts::CSignedNumber< TInternalScalar > ));
    BOOST_STATIC_ASSERT(( TSpace::dimension == 3 ));
    BOOST_STATIC_ASSERT(( TInputPoint::dimension == 3 ));
 
    // ----------------------- public types ------------------------------
  public:
    typedef TSpace Space;
    typedef TInputPoint InputPoint;
    typedef TInternalScalar InternalScalar;
    typedef InputPoint InputVector;
    typedef typename InputVector::Component Component;
    typedef typename InputPoint::Coordinate Coordinate;
    typedef PointVector<3,InternalScalar> InternalVector;
 
    typedef std::set< InputPoint > InputPointSet;
    typedef typename InputPointSet::size_type Size;
    typedef typename InputPointSet::const_iterator ConstIterator;
    typedef typename InputPointSet::iterator Iterator;
    typedef typename Space::Point Point;
    typedef ParallelStrip<Space,true,true> Primitive;
 
    // ----------------------- std public types ------------------------------
  public:
    typedef typename InputPointSet::const_iterator const_iterator;
    typedef typename InputPointSet::const_pointer const_pointer;
    typedef typename InputPointSet::const_reference const_reference;
    typedef typename InputPointSet::value_type value_type;
    typedef typename InputPointSet::difference_type difference_type;
    typedef typename InputPointSet::size_type size_type;
 
    // ----------------------- internal types ------------------------------
  private:
    struct State {
      InternalScalar height;   
      InputPoint A,B,C;        
      InternalVector N;        
      InputPoint ptMax;        
      InputPoint ptMin;        
      InternalScalar max;      
      InternalScalar min;      
      unsigned int nbValid;    
    };
 
    // ----------------------- Standard services ------------------------------
  public:
 
    ~ChordNaivePlaneComputer();
 
    ChordNaivePlaneComputer();
 
    ChordNaivePlaneComputer ( const ChordNaivePlaneComputer & other );
 
    ChordNaivePlaneComputer & operator= ( const ChordNaivePlaneComputer & other );
 
    void clear();
 
    void init( Dimension axis,
               InternalScalar widthNumerator = NumberTraits< InternalScalar >::ONE,
               InternalScalar widthDenominator = NumberTraits< InternalScalar >::ONE );
 
    template <typename TInputIterator>
    static
    std::pair<InternalScalar, InternalScalar>
    computeAxisWidth( Dimension axis, TInputIterator itB, TInputIterator itE );
 
 
    //-------------------- model of ForwardContainer -----------------------------
  public:
 
    Size size() const;
 
    bool empty() const;
 
    ConstIterator begin() const;
 
    ConstIterator end() const;
 
    Size max_size() const;
 
    Size maxSize() const;
 
 
    //-------------------- model of concepts::CPointPredicate -----------------------------
  public:
 
    bool operator()( const Point & p ) const;
 
    //-------------------- model of CIncrementalPrimitiveComputer -----------------------------
  public:
 
    bool extendAsIs( const InputPoint & p );
 
    bool extend( const InputPoint & p );
 
    bool isExtendable( const InputPoint & p ) const;
 
    //-------------------- model of CAdditivePrimitiveComputer -----------------------------
  public:
 
    template <typename TInputIterator>
    bool extend( TInputIterator itB, TInputIterator itE );
 
    template <typename TInputIterator>
    bool isExtendable( TInputIterator itB, TInputIterator itE ) const;
 
    template <typename TInputIterator>
    bool
    satisfies( TInputIterator itB, TInputIterator itE ) const;
 
    template <typename TInputIterator>
    std::pair<InternalScalar, InternalScalar>
    axisWidth( TInputIterator itB, TInputIterator itE ) const;
 
    //-------------------- Primitive services -----------------------------
  public:
 
    Primitive primitive() const;
 
    //-------------------- Parameters services -----------------------------
  public:
 
    template <typename Vector3D>
    void getNormal( Vector3D & normal ) const;
 
    const InternalVector & exactNormal() const;
 
    template <typename Vector3D>
    void getUnitNormal( Vector3D & normal ) const;
 
    void getBounds( double & min, double & max ) const;
 
    const InputPoint & minimalPoint() const;
 
    const InputPoint & maximalPoint() const;
 
    // ----------------------- Utilities --------------------------------------
  public:
 
    template <typename TVector1, typename TVector2>
    static InternalScalar internalDot( const TVector1 & u, const TVector2 & v );
 
    template <typename TVector1, typename TVector2>
    static void internalCross( InternalVector & n, const TVector1 & u, const TVector2 & v );
 
    // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
    // ------------------------- Private Datas --------------------------------
  private:
    Dimension z;               
    Dimension x,y;             
    InternalScalar myWidth0;   
    InternalScalar myWidth1;   
    mutable InputPointSet myPointSet;
    State myState;             
    mutable State _state;      
    mutable InternalScalar _d; 
    // ------------------------- Hidden services ------------------------------
  protected:
 
 
    // ------------------------- Internals ------------------------------------
  private:
 
    template <typename TInputIterator>
    void computeMinMax( State & state, TInputIterator itB, TInputIterator itE ) const;
 
    template <typename TInputIterator>
    bool updateMinMax( State & state, TInputIterator itB, TInputIterator itE ) const;
 
    bool checkPlaneWidth( const State & state ) const;
 
    bool setUp1( const InputPoint & p1 );
 
    void setUpNormal( State & state ) const;
 
    void setUpNormal1( State & state ) const;
 
    void setUpNormal2( State & state ) const;
 
    void setUpNormal3( State & state ) const;
 
    int signDelta( const InputPoint & A, const InputPoint & B, const InputPoint & C ) const;
 
    int signDelta( const InputPoint & A, const InputPoint & C ) const;
 
    bool newCurrentTriangle( State & state, const InputPoint & M ) const;
 
    template <typename TInputIterator>
    unsigned int findTriangle( State & state, TInputIterator itB, TInputIterator itE ) const;
 
    template <typename TInputIterator>
    unsigned int findTriangle1( State & state, TInputIterator itB, TInputIterator itE ) const;
    template <typename TInputIterator>
    unsigned int findTriangle2( State & state, TInputIterator itB, TInputIterator itE ) const;
    template <typename TInputIterator1, typename TInputIterator2>
    unsigned int findMixedTriangle( State & state,
                                    TInputIterator1 itB1, TInputIterator1 itE1,
                                    TInputIterator2 itB2, TInputIterator2 itE2 ) const;
 
    bool alignedAlongAxis( const InputPoint & p1, const InputPoint & p2 ) const;
 
    void computeHeight( State & state ) const;
 
    void computeNormal( State & state ) const;
 
    template <typename TInputIterator>
    bool
    satisfies( State & state, TInputIterator itB, TInputIterator itE ) const;
 
 
    template <typename TInputIterator>
    std::pair<InternalScalar, InternalScalar>
    axisWidth( State & state, TInputIterator itB, TInputIterator itE ) const;
 
 
    bool checkWidth( const State & state ) const;
 
    void selfDisplay ( std::ostream & out, const State & state ) const;
 
 
  }; // end of class ChordNaivePlaneComputer
 
 
  template <typename TSpace, typename TInputPoint, typename TInternalScalar>
  std::ostream&
  operator<< ( std::ostream & out,
               const ChordNaivePlaneComputer<TSpace, TInputPoint, TInternalScalar> & object );
 
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/geometry/surfaces/ChordNaivePlaneComputer.ih"
 
//                                                                           //
 
#endif // !defined ChordNaivePlaneComputer_h
 
#undef ChordNaivePlaneComputer_RECURSES
#endif // else defined(ChordNaivePlaneComputer_RECURSES)