#include <itkTimeVaryingBSplineVelocityFieldImageRegistrationMethod.h>

Inheritance diagram for itk::TimeVaryingBSplineVelocityFieldImageRegistrationMethod< TFixedImage, TMovingImage, TOutputTransform >:

Collaboration diagram for itk::TimeVaryingBSplineVelocityFieldImageRegistrationMethod< TFixedImage, TMovingImage, TOutputTransform >:

List of all members.

Public Types
typedef Superclass::CompositeTransformType	CompositeTransformType
typedef SmartPointer< const Self >	ConstPointer
typedef DecoratedOutputTransformType::Pointer	DecoratedOutputTransformPointer
typedef Superclass::DecoratedOutputTransformType	DecoratedOutputTransformType
typedef OutputTransformType::DerivativeType	DerivativeType
typedef DerivativeType::ValueType	DerivativeValueType
typedef OutputTransformType::DisplacementFieldType	DisplacementFieldType
typedef TimeVaryingVelocityFieldControlPointLatticeType::PixelType	DisplacementVectorType
typedef FixedImageType::Pointer	FixedImagePointer
typedef TFixedImage	FixedImageType
typedef MetricType::MeasureType	MeasureType
typedef MetricType::Pointer	MetricPointer
typedef Superclass::MetricType	MetricType
typedef MovingImageType::Pointer	MovingImagePointer
typedef TMovingImage	MovingImageType
typedef Array< SizeValueType >	NumberOfIterationsArrayType
typedef OutputTransformType::Pointer	OutputTransformPointer
typedef TOutputTransform	OutputTransformType
typedef SmartPointer< Self >	Pointer
typedef OutputTransformType::ScalarType	RealType
typedef TimeVaryingBSplineVelocityFieldImageRegistrationMethod	Self
typedef ImageRegistrationMethodv4 < TFixedImage, TMovingImage, TOutputTransform >	Superclass
typedef OutputTransformType::TimeVaryingVelocityFieldControlPointLatticePointer	TimeVaryingVelocityFieldControlPointLatticePointer
typedef OutputTransformType::TimeVaryingVelocityFieldControlPointLatticeType	TimeVaryingVelocityFieldControlPointLatticeType
typedef OutputTransformType::TimeVaryingVelocityFieldControlPointLatticePointer	TimeVaryingVelocityFieldPointer
typedef OutputTransformType::TimeVaryingVelocityFieldControlPointLatticeType	TimeVaryingVelocityFieldType
typedef MetricType::VirtualImageType	VirtualImageType
Public Member Functions
virtual ::itk::LightObject::Pointer	CreateAnother (void) const
virtual const char *	GetNameOfClass () const

virtual void	SetLearningRate (RealType _arg)
virtual RealType	GetLearningRate () const

virtual void	SetNumberOfIterationsPerLevel (NumberOfIterationsArrayType _arg)
virtual NumberOfIterationsArrayType	GetNumberOfIterationsPerLevel () const

virtual void	SetConvergenceThreshold (RealType _arg)
virtual RealType	GetConvergenceThreshold () const

virtual void	SetConvergenceWindowSize (unsigned int _arg)
virtual unsigned int	GetConvergenceWindowSize () const

virtual void	SetNumberOfTimePointSamples (SizeValueType _arg)
virtual SizeValueType	GetNumberOfTimePointSamples () const
Static Public Member Functions
static Pointer	New ()
Static Public Attributes
static const unsigned int	ImageDimension = TFixedImage::ImageDimension
Protected Member Functions
virtual void	GenerateData ()
virtual void	PrintSelf (std::ostream &os, Indent indent) const
virtual void	StartOptimization ()
	TimeVaryingBSplineVelocityFieldImageRegistrationMethod ()
virtual	~TimeVaryingBSplineVelocityFieldImageRegistrationMethod ()
Private Member Functions
void	operator= (const Self &)
	TimeVaryingBSplineVelocityFieldImageRegistrationMethod (const Self &)
Private Attributes
RealType	m_ConvergenceThreshold
unsigned int	m_ConvergenceWindowSize
RealType	m_LearningRate
NumberOfIterationsArrayType	m_NumberOfIterationsPerLevel
SizeValueType	m_NumberOfTimePointSamples

Detailed Description

template<typename TFixedImage, typename TMovingImage, typename TOutputTransform = TimeVaryingBSplineVelocityFieldTransform<double, GetImageDimension<TFixedImage>::ImageDimension>>
class itk::TimeVaryingBSplineVelocityFieldImageRegistrationMethod< TFixedImage, TMovingImage, TOutputTransform >

Interface method for the current registration framework using the time varying velocity field transform.

Output: The output is the updated transform which has been added to the composite transform.

This derived class from the SimpleImageRegistrationMethod class is specialized to handle time-varying velocity field transforms of which there are currently 2:

TimeVaryingDisplacementFieldTransform
GaussianSmoothingOnUpdateTimeVaryingDisplacementFieldTransform

The latter is derived from the former and performs an optional spatial and temporal smoothing on the update and total velocity fields. Integration of the velocity field is performed using 4th order Runge Kutta and is performed using the class itkTimeVaryingBSplineVelocityFieldIntegrationImageFilter.

Optimization occurs in an iterative fashion where for each sample time point, t, in the velocity field, we integrate the velocity field in the range [0, t] to yield the displacement field which warps fixed image to time point t. Simultaneously, we integrate the velocity field in the range [t, 1] to yield the displacement field transform which warps the moving image to time point t. The metric derivative for each time point of the velocity field calculated in this way produces the normalized update field (or gradient) of the velocity field to be added to the total field at each iteration after being multiplied by the learning rate and optionally smoothed. Mathematically, this can be described as

$v_{total} = G_1( v_{total} + \lambda * G_2( v_{update} ) )$

where

$ G_1 = $ gaussian smoothing on the total field $ G_2 = $ gaussian smoothing on the update field $\lambda =$ learning rate $v_{update} =$ the normalized velocity field where we normalize the velocity field at each time point separately by the max norm of the field at that time point. This is done due to a weakly necessary (but not sufficient) condition being that the velocity field have a constant norm for all time points.