Research Article

Diffeomorphic Registration of Images with Variable Contrast Enhancement

Figure 1

Comparison between additive and compositive field accumulations. Warping is implemented using linear interpolation. In (a) and (b), two different displacement fields are defined on the plane (for visual clarity). In (c), the field 𝐷 1 warped by 𝐷 2 , that is, 𝐷 1 𝐷 2 = 𝐷 1 Δ 2 . In (d), the field 𝐷 2 is applied on the pixel grid. In (e), the grid is warped by the field resulting from an addition-based accumulation of 𝐷 1 and 𝐷 2 . In (f), the grid is warped by the displacement field 𝐷 1 𝐷 2 arising by the composition of Δ 1 and Δ 2 , which is the sum of the dark blue and gray arrows (given by Δ 1 Δ 2 I d ). This composition is really the accumulation that matters since it corresponds to the way displacement fields are iteratively applied to an image (see Section 3). Since 𝐷 1 𝐷 2 = 𝐷 2 + 𝐷 1 𝐷 2 , the summed vectors in (f) correspond to the vectors in (a) and (c). In (g), a moving object 𝑚 , divided in 4 colors (regions between pixel centers). In (h), the result of the warping of 𝑚 by the sum of the fields. Clearly, the surfaces are inverted (mirror effect, visible because of the inversion of colors), leading to nonphysical deformations (negative Jacobians). In (i), the result of the warping of 𝑚 by the composition of the fields. One can notice that in spite of the deformation of the shape of the object, the location of the colors is conserved.
891585.fig.001a
(a) 𝐷 2
891585.fig.001b
(b) 𝐷 1
891585.fig.001c
(c) 𝐷 1 𝐷 2
891585.fig.001d
(d) 𝑥 + 𝐷 2
891585.fig.001e
(e) 𝑥 + ( 𝐷 1 + 𝐷 2 )
891585.fig.001f
(f) 𝑥 + ( 𝐷 1 𝐷 2 )
891585.fig.001g
(g) moving image 𝑚
891585.fig.001h
(h) 𝑚 ( 𝐷 1 + 𝐷 2 )
891585.fig.001i
(i) 𝑚 ( 𝐷 1 𝐷 2 )