$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 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 
58  // Tech. note: using n (instead of V::dim) above helps g++-2.95.
59 
60  float operator()(const V& v) const;
61 
63 
64  mean_t mean() const;
65 
66  protected:
69  };
70 
71 
72 # ifndef MLN_INCLUDE_ONLY
73 
74  template <typename V>
75  inline
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