Iglesias, Juan Eugenio3; Sabuncu, Mert Rory3; Van Leemput, Koen1
1 Department of Applied Mathematics and Computer Science, Technical University of Denmark2 Image Analysis & Computer Graphics, Department of Applied Mathematics and Computer Science, Technical University of Denmark3 Massachusetts General Hospital
Many segmentation algorithms in medical image analysis use Bayesian modeling to augment local image appearance with prior anatomical knowledge. Such methods often contain a large number of free parameters that are first estimated and then kept fixed during the actual segmentation process. However, a faithful Bayesian analysis would marginalize over such parameters, accounting for their uncertainty by considering all possible values they may take. Here we propose to incorporate this uncertainty into Bayesian segmentation methods in order to improve the inference process. In particular, we approximate the required marginalization over model parameters using computationally efficient Markov chain Monte Carlo techniques. We illustrate the proposed approach using a recently developed Bayesian method for the segmentation of hippocampal subfields in brain MRI scans, showing a significant improvement in an Alzheimer’s disease classification task. As an additional benefit, the technique also allows one to compute informative “error bars” on the volume estimates of individual structures.
Medical Image Analysis, 2013, Vol 17, Issue 7, p. 766-778
Bayesian modeling; Segmentation; Monte Carlo sampling; Hippocampal subfields
Main Research Area:
15th International Conference on Medical Image Computing and Computer Assisted InterventionMedical Image Computing and Computer Assisted Intervention, 2013