itkMovingHistogramImageFilterBase.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkMovingHistogramImageFilterBase.h,v $
  Language:  C++
  Date:      $Date: 2009-04-23 03:43:42 $
  Version:   $Revision: 1.7 $

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

#include "itkKernelImageFilter.h"
#include <list>
#include <map>
#include <set>

namespace itk {

template<class TInputImage, class TOutputImage, class TKernel >
class ITK_EXPORT MovingHistogramImageFilterBase :
    public KernelImageFilter<TInputImage, TOutputImage, TKernel>
{
public:
  typedef MovingHistogramImageFilterBase                         Self;
  typedef KernelImageFilter<TInputImage, TOutputImage, TKernel>  Superclass;
  typedef SmartPointer<Self>                                     Pointer;
  typedef SmartPointer<const Self>                               ConstPointer;

  itkNewMacro(Self);  

  itkTypeMacro(MovingHistogramImageFilterBase, 
               KernelImageFilter);

  typedef TInputImage                                            InputImageType;
  typedef TOutputImage                                           OutputImageType;
  typedef typename TInputImage::RegionType                       RegionType;
  typedef typename TInputImage::SizeType                         SizeType;
  typedef typename TInputImage::IndexType                        IndexType;
  typedef typename TInputImage::PixelType                        PixelType;
  typedef typename TInputImage::OffsetType                       OffsetType;
  typedef typename Superclass::OutputImageRegionType             OutputImageRegionType;
  typedef typename TOutputImage::PixelType                       OutputPixelType;

  itkStaticConstMacro(ImageDimension, unsigned int,
                      TInputImage::ImageDimension);

  typedef TKernel                                                KernelType;

  typedef typename KernelType::ConstIterator                     KernelIteratorType;

  typedef typename KernelType::SizeType                          RadiusType;

  typedef typename std::list< OffsetType >                       OffsetListType;

  typedef typename std::map< OffsetType, OffsetListType, typename OffsetType::LexicographicCompare >          OffsetMapType;

  void SetKernel( const KernelType& kernel );

  itkGetConstMacro(PixelsPerTranslation, unsigned long);
  
protected:
  MovingHistogramImageFilterBase();
  ~MovingHistogramImageFilterBase() {};

  void PrintSelf(std::ostream& os, Indent indent) const;
  
  void GetDirAndOffset(const IndexType LineStart, 
                      const IndexType PrevLineStart,
                      OffsetType &LineOffset,
                      OffsetType &Changes,
                      int &LineDirection);

  // store the added and removed pixel offset in a list
  OffsetMapType m_AddedOffsets;
  OffsetMapType m_RemovedOffsets;

  // store the offset of the kernel to initialize the histogram
  OffsetListType m_KernelOffsets;

  FixedArray< int, itkGetStaticConstMacro(ImageDimension) > m_Axes;

  unsigned long m_PixelsPerTranslation;


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

  class DirectionCost
    {
    public :
    DirectionCost( int dimension, int count )
      {
      m_Dimension = dimension;
      m_Count = count;
      }
    
    inline bool operator< ( const DirectionCost &dc ) const
      {
      if( m_Count > dc.m_Count )
        { return true; }
      else if( m_Count < dc.m_Count )
        { return false; }
      else //if (m_Count == dc.m_Count) 
        { return m_Dimension > dc.m_Dimension; }
      }

    int m_Dimension;
    int m_Count;
    };

}; // end of class

} // end namespace itk
  
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkMovingHistogramImageFilterBase.txx"
#endif

#endif