testDigitalShapesDecorator.cpp

#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/shapes/GaussDigitizer.h"

#include "DGtal/shapes/DigitalShapesDecorator.h"

#include "DGtal/shapes/parametric/Ball2D.h"


using namespace std;
using namespace DGtal;

// Functions for testing class DigitalShapesDecorator.

bool testDigitalShapesDecorator()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;

  trace.beginBlock ( "Testing Unary operation on digital shapes ..." );

  typedef Ball2D< Z2i::Space > ShapeA;
  typedef GaussDigitizer< Z2i::Space, ShapeA > MyGaussDigitizerA;

  typedef ShapeA::Point Point;

  double h = 0.5;

  ShapeA shapeA(-2.501, 0.0, 2.5013);
  ShapeA shapeB(2, 0.0, 2.5013);
  ShapeA shapeC(0.0, 0.0, 2.5);

  MyGaussDigitizerA digShapeA;
  digShapeA.attach( shapeA );
  digShapeA.init( shapeA.getLowerBound(), shapeA.getUpperBound(), h );

  MyGaussDigitizerA digShapeB;
  digShapeB.attach( shapeB );
  digShapeB.init( shapeB.getLowerBound(), shapeB.getUpperBound(), h );

  MyGaussDigitizerA digShapeC;
  digShapeC.attach( shapeC );
  digShapeC.init( shapeC.getLowerBound(), shapeC.getUpperBound(), h );


  typedef DigitalShapesCSG< MyGaussDigitizerA, MyGaussDigitizerA > CSG;
  CSG s_union ( digShapeA );
  s_union.plus( digShapeB );

  CSG s_intersec ( digShapeA );
  s_intersec.plus( digShapeB );
  s_intersec.intersection( digShapeC );

  CSG s_minus ( digShapeA );
  s_minus.minus( digShapeC );


  nbok += (( s_union.orientation( Point( -12, 0 )) == INSIDE )
          || ( s_union.orientation( Point( -12, 0 )) == ON  ))? 0 : 1;
  nbok += (( s_union.orientation( Point( 0, 0 )) == INSIDE )
          || ( s_union.orientation( Point( 0, 0 )) == ON  ))? 1 : 0;
  nbok += (( s_union.orientation( Point( 10, 0 )) == INSIDE )
          || ( s_union.orientation( Point( 10, 0 )) == ON  ))? 0 : 1;
  nbok += (( s_union.orientation( Point( 9, 0 )) == INSIDE )
          || ( s_union.orientation( Point( 9, 0 )) == ON  ))? 1 : 0;
  nbok += (( s_union.orientation( Point( -10, 0 )) == INSIDE )
          || ( s_union.orientation( Point( -10, 0 )) == ON  ))? 1 : 0;


  nbok += (( s_intersec.orientation( Point( -6, 0 )) == INSIDE )
          || ( s_intersec.orientation( Point( -6, 0 )) == ON  ))? 0 : 1;
  nbok += (( s_intersec.orientation( Point( 4, 0 )) == INSIDE )
          || ( s_intersec.orientation( Point( 4, 0 )) == ON  ))? 1 : 0;
  nbok += (( s_intersec.orientation( Point( 6, 0 )) == INSIDE )
          || ( s_intersec.orientation( Point( 6, 0 )) == ON  ))? 0 : 1;
  nbok += (( s_intersec.orientation( Point( 0, 5 )) == INSIDE )
          || ( s_intersec.orientation( Point( 0, 5 )) == ON  ))? 0 : 1;


  nbok += (( s_minus.orientation( Point( -9, 0 )) == INSIDE )
          || ( s_minus.orientation( Point( -9, 0 )) == ON  ))? 1 : 0;
  nbok += (( s_minus.orientation( Point( -2, 0 )) == INSIDE )
          || ( s_minus.orientation( Point( -2, 0 )) == ON  ))? 0 : 1;
  nbok += (( s_minus.orientation( Point( 0, 2 )) == INSIDE )
          || ( s_minus.orientation( Point( 0, 2 )) == ON  ))? 0 : 1;
  nbok += (( s_minus.orientation( Point( -8, 0 )) == INSIDE )
          || ( s_minus.orientation( Point( -8, 0 )) == ON  ))? 1 : 0;

  nb = 13;


  trace.info() << "(" << nbok << "/" << nb << ") "
               << "true == true" << std::endl;
  trace.endBlock();
  return  nbok == nb;
}

// Standard services - public :

int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class DigitalShapesDecorator" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << std::endl;

  bool res = testDigitalShapesDecorator(); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << std::endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //