summaryrefslogtreecommitdiff
path: root/2005/i/robert/src/ovision/see/hotelling.hh
blob: b3c4081aac785a82644208caa0019ebdf68677e2 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
// robert - programme du robot 2005 
// 
// Copyright (C) 2005 Olivier Gaillard
//
// Robot APB Team/Efrei 2005.
//        Web: http://assos.efrei.fr/robot/
//      Email: robot AT efrei DOT fr
//
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
// 
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.

#ifndef hotelling_h
#define hotelling_h

#include <vector>

#include "oconfig.hh"
#include "hpoint.hh"

/// Detecte l'orientation principale d'une masse de donn�es
class Hotelling
{
    /// Classe oconfig
    OConfig *oconfig_;
    /// Matrice de covariance
    double covMatrix_[2][2];
    /// Donn�es � trait�es ajust�es
    Hpoint *data_;
    /// Nombre d'�chantillons de donn�es
    unsigned dataSize_;
    /// Moyenne de la premi�re et deuxi�me composante
    double meanX_, meanY_;
    /// EigenValues
    double val1_, val2_; 
    /// EigenVectors
    Hpoint e1_, e2_;

  public:
    /// Constructeur
    Hotelling (const std::vector<Hpoint> &dataList);
    /// Destructeur
    ~Hotelling (void);
    /// Affiche la matrice de covariance
    void showCovMatrix () const;
    /// Affiche les eigenvectors
    void show () const;
    /// Calcul des eigenvectors
    void eigenVectors ();
    /// Accessors
    void getPC (Hpoint &pca);
    void getPC (double &x, double &y);
    void get (double &x, double &y, double &x2, double &y2);

  private:
    /// Calcul des eigenvalues
    void eigenValues ();
    /// Calcul de la variance
    double cov (const unsigned i, const unsigned j);
    /// Calcul de la matrice de covariance
    void covMatrix ();
    /// Calcul de la moyenne
    void mean (const std::vector<Hpoint> &dataList);
};


#endif // hotelling_h