Watershed-based segmentation of 3D MR data for volume quantization

The aim of this work is the development of a semiautomatic segmentation technique for efficient and accurate volume quantization of Magnetic Resonance (MR) data. The proposed technique uses a 3D variant of Vincent and Soilles immersion-based watershed algorithm that is applied to the gradient magnitude of the MR data and that produces small volume primitives. The known drawback of the watershed algorithm, oversegmentation, is strongly reduced by a priori application of a 3D adaptive anisotropic diffusion filter to the MR data. Furthermore, oversegmentation is a posteriori reduced by properly merging small volume primitives that have similar gray level distributions. The outcome of the preceding image processing steps is presented to the user for manual segmentation. Through selection of volume primitives, the user quickly segments the first slice, which contains the object of interest. Afterwards, the subsequent slices are automatically segmented by extrapolation. Segmentation results are contingently manually corrected. The proposed segmentation technique is tested on phantom objects, where segmentation errors less than 2% are observed. In addition, the technique is demonstrated on 3D MR data of the mouse head from which the cerebellum is extracted. Volumes of the mouse cerebellum and the mouse brains in toto are calculated.