$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
mahalanobis.hh
1
// Copyright (C) 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_STAT_MAHALANOBIS_HH
27
# define MLN_FUN_STAT_MAHALANOBIS_HH
28
32
33
# include <cmath>
34
# include <mln/core/concept/function.hh>
35
# include <mln/algebra/vec.hh>
36
# include <mln/algebra/mat.hh>
37
38
39
namespace
mln
40
{
41
42
namespace
fun
43
{
44
45
namespace
stat
46
{
47
48
template
<
typename
V>
49
struct
mahalanobis
50
:
public
Function_v2v
< mahalanobis<V> >,
51
private
metal::equal< V, algebra::vec<V::dim,float> >::check_t
52
{
53
enum
{
n
= V::dim };
54
typedef
float
result
;
55
56
mahalanobis
(
const
algebra::mat<n,n,float>
& var,
57
const
algebra::vec<n,float>
&
mean
);
58
// Tech. note: using n (instead of V::dim) above helps g++-2.95.
59
60
float
operator()
(
const
V& v)
const
;
61
62
typedef
algebra::vec<n,float>
mean_t
;
63
64
mean_t
mean
()
const
;
65
66
protected
:
67
algebra::mat<n,n,float>
var_1_
;
68
algebra::vec<n,float>
mean_
;
69
};
70
71
72
# ifndef MLN_INCLUDE_ONLY
73
74
template
<
typename
V>
75
inline
76
mahalanobis<V>::mahalanobis
(
const
algebra::mat<n,n,float>
& var,
77
const
algebra::vec<n,float>
& mean)
78
{
79
var_1_ = var.
_1
();
80
mean_ = mean;
81
}
82
83
template
<
typename
V>
84
inline
85
float
86
mahalanobis<V>::operator()
(
const
V& v)
const
87
{
88
return
std::sqrt
((v - mean_).t() * var_1_ * (v - mean_));
89
}
90
91
template
<
typename
V>
92
inline
93
typename
mahalanobis<V>::mean_t
94
mahalanobis<V>::mean
()
const
95
{
96
return
mean_;
97
}
98
99
# endif // ! MLN_INCLUDE_ONLY
100
101
}
// end of namespace mln::fun::stat
102
103
}
// end of namespace mln::fun
104
105
}
// end of namespace mln
106
107
108
#endif // ! MLN_FUN_STAT_MAHALANOBIS_HH
mln
fun
stat
mahalanobis.hh
Copyright (C) 2012 EPITA Research and Development Laboratory (LRDE)