DGtal  1.5.beta
IntegralInvariantVolumeEstimator.h
1 
17 #pragma once
18 
34 #if defined(IntegralInvariantVolumeEstimator_RECURSES)
35 #error Recursive header files inclusion detected in IntegralInvariantVolumeEstimator.h
36 #else // defined(IntegralInvariantVolumeEstimator_RECURSES)
38 #define IntegralInvariantVolumeEstimator_RECURSES
39 
40 #if !defined IntegralInvariantVolumeEstimator_h
42 #define IntegralInvariantVolumeEstimator_h
43 
45 // Inclusions
46 #include <iostream>
47 #include "DGtal/base/Common.h"
48 
49 #include "DGtal/kernel/BasicPointFunctors.h"
50 #include "DGtal/kernel/CPointPredicate.h"
51 #include "DGtal/topology/CCellFunctor.h"
52 #include "DGtal/topology/CCellularGridSpaceND.h"
53 #include "DGtal/geometry/surfaces/FunctorOnCells.h"
54 
55 #include "DGtal/shapes/GaussDigitizer.h"
56 #include "DGtal/shapes/Shapes.h"
57 
58 #include "DGtal/geometry/surfaces/DigitalSurfaceConvolver.h"
59 #include "DGtal/geometry/surfaces/estimation/IIGeometricFunctors.h"
60 #include "DGtal/shapes/EuclideanShapesDecorator.h"
61 
62 #include "DGtal/shapes/implicit/ImplicitBall.h"
64 
65 
66 namespace DGtal
67 {
68 
70 // template class IntegralInvariantVolumeEstimator
112 template <typename TKSpace, typename TPointPredicate, typename TVolumeFunctor>
114 {
115 public:
117  typedef TKSpace KSpace;
118  typedef TPointPredicate PointPredicate;
119  typedef TVolumeFunctor VolumeFunctor;
120 
123 
124  typedef typename KSpace::Space Space;
126  typedef typename Space::Point Point;
127  typedef typename Space::RealPoint RealPoint;
128  typedef typename Space::RealVector RealVector;
130  typedef typename KSpace::SCell Spel;
131  typedef typename KSpace::Surfel Surfel;
132  typedef typename KSpace::SurfelSet SurfelSet;
133  typedef typename SurfelSet::const_iterator ConstIteratorKernel;
134 
136  typedef typename VolumeFunctor::Quantity Quantity;
138  typedef int Value;
139 
146 
147 
153 
158  typedef typename Matrix::Component Component;
159  typedef double Scalar;
163  typename VolumeFunctor::Argument >::value ));
164 
165 
166  // ----------------------- Standard services ------------------------------
167 public:
168 
177 
188  ConstAlias< PointPredicate > aPointPredicate,
189  VolumeFunctor fct = VolumeFunctor() );
190 
195 
201 
207  Self& operator= ( const Self& other );
208 
212  void clear();
213 
214  // ----------------------- Interface --------------------------------------
215 public:
216 
218  Scalar h() const;
219 
228  ConstAlias<PointPredicate> aPointPredicate );
229 
235  void setParams( const double dRadius );
236 
245  template <typename SurfelConstIterator>
246  void init( const double _h, SurfelConstIterator itb, SurfelConstIterator ite );
247 
262  template< typename SurfelConstIterator >
263  Quantity eval ( SurfelConstIterator it ) const;
264 
265 
286  template <typename OutputIterator, typename SurfelConstIterator>
287  OutputIterator eval( SurfelConstIterator itb,
288  SurfelConstIterator ite,
289  OutputIterator result ) const;
290 
295  void selfDisplay ( std::ostream & out ) const;
296 
301  bool isValid() const;
302 
303  // ------------------------- Private Datas --------------------------------
304 private:
305 
308  std::vector< PairIterators > myKernels;
309  std::vector< DigitalSet * > myKernelsSet;
319 
320 private:
321 
322 
323 }; // end of class IntegralInvariantVolumeEstimator
324 
331  template <typename TKSpace, typename TPointPredicate, typename TVolumeFunctor>
332  std::ostream&
333  operator<< ( std::ostream & out,
335 
336 } // namespace DGtal
337 
338 
340 // Includes inline functions.
341 #include "DGtal/geometry/surfaces/estimation/IntegralInvariantVolumeEstimator.ih"
342 
343 // //
345 
346 #endif // !defined IntegralInvariantVolumeEstimator_h
347 
348 #undef IntegralInvariantVolumeEstimator_RECURSES
349 #endif // else defined(IntegralInvariantVolumeEstimator_RECURSES)
Aim: This class encapsulates its parameter class so that to indicate to the user that the object/poin...
Definition: ConstAlias.h:187
Aim: Smart pointer based on reference counts.
Definition: CountedPtr.h:80
Aim: A wrapper class around a STL associative container for storing sets of digital points within som...
std::pair< KernelConstIterator, KernelConstIterator > PairIterators
Aim: Constructive Solid Geometry (CSG) between models of CEuclideanBoundedShape and CEuclideanOriente...
Aim: Convert a functor on Digital Point to a Functor on Khalimsky Cell.
Aim: A class for computing the Gauss digitization of some Euclidean shape, i.e. its intersection with...
Aim: model of CEuclideanOrientedShape and CEuclideanBoundedShape concepts to create a ball in nD....
Definition: ImplicitBall.h:65
Aim: This class implement an Integral Invariant estimator which computes for each surfel the volume o...
Self & operator=(const Self &other)
CountedConstPtrOrConstPtr< PointPredicate > myPointPredicate
Smart pointer (if required) on a point predicate.
const KernelSpelFunctor myKernelFunctor
Kernel functor (on Spel)
IntegralInvariantVolumeEstimator< TKSpace, TPointPredicate, TVolumeFunctor > Self
VolumeFunctor myFct
The volume functor that transforms the volume into a quantity.
std::vector< DigitalSet * > myKernelsSet
Array of shifting masks. Size = 9 for each shifting (0-adjacent and full kernel included)
BOOST_CONCEPT_ASSERT((concepts::CCellFunctor< ShapeSpelFunctor >))
FunctorOnCells< ShapePointFunctor, KSpace > ShapeSpelFunctor
IntegralInvariantVolumeEstimator(VolumeFunctor fct=VolumeFunctor())
GaussDigitizer< Space, EuclideanMinus > DigitalShape
DigitalSurfaceConvolver< ShapeSpelFunctor, KernelSpelFunctor, KSpace, DigitalShapeKernel > Convolver
functors::ConstValueCell< Value, Spel > KernelSpelFunctor
void setParams(const double dRadius)
BOOST_CONCEPT_ASSERT((concepts::CCellularGridSpaceND< KSpace >))
void attach(ConstAlias< KSpace > K, ConstAlias< PointPredicate > aPointPredicate)
Quantity eval(SurfelConstIterator it) const
OutputIterator eval(SurfelConstIterator itb, SurfelConstIterator ite, OutputIterator result) const
VolumeFunctor::Quantity Quantity
The returned type of the estimator, depends on the functor.
void init(const double _h, SurfelConstIterator itb, SurfelConstIterator ite)
void selfDisplay(std::ostream &out) const
CountedPtr< ShapeSpelFunctor > myShapeSpelFunctor
Smart pointer on functor spel -> {0,1}.
IntegralInvariantVolumeEstimator(const Self &other)
BOOST_STATIC_ASSERT((concepts::ConceptUtils::SameType< typename Convolver::Quantity, typename VolumeFunctor::Argument >::value))
CountedPtr< DigitalShapeKernel > myDigKernel
Digital kernel.
std::vector< PairIterators > myKernels
array of begin/end iterator of shifting masks.
CountedPtr< ShapePointFunctor > myShapePointFunctor
Smart pointer on functor point -> {0,1}.
GaussDigitizer< Space, KernelSupport > DigitalShapeKernel
CountedPtr< KernelSupport > myKernel
Euclidean kernel.
Scalar myRadius
"digital" radius of the kernel (buy may be non integer).
BOOST_CONCEPT_ASSERT((concepts::CPointPredicate< PointPredicate >))
BOOST_CONCEPT_ASSERT((concepts::CUnaryFunctor< VolumeFunctor, Component, Quantity >))
DigitalSetSelector< Domain, BIG_DS+HIGH_VAR_DS >::Type DigitalSet
CountedPtr< Domain > myShapeDomain
Smart pointer on domain
functors::PointFunctorFromPointPredicateAndDomain< PointPredicate, Domain, unsigned int > ShapePointFunctor
IntegralInvariantVolumeEstimator(ConstAlias< KSpace > K, ConstAlias< PointPredicate > aPointPredicate, VolumeFunctor fct=VolumeFunctor())
EuclideanShapesCSG< KernelSupport, KernelSupport > EuclideanMinus
std::set< SCell > SurfelSet
Preferred type for defining a set of surfels (always signed cells).
Aim: implements basic MxN Matrix services (M,N>=1).
Definition: SimpleMatrix.h:76
TComponent Component
Definition: SimpleMatrix.h:79
Aim: Define a simple functor that returns a constant quantity (0 by default).
DGtal is the top-level namespace which contains all DGtal functions and types.
std::ostream & operator<<(std::ostream &out, const ATu0v1< TKSpace, TLinearAlgebra > &object)
Represents a signed cell in a cellular grid space by its Khalimsky coordinates and a boolean value.
Aim: Defines a functor on cells.
Definition: CCellFunctor.h:90
Aim: This concept describes a cellular grid space in nD. In these spaces obtained by cartesian produc...
Aim: Defines a predicate on a point.
Aim: Defines a unary functor, which associates arguments to results.
Definition: CUnaryFunctor.h:90
Create a point functor from a point predicate and a domain.
KSpace K