Automatic segmentation of rodent spinal cord diffusion MR images

MRI, is a key tool for noninvasive spinal cord lesion analysis; however, accurate, quantitative methods for this analysis are lacking. A new, multistep, multidimensional approach, utilizing the classification expectation maximization algorithm, is proposed for MRI segmentation of spinal cord tissues. Diffusion tensor imaging is used to generate multiple images of each spinal slice, with different diffusion direction weightings. The maximum likelihood tissue classifications are then jointly estimated to produce a binary classification image, corresponding to voxels containing either spinal cord or background. Edge detection is employed to find a nonparametric curve encapsulating the entire spinal cord. The algorithm is evaluated using data from in vivo diffusion tensor imaging of control and injured mouse spinal cords. The algorithm is shown to remain accurate for whole spinal cord, white matter, and hemorrhage segmentation in the presence of significant injury. The results of the method are shown to be at least on par with expert manual segmentation.