A SAEM algorithm for the estimation of template and deformation parameters in medical image sequences

This paper is about object deformations observed throughout a sequence of images. We present a statistical framework in which the observed images are defined as noisy realizations of a randomly deformed template image. In this framework, we focus on the problem of the estimation of parameters relate...

Full description

Saved in:
Bibliographic Details
Published in:Statistics and computing 2009-12, Vol.19 (4), p.465-478
Main Authors: Richard, Frédéric J. P., Samson, Adeline M. M., Cuénod, Charles A.
Format: Article
Language:English
Subjects:
Algorithms
Artificial Intelligence
Computer Science
Computer simulation
Construction
Deformation
Engineering Sciences
Image Processing
Mathematical models
Mathematics and Statistics
Medical Imaging
Methodology
Monte Carlo methods
Probability and Statistics in Computer Science
Representations
Signal and Image processing
Statistical Theory and Methods
Statistics
Statistics and Computing/Statistics Programs
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:This paper is about object deformations observed throughout a sequence of images. We present a statistical framework in which the observed images are defined as noisy realizations of a randomly deformed template image. In this framework, we focus on the problem of the estimation of parameters related to the template and deformations. Our main motivation is the construction of estimation framework and algorithm which can be applied to short sequences of complex and highly-dimensional images. The originality of our approach lies in the representations of the template and deformations, which are defined on a common triangulated domain, adapted to the geometry of the observed images. In this way, we have joint representations of the template and deformations which are compact and parsimonious. Using such representations, we are able to drastically reduce the number of parameters in the model. Besides, we adapt to our framework the Stochastic Approximation EM algorithm combined with a Markov Chain Monte Carlo procedure which was proposed in 2004 by Kuhn and Lavielle. Our implementation of this algorithm takes advantage of some properties which are specific to our framework. More precisely, we use the Markovian properties of deformations to build an efficient simulation strategy based on a Metropolis-Hasting-Within-Gibbs sampler. Finally, we present some experiments on sequences of medical images and synthetic data.
ISSN:0960-3174
1573-1375
DOI:10.1007/s11222-008-9106-7