Improved inference in Bayesian segmentation using Monte Carlo sampling: application to hippocampal subfield volumetry

Med Image Anal. 2013 Oct;17(7):766-78. doi: 10.1016/j.media.2013.04.005. Epub 2013 May 22.

Abstract

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.

Keywords: Bayesian modeling; Hippocampal subfields; Monte Carlo sampling; Segmentation.

Publication types

  • Research Support, N.I.H., Extramural
  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms*
  • Alzheimer Disease / pathology*
  • Artificial Intelligence
  • Bayes Theorem
  • Data Interpretation, Statistical*
  • Hippocampus / pathology*
  • Humans
  • Image Enhancement / methods
  • Image Interpretation, Computer-Assisted / methods*
  • Imaging, Three-Dimensional / methods*
  • Monte Carlo Method
  • Organ Size
  • Pattern Recognition, Automated / methods*
  • Reproducibility of Results
  • Sample Size
  • Sensitivity and Specificity