testDigitalShapesDecorator.cpp

Go to the documentation of this file.

 
#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 Z2i::Space Space;
  typedef Ball2D< Space > ShapeA;
  typedef GaussDigitizer< Space, ShapeA > MyGaussDigitizerA;
 
  typedef Space::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;
}
//                                                                           //