$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
closest_point.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_X2P_CLOSEST_POINT_HH
27 # define MLN_FUN_X2P_CLOSEST_POINT_HH
28 
29 # include <mln/algebra/vec.hh>
30 # include <mln/norm/l2.hh>
31 # include <mln/core/site_set/p_array.hh>
32 # include <mln/core/site_set/box.hh>
33 
34 namespace mln
35 {
36 
37  namespace fun
38  {
39 
40  namespace x2p
41  {
42 
44  template <typename P>
46  {
48  typedef P result;
49 
51  : box_(box), X(X)
52  , log_functor_call(0)
53 
54  { }
55 
56  result
57  //inline
58  operator () (const input& Ck) const
59  {
61 
63  algebra::vec<P::dim,float> best_x = convert::to< algebra::vec<P::dim,float> >(X[0]);
64  float best_d = norm::l2(Cki - best_x);
65  for (unsigned j = 1; j < X.nsites(); ++j)
66  {
67  algebra::vec<P::dim,float> Xj = convert::to< algebra::vec<P::dim,float> >(X[j]);
68  float d = norm::l2(Cki - Xj);
69  if (d < best_d)
70  {
71  best_d = d;
72  best_x = Xj;
73  }
74  }
75  return convert::to<P>(best_x);
76  }
77 
78  const box<P>& domain() const
79  {
80  return box_;
81  }
82 
83  const box<P> box_;
84  const p_array<P>& X;
85 
86  // log call to the functor
87  mutable unsigned log_functor_call;
88  };
89 
90  } // end of namespace mln::fun::x2p
91 
92  } // end of namespace mln::fun
93 
94 } // end of namespace mln
95 
96 #endif // ! MLN_FUN_X2P_CLOSEST_POINT_HH