DGtal  0.9.2
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)
37 
38 #define IntegralInvariantVolumeEstimator_RECURSES
39 
40 #if !defined IntegralInvariantVolumeEstimator_h
41 
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 
155  KSpace, DigitalShapeKernel> Convolver;
158  typedef typename Matrix::Component Component;
159  typedef double Scalar;
163  typename VolumeFunctor::Argument >::value ));
164 
165 
166  // ----------------------- Standard services ------------------------------
167 public:
168 
176  IntegralInvariantVolumeEstimator( VolumeFunctor fct = VolumeFunctor() );
177 
188  ConstAlias< PointPredicate > aPointPredicate,
189  VolumeFunctor fct = VolumeFunctor() );
190 
195 
200  IntegralInvariantVolumeEstimator ( const Self& other );
201 
207  Self& operator= ( const Self& other );
208 
212  void clear();
213 
214  // ----------------------- Interface --------------------------------------
215 public:
216 
218  Scalar h() const;
219 
227  void attach( ConstAlias< KSpace > K,
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 
306  VolumeFunctor myFct;
307  const KernelSpelFunctor myKernelFunctor;
308  std::vector< PairIterators > myKernels;
309  std::vector< DigitalSet * > myKernelsSet;
317  Scalar myH;
318  Scalar myRadius;
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)
void selfDisplay(std::ostream &out) const
GaussDigitizer< Space, KernelSupport > DigitalShapeKernel
std::pair< KernelConstIterator, KernelConstIterator > PairIterators
Aim: Smart pointer based on reference counts.
Definition: CountedPtr.h:79
Aim: This class encapsulates its parameter class so that to indicate to the user that the object/poin...
Definition: ConstAlias.h:186
Aim: Defines a functor on cells.
Definition: CCellFunctor.h:89
CountedPtr< Domain > myShapeDomain
Smart pointer on domain.
functors::PointFunctorFromPointPredicateAndDomain< PointPredicate, Domain, unsigned int > ShapePointFunctor
IntegralInvariantVolumeEstimator(VolumeFunctor fct=VolumeFunctor())
EuclideanShapesCSG< KernelSupport, KernelSupport > EuclideanMinus
FunctorOnCells< ShapePointFunctor, KSpace > ShapeSpelFunctor
GaussDigitizer< Space, EuclideanMinus > DigitalShape
Aim: This class implement an Integral Invariant estimator which computes for each surfel the volume o...
void attach(ConstAlias< KSpace > K, ConstAlias< PointPredicate > aPointPredicate)
VolumeFunctor::Quantity Quantity
The returned type of the estimator, depends on the functor.
Aim: Defines a unary functor, which associates arguments to results.
Definition: CUnaryFunctor.h:89
Aim: implements basic MxN Matrix services (M,N>=1).
Definition: SimpleMatrix.h:75
Aim: Convert a functor on Digital Point to a Functor on Khalimsky Cell.
TComponent Component
Definition: SimpleMatrix.h:79
Self & operator=(const Self &other)
Aim: Constructive Solid Geometry (CSG) between models of CEuclideanBoundedShape and CEuclideanOriente...
Aim: Defines a predicate on a point.
Scalar myRadius
"digital" radius of the kernel (buy may be non integer).
Create a point functor from a point predicate and a domain.
Quantity eval(SurfelConstIterator it) const
Aim: A class for computing the Gauss digitization of some Euclidean shape, i.e. its intersection with...
std::ostream & operator<<(std::ostream &out, const ClosedIntegerHalfPlane< TSpace > &object)
CountedPtr< ShapeSpelFunctor > myShapeSpelFunctor
Smart pointer on functor spel -> {0,1}.
VolumeFunctor myFct
The volume functor that transforms the volume into a quantity.
Aim: model of CEuclideanOrientedShape and CEuclideanBoundedShape concepts to create a ball in nD...
Definition: ImplicitBall.h:64
CountedConstPtrOrConstPtr< PointPredicate > myPointPredicate
Smart pointer (if required) on a point predicate.
DigitalSurfaceConvolver< ShapeSpelFunctor, KernelSpelFunctor, KSpace, DigitalShapeKernel > Convolver
CountedPtr< KernelSupport > myKernel
Euclidean kernel.
DGtal is the top-level namespace which contains all DGtal functions and types.
Aim: A wrapper class around a STL associative container for storing sets of digital points within som...
DigitalSetSelector< Domain, BIG_DS+HIGH_VAR_DS >::Type DigitalSet
BOOST_STATIC_ASSERT((concepts::ConceptUtils::SameType< typename Convolver::Quantity, typename VolumeFunctor::Argument >::value))
const KernelSpelFunctor myKernelFunctor
Kernel functor (on Spel)
Aim: This concept describes a cellular grid space in nD. In these spaces obtained by cartesian produc...
void init(const double _h, SurfelConstIterator itb, SurfelConstIterator ite)
void setParams(const double dRadius)
IntegralInvariantVolumeEstimator< TKSpace, TPointPredicate, TVolumeFunctor > Self
functors::ConstValueCell< Value, Spel > KernelSpelFunctor
BOOST_CONCEPT_ASSERT((concepts::CCellularGridSpaceND< KSpace >))
Aim: Define a simple functor that returns a constant quantity (0 by default).
CountedPtr< DigitalShapeKernel > myDigKernel
Digital kernel.
std::vector< DigitalSet * > myKernelsSet
Array of shifting masks. Size = 9 for each shifting (0-adjacent and full kernel included) ...
std::vector< PairIterators > myKernels
array of begin/end iterator of shifting masks.
CountedPtr< ShapePointFunctor > myShapePointFunctor
Smart pointer on functor point -> {0,1}.