ChordGenericStandardPlaneComputer.h

 
#pragma once
 
#if defined(ChordGenericStandardPlaneComputer_RECURSES)
#error Recursive header files inclusion detected in ChordGenericStandardPlaneComputer.h
#else // defined(ChordGenericStandardPlaneComputer_RECURSES)
#define ChordGenericStandardPlaneComputer_RECURSES
 
#if !defined ChordGenericStandardPlaneComputer_h
#define ChordGenericStandardPlaneComputer_h
 
// Inclusions
#include <iostream>
#include <set>
#include "DGtal/base/Common.h"
#include "DGtal/base/IteratorAdapter.h"
#include "DGtal/kernel/CInteger.h"
#include "DGtal/kernel/CSpace.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/arithmetic/IntegerComputer.h"
#include "DGtal/geometry/surfaces/ChordNaivePlaneComputer.h"
 
namespace DGtal
{
 
  // template class ChordGenericStandardPlaneComputer
  template < typename TSpace,
             typename TInputPoint,
             typename TInternalScalar >
  class ChordGenericStandardPlaneComputer
  {
    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 ChordNaivePlaneComputer< Space, InputPoint, InternalScalar > ChordComputer;
    typedef typename ChordComputer::Component Component;
    typedef typename ChordComputer::Coordinate Coordinate;
    typedef typename ChordComputer::InternalVector InternalVector;
    typedef typename ChordComputer::Point Point;
    typedef typename ChordComputer::Primitive Primitive;
    typedef typename ChordComputer::InputPointSet InputPointSet;
    typedef typename ChordComputer::Size Size;
    typedef typename ChordComputer::ConstIterator InputPointSetConstIterator;
    typedef typename ChordComputer::Iterator InputPointSetIterator;
 
    // ----------------------- internal types ------------------------------
  private:
    typedef std::vector<Dimension>::iterator OrthantIterator;
    typedef std::vector<Dimension>::const_iterator OrthantConstIterator;
 
 
    // --------- adapter types for standard to naive transformation --------
  public:
 
    struct Transform {
      typedef InputPoint Value;
      typedef InputPoint value_type;
      bool _posX; 
      bool _posY; 
      inline Transform() {}
      inline Transform( bool posX, bool posY ) : _posX( posX ), _posY( posY ) {}
      inline InputPoint operator()( const InputPoint & p ) const
      { return InputPoint( p[ 0 ] + ( _posX ? p[ 2 ] : -p[ 2 ] ),
                           p[ 1 ] + ( _posY ? p[ 2 ] : -p[ 2 ] ),
                           p[ 2 ] ); }
      inline bool isPlusZOnX() const { return _posX; }
      inline bool isPlusZOnY() const { return _posY; }
      inline InputPoint inverse( const InputPoint & p ) const
      { return InputPoint( p[ 0 ] + ( _posX ? -p[ 2 ] : p[ 2 ] ),
                           p[ 1 ] + ( _posY ? -p[ 2 ] : p[ 2 ] ),
                           p[ 2 ] ); }
      template <typename Vector3D>
      void transformBack( Vector3D & v )
      {
        v[ 2 ] += ( _posX ? v[ 0 ] : -v[ 0 ] );
        v[ 2 ] += ( _posY ? v[ 1 ] : -v[ 1 ] );
      }
      template <typename Vector3D>
      void transform( Vector3D & v )
      {
        v[ 2 ] += ( _posX ? -v[ 0 ] : v[ 0 ] );
        v[ 2 ] += ( _posY ? -v[ 1 ] : v[ 1 ] );
      }
    };
 
    typedef boost::transform_iterator<Transform,InputPointSetConstIterator,
                                      InputPoint,InputPoint>
    ConstIterator;
 
    // ----------------------- std public types ------------------------------
  public:
    typedef ConstIterator const_iterator;
    typedef const InputPoint* const_pointer;
    typedef const InputPoint& const_reference;
    typedef InputPoint value_type;
    typedef typename InputPointSet::difference_type difference_type;
    typedef typename InputPointSet::size_type size_type;
 
    // ----------------------- Standard services ------------------------------
  public:
 
    ~ChordGenericStandardPlaneComputer();
 
    ChordGenericStandardPlaneComputer();
 
    ChordGenericStandardPlaneComputer ( const ChordGenericStandardPlaneComputer & other );
 
    ChordGenericStandardPlaneComputer & operator= ( const ChordGenericStandardPlaneComputer & other );
 
    Dimension active() const;
 
    void clear();
 
    void init( InternalScalar widthNumerator = NumberTraits< InternalScalar >::ONE,
               InternalScalar widthDenominator = NumberTraits< InternalScalar >::ONE );
 
    //-------------------- 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 it, TInputIterator itE );
 
    template <typename TInputIterator>
    bool isExtendable( TInputIterator it, TInputIterator itE ) const;
 
    //-------------------- Primitive services -----------------------------
  public:
 
    Primitive primitive() const;
 
    //-------------------- Parameters services -----------------------------
  public:
 
    template <typename Vector3D>
    void getNormal( Vector3D & normal ) const;
 
    template <typename Vector3D>
    void getUnitNormal( Vector3D & normal ) const;
 
    void getBounds( double & min, double & max ) const;
 
    InputPoint minimalPoint() const;
 
    InputPoint maximalPoint() const;
 
    void getCharacteristics( InternalVector & n,
                             InternalScalar & imin, InternalScalar & imax,
                             InputPoint & p_min, InputPoint & p_max ) const;
 
    // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
    // ------------------------- Private Datas --------------------------------
  private:
    std::vector<Dimension> myOrthants; 
    ChordComputer myComputers[ 4 ];    
    static Transform myTransforms[ 4 ];
    mutable std::vector<Dimension> _orthantsToErase; 
    // ------------------------- Hidden services ------------------------------
  protected:
    static Transform t( Dimension orthant );
    static Transform invT( Dimension orthant );
 
    // ------------------------- Internals ------------------------------------
  private:
  }; // end of class ChordGenericStandardPlaneComputer
 
 
  template <typename TSpace, typename TInputPoint, typename TInternalScalar>
  std::ostream&
  operator<< ( std::ostream & out,
               const ChordGenericStandardPlaneComputer<TSpace, TInputPoint, TInternalScalar> & object );
 
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/geometry/surfaces/ChordGenericStandardPlaneComputer.ih"
 
//                                                                           //
 
#endif // !defined ChordGenericStandardPlaneComputer_h
 
#undef ChordGenericStandardPlaneComputer_RECURSES
#endif // else defined(ChordGenericStandardPlaneComputer_RECURSES)