itkGaussianGoodnessOfFitComponent.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkGaussianGoodnessOfFitComponent.h,v $
  Language:  C++
  Date:      $Date: 2005-07-26 15:54:55 $
  Version:   $Revision: 1.9 $

  Copyright (c) Insight Software 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 __itkGaussianGoodnessOfFitComponent_h
#define __itkGaussianGoodnessOfFitComponent_h

#include "itkGoodnessOfFitComponentBase.h"
#include "itkGaussianDensityFunction.h"
#include "itkFunctionBase.h"
#include "itkWeightedCovarianceCalculator.h"
#include "itkSymmetricEigenAnalysis.h"

namespace itk{ 
namespace Statistics{

template< class TInputSample >
class GaussianGoodnessOfFitComponent 
  : public GoodnessOfFitComponentBase< TInputSample > 
{
public:
  typedef GaussianGoodnessOfFitComponent Self;
  typedef GoodnessOfFitComponentBase< TInputSample > Superclass;
  typedef SmartPointer< Self > Pointer;
  typedef SmartPointer< const Self > ConstPointer;

  itkTypeMacro(GaussianGoodnessOfFitComponent, 
               GoodnessOfFitComponentBase) ;

  itkNewMacro(Self) ;

  
  typedef typename TInputSample::MeasurementType MeasurementType ;
  typedef typename TInputSample::MeasurementVectorType MeasurementVectorType ;

  typedef typename Superclass::CenterType CenterType ;
  typedef typename Superclass::RadiusType RadiusType ;
  typedef typename Superclass::MeanType MeanType ;
  typedef typename Superclass::StandardDeviationType StandardDeviationType ;
  typedef typename Superclass::ResampledSampleType ResampledSampleType ;
  typedef typename Superclass::ProjectionAxisArrayType ProjectionAxisArrayType;
  typedef typename Superclass::MeasurementVectorSizeType MeasurementVectorSizeType;
  typedef Array< double > ParametersType ;

  typedef GaussianDensityFunction< MeasurementVectorType > 
  ProbabilityDensityFunctionType ;

  typedef typename ProbabilityDensityFunctionType::CovarianceType CovarianceType ;

  typedef WeightedCovarianceCalculator< ResampledSampleType > 
  CovarianceCalculatorType ;

  typedef Array< double > EigenValuesArrayType;
  typedef SymmetricEigenAnalysis< ProjectionAxisArrayType, EigenValuesArrayType >
          ProjectionAxisCalculatorType;

  unsigned int GetNumberOfParameters() const
  { return (unsigned int)(this->GetMeasurementVectorSize() + 1) ; }

  void SetParameters(const ParametersType &parameter) ;

  CenterType* GetCenter() ;

  RadiusType* GetRadius() ;

  MeanType* GetMean() ;

  StandardDeviationType* GetStandardDeviation() ;

  double GetCumulativeProbability(double x) const ;

  double GetProbabilityDensity(MeasurementVectorType &measurements) const ;

  void PrintParameters(std::ostream &os) const ;

  ParametersType GetFullParameters() const ;

  virtual void SetInputSample( const TInputSample* sample );

protected:
  GaussianGoodnessOfFitComponent() ;
  virtual ~GaussianGoodnessOfFitComponent() ;
  virtual void PrintSelf(std::ostream& os, Indent indent) const ;

  virtual void CalculateProjectionAxes() ;

private:
  typename ProbabilityDensityFunctionType::Pointer 
  m_ProbabilityDensityFunction ;
  typename CovarianceCalculatorType::Pointer m_CovarianceCalculator ;
  ProjectionAxisCalculatorType * m_ProjectionAxisCalculator ;

  MeanType m_Mean ;
  CenterType m_Center ;
  RadiusType m_Radius ;
  StandardDeviationType m_StandardDeviation ;
  CovarianceType m_Covariance ;

  unsigned int m_NumberOfParameters ;

  int m_LongestAxisIndex ;
  double m_LargestEigenValue ;
} ; // end of class

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

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

#endif