itkGaussianGoodnessOfFitComponent.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkGaussianGoodnessOfFitComponent.h,v $
  Language:  C++
  Date:      $Date: 2009-03-04 15:23:48 $
  Version:   $Revision: 1.10 $

  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