Numerics/Statistics/itkHistogramToImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkHistogramToImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2009-04-06 11:51:06 $
  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 __itkHistogramToImageFilter_h
#define __itkHistogramToImageFilter_h

#include "itkImageSource.h"
#include "itkConceptChecking.h"
#include "itkHistogram.h"
#include "itkImage.h"
#include "itkImageRegionIteratorWithIndex.h"
#include "itkSimpleDataObjectDecorator.h"


namespace itk
{

template <class THistogram, class TFunction>
class ITK_EXPORT HistogramToImageFilter :
  public ImageSource<Image< typename TFunction::OutputPixelType , 
  ::itk::Statistics::GetHistogramDimension<THistogram>::HistogramDimension> >
{
public:

  typedef TFunction                             FunctorType;
  typedef typename FunctorType::OutputPixelType OutputPixelType;
  typedef HistogramToImageFilter                Self;
  typedef ImageSource< Image<OutputPixelType, 
    ::itk::Statistics::GetHistogramDimension<THistogram>::HistogramDimension> >
    Superclass;
  typedef SmartPointer<Self>                   Pointer;
  typedef SmartPointer<const Self>             ConstPointer;

  typedef Image<OutputPixelType, 
    ::itk::Statistics::GetHistogramDimension<THistogram>::HistogramDimension >
    OutputImageType; 
  typedef typename Superclass::Pointer    OutputImagePointer;
  typedef typename OutputImageType::SpacingType SpacingType;
  typedef typename OutputImageType::PointType   PointType;

  
  // Define an iterator to iterate through the image
  typedef itk::ImageRegionIteratorWithIndex< Image<OutputPixelType, 
    ::itk::Statistics::GetHistogramDimension<THistogram>::HistogramDimension> >
                                                             ImageIteratorType;

  itkNewMacro(Self);

  itkTypeMacro(HistogramToImageFilter,ImageSource);

  typedef typename Superclass::OutputImageRegionType OutputImageRegionType;

  typedef THistogram                                     HistogramType;
  typedef typename HistogramType::MeasurementVectorType  MeasurementVectorType;
  typedef typename HistogramType::SizeType               HistogramSizeType;
  typedef typename HistogramType::SizeType               SizeType;

  typedef SimpleDataObjectDecorator< HistogramType* > InputHistogramObjectType;

  itkStaticConstMacro(ImageDimension, unsigned int,
    ::itk::Statistics::GetHistogramDimension<THistogram>::HistogramDimension );

  virtual void SetInput( const HistogramType *histogram);
  virtual void SetInput( const InputHistogramObjectType *inputObject);
  const InputHistogramObjectType* GetInput(void);

                      
  itkSetMacro(Spacing,SpacingType);
  virtual void SetSpacing(const double* values);


  itkGetConstReferenceMacro(Spacing,SpacingType);

  itkSetMacro(Origin,PointType);
  virtual void SetOrigin(const double* values);

  itkGetConstReferenceMacro(Origin,PointType);

  itkGetMacro(Size,HistogramSizeType);


  void SetFunctor(const FunctorType& functor)
    {
    m_Functor = functor;
    this->Modified();
    }

  FunctorType& GetFunctor() { return m_Functor; }
  const FunctorType& GetFunctor() const { return m_Functor; }

  void SetTotalFrequency( unsigned long n );
  

  
protected:
  HistogramToImageFilter();
  ~HistogramToImageFilter();

  virtual void GenerateOutputInformation();
  virtual void GenerateData();
  
  FunctorType m_Functor;
  
  SizeType        m_Size;
  SpacingType     m_Spacing;
  PointType       m_Origin;

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

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


};

} // end namespace itk

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

#endif