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:
00059
typedef NCCRegistrationFunction
Self;
00060 typedef PDEDeformableRegistrationFunction< TFixedImage,
00061 TMovingImage, TDeformationField >
Superclass;
00062 typedef SmartPointer<Self> Pointer;
00063 typedef SmartPointer<const Self> ConstPointer;
00064
00066
itkNewMacro(
Self);
00067
00069
itkTypeMacro( NCCRegistrationFunction,
00070
PDEDeformableRegistrationFunction );
00071
00073
typedef typename Superclass::MovingImageType
MovingImageType;
00074 typedef typename Superclass::MovingImagePointer
MovingImagePointer;
00075
00077
typedef typename Superclass::FixedImageType
FixedImageType;
00078 typedef typename Superclass::FixedImagePointer
FixedImagePointer;
00079 typedef typename FixedImageType::IndexType
IndexType;
00080 typedef typename FixedImageType::SizeType
SizeType;
00081
00083
typedef typename Superclass::DeformationFieldType
DeformationFieldType;
00084 typedef typename Superclass::DeformationFieldTypePointer
00085
DeformationFieldTypePointer;
00086
00088
itkStaticConstMacro(ImageDimension,
unsigned int,Superclass::ImageDimension);
00089
00091
typedef typename Superclass::PixelType
PixelType;
00092 typedef typename Superclass::RadiusType
RadiusType;
00093 typedef typename Superclass::NeighborhoodType
NeighborhoodType;
00094
00095
typedef typename Superclass::FloatOffsetType
FloatOffsetType;
00096 typedef typename Superclass::TimeStepType
TimeStepType;
00097
00099
typedef double CoordRepType;
00100 typedef InterpolateImageFunction<MovingImageType,CoordRepType> InterpolatorType;
00101 typedef typename InterpolatorType::Pointer
InterpolatorPointer;
00102 typedef typename InterpolatorType::PointType
PointType;
00103 typedef LinearInterpolateImageFunction<MovingImageType,CoordRepType>
00104
DefaultInterpolatorType;
00105
00107
typedef CovariantVector<double,itkGetStaticConstMacro(ImageDimension)> CovariantVectorType;
00108
00110
typedef CentralDifferenceImageFunction<FixedImageType> GradientCalculatorType;
00111 typedef typename GradientCalculatorType::Pointer
GradientCalculatorPointer;
00112
00114
void SetMovingImageInterpolator(
InterpolatorType * ptr )
00115 { m_MovingImageInterpolator = ptr; }
00116
00118 InterpolatorType * GetMovingImageInterpolator(
void)
00119 {
return m_MovingImageInterpolator; }
00120
00122
virtual TimeStepType ComputeGlobalTimeStep(
void *
itkNotUsed(GlobalData))
const
00123 {
return m_TimeStep; }
00124
00127
virtual void *GetGlobalDataPointer()
const
00128 {
00129
GlobalDataStruct *global =
new GlobalDataStruct();
00130
return global;
00131 }
00132
00134
virtual void ReleaseGlobalDataPointer(
void *GlobalData )
const
00135
{
delete (GlobalDataStruct *) GlobalData; }
00136
00138
virtual void InitializeIteration();
00139
00142
virtual PixelType ComputeUpdate(
const NeighborhoodType &neighborhood,
00143
void *globalData,
00144
const FloatOffsetType &offset =
FloatOffsetType(0.0));
00145
00146
protected:
00147 NCCRegistrationFunction();
00148 ~NCCRegistrationFunction() {}
00149
void PrintSelf(std::ostream& os,
Indent indent)
const;
00150
00152
typedef ConstNeighborhoodIterator<FixedImageType> FixedImageNeighborhoodIteratorType;
00153
00156
struct GlobalDataStruct
00157 {
00158 FixedImageNeighborhoodIteratorType m_FixedImageIterator;
00159 };
00160
00161
00162
private:
00163 NCCRegistrationFunction(
const Self&);
00164
void operator=(
const Self&);
00165
00167
const double * m_FixedImageSpacing;
00168
const double * m_FixedImageOrigin;
00169
00171
GradientCalculatorPointer m_FixedImageGradientCalculator;
00172
00173
00175
InterpolatorPointer m_MovingImageInterpolator;
00176
00178
TimeStepType m_TimeStep;
00179
00181
double m_DenominatorThreshold;
00182
00184
double m_IntensityDifferenceThreshold;
00185
00186
mutable double m_MetricTotal;
00187
00188 };
00189
00190
00191 }
00192
00193
#ifndef ITK_MANUAL_INSTANTIATION
00194
#include "itkNCCRegistrationFunction.txx"
00195
#endif
00196
00197
#endif