itkWeightedCovarianceCalculator.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkWeightedCovarianceCalculator.h,v $
  Language:  C++
  Date:      $Date: 2009-03-04 19:29:54 $
  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 __itkWeightedCovarianceCalculator_h
#define __itkWeightedCovarianceCalculator_h

#include "itkArray.h"
#include "itkVariableSizeMatrix.h"
#include "itkSampleAlgorithmBase.h"

namespace itk { 
namespace Statistics {
  
template< class TSample >
class WeightedCovarianceCalculator :
      public SampleAlgorithmBase< TSample >
{
public:
  typedef WeightedCovarianceCalculator   Self;
  typedef SampleAlgorithmBase< TSample > Superclass;
  typedef SmartPointer<Self>             Pointer;
  typedef SmartPointer<const Self>       ConstPointer;

  itkTypeMacro(WeightedCovarianceCalculator, SampleAlgorithmBase);
  itkNewMacro(Self);

  typedef typename Superclass::MeasurementVectorSizeType MeasurementVectorSizeType;

  typedef typename TSample::MeasurementVectorType MeasurementVectorType;

  typedef FunctionBase< MeasurementVectorType, double > WeightFunctionType;

  typedef Array< double >                               MeanType;

  typedef VariableSizeMatrix< double >                  OutputType;

  typedef Array< double > WeightArrayType;

  void SetWeights(WeightArrayType* array);

  WeightArrayType* GetWeights( void );

  void SetWeightFunction(WeightFunctionType* func);

  WeightFunctionType* GetWeightFunction( void );

  void SetMean(MeanType* mean);

  MeanType* GetMean( void );

  OutputType* GetOutput( void );

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

  void GenerateData( void );

  void ComputeCovarianceWithGivenMean( void );
  
  void ComputeCovarianceWithoutGivenMean( void );

private:
  OutputType*         m_Output;
  MeanType*           m_Mean;
  MeanType*           m_InternalMean;
  WeightArrayType*    m_Weights;
  WeightFunctionType* m_WeightFunction;
}; // end of class
    
} // end of namespace Statistics 
} // end of namespace itk 

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

#endif