PCA.hpp

Go to the documentation of this file.

//
// This file is part of MorphoGraphX - http://www.MorphoGraphX.org
// Copyright (C) 2012-2015 Richard S. Smith and collaborators.
//
// If you use MorphoGraphX in your work, please cite:
//   http://dx.doi.org/10.7554/eLife.05864
//
// MorphoGraphX is free software, and is licensed under under the terms of the 
// GNU General (GPL) Public License version 2.0, http://www.gnu.org/licenses.
// 
#ifndef PCA_HPP
#define PCA_HPP
 
#include <Process.hpp>
 
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_vector.h>
#include <gsl/gsl_eigen.h>
 
namespace mgx 
{
 
  // Class to calculate principle directions of cell to find plane of projection
  class mgx_EXPORT PCA
  {
 
    public:
 
    Point3f p1, p2, p3;
    Point3f ev;
    gsl_matrix *mat;
    gsl_vector *eval;
    gsl_matrix *evec;
    gsl_eigen_symmv_workspace *w;
 
    PCA() : mat(gsl_matrix_alloc(3, 3)), eval(gsl_vector_alloc(3)),
            evec(gsl_matrix_alloc(3, 3)), w(gsl_eigen_symmv_alloc(3)) {} 
 
    ~PCA()
    {
      gsl_eigen_symmv_free(w);
      gsl_vector_free(eval);
      gsl_matrix_free(evec);
      gsl_matrix_free(mat);
    }
 
    bool operator()(const Matrix3f& corr);
    
  };
 
 
}
#endif