ATLAS CONSTRUCTION FROM HIGH ANGULAR RESOLUTION DIFFUSION IMAGING DATA REPRESENTED BY GAUSSIAN MIXTURE FIELDS

Reference-free brain template construction with population symmetric registration

Population registration has been proposed for normalizing a large group of images into a common space, which is important in many clinical and research studies, such as brain development, aging, and atlas construction. Different from pairwise registration problem that aligns the target image to the reference directly, determining the reference or the hidden common space with the least bias is important in population registration. In order to decrease this bias, a lot of work takes the arithmetic mean image as the reference. However, the arithmetic mean image is usually too smooth to guide the population registration. This work presents an efficient symmetric population registration strategy for brain template construction, which defines the symmetric population center guiding population registration. This is important because the population registration problem can be translated into a series of pairwise registration problem which is easier to optimize and implement. Another prominent merit of proposed population registration algorithm is reference-free, which eliminates the reference dependency-related bias in population registration. Based on this symmetric population registration, the brain template is constructed by approximating both the population's intensity and gradient information. In addition, we also present a new measurement named with average bias for evaluating the unbiasedness of brain template. Experiments were first carried out on four synthetic images created with controllable transforms, which aim at comparing the difference between conventional method and proposed method. Further experiment is designed for reference-free validation. Finally, in real inter-subject brain data, twenty MRI T1 volumes with size 256 × 256 × 176 are used to construct a symmetric brain template with proposed population registration method. The constructed brain template has a small bias and clear brain details comparing with DARTEL.

A Method for Automated Classification of Parkinson's Disease Diagnosis Using an Ensemble Average Propagator Template Brain Map Estimated from Diffusion MRI

Parkinson's disease (PD) is a common and debilitating neurodegenerative disorder that affects patients in all countries and of all nationalities. Magnetic resonance imaging (MRI) is currently one of the most widely used diagnostic imaging techniques utilized for detection of neurologic diseases. Changes in structural biomarkers will likely play an important future role in assessing progression of many neurological diseases inclusive of PD. In this paper, we derived structural biomarkers from diffusion MRI (dMRI), a structural modality that allows for non-invasive inference of neuronal fiber connectivity patterns. The structural biomarker we use is the ensemble average propagator (EAP), a probability density function fully characterizing the diffusion locally at a voxel level. To assess changes with respect to a normal anatomy, we construct an unbiased template brain map from the EAP fields of a control population. Use of an EAP captures both orientation and shape information of the diffusion process at each voxel in the dMRI data, and this feature can be a powerful representation to achieve enhanced PD brain mapping. This template brain map construction method is applicable to small animal models as well as to human brains. The differences between the control template brain map and novel patient data can then be assessed via a nonrigid warping algorithm that transforms the novel data into correspondence with the template brain map, thereby capturing the amount of elastic deformation needed to achieve this correspondence. We present the use of a manifold-valued feature called the Cauchy deformation tensor (CDT), which facilitates morphometric analysis and automated classification of a PD versus a control population. Finally, we present preliminary results of automated discrimination between a group of 22 controls and 46 PD patients using CDT. This method may be possibly applied to larger population sizes and other parkinsonian syndromes in the near future.

Tractography from HARDI using an intrinsic unscented Kalman filter

A novel adaptation of the unscented Kalman filter (UKF) was recently introduced in literature for simultaneous multitensor estimation and fiber tractography from diffusion MRI. This technique has the advantage over other tractography methods in terms of computational efficiency, due to the fact that the UKF simultaneously estimates the diffusion tensors and propagates the most consistent direction to track along. This UKF and its variants reported later in literature however are not intrinsic to the space of diffusion tensors. Lack of this key property can possibly lead to inaccuracies in the multitensor estimation as well as in the tractography. In this paper, we propose a novel intrinsic unscented Kalman filter (IUKF) in the space of diffusion tensors which are symmetric positive definite matrices, that can be used for simultaneous recursive estimation of multitensors and propagation of directional information for use in fiber tractography from diffusion weighted MR data. In addition to being more accurate, IUKF retains all the advantages of UKF mentioned above. We demonstrate the accuracy and effectiveness of the proposed method via experiments publicly available phantom data from the fiber cup-challenge (MICCAI 2009) and diffusion weighted MR scans acquired from human brains and rat spinal cords.

iPGA: incremental principal geodesic analysis with applications to movement disorder classification

The nonlinear version of the well known PCA called the Prinicipal Geodesic Analysis (PGA) was introduced in the past decade for statistical analysis of shapes as well as diffusion tensors. PGA of diffusion tensor fields or any other manifold-valued fields can be a computationally demanding task due to the dimensionality of the problem and thus establishing motivation for an incremental PGA (iPGA) algorithm. In this paper, we present a novel iPGA algorithm that incrementally updates the current Karcher mean and the principal sub-manifolds with any newly introduced data into the pool without having to recompute the PGA from scratch. We demonstrate substantial computational and memory savings of iPGA over the batch mode PGA for diffusion tensor fields via synthetic and real data examples. Further, we use the iPGA derived representation in an NN classifier to automatically discriminate between controls, Parkinson's Disease and Essential Tremor patients, given their HARDI brain scans.

A NOVEL INTRINSIC UNSCENTED KALMAN FILTER FOR TRACTOGRAPHY FROM HARDI*

The unscented Kalman filter (UKF) was recently introduced in literature for simultaneous multi-tensor estimation and tractography. This UKF however was not intrinsic to the space of diffusion tensors. Lack of this key property leads to inaccuracies in the multi-tensor estimation as well as in tractography. In this paper, we propose an novel intrinsic unscented Kalman filter (IUKF) in the space of symmetric positive definite matrices, which can be used for simultaneous recursive estimation of multi-tensors and tractography from diffusion weighted MR data. In addition to being more accurate, IUKF retains all the advantages of UKF for instance, multi-tensor estimation is only performed in the places where it is needed for tractography, which would be much more efficient than the two stage process involved in methods that do tracking post diffusion tensor estimation. The accuracy and effectiveness of the proposed method is demonstrated via real data experiments.