Proposals:Oriented Image Registration
itkOrientedImage should be the preferred image to use for performing image registration.
Currently the direction cosines contained in this image are not used consistently by all the classes involved in image registration.
This page proposes a comprehensive approach for fixing this problem.
Most image metrics internally compute image gradients as part of the process of computing the derivatives of the Metric. Such image gradients are currently computed in the space of the image grid, taking into account the pixel spacing but not taking into account the image direction
This issue has been logged as Bug Entry: 5081
and has been discussed multiple times in the developers list (citation needed)
The Proposed Solution
The prososed solution involves modifications of the toolkit at four different levels
- Adding a methods to itk::ImageBase for transforming index space Gradients into physical space Gradients and back
- These methods will be overloaded in the oriented image.
- The methods will be implemented using Template Metaprogramming in order to reduce their impact on performance
- Adding boolean flags to all classes that compute image gradients
- When the boolean flags are ON the image direction will be taken into account
- The boolean flags will be initialized to OFF by default in order to provide backward-compatibility
- The ImageMetrics will turn ON the boolean flags of helper classes used for computing gradients
- An advanced CMake variable will enable the use of all the mechanisms defined above
- The variable will be OFF by default in order to provide backward-compatibility
Details on each one of these items follow
Adding Gradient Transformation Methods to the itk::ImageBase
Following the naming of the existing methods
we propose to add the methods
or we could be more specific and say "IndexGradient"
to avoid the ambiguity of whether we are actually transforming indices into gradients, or transforming gradients that were computed in index space into gradient that are computed in physical space.
Note that strictly speaking, currently the gradients are not computed in index space, but in a space that takes the spacing but not the direction into account. This is not quite index space, but it is not physical space either.
We could acknowledge that fact and use a naming such as one of the following