testImplicitShape.cpp

Go to the documentation of this file.

 
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/shapes/Shapes.h"
#include "DGtal/shapes/ShapeFactory.h"
#include "DGtal/io/boards/Board2D.h"
#include "DGtal/io/colormaps/GrayscaleColorMap.h"
#include "DGtal/io/writers/VolWriter.h"
#include "DGtal/images/ImageContainerBySTLVector.h"
 
using namespace std;
using namespace DGtal;
 
// Functions for testing class ImplicitShape.
 
bool testImplicitShape()
{
  trace.beginBlock ( "Testing implicit shaper ..." );
  Z2i::Point a(0,0);
  Z2i::Point b(64,64);
  Z2i::Point c(32,32);
  
  Board2D board;
  
  Z2i::Domain domain(a,b);
  Z2i::DigitalSet set(domain);
  
  Shapes<Z2i::Domain>::euclideanShaper( set,
             ImplicitBall<Z2i::Space>( c, 10));
  board << set;
  board.saveSVG("implicitball.svg");
 
  set.clear();
  board.clear();
  Shapes<Z2i::Domain>::euclideanShaper( set,
             ImplicitHyperCube<Z2i::Space>( c, 10));
  board << set;
  board.saveSVG("implicitcube.svg");
  
 
  set.clear();
  board.clear();
  Shapes<Z2i::Domain>::euclideanShaper( set,
             ImplicitNorm1Ball<Z2i::Space>( c, 10));
  board << set;
  board.saveSVG("implicitlosange.svg");
  
  set.clear();
  board.clear();
  Shapes<Z2i::Domain>::euclideanShaper( set,
             ImplicitRoundedHyperCube<Z2i::Space>( c, 10, 1));
  board << set;
  board.saveSVG("implicitrounded-1.svg");
  
 
  set.clear();
  board.clear();
  Shapes<Z2i::Domain>::euclideanShaper( set,
             ImplicitRoundedHyperCube<Z2i::Space>( c, 10, 2.5));
  board << set;
  board.saveSVG("implicitrounded-2.5.svg");
  
 
  return true;
}
bool testImplicitShape3D()
{
  trace.beginBlock ( "Testing implicit shaper in 3D..." );
  Z3i::Point a(0,0);
  Z3i::Point b(64,64,64);
  Z3i::Point c(32,32,32);
  
  Z3i::Domain domain(a,b);
 
  typedef ImageContainerBySTLVector<Z3i::Domain, DGtal::uint8_t> Image;
 
  Image image(domain);
  Z3i::DigitalSet set(domain);
 
  Shapes<Z3i::Domain>::euclideanShaper( set,
             ImplicitRoundedHyperCube<Z3i::Space>( c, 10, 2.5));
 
  for(Z3i::DigitalSet::ConstIterator it=set.begin(), itend=set.end();
      it != itend;
      ++it)
    {
      image.setValue((*it), 128);
    }
 
  VolWriter<Image>::exportVol("implicitrounded.vol",image);
  
  return true;
}
 
// Standard services - public :
 
int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class ImplicitShape" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;
 
  bool res = testImplicitShape() && testImplicitShape3D(); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //