Computing large deformation metric mappings via geodesic flows of diffeomorphisms

Abstract

This paper examine the Euler-Lagrange equations for the solution of the large deformation diffeomorphic metric mapping problem studied in Dupuis et al. (1998) and Trouvé (1995) in which two images I0, I1 are given and connected via the diffeomorphic change of coordinates I 0 o φ-1 = I1 where φ = φ1 is the end point at t = 1 of curve φt, t ∈ [0, 1] satisfying φt, = vt(φt), t ∈ [0, 1] with φ0 = id. The variational problem takes the form argmin ν: φ̇t=ν1(φ1) (∫01 ∥νt∥V2dt + ∥I0o φ1-1 - I 1∥L22 where ∥νt∥V is an appropriate Sobolev norm on the velocity field νt(·), and the second term enforces matching of the images with ∥·∥L2 representing the squared-error norm. In this paper we derive the Euler-Lagrange equations characterizing the minimizing vector fields νt, t ∈ [0, 1] assuming sufficient smoothness of the norm to guarantee existence of solutions in the space of diffeomorphisms. We describe the implementation of the Euler equations using semi-lagrangian method of computing particle flows and show the solutions for various examples. As well, we compute the metric distance on several anatomical configurations as measured by ∫01 ∥νt∥V dt on the geodesic shortest paths.