Deep learning-based parameter mapping for joint relaxation and diffusion tensor MR Fingerprinting

Magnetic Resonance Fingerprinting (MRF) enables the simultaneous quantification of multiple properties of biological tissues. It relies on a pseudo-random acquisition and the matching of acquired signal evolutions to a precomputed dictionary. However, the dictionary is not scalable to higher-parametric spaces, limiting MRF to the simultaneous mapping of only a small number of parameters (proton density, T1 and T2 in general). Inspired by diffusion-weighted SSFP imaging, we present a proof-of-concept of a novel MRF sequence with embedded diffusion-encoding gradients along all three axes to efficiently encode orientational diffusion and T1 and T2 relaxation. We take advantage of a convolutional neural network (CNN) to reconstruct multiple quantitative maps from this single, highly undersampled acquisition. We bypass expensive dictionary matching by learning the implicit physical relationships between the spatiotemporal MRF data and the T1, T2 and diffusion tensor parameters. The predicted parameter maps and the derived scalar diffusion metrics agree well with state-of-the-art reference protocols. Orientational diffusion information is captured as seen from the estimated primary diffusion directions. In addition to this, the joint acquisition and reconstruction framework proves capable of preserving tissue abnormalities in multiple sclerosis lesions

A model-based method with joint sparsity constant for direct diffusion tensor estimation

Diffusion tensor imaging (DTI) has been widely used for nondestructive characterization of microstructures of myocardium or brain connectivity. It requires repeated acquisition with different diffusion gradients. The long acquisition time greatly limits the clinical application of DTI. In this paper, a novel method, named model-based method with joint sparsity constraint (MB-JSC), effectively incorporates the prior information on the joint sparsity of different diffusion-weighted images in direct estimation of the diffusion tensor from highly undersampled k-space data. Experimental results demonstrate that the proposed method is able to estimate the diffusion tensors more accurately than the existing method when a high net reduction factor is used.published_or_final_versionThe 9th IEEE International Symposium on Biomedical Imaging (ISBI 2012), Barcelona, Spain, 2-5 May 2012. In Proceedings of the 9th ISBI, 2012, p. 510-51

Doctor of Philosophy

dissertationDiffusion tensor MRI (DT-MRI or DTI) has been proven useful for characterizing biological tissue microstructure, with the majority of DTI studies having been performed previously in the brain. Other studies have shown that changes in DTI parameters are detectable in the presence of cardiac pathology, recovery, and development, and provide insight into the microstructural mechanisms of these processes. However, the technical challenges of implementing cardiac DTI in vivo, including prohibitive scan times inherent to DTI and measuring small-scale diffusion in the beating heart, have limited its widespread usage. This research aims to address these technical challenges by: (1) formulating a model-based reconstruction algorithm to accurately estimate DTI parameters directly from fewer MRI measurements and (2) designing novel diffusion encoding MRI pulse sequences that compensate for the higher-order motion of the beating heart. The model-based reconstruction method was tested on undersampled DTI data and its performance was compared against other state-of-the-art reconstruction algorithms. Model-based reconstruction was shown to produce DTI parameter maps with less blurring and noise and to estimate global DTI parameters more accurately than alternative methods. Through numerical simulations and experimental demonstrations in live rats, higher-order motion compensated diffusion-encoding was shown to successfully eliminate signal loss due to motion, which in turn produced data of sufficient quality to accurately estimate DTI parameters, such as fiber helix angle. Ultimately, the model-based reconstruction and higher-order motion compensation methods were combined to characterize changes in the cardiac microstructure in a rat model with inducible arterial hypertension in order to demonstrate the ability of cardiac DTI to detect pathological changes in living myocardium

Simultaneous Multislice (SMS) Imaging Techniques

Item does not contain fulltex

Denoising and fast diffusion imaging with physically constrained sparse dictionary learning

International audienceDiffusion-weighted imaging (DWI) allows imaging the geometry of water diffusion in biological tissues. However, DW images are noisy at high b-values and acquisitions are slow when using a large number of measurements, such as in Diffusion Spectrum Imaging (DSI). This work aims to denoise DWI and reduce the number of required measurements, while maintaining data quality. To capture the structure of DWI data, we use sparse dictionary learning constrained by the physical properties of the signal: symmetry and positivity. The method learns a dictionary of diffusion profiles on all the DW images at the same time and then scales to full brain data. Its performance is investigated with simulations and two real DSI datasets. We obtain better signal estimates from noisy measurements than by applying mirror symmetry through the q-space origin, Gaussian denoising or state-of- the-art non-local means denoising. Using a high-resolution dictionary learnt on another subject, we show that we can reduce the number of images acquired while still generating high resolution DSI data. Using dictionary learning, one can denoise DW images effectively and perform faster acquisitions. Higher b-value acquisitions and DSI techniques are possible with approximately 40 measurements. This opens important perspectives for the connectomics community using DSI

Spherical deconvolution of multichannel diffusion MRI data with non-Gaussian noise models and spatial regularization

Spherical deconvolution (SD) methods are widely used to estimate the intra-voxel white-matter fiber orientations from diffusion MRI data. However, while some of these methods assume a zero-mean Gaussian distribution for the underlying noise, its real distribution is known to be non-Gaussian and to depend on the methodology used to combine multichannel signals. Indeed, the two prevailing methods for multichannel signal combination lead to Rician and noncentral Chi noise distributions. Here we develop a Robust and Unbiased Model-BAsed Spherical Deconvolution (RUMBA-SD) technique, intended to deal with realistic MRI noise, based on a Richardson-Lucy (RL) algorithm adapted to Rician and noncentral Chi likelihood models. To quantify the benefits of using proper noise models, RUMBA-SD was compared with dRL-SD, a well-established method based on the RL algorithm for Gaussian noise. Another aim of the study was to quantify the impact of including a total variation (TV) spatial regularization term in the estimation framework. To do this, we developed TV spatially-regularized versions of both RUMBA-SD and dRL-SD algorithms. The evaluation was performed by comparing various quality metrics on 132 three-dimensional synthetic phantoms involving different inter-fiber angles and volume fractions, which were contaminated with noise mimicking patterns generated by data processing in multichannel scanners. The results demonstrate that the inclusion of proper likelihood models leads to an increased ability to resolve fiber crossings with smaller inter-fiber angles and to better detect non-dominant fibers. The inclusion of TV regularization dramatically improved the resolution power of both techniques. The above findings were also verified in brain data

Fast diffusion MRI based on sparse acquisition and reconstruction for long-term population imaging

Diffusion weighted magnetic resonance imaging (dMRI) is a unique MRI modality to probe the diffusive molecular transport in biological tissue. Due to its noninvasiveness and its ability to investigate the living human brain at submillimeter scale, dMRI is frequently performed in clinical and biomedical research to study the brain’s complex microstructural architecture. Over the last decades large prospective cohort studies have been set up with the aim to gain new insights into the development and progression of brain diseases across the life span and to discover biomarkers for disease prediction and potentially prevention. To allow for diverse brain imaging using different MRI modalities, stringent scan time limits are typically imposed in population imaging. Nevertheless, population studies aim to apply advanced and thereby time consuming dMRI protocols that deliver high quality data with great potential for future analysis. To allow for time-efficient but also versatile diffusion imaging, this thesis contributes to the investigation of accelerating diffusion spectrum imaging (DSI), an advanced dMRI technique that acquires imaging data with high intra-voxel resolution of tissue microstructure. Combining state-of-the-art parallel imaging and the theory of compressed sensing (CS) enables the acceleration of spatial encoding and diffusion encoding in dMRI. In this way, the otherwise long acquisition times in DSI can be reduced significantly. In this thesis, first, suitable q-space sampling strategies and basis functions are explored that fulfill the requirements of CS theory for accurate sparse DSI reconstruction. Novel 3D q-space sample distributions are investigated for CS-DSI. Moreover, conventional CS-DSI based on the discrete Fourier transform is compared for the first time to CS-DSI based on the continuous SHORE (simple harmonic oscillator based reconstruction and estimation) basis functions. Based on these findings, a CS-DSI protocol is proposed for application in a prospective cohort study, the Rhineland Study. A pilot study was designed and conducted to evaluate the CS-DSI protocol in comparison with state-of-the-art 3-shell dMRI and dedicated protocols for diffusion tensor imaging (DTI) and for the combined hindered and restricted model of diffusion (CHARMED). Population imaging requires processing techniques preferably with low computational cost to process and analyze the acquired big data within a reasonable time frame. Therefore, a pipeline for automated processing of CS-DSI acquisitions was implemented including both in-house developed and existing state-of-the-art processing tools. The last contribution of this thesis is a novel method for automatic detection and imputation of signal dropout due to fast bulk motion during the diffusion encoding in dMRI. Subject motion is a common source of artifacts, especially when conducting clinical or population studies with children, the elderly or patients. Related artifacts degrade image quality and adversely affect data analysis. It is, thus, highly desired to detect and then exclude or potentially impute defective measurements prior to dMRI analysis. Our proposed method applies dMRI signal modeling in the SHORE basis and determines outliers based on the weighted model residuals. Signal imputation reconstructs corrupted and therefore discarded measurements from the sparse set of inliers. This approach allows for fast and robust correction of imaging artifacts in dMRI which is essential to estimate accurate and precise model parameters that reflect the diffusive transport of water molecules and the underlying microstructural environment in brain tissue.Die diffusionsgewichtete Magnetresonanztomographie (dMRT) ist ein einzigartiges MRTBildgebungsverfahren, um die Diffusionsbewegung von Wassermolekülen in biologischem Gewebe zu messen. Aufgrund der Möglichkeit Schichtbilder nicht invasiv aufzunehmen und das lebende menschliche Gehirn im Submillimeter-Bereich zu untersuchen, ist die dMRT ein häufig verwendetes Bildgebungsverfahren in klinischen und biomedizinischen Studien zur Erforschung der komplexen mikrostrukturellen Architektur des Gehirns. In den letzten Jahrzehnten wurden große prospektive Kohortenstudien angelegt, um neue Einblicke in die Entwicklung und den Verlauf von Gehirnkrankheiten über die Lebenspanne zu erhalten und um Biomarker zur Krankheitserkennung und -vorbeugung zu bestimmen. Um durch die Verwendung unterschiedlicher MRT-Verfahren verschiedenartige Schichtbildaufnahmen des Gehirns zu ermöglich, müssen Scanzeiten typischerweise stark begrenzt werden. Dennoch streben Populationsstudien die Anwendung von fortschrittlichen und daher zeitintensiven dMRT-Protokollen an, um Bilddaten in hoher Qualität und mit großem Potential für zukünftige Analysen zu akquirieren. Um eine zeiteffizente und gleichzeitig vielseitige Diffusionsbildgebung zu ermöglichen, leistet diese Dissertation Beiträge zur Untersuchung von Beschleunigungsverfahren für die Bildgebung mittels diffusion spectrum imaging (DSI). DSI ist ein fortschrittliches dMRT-Verfahren, das Bilddaten mit hoher intra-voxel Auflösung der Gewebestruktur erhebt. Werden modernste Verfahren zur parallelen MRT-Bildgebung mit der compressed sensing (CS) Theorie kombiniert, ermöglicht dies eine Beschleunigung der räumliche Kodierung und der Diffusionskodierung in der dMRT. Dadurch können die ansonsten langen Aufnahmezeiten für DSI erheblich reduziert werden. In dieser Arbeit werden zuerst geeigenete Strategien zur Abtastung des q-space sowie Basisfunktionen untersucht, welche die Anforderungen der CS-Theorie für eine korrekte Signalrekonstruktion der dünnbesetzten DSI-Daten erfüllen. Neue 3D-Verteilungen von Messpunkten im q-space werden für die Verwendung in CS-DSI untersucht. Außerdem wird konventionell auf der diskreten Fourier-Transformation basierendes CS-DSI zum ersten Mal mit einem CS-DSI Verfahren verglichen, welches kontinuierliche SHORE (simple harmonic oscillator based reconstruction and estimation) Basisfunktionen verwendet. Aufbauend auf diesen Ergebnissen wird ein CS-DSI-Protokoll zur Anwendung in einer prospektiven Kohortenstudie, der Rheinland Studie, vorgestellt. Eine Pilotstudie wurde entworfen und durchgeführt, um das CS-DSI-Protokoll im Vergleich mit modernster 3-shell-dMRT und mit dedizierten Protokollen für diffusion tensor imaging (DTI) und für das combined hindered and restricted model of diffusion (CHARMED) zu evaluieren. Populationsbildgebung erfordert Prozessierungsverfahren mit möglichst geringem Rechenaufwand, um große akquirierte Datenmengen in einem angemessenen Zeitrahmen zu verarbeiten und zu analysieren. Dafür wurde eine Pipeline zur automatisierten Verarbeitung von CS-DSI-Daten implementiert, welche sowohl eigenentwickelte als auch bereits existierende moderene Verarbeitungsprogramme enthält. Der letzte Beitrag dieser Arbeit ist eine neue Methode zur automatischen Detektion und Imputation von Signalabfall, welcher durch schnelle Bewegungen während der Diffusionskodierung in der dMRT entsteht. Bewegungen der Probanden während der dMRT-Aufnahme sind eine häufige Ursache für Bildfehler, vor allem in klinischen oder Populationsstudien mit Kindern, alten Menschen oder Patienten. Diese Artefakte vermindern die Datenqualität und haben einen negativen Einfluss auf die Datenanalyse. Daher ist es das Ziel, fehlerhafte Messungen vor der dMRI-Analyse zu erkennen und dann auszuschließen oder wenn möglich zu ersetzen. Die vorgestellte Methode verwendet die SHORE-Basis zur dMRT-Signalmodellierung und bestimmt Ausreißer mit Hilfe von gewichteten Modellresidualen. Die Datenimputation rekonstruiert die unbrauchbaren und daher verworfenen Messungen mit Hilfe der verbleibenden, dünnbesetzten Menge an Messungen. Dieser Ansatz ermöglicht eine schnelle und robuste Korrektur von Bildartefakten in der dMRT, welche erforderlich ist, um korrekte und präzise Modellparameter zu schätzen, die die Diffusionsbewegung von Wassermolekülen und die zugrundeliegende Mikrostruktur des Gehirngewebes reflektieren

Study of Human Muscle Structure and Function with Velocity Encoded Phase Contrast and Diffusion Tensor Magnetic Resonance Imaging Techniques

The disproportionate loss of muscle force with aging and disuse atrophy compared to the loss of muscle mass is not yet completely understood. In addition to well-established neural and contractile determinants of force loss, remodeling of the extracellular matrix (ECM) has been recently shown in animal models to be another important contributor. In-vivo human studies exploring the structural remodeling of the ECM and its functional consequences are lacking due to the paucity of appropriate imaging techniques. This study focuses on the development and application of advanced Magnetic Resonance Imaging (MRI) methods to elucidate the mechanisms of loss of force with aging and disuse atrophy with the focus on ECM. Functional changes are investigated by strain and strain rate tensor mapping of muscle under different contraction paradigms using Velocity Encoded Phase-Contrast MRI. Methodological advances include improvements in hardware and software control of the dynamic studies. To overcome the limitation of long scan times, compressed sensing MR acquisition and reconstruction framework to reduce scan times to under a minute were developed. A multi-step automated analysis pipeline to extract 3D strain/strain rate tensors from the velocity images was implemented to process the large dynamic volumes. Strain indices reflecting the material properties of the ECM were shown to correlate with force loss leading to a hypothesis that shear strain may serve as a surrogate marker for lateral transmission of force. Diffusion tensor imaging has been applied previously to study skeletal muscle fiber architecture. The resolution of the images precludes direct inferences to be made about the microstructure. To address this limitation, bicompartmental and Random Permeable Barrier models of diffusion were applied to the diffusion data obtained with spin-echo and custom-developed stimulated echo echo-planar-imaging sequences respectively. Model derived parameters (fiber diameter, wall permeability) obtained from fitting time-dependent diffusion data were in physiologically reasonable range, with potential for tracking age related changes in muscle microstructure. The developed imaging and modeling techniques were applied to a cohort of young/senior subjects and to longitudinal tracking of disuse atrophy induced by Unilateral Limb Suspension. These studies may potentially provide the causal link between age- and disuse-related structural remodeling and its functional consequences

Hybrid-space reconstruction with add-on distortion correction for simultaneous multi-slab diffusion MRI

Purpose: This study aims to propose a model-based reconstruction algorithm for simultaneous multi-slab diffusion MRI acquired with blipped-CAIPI gradients (blipped-SMSlab), which can also incorporate distortion correction. Methods: We formulate blipped-SMSlab in a 4D k-space with kz gradients for the intra-slab slice encoding and km (blipped-CAIPI) gradients for the inter-slab encoding. Because kz and km gradients share the same physical axis, the blipped-CAIPI gradients introduce phase interference in the z-km domain while motion induces phase variations in the kz-m domain. Thus, our previous k-space-based reconstruction would need multiple steps to transform data back and forth between k-space and image space for phase correction. Here we propose a model-based hybrid-space reconstruction algorithm to correct the phase errors simultaneously. Moreover, the proposed algorithm is combined with distortion correction, and jointly reconstructs data acquired with the blip-up/down acquisition to reduce the g-factor penalty. Results: The blipped-CAIPI-induced phase interference is corrected by the hybrid-space reconstruction. Blipped-CAIPI can reduce the g-factor penalty compared to the non-blipped acquisition in the basic reconstruction. Additionally, the joint reconstruction simultaneously corrects the image distortions and improves the 1/g-factors by around 50%. Furthermore, through the joint reconstruction, SMSlab acquisitions without the blipped-CAIPI gradients also show comparable correction performance with blipped-SMSlab. Conclusion: The proposed model-based hybrid-space reconstruction can reconstruct blipped-SMSlab diffusion MRI successfully. Its extension to a joint reconstruction of the blip-up/down acquisition can correct EPI distortions and further reduce the g-factor penalty compared with the separate reconstruction.Comment: 10 figures+tables, 8 supplementary figure

Development of Methodologies for Diffusion-weighted Magnetic Resonance Imaging at High Field Strength

Diffusion-weighted imaging of small animals at high field strengths is a challenging prospect due to its extreme sensitivity to motion. Periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) was introduced at 9.4T as an imaging method that is robust to motion and distortion. Proton density (PD)-weighted and T2-weighted PROPELLER data were generally superior to that acquired with single-shot, Cartesian and echo planar imaging-based methods in terms of signal-to-noise ratio (SNR), contrast-to-noise ratio and resistance to artifacts. Simulations and experiments revealed that PROPELLER image quality was dependent on the field strength and echo times specified. In particular, PD-weighted imaging at high field led to artifacts that reduced image contrast. In PROPELLER, data are acquired in progressively rotated blades in k-space and combined on a Cartesian grid. PROPELLER with echo truncation at low spatial frequencies (PETALS) was conceived as a postprocessing method that improved contrast by reducing the overlap of k-space data from different blades with different echo times. Where the addition of diffusion weighting gradients typically leads to catastrophic motion artifacts in multi-shot sequences, diffusion-weighted PROPELLER enabled the acquisition of high quality, motion-robust data. Applications in the healthy mouse brain and abdomen at 9.4T and in stroke patients at 3T are presented. PROPELLER increases the minimum scan time by approximately 50%. Consequently, methods were explored to reduce the acquisition time. Two k-space undersampling regimes were investigated by examining image fidelity as a function of degree of undersampling. Undersampling by acquiring fewer k-space blades was shown to be more robust to motion and artifacts than undersampling by expanding the distance between successive phase encoding steps. To improve the consistency of undersampled data, the non-uniform fast Fourier transform was employed. It was found that acceleration factors of up to two could be used with minimal visual impact on image fidelity. To reduce the number of scans required for isotropic diffusion weighting, the use of rotating diffusion gradients was investigated, exploiting the rotational symmetry of the PROPELLER acquisition. Fixing the diffusion weighting direction to the individual rotating blades yielded geometry and anisotropy-dependent diffusion measurements. However, alternating the orientations of diffusion weighting with successive blades led to more accurate measurements of the apparent diffusion coefficient while halving the overall acquisition time. Optimized strategies are proposed for the use of PROPELLER in rapid high resolution imaging at high field strength