ITK
4.4.0
Insight Segmentation and Registration Toolkit
|
#include <itkPDEDeformableRegistrationFilter.h>
Deformably register two images using a PDE algorithm.
PDEDeformableRegistrationFilter is a base case for filter implementing a PDE deformable algorithm that register two images by computing the displacement field which will map a moving image onto a fixed image.
A displacement 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 displacement Field type.
The input fixed and moving images are set via methods SetFixedImage and SetMovingImage respectively. An initial displacement field maybe set via SetInitialDisplacementField or SetInput. If no initial field is set, a zero field is used as the initial condition.
The output displacement field can be obtained via methods GetOutput or GetDisplacementField.
The PDE algorithm is run for a user defined number of iterations. Typically the PDE algorithm requires period Gaussin smoothing of the displacement field to enforce an elastic-like condition. The amount of smoothing is governed by a set of user defined standard deviations (one for each dimension).
In terms of memory, this filter keeps two internal buffers: one for storing the intermediate updates to the field and one for double-buffering when smoothing the displacement field. Both buffers are the same type and size as the output displacement field.
This class make use of the finite difference solver hierarchy. Update for each iteration is computed using a PDEDeformableRegistrationFunction.
Definition at line 73 of file itkPDEDeformableRegistrationFilter.h.
Static Public Member Functions | |
static Pointer | New () |
Static Public Attributes | |
static const unsigned int | ImageDimension = Superclass::ImageDimension |
Static Public Attributes inherited from itk::DenseFiniteDifferenceImageFilter< TDisplacementField, TDisplacementField > | |
static const unsigned int | ImageDimension |
Static Public Attributes inherited from itk::FiniteDifferenceImageFilter< TDisplacementField, TDisplacementField > | |
static const unsigned int | ImageDimension |
Static Public Attributes inherited from itk::InPlaceImageFilter< TDisplacementField, TDisplacementField > | |
static const unsigned int | InputImageDimension |
static const unsigned int | OutputImageDimension |
Static Public Attributes inherited from itk::ImageToImageFilter< TDisplacementField, TDisplacementField > | |
static const unsigned int | InputImageDimension |
static const unsigned int | OutputImageDimension |
Static Public Attributes inherited from itk::ImageSource< TDisplacementField > | |
static const unsigned int | OutputImageDimension |
Private Member Functions | |
void | operator= (const Self &) |
PDEDeformableRegistrationFilter (const Self &) | |
Private Attributes | |
double | m_MaximumError |
unsigned int | m_MaximumKernelWidth |
bool | m_SmoothDisplacementField |
bool | m_SmoothUpdateField |
StandardDeviationsType | m_StandardDeviations |
bool | m_StopRegistrationFlag |
DisplacementFieldPointer | m_TempField |
StandardDeviationsType | m_UpdateFieldStandardDeviations |
Additional Inherited Members | |
Protected Types inherited from itk::DenseFiniteDifferenceImageFilter< TDisplacementField, TDisplacementField > | |
typedef UpdateBufferType::RegionType | ThreadRegionType |
Static Protected Member Functions inherited from itk::ImageSource< TDisplacementField > | |
static const ImageRegionSplitterBase * | GetGlobalDefaultSplitter () |
static ITK_THREAD_RETURN_TYPE | ThreaderCallback (void *arg) |
Protected Attributes inherited from itk::FiniteDifferenceImageFilter< TDisplacementField, TDisplacementField > | |
IdentifierType | m_ElapsedIterations |
bool | m_IsInitialized |
bool | m_ManualReinitialization |
double | m_MaximumRMSError |
IdentifierType | m_NumberOfIterations |
double | m_RMSChange |
typedef SmartPointer< const Self > itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::ConstPointer |
Definition at line 82 of file itkPDEDeformableRegistrationFilter.h.
typedef DisplacementFieldType::Pointer itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::DisplacementFieldPointer |
Definition at line 103 of file itkPDEDeformableRegistrationFilter.h.
typedef TDisplacementField itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::DisplacementFieldType |
Deformation field type.
Definition at line 102 of file itkPDEDeformableRegistrationFilter.h.
typedef Superclass::FiniteDifferenceFunctionType itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::FiniteDifferenceFunctionType |
FiniteDifferenceFunction type.
Definition at line 115 of file itkPDEDeformableRegistrationFilter.h.
typedef FixedImageType::ConstPointer itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::FixedImageConstPointer |
Definition at line 94 of file itkPDEDeformableRegistrationFilter.h.
typedef FixedImageType::Pointer itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::FixedImagePointer |
Definition at line 93 of file itkPDEDeformableRegistrationFilter.h.
typedef TFixedImage itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::FixedImageType |
FixedImage image type.
Definition at line 89 of file itkPDEDeformableRegistrationFilter.h.
typedef MovingImageType::ConstPointer itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::MovingImageConstPointer |
Definition at line 99 of file itkPDEDeformableRegistrationFilter.h.
typedef MovingImageType::Pointer itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::MovingImagePointer |
Definition at line 98 of file itkPDEDeformableRegistrationFilter.h.
typedef TMovingImage itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::MovingImageType |
MovingImage image type.
Definition at line 97 of file itkPDEDeformableRegistrationFilter.h.
typedef Superclass::OutputImageType itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::OutputImageType |
Types inherithed from the superclass
Definition at line 111 of file itkPDEDeformableRegistrationFilter.h.
typedef PDEDeformableRegistrationFunction< FixedImageType, MovingImageType, DisplacementFieldType > itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::PDEDeformableRegistrationFunctionType |
PDEDeformableRegistrationFilterFunction type.
Definition at line 119 of file itkPDEDeformableRegistrationFilter.h.
typedef SmartPointer< Self > itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::Pointer |
Definition at line 81 of file itkPDEDeformableRegistrationFilter.h.
typedef PDEDeformableRegistrationFilter itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::Self |
Standard class typedefs.
Definition at line 79 of file itkPDEDeformableRegistrationFilter.h.
typedef FixedArray< double, ImageDimension > itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::StandardDeviationsType |
Definition at line 172 of file itkPDEDeformableRegistrationFilter.h.
typedef DenseFiniteDifferenceImageFilter< TDisplacementField, TDisplacementField > itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::Superclass |
Definition at line 80 of file itkPDEDeformableRegistrationFilter.h.
|
protected |
|
inlineprotected |
Definition at line 245 of file itkPDEDeformableRegistrationFilter.h.
|
private |
|
protectedvirtual |
A simple method to copy the data from the input to the output. If the input does not exist, a zero field is written to the output.
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.
Reimplemented in itk::DiffeomorphicDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, itk::FastSymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, and itk::SymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >.
|
protectedvirtual |
It is difficult to compute in advance the input moving image region required to compute the requested output region. Thus the safest thing to do is to request for the whole moving image.
For the fixed image and displacement field, the input requested region set to be the same as that of the output requested region.
Reimplemented from itk::FiniteDifferenceImageFilter< TDisplacementField, TDisplacementField >.
|
protectedvirtual |
By default the output displacement field has the same Spacing, Origin and LargestPossibleRegion as the input/initial displacement field. If the initial displacement field is not set, the output information is copied from the fixed image.
Reimplemented from itk::ProcessObject.
|
inline |
Get output displacement field.
Definition at line 142 of file itkPDEDeformableRegistrationFilter.h.
const FixedImageType* itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::GetFixedImage | ( | void | ) | const |
Get the fixed image.
|
virtual |
Set/Get the desired maximum error of the Guassian kernel approximate.
|
virtual |
Set/Get the desired limits of the Gaussian kernel width.
const MovingImageType* itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::GetMovingImage | ( | void | ) | const |
Get the moving image.
|
virtual |
Run-time type information (and related methods)
Reimplemented from itk::DenseFiniteDifferenceImageFilter< TDisplacementField, TDisplacementField >.
Reimplemented in itk::CurvatureRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField, TImageForceFunction >, itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, itk::DiffeomorphicDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, itk::FastSymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, itk::SymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, and itk::DemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >.
|
virtual |
Get the number of valid inputs. For PDEDeformableRegistration, this checks whether the fixed and moving images have been set. While PDEDeformableRegistration can take a third input as an initial displacement field, this input is not a required input.
Reimplemented from itk::ProcessObject.
|
virtual |
Set/Get whether the displacement field is smoothed (regularized). Smoothing the displacement yields a solution elastic in nature. If SmoothDisplacementField is on, then the displacement field is smoothed with a Gaussian whose standard deviations are specified with SetStandardDeviations()
|
virtual |
Set/Get whether the update field is smoothed (regularized). Smoothing the update field yields a solution viscous in nature. If SmoothUpdateField is on, then the update field is smoothed with a Gaussian whose standard deviations are specified with SetUpdateFieldStandardDeviations()
|
virtual |
Get the Gaussian smoothing standard deviations use for smoothing the displacement field.
|
virtual |
Get the Gaussian smoothing standard deviations used for smoothing the update field.
|
inlineprotectedvirtual |
Supplies the halting criteria for this class of filters. The algorithm will stop after a user-specified number of iterations.
Reimplemented from itk::FiniteDifferenceImageFilter< TDisplacementField, TDisplacementField >.
Reimplemented in itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >.
Definition at line 250 of file itkPDEDeformableRegistrationFilter.h.
|
protectedvirtual |
This method is called before iterating the solution.
Reimplemented from itk::FiniteDifferenceImageFilter< TDisplacementField, TDisplacementField >.
Reimplemented in itk::CurvatureRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField, TImageForceFunction >.
|
protectedvirtual |
Initialize the state of filter and equation before each iteration. Progress feeback is implemented as part of this method.
Reimplemented from itk::FiniteDifferenceImageFilter< TDisplacementField, TDisplacementField >.
Reimplemented in itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, itk::DiffeomorphicDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, itk::FastSymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, itk::DemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, and itk::SymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >.
|
static |
Method for creation through the object factory.
|
private |
|
protectedvirtual |
This method is called after the solution has been generated. In this case, the filter release the memory of the internal buffers.
Reimplemented from itk::FiniteDifferenceImageFilter< TDisplacementField, TDisplacementField >.
|
protectedvirtual |
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 >.
Reimplemented in itk::DiffeomorphicDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, itk::FastSymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >, and itk::SymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >.
void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::SetFixedImage | ( | const FixedImageType * | ptr | ) |
Set the fixed image.
|
inline |
Set initial displacement field.
Definition at line 138 of file itkPDEDeformableRegistrationFilter.h.
|
virtual |
Set/Get the desired maximum error of the Guassian kernel approximate.
|
virtual |
Set/Get the desired limits of the Gaussian kernel width.
void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField >::SetMovingImage | ( | const MovingImageType * | ptr | ) |
Set the moving image.
|
virtual |
Set/Get whether the displacement field is smoothed (regularized). Smoothing the displacement yields a solution elastic in nature. If SmoothDisplacementField is on, then the displacement field is smoothed with a Gaussian whose standard deviations are specified with SetStandardDeviations()
|
virtual |
Set/Get whether the update field is smoothed (regularized). Smoothing the update field yields a solution viscous in nature. If SmoothUpdateField is on, then the update field is smoothed with a Gaussian whose standard deviations are specified with SetUpdateFieldStandardDeviations()
|
virtual |
Set the Gaussian smoothing standard deviations for the displacement field. The values are set with respect to pixel coordinates.
|
virtual |
Set the Gaussian smoothing standard deviations for the displacement field. The values are set with respect to pixel coordinates.
|
virtual |
Set the Gaussian smoothing standard deviations for the update field. The values are set with respect to pixel coordinates.
|
virtual |
Set the Gaussian smoothing standard deviations for the update field. The values are set with respect to pixel coordinates.
|
protectedvirtual |
Utility to smooth the displacement field (represented in the Output) using a Guassian operator. The amount of smoothing can be specified by setting the StandardDeviations.
|
virtual |
Set/Get whether the displacement field is smoothed (regularized). Smoothing the displacement yields a solution elastic in nature. If SmoothDisplacementField is on, then the displacement field is smoothed with a Gaussian whose standard deviations are specified with SetStandardDeviations()
|
virtual |
Set/Get whether the displacement field is smoothed (regularized). Smoothing the displacement yields a solution elastic in nature. If SmoothDisplacementField is on, then the displacement field is smoothed with a Gaussian whose standard deviations are specified with SetStandardDeviations()
|
protectedvirtual |
Utility to smooth the UpdateBuffer using a Gaussian operator. The amount of smoothing can be specified by setting the UpdateFieldStandardDeviations.
|
virtual |
Set/Get whether the update field is smoothed (regularized). Smoothing the update field yields a solution viscous in nature. If SmoothUpdateField is on, then the update field is smoothed with a Gaussian whose standard deviations are specified with SetUpdateFieldStandardDeviations()
|
virtual |
Set/Get whether the update field is smoothed (regularized). Smoothing the update field yields a solution viscous in nature. If SmoothUpdateField is on, then the update field is smoothed with a Gaussian whose standard deviations are specified with SetUpdateFieldStandardDeviations()
|
inlinevirtual |
Stop the registration after the current iteration.
Definition at line 228 of file itkPDEDeformableRegistrationFilter.h.
|
static |
Inherit some enums and typedefs from the superclass.
Definition at line 123 of file itkPDEDeformableRegistrationFilter.h.
|
private |
Maximum error for Gaussian operator approximation.
Definition at line 323 of file itkPDEDeformableRegistrationFilter.h.
|
private |
Limits of Guassian kernel width.
Definition at line 326 of file itkPDEDeformableRegistrationFilter.h.
|
private |
Modes to control smoothing of the update and displacement fields
Definition at line 314 of file itkPDEDeformableRegistrationFilter.h.
|
private |
Definition at line 315 of file itkPDEDeformableRegistrationFilter.h.
|
private |
Standard deviation for Gaussian smoothing
Definition at line 310 of file itkPDEDeformableRegistrationFilter.h.
|
private |
Flag to indicate user stop registration request.
Definition at line 329 of file itkPDEDeformableRegistrationFilter.h.
|
private |
Temporary displacement field use for smoothing the the displacement field.
Definition at line 319 of file itkPDEDeformableRegistrationFilter.h.
|
private |
Definition at line 311 of file itkPDEDeformableRegistrationFilter.h.