DGtal  0.9.2
TrueLocalEstimatorOnPoints.ih
1 /**
2  * This program is free software: you can redistribute it and/or modify
3  * it under the terms of the GNU Lesser General Public License as
4  * published by the Free Software Foundation, either version 3 of the
5  * License, or (at your option) any later version.
6  *
7  * This program is distributed in the hope that it will be useful,
8  * but WITHOUT ANY WARRANTY; without even the implied warranty of
9  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
10  * GNU General Public License for more details.
11  *
12  * You should have received a copy of the GNU General Public License
13  * along with this program. If not, see <http://www.gnu.org/licenses/>.
14  *
15  **/
16 
17 /**
18  * @file TrueLocalEstimatorOnPoints.ih
19  * @author David Coeurjolly (\c david.coeurjolly@liris.cnrs.fr )
20  * Laboratoire d'InfoRmatique en Image et Systèmes d'information - LIRIS (CNRS, UMR 5205), CNRS, France
21  *
22  * @date 2011/06/27
23  *
24  * Implementation of inline methods defined in TrueLocalEstimatorOnPoints.h
25  *
26  * This file is part of the DGtal library.
27  */
28 
29 
30 //////////////////////////////////////////////////////////////////////////////
31 #include <cstdlib>
32 //////////////////////////////////////////////////////////////////////////////
33 
34 ///////////////////////////////////////////////////////////////////////////////
35 // IMPLEMENTATION of inline methods.
36 ///////////////////////////////////////////////////////////////////////////////
37 
38 ///////////////////////////////////////////////////////////////////////////////
39 // ----------------------- Standard services ------------------------------
40 
41 
42 // ------------------------------------------------------------------------
43 template <typename CIt, typename PShape, typename PShapeFunctor>
44 inline
45 DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::TrueLocalEstimatorOnPoints()
46  :myH(0.0), myBegin(), myEnd(), myFunctorPtr(0)
47 {
48 }
49 
50 // ------------------------------------------------------------------------
51 template <typename CIt, typename PShape, typename PShapeFunctor>
52 inline
53 bool
54 DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::isValid() const
55 {
56  return ( (myH > 0)&&(myFunctorPtr != 0) );
57 }
58 
59 // ------------------------------------------------------------------------
60 template <typename CIt, typename PShape, typename PShapeFunctor>
61 inline
62 void
63 DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>
64 ::init(const double h,
65  const ConstIterator& itb,
66  const ConstIterator& ite)
67 {
68  myH = h;
69  ASSERT(myH > 0);
70  myBegin = itb;
71  myEnd = ite;
72 }
73 
74 // ------------------------------------------------------------------------
75 template <typename CIt, typename PShape, typename PShapeFunctor>
76 inline
77 void
78 DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>
79 ::attach(ParametricShape* aShapePtr)
80 {
81  myFunctorPtr = new ParametricShapeFunctor(aShapePtr);
82 }
83 
84 // ------------------------------------------------------------------------
85 template <typename CIt, typename PShape, typename PShapeFunctor>
86 inline
87 DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::~TrueLocalEstimatorOnPoints()
88 {
89  delete myFunctorPtr;
90 }
91 
92 // ------------------------------------------------------------------------
93 template <typename CIt, typename PShape, typename PShapeFunctor>
94 inline
95 typename DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::Quantity
96 DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>
97 ::eval(const ConstIterator& it) const
98 {
99  ASSERT( isValid() );
100  RealPoint p( *it );
101  p *= myH;
102  return myFunctorPtr->operator()(p);
103 }
104 
105 // ------------------------------------------------------------------------
106 template <typename CIt, typename PShape, typename PShapeFunctor>
107 template <typename OutputIterator>
108 inline
109 OutputIterator
110 DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>
111 ::eval(const ConstIterator& itb,
112  const ConstIterator& ite,
113  OutputIterator result) const
114 {
115  ASSERT( isValid() );
116  // do-while loop to deal with the case of a whole circular range
117  if (isNotEmpty(itb, ite))
118  {
119  ConstIterator it = itb;
120  do
121  {
122  *result++ = eval( it );
123  ++it;
124  } while (it != ite);
125  }
126  return result;
127 }