#include <itkVariationalRegistrationNCCFunction.h>

Inheritance diagram for itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >:

Collaboration diagram for itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >:

Detailed Description

template<class TFixedImage, class TMovingImage, class TDisplacementField>
class itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >

This class computes NCC forces in the variational registration framework.

This class implements NCC forces as given in Hermosillo, Chefd'Hotel, and Faugeras. "Variational methods for multimodal image matching." IJCV 50(3), 2002: 329-343 and Avants et al. "Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain." Medical image analysis 12(1), 2008: 26-41 (except Jacobian term). We define the derivative of NCC between two images as:

$f^{NCC}(x)=\tau\kappa\frac{2\sum_w (F-\bar{F})(M-\bar{M})}{\sum_w (F-\bar{F})^2\ \sum_w (M-\bar{M})^2}\left((M-\bar{M}) - \frac{\sum_w (F-\bar{F})(M-\bar{M})}{\sum_w (F-\bar{F})^2}(F-\bar{F})\right) \nabla M(x+u(x))$

$\tau$ is the step size and $\kappa$ is the mean squared spacing. Use SetRadius() (see FiniteDifferenceFunction) to set the size of the neighbourhood to compute local mean values $\bar{F}$ and $\bar{M}$ and the local sums. Alternative, the classical gradient $\nabla M(x+u(x))$ can be replaced by $\nabla F(x)$ or $\frac{\nabla F(x) + \nabla M(x+u(x))}{2}$ .

See Also: VariationalRegistrationFilter; VariationalRegistrationFunction

Note: This class was developed with funding from the German Research Foundation (DFG: EH 224/3-1 and HA 235/9-1).

Author: Alexander Schmidt-Richberg; Rene Werner; Jan Ehrhardt

Definition at line 64 of file itkVariationalRegistrationNCCFunction.h.

Public Types
typedef SmartPointer< const Self >	ConstPointer

typedef Superclass::DisplacementFieldType	DisplacementFieldType

typedef Superclass::DisplacementFieldTypePointer	DisplacementFieldTypePointer

typedef Superclass::FixedImagePointer	FixedImagePointer

typedef Superclass::FixedImageType	FixedImageType

typedef Superclass::FloatOffsetType	FloatOffsetType

typedef GradientCalculatorType::Pointer	GradientCalculatorPointer

typedef CentralDifferenceImageFunction < FixedImageType >	GradientCalculatorType

typedef FixedImageType::IndexType	IndexType

typedef Superclass::MaskImagePointer	MaskImagePointer

typedef Superclass::MaskImageType	MaskImageType

typedef Superclass::MovingImagePointer	MovingImagePointer

typedef Superclass::MovingImageType	MovingImageType

typedef Superclass::NeighborhoodType	NeighborhoodType

typedef Superclass::PixelType	PixelType

typedef SmartPointer< Self >	Pointer

typedef Superclass::RadiusType	RadiusType

typedef VariationalRegistrationNCCFunction	Self

typedef FixedImageType::SizeType	SizeType

typedef FixedImageType::SpacingType	SpacingType

typedef VariationalRegistrationFunction < TFixedImage, TMovingImage, TDisplacementField >	Superclass

Public Types inherited from itk::VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >
typedef SmartPointer< const Self >	ConstPointer

typedef TDisplacementField	DisplacementFieldType

typedef DisplacementFieldType::ConstPointer	DisplacementFieldTypePointer

typedef FixedImageType::ConstPointer	FixedImagePointer

typedef TFixedImage	FixedImageType

typedef unsigned char	MaskImagePixelType

typedef MaskImageType::ConstPointer	MaskImagePointer

typedef Image < MaskImagePixelType, ImageDimension >	MaskImageType

typedef MovingImageType::ConstPointer	MovingImagePointer

typedef TMovingImage	MovingImageType

typedef MovingImageWarperType::Pointer	MovingImageWarperPointer

typedef itk::ContinuousBorderWarpImageFilter < FixedImageType, WarpedImageType, DisplacementFieldType >	MovingImageWarperType

typedef SmartPointer< Self >	Pointer

typedef VariationalRegistrationFunction	Self

typedef FiniteDifferenceFunction < TDisplacementField >	Superclass

typedef Superclass::TimeStepType	TimeStepType

typedef WarpedImageType::ConstPointer	WarpedImagePointer

typedef TFixedImage	WarpedImageType

Public Types inherited from itk::FiniteDifferenceFunction< TDisplacementField >
typedef SmartPointer< const Self >	ConstPointer

typedef ZeroFluxNeumannBoundaryCondition < ImageType >	DefaultBoundaryConditionType

typedef Vector< float, itkGetStaticConstMacro(ImageDimension) >	FloatOffsetType

typedef TDisplacementField	ImageType

typedef Vector< PixelRealType, itkGetStaticConstMacro(ImageDimension) >	NeighborhoodScalesType

typedef ConstNeighborhoodIterator < TDisplacementField, DefaultBoundaryConditionType >	NeighborhoodType

typedef double	PixelRealType

typedef ImageType::PixelType	PixelType

typedef SmartPointer< Self >	Pointer

typedef ConstNeighborhoodIterator < TDisplacementField > ::RadiusType	RadiusType

typedef FiniteDifferenceFunction	Self

typedef LightObject	Superclass

typedef double	TimeStepType

Public Types inherited from itk::LightObject
typedef SmartPointer< const Self >	ConstPointer

typedef SmartPointer< Self >	Pointer

typedef LightObject	Self

Public Member Functions
virtual Superclass::TimeStepType	ComputeGlobalTimeStep (void *gd) const

virtual PixelType	ComputeUpdate (const NeighborhoodType &neighborhood, void *globalData, const FloatOffsetType &offset=FloatOffsetType(0.0))

virtual ::itk::LightObject::Pointer	CreateAnother () const

virtual const char *	GetNameOfClass () const

virtual void	InitializeIteration ()

virtual void	SetGradientTypeToFixedImage ()

virtual void	SetGradientTypeToSymmetric ()

virtual void	SetGradientTypeToWarpedMovingImage ()

Public Member Functions inherited from itk::VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >
virtual const DisplacementFieldType *	GetDisplacementField (void) const

virtual const FixedImageType *	GetFixedImage (void) const

virtual void *	GetGlobalDataPointer () const

virtual MaskImagePixelType	GetMaskBackgroundThreshold (void) const

virtual const MaskImageType *	GetMaskImage (void) const

virtual double	GetMetric () const

virtual const MovingImageType *	GetMovingImage (void) const

virtual const MovingImageWarperType *	GetMovingImageWarper (void) const

virtual double	GetRMSChange () const

virtual const TimeStepType	GetTimeStep (void) const

virtual void	ReleaseGlobalDataPointer (void *GlobalData) const

virtual void	SetDisplacementField (DisplacementFieldType *ptr)

virtual void	SetFixedImage (const FixedImageType *ptr)

virtual void	SetMaskBackgroundThreshold (MaskImagePixelType threshold)

virtual void	SetMaskImage (const MaskImageType *ptr)

virtual void	SetMovingImage (const MovingImageType *ptr)

virtual void	SetMovingImageWarper (MovingImageWarperType *ptr)

virtual void	SetTimeStep (TimeStepType timeStep)

Public Member Functions inherited from itk::FiniteDifferenceFunction< TDisplacementField >
const NeighborhoodScalesType	ComputeNeighborhoodScales () const

const RadiusType &	GetRadius () const

void	GetScaleCoefficients (PixelRealType vals[ImageDimension]) const

void	SetRadius (const RadiusType &r)

void	SetScaleCoefficients (PixelRealType vals[ImageDimension])

Public Member Functions inherited from itk::LightObject
virtual void	Delete ()

virtual int	GetReferenceCount () const

	itkCloneMacro (Self)

void	Print (std::ostream &os, Indent indent=0) const

virtual void	Register () const

virtual void	SetReferenceCount (int)

virtual void	UnRegister () const noexcept

Static Public Member Functions
static Pointer	New ()

Static Public Member Functions inherited from itk::LightObject
static void	BreakOnError ()

static Pointer	New ()

Static Public Attributes
static const unsigned int	ImageDimension = Superclass::ImageDimension

Static Public Attributes inherited from itk::VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >
static const unsigned int	ImageDimension = Superclass::ImageDimension

Static Public Attributes inherited from itk::FiniteDifferenceFunction< TDisplacementField >
static const unsigned int	ImageDimension

Protected Types
typedef ConstNeighborhoodIterator < FixedImageType >	FixedImageNeighborhoodIteratorType

typedef Superclass::GlobalDataStruct	GlobalDataStruct

enum	GradientType { GRADIENT_TYPE_WARPED = 0, GRADIENT_TYPE_FIXED = 1, GRADIENT_TYPE_SYMMETRIC = 2 }

Protected Member Functions
virtual void	PrintSelf (std::ostream &os, Indent indent) const

	VariationalRegistrationNCCFunction ()

	~VariationalRegistrationNCCFunction ()

Protected Member Functions inherited from itk::VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >
virtual const WarpedImagePointer	GetWarpedImage (void) const

	VariationalRegistrationFunction ()

virtual void	WarpMovingImage (void)

	~VariationalRegistrationFunction ()

Protected Member Functions inherited from itk::FiniteDifferenceFunction< TDisplacementField >
	FiniteDifferenceFunction ()

	~FiniteDifferenceFunction ()

Protected Member Functions inherited from itk::LightObject
virtual LightObject::Pointer	InternalClone () const

	LightObject ()

virtual void	PrintHeader (std::ostream &os, Indent indent) const

virtual void	PrintTrailer (std::ostream &os, Indent indent) const

virtual	~LightObject ()

Private Member Functions
void	operator= (const Self &)

	VariationalRegistrationNCCFunction (const Self &)

Private Attributes
GradientCalculatorPointer	m_FixedImageGradientCalculator

GradientType	m_GradientType

double	m_Normalizer

GradientCalculatorPointer	m_WarpedImageGradientCalculator

Additional Inherited Members
Protected Attributes inherited from itk::FiniteDifferenceFunction< TDisplacementField >
RadiusType	m_Radius

PixelRealType	m_ScaleCoefficients [ImageDimension]

Protected Attributes inherited from itk::LightObject
AtomicInt< int >	m_ReferenceCount

Member Typedef Documentation

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef SmartPointer< const Self > itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::ConstPointer

Definition at line 73 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::DisplacementFieldType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::DisplacementFieldType

Deformation field type.

Definition at line 101 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::DisplacementFieldTypePointer itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::DisplacementFieldTypePointer

Definition at line 103 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef ConstNeighborhoodIterator<FixedImageType> itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::FixedImageNeighborhoodIteratorType

protected

FixedImage image neighborhood iterator type.

Definition at line 151 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::FixedImagePointer itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::FixedImagePointer

Definition at line 90 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::FixedImageType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::FixedImageType

FixedImage image type.

Definition at line 89 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::FloatOffsetType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::FloatOffsetType

Definition at line 109 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::GlobalDataStruct itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::GlobalDataStruct

protected

Definition at line 153 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef GradientCalculatorType::Pointer itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::GradientCalculatorPointer

Definition at line 114 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef CentralDifferenceImageFunction<FixedImageType> itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::GradientCalculatorType

Gradient calculator type.

Definition at line 113 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef FixedImageType::IndexType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::IndexType

Definition at line 96 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::MaskImagePointer itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::MaskImagePointer

Definition at line 94 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::MaskImageType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::MaskImageType

MaskImage image type.

Definition at line 93 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::MovingImagePointer itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::MovingImagePointer

Definition at line 86 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::MovingImageType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::MovingImageType

MovingImage image type.

Definition at line 85 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::NeighborhoodType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::NeighborhoodType

Definition at line 108 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::PixelType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::PixelType

Inherit some types from the superclass.

Definition at line 106 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef SmartPointer< Self > itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::Pointer

Definition at line 72 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef Superclass::RadiusType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::RadiusType

Definition at line 107 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef VariationalRegistrationNCCFunction itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::Self

Standard class typedefs.

Definition at line 69 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef FixedImageType::SizeType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::SizeType

Definition at line 97 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef FixedImageType::SpacingType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::SpacingType

Definition at line 98 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

typedef VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField > itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::Superclass

Definition at line 71 of file itkVariationalRegistrationNCCFunction.h.

Member Enumeration Documentation

template<class TFixedImage , class TMovingImage , class TDisplacementField >

enum itk::VariationalRegistrationNCCFunction::GradientType

protected

Type of available image forces

Enumerator
GRADIENT_TYPE_WARPED
GRADIENT_TYPE_FIXED
GRADIENT_TYPE_SYMMETRIC

Definition at line 156 of file itkVariationalRegistrationNCCFunction.h.

Constructor & Destructor Documentation

template<class TFixedImage , class TMovingImage , class TDisplacementField >

itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::VariationalRegistrationNCCFunction ( )

protected

template<class TFixedImage , class TMovingImage , class TDisplacementField >

itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::~VariationalRegistrationNCCFunction ( )

inlineprotected

Definition at line 145 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::VariationalRegistrationNCCFunction ( const Self & )

private

Member Function Documentation

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual Superclass::TimeStepType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::ComputeGlobalTimeStep ( void * gd ) const

virtual

Computes the time step for an update. Returns the constant time step scaled with the mean squared spacing.

See Also: SetTimeStep()

Reimplemented from itk::VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual PixelType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::ComputeUpdate	(	const NeighborhoodType &	neighborhood,
		void *	globalData,
		const FloatOffsetType &	offset = `FloatOffsetType(0.0)`
	)

virtual

This method is called by a finite difference solver image filter at each pixel that does not lie on a data set boundary

Implements itk::FiniteDifferenceFunction< TDisplacementField >.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual::itk::LightObject::Pointer itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::CreateAnother ( ) const

virtual

Create an object from an instance, potentially deferring to a factory. This method allows you to create an instance of an object that is exactly the same type as the referring object. This is useful in cases where an object has been cast back to a base class.

Reimplemented from itk::LightObject.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual const char* itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::GetNameOfClass ( ) const

virtual

Run-time type information (and related methods).

Reimplemented from itk::VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual void itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::InitializeIteration ( )

virtual

Set the object's state before each iteration.

Reimplemented from itk::VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

static Pointer itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::New ( )

static

Method for creation through the object factory.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

void itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::operator= ( const Self & )

private

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual void itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::PrintSelf	(	std::ostream &	os,
		Indent	indent
	)		const

protectedvirtual

Print information about the filter.

Reimplemented from itk::VariationalRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual void itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::SetGradientTypeToFixedImage ( )

inlinevirtual

Select that the fixed image gradient is used for computing the forces.

Definition at line 126 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual void itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::SetGradientTypeToSymmetric ( )

inlinevirtual

Select that fixed and warped image gradients are used for computing the forces.

Definition at line 135 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

virtual void itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::SetGradientTypeToWarpedMovingImage ( )

inlinevirtual

Select that the warped image gradient is used for computing the forces.

Definition at line 130 of file itkVariationalRegistrationNCCFunction.h.

Member Data Documentation

template<class TFixedImage , class TMovingImage , class TDisplacementField >

const unsigned int itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::ImageDimension = Superclass::ImageDimension

static

Get image dimension.

Definition at line 82 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

GradientCalculatorPointer itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::m_FixedImageGradientCalculator

private

Function to compute derivatives of the fixed image.

Definition at line 167 of file itkVariationalRegistrationNCCFunction.h.

template<class TFixedImage , class TMovingImage , class TDisplacementField >

GradientType itk::VariationalRegistrationNCCFunction< TFixedImage, TMovingImage, TDisplacementField >::m_GradientType

private

Set if warped or fixed image gradient is used for force computation.