Bayesian template estimation in computational anatomy

Templates play a fundamental role in Computational Anatomy. In this paper, we present a Bayesian model for template estimation. It is assumed that observed images I 1, I 2, …, I N are generated by shooting the template J through Gaussian distributed random initial momenta θ 1, θ 2, …, θ N. The template is J modeled as a deformation from a given hypertemplate J 0 with initial momentum μ, which has a Gaussian prior. We apply a mode approximation of the EM (MAEM) procedure, where the conditional expectation is replaced by a Dirac measure at the mode. This leads us to an image matching problem with a Jacobian weight term, and we solve it by deriving the weighted Euler–Lagrange equation. The results of template estimation for hippocampus and cardiac data are presented.

Authors

• Jun Ma University of Electronic Science and Technology of China
• Michael I. Miller
• Alain Trouvé
• Laurent Younes

Keywords

• Template estimation
• Computational anatomy
• Bayesian
• Weighted Euler–Lagrange equation