PCA.hh 3.54 KB
Newer Older
Jan Möbius's avatar
 
Jan Möbius committed
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
//=============================================================================
//
//                               OpenFlipper
//        Copyright (C) 2008 by Computer Graphics Group, RWTH Aachen
//                           www.openflipper.org
//
//-----------------------------------------------------------------------------
//
//                                License
//
//  OpenFlipper is free software: you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or
//  (at your option) any later version.
// 
//  OpenFlipper 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 Lesser General Public License for more details.
// 
//  You should have received a copy of the GNU Lesser General Public License
//  along with OpenFlipper.  If not, see <http://www.gnu.org/licenses/>.
//
//-----------------------------------------------------------------------------
//
//   $Revision$
//   $Author$
//   $Date$
//
//=============================================================================




//=============================================================================
//
//  CLASS PCA
//
//=============================================================================


#ifndef PCA_HH
#define PCA_HH

/*! \file PCA.hh
    \brief Classes for doing Principal component analysis
    
    Classes for doing Principal component analysis
*/

//== INCLUDES =================================================================

Jan Möbius's avatar
Jan Möbius committed
53
#include <gmm/gmm.h>
Jan Möbius's avatar
 
Jan Möbius committed
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121

//== FORWARDDECLARATIONS ======================================================

//== NAMESPACES ===============================================================

/// Namespace for principal Component Analysis
namespace Pca {
 
//== CLASS DEFINITION =========================================================
	      
/** \class PCA PCA.hh </PCA.hh>
    \brief Class for principal component Analysis

    Class for principal component Analysis
*/

template < typename VectorT >
class PCA
{
  typedef gmm::dense_matrix<double>  Matrix;
  typedef std::vector< double >      Vector;

public:
   
  /// Constructor
  PCA() {}
  
  /** Extended constructor
    Parameters : see pca()
  */
  PCA( std::vector< VectorT >& _points, VectorT& _first , VectorT& _second , VectorT& _third);
 
  /// Destructor
  ~PCA();
  
  /**
    Compute center of gravity for a vector of points
    @param _points set of points
  */
  inline VectorT center_of_gravity(const std::vector< VectorT >& _points );
  
  /**
    Compute the principal component analysys for a vector of points
    @param _points set of points
    @param _first  First main axis
    @param _second second main axis
    @param _third  third main axis
  */
  void pca(std::vector< VectorT >& _points , VectorT& _first , VectorT& _second , VectorT& _third);

  bool SymRightEigenproblem( Matrix &_mat_A, Matrix & _mat_VR,
                             std::vector< double > & _vec_EV );
private:

};


//=============================================================================
} //namespace Pca
//=============================================================================
#if defined(INCLUDE_TEMPLATES) && !defined(PCA_C)
#define PCA_TEMPLATES
#include "PCA.cc"
#endif
//=============================================================================
#endif // PCA_HH defined
//=============================================================================