itkSampleToHistogramProjectionFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkSampleToHistogramProjectionFilter.h,v $
  Language:  C++
  Date:      $Date: 2009-03-04 19:29:54 $
  Version:   $Revision: 1.6 $

  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