itkComplexBSplineInterpolateImageFunction.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkComplexBSplineInterpolateImageFunction.h,v $
  Language:  C++
  Date:      $Date: 2009-05-15 20:10:02 $
  Version:   $Revision: 1.2 $

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

#include "itkInterpolateImageFunction.h"
#include "itkBSplineInterpolateImageFunction.h"
#include "itkComplexToRealImageFilter.h"
#include "itkComplexToImaginaryImageFilter.h"

namespace itk {
 template< class TImageType, class TCoordRep = double, class TCoefficientType = double >
 class ITK_EXPORT ComplexBSplineInterpolateImageFunction : 
  public InterpolateImageFunction< TImageType, TCoordRep > 
{
public:

  typedef ComplexBSplineInterpolateImageFunction       Self;

  typedef InterpolateImageFunction< TImageType, TCoordRep >  Superclass;

  typedef SmartPointer< Self >                    Pointer;

  typedef SmartPointer< const Self >              ConstPointer;

    itkTypeMacro( ComplexBSplineInterpolateImageFunction, InterpolateImageFunction );

 
  itkNewMacro( Self );

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

  typedef typename Superclass::OutputType OutputType;

  typedef typename Superclass::InputImageType InputImageType;

  typedef typename Superclass::IndexType IndexType;

  typedef typename Superclass::ContinuousIndexType ContinuousIndexType;

  typedef typename Superclass::PointType PointType;

  typedef Image< double, itkGetStaticConstMacro( ImageDimension ) > InternalImageType;

  typedef ComplexToRealImageFilter< InputImageType, InternalImageType > RealFilterType;
  typedef ComplexToImaginaryImageFilter< InputImageType, InternalImageType > ImaginaryFilterType;

  typedef BSplineInterpolateImageFunction< InternalImageType, TCoordRep, TCoefficientType > InterpolatorType;


  
  virtual OutputType EvaluateAtContinuousIndex( const ContinuousIndexType & index ) const; 

/*  typedef CovariantVector< OutputType, itkGetStaticConstMacro( ImageDimension ) > CovariantVectorType;

  CovariantVectorType EvaluateDerivative( const PointType & point ) const
  {
   ContinuousIndexType index;
   this->GetInputImage()->TransformPhysicalPointToContinuousIndex( point, index );
   return ( this->EvaluateDerivativeAtContinuousIndex( index ) );
  } 

  CovariantVectorType EvaluateDerivativeAtContinuousIndex( const ContinuousIndexType & x ) const;
*/

  void SetSplineOrder( unsigned int SplineOrder );
  itkGetConstMacro( SplineOrder, int );

  virtual void SetInputImage(const TImageType * inputData);

  
protected:
  ComplexBSplineInterpolateImageFunction();
  virtual ~ComplexBSplineInterpolateImageFunction() {};
  void operator=( const Self& ); //purposely not implemented
  void PrintSelf(std::ostream& os, Indent indent) const;
  
  
private:
  ComplexBSplineInterpolateImageFunction( const Self& ); //purposely not implemented
  
  
  unsigned int                          m_SplineOrder;
  typename InterpolatorType::Pointer    m_RealInterpolator;
  typename InterpolatorType::Pointer    m_ImaginaryInterpolator;
  typename RealFilterType::Pointer      m_RealFilter;
  typename ImaginaryFilterType::Pointer m_ImaginaryFilter;
  
}; // class


} // namespace


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

#endif