itkCurvatureRegistrationFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkCurvatureRegistrationFilter.h,v $
  Language:  C++
  Date:      $Date: 2008-10-23 16:15:23 $
  Version:   $Revision: 1.4 $

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

#include "itkPDEDeformableRegistrationFilter.h"
#include "itkMeanSquareRegistrationFunction.h"

#if defined(USE_FFTWF) || defined(USE_FFTWD)
#include "fftw3.h"

namespace itk {

template<class TFixedImage, class TMovingImage, class TDeformationField,
  class TImageForceFunction =
  MeanSquareRegistrationFunction<TFixedImage,TMovingImage,TDeformationField> >
class ITK_EXPORT CurvatureRegistrationFilter :
    public PDEDeformableRegistrationFilter< TFixedImage, TMovingImage,
                                            TDeformationField>
{
public:
  typedef CurvatureRegistrationFilter             Self;
  typedef PDEDeformableRegistrationFilter<
    TFixedImage, TMovingImage,TDeformationField>  Superclass;
  typedef SmartPointer<Self>                      Pointer;
  typedef SmartPointer<const Self>                ConstPointer;

  itkNewMacro(Self);

  itkTypeMacro( CurvatureRegistrationFilter, 
                PDEDeformableRegistrationFilter );

  typedef typename Superclass::TimeStepType  TimeStepType;

  typedef typename Superclass::FixedImageType     FixedImageType;
  typedef typename Superclass::FixedImagePointer  FixedImagePointer;
  itkStaticConstMacro(ImageDimension, unsigned int, FixedImageType::ImageDimension );

  typedef typename Superclass::MovingImageType     MovingImageType;
  typedef typename Superclass::MovingImagePointer  MovingImagePointer;

  typedef typename Superclass::DeformationFieldType 
                                                   DeformationFieldType;
  typedef typename Superclass::DeformationFieldPointer  
                                                   DeformationFieldPointer;

  typedef typename TDeformationField::PixelType         DeformationFieldPixelType;
  typedef typename DeformationFieldPixelType::ValueType DeformationFieldComponentType;
  itkStaticConstMacro(DeformationVectorDimension, unsigned int, DeformationFieldPixelType::Dimension );

  #if defined(USE_FFTWD)
  //Prefer to use double precision
  typedef double RealTypeDFT;
  #else
    #if defined(USE_FFTWF)
      //Allow to use single precision
      #warning "Using single precision for FFT computations!"
      typedef double RealTypeDFT;
    #endif
  #endif

  typedef Image<RealTypeDFT,TDeformationField::ImageDimension> DeformationFieldComponentImageType;
  typedef typename DeformationFieldComponentImageType::Pointer DeformationFieldComponentImagePointer;

  typedef typename Superclass::FiniteDifferenceFunctionType
  FiniteDifferenceFunctionType;

  typedef TImageForceFunction RegistrationFunctionType;

  void SetConstraintWeight( const float w ) { m_ConstraintWeight = w; }

  void SetTimeStep( const TimeStepType ts ) { m_TimeStep = ts; }

  virtual double GetMetric() const;

protected:
  CurvatureRegistrationFilter();
  ~CurvatureRegistrationFilter();
  void PrintSelf(std::ostream& os, Indent indent) const;

  virtual void Initialize();

  virtual void ApplyUpdate(TimeStepType dt);

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

  int m_FixedImageDimensions[ImageDimension];

  RealTypeDFT* m_DeformationFieldComponentImage;
  RealTypeDFT* m_DeformationFieldComponentImageDCT;

  float m_ConstraintWeight;

  fftw_plan m_PlanForwardDCT;
  fftw_plan m_PlanBackwardDCT;

  TimeStepType m_TimeStep;

  RealTypeDFT* m_DiagonalElements[ImageDimension];
};


} // end namespace itk

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

#endif //defined(USE_FFTWF) || defined(USE_FFTWD)

#endif