exampleQuickHull3D.cpp

Go to the documentation of this file.

 
#include "DGtal/base/Common.h"
#include "DGtal/geometry/volumes/ConvexityHelper.h"
#include "DGtal/shapes/SurfaceMesh.h"
#include "DGtal/io/writers/SurfaceMeshWriter.h"
#include "ConfigExamples.h"
 
double rand01() { return (double) rand() / (double) RAND_MAX; }
 
using namespace DGtal;
using namespace DGtal::Z3i;
int main( int argc, char* argv[] )
{
  int precision = argc > 1 ? atoi( argv[ 1 ] ) : 100;  // precision
  std::string inputFilename = argc > 2
    ? std::string( argv[ 2 ] )
    : examplesPath + "samples/bunny.dat" ;
  std::vector< RealPoint > points;
  std::ifstream finput( inputFilename.c_str() );
  std::string linestr;
  while ( std::getline( finput, linestr ) )
    {
      std::istringstream iss( linestr );
      double a, b, c;
      if ( ! (iss >> a >> b >> c) ) break;
      points.push_back( RealPoint( a, b, c ) );
    }
  trace.info() << "Read " << points.size() << " 3D points." << std::endl;
 
  // Build rational polytope
  typedef ConvexityHelper< 3 > Helper;
  const auto polytope
    = Helper::computeRationalPolytope( points, precision );
  trace.info() << polytope << std::endl;
 
  // Build the boundary of the convex hull as a surface mesh
  SurfaceMesh< RealPoint, RealVector > mesh;
  bool ok = Helper::computeConvexHullBoundary( mesh, points, precision );
  trace.info() << mesh << std::endl;
  std::ofstream out( "qhull-mesh.obj" );
  SurfaceMeshWriter< RealPoint, RealVector >::writeOBJ( out, mesh );
  out.close();
 
  // Build the boundary of the convex hull as a polygonal surface
  PolygonalSurface< RealPoint > polysurf;
  bool ok2 = Helper::computeConvexHullBoundary( polysurf, points, precision );
  trace.info() << polysurf << std::endl;
 
  // Build the convex hull as a convex cell complex.
  ConvexCellComplex< RealPoint > cvx_complex;
  bool ok3 = Helper::computeConvexHullCellComplex( cvx_complex, points, precision );
  trace.info() << cvx_complex << std::endl;
  
  return ( ok && ok2 && ok3 ) ? 0 : 1;
}