itkSampleToHistogramProjectionFilter.h

Go to the documentation of this file.

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

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

  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 __itkSampleToHistogramProjectionFilter_h
#define __itkSampleToHistogramProjectionFilter_h

#include "itkMacro.h"
#include "itkObject.h"
#include "itkSample.h"
#include "itkSubsample.h"
#include "itkHistogram.h"
#include "itkFunctionBase.h"
#include "itkSampleAlgorithmBase.h"
#include "itkArray.h"

namespace itk{ 
namespace Statistics{

template < class TInputSample, class THistogramMeasurement >
class ITK_EXPORT SampleToHistogramProjectionFilter : 
    public SampleAlgorithmBase< TInputSample >
{
public:
  typedef SampleToHistogramProjectionFilter  Self;
  typedef SampleAlgorithmBase< TInputSample > Superclass ;
  typedef SmartPointer< Self > Pointer ;
  typedef const SmartPointer< Self > ConstPointer ;

  itkTypeMacro(SampleToHistogramProjectionFilter, SampleAlgorithmBase);

  itkNewMacro(Self) ;


  typedef typename TInputSample::MeasurementVectorType MeasurementVectorType ;
  typedef typename TInputSample::MeasurementType MeasurementType ;
  typedef typename TInputSample::FrequencyType FrequencyType ;
  typedef typename TInputSample::InstanceIdentifier InstanceIdentifier ;

  typedef typename Superclass::InputSampleType InputSampleType ;
  typedef typename Superclass::MeasurementVectorSizeType MeasurementVectorSizeType;

  typedef   Array< double > ArrayType;

  typedef Array< double > MeanType;

  typedef Histogram< THistogramMeasurement, 1 > HistogramType ;

  virtual void SetHistogram(HistogramType* histogram) ;

  void SetMean(MeanType* center) ;

  MeanType* GetMean() ;

  void SetStandardDeviation(double* value) ;

  double* GetStandardDeviation() ;

  void SetProjectionAxis(ArrayType* axis) ;

  ArrayType* GetProjectionAxis() ;

  void SetHistogramBinOverlap(double overlap) ;

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

  float CalculateOverlap(int binIndex,
                         float dotProduct, 
                         float scale,
                         float marginalDistance,
                         bool firstHalf) ;

  void GenerateData() ;

private:
  SampleToHistogramProjectionFilter(const Self&) ; //purposely not implemented
  void operator=(const Self&) ; //purposely not implemented

  bool m_HistogramUseEquiProbableBins ;
  double m_HistogramBinOverlap ;

  MeanType* m_Mean ;
  ArrayType* m_ProjectionAxis ;
  double* m_StandardDeviation ;
  HistogramType* m_Histogram ;
  FrequencyType m_MinimumFrequency ;
} ; // end of class

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

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

#endif