Main Page   Groups   Namespace List   Class Hierarchy   Alphabetical List   Compound List   File List   Namespace Members   Compound Members   File Members   Concepts

itkMeanSquareRegistrationFunction.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkMeanSquareRegistrationFunction.h,v $
00005   Language:  C++
00006   Date:      $Date: 2009-01-24 20:02:59 $
00007   Version:   $Revision: 1.8 $
00008 
00009   Copyright (c) Insight Software Consortium. All rights reserved.
00010   See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
00011 
00012      This software is distributed WITHOUT ANY WARRANTY; without even 
00013      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
00014      PURPOSE.  See the above copyright notices for more information.
00015 
00016 =========================================================================*/
00017 #ifndef __itkMeanSquareRegistrationFunction_h
00018 #define __itkMeanSquareRegistrationFunction_h
00019 
00020 #include "itkPDEDeformableRegistrationFunction.h"
00021 #include "itkPoint.h"
00022 #include "itkCovariantVector.h"
00023 #include "itkInterpolateImageFunction.h"
00024 #include "itkLinearInterpolateImageFunction.h"
00025 #include "itkCentralDifferenceImageFunction.h"
00026 
00027 namespace itk {
00028 
00052 template<class TFixedImage, class TMovingImage, class TDeformationField>
00053 class ITK_EXPORT MeanSquareRegistrationFunction : 
00054   public PDEDeformableRegistrationFunction< TFixedImage,
00055     TMovingImage, TDeformationField>
00056 {
00057 public:
00058 
00060   typedef MeanSquareRegistrationFunction    Self;
00061   typedef PDEDeformableRegistrationFunction< TFixedImage,
00062     TMovingImage, TDeformationField >       Superclass;
00063   typedef SmartPointer<Self>                Pointer;
00064   typedef SmartPointer<const Self>          ConstPointer;
00065 
00067   itkNewMacro(Self);
00068 
00070   itkTypeMacro( MeanSquareRegistrationFunction, 
00071     PDEDeformableRegistrationFunction );
00072 
00074   typedef typename Superclass::MovingImageType     MovingImageType;
00075   typedef typename Superclass::MovingImagePointer  MovingImagePointer;
00076 
00078   typedef typename Superclass::FixedImageType     FixedImageType;
00079   typedef typename Superclass::FixedImagePointer  FixedImagePointer;
00080   typedef typename FixedImageType::IndexType      IndexType;
00081   typedef typename FixedImageType::SizeType       SizeType;
00082   typedef typename FixedImageType::SpacingType    SpacingType;
00083 
00085   typedef typename Superclass::DeformationFieldType    DeformationFieldType;
00086   typedef typename DeformationFieldType::PixelType     DeformationFieldPixelType;
00087   typedef typename Superclass::DeformationFieldTypePointer   
00088     DeformationFieldTypePointer;
00089 
00091   itkStaticConstMacro(ImageDimension, unsigned int,Superclass::ImageDimension);
00092 
00094   typedef typename Superclass::PixelType        PixelType;
00095   typedef typename Superclass::RadiusType       RadiusType;
00096   typedef typename Superclass::NeighborhoodType NeighborhoodType;
00097   typedef typename Superclass::FloatOffsetType  FloatOffsetType;
00098   typedef typename Superclass::TimeStepType     TimeStepType;
00099 
00101   typedef double                                     CoordRepType;
00102   typedef InterpolateImageFunction<MovingImageType,CoordRepType>
00103                                                      InterpolatorType;
00104   typedef typename InterpolatorType::Pointer         InterpolatorPointer;
00105   typedef typename InterpolatorType::PointType       PointType;
00106   typedef LinearInterpolateImageFunction<MovingImageType,CoordRepType>
00107                                                      DefaultInterpolatorType;
00108 
00110   typedef CovariantVector<double,itkGetStaticConstMacro(ImageDimension)> CovariantVectorType;
00111 
00113   typedef CentralDifferenceImageFunction<FixedImageType> GradientCalculatorType;
00114   typedef typename GradientCalculatorType::Pointer   GradientCalculatorPointer;
00115 
00117   void SetMovingImageInterpolator( InterpolatorType * ptr )
00118     { m_MovingImageInterpolator = ptr; }
00119 
00121   InterpolatorType * GetMovingImageInterpolator(void)
00122     { return m_MovingImageInterpolator; }
00123 
00125   virtual TimeStepType ComputeGlobalTimeStep(void *itkNotUsed(GlobalData)) const
00126     { return m_TimeStep; }
00127 
00130   virtual void *GetGlobalDataPointer() const
00131     {
00132     GlobalDataStruct *global = new GlobalDataStruct();
00133     return global;
00134     }
00136 
00138   virtual void ReleaseGlobalDataPointer( void *GlobalData ) const
00139     { delete (GlobalDataStruct *) GlobalData;  }
00140 
00142   virtual void InitializeIteration();
00143 
00146   virtual PixelType  ComputeUpdate(const NeighborhoodType &neighborhood,
00147                      void *globalData,
00148                      const FloatOffsetType &offset = FloatOffsetType(0.0));
00149 
00150 
00151 protected:
00152   MeanSquareRegistrationFunction();
00153   ~MeanSquareRegistrationFunction() {}
00154   void PrintSelf(std::ostream& os, Indent indent) const;
00155 
00157   typedef ConstNeighborhoodIterator<FixedImageType> FixedImageNeighborhoodIteratorType;
00158 
00161   struct GlobalDataStruct
00162     {
00163     FixedImageNeighborhoodIteratorType   m_FixedImageIterator;
00164     };
00165 
00166 private:
00167   MeanSquareRegistrationFunction(const Self&); //purposely not implemented
00168   void operator=(const Self&); //purposely not implemented
00169 
00171   SpacingType                     m_FixedImageSpacing;
00172 
00174   GradientCalculatorPointer       m_FixedImageGradientCalculator;
00175 
00177   InterpolatorPointer             m_MovingImageInterpolator;
00178 
00180   TimeStepType                    m_TimeStep;
00181 
00183   double                          m_DenominatorThreshold;
00184 
00186   double                          m_IntensityDifferenceThreshold;
00187 };
00188 
00189 
00190 } // end namespace itk
00191 
00192 #ifndef ITK_MANUAL_INSTANTIATION
00193 #include "itkMeanSquareRegistrationFunction.txx"
00194 #endif
00195 
00196 #endif
00197 

Generated at Sat Feb 28 12:59:24 2009 for ITK by doxygen 1.5.6 written by Dimitri van Heesch, © 1997-2000