Bayesian approach to segmentation of statistical parametric maps

A contextual segmentation technique to detect brain activation from functional brain images is presented in the Bayesian framework. Unlike earlier similar approaches [Holmes and Ford (1993) and Descombes et al. (1998)], a Markov random field (MRF) is used to represent configurations of activated brain voxels, and likelihoods given by statistical parametric maps (SPM's) are directly used to find the maximum a posteriori (MAP) estimation of segmentation. The iterative segmentation algorithm, which is based on a simulated annealing scheme, is fully data-driven and capable of analyzing experiments involving multiple-input stimuli. Simulation results and comparisons with the simple thresholding and the statistical parametric mapping (SPM) approaches are presented with synthetic images, and functional MR images acquired in memory retrieval and event-related working memory tasks. The experiments show that an MRF is a valid representation of the activation patterns obtained in functional brain images, and the present technique renders a superior segmentation scheme to the context-free approach and the SPM approach.