testMeasureSet.cpp

#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/shapes/Shapes.h"
#include "DGtal/shapes/ShapeFactory.h"
 
#include "DGtal/geometry/volumes/estimation/Measure.h"

 
using namespace std;
using namespace DGtal;

// Functions for testing class Measure.
 
bool testConcept()
{
  
  return true;
}
 
 

bool testMeasure()
{
  trace.beginBlock ( "Testing Measure ..." );
 
 
  Z3i::Point a(0,0);
  Z3i::Point b(64,64,64);
  Z3i::Point c(32,32,32);
  Z3i::Domain domain(a,b);
  Z3i::DigitalSet set(domain);
  
  Shapes<Z3i::Domain>::euclideanShaper( set,
                                        ImplicitBall<Z3i::Space>( c, 10));
  
  Measure< Z3i::DigitalSet> measure;
 
  trace.info() << "Input set= "<<set<<std::endl;
 
  trace.info()<<measure<<std::endl;
 
  
  measure.init(10, set);
  trace.info() << "Volume (h=10) "<<measure.eval()<<std::endl;
  measure.init(100, set);
  trace.info() << "Volume (h=100) "<<measure.eval()<<std::endl;
  
  trace.info()<<measure<<std::endl;
 
  trace.endBlock();
  
  return true;
}

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