Proyecto visualizador de imágenes cerebrales a través de tractografías

En este trabajo se presenta un proyecto de investigación con el fin de ayudar a estandarizar el problema de la construcción de tractografías a partir de las imágenes ponderada por difusión (DWI) el cual normalmente es un problema al que se enfrentan los neurólogos y neurocirujanos para identificar tractos cerebrales. Con el fin de solucionar el problema, la investigación se plantea desde la identificación y comparación de las herramientas y algoritmos utilizados hasta la elaboración de una herramienta web de prueba la cual permite la reconstrucción completa desde la imagen DWI hasta la tractografía

Imagerie de diffusion en temps-réel (correction du bruit et inférence de la connectivité cérébrale)

La plupart des constructeurs de systèmes d'imagerie par résonance magnétique (IRM) proposent un large choix d'applications de post-traitement sur les données IRM reconstruites a posteriori, mais très peu de ces applications peuvent être exécutées en temps réel pendant l'examen. Mises à part certaines solutions dédiées à l'IRM fonctionnelle permettant des expériences relativement simples ainsi que d'autres solutions pour l'IRM interventionnelle produisant des scans anatomiques pendant un acte de chirurgie, aucun outil n'a été développé pour l'IRM pondérée en diffusion (IRMd). Cependant, comme les examens d'IRMd sont extrêmement sensibles à des perturbations du système hardware ou à des perturbations provoquées par le sujet et qui induisent des données corrompues, il peut être intéressant d'investiguer la possibilité de reconstruire les données d'IRMd directement lors de l'examen. Cette thèse est dédiée à ce projet innovant. La contribution majeure de cette thèse a consisté en des solutions de débruitage des données d'IRMd en temps réel. En effet, le signal pondéré en diffusion peut être corrompu par un niveau élevé de bruit qui n'est plus gaussien, mais ricien ou chi non centré. Après avoir réalisé un état de l'art détaillé de la littérature sur le bruit en IRM, nous avons étendu l'estimateur linéaire qui minimise l'erreur quadratique moyenne (LMMSE) et nous l'avons adapté à notre cadre de temps réel réalisé avec un filtre de Kalman. Nous avons comparé les performances de cette solution à celles d'un filtrage gaussien standard, difficile à implémenter car il nécessite une modification de la chaîne de reconstruction pour y être inséré immédiatement après la démodulation du signal acquis dans l'espace de Fourier. Nous avons aussi développé un filtre de Kalman parallèle qui permet d'appréhender toute distribution de bruit et nous avons montré que ses performances étaient comparables à celles de notre méthode précédente utilisant un filtre de Kalman non parallèle. Enfin, nous avons investigué la faisabilité de réaliser une tractographie en temps-réel pour déterminer la connectivité structurelle en direct, pendant l'examen. Nous espérons que ce panel de développements méthodologiques permettra d'améliorer et d'accélérer le diagnostic en cas d'urgence pour vérifier l'état des faisceaux de fibres de la substance blanche.Most magnetic resonance imaging (MRI) system manufacturers propose a huge set of software applications to post-process the reconstructed MRI data a posteriori, but few of them can run in real-time during the ongoing scan. To our knowledge, apart from solutions dedicated to functional MRI allowing relatively simple experiments or for interventional MRI to perform anatomical scans during surgery, no tool has been developed in the field of diffusion-weighted MRI (dMRI). However, because dMRI scans are extremely sensitive to lots of hardware or subject-based perturbations inducing corrupted data, it can be interesting to investigate the possibility of processing dMRI data directly during the ongoing scan and this thesis is dedicated to this challenging topic. The major contribution of this thesis aimed at providing solutions to denoise dMRI data in real-time. Indeed, the diffusion-weighted signal may be corrupted by a significant level of noise which is not Gaussian anymore, but Rician or noncentral chi. After making a detailed review of the literature, we extended the linear minimum mean square error (LMMSE) estimator and adapted it to our real-time framework with a Kalman filter. We compared its efficiency to the standard Gaussian filtering, difficult to implement, as it requires a modification of the reconstruction pipeline to insert the filter immediately after the demodulation of the acquired signal in the Fourier space. We also developed a parallel Kalman filter to deal with any noise distribution and we showed that its efficiency was quite comparable to the non parallel Kalman filter approach. Last, we addressed the feasibility of performing tractography in real-time in order to infer the structural connectivity online. We hope that this set of methodological developments will help improving and accelerating a diagnosis in case of emergency to check the integrity of white matter fiber bundles.PARIS11-SCD-Bib. électronique (914719901) / SudocSudocFranceF

Quantitative evaluation of 10 tractography algorithms on a realistic diffusion MR phantom.

International audienceAs it provides the only method for mapping white matter fibers in vivo, diffusion MRI tractography is gaining importance in clinical and neuroscience research. However, despite the increasing availability of different diffusion models and tractography algorithms, it remains unclear how to select the optimal fiber reconstruction method, given certain imaging parameters. Consequently, it is of utmost importance to have a quantitative comparison of these models and algorithms and a deeper understanding of the corresponding strengths and weaknesses. In this work, we use a common dataset with known ground truth and a reproducible methodology to quantitatively evaluate the performance of various diffusion models and tractography algorithms. To examine a wide range of methods, the dataset, but not the ground truth, was released to the public for evaluation in a contest, the "Fiber Cup". 10 fiber reconstruction methods were evaluated. The results provide evidence that: 1. For high SNR datasets, diffusion models such as (fiber) orientation distribution functions correctly model the underlying fiber distribution and can be used in conjunction with streamline tractography, and 2. For medium or low SNR datasets, a prior on the spatial smoothness of either the diffusion model or the fibers is recommended for correct modelling of the fiber distribution and proper tractography results. The phantom dataset, the ground truth fibers, the evaluation methodology and the results obtained so far will remain publicly available on: http://www.lnao.fr/spip.php?rubrique79 to serve as a comparison basis for existing or new tractography methods. New results can be submitted to [email protected] and updates will be published on the webpage

Higher-Order Tensors and Differential Topology in Diffusion MRI Modeling and Visualization

Diffusion Weighted Magnetic Resonance Imaging (DW-MRI) is a noninvasive method for creating three-dimensional scans of the human brain. It originated mostly in the 1970s and started its use in clinical applications in the 1980s. Due to its low risk and relatively high image quality it proved to be an indispensable tool for studying medical conditions as well as for general scientific research. For example, it allows to map fiber bundles, the major neuronal pathways through the brain. But all evaluation of scanned data depends on mathematical signal models that describe the raw signal output and map it to biologically more meaningful values. And here we find the most potential for improvement. In this thesis we first present a new multi-tensor kurtosis signal model for DW-MRI. That means it can detect multiple overlapping fiber bundles and map them to a set of tensors. Compared to other already widely used multi-tensor models, we also add higher order kurtosis terms to each fiber. This gives a more detailed quantification of fibers. These additional values can also be estimated by the Diffusion Kurtosis Imaging (DKI) method, but we show that these values are drastically affected by fiber crossings in DKI, whereas our model handles them as intrinsic properties of fiber bundles. This reduces the effects of fiber crossings and allows a more direct examination of fibers. Next, we take a closer look at spherical deconvolution. It can be seen as a generalization of multi-fiber signal models to a continuous distribution of fiber directions. To this approach we introduce a novel mathematical constraint. We show, that state-of-the-art methods for estimating the fiber distribution become more robust and gain accuracy when enforcing our constraint. Additionally, in the context of our own deconvolution scheme, it is algebraically equivalent to enforcing that the signal can be decomposed into fibers. This means, tractography and other methods that depend on identifying a discrete set of fiber directions greatly benefit from our constraint. Our third major contribution to DW-MRI deals with macroscopic structures of fiber bundle geometry. In recent years the question emerged, whether or not, crossing bundles form two-dimensional surfaces inside the brain. Although not completely obvious, there is a mathematical obstacle coming from differential topology, that prevents general tangential planes spanned by fiber directions at each point to be connected into consistent surfaces. Research into how well this constraint is fulfilled in our brain is hindered by the high precision and complexity needed by previous evaluation methods. This is why we present a drastically simpler method that negates the need for precisely finding fiber directions and instead only depends on the simple diffusion tensor method (DTI). We then use our new method to explore and improve streamsurface visualization.<br /

Spatially Regularized Reconstruction of Fibre Orientation Distributions in the Presence of Isotropic Diffusion

The connectivity and structural integrity of the white matter of the brain is known to be implicated in a wide range of brain-related diseases and injuries. However, it is only since the advent of diffusion magnetic resonance imaging (dMRI) that researchers have been able to probe the miscrostructure of white matter in vivo. Presently, among a range of methods of dMRI, high angular resolution diffusion imaging (HARDI) is known to excel in its ability to provide reliable information about the local orientations of neural fasciculi (aka fibre tracts). It preserves the high angular resolution property of diffusion spectrum imaging (DSI) but requires less measurements. Meanwhile, as opposed to the more traditional diffusion tensor imaging (DTI), HARDI is capable of distinguishing the orientations of multiple fibres passing through a given spatial voxel. Unfortunately, the ability of HARDI to discriminate neural fibres that cross each other at acute angles is always limited. The limitation becomes the motivation to develop numerous post-processing tools, aiming at the improvement of the angular resolution of HARDI. Among such methods, spherical deconvolution (SD) is the one which attracts the most attentions. Due to its ill-posed nature, however, standard SD relies on a number of a priori assumptions needed to render its results unique and stable. In the present thesis, we introduce a novel approach to the problem of non-blind SD of HARDI signals, which does not only consider the existence of anisotropic diffusion component of HARDI signal but also explicitly take the isotropic diffusion component into account. As a result of that, in addition to reconstruction of fODFs, our algorithm can also yield a useful estimation of its related IDM, which quantifies a relative contribution of the isotropic diffusion component as well as its spatial pattern. Moreover, one of the principal contributions is to demonstrate the effectiveness of exploiting different prior models for regularization of the spatial-domain behaviours of the reconstructed fODFs and IDMs. Specifically, the fibre continuity model has been used to force the local maxima of the fODFs to vary consistently throughout the brain, whereas the bounded variation model has helped us to achieve piecewise smooth reconstruction of the IDMs. The proposed algorithm is formulated as a convex minimization problem, which admits a unique and stable minimizer. Moreover, using ADMM, we have been able to find the optimal solution via a sequence of simpler optimization problems, which are both computationally efficient and amenable to parallel computations. In a series of both in silico and in vivo experiments, we demonstrate how the proposed solution can be used to successfully overcome the effect of partial voluming, while preserving the spatial coherency of cerebral diffusion at moderate to severe noise levels. The performance of the proposed method is compared with that of several available alternatives, with the comparative results clearly supporting the viability and usefulness of our approach. Moreover, the results illustrate the power of applied spatial regularization terms