Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkNCCRegistrationFunction_h
00018 #define __itkNCCRegistrationFunction_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 NCCRegistrationFunction :
00054 public PDEDeformableRegistrationFunction< TFixedImage,
00055 TMovingImage, TDeformationField>
00056 {
00057 public:
00058
00060 typedef NCCRegistrationFunction 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( NCCRegistrationFunction,
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 Superclass::DeformationFieldTypePointer
00087 DeformationFieldTypePointer;
00088
00090 itkStaticConstMacro(ImageDimension, unsigned int,Superclass::ImageDimension);
00091
00093 typedef typename Superclass::PixelType PixelType;
00094 typedef typename Superclass::RadiusType RadiusType;
00095 typedef typename Superclass::NeighborhoodType NeighborhoodType;
00096
00097 typedef typename Superclass::FloatOffsetType FloatOffsetType;
00098 typedef typename Superclass::TimeStepType TimeStepType;
00099
00101 typedef double CoordRepType;
00102 typedef InterpolateImageFunction<MovingImageType,CoordRepType> InterpolatorType;
00103
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 protected:
00151 NCCRegistrationFunction();
00152 ~NCCRegistrationFunction() {}
00153 void PrintSelf(std::ostream& os, Indent indent) const;
00154
00156 typedef ConstNeighborhoodIterator<FixedImageType> FixedImageNeighborhoodIteratorType;
00157
00160 struct GlobalDataStruct
00161 {
00162 FixedImageNeighborhoodIteratorType m_FixedImageIterator;
00163 };
00164
00165 private:
00166 NCCRegistrationFunction(const Self&);
00167 void operator=(const Self&);
00168
00170 SpacingType m_FixedImageSpacing;
00171
00173 GradientCalculatorPointer m_FixedImageGradientCalculator;
00174
00175
00177 InterpolatorPointer m_MovingImageInterpolator;
00178
00180 TimeStepType m_TimeStep;
00181
00183 double m_DenominatorThreshold;
00184
00186 double m_IntensityDifferenceThreshold;
00187
00188 mutable double m_MetricTotal;
00189
00190 };
00191
00192
00193 }
00194
00195 #ifndef ITK_MANUAL_INSTANTIATION
00196 #include "itkNCCRegistrationFunction.txx"
00197 #endif
00198
00199 #endif
00200