testVoronoiCovarianceMeasure.cpp File Reference

#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"

Go to the source code of this file.

Functions
bool	testVoronoiCovarianceMeasure ()
int	main (int, char **)

Detailed Description

This program is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Author

Jacques-Olivier Lachaud (jacqu.nosp@m.es-o.nosp@m.livie.nosp@m.r.la.nosp@m.chaud.nosp@m.@uni.nosp@m.v-sav.nosp@m.oie..nosp@m.fr ) Laboratory of Mathematics (CNRS, UMR 5127), University of Savoie, France

Date

2014/02/09

Functions for testing class VoronoiCovarianceMeasure, and also IntegralInvariantEstimator.

This file is part of the DGtal library.

This program is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Author

Jacques-Olivier Lachaud (jacqu.nosp@m.es-o.nosp@m.livie.nosp@m.r.la.nosp@m.chaud.nosp@m.@uni.nosp@m.v-sav.nosp@m.oie..nosp@m.fr ) Laboratory of Mathematics (CNRS, UMR 5127), University of Savoie, France

Date

2014/02/09

Functions for testing class VoronoiCovarianceMeasure.

This file is part of the DGtal library.

Definition in file testVoronoiCovarianceMeasure.cpp.

Function Documentation

main()

int main	(	int	,
		char **
	)

Definition at line 112 of file testVoronoiCovarianceMeasure.cpp.

{
  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;
}

References DGtal::Trace::beginBlock(), DGtal::Trace::emphase(), DGtal::Trace::endBlock(), and testVoronoiCovarianceMeasure().

testVoronoiCovarianceMeasure()

bool testVoronoiCovarianceMeasure

(

)

Example of a test. To be completed.

Definition at line 48 of file testVoronoiCovarianceMeasure.cpp.

{
  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;
}

References DGtal::HyperRectDomain< TSpace >::begin(), DGtal::Trace::beginBlock(), domain, DGtal::HyperRectDomain< TSpace >::end(), DGtal::Trace::endBlock(), DGtal::Trace::info(), and DGtal::HyperRectDomain< TSpace >::size().

Referenced by main().