Inverse consistency

From HandWiki

In image registration, inverse consistency measures the consistency of mappings between images produced by a registration algorithm. The inverse consistency error, introduced by Christiansen and Johnson in 2001, quantifies the distance between the composition of the mappings from each image to the other, produced by the registration procedure, and the identity function, and is used as a regularisation constraint in the loss function of many registration algorithms to enforce consistent mappings.[1][2][3][4] Inverse consistency is necessary for good image registration but it is not sufficient, since a mapping can be perfectly consistent but not register the images at all.[5]

Definition

Image registration is the process of establishing a common coordinate system between two images, and given two images

[math]\displaystyle{ \begin{align} I_1: \Omega_1 \to \mathbb{R} \\ I_2: \Omega_2 \to \mathbb{R} \end{align} }[/math]

registering a source image [math]\displaystyle{ I_1 }[/math] to a target image [math]\displaystyle{ I_2 }[/math] consists of determining a transformation [math]\displaystyle{ f_1: \Omega_2 \to \Omega_1 }[/math] that maps points from the target space to the source space.[6] An ideal registration algorithm should not be sensitive to which image in the pair is used as source or target, and the registration operator should be antisymmetric such that the mappings

[math]\displaystyle{ \begin{align} f_1: \Omega_2 \to \Omega_1 \\ f_2: \Omega_1 \to \Omega_2 \end{align} }[/math]

produced when registering [math]\displaystyle{ I_1 }[/math] to [math]\displaystyle{ I_2 }[/math] and [math]\displaystyle{ I_2 }[/math] to [math]\displaystyle{ I_1 }[/math] respectively should be the inverse of each other, i.e. [math]\displaystyle{ f_2 = f_1^{-1} }[/math] and [math]\displaystyle{ f_1 = f_2^{-1} }[/math] or, equivalently, [math]\displaystyle{ f_2 \circ f_1 = \operatorname{id}_{\Omega_2} }[/math] and [math]\displaystyle{ f_1 \circ f_2 = \operatorname{id}_{\Omega_1} }[/math], where [math]\displaystyle{ \circ }[/math] denotes the function composition operator.[7][8]

Real algorithms are not perfect, and when swapping the role of source and target image in a registration problem the so obtained transformations are not the inverse of each other. Inverse consistency can be enforced by adding to the loss function of the registration a symmetric regularisation term that penalises inconsistent transformations[1]

[math]\displaystyle{ \int_{\Omega_2} \left\Vert f_2(f_1(x)) - x \right\Vert^2 \mathrm{d}x + \int_{\Omega_1} \left\Vert f_1(f_2(x)) - x \right\Vert^2 \mathrm{d}x . }[/math]

Inverse consistency can be used as a quality metric to evaluate image registration results. The inverse consistency error ([math]\displaystyle{ ICE }[/math]) measures the distance between the composition of the two transforms and the identity function, and it can be formulated in terms of both average ([math]\displaystyle{ ICE_a }[/math]) or maximum ([math]\displaystyle{ ICE_m }[/math]) over a region of interest [math]\displaystyle{ \Omega }[/math] of the image:[9]

[math]\displaystyle{ \begin{align} ICE_a &= \frac{1}{\int_{\Omega} \mathrm{d}x} \int_{\Omega} \left\Vert f_2(f_1(x)) - x \right\Vert \mathrm{d}x \\ ICE_m &= \max_{x \in \Omega} \left\Vert f_2(f_1(x)) - x \right\Vert . \end{align} }[/math]

While inverse consistency is a necessary property of good registration algorithms, inverse consistency error alone is not a sufficient metric to evaluate the quality of image registration results, since a perfectly consistent mapping, with no other constraint, may be not even close to correctly register a pair of images.[5]

References

  1. 1.0 1.1 Christiansen and Johnson (2001)
  2. Bhatia et al. (2004)
  3. Johnson and Christensen (2002)
  4. Arganda-Carreras et al. (2006)
  5. 5.0 5.1 Rohlfing (2012)
  6. The transformation goes in this direction in order to allow resampling the source image in target space.
  7. Bender and Tomé (2009)
  8. Leow et al. (2005)
  9. Beg and Khan (2007)

External links