itkRecursiveSeparableImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkRecursiveSeparableImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2003/02/26 23:20:31 $
  Version:   $Revision: 1.16 $

  Copyright (c) 2002 Insight 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 __itkRecursiveSeparableImageFilter_h
#define __itkRecursiveSeparableImageFilter_h

#include "itkImageToImageFilter.h"
#include "itkNumericTraits.h"

namespace itk
{
  
template <typename TInputImage, typename TOutputImage=TInputImage>
class ITK_EXPORT RecursiveSeparableImageFilter :
    public ImageToImageFilter<TInputImage,TOutputImage> 
{
public:
  typedef RecursiveSeparableImageFilter  Self;
  typedef ImageToImageFilter<TInputImage,TOutputImage>   Superclass;
  typedef SmartPointer<Self>   Pointer;
  typedef SmartPointer<const Self>  ConstPointer;

  itkTypeMacro( RecursiveSeparableImageFilter, ImageToImageFilter );

  typedef typename TInputImage::Pointer  InputImagePointer;
  typedef typename TInputImage::ConstPointer  InputImageConstPointer;

  typedef typename TInputImage::PixelType                   InputPixelType;
  typedef typename NumericTraits<InputPixelType>::RealType  RealType;

  typedef TInputImage      InputImageType;

  typedef TOutputImage      OutputImageType;

  itkGetMacro(Direction, unsigned int);

  itkSetMacro(Direction, unsigned int);

  void SetInputImage( const TInputImage * );
    
  const TInputImage * GetInputImage( void );

  virtual void GenerateInputRequestedRegion() throw(InvalidRequestedRegionError);

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

  void GenerateData(void);

  // Override since the filter produces the entire dataset
  void EnlargeOutputRequestedRegion(DataObject *output);

  virtual void SetUp(void) = 0;

  void ApplyRecursiveFilter(unsigned int dimension);

  virtual void ComputeFilterCoefficients(bool symmetric) = 0;

  void FilterDataArray(RealType *outs,
                 const RealType *data, unsigned int ln);

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

  unsigned int m_Direction;

protected:
  RealType m_K;                       

  RealType m_Spacing;

  RealType m_A0;
  RealType m_A1;
  RealType m_B0;
  RealType m_B1;
  RealType m_C0;
  RealType m_C1;
  RealType m_W0;
  RealType m_W1; 

  RealType m_N00;
  RealType m_N11;
  RealType m_N22;
  RealType m_N33; 
  
  RealType m_D11;
  RealType m_D22;
  RealType m_D33;
  RealType m_D44; 
  
  RealType m_M11;
  RealType m_M22;
  RealType m_M33;
  RealType m_M44; 

  RealType m_BN1;
  RealType m_BN2;
  RealType m_BN3;
  RealType m_BN4; 
 
  RealType m_BM1;
  RealType m_BM2;
  RealType m_BM3;
  RealType m_BM4; 
 
};


} // end namespace itk

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


#endif

---