COBAGenericStandardPlaneComputer.h

 
#pragma once
 
#if defined(COBAGenericStandardPlaneComputer_RECURSES)
#error Recursive header files inclusion detected in COBAGenericStandardPlaneComputer.h
#else // defined(COBAGenericStandardPlaneComputer_RECURSES)
#define COBAGenericStandardPlaneComputer_RECURSES
 
#if !defined COBAGenericStandardPlaneComputer_h
#define COBAGenericStandardPlaneComputer_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/COBANaivePlaneComputer.h"
 
namespace DGtal
{
 
  // template class COBAGenericStandardPlaneComputer
  template < typename TSpace,
             typename TInternalInteger >
  class COBAGenericStandardPlaneComputer
  {
 
    BOOST_CONCEPT_ASSERT(( concepts::CSpace< TSpace > ));
    BOOST_CONCEPT_ASSERT(( concepts::CInteger< TInternalInteger > ));
    BOOST_STATIC_ASSERT(( TSpace::dimension == 3 ));
 
    // ----------------------- public types ------------------------------
  public:
    typedef TSpace Space;
    typedef typename Space::Point Point;
    typedef std::set< Point > PointSet;
    typedef typename PointSet::size_type Size;
    typedef typename PointSet::const_iterator PointSetConstIterator;
    typedef typename PointSet::iterator PointSetIterator;
    typedef TInternalInteger InternalInteger;
    typedef IntegerComputer< InternalInteger > MyIntegerComputer;
    typedef COBANaivePlaneComputer< Space, InternalInteger > COBAComputer;
    typedef typename COBAComputer::Primitive Primitive;
    typedef typename COBAComputer::IntegerVector3 IntegerVector3;
 
    // ----------------------- 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 Point Value;
      typedef Point value_type;
      bool _posX; 
      bool _posY; 
      inline Transform() {}
      inline Transform( bool posX, bool posY ) : _posX( posX ), _posY( posY ) {}
      inline Point operator()( const Point & p ) const
      {
        if ( _posX ) {
          if ( _posY ) return Point( p[ 0 ] + p[ 2 ], p[ 1 ] + p[ 2 ], p[ 2 ] );
          else         return Point( p[ 0 ] + p[ 2 ], p[ 1 ] - p[ 2 ], p[ 2 ] );
        } else {
          if ( _posY ) return Point( p[ 0 ] - p[ 2 ], p[ 1 ] + p[ 2 ], p[ 2 ] );
          else         return Point( p[ 0 ] - p[ 2 ], p[ 1 ] - p[ 2 ], p[ 2 ] );
        }
      }
      inline bool isPlusZOnX() const { return _posX; }
      inline bool isPlusZOnY() const { return _posY; }
      inline Point inverse( const Point & p ) const
      {
        if ( _posX ) {
          if ( _posY ) return Point( p[ 0 ] - p[ 2 ], p[ 1 ] - p[ 2 ], p[ 2 ] );
          else         return Point( p[ 0 ] - p[ 2 ], p[ 1 ] + p[ 2 ], p[ 2 ] );
        } else {
          if ( _posY ) return Point( p[ 0 ] + p[ 2 ], p[ 1 ] - p[ 2 ], p[ 2 ] );
          else         return Point( p[ 0 ] + p[ 2 ], p[ 1 ] + p[ 2 ], p[ 2 ] );
        }
      }
      template <typename Vector3D>
      void transformBack( Vector3D & v )
      {
        if ( _posX ) v[ 2 ] += v[ 0 ];
        else         v[ 2 ] -= v[ 0 ];
        if ( _posY ) v[ 2 ] += v[ 1 ];
        else         v[ 2 ] -= v[ 1 ];
      }
      template <typename Vector3D>
      void transform( Vector3D & v )
      {
        if ( _posX ) v[ 2 ] -= v[ 0 ];
        else         v[ 2 ] += v[ 0 ];
        if ( _posY ) v[ 2 ] -= v[ 1 ];
        else         v[ 2 ] += v[ 1 ];
      }
    };
 
    typedef boost::transform_iterator<Transform,PointSetConstIterator,Point,Point>
    ConstIterator;
 
    // ----------------------- std public types ------------------------------
  public:
    typedef ConstIterator const_iterator;
    typedef const Point* const_pointer;
    typedef const Point& const_reference;
    typedef Point value_type;
    typedef typename PointSet::difference_type difference_type;
    typedef typename PointSet::size_type size_type;
 
    // ----------------------- Standard services ------------------------------
  public:
 
    ~COBAGenericStandardPlaneComputer();
 
    COBAGenericStandardPlaneComputer();
 
    COBAGenericStandardPlaneComputer ( const COBAGenericStandardPlaneComputer & other );
 
    COBAGenericStandardPlaneComputer & operator= ( const COBAGenericStandardPlaneComputer & other );
 
    MyIntegerComputer & ic() const;
 
    Dimension active() const;
 
    void clear();
 
    void init( InternalInteger diameter,
               InternalInteger widthNumerator = NumberTraits< InternalInteger >::ONE,
               InternalInteger widthDenominator = NumberTraits< InternalInteger >::ONE );
 
    Size complexity() const;
 
    //-------------------- 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;
 
    bool extendAsIs( const Point & p );
 
    //-------------------- model of CIncrementalPrimitiveComputer -----------------------------
  public:
 
    bool extend( const Point & p );
 
    bool isExtendable( const Point & 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;
 
    Point minimalPoint() const;
 
    Point maximalPoint() const;
 
    void getCharacteristics( IntegerVector3 & n,
                             InternalInteger & imin, InternalInteger & imax,
                             Point & p_min, Point & p_max ) const;
 
    // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
    // ------------------------- Private Datas --------------------------------
  private:
    std::vector<Dimension> myOrthants; 
    COBAComputer 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 COBAGenericStandardPlaneComputer
 
 
  template <typename TSpace, typename TInternalInteger>
  std::ostream&
  operator<< ( std::ostream & out, const COBAGenericStandardPlaneComputer<TSpace, TInternalInteger> & object );
 
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/geometry/surfaces/COBAGenericStandardPlaneComputer.ih"
 
//                                                                           //
 
#endif // !defined COBAGenericStandardPlaneComputer_h
 
#undef COBAGenericStandardPlaneComputer_RECURSES
#endif // else defined(COBAGenericStandardPlaneComputer_RECURSES)