Incorporating Relaxivities to More Accurately Reconstruct MR Images

Purpose To develop a mathematical model that incorporates the magnetic resonance relaxivities into the image reconstruction process in a single step. Materials and methods In magnetic resonance imaging, the complex-valued measurements of the acquired signal at each point in frequency space are expressed as a Fourier transformation of the proton spin density weighted by Fourier encoding anomalies: T2⁎, T1, and a phase determined by magnetic field inhomogeneity (∆B) according to the MR signal equation. Such anomalies alter the expected symmetry and the signal strength of the k-space observations, resulting in images distorted by image warping, blurring, and loss in image intensity. Although T1 on tissue relaxation time provides valuable quantitative information on tissue characteristics, the T1 recovery term is typically neglected by assuming a long repetition time. In this study, the linear framework presented in the work of Rowe et al., 2007, and of Nencka et al., 2009 is extended to develop a Fourier reconstruction operation in terms of a real-valued isomorphism that incorporates the effects of T2⁎, ∆B, and T1. This framework provides a way to precisely quantify the statistical properties of the corrected image-space data by offering a linear relationship between the observed frequency space measurements and reconstructed corrected image-space measurements. The model is illustrated both on theoretical data generated by considering T2⁎, T1, and/or ∆B effects, and on experimentally acquired fMRI data by focusing on the incorporation of T1. A comparison is also made between the activation statistics computed from the reconstructed data with and without the incorporation of T1 effects. Result Accounting for T1 effects in image reconstruction is shown to recover image contrast that exists prior to T1 equilibrium. The incorporation of T1 is also shown to induce negligible correlation in reconstructed images and preserve functional activations. Conclusion With the use of the proposed method, the effects of T2⁎ and ∆B can be corrected, and T1 can be incorporated into the time series image-space data during image reconstruction in a single step. Incorporation of T1 provides improved tissue segmentation over the course of time series and therefore can improve the precision of motion correction and image registration

Fat fraction mapping using bSSFP Signal Profile Asymmetries for Robust multi-Compartment Quantification (SPARCQ)

Purpose: To develop a novel quantitative method for detection of different tissue compartments based on bSSFP signal profile asymmetries (SPARCQ) and to provide a validation and proof-of-concept for voxel-wise water-fat separation and fat fraction mapping. Methods: The SPARCQ framework uses phase-cycled bSSFP acquisitions to obtain bSSFP signal profiles. For each voxel, the profile is decomposed into a weighted sum of simulated profiles with specific off-resonance and relaxation time ratios. From the obtained set of weights, voxel-wise estimations of the fractions of the different components and their equilibrium magnetization are extracted. For the entire image volume, component-specific quantitative maps as well as banding-artifact-free images are generated. A SPARCQ proof-of-concept was provided for water-fat separation and fat fraction mapping. Noise robustness was assessed using simulations. A dedicated water-fat phantom was used to validate fat fractions estimated with SPARCQ against gold-standard 1H MRS. Quantitative maps were obtained in knees of six healthy volunteers, and SPARCQ repeatability was evaluated in scan rescan experiments. Results: Simulations showed that fat fraction estimations are accurate and robust for signal-to-noise ratios above 20. Phantom experiments showed good agreement between SPARCQ and gold-standard (GS) fat fractions (fF(SPARCQ) = 1.02*fF(GS) + 0.00235). In volunteers, quantitative maps and banding-artifact-free water-fat-separated images obtained with SPARCQ demonstrated the expected contrast between fatty and non-fatty tissues. The coefficient of repeatability of SPARCQ fat fraction was 0.0512. Conclusion: The SPARCQ framework was proposed as a novel quantitative mapping technique for detecting different tissue compartments, and its potential was demonstrated for quantitative water-fat separation.Comment: 20 pages, 7 figures, submitted to Magnetic Resonance in Medicin

A feasible and automatic free tool for T1 and ECV mapping

Purpose: Cardiac magnetic resonance (CMR) is a useful non-invasive tool for characterizing tissues and detecting myocardial fibrosis and edema. Estimation of extracellular volume fraction (ECV) using T1 sequences is emerging as an accurate biomarker in cardiac diseases associated with diffuse fibrosis. In this study, automatic software for T1 and ECV map generation consisting of an executable file was developed and validated using phantom and human data. Methods: T1 mapping was performed in phantoms and 30 subjects (22 patients and 8 healthy subjects) on a 1.5T MR scanner using the modified Look-Locker inversion-recovery (MOLLI) sequence prototype before and 15 min after contrast agent administration. T1 maps were generated using a Fast Nonlinear Least Squares algorithm. Myocardial ECV maps were generated using both pre- and post-contrast T1 image registration and automatic extraction of blood relaxation rates. Results: Using our software, pre- and post-contrast T1 maps were obtained in phantoms and healthy subjects resulting in a robust and reliable quantification as compared to reference software. Coregistration of pre- and post-contrast images improved the quality of ECV maps. Mean ECV value in healthy subjects was 24.5% ± 2.5%. Conclusions: This study demonstrated that it is possible to obtain accurate T1 maps and informative ECV maps using our software. Pixel-wise ECV maps obtained with this automatic software made it possible to visualize and evaluate the extent and severity of ECV alterations

HYDRA: Hybrid Deep Magnetic Resonance Fingerprinting

Purpose: Magnetic resonance fingerprinting (MRF) methods typically rely on dictio-nary matching to map the temporal MRF signals to quantitative tissue parameters. Such approaches suffer from inherent discretization errors, as well as high computational complexity as the dictionary size grows. To alleviate these issues, we propose a HYbrid Deep magnetic ResonAnce fingerprinting approach, referred to as HYDRA. Methods: HYDRA involves two stages: a model-based signature restoration phase and a learning-based parameter restoration phase. Signal restoration is implemented using low-rank based de-aliasing techniques while parameter restoration is performed using a deep nonlocal residual convolutional neural network. The designed network is trained on synthesized MRF data simulated with the Bloch equations and fast imaging with steady state precession (FISP) sequences. In test mode, it takes a temporal MRF signal as input and produces the corresponding tissue parameters. Results: We validated our approach on both synthetic data and anatomical data generated from a healthy subject. The results demonstrate that, in contrast to conventional dictionary-matching based MRF techniques, our approach significantly improves inference speed by eliminating the time-consuming dictionary matching operation, and alleviates discretization errors by outputting continuous-valued parameters. We further avoid the need to store a large dictionary, thus reducing memory requirements. Conclusions: Our approach demonstrates advantages in terms of inference speed, accuracy and storage requirements over competing MRF method

MRI overview for fat quantification in non-alcoholic fatty liver disease in the clinical and research settings

The general purpose of this master’s thesis is to describe the MRI techniques used in scanning and post processing for quantifying liver fat percentages for the purpose of diagnosis and research. At the onset we will look at epidemiological data regarding nonalcoholic fatty liver disease, which is often called by the name of hepatic steatosis. Based on the prevalence of this disease it is worthwhile to fully understand non-invasive (MRI) analysis, and its use in the clinical and research setting. Following an introductory section regarding the basis of magnetic resonance imaging, we will take a more in-depth look at current methods utilized for liver fat quantification. Due to the massive population of those of suffer from this disease worldwide it is prudent to analyze current methods, as well as the implications that such research has and will have on the pharmaceutical approach to treating this disease. The purpose of this thesis is to elucidate the MRI techniques utilized for liver fat quantification and provide a comprehensive view of how these techniques are used for diagnosis in the clinical setting, and longitudinal studies in the research setting to measure liver fat levels and how they react to various treatment approaches

Towards in vivo g-ratio mapping using MRI: unifying myelin and diffusion imaging

The g-ratio, quantifying the comparative thickness of the myelin sheath encasing an axon, is a geometrical invariant that has high functional relevance because of its importance in determining neuronal conduction velocity. Advances in MRI data acquisition and signal modelling have put in vivo mapping of the g-ratio, across the entire white matter, within our reach. This capacity would greatly increase our knowledge of the nervous system: how it functions, and how it is impacted by disease. This is the second review on the topic of g-ratio mapping using MRI. As such, it summarizes the most recent developments in the field, while also providing methodological background pertinent to aggregate g-ratio weighted mapping, and discussing pitfalls associated with these approaches. Using simulations based on recently published data, this review demonstrates the relevance of the calibration step for three myelin-markers (macromolecular tissue volume, myelin water fraction, and bound pool fraction). It highlights the need to estimate both the slope and offset of the relationship between these MRI-based markers and the true myelin volume fraction if we are really to achieve the goal of precise, high sensitivity g-ratio mapping in vivo. Other challenges discussed in this review further evidence the need for gold standard measurements of human brain tissue from ex vivo histology. We conclude that the quest to find the most appropriate MRI biomarkers to enable in vivo g-ratio mapping is ongoing, with the potential of many novel techniques yet to be investigated.Comment: Will be published as a review article in Journal of Neuroscience Methods as parf of the Special Issue with Hu Cheng and Vince Calhoun as Guest Editor

Characterization of nanoparticle transport through quartz and dolomite gravels by magnetic resonance imaging

Magnetic resonance imaging (MRI) has tremendous potential for revealing transport processes in engineered and geologic systems. Here, we utilize MRI to image nanoparticle (NP) transport through a saturated coarse-grained system. Commercially available paramagnetically tagged NPs are used; the paramagnetic tag making the NP visible to MRI. NP transport was imaged as NPs migrated through packed columns of quartz and dolomite gravel. Changes in T2-weighted image intensity were calibrated to provide fully quantitative maps of NP concentration at regular time intervals (T 2 being the spin–spin relaxation time of 1H nuclei). Transport of nanoparticles was significantly retarded in dolomite compared to quartz due to electrostatic attraction between nanoparticle and dolomite surfaces. NP concentration profiles were evaluated with the CXTFIT computer package to estimate nanoparticle transport parameters at multiple points along the length of the column. This provided temporally resolved parameters that standard breakthrough curve analysis cannot provide. Particle–surface interaction energy profiles were described through Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. While dispersion coefficients and fast deposition rate constant (k fast) were found to increase with distance, deposition rate constant (k) and collision efficiency (α) were found to decrease with distance. These length-dependant variations have significant scaling-up implications for transport models used to predict NP transport in natural and engineered coarse-grained systems, such as sustainable urban drainage systems and river beds