EuclideanShapesDecorator.h

 
#pragma once
 
#if defined(EuclideanShapesDecorator_RECURSES)
#error Recursive header files inclusion detected in EuclideanShapesDecorator.h
#else // defined(EuclideanShapesDecorator_RECURSES)
#define EuclideanShapesDecorator_RECURSES
 
#if !defined EuclideanShapesDecorator_h
#define EuclideanShapesDecorator_h
 
// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/base/ConstAlias.h"
 
#include "DGtal/shapes/CEuclideanBoundedShape.h"
#include "DGtal/shapes/CEuclideanOrientedShape.h"
 
namespace DGtal
{
 
  // template class EuclideanShapesCSG
  template <typename ShapeA, typename ShapeB>
  class EuclideanShapesCSG
  {
  protected:
    enum e_operator
    {
      e_plus,
      e_intersection,
      e_minus
    };
 
  public:
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeA > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeA > ));
 
    typedef typename ShapeA::Space Space;
    typedef typename ShapeA::RealPoint RealPoint;
 
    EuclideanShapesCSG( )
      : bIsValid(false)
    {}
 
    EuclideanShapesCSG ( const EuclideanShapesCSG & other )
      : myShapeA(other.myShapeA), v_shapes(other.v_shapes),
        myLowerBound(other.myLowerBound), myUpperBound(other.myUpperBound),
        bIsValid(other.bIsValid)
    {}
 
    EuclideanShapesCSG( ConstAlias<ShapeA> a )
      : myShapeA( a )
    {
      myLowerBound = myShapeA->getLowerBound();
      myUpperBound = myShapeA->getUpperBound();
 
      bIsValid = true;
    }
 
    EuclideanShapesCSG & operator= ( const EuclideanShapesCSG & other )
    {
      myShapeA = other.myShapeA;
      v_shapes = other.v_shapes;
 
      myLowerBound = other.myLowerBound;
      myUpperBound = other.myUpperBound;
 
      bIsValid = other.bIsValid;
 
      return *this;
    }
 
    void setParams( ConstAlias<ShapeA> a )
    {
      myShapeA = a;
 
      myLowerBound = myShapeA->getLowerBound();
      myUpperBound = myShapeA->getUpperBound();
 
      bIsValid = true;
    }
 
    void plus( ConstAlias<ShapeB> b )
    {
      BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeB > ));
      BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeB > ));
 
      FATAL_ERROR_MSG( isValid(), "Operation invalid. Maybe you don't set a ShapeA object." );
 
      std::pair<e_operator, CountedConstPtrOrConstPtr< ShapeB > > shape( e_plus, b );
 
      for(Dimension i =0; i < Space::dimension; ++i)
      {
        myLowerBound[i] = std::min(myLowerBound[i], b->getLowerBound()[i]);
        myUpperBound[i] = std::max(myUpperBound[i], b->getUpperBound()[i]);
      }
 
      v_shapes.push_back(shape); 
    }
 
    void intersection( ConstAlias<ShapeB> b )
    {
      BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeB > ));
      BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeB > ));
 
      FATAL_ERROR_MSG( isValid(), "Operation invalid. Maybe you don't set a ShapeA object." );
 
      std::pair<e_operator, CountedConstPtrOrConstPtr< ShapeB > > shape( e_intersection, b );
 
      for(Dimension i=0; i < Space::dimension; ++i)
      {
        myLowerBound[i] = std::max(myLowerBound[i], b->getLowerBound()[i]);
        myUpperBound[i] = std::min(myUpperBound[i], b->getUpperBound()[i]);
      }
 
      v_shapes.push_back(shape); 
    }
 
    void minus( ConstAlias<ShapeB> b )
    {
      BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeB > ));
      BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeB > ));
 
      FATAL_ERROR_MSG( isValid(), "Operation invalid. Maybe you don't set a ShapeA object." );
 
      std::pair<e_operator, CountedConstPtrOrConstPtr< ShapeB > > shape( e_minus, b );
 
      v_shapes.push_back(shape);
    }
 
    RealPoint getLowerBound() const
    {
      FATAL_ERROR_MSG( isValid(), "Operation invalid. Maybe you don't set a ShapeA object." );
 
      return myLowerBound;
    }
 
    RealPoint getUpperBound() const
    {
      FATAL_ERROR_MSG( isValid(), "Operation invalid. Maybe you don't set a ShapeA object." );
 
      return myUpperBound;
    }
 
    Orientation orientation( const RealPoint & p ) const
    {
      FATAL_ERROR_MSG( isValid(), "Operation invalid. Maybe you don't set a ShapeA object." );
 
      Orientation orient = myShapeA->orientation( p );
 
      for(unsigned int i = 0; i < v_shapes.size(); ++i)
      {
        if( v_shapes[i].first == e_minus )
        {
          if (( v_shapes[i].second->orientation( p ) == INSIDE ) || ( v_shapes[i].second->orientation( p ) == ON ))
          {
            orient = OUTSIDE;
          }
        }
        else if( v_shapes[i].first == e_intersection )
        {
          if (( orient == ON ) && ( v_shapes[i].second->orientation( p ) != OUTSIDE ))
          {
            orient = ON;
          }
          else if (( v_shapes[i].second->orientation( p ) == ON ) && ( orient != OUTSIDE ))
          {
            orient = ON;
          }
          else if (( orient == INSIDE ) && ( v_shapes[i].second->orientation( p ) == INSIDE ))
          {
            orient = INSIDE;
          }
          else
          {
            orient = OUTSIDE;
          }
        }
        else 
        {
          if (( orient == INSIDE ) || ( v_shapes[i].second->orientation( p ) == INSIDE ))
          {
            orient = INSIDE;
          }
          else if (( orient == ON ) || ( v_shapes[i].second->orientation( p ) == ON ))
          {
            orient = ON;
          }
          else
          {
            orient = OUTSIDE;
          }
        }
      }
 
      return orient;
    }
 
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const
    {
      return bIsValid;
    }
 
    // ------------------------- Internals ------------------------------------
  private:
 
    CountedConstPtrOrConstPtr< ShapeA > myShapeA;
 
    std::vector< std::pair<e_operator,  CountedConstPtrOrConstPtr< ShapeB > > > v_shapes;
 
    RealPoint myLowerBound;
 
    RealPoint myUpperBound;
 
    bool bIsValid;
 
  };
 
namespace deprecated
{
// template class EuclideanShapesDecorator
  template <typename ShapeA, typename ShapeB>
  class EuclideanShapesUnion
  {
    // ----------------------- Standard services ------------------------------
  public:
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeA > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeA > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeB > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeB > ));
 
    typedef typename ShapeA::Space Space;
    typedef typename ShapeA::RealPoint RealPoint;
 
    EuclideanShapesUnion( ConstAlias< ShapeA > a, ConstAlias< ShapeB > b )
      : myShapeA( a ),
        myShapeB( b )
    {
      RealPoint shapeALowerBoundary = myShapeA.getLowerBound();
      RealPoint shapeBLowerBoundary = myShapeB.getLowerBound();
      RealPoint shapeAUpperBoundary = myShapeA.getUpperBound();
      RealPoint shapeBUpperBoundary = myShapeB.getUpperBound();
      for ( unsigned int i = 0; i < myLowerBound.size(); ++i )
      {
        myLowerBound[ i ] = std::min( shapeALowerBoundary[ i ], shapeBLowerBoundary[ i ] );
        myUpperBound[ i ] = std::max( shapeAUpperBoundary[ i ], shapeBUpperBoundary[ i ] );
      }
    }
 
    RealPoint getLowerBound() const
    {
      return myLowerBound;
    }
 
    RealPoint getUpperBound() const
    {
      return myUpperBound;
    }
 
    Orientation orientation( const RealPoint & p ) const
    {
      if (( myShapeA.orientation( p ) == INSIDE ) || ( myShapeB.orientation( p ) == INSIDE ))
        {
            return INSIDE;
        }
      else if (( myShapeA.orientation( p ) == ON ) || ( myShapeB.orientation( p ) == ON ))
        {
            return ON;
        }
        return OUTSIDE;
    }
 
    ~EuclideanShapesUnion(){}
 
    // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
    // ------------------------- Hidden services ------------------------------
  protected:
 
    EuclideanShapesUnion();
 
  private:
 
    EuclideanShapesUnion ( const EuclideanShapesUnion & other );
 
    EuclideanShapesUnion & operator= ( const EuclideanShapesUnion & other );
 
    // ------------------------- Internals ------------------------------------
  private:
    const ShapeA & myShapeA;
    const ShapeB & myShapeB;
 
    RealPoint myLowerBound;
    RealPoint myUpperBound;
 
  }; // end of class EuclideanShapesUnion
 
  // template class EuclideanShapesIntersection
  template <typename ShapeA, typename ShapeB>
  class EuclideanShapesIntersection
  {
    // ----------------------- Standard services ------------------------------
  public:
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeA > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeA > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeB > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeB > ));
 
    typedef typename ShapeA::Space Space;
    typedef typename ShapeA::RealPoint RealPoint;
 
    EuclideanShapesIntersection( ConstAlias< ShapeA > a, ConstAlias< ShapeB > b )
      : myShapeA( a ),
        myShapeB( b )
    {
      RealPoint shapeALowerBoundary = myShapeA.getLowerBound();
      RealPoint shapeBLowerBoundary = myShapeB.getLowerBound();
      RealPoint shapeAUpperBoundary = myShapeA.getUpperBound();
      RealPoint shapeBUpperBoundary = myShapeB.getUpperBound();
      for ( unsigned int i = 0; i < myLowerBound.size(); ++i )
      {
        myLowerBound[ i ] = std::min( shapeALowerBoundary[ i ], shapeBLowerBoundary[ i ] );
        myUpperBound[ i ] = std::max( shapeAUpperBoundary[ i ], shapeBUpperBoundary[ i ] );
      }
    }
 
 
    RealPoint getLowerBound() const
    {
      return myLowerBound;
    }
 
    RealPoint getUpperBound() const
    {
      return myUpperBound;
    }
 
    Orientation orientation( const RealPoint & p ) const
    {
      if (( myShapeA.orientation( p ) == ON ) && ( myShapeB.orientation( p ) != OUTSIDE ))
      {
        return ON;
      }
      else if (( myShapeB.orientation( p ) == ON ) && ( myShapeA.orientation( p ) != OUTSIDE ))
      {
        return ON;
      }
      else if (( myShapeA.orientation( p ) == INSIDE ) && ( myShapeB.orientation( p ) == INSIDE ))
      {
        return INSIDE;
      }
 
      return OUTSIDE;
    }
 
 
    ~EuclideanShapesIntersection(){}
 
    // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
    // ------------------------- Hidden services ------------------------------
  protected:
 
    EuclideanShapesIntersection();
 
  private:
 
    EuclideanShapesIntersection ( const EuclideanShapesIntersection & other );
 
    EuclideanShapesIntersection & operator= ( const EuclideanShapesIntersection & other );
 
    // ------------------------- Internals ------------------------------------
  private:
    const ShapeA & myShapeA;
    const ShapeB & myShapeB;
 
    RealPoint myLowerBound;
    RealPoint myUpperBound;
 
  }; // end of class EuclideanShapesIntersection
 
  // template class EuclideanShapesMinus
  template <typename ShapeA, typename ShapeB>
  class EuclideanShapesMinus
  {
    // ----------------------- Standard services ------------------------------
  public:
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeA > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeA > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanBoundedShape< ShapeB > ));
    BOOST_CONCEPT_ASSERT (( concepts::CEuclideanOrientedShape< ShapeB > ));
 
    typedef typename ShapeA::Space Space;
    typedef typename ShapeA::RealPoint RealPoint;
 
    EuclideanShapesMinus( ConstAlias< ShapeA > a, ConstAlias< ShapeB > b )
      : myShapeA( a ),
        myShapeB( b )
    {
      RealPoint shapeALowerBoundary = myShapeA.getLowerBound();
      RealPoint shapeBLowerBoundary = myShapeB.getLowerBound();
      RealPoint shapeAUpperBoundary = myShapeA.getUpperBound();
      RealPoint shapeBUpperBoundary = myShapeB.getUpperBound();
      for ( unsigned int i = 0; i < myLowerBound.size(); ++i )
      {
        myLowerBound[ i ] = std::min( shapeALowerBoundary[ i ], shapeBLowerBoundary[ i ] );
        myUpperBound[ i ] = std::max( shapeAUpperBoundary[ i ], shapeBUpperBoundary[ i ] );
      }
    }
 
    RealPoint getLowerBound() const
    {
      return myLowerBound;
    }
 
    RealPoint getUpperBound() const
    {
      return myUpperBound;
    }
 
    Orientation orientation( const RealPoint & p ) const
    {
      if (( myShapeB.orientation( p ) == INSIDE ) || ( myShapeB.orientation( p ) == ON ))
      {
        return OUTSIDE;
      }
      return myShapeA.orientation( p );
    }
 
 
    ~EuclideanShapesMinus(){}
 
    // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
    // ------------------------- Hidden services ------------------------------
  protected:
 
    EuclideanShapesMinus();
 
  private:
 
    EuclideanShapesMinus ( const EuclideanShapesMinus & other );
 
    EuclideanShapesMinus & operator= ( const EuclideanShapesMinus & other );
 
    // ------------------------- Internals ------------------------------------
  private:
    const ShapeA & myShapeA;
    const ShapeB & myShapeB;
 
    RealPoint myLowerBound;
    RealPoint myUpperBound;
 
  }; // end of class EuclideanShapesMinus
 
 
 }
 
 
  template <typename ShapeA, typename ShapeB>
  std::ostream&
  operator<< ( std::ostream & out, const deprecated::EuclideanShapesUnion<ShapeA, ShapeB> & object );
 
  template <typename ShapeA, typename ShapeB>
  std::ostream&
  operator<< ( std::ostream & out, const deprecated::EuclideanShapesIntersection<ShapeA, ShapeB> & object );
 
  template <typename ShapeA, typename ShapeB>
  std::ostream&
  operator<< ( std::ostream & out, const deprecated::EuclideanShapesMinus<ShapeA, ShapeB> & object );
 
} // namespace DGtal
 
 
//                                                                           //
 
#endif // !defined EuclideanShapesDecorator_h
 
#undef EuclideanShapesDecorator_RECURSES
#endif // else defined(EuclideanShapesDecorator_RECURSES)