template<typename TFixedImage, typename TMovingImage, typename TDisplacementField, typename TImageForceFunction = MeanSquareRegistrationFunction< TFixedImage, TMovingImage, TDisplacementField >>
class itk::CurvatureRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField, TImageForceFunction >
Deformably register two images using the fast curvature algorithm.
CurvatureRegistrationFilter implements the fast (i.e., O(n log n) ) registration method described in B. Fischer and J. Modersitzki, "A unified approach to fast image registration and a new curvature
based registration technique," Linear Algebra and its Applications, vol. 380, pp. 107-124, 2004.
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, the deformation field type, and the images forces function. Fundamentally, any image force function in the PDE framework should work (e.g., mean squares, demons, correlation, mutual information, etc.).
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 output deformation field can be obtained via methods GetOutput or GetDisplacementField.
This class makes use of the finite difference solver hierarchy. Update for each iteration is computed in the function class defined by the template parameter TImageForceFunction.
The algorithm has three parameters: the number of iteration to be performed, the time step, and the weight of the curvature regularization term in the overall cost function.
- Warning
- The algorithm parameters (time step and constraint weight) as given by F&M (dt=100; alpha=0.01) do not seem to work with this implementation. In fact, they seem to be off by orders of magnitude. This may be an effect of different normalization of the image forces, or this implementation may be incorrect. Be aware.
-
This filter assumes that the fixed image type, moving image type and deformation field type all have the same number of dimensions.
-
There is something sketchy going on with the DCT (see FFTW documentation) regarding numerical stability of the "ordinary" DCT. Also, F&M are slightly ambiguous in their use of the DCT operator. So it is quite possible that the implementation of their algorithm is not quite correct. Ultimately, they are saying that the curvature smoother is equivalent to a special low-pass filter that is applied to the DCT coefficients.
- Note
- This class was developed with funding from:
"CNS Deficits: Interaction of Age and Alcoholism" NIAAA AA05965, PI: A. Pfefferbaum
"INIA: Imaging Core" NIAAA AA13521, PI: A. Pfefferbaum
- See Also
- CurvatureRegistrationFunction
- Author
- Torsten Rohlfing, SRI International, Neuroscience Program
- Examples:
- Examples/RegistrationITKv3/DeformableRegistration10.cxx, Examples/RegistrationITKv3/DeformableRegistration9.cxx, Examples/RegistrationITKv4/DeformableRegistration10.cxx, and Examples/RegistrationITKv4/DeformableRegistration9.cxx.
Definition at line 100 of file itkCurvatureRegistrationFilter.h.
|
using | ConstPointer = SmartPointer< const Self > |
|
using | DisplacementFieldComponentImagePointer = typename DisplacementFieldComponentImageType::Pointer |
|
using | DisplacementFieldComponentImageType = Image< RealTypeDFT, TDisplacementField::ImageDimension > |
|
using | DisplacementFieldComponentType = typename DisplacementFieldPixelType::ValueType |
|
using | DisplacementFieldPixelType = typename TDisplacementField::PixelType |
|
using | DisplacementFieldPointer = typename Superclass::DisplacementFieldPointer |
|
using | DisplacementFieldType = typename Superclass::DisplacementFieldType |
|
using | FiniteDifferenceFunctionType = typename Superclass::FiniteDifferenceFunctionType |
|
using | FixedImagePointer = typename Superclass::FixedImagePointer |
|
using | FixedImageType = typename Superclass::FixedImageType |
|
using | MovingImagePointer = typename Superclass::MovingImagePointer |
|
using | MovingImageType = typename Superclass::MovingImageType |
|
using | Pointer = SmartPointer< Self > |
|
using | RealTypeDFT = double |
|
using | RegistrationFunctionType = TImageForceFunction |
|
using | Self = CurvatureRegistrationFilter |
|
using | Superclass = PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDisplacementField > |
|
using | TimeStepType = typename Superclass::TimeStepType |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | DisplacementFieldPointer = typename DisplacementFieldType::Pointer |
|
using | DisplacementFieldType = TDisplacementField |
|
using | FiniteDifferenceFunctionType = typename Superclass::FiniteDifferenceFunctionType |
|
using | FixedImageConstPointer = typename FixedImageType::ConstPointer |
|
using | FixedImagePointer = typename FixedImageType::Pointer |
|
using | FixedImageType = TFixedImage |
|
using | MovingImageConstPointer = typename MovingImageType::ConstPointer |
|
using | MovingImagePointer = typename MovingImageType::Pointer |
|
using | MovingImageType = TMovingImage |
|
using | OutputImageType = typename Superclass::OutputImageType |
|
using | PDEDeformableRegistrationFunctionType = PDEDeformableRegistrationFunction< FixedImageType, MovingImageType, DisplacementFieldType > |
|
using | Pointer = SmartPointer< Self > |
|
using | Self = PDEDeformableRegistrationFilter |
|
using | StandardDeviationsType = FixedArray< double, ImageDimension > |
|
using | Superclass = DenseFiniteDifferenceImageFilter< TDisplacementField, TDisplacementField > |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | FiniteDifferenceFunctionType = typename Superclass::FiniteDifferenceFunctionType |
|
using | InputImageType = typename Superclass::InputImageType |
|
using | OutputImageType = typename Superclass::OutputImageType |
|
using | PixelType = typename Superclass::PixelType |
|
using | Pointer = SmartPointer< Self > |
|
using | Self = DenseFiniteDifferenceImageFilter |
|
using | Superclass = FiniteDifferenceImageFilter< TDisplacementField, TDisplacementField > |
|
using | TimeStepType = typename Superclass::TimeStepType |
|
using | UpdateBufferType = OutputImageType |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | FiniteDifferenceFunctionType = FiniteDifferenceFunction< TDisplacementField > |
|
using | InputImageType = TDisplacementField |
|
using | InputPixelType = typename TDisplacementField::PixelType |
|
using | InputPixelValueType = typename NumericTraits< InputPixelType >::ValueType |
|
using | NeighborhoodScalesType = typename FiniteDifferenceFunctionType::NeighborhoodScalesType |
|
using | OutputImageType = TDisplacementField |
|
using | OutputPixelType = typename TDisplacementField::PixelType |
|
using | OutputPixelValueType = typename NumericTraits< OutputPixelType >::ValueType |
|
using | PixelType = OutputPixelType |
|
using | Pointer = SmartPointer< Self > |
|
using | RadiusType = typename FiniteDifferenceFunctionType::RadiusType |
|
using | Self = FiniteDifferenceImageFilter |
|
using | Superclass = InPlaceImageFilter< TDisplacementField, TDisplacementField > |
|
using | TimeStepType = typename FiniteDifferenceFunctionType::TimeStepType |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | InputImageConstPointer = typename InputImageType::ConstPointer |
|
using | InputImagePixelType = typename InputImageType::PixelType |
|
using | InputImagePointer = typename InputImageType::Pointer |
|
using | InputImageRegionType = typename InputImageType::RegionType |
|
using | InputImageType = TDisplacementField |
|
using | OutputImagePixelType = typename Superclass::OutputImagePixelType |
|
using | OutputImagePointer = typename Superclass::OutputImagePointer |
|
using | OutputImageRegionType = typename Superclass::OutputImageRegionType |
|
using | OutputImageType = typename Superclass::OutputImageType |
|
using | Pointer = SmartPointer< Self > |
|
using | Self = InPlaceImageFilter |
|
using | Superclass = ImageToImageFilter< TDisplacementField, TDisplacementField > |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | InputImageConstPointer = typename InputImageType::ConstPointer |
|
using | InputImagePixelType = typename InputImageType::PixelType |
|
using | InputImagePointer = typename InputImageType::Pointer |
|
using | InputImageRegionType = typename InputImageType::RegionType |
|
using | InputImageType = TDisplacementField |
|
using | OutputImagePixelType = typename Superclass::OutputImagePixelType |
|
using | OutputImageRegionType = typename Superclass::OutputImageRegionType |
|
using | Pointer = SmartPointer< Self > |
|
using | Self = ImageToImageFilter |
|
using | Superclass = ImageSource< TDisplacementField > |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | DataObjectIdentifierType = Superclass::DataObjectIdentifierType |
|
using | DataObjectPointer = DataObject::Pointer |
|
using | DataObjectPointerArraySizeType = Superclass::DataObjectPointerArraySizeType |
|
using | OutputImagePixelType = typename OutputImageType::PixelType |
|
using | OutputImagePointer = typename OutputImageType::Pointer |
|
using | OutputImageRegionType = typename OutputImageType::RegionType |
|
using | OutputImageType = TDisplacementField |
|
using | Pointer = SmartPointer< Self > |
|
using | Self = ImageSource |
|
using | Superclass = ProcessObject |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | DataObjectIdentifierType = DataObject::DataObjectIdentifierType |
|
using | DataObjectPointer = DataObject::Pointer |
|
using | DataObjectPointerArray = std::vector< DataObjectPointer > |
|
using | DataObjectPointerArraySizeType = DataObjectPointerArray::size_type |
|
using | MultiThreaderType = MultiThreaderBase |
|
using | NameArray = std::vector< DataObjectIdentifierType > |
|
using | Pointer = SmartPointer< Self > |
|
using | Self = ProcessObject |
|
using | Superclass = Object |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | Pointer = SmartPointer< Self > |
|
using | Self = Object |
|
using | Superclass = LightObject |
|
using | ConstPointer = SmartPointer< const Self > |
|
using | Pointer = SmartPointer< Self > |
|
using | Self = LightObject |
|
|
void | ApplyUpdate (const TimeStepType &dt) override |
|
| CurvatureRegistrationFilter () |
|
void | Initialize () override |
|
void | PrintSelf (std::ostream &os, Indent indent) const override |
|
| ~CurvatureRegistrationFilter () override |
|
void | CopyInputToOutput () override |
|
void | GenerateInputRequestedRegion () override |
|
void | GenerateOutputInformation () override |
|
bool | Halt () override |
|
void | Initialize () override |
|
void | InitializeIteration () override |
|
| PDEDeformableRegistrationFilter () |
|
void | PostProcessOutput () override |
|
void | PrintSelf (std::ostream &os, Indent indent) const override |
|
virtual void | SmoothDisplacementField () |
|
virtual void | SmoothUpdateField () |
|
| ~PDEDeformableRegistrationFilter () override=default |
|
void | AllocateUpdateBuffer () override |
|
void | ApplyUpdate (const TimeStepType &dt) override |
|
TimeStepType | CalculateChange () override |
|
void | CopyInputToOutput () override |
|
| DenseFiniteDifferenceImageFilter () |
|
virtual UpdateBufferType * | GetUpdateBuffer () |
|
void | PrintSelf (std::ostream &os, Indent indent) const override |
|
virtual void | ThreadedApplyUpdate (const TimeStepType &dt, const ThreadRegionType ®ionToProcess, ThreadIdType threadId) |
|
virtual TimeStepType | ThreadedCalculateChange (const ThreadRegionType ®ionToProcess, ThreadIdType threadId) |
|
| ~DenseFiniteDifferenceImageFilter () override=default |
|
virtual void | ApplyUpdate (const TimeStepType &dt)=0 |
|
| FiniteDifferenceImageFilter () |
|
void | GenerateData () override |
|
void | GenerateInputRequestedRegion () override |
|
void | PrintSelf (std::ostream &os, Indent indent) const override |
|
virtual TimeStepType | ResolveTimeStep (const std::vector< TimeStepType > &timeStepList, const std::vector< bool > &valid) const |
|
virtual void | SetElapsedIterations (IdentifierType _arg) |
|
virtual bool | ThreadedHalt (void *) |
|
| ~FiniteDifferenceImageFilter () override=default |
|
void | AllocateOutputs () override |
|
virtual bool | GetRunningInPlace () const |
|
| InPlaceImageFilter ()=default |
|
void | ReleaseInputs () override |
|
| ~InPlaceImageFilter () override=default |
|
virtual void | CallCopyInputRegionToOutputRegion (OutputImageRegionType &destRegion, const InputImageRegionType &srcRegion) |
|
virtual void | CallCopyOutputRegionToInputRegion (InputImageRegionType &destRegion, const OutputImageRegionType &srcRegion) |
|
| ImageToImageFilter () |
|
void | VerifyInputInformation () ITKv5_CONST override |
|
| ~ImageToImageFilter () override |
|
void | PushBackInput (const DataObject *input) override |
|
void | PushFrontInput (const DataObject *input) override |
|
virtual void | AfterThreadedGenerateData () |
|
virtual void | BeforeThreadedGenerateData () |
|
void | ClassicMultiThread (ThreadFunctionType callbackFunction) |
|
virtual const
ImageRegionSplitterBase * | GetImageRegionSplitter () const |
|
| ImageSource () |
|
virtual unsigned int | SplitRequestedRegion (unsigned int i, unsigned int pieces, OutputImageRegionType &splitRegion) |
|
| ~ImageSource () override=default |
|
virtual void | ThreadedGenerateData (const OutputImageRegionType &outputRegionForThread, ThreadIdType threadId) |
|
virtual void | DynamicThreadedGenerateData (const OutputImageRegionType &outputRegionForThread) |
|
virtual bool | GetDynamicMultiThreading () const |
|
virtual void | SetDynamicMultiThreading (bool _arg) |
|
virtual void | DynamicMultiThreadingOn () |
|
virtual void | DynamicMultiThreadingOff () |
|
virtual void | AddInput (DataObject *input) |
|
void | AddOptionalInputName (const DataObjectIdentifierType &) |
|
void | AddOptionalInputName (const DataObjectIdentifierType &, DataObjectPointerArraySizeType idx) |
|
virtual void | AddOutput (DataObject *output) |
|
bool | AddRequiredInputName (const DataObjectIdentifierType &) |
|
bool | AddRequiredInputName (const DataObjectIdentifierType &, DataObjectPointerArraySizeType idx) |
|
virtual void | CacheInputReleaseDataFlags () |
|
virtual void | GenerateOutputRequestedRegion (DataObject *output) |
|
DataObject * | GetInput (const DataObjectIdentifierType &key) |
|
const DataObject * | GetInput (const DataObjectIdentifierType &key) const |
|
virtual const
DataObjectPointerArraySizeType & | GetNumberOfRequiredInputs () const |
|
virtual const
DataObjectPointerArraySizeType & | GetNumberOfRequiredOutputs () const |
|
bool | IsIndexedInputName (const DataObjectIdentifierType &) const |
|
bool | IsIndexedOutputName (const DataObjectIdentifierType &) const |
|
bool | IsRequiredInputName (const DataObjectIdentifierType &) const |
|
DataObjectPointerArraySizeType | MakeIndexFromInputName (const DataObjectIdentifierType &name) const |
|
DataObjectPointerArraySizeType | MakeIndexFromOutputName (const DataObjectIdentifierType &name) const |
|
DataObjectIdentifierType | MakeNameFromInputIndex (DataObjectPointerArraySizeType idx) const |
|
DataObjectIdentifierType | MakeNameFromOutputIndex (DataObjectPointerArraySizeType idx) const |
|
| ProcessObject () |
|
virtual void | PropagateResetPipeline () |
|
virtual void | RemoveInput (const DataObjectIdentifierType &key) |
|
virtual void | RemoveInput (DataObjectPointerArraySizeType) |
|
virtual void | RemoveOutput (const DataObjectIdentifierType &key) |
|
virtual void | RemoveOutput (DataObjectPointerArraySizeType idx) |
|
bool | RemoveRequiredInputName (const DataObjectIdentifierType &) |
|
virtual void | RestoreInputReleaseDataFlags () |
|
virtual void | SetInput (const DataObjectIdentifierType &key, DataObject *input) |
|
virtual void | SetNthInput (DataObjectPointerArraySizeType num, DataObject *input) |
|
virtual void | SetNthOutput (DataObjectPointerArraySizeType num, DataObject *output) |
|
void | SetNumberOfIndexedInputs (DataObjectPointerArraySizeType num) |
|
void | SetNumberOfIndexedOutputs (DataObjectPointerArraySizeType num) |
|
virtual void | SetNumberOfRequiredInputs (DataObjectPointerArraySizeType) |
|
virtual void | SetNumberOfRequiredOutputs (DataObjectPointerArraySizeType _arg) |
|
virtual void | SetOutput (const DataObjectIdentifierType &key, DataObject *output) |
|
virtual void | SetPrimaryInput (DataObject *input) |
|
virtual void | SetPrimaryOutput (DataObject *output) |
|
void | SetRequiredInputNames (const NameArray &) |
|
virtual void | VerifyPreconditions () ITKv5_CONST |
|
| ~ProcessObject () override |
|
DataObject * | GetInput (DataObjectPointerArraySizeType idx) |
|
const DataObject * | GetInput (DataObjectPointerArraySizeType idx) const |
|
DataObject * | GetPrimaryInput () |
|
const DataObject * | GetPrimaryInput () const |
|
virtual void | SetPrimaryInputName (const DataObjectIdentifierType &key) |
|
virtual const char * | GetPrimaryInputName () const |
|
DataObject * | GetOutput (const DataObjectIdentifierType &key) |
|
const DataObject * | GetOutput (const DataObjectIdentifierType &key) const |
|
virtual void | SetPrimaryOutputName (const DataObjectIdentifierType &key) |
|
virtual const char * | GetPrimaryOutputName () const |
|
DataObject * | GetOutput (DataObjectPointerArraySizeType idx) |
|
const DataObject * | GetOutput (DataObjectPointerArraySizeType idx) const |
|
DataObject * | GetPrimaryOutput () |
|
const DataObject * | GetPrimaryOutput () const |
|
| Object () |
|
bool | PrintObservers (std::ostream &os, Indent indent) const |
|
virtual void | SetTimeStamp (const TimeStamp &time) |
|
| ~Object () override |
|
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 () |
|