testVoronoiCovarianceMeasure.cpp

Go to the documentation of this file.

 
#include <iostream>
#include <vector>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/geometry/volumes/estimation/VoronoiCovarianceMeasure.h"
#include "DGtal/geometry/tools/SpatialCubicalSubdivision.h"
 
 
// Functions for testing class VoronoiCovarianceMeasure
 
bool testVoronoiCovarianceMeasure()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;
  
  using namespace DGtal;
  using namespace DGtal::Z3i; // gets Space, Point, Domain
  trace.beginBlock ( "testVoronoiCovarianceMeasure" );
  typedef ExactPredicateLpSeparableMetric<Space,2> Metric;
  typedef VoronoiCovarianceMeasure<Space, Metric> VCM;
  typedef VCM::MatrixNN Matrix;
 
  Point a(0,0);
  Point c(32,32,32);
  Domain domain(a,c);
  std::vector<Point> pts;
  pts.push_back( Point( 10,10,10 ) );
  pts.push_back( Point( 10,10,11 ) );
  pts.push_back( Point( 10,11,11 ) );
  pts.push_back( Point( 11,11,12 ) );
  pts.push_back( Point( 11,12,12 ) );
  pts.push_back( Point( 11,12,12 ) );
  pts.push_back( Point( 12,12,12 ) );
  pts.push_back( Point( 20,20,15 ) );
  pts.push_back( Point( 30,10,25 ) );
  pts.push_back( Point( 25,25,20 ) );
  Metric l2;
  VCM vcm( 5.0, 4.0, l2, true );
  vcm.init( pts.begin(), pts.end() );
  Domain d = vcm.domain();
  double sum_dist = 0.0;
  std::map<Point,int> sizeCells;
  for ( Domain::ConstIterator it = d.begin(), itE = d.end(); it != itE; ++it )
    {
      sum_dist += l2( *it, vcm.voronoiMap()( *it ) );
      sizeCells[ vcm.voronoiMap()( *it ) ] += 1;
    }
  double avg_dist = ( sum_dist / d.size() );
  trace.info() << "Average distance is " << avg_dist << std::endl;
  nbok += ( 9 <= avg_dist && avg_dist <= 11 ) ? 1 : 0; 
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
               << "9 <= " << avg_dist << " <= 11 " << std::endl;
  for ( std::map<Point,int>::const_iterator it = sizeCells.begin(), itE = sizeCells.end();
        it != itE; ++it )
    trace.info() << "Nb(" << it->first << ")=" << it->second << std::endl;
  nbok += sizeCells.size() == 9 ? 1 : 0;
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
               << "sizeCells.size() == 10" << std::endl;
 
  functors::HatPointFunction< Point, double > chi_r( 1.0, 4.0 );
  Matrix vcm_r = vcm.measure( chi_r, Point( 10,10,10 ) );
  trace.info() << "- vcm_r.row(0) = " << vcm_r.row( 0 ) << std::endl;
  trace.info() << "- vcm_r.row(1) = " << vcm_r.row( 1 ) << std::endl;
  trace.info() << "- vcm_r.row(2) = " << vcm_r.row( 2 ) << std::endl;
  trace.endBlock();
  
  return nbok == nb;
}
 
// Standard services - public :
 
int main( int /* argc */, char** /* argv */ )
{
  using namespace std;
  using namespace DGtal;
  trace.beginBlock ( "Testing VoronoiCovarianceMeasure ..." );
  bool res = testVoronoiCovarianceMeasure();
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //