itkSpatialObjectToImageStatisticsCalculator.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkSpatialObjectToImageStatisticsCalculator.h,v $
  Language:  C++
  Date:      $Date: 2008-10-18 16:11:14 $
  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 __itkSpatialObjectToImageStatisticsCalculator_h
#define __itkSpatialObjectToImageStatisticsCalculator_h

#include "itkObject.h"
#include "itkObjectFactory.h"
#include "itkFloodFilledSpatialFunctionConditionalConstIterator.h"
#include "itkMatrix.h"
#include "itkNumericTraits.h"
#include "itkListSample.h"
#include "itkVector.h"

namespace itk
{

template <class TInputImage, class TInputSpatialObject, unsigned int TSampleDimension=1>
class ITK_EXPORT SpatialObjectToImageStatisticsCalculator : public Object 
{
public:
  typedef SpatialObjectToImageStatisticsCalculator Self;
  typedef Object                                   Superclass;
  typedef SmartPointer<Self>                       Pointer;
  typedef SmartPointer<const Self>                 ConstPointer;

  itkNewMacro(Self);

  itkTypeMacro(SpatialObjectToImageStatisticsCalculator, Object);

  typedef TInputImage                                  ImageType;
  typedef typename TInputImage::Pointer                ImagePointer;
  typedef typename TInputImage::ConstPointer           ImageConstPointer;
  typedef typename TInputImage::PixelType              PixelType;
  typedef typename TInputImage::IndexType              IndexType;

  typedef  typename NumericTraits< PixelType >::AccumulateType AccumulateType;

  itkStaticConstMacro(ImageDimension, unsigned int,
                      ImageType::ImageDimension);
  
  itkStaticConstMacro(SampleDimension, unsigned int,
                      TSampleDimension);

  itkStaticConstMacro(ObjectDimension, unsigned int,
                      TInputSpatialObject::ObjectDimension);


  typedef TInputSpatialObject SpatialObjectType;  
  typedef typename SpatialObjectType::Pointer  SpatialObjectPointer;
  typedef typename SpatialObjectType::ConstPointer SpatialObjectConstPointer;

  typedef itk::FloodFilledSpatialFunctionConditionalConstIterator<ImageType,
                                                 SpatialObjectType> IteratorType;

  typedef Vector< double, TSampleDimension>                    VectorType;
  typedef Matrix< double, TSampleDimension, TSampleDimension > MatrixType;

  typedef itk::Statistics::ListSample< VectorType > SampleType;

  itkSetMacro(SampleDirection,unsigned int);
  itkGetMacro(SampleDirection,unsigned int);

  itkSetConstObjectMacro(Image,ImageType);

  itkSetObjectMacro(SpatialObject,SpatialObjectType);

  const VectorType & GetMean() const {return m_Mean;}

  const MatrixType & GetCovarianceMatrix() const {return m_CovarianceMatrix;}

  AccumulateType GetSum() const {return m_Sum;}

  itkGetConstMacro(NumberOfPixels,unsigned long);

  void Update(void);

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

  bool ComputeStatistics();

private:
  SpatialObjectToImageStatisticsCalculator(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented
  
  ImageConstPointer            m_Image;
  SpatialObjectPointer         m_SpatialObject;
  VectorType                   m_Mean;
  AccumulateType               m_Sum;
  unsigned long                m_NumberOfPixels;
  MatrixType                   m_CovarianceMatrix;
  unsigned int                 m_SampleDirection;
  unsigned long                m_InternalImageTime;
  unsigned long                m_InternalSpatialObjectTime;
  TimeStamp                    m_ModifiedTime;
  typename SampleType::Pointer m_Sample;
};

} // end namespace itk


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

#endif /* __itkSpatialObjectToImageStatisticsCalculator_h */