$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
x2v/linear.hh
1 // Copyright (C) 2008, 2009 EPITA Research and Development Laboratory (LRDE)
2 //
3 // This file is part of Olena.
4 //
5 // Olena is free software: you can redistribute it and/or modify it under
6 // the terms of the GNU General Public License as published by the Free
7 // Software Foundation, version 2 of the License.
8 //
9 // Olena is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // General Public License for more details.
13 //
14 // You should have received a copy of the GNU General Public License
15 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
16 //
17 // As a special exception, you may use this file as part of a free
18 // software project without restriction. Specifically, if other files
19 // instantiate templates or use macros or inline functions from this
20 // file, or you compile this file and link it with other files to produce
21 // an executable, this file does not by itself cause the resulting
22 // executable to be covered by the GNU General Public License. This
23 // exception does not however invalidate any other reasons why the
24 // executable file might be covered by the GNU General Public License.
25 
26 #ifndef MLN_FUN_X2V_LINEAR_HH
27 # define MLN_FUN_X2V_LINEAR_HH
28 
34 # include <mln/core/image/image1d.hh>
35 # include <mln/core/concept/function.hh>
36 # include <mln/fun/internal/selector.hh>
37 # include <mln/convert/to.hh>
38 # include <mln/algebra/vec.hh>
39 
40 namespace mln
41 {
42 
43  namespace fun
44  {
45 
46  namespace x2x
47  {
48 
52  template < typename I >
53  struct linear
54  : public fun::internal::selector_<const algebra::vec<1,float>,
55  // float is a dummy parameter (real is C)
56  mln_value(I), linear<I> >::ret
57  {
58  typedef mln_value(I) result;
59 
61  linear(const I& ima);
62 
65  template <typename C>
66  mln_value(I)
67  operator()(const algebra::vec<1,C>& v) const;
68 
70  const I& ima;
71  };
72 
73 
74 # ifndef MLN_INCLUDE_ONLY
75 
76  template <typename I>
77  linear<I>::linear(const I& ima) : ima(ima)
78  {
79  mlc_bool(I::psite::dim == 1)::check();
80  }
81 
82  template <typename I>
83  template <typename C>
84  mln_value(I)
85  linear<I>::operator()(const algebra::vec<1,C>& v) const
86  {
87  typedef mln_sum(mln_value(I)) vsum;
88 
89  // looking for img(x);
90  double x = v[0];
91 
92  // p1
93  double xa = mln_point(I)::coord(v[0]);
94  vsum ya = ima(point1d(xa));
95 
96  // x makes sens in img
97  if (x == xa)
98  return ima(xa);
99 
100  // p2
101  double xb = mln_point(I)::coord(v[0] + 1);
102  vsum yb = ima(point1d(xb));
103 
104  // Taylor-young
105  return convert::to<mln_value(I)>
106  (ya + (x - xa) * (yb - ya) / (xb - xa));
107  }
108 
109 # endif // ! MLN_INCLUDE_ONLY
110 
111  } // end of namespace mln::fun::x2x
112 
113  } // end of namespace mln::fun
114 
115 } // end of namespace mln
116 
117 
118 #endif // ! MLN_FUN_X2V_LINEAR_HH