ITK
4.13.0
Insight Segmentation and Registration Toolkit
|
#include <itkLevelSetMotionRegistrationFilter.h>
Deformably register two images using level set motion.
LevelSetMotionFilter implements a deformable registration algorithm that aligns a fixed and a moving image under level set motion. The equations of motion are similar to those of the DemonsRegistrationFilter. The main differences are: (1) Gradients of the moving image are calculated on a smoothed image while intensity difference are measured on the original images (2) Magnitude of the motion vector is a function of the differences in intensity between the fixed and moving pixel. An adaptive timestep is calculated based on the maximum motion vector over the entire field to ensure stability. The timestep also implictly converts the motion vector measured in units of intensity to a vector measured in physical units. Demons, on the other hand, defines its motion vectors as function of both the intensity differences and gradient magnitude at each respective pixel. Consider two separate pixels with the same intensity differences between the corresponding fixed and moving pixel pairs. In demons, the motion vector of the pixel over a low gradient region will be larger than the motion vector of the pixel over a large gradient region. This leads to an unstable vector field. In the levelset approach, the motion vectors will be proportional to the gradients, scaled by the maximum gradient over the entire field. The pixel with at the lower gradient position will more less than the pixel at the higher gradient position. (3) Gradients are calculated using minmod finite difference instead of using central differences.
A deformation field is represented as a image whose pixel type is some vector type with at least N elements, where N is the dimension of the fixed image. The vector type must support element access via operator []. It is assumed that the vector elements behave like floating point scalars.
This class is templated over the fixed image type, moving image type and the deformation field type.
The input fixed and moving images are set via methods SetFixedImage and SetMovingImage respectively. An initial deformation field maybe set via SetInitialDisplacementField or SetInput. If no initial field is set, a zero field is used as the initial condition.
The algorithm has one parameters: the number of iteration to be performed.
The output deformation field can be obtained via methods GetOutput or GetDisplacementField.
This class make use of the finite difference solver hierarchy. Update for each iteration is computed in LevelSetMotionFunction.
Ref: B.C. Vemuri, J. Ye, Y. Chen, C.M. Leonard. "Image registration via level-set motion: applications to atlas-based segmentation". Medical Image Analysis. Vol. 7. pp. 1-20. 2003.
Definition at line 91 of file itkLevelSetMotionRegistrationFilter.h.
Static Public Member Functions | |
static Pointer | New () |
Static Public Member Functions inherited from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField > | |
static Pointer | New () |
Static Public Member Functions inherited from itk::Object | |
static bool | GetGlobalWarningDisplay () |
static void | GlobalWarningDisplayOff () |
static void | GlobalWarningDisplayOn () |
static Pointer | New () |
static void | SetGlobalWarningDisplay (bool flag) |
Static Public Member Functions inherited from itk::LightObject | |
static void | BreakOnError () |
static Pointer | New () |
typedef SmartPointer< const Self > itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::ConstPointer |
Definition at line 100 of file itkLevelSetMotionRegistrationFilter.h.
typedef Superclass::DisplacementFieldPointer itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::DisplacementFieldPointer |
Definition at line 124 of file itkLevelSetMotionRegistrationFilter.h.
typedef Superclass::DisplacementFieldType itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::DisplacementFieldType |
Deformation field type.
Definition at line 122 of file itkLevelSetMotionRegistrationFilter.h.
typedef Superclass::FiniteDifferenceFunctionType itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::FiniteDifferenceFunctionType |
FiniteDifferenceFunction type.
Definition at line 128 of file itkLevelSetMotionRegistrationFilter.h.
typedef Superclass::FixedImagePointer itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::FixedImagePointer |
Definition at line 114 of file itkLevelSetMotionRegistrationFilter.h.
typedef Superclass::FixedImageType itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::FixedImageType |
FixedImage image type.
Definition at line 113 of file itkLevelSetMotionRegistrationFilter.h.
typedef LevelSetMotionRegistrationFunction< FixedImageType, MovingImageType, DisplacementFieldType > itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::LevelSetMotionFunctionType |
LevelSetMotionFilterFunction type.
Definition at line 132 of file itkLevelSetMotionRegistrationFilter.h.
typedef Superclass::MovingImagePointer itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::MovingImagePointer |
Definition at line 118 of file itkLevelSetMotionRegistrationFilter.h.
typedef Superclass::MovingImageType itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::MovingImageType |
MovingImage image type.
Definition at line 117 of file itkLevelSetMotionRegistrationFilter.h.
typedef SmartPointer< Self > itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::Pointer |
Definition at line 99 of file itkLevelSetMotionRegistrationFilter.h.
typedef LevelSetMotionRegistrationFilter itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::Self |
Standard class typedefs.
Definition at line 97 of file itkLevelSetMotionRegistrationFilter.h.
typedef PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField > itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::Superclass |
Definition at line 98 of file itkLevelSetMotionRegistrationFilter.h.
typedef Superclass::TimeStepType itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::TimeStepType |
Inherit types from superclass.
Definition at line 107 of file itkLevelSetMotionRegistrationFilter.h.
|
protected |
|
inlineoverrideprotected |
Definition at line 183 of file itkLevelSetMotionRegistrationFilter.h.
|
overrideprotectedvirtual |
Apply update.
Reimplemented from itk::DenseFiniteDifferenceImageFilter< TDisplacementField, TDisplacementField >.
|
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::Object.
|
virtual |
|
virtual |
|
virtual |
|
virtual |
|
virtual |
Get the metric value. The metric value is the mean square difference in intensity between the fixed image and transforming moving image computed over the the overlapping region between the two images. This is value is only available for the previous iteration and NOT the current iteration.
|
virtual |
Run-time type information (and related methods).
Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >.
|
overrideprotectedvirtual |
This method returns true when the current iterative solution of the equation has met the criteria to stop solving. This version calls the superclass' version but also Halts if the RMSChange is zero.
Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >.
|
overrideprotectedvirtual |
Initialize the state of filter and equation before each iteration.
Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >.
|
static |
Method for creation through the object factory.
|
overrideprotectedvirtual |
Methods invoked by Print() to print information about the object including superclasses. Typically not called by the user (use Print() instead) but used in the hierarchical print process to combine the output of several classes.
Reimplemented from itk::InPlaceImageFilter< TDisplacementField, TDisplacementField >.
|
virtual |
Set/Get the parameter alpha. Alpha is added to the calculated gradient magnitude prior to normalizing the gradient to protect against numerical instability as the gradient magnitude approaches zero. This should be set as a small fraction of the intensity dynamic range, for instance 0.04%. Default is the absolute (not percentage) value of 0.1.
|
virtual |
Set/Get the threshold below which the gradient magnitude is considered the zero vector. Default is 1e-9.
|
virtual |
Set/Get the standard deviation used for smoothing the moving image prior to calculating gradients. The standard deviation is measured in physical units (for instance mm). Note that this smoothing value is not to be confused with the PDEDeformableRegistrationFilter::SetStandardDeviations() method. The method in PDEDeformableRegistrationFilter is for setting the smoothing parameters for regularizing the deformation field between interations. Those smoothing parameters are set in pixel units not physical units. Deformation field smoothing is not done by default in LevelSetMotionRegistration. This smoothing parameter is to condition the gradient calculation and parameter is specified in physical units.
|
virtual |
Set/Get the threshold below which the absolute difference of intensity yields a match. When the intensities match between a moving and fixed image pixel, the update vector (for that iteration) will be the zero vector. Default is 0.001.