itkMahalanobisDistanceMembershipFunction.h

Go to the documentation of this file.

/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkMahalanobisDistanceMembershipFunction.h,v $
  Language:  C++
  Date:      $Date: 2002/11/13 00:52:54 $
  Version:   $Revision: 1.3 $

  Copyright (c) 2002 Insight Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#ifndef __itkMahalanobisDistanceMembershipFunction_h
#define __itkMahalanobisDistanceMembershipFunction_h

#include <vnl/vnl_vector.h>
#include <vnl/vnl_transpose.h>
#include <vnl/vnl_matrix.h>
#include <vnl/algo/vnl_matrix_inverse.h>
#include <vnl/algo/vnl_determinant.h>


#include "itkMembershipFunctionBase.h"

#define PI 3.141592

namespace itk{ 
  namespace Statistics{

template< class TVector >
class ITK_EXPORT MahalanobisDistanceMembershipFunction :
      public MembershipFunctionBase< TVector >
{
public:
  typedef MahalanobisDistanceMembershipFunction Self;
  typedef MembershipFunctionBase< TVector > Superclass ;
  typedef SmartPointer<Self> Pointer;

  itkTypeMacro(MahalanobisDistanceMembershipFunction, MembershipFunctionBase);
  itkNewMacro(Self);
  
  typedef TVector MeasurementVectorType ;

  void SetMean(const vnl_vector<double> &mean) ;
  
  const vnl_vector<double> & GetMean() const;

  void SetCovariance(const vnl_matrix<double> &cov); 
  
  const vnl_matrix<double> & GetCovariance() const;

  void SetInverseCovariance(const vnl_matrix<double> &invcov); 
  
  const vnl_matrix<double> & GetInverseCovariance() const;

  itkSetMacro( NumberOfSamples, double );

  itkGetMacro( NumberOfSamples, double );

  double Evaluate(const MeasurementVectorType &measurement) const;
  
protected:
  MahalanobisDistanceMembershipFunction(void) ;
  virtual ~MahalanobisDistanceMembershipFunction(void) {}
  void PrintSelf(std::ostream& os, Indent indent) const;

private:
  vnl_vector< double >  m_Mean;              // mean
  vnl_matrix< double >  m_Covariance;         // covariance matrix

  // inverse covariance matrix which is automatically calculated 
  // when covariace matirx is set.  This speed up the GetProbability()
  vnl_matrix< double >  m_InverseCovariance; 

  // Number of samples defining this density
  double m_NumberOfSamples;
  // pre_factor which is automatically calculated 
  // when covariace matirx is set.  This speeds up the GetProbability()  
  double       m_PreFactor;
  double       m_Epsilon;
  double       m_DoubleMax;
 
  bool         m_ValidInverseCovarianceFlag;

  itkStaticConstMacro(VectorDimension, unsigned int, TVector::Dimension);
  typedef vnl_matrix_fixed<double,1,itkGetStaticConstMacro(VectorDimension)> ColumnVectorType;

  mutable ColumnVectorType      m_TempVec;
  mutable ColumnVectorType      m_TempMat;

  void CalculateInverseCovariance(); 
};

  } // end of namespace Statistics
} // end namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#include "itkMahalanobisDistanceMembershipFunction.txx"
#endif

#endif

---