A Riemannian approach to diffusion tensor images segmentation

We address the problem of the segmentation of cerebral white matter structures from diffusion tensor images. Our approach is grounded on the theoretically well-founded differential geometrical properties of the space of multivariate normal distributions. We introduce a variational formulation, in the level set framework, to estimate the optimal segmentation according to the following hypothesis: Diffusion tensors exhibit a Gaussian distribution in the different partitions. Moreover, we must respect the geometric constraints imposed by the interfaces existing among the cerebral structures and detected by the gradient of the diffusion tensor image. We validate our algorithm on synthetic data and report interesting results on real datasets. We focus on two structures of the white matter with different properties and respectively known as the corpus callosum and the corticospinal tract.