Search

Conventional magnetic resonance imaging (MRI) is hampered by long scan times and only qualitative image contrasts that prohibit a direct comparison between different systems. To address these limitations, model-based reconstructions explicitly model the physical laws that govern the MRI signal generation. By formulating image reconstruction as an inverse problem, quantitative maps of the underlying physical parameters can then be extracted directly from efficiently acquired k-space signals without intermediate image reconstruction-addressing both shortcomings of conventional MRI at the same time. This review will discuss basic concepts of model-based reconstructions and report on our experience in developing several model-based methods over the last decade using selected examples that are provided complete with data and code. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.

A novel method is developed that adaptively generates tiny dictionaries for joint T1-T2 mapping in magnetic resonance imaging. This work breaks the bond between dictionary size and representation accuracy (i) by approximating the Bloch-response manifold by piece-wise linear functions and (ii) by adaptively refining the sampling grid depending on the locally-linear approximation error. Data acquisition is accomplished with use of an 2D radially sampled Inversion-Recovery Hybrid-State Free Precession sequence. Adaptive dictionaries are generated with different error tolerances and compared to a heuristically designed dictionary. Based on simulation results, tiny dictionaries were used for T1-T2 mapping in phantom and in vivo studies. Reconstruction and parameter mapping were performed entirely in subspace. All experiments demonstrated excellent agreement between the proposed mapping technique and template matching using heuristic dictionaries. Adaptive dictionaries in combination with manifold projection allow to reduce the necessary dictionary sizes by one to two orders of magnitude.

Purpose: A novel subspace-based reconstruction method for frequency-modulated balanced steady-state free precession (fmSSFP) MRI is presented. In this work, suitable data acquisition schemes, subspace sizes, and efficiencies for banding removal are investigated., Theory and Methods: By combining a fmSSFP MRI sequence with a 3D stack-of-stars trajectory, scan efficiency is maximized as spectral information is obtained without intermediate preparation phases. A memory-efficient reconstruction routine is implemented by introducing the low-frequency Fourier transform as a subspace which allows for the formulation of a convex reconstruction problem. The removal of banding artifacts is investigated by comparing the proposed acquisition and reconstruction technique to phase-cycled bSSFP MRI. Aliasing properties of different undersampling schemes are analyzed and water/fat separation is demonstrated by reweighting the reconstructed subspace coefficients to generate virtual spectral responses in a post-processing step., Results: A simple root-of-sum-of-squares combination of the reconstructed subspace coefficients yields high-SNR images with the characteristic bSSFP contrast but without banding artifacts. Compared to Golden-Angle trajectories, turn-based sampling schemes were superior in minimizing aliasing across reconstructed subspace coefficients. Water/fat separated images of the human knee were obtained by reweighting subspace coefficients., Conclusions: The novel subspace-based fmSSFP MRI technique emerges as a time-efficient alternative to phase-cycled bSFFP. The method does not need intermediate preparation phases, offers high SNR and avoids banding artifacts. Reweighting of the reconstructed subspace coefficients allows for generating virtual spectral responses with applications to water/fat separation., (© 2018 International Society for Magnetic Resonance in Medicine.)

Purpose: To develop a high-speed multislice T1 mapping method based on a single-shot inversion-recovery (IR) radial FLASH acquisition and a regularized model-based reconstruction., Methods: Multislice radial k-space data are continuously acquired after a single nonselective inversion pulse using a golden-angle sampling scheme in a spoke-interleaved manner with optimized flip angles. Parameter maps and coil sensitivities of each slice are estimated directly from highly undersampled radial k-space data using a model-based nonlinear inverse reconstruction in conjunction with joint sparsity constraints. The performance of the method has been validated using a numerical and experimental T1 phantom as well as demonstrated for studies of the human brain and liver at 3T., Results: The proposed method allows for 7 simultaneous T1 maps of the brain at 0.5 × 0.5 × 4 mm 3 resolution within a single IR experiment of 4 s duration. Phantom studies confirm similar accuracy and precision as obtained for a single-slice acquisition. For abdominal applications, the proposed method yields three simultaneous T1 maps at 1.25 × 1.25 × 6 mm 3 resolution within a 4 s breath hold., Conclusion: Rapid, robust, accurate, and precise multislice T1 mapping may be achieved by combining the advantages of a model-based nonlinear inverse reconstruction, radial sampling, parallel imaging, and compressed sensing.

Purpose: To develop a model-based reconstruction technique for single-shot T 1 mapping with high spatial resolution, accuracy, and precision using an inversion-recovery (IR) fast low-angle shot (FLASH) acquisition with radial encoding., Methods: The proposed model-based reconstruction jointly estimates all model parameters, that is, the equilibrium magnetization, steady-state magnetization, 1/ T1*, and all coil sensitivities from the data of a single-shot IR FLASH acquisition with a small golden-angle radial trajectory. Joint sparsity constraints on the parameter maps are exploited to improve the performance of the iteratively regularized Gauss-Newton method chosen for solving the nonlinear inverse problem. Validations include both a numerical and experimental T 1 phantom, as well as in vivo studies of the human brain and liver at 3 T., Results: In comparison to previous reconstruction methods for single-shot T 1 mapping, which are based on real-time MRI with pixel-wise fitting and a model-based approach with a predetermination of coil sensitivities, the proposed method presents with improved robustness against phase errors and numerical precision in both phantom and in vivo studies., Conclusion: The comprehensive model-based reconstruction with L1 regularization offers rapid and robust T 1 mapping with high accuracy and precision. The method warrants accelerated computing and online implementation for extended clinical trials. Magn Reson Med 79:730-740, 2018. © 2017 International Society for Magnetic Resonance in Medicine., (© 2017 International Society for Magnetic Resonance in Medicine.)

Purpose: To develop a model-based reconstruction technique for real-time phase-contrast flow MRI with improved spatiotemporal accuracy in comparison to methods using phase differences of two separately reconstructed images with differential flow encodings., Methods: The proposed method jointly computes a common image, a phase-contrast map, and a set of coil sensitivities from every pair of flow-compensated and flow-encoded datasets obtained by highly undersampled radial FLASH. Real-time acquisitions with five and seven radial spokes per image resulted in 25.6 and 35.7 ms measuring time per phase-contrast map, respectively. The signal model for phase-contrast flow MRI requires the solution of a nonlinear inverse problem, which is accomplished by an iteratively regularized Gauss-Newton method. Aspects of regularization and scaling are discussed. The model-based reconstruction was validated for a numerical and experimental flow phantom and applied to real-time phase-contrast MRI of the human aorta for 10 healthy subjects and 2 patients., Results: Under all conditions, and compared with a previously developed real-time flow MRI method, the proposed method yields quantitatively accurate phase-contrast maps (i.e., flow velocities) with improved spatial acuity, reduced phase noise and reduced streaking artifacts., Conclusion: This novel model-based reconstruction technique may become a new tool for clinical flow MRI in real time. Magn Reson Med 77:1082-1093, 2017. © 2016 International Society for Magnetic Resonance in Medicine., (© 2016 International Society for Magnetic Resonance in Medicine.)

Purpose: Together with the development of MRI contrasts that are inherently small in their magnitude, increased magnetic field accuracy is also required. Hence, mapping of the static magnetic field (B 0 ) and the excitation field (B 1 ) is not only important to feedback shim algorithms, but also for postprocess contrast-correction procedures., Methods: A novel field-inhomogeneity mapping method is presented that allows simultaneous mapping of the water shift and B 1 (WASABI) using an off-resonant rectangular preparation pulse. The induced Rabi oscillations lead to a sinc-like spectrum in the frequency-offset dimension and allow for determination of B 0 by its symmetry axis and of B 1 by its oscillation frequency., Results: Stability of the WASABI method with regard to the influences of T 1 , T 2 , magnetization transfer, and repetition time was investigated and its convergence interval was verified. B 0 and B 1 maps obtained simultaneously by means of WASABI in the human brain at 3 T and 7 T can compete well with maps obtained by standard methods. Finally, the method was applied successfully for B 0 and B 1 correction of chemical exchange saturation transfer MRI (CEST-MRI) data of the human brain., Conclusion: The proposed WASABI method yields a novel simultaneous B 0 and B 1 mapping within 1 min that is robust and easy to implement. Magn Reson Med 77:571-580, 2017. © 2016 International Society for Magnetic Resonance in Medicine., (© 2016 International Society for Magnetic Resonance in Medicine.)

Purpose: To develop a method for spoiling transverse magnetizations without additional gradients to minimize repetition times for radial fast low angle shot (FLASH) MRI., Methods: Residual steady state transverse magnetizations and corresponding image artifacts were analyzed for radial gradient echo sequences with constant and randomized RF phases in comparison with a sequence with refocused frequency-encoding gradients, constant spoiler gradient, and conventional RF spoiling (gold standard). The spoiling performance was assessed for different radial trajectories using numerical simulations, phantom experiments, and in vivo MRI studies of the human brain., Results: Simulations as well as phantom and in vivo measurements reveal a highly efficient spoiling capacity for randomized RF phases and radial FLASH sequences without the need for gradient rewinding and spoiler gradients. The data also demonstrate a strong dependence of the spoiling performance on the chosen radial trajectory (ie, the azimuthal angular increment between successive projections) with excellent results for an interleaved multiturn scheme., Conclusion: Effective spoiling of transverse magnetizations in radial FLASH MRI may be achieved by randomized RF phases without additional spoiler gradients. The technique allows for short repetition times as required for high-speed real-time MRI., (© 2015 Wiley Periodicals, Inc.)

Purpose: To provide multidimensional velocity compensation for real-time phase-contrast flow MRI., Methods: The proposed method introduces asymmetric gradient echoes for highly undersampled radial FLASH MRI with phase-sensitive image reconstruction by regularized nonlinear inversion (NLINV). Using an adapted gradient delay correction the resulting image quality was analyzed by simulations and experimentally validated at 3 Tesla. For real-time flow MRI the reduced gradient-echo timing allowed for the incorporation of velocity-compensating waveforms for all imaging gradients at even shorter repetition times., Results: The results reveal a usable degree of 20% asymmetry. Real-time flow MRI with full velocity compensation eliminated signal void in a flow phantom, confirmed flow parameters in healthy subjects and demonstrated signal recovery and phase conservation in a patient with aortic valve insufficiency and stenosis. Exemplary protocols at 1.4-1.5 mm resolution and 6 mm slice thickness achieved total acquisition times of 33.3-35.7 ms for two images (7 spokes each) with and without flow-encoding gradient., Conclusion: Asymmetric gradient echoes were successfully implemented for highly undersampled radial trajectories. The resulting temporal gain offers full velocity compensation for real-time phase-contrast flow MRI which minimizes false-positive contributions from complex flow and further enhances the temporal resolution compared with acquisitions with symmetric echoes., (© 2015 Wiley Periodicals, Inc.)

Amide protons of peptide bonds induce an important chemical exchange saturation transfer (CEST) contrast in vivo. As a simple in vitro model for a peptide amide proton CEST effect, we suggest herein the dipeptide carnosine. We show that the metabolite carnosine creates a CEST effect and we study the properties of the exchange of the amide proton (-NH) of the carnosine peptide bond (NHCPB) in model solutions for a pH range from 6 to 8.3 and a temperature range from T = 5 °C to 43 °C by means of CEST and water exchange spectroscopy (WEX) experiments on a 3 T whole-body MR tomograph. The dependence of the NHCPB chemical exchange rate k(sw) on pH and temperature T was determined using WEX. For physiological conditions (T = 37 °C, pH = 7.10) we obtained k(sw) = (47.07 ± 7.90)/s. With similar chemical shift and exchange properties to amide protons in vivo, carnosine forms a simple model system for optimization of CEST pulse sequences in vitro. The potential for direct detection of the metabolite carnosine in vivo is discussed., (Copyright © 2015 John Wiley & Sons, Ltd.)

Using diffusion-tensor magnetic resonance imaging and fiber tractography the topographic organization of the human corpus callosum (CC) has been described to comprise five segments with fibers projecting into prefrontal (I), premotor and supplementary motor (II), primary motor (III), and primary sensory areas (IV), as well as into parietal, temporal, and occipital cortical areas (V). In order to more rapidly characterize the underlying anatomy of these segments, this study used a novel single-shot T1 mapping method to quantitatively determine T1 relaxation times in the human CC. A region-of-interest analysis revealed a tendency for the lowest T1 relaxation times in the genu and the highest T1 relaxation times in the somatomotor region of the CC. This observation separates regions dominated by myelinated fibers with large diameters (somatomotor area) from densely packed smaller axonal bundles (genu) with less myelin. The results indicate that characteristic T1 relaxation times in callosal profiles provide an additional means to monitor differences in fiber anatomy, fiber density, and gray matter in respective neocortical areas. In conclusion, rapid T1 mapping allows for a characterization of the axonal architecture in an individual CC in less than 10 s. The approach emerges as a valuable means for studying neocortical brain anatomy with possible implications for the diagnosis of neurodegenerative processes.

Off-resonant spinlock (SL) enables an NMR imaging technique that can detect dilute metabolites similar to chemical exchange saturation transfer. However, in clinical MR scanners, RF pulse widths are restricted due to recommended specific absorption rate limits. Therefore, trains of short RF pulses that provide effective saturation during the required irradiation period are commonly employed. Quantitative evaluation of spectra obtained by pulsed saturation schemes is harder to achieve, since the theory of continuous wave saturation cannot be applied directly. In this paper we demonstrate the general feasibility of quantifying proton exchange rates from data obtained in pulsed SL experiments on a clinical 3 T MR scanner. We also propose a theoretical treatment of pulsed SL in the presence of chemical exchange using an interleaved saturation-relaxation approach. We show that modeling magnetization transfer during the pauses between the RF pulses is crucial, especially in the case of exchange rates that are small with respect to the delay times. The dynamics is still governed by a monoexponential decay towards steady state, for which we give the effective rate constant. The derived analytical model agrees well with the full numerical simulation of the Bloch-McConnell equations for a broad range of values of the system parameters., (Copyright © 2014 John Wiley & Sons, Ltd.)