CGAL::Poisson_reconstruction_function<GeomTraits>

Definition

Given a set of 3D points with oriented normals sampled on the boundary of a 3D solid, the Poisson Surface Reconstruction method [KBH05] solves for an approximate indicator function of the inferred solid, whose gradient best matches the input normals. The output scalar function, represented in an adaptive octree, is then iso-contoured using an adaptive marching cubes.

Poisson_reconstruction_function implements a variant of this algorithm which solves for a piecewise linear function on a 3D Delaunay triangulation instead of an adaptive octree and uses the TAUCS sparse linear solver.

#include <CGAL/Poisson_reconstruction_function.h>

Parameters

template<
class Gt>
class Poisson_reconstruction_function;

Parameters


Gt: Geometric traits class.

Is Model for the Concepts

Model of the ImplicitFunction concept.

Types

Poisson_reconstruction_function<GeomTraits>::Geom_traits
Geometric traits class.

Poisson_reconstruction_function<GeomTraits>::FT
typedef to Geom_traits::FT

Poisson_reconstruction_function<GeomTraits>::Point
typedef to Geom_traits::Point_3

Poisson_reconstruction_function<GeomTraits>::Vector
typedef to Geom_traits::Vector_3

Poisson_reconstruction_function<GeomTraits>::Sphere
typedef to Geom_traits::Sphere_3

Creation

template<typename InputIterator, typename PointPMap, typename NormalPMap>
Poisson_reconstruction_function<GeomTraits> fct ( InputIterator first,
InputIterator beyond,
PointPMap point_pmap,
NormalPMap normal_pmap);
Creates a Poisson implicit function from the [first, beyond) range of points.
Template Parameters: 
InputIterator: iterator over input points. PointPMap: is a model of boost::ReadablePropertyMap with a value_type = Point_3. It can be omitted if InputIterator value_type is convertible to Point_3. NormalPMap: is a model of boost::ReadablePropertyMap with a value_type = Vector_3.
Parameters: 
first: iterator over the first input point. beyond: past-the-end iterator over the input points. point_pmap: property map to access the position of an input point. normal_pmap: property map to access the oriented normal of an input point.

Operations

Sphere fct.bounding_sphere () const Returns a sphere bounding the inferred surface.
bool fct.compute_implicit_function () The function compute_implicit_function() must be called after each insertion of oriented points. It computes the piecewise linear scalar function operator() by: applying Delaunay refinement, solving for operator() at each vertex of the triangulation with a sparse linear solver, and shifting and orienting operator() such that it is 0 at all input points and negative inside the inferred surface.
Returns false if the linear solver fails.
FT fct.operator() ( const Point& p) const
ImplicitFunction interface: evaluates the implicit function at a given 3D query point.
Point fct.get_inner_point () const Returns a point located inside the inferred surface.

Example

See poisson_reconstruction_example.cpp.