testBinomialConvolver.cpp

#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/geometry/curves/BinomialConvolver.h"

using namespace std;
using namespace DGtal;

// Functions for testing class Signal.

bool testBinomialConvolver()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;
  
  trace.beginBlock ( "Testing block ..." );
  typedef PointVector<2, double> RealPoint;
  std::vector< RealPoint > points;
#ifdef CPP11_INITIALIZER_LIST
  points.push_back( RealPoint( { 0.0, 0.0 } ) ); 
  points.push_back( RealPoint( { 1.0, 0.0 } ) ); 
  points.push_back( RealPoint( { 2.0, 0.0 } ) ); 
  points.push_back( RealPoint( { 2.0, 1.0 } ) ); 
  points.push_back( RealPoint( { 2.0, 2.0 } ) ); 
  points.push_back( RealPoint( { 1.0, 2.0 } ) ); 
  points.push_back( RealPoint( { 0.0, 2.0 } ) ); 
  points.push_back( RealPoint( { 0.0, 1.0 } ) ); 
#else
   points.push_back( RealPoint(  0.0, 0.0  ) ); 
  points.push_back( RealPoint(  1.0, 0.0  ) ); 
  points.push_back( RealPoint(  2.0, 0.0  ) ); 
  points.push_back( RealPoint(  2.0, 1.0  ) ); 
  points.push_back( RealPoint( 2.0, 2.0  ) ); 
  points.push_back( RealPoint( 1.0, 2.0  ) ); 
  points.push_back( RealPoint(  0.0, 2.0  ) ); 
  points.push_back( RealPoint(  0.0, 1.0  ) ); 
#endif
  
  typedef std::vector< RealPoint >::const_iterator ConstIteratorOnPoints;
  typedef BinomialConvolver<ConstIteratorOnPoints, double> MyBinomialConvolver;
  
  for ( unsigned int n = 1; n < 10; ++n )
    {
      trace.info() << "Binomial convolver n=" << n << std::endl;
      MyBinomialConvolver bcc( n );
      bcc.init( 1.0, points.begin(), points.end(), true );
      for ( unsigned int i = 0; i < 8; ++i )
        std::cout << i
            << " " << bcc.x( i ).first
            << " " << bcc.x( i ).second
            << " " << bcc.tangent( i ).first
            << " " << bcc.tangent( i ).second
            << " " << bcc.curvature( i )
            << std::endl;
    }
  unsigned int n = MyBinomialConvolver::suggestedSize( 1.0, points.begin(), points.end() );
  trace.info() << "Binomial convolver suggested n=" 
         << n
         << std::endl;

  typedef TangentFromBinomialConvolverFunctor< MyBinomialConvolver, RealPoint >
    TangentBCFct;
  typedef CurvatureFromBinomialConvolverFunctor< MyBinomialConvolver, double >
    CurvatureBCFct;
  BinomialConvolverEstimator< MyBinomialConvolver, TangentBCFct> tgtEstimator;
  BinomialConvolverEstimator< MyBinomialConvolver, CurvatureBCFct> curvEstimator;
  tgtEstimator.init( 1.0, points.begin(), points.end(), true );
  curvEstimator.init( 1.0, points.begin(), points.end(), true );
  for ( ConstIteratorOnPoints it = points.begin(), it_end = points.end();
  it != it_end; ++it )
    {
      std::cout << *it 
    << " " << tgtEstimator.eval( it ) 
    << " " << curvEstimator.eval( it ) 
    << std::endl;
    }
  nbok += true ? 1 : 0; 
  nb++;
  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 BinomialConvolver" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;

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