Aim: model of a CEuclideanOrientedShape from an implicit function from an oriented point cloud. The implicit function is given by the generalized winding number of the oriented point cloud [10] . We use the libIGL implementation. More...

#include <DGtal/shapes/WindingNumbersShape.h>

Public Types
using	Space = TSpace
	Types. More...

using	RealPoint = typename Space::RealPoint

using	RealVector = typename Space::RealVector

using	Orientation = DGtal::Orientation

Public Member Functions
	WindingNumbersShape ()=delete

	WindingNumbersShape (ConstAlias< Eigen::MatrixXd > points, ConstAlias< Eigen::MatrixXd > normals, bool skipPointAreas=false)

	WindingNumbersShape (ConstAlias< Eigen::MatrixXd > points, ConstAlias< Eigen::MatrixXd > normals, const Eigen::VectorXd &areas)

void	setPointAreas (ConstAlias< Eigen::VectorXd > areas)

Orientation	orientation (const RealPoint aPoint, const double threshold=0.3) const

std::vector< Orientation >	orientationBatch (const Eigen::MatrixXd &queries, const double threshold=0.3) const

Data Fields
CountedConstPtrOrConstPtr< Eigen::MatrixXd >	myPoints
	Const alias to the points. More...

CountedConstPtrOrConstPtr< Eigen::MatrixXd >	myNormals
	Const alias to the normals. More...

Eigen::VectorXd	myPointAreas
	Const alias to point area measure. More...

std::vector< std::vector< int > >	myO_PI
	libIGL octree for fast queries data structure More...

Eigen::MatrixXi	myO_CH
	libIGL octree for fast queries data structure More...

Eigen::MatrixXd	myO_CN
	libIGL octree for fast queries data structure More...

Eigen::VectorXd	myO_W
	libIGL octree for fast queries data structure More...

Detailed Description

template<typename TSpace>
struct DGtal::WindingNumbersShape< TSpace >

Aim: model of a CEuclideanOrientedShape from an implicit function from an oriented point cloud. The implicit function is given by the generalized winding number of the oriented point cloud [10] . We use the libIGL implementation.

Description of template class 'WindingNumbersShape'

See also: testWindingNumberShape, windingNumberShape

Template Parameters

TSPace the digital space type (a model of CSpace)

Definition at line 66 of file WindingNumbersShape.h.

Member Typedef Documentation

◆ Orientation

template<typename TSpace >

using DGtal::WindingNumbersShape< TSpace >::Orientation = DGtal::Orientation

Definition at line 72 of file WindingNumbersShape.h.

◆ RealPoint

template<typename TSpace >

using DGtal::WindingNumbersShape< TSpace >::RealPoint = typename Space::RealPoint

Definition at line 70 of file WindingNumbersShape.h.

◆ RealVector

template<typename TSpace >

using DGtal::WindingNumbersShape< TSpace >::RealVector = typename Space::RealVector

Definition at line 71 of file WindingNumbersShape.h.

◆ Space

template<typename TSpace >

using DGtal::WindingNumbersShape< TSpace >::Space = TSpace

Types.

Definition at line 69 of file WindingNumbersShape.h.

Constructor & Destructor Documentation

◆ WindingNumbersShape() [1/3]

template<typename TSpace >

DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape ( )

delete

◆ WindingNumbersShape() [2/3]

template<typename TSpace >

DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape	(	ConstAlias< Eigen::MatrixXd >	points,
		ConstAlias< Eigen::MatrixXd >	normals,
		bool	skipPointAreas = `false`
	)

inline

Construct a WindingNumberShape Euclidean shape from an oriented point cloud. This constructor estimates the area of each point using CGAL.

If the number of points is greater than 20, CGAL is used to estimate the per-sample area (if skipPointAreas is set to false).

Parameters

points	a "nx3" matrix with the sample coordinates.
normals	a "nx3" matrix for the normal vectors.
skipPointAreas	if true, we do not estimate the point areas using CGAL (default: false)

Definition at line 86 of file WindingNumbersShape.h.

     {
       myPoints  = points;
       myNormals = normals;
       myPointAreas = Eigen::VectorXd::Ones(myPoints->rows());
       // Build octree, from libIGL tutorials
       igl::octree(*myPoints,myO_PI,myO_CH,myO_CN,myO_W);
       if (skipPointAreas)
         trace.warning()<<"[WindingNumberShape] Skipping the CGAL point area estimation. By default, point areas are set to 1.0"<<std::endl;
       else
         if (points->rows()> 20)
         {
           Eigen::MatrixXi I;
           igl::knn(*myPoints,(int)points->rows(),myO_PI,myO_CH,myO_CN,myO_W,I);
           // CGAL is only used to help get point areas
           igl::copyleft::cgal::point_areas(*myPoints,I,*myNormals,myPointAreas);
           trace.info()<<"[WindingNumberShape] Min/max point area : "<<myPointAreas.minCoeff()<<" -- "<<myPointAreas.maxCoeff()<<std::endl;
         }
         else
         {
           trace.warning()<<"[WindingNumberShape] Too few points to use CGAL point_areas. Using the constant area setting."<<std::endl;
         }
     }

References DGtal::Trace::info(), DGtal::WindingNumbersShape< TSpace >::myNormals, DGtal::WindingNumbersShape< TSpace >::myO_CH, DGtal::WindingNumbersShape< TSpace >::myO_CN, DGtal::WindingNumbersShape< TSpace >::myO_PI, DGtal::WindingNumbersShape< TSpace >::myO_W, DGtal::WindingNumbersShape< TSpace >::myPointAreas, DGtal::WindingNumbersShape< TSpace >::myPoints, DGtal::trace, and DGtal::Trace::warning().

◆ WindingNumbersShape() [3/3]

template<typename TSpace >

DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape	(	ConstAlias< Eigen::MatrixXd >	points,
		ConstAlias< Eigen::MatrixXd >	normals,
		const Eigen::VectorXd &	areas
	)

inline

Construct a WindingNumberShape Euclidean shape from an oriented point cloud. For this constructor, the area of each point is given by the user.

Parameters

points	a "nx3" matrix with the sample coordinates.
normals	a "nx3" matrix for the normal vectors.
areas	a "n" vector with the area of each point.

Definition at line 118 of file WindingNumbersShape.h.

     {
       myPoints  = points;
       myNormals = normals;
       myPointAreas = areas;
       igl::octree(*myPoints,myO_PI,myO_CH,myO_CN,myO_W);
     }

References DGtal::WindingNumbersShape< TSpace >::myNormals, DGtal::WindingNumbersShape< TSpace >::myO_CH, DGtal::WindingNumbersShape< TSpace >::myO_CN, DGtal::WindingNumbersShape< TSpace >::myO_PI, DGtal::WindingNumbersShape< TSpace >::myO_W, DGtal::WindingNumbersShape< TSpace >::myPointAreas, and DGtal::WindingNumbersShape< TSpace >::myPoints.

Member Function Documentation

◆ orientation()

template<typename TSpace >

Orientation DGtal::WindingNumbersShape< TSpace >::orientation	(	const RealPoint	aPoint,
		const double	threshold = `0.3`
	)		const

inline

Orientation of a point using the winding number value from an oriented pointcloud.

Note: For multiple queries, orientationBatch() should be used.

Parameters

aPoint	[in] a point in space
threshold	[in] the iso-value of the surface of the winding number implicit map (default = 0.3).

Returns: a DGtal::Orientation value

Definition at line 144 of file WindingNumbersShape.h.

     {
       Eigen::MatrixXd queries(1,3);
       queries << aPoint(0) , aPoint(1) , aPoint(2);
       auto singlePoint = orientationBatch(queries, threshold);
       return singlePoint[0];
     }

References aPoint(), and DGtal::WindingNumbersShape< TSpace >::orientationBatch().

◆ orientationBatch()

template<typename TSpace >

std::vector<Orientation> DGtal::WindingNumbersShape< TSpace >::orientationBatch	(	const Eigen::MatrixXd &	queries,
		const double	threshold = `0.3`
	)		const

inline

Orientation of a set of points (queries) using the winding number value from an oriented pointcloud.

Parameters

queries	[in] a "nx3" matrix with the query points in space.
threshold	[in] the iso-value of the surface of the winding number implicit map (default = 0.3).

Returns: a DGtal::Orientation value vector for each query point.

Definition at line 159 of file WindingNumbersShape.h.

     {
       Eigen::VectorXd W;
       std::vector<Orientation> results( queries.rows(), DGtal::OUTSIDE );
       Eigen::MatrixXd O_CM;
       Eigen::VectorXd O_R;
       Eigen::MatrixXd O_EC;
     
       //Checking if the areas
       igl::fast_winding_number(*myPoints,*myNormals,myPointAreas,myO_PI,myO_CH,2,O_CM,O_R,O_EC);
       igl::fast_winding_number(*myPoints,*myNormals,myPointAreas,myO_PI,myO_CH,O_CM,O_R,O_EC,queries,2,W);
       
       //Reformating the output
       for(auto i=0u; i < queries.rows(); ++i)
       {
         if (std::abs(W(i)) < threshold )
           results[i] = DGtal::OUTSIDE;
         else
           if (std::abs(W(i)) > threshold)
             results[i] = DGtal::INSIDE;
           else
             results[i] = DGtal::ON;
       }
       return results;
     }

References DGtal::INSIDE, DGtal::WindingNumbersShape< TSpace >::myNormals, DGtal::WindingNumbersShape< TSpace >::myO_CH, DGtal::WindingNumbersShape< TSpace >::myO_PI, DGtal::WindingNumbersShape< TSpace >::myPointAreas, DGtal::WindingNumbersShape< TSpace >::myPoints, DGtal::ON, and DGtal::OUTSIDE.

Referenced by DGtal::WindingNumbersShape< TSpace >::orientation().

◆ setPointAreas()

template<typename TSpace >

void DGtal::WindingNumbersShape< TSpace >::setPointAreas ( ConstAlias< Eigen::VectorXd > areas )

inline

Set the area map for each point.

Parameters

areas a Eigen vector of estimated area for each input point.

Definition at line 131 of file WindingNumbersShape.h.

     {
       myPointAreas = areas;
     }

References DGtal::WindingNumbersShape< TSpace >::myPointAreas.

Field Documentation

◆ myNormals

template<typename TSpace >

CountedConstPtrOrConstPtr<Eigen::MatrixXd> DGtal::WindingNumbersShape< TSpace >::myNormals

Const alias to the normals.

Definition at line 189 of file WindingNumbersShape.h.

Referenced by DGtal::WindingNumbersShape< TSpace >::orientationBatch(), and DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape().

◆ myO_CH

template<typename TSpace >

Eigen::MatrixXi DGtal::WindingNumbersShape< TSpace >::myO_CH

libIGL octree for fast queries data structure

Definition at line 196 of file WindingNumbersShape.h.

Referenced by DGtal::WindingNumbersShape< TSpace >::orientationBatch(), and DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape().

◆ myO_CN

template<typename TSpace >

Eigen::MatrixXd DGtal::WindingNumbersShape< TSpace >::myO_CN

libIGL octree for fast queries data structure

Definition at line 198 of file WindingNumbersShape.h.

Referenced by DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape().

◆ myO_PI

template<typename TSpace >

std::vector<std::vector<int > > DGtal::WindingNumbersShape< TSpace >::myO_PI

libIGL octree for fast queries data structure

Definition at line 194 of file WindingNumbersShape.h.

Referenced by DGtal::WindingNumbersShape< TSpace >::orientationBatch(), and DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape().

◆ myO_W

template<typename TSpace >

Eigen::VectorXd DGtal::WindingNumbersShape< TSpace >::myO_W

libIGL octree for fast queries data structure

Definition at line 200 of file WindingNumbersShape.h.

Referenced by DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape().

◆ myPointAreas

template<typename TSpace >

Eigen::VectorXd DGtal::WindingNumbersShape< TSpace >::myPointAreas

Const alias to point area measure.

Definition at line 191 of file WindingNumbersShape.h.

Referenced by DGtal::WindingNumbersShape< TSpace >::orientationBatch(), DGtal::WindingNumbersShape< TSpace >::setPointAreas(), and DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape().

◆ myPoints

template<typename TSpace >

CountedConstPtrOrConstPtr<Eigen::MatrixXd> DGtal::WindingNumbersShape< TSpace >::myPoints

Const alias to the points.

Definition at line 187 of file WindingNumbersShape.h.

Referenced by DGtal::WindingNumbersShape< TSpace >::orientationBatch(), and DGtal::WindingNumbersShape< TSpace >::WindingNumbersShape().

The documentation for this struct was generated from the following file:

WindingNumbersShape.h

Public Types

Public Member Functions

Data Fields

Detailed Description

template<typename TSpace> struct DGtal::WindingNumbersShape< TSpace >

Member Typedef Documentation

◆ Orientation

◆ RealPoint

◆ RealVector

◆ Space

Constructor & Destructor Documentation

◆ WindingNumbersShape() [1/3]

◆ WindingNumbersShape() [2/3]

◆ WindingNumbersShape() [3/3]

Member Function Documentation

◆ orientation()

◆ orientationBatch()

◆ setPointAreas()

Field Documentation

◆ myNormals

◆ myO_CH

◆ myO_CN

◆ myO_PI

◆ myO_W

◆ myPointAreas

◆ myPoints

template<typename TSpace>
struct DGtal::WindingNumbersShape< TSpace >