$extrastylesheet
Olena
User documentation 2.1
An Image Processing Platform
Milena
Getting started
API Reference Manual
All Classes
Examples
Demos
Publications
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>
45
struct
closest_point
46
{
47
typedef
algebra::vec<P::dim, float>
input
;
48
typedef
P
result
;
49
50
closest_point
(
const
p_array<P>
&
X
,
const
box<P>
&
box
)
51
:
box_
(box), X(X)
52
,
log_functor_call
(0)
53
54
{ }
55
56
result
57
//inline
58
operator ()
(
const
input
& Ck)
const
59
{
60
++
log_functor_call
;
61
62
algebra::vec<P::dim,float>
Cki = Ck;
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
mln
fun
x2p
closest_point.hh
Copyright (C) 2012 EPITA Research and Development Laboratory (LRDE)