FAST AND CLOSED-FORM ENSEMBLE-AVERAGE-PROPAGATOR APPROXIMATION FROM THE 4TH-ORDER DIFFUSION TENSOR

International audienceGeneralized Diffusion Tensor Imaging (GDTI) was developed to model complex Apparent Diffusivity Coefficient (ADC) using Higher Order Tensors (HOT) and to overcome the inherent single-peak shortcoming of DTI. However, the geometry of a complex ADC profile doesn't correspond to the underlying structure of fibers. This tissue geometry can be inferred from the shape of the Ensemble Average Propagator (EAP). Though interesting methods for estimating a positive ADC using 4th order diffusion tensors were developed, GDTI in general was overtaken by other approaches, e.g. the Orientation Distribution Function (ODF), since it is considerably difficult to recuperate the EAP from a HOT model of the ADC in GDTI. In this paper we present a novel closed-form approximation of the EAP using Hermite Polynomials from a modified HOT model of the original GDTI-ADC. Since the solution is analytical, it is fast, differentiable, and the approximation converges well to the true EAP. This method also makes the effort of computing a positive ADC worthwhile, since now both the ADC and the EAP can be used and have closed forms. We demonstrate on 4th order diffusion tensors

Diffeomorphic Metric Mapping and Probabilistic Atlas Generation of Hybrid Diffusion Imaging based on BFOR Signal Basis

We propose a large deformation diffeomorphic metric mapping algorithm to align multiple b-value diffusion weighted imaging (mDWI) data, specifically acquired via hybrid diffusion imaging (HYDI), denoted as LDDMM-HYDI. We then propose a Bayesian model for estimating the white matter atlas from HYDIs. We adopt the work given in Hosseinbor et al. (2012) and represent the q-space diffusion signal with the Bessel Fourier orientation reconstruction (BFOR) signal basis. The BFOR framework provides the representation of mDWI in the q-space and thus reduces memory requirement. In addition, since the BFOR signal basis is orthonormal, the L2 norm that quantifies the differences in the q-space signals of any two mDWI datasets can be easily computed as the sum of the squared differences in the BFOR expansion coefficients. In this work, we show that the reorientation of the $q$-space signal due to spatial transformation can be easily defined on the BFOR signal basis. We incorporate the BFOR signal basis into the LDDMM framework and derive the gradient descent algorithm for LDDMM-HYDI with explicit orientation optimization. Additionally, we extend the previous Bayesian atlas estimation framework for scalar-valued images to HYDIs and derive the expectation-maximization algorithm for solving the HYDI atlas estimation problem. Using real HYDI datasets, we show the Bayesian model generates the white matter atlas with anatomical details. Moreover, we show that it is important to consider the variation of mDWI reorientation due to a small change in diffeomorphic transformation in the LDDMM-HYDI optimization and to incorporate the full information of HYDI for aligning mDWI

Extraction of Structural Metrics from Crossing Fiber Models

Diffusion MRI (dMRI) measurements allow us to infer the microstructural properties of white matter and to reconstruct fiber pathways in-vivo. High angular diffusion imaging (HARDI) allows for the creation of more and more complex local models connecting the microstructure to the measured signal. One of the challenges is the derivation of meaningful metrics describing the underlying structure from the local models. The aim hereby is to increase the specificity of the widely used metric fractional anisotropy (FA) by using the additional information contained within the HARDI data. A local model which is connected directly to the underlying microstructure through the model of a single fiber population is spherical deconvolution. It produces a fiber orientation density function (fODF), which can often be interpreted as superposition of multiple peaks, each associated to one relatively coherent fiber population (bundle). Parameterizing these peaks one is able to disentangle and characterize these bundles. In this work, the fODF peaks are approximated by Bingham distributions, capturing first and second order statistics of the fiber orientations, from which metrics for the parametric quantification of fiber bundles are derived. Meaningful relationships between these measures and the underlying microstructural properties are proposed. The focus lies on metrics derived directly from properties of the Bingham distribution, such as peak length, peak direction, peak spread, integral over the peak, as well as a metric derived from the comparison of the largest peaks, which probes the complexity of the underlying microstructure. These metrics are compared to the conventionally used fractional anisotropy (FA) and it is shown how they may help to increase the specificity of the characterization of microstructural properties. Visualization of the micro-structural arrangement is another application of dMRI. This is done by using tractography to propagate the fiber layout, extracted from the local model, in each voxel. In practice most tractography algorithms use little of the additional information gained from HARDI based local models aside from the reconstructed fiber bundle directions. In this work an approach to tractography based on the Bingham parameterization of the fODF is introduced. For each of the fiber populations present in a voxel the diffusion signal and tensor are computed. Then tensor deflection tractography is performed. This allows incorporating the complete bundle information, performing local interpolation as well as using multiple directions per voxel for generating tracts. Another aspect of this work is the investigation of the spherical harmonic representation which is used most commonly for the fODF by means of the parameters derived from the Bingham distribution fit. Here a strong connection between the approximation errors in the spherical representation of the Dirac delta function and the distribution of crossing angles recovered from the fODF was discovered. The final aspect of this work is the application of the metrics derived from the Bingham fit to a number of fetal datasets for quantifying the brain’s development. This is done by introducing the Gini-coefficient as a metric describing the brain’s age

Estimation and Processing of Ensemble Average Propagator and Its Features in Diffusion MRI

Diffusion MRI (dMRI) is the unique technique to infer the microstructure of the white matter in vivo and noninvasively, by modeling the diffusion of water molecules. Ensemble Average Propagator (EAP) and Orientation Distribution Function (ODF) are two important Probability Density Functions (PDFs) which reflect the water diffusion. Estimation and processing of EAP and ODF is the central problem in dMRI, and is also the first step for tractography. Diffusion Tensor Imaging (DTI) is the most widely used estimation method which assumes EAP as a Gaussian distribution parameterized by a tensor. Riemannian framework for tensors has been proposed successfully in tensor estimation and processing. However, since the Gaussian EAP assumption is oversimplified, DTI can not reflect complex microstructure like fiber crossing. High Angular Resolution Diffusion Imaging (HARDI) is a category of methods proposed to avoid the limitations of DTI. Most HARDI methods like Q-Ball Imaging (QBI) need some assumptions and only can handle the data from single shell (single $b$ value), which are called as single shell HARDI (sHARDI) methods. However, with the development of scanners and acquisition methods, multiple shell data becomes more and more practical and popular. This thesis focuses on the estimation and processing methods in multiple shell HARDI (mHARDI) which can handle the diffusion data from arbitrary sampling scheme. There are many original contributions in this thesis. -First, we develop the analytical Spherical Polar Fourier Imaging (SPFI), which represents the signal using SPF basis and obtains EAP and its various features including ODFs and some scalar indices like Generalized Fractional Anisotropy (GFA) from analytical linear transforms. In the implementation of SPFI, we present two ways for scale estimation and propose to consider the prior $E(0)=1$ in estimation process. -Second, a novel Analytical Fourier Transform in Spherical Coordinate (AFT-SC) framework is proposed to incorporate many sHARDI and mHARDI methods, explore their relation and devise new analytical EAP/ODF estimation methods. -Third, we present some important criteria to compare different HARDI methods and illustrate their advantages and limitations. -Fourth, we propose a novel diffeomorphism invariant Riemannian framework for ODF and EAP processing, which is a natural generalization of previous Riemannian framework for tensors, and can be used for general PDF computing by representing the square root of the PDF called wavefunction with orthonormal basis. In this Riemannian framework, the exponential map, logarithmic map and geodesic have closed forms, the weighted Riemannian mean and median uniquely exist and can be estimated from an efficient gradient descent. Log-Euclidean framework and Affine-Euclidean framework are developed for fast data processing. -Fifth, we theoretically and experimentally compare the Euclidean metric and Riemannian metric for tensors, ODFs and EAPs. -Finally, we propose the Geodesic Anisotropy (GA) to measure the anisotropy of EAPs, Square Root Parameterized Estimation (SRPE) for nonnegative definite ODF/EAP estimation, weighted Riemannian mean/median for ODF/EAP interpolation, smoothing, atlas estimation. The concept of \emph{reasonable mean value interpolation} is presented for interpolation of general PDF data.L'IRM de diffusion est a ce jour la seule technique a meme d'observer in vivo et de fac¸on non-invasive les structures fines de la mati'ere blanche, en modelisant la diffusion des molecules d'eau. Le propagateur moyen (EAP pour Ensemble average Propagator en anglais) et la fonction de distribution d'orientation (ODF pour Orientation Distribution Function en anglais) sont les deux fonctions de probabilites d'int'erˆet pour caracteriser la diffusion des molecules d'eau. Le probleme central en IRM de diffusion est la reconstruction et le traitement de ces fonctions (EAP et ODF); c'est aussi le point de depart pour la tractographie des fibres de la mati'ere blanche. Le formalisme du tenseur de diffusion (DTI pour Diffusion Tensor Imaging en anglais) est le modele le plus couramment utilise, et se base sur une hypothese de diffusion gaussienne. Il existe un cadre riemannien qui permet d'estimer et de traiter correctement les images de tenseur de diffusion. Cependant, l'hypothese d'une diffusion gaussienne est une simplification, qui ne permet pas de d'écrire les cas ou la structure microscopique sous-jacente est complexe, tels que les croisements de faisceaux de fibres. L'imagerie 'a haute resolution angulaire (HARDI pour High Angular Resolution Diffusion Imaging en anglais) est un ensemble de methodes qui permettent de contourner les limites du modele tensoriel. La plupart des m'ethodes HARDI 'a ce jour, telles que l'imagerie spherique de l'espace de Fourier (QBI pour Q-Ball Imaging en anglais) se basent sur des hypoth'eses reductrices, et prennent en compte des acquisitions qui ne se font que sur une seule sphere dans l'espace de Fourier (sHARDI pour single-shell HARDI en anglais), c'est-a-dire une seule valeur du coefficient de ponderation b. Cependant, avec le developpement des scanners IRM et des techniques d'acquisition, il devient plus facile d'acquerir des donn'ees sur plusieurs sph'eres concentriques. Cette th'ese porte sur les methodes d'estimation et de traitement de donnees sur plusieurs spheres (mHARDI pour multiple-shell HARDI en anglais), et de facon generale sur les methodes de reconstruction independantes du schema d'echantillonnage. Cette these presente plusieurs contributions originales. En premier lieu, nous developpons l'imagerie par transformee de Fourier en coordonnees spheriques (SPFI pour Spherical Polar Fourier Imaging en anglais), qui se base sur une representation du signal dans une base de fonctions a parties radiale et angulaire separables (SPF basis pour Spherical Polar Fourier en anglais). Nous obtenons, de fac¸on analytique et par transformations lineaires, l'EAP ainsi que ses caracteristiques importantes : l'ODF, et des indices scalaires tels que l'anisotropie fractionnelle generalisee (GFA pour Generalized Fractional Anisotropy en anglais). En ce qui concerne l'implementation de SPFI, nous presentons deux methodes pour determiner le facteur d'echelle, et nous prenons en compte le fait que E(0) = 1 dans l'estimation. En second lieu, nous presentons un nouveau cadre pour une transformee de Fourier analytique en coordonnees spheriques (AFT-SC pour Analytical Fourier Transform in Spherical Coordinate en anglais), ce qui permet de considerer aussi bien les methodes mHARDI que sHARDI, d'explorer les relations entre ces methodes, et de developper de nouvelles techniques d'estimation de l'EAP et de l'ODF. Nous presentons en troisieme lieu d'importants crit'eres de comparaison des differentes methodes HARDI, ce qui permet de mettre en lumiere leurs avantages et leurs limites. Dans une quatrieme partie, nous proposons un nouveau cadre riemannien invariant par diffeomorphisme pour le traitement de l'EAP et de l'ODF. Ce cadre est une generalisation de la m'ethode riemannienne precedemment appliquee au tenseur de diffusion. Il peut etre utilise pour l'estimation d'une fonction de probabilite representee par sa racine carree, appelee fonction d'onde, dans une base de fonctions orthonormale. Dans ce cadre riemannien, les applications exponentielle et logarithmique, ainsi que les geodesiques ont une forme analytique. La moyenne riemannienne ponderee ainsi que la mediane existent et sont uniques, et peuvent etre calculees de facon efficace par descente de gradient. Nous developpons egalement un cadre log-euclidien et un cadre affine-euclidien pour un traitement rapide des donnees. En cinquieme partie, nous comparons, theoriquement et sur un plan exp'erimental, les metriques euclidiennes et riemanniennes pour les tenseurs, l'ODF et l'EAP. Finalement, nous proposons l'anisotropie geodesique (GA pour Geodesic Anisotropy en anglais) pour mesurer l'anisotropie de l'EAP; une parametrisation par la racine carrée (SRPE pour Square-Root Parameterized Estimation en anglais) pour l'estimation d'un EAP et d'une ODF positifs; la mediane et la moyenne riemanniennes ponderees pour l'interpolation, le lissage et la construction d'atlas bas'es sur l'ODF et de l'EAP. Nous introduisons la notion de valeur moyenne raisonnable pour l'interpolation de fonction de probabilites en general

Spatially Regularizing High Angular Resolution Diffusion Imaging

Many recent high angular resolution diffusion imaging (HARDI) reconstruction techniques have been introduced to infer ensemble average propagator (EAP),describing the three-dimensional (3D) average diffusion process of water molecules or the angular structure information contained in EAP, orientation distribution function (ODF). Most of these methods perform reconstruction independently at each voxel, which essentially ignoring the functional nature of the HARDI data at different voxels in space. The aim of my thesis is to develop methods which can spatially and adaptively infer the EAP, or ODF of water diffusion in regions with complex fiber configurations. In Chapter 3, we propose a penalized multi-scale adaptive regression model (PMARM) framework to spatially and adaptively infer the ODF of water diffusion in regions with complex fiber configurations. We first represent DW-MRI signals using Spherical Harmonic (SH) basis, then apply PMARM on advanced statistical methods to calculate the coefficients of SH representation, from which ODF representation is calculated using Funk-Radon transformation. PMARM reconstructs the ODF at each voxel by adaptively borrowing the spatial information from the neighboring voxels. We show in the real and simulated data sets that PMARM can substantially reduce the noise level, while improving the ODF reconstruction. In Chapter 4, we propose a robust multi-scale adaptive and sequential smoothing (MASS) method framework to robustly, spatially and adaptively infer the EAP of water diffusion in regions with complex fiber configurations. We first calculate spherical polar Fourier basis representation of the DW-MRI signals, and then apply MASS adaptively and sequentially updating SPF representation by borrowing the spatial information from the neighboring voxels. We show in the real and simulated data sets that MASS can reduce the angle detection errors on fiber crossing area and provides more accurate reconstructions than standard voxel-wise methods and robust MASS performs very well with the presence of outliers. In Chapter 5, we extend multi-scale adaptive method framework to dictionary learning methods, and show that by adding smoothing technique, we can significantly improve the accuracy of EAP reconstruction and reduce the angle detection errors on fiber crossing, even in very low signal-to-noise ratio situation.Doctor of Philosoph

Modélisation locale en imagerie par résonance magnétique de diffusion : de l'acquisition comprimée au connectome

L’imagerie par résonance magnétique pondérée en diffusion est une modalité d’imagerie médicale non invasive qui permet de mesurer les déplacements microscopiques des molécules d’eau dans les tissus biologiques. Il est possible d’utiliser cette information pour inférer la structure du cerveau. Les techniques de modélisation locale de la diffusion permettent de calculer l’orientation et la géométrie des tissus de la matière blanche. Cette thèse s’intéresse à l’optimisation des métaparamètres utilisés par les modèles locaux. Nous dérivons des paramètres optimaux qui améliorent la qualité des métriques de diffusion locale, de la tractographie de la matière blanche et de la connectivité globale. L’échantillonnage de l’espace-q est un des paramètres principaux qui limitent les types de modèle et d’inférence applicable sur des données acquises en clinique. Dans cette thèse, nous développons une technique d’échantillonnage de l’espace-q permettant d’utiliser l’acquisition comprimée pour réduire le temps d’acquisition nécessaire

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: algorithms and result

Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computational approaches have been proposed to harmonize diffusion MRI data or remove scanner-specific differences. To date, these methods have mostly been developed for or evaluated on single b-value diffusion MRI data. In this work, we present the evaluation results of 19 algorithms that are developed to harmonize the cross-scanner and cross-protocol variability of multi-shell diffusion MRI using a benchmark database. The proposed algorithms rely on various signal representation approaches and computational tools, such as rotational invariant spherical harmonics, deep neural networks and hybrid biophysical and statistical approaches. The benchmark database consists of data acquired from the same subjects on two scanners with different maximum gradient strength (80 and 300 ​mT/m) and with two protocols. We evaluated the performance of these algorithms for mapping multi-shell diffusion MRI data across scanners and across protocols using several state-of-the-art imaging measures. The results show that data harmonization algorithms can reduce the cross-scanner and cross-protocol variabilities to a similar level as scan-rescan variability using the same scanner and protocol. In particular, the LinearRISH algorithm based on adaptive linear mapping of rotational invariant spherical harmonics features yields the lowest variability for our data in predicting the fractional anisotropy (FA), mean diffusivity (MD), mean kurtosis (MK) and the rotationally invariant spherical harmonic (RISH) features. But other algorithms, such as DIAMOND, SHResNet, DIQT, CMResNet show further improvement in harmonizing the return-to-origin probability (RTOP). The performance of different approaches provides useful guidelines on data harmonization in future multi-site studies