Aim: Estimates normal vector using CGAL Jet Fitting and Monge Form. More...

#include <DGtal/geometry/surfaces/estimation/estimationFunctors/MongeJetFittingNormalVectorEstimator.h>

Public Types
typedef TSurfel	Surfel

typedef TEmbedder	SCellEmbedder

typedef SCellEmbedder::RealPoint	RealPoint

typedef RealPoint	Quantity

typedef CGAL::Cartesian< double >	CGALKernel

typedef CGALKernel::Point_3	CGALPoint

typedef CGALKernel::Vector_3	CGALVector

typedef CGAL::Monge_via_jet_fitting< CGALKernel >	CGALMongeViaJet

typedef CGALMongeViaJet::Monge_form	CGALMongeForm

Public Member Functions
	MongeJetFittingNormalVectorEstimator (ConstAlias< SCellEmbedder > anEmbedder, const double h, unsigned int d=4)

void	pushSurfel (const Surfel &aSurf, const double aDistance)

Quantity	eval ()

void	reset ()

Private Attributes
const SCellEmbedder *	myEmbedder
	Alias of the geometrical embedder.

std::vector< CGALPoint >	myPoints
	Array of CGAL points.

double	myH

unsigned int	myD
	Degree of the polynomial surface to fit.

Detailed Description

template<typename TSurfel, typename TEmbedder>
class DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >

Aim: Estimates normal vector using CGAL Jet Fitting and Monge Form.

Description of template class 'MongeJetFittingNormalVectorEstimator'

As discussed in MongeJetFittingMeanCurvatureEstimator, only the estimated normal vector is given up to a sign. Hence, depending on the local orientation of the surface, you would problably have to reverse the estimated vector \( \vec{n} \) to \( -\vec{n}\).

model of concepts::CLocalEstimatorFromSurfelFunctor

Template Parameters

TSurfel	type of surfels
TEmbedder	type of functors which embed surfel to \( \mathbb{R}^3\)

Definition at line 81 of file MongeJetFittingNormalVectorEstimator.h.

Member Typedef Documentation

◆ CGALKernel

template<typename TSurfel , typename TEmbedder >

typedef CGAL::Cartesian<double> DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::CGALKernel

Definition at line 90 of file MongeJetFittingNormalVectorEstimator.h.

◆ CGALMongeForm

template<typename TSurfel , typename TEmbedder >

typedef CGALMongeViaJet::Monge_form DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::CGALMongeForm

Definition at line 94 of file MongeJetFittingNormalVectorEstimator.h.

◆ CGALMongeViaJet

template<typename TSurfel , typename TEmbedder >

typedef CGAL::Monge_via_jet_fitting<CGALKernel> DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::CGALMongeViaJet

Definition at line 93 of file MongeJetFittingNormalVectorEstimator.h.

◆ CGALPoint

template<typename TSurfel , typename TEmbedder >

typedef CGALKernel::Point_3 DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::CGALPoint

Definition at line 91 of file MongeJetFittingNormalVectorEstimator.h.

◆ CGALVector

template<typename TSurfel , typename TEmbedder >

typedef CGALKernel::Vector_3 DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::CGALVector

Definition at line 92 of file MongeJetFittingNormalVectorEstimator.h.

◆ Quantity

template<typename TSurfel , typename TEmbedder >

typedef RealPoint DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::Quantity

Definition at line 88 of file MongeJetFittingNormalVectorEstimator.h.

◆ RealPoint

template<typename TSurfel , typename TEmbedder >

typedef SCellEmbedder::RealPoint DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::RealPoint

Definition at line 87 of file MongeJetFittingNormalVectorEstimator.h.

◆ SCellEmbedder

template<typename TSurfel , typename TEmbedder >

typedef TEmbedder DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::SCellEmbedder

Definition at line 86 of file MongeJetFittingNormalVectorEstimator.h.

◆ Surfel

template<typename TSurfel , typename TEmbedder >

typedef TSurfel DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::Surfel

Definition at line 85 of file MongeJetFittingNormalVectorEstimator.h.

Constructor & Destructor Documentation

◆ MongeJetFittingNormalVectorEstimator()

template<typename TSurfel , typename TEmbedder >

DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::MongeJetFittingNormalVectorEstimator	(	ConstAlias< SCellEmbedder >	anEmbedder,
		const double	h,
		unsigned int	d = 4 )

inline

Constructor.

Parameters

anEmbedder	embedder to map surfel to R^n.
h	grid step
d	degree of the polynomial surface to fit (default d=4).

Definition at line 103 of file MongeJetFittingNormalVectorEstimator.h.

                                                                            :
      myEmbedder(&anEmbedder), myH(h), myD(d)
    {
      FATAL_ERROR_MSG(d>=2, "Polynomial surface degree must be greater than 2");
    }

Member Function Documentation

◆ eval()

template<typename TSurfel , typename TEmbedder >

Quantity DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::eval ( )

inline

Evaluate the normal vector from Monge form.

Returns: the mean curvature

Definition at line 131 of file MongeJetFittingNormalVectorEstimator.h.

    {
      CGALMongeForm monge_form;
      CGALMongeViaJet monge_fit;
 
      monge_form = monge_fit(myPoints.begin() , myPoints.end(), myD, (2<myD)? myD : 2);
 
      CGALVector v= monge_form.normal_direction (  );
 
      return RealPoint(v.x(),v.y(),v.z());
    }

References DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::myD, and DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::myPoints.

◆ pushSurfel()

template<typename TSurfel , typename TEmbedder >

void DGtal::functors::MongeJetFittingNormalVectorEstimator< TSurfel, TEmbedder >::pushSurfel	(	const Surfel &	aSurf,
		const double	aDistance )

inline

Add the geometrical embedding of a surfel to the point list

Parameters

aSurf	a surfel to add
aDistance	distance of aSurf to the neighborhood boundary

Definition at line 116 of file MongeJetFittingNormalVectorEstimator.h.

    {
      BOOST_VERIFY(aDistance == aDistance);
 
      RealPoint p = myEmbedder->operator()(aSurf);
      CGALPoint pp(p[0]*myH,p[1]*myH,p[2]*myH);
      myPoints.push_back(pp);
    }