Stratified mixture modeling for segmentation of white-matter lesions in brain MR images

Accurate characterization of white-matter lesions from magnetic resonance (MR) images has increasing importance for diagnosis and management of treatment of certain neurological diseases, and can be performed in an objective and effective way by automated lesion segmentation. This usually involves modeling the whole-brain MR intensity distribution, however, capturing various sources of MR intensity variability and lesion heterogeneity results in highly complex whole-brain MR intensity models, thus their robust estimation on a large set of MR images presents a huge challenge. We propose a novel approach employing stratified mixture modeling, where the main premise is that the otherwise complex whole-brain model can be reduced to a tractable parametric form in small brain subregions. We show on MR images of multiple sclerosis (MS) patients with different lesion loads that robust estimators enable accurate mixture modeling of MR intensity in small brain subregions even in the presence of lesions. Recombination of the mixture models across strata provided an accurate whole-brain MR intensity model. Increasing the number of subregions and, thereby, the model complexity, consistently improved the accuracy of whole-brain MR intensity modeling and segmentation of normal structures. The proposed approach was incorporated into three unsupervised lesion segmentation methods and, compared to original and three other state-of-the-art methods, the proposed modeling approach significantly improved lesion segmentation according to increased Dice similarity indices and lower number of false positives on real MR images of 30 patients with MS.