testBallQuadViewer.cpp

 
#include <iostream>
#include "DGtal/base/Common.h"
#include "ConfigTest.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/topology/ImplicitDigitalSurface.h"
#include "DGtal/io/viewers/Viewer3D.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;
};
 
 
bool testBallQuad(int argc, char **argv)
{
  unsigned int nbok = 0;
  unsigned int nb = 0;
 
  QApplication application(argc, argv);
 
  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;
 
  Viewer3D<Space,KSpace> viewer(K);
  viewer.setWindowTitle("simpleViewer");
  viewer.show();
 
 
  for ( ConstIterator it = boundary.begin(), it_end = boundary.end();
        it != it_end; ++it )
    {
      ++nbsurfels;
 
      Display3DFactory<>::drawOrientedSurfelWithNormal(viewer,
                                                       *it,
                                                       viewer.embedKS(*it).getNormalized());
    }
 
  trace.info() << nbsurfels << " surfels found." << std::endl;
  viewer  << Display3D<Space, KSpace>::updateDisplay;
 
  bool res = application.exec();
 
  return res;
}
 
// 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;
 
  bool res = testBallQuad(argc,argv); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //