testBallQuadViewer.cpp

Go to the documentation of this file.

#include <iostream>
#include "DGtal/base/Common.h"
#include "ConfigTest.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/topology/ImplicitDigitalSurface.h"
#include "DGtal/io/viewers/PolyscopeViewer.h"
 
using namespace std;
using namespace DGtal;

// Functions for testing class BallQuad.

template <typename TPoint3>
struct ImplicitDigitalBall3 {
  typedef TPoint3 Point;
  inline
  ImplicitDigitalBall3( double r )
    : myR( r )
  {}
  inline
  bool operator()( const TPoint3 & p ) const
  {
    double x = ( (double) p[ 0 ]  );
    double y = ( (double) p[ 1 ]  );
    double z = ( (double) p[ 2 ]  );
    return ( x*x + y*y + z*z  -myR*myR) <= 0.0;
  }
  double myR;
};
 
 
void testBallQuad()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;
 
  trace.beginBlock ( "Testing... Ball with quadnormal");
  using namespace Z3i;
  typedef ImplicitDigitalBall3<Point> ImplicitDigitalBall;
  typedef ImplicitDigitalSurface<KSpace,ImplicitDigitalBall> Boundary;
  typedef Boundary::SurfelConstIterator ConstIterator;
  typedef Boundary::Surfel Surfel;
  Point p1( -100, -100, -100 );
  Point p2( 100, 100, 100 );
  KSpace K;
  nbok += K.init( p1, p2, true ) ? 1 : 0;
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
               << "K.init() is ok" << std::endl;
  ImplicitDigitalBall ball( 60.0 );
  Surfel bel = Surfaces<KSpace>::findABel( K, ball, 10000 );
  Boundary boundary( K, ball,
                     SurfelAdjacency<KSpace::dimension>( true ), bel );
  unsigned int nbsurfels = 0;
 
  PolyscopeViewer<Space,KSpace> viewer(K);
  viewer.allowReuseList = true;
 
  for ( ConstIterator it = boundary.begin(), it_end = boundary.end();
        it != it_end; ++it )
  {
    ++nbsurfels;
    viewer << *it;
  }
 
  trace.info() << nbsurfels << " surfels found." << std::endl;
 
  viewer.show();
}

// Standard services - public :
 
int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class BallQuad" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;
 
  testBallQuad(); // && ... other tests
  trace.emphase() << "Passed." << endl;
  trace.endBlock();
  return 0;
}
//                                                                           //