Fast Interleaved Multislice T1 Mapping: Model-Based Reconstruction of Single-Shot Inversion-Recovery Radial FLASH
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 mm3 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 mm3 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.
Quantitative T1 mapping is of great interest for application in both research and diagnostics . Fast T1 mapping usually employs an inversion-recovery (IR) look-locker (LL) sequence where RF pulses are continuously applied after inversion with T1 determination in a postprocessing step [2, 3]. Based on the IR LL sequence, recent advances in real-time MRI [4, 5] and model-based reconstructions [6–8] enabled single-shot T1 mapping within a few seconds. The first method  defines the IR readout by a series of highly undersampled radial fast low-angle shot (FLASH) images with nonlinear inversion (NLINV) reconstruction , while the latter methods reconstruct all parameter maps directly from highly undersampled k-space data, bypassing the intermediate step of image reconstruction completely. Compared to image-space methods such as NLINV, model-based reconstructions  allow for a more efficient exploitation of redundancy in the data. Thus, model-based approaches have also been investigated for accelerated T2 mapping [11–13], diffusion , flow , and MR fingerprinting .
In general, quantitative parametric mapping so far focused on single slices, whereas in clinical practice, multislice mapping is highly desirable. For example, it has been recommended to perform myocardial T1 mapping in at least three short-axis sections to capture potential heterogeneity across the left ventricular wall . Previous studies have used a variable flip angle (VFA) approach  to acquire volume T1 maps within a clinically acceptable time. While efficient, this method is sensitive to B1 field inhomogeneity and motion might be a problem for abdominal applications [19, 20]. Compared to the VFA approach, the IR LL technique is more tolerant to B1 field inhomogeneities [2, 3]. Early efforts based on the IR LL sequence usually employed segmented Cartesian acquisition schemes which involve a waiting period between the steady state of one segment and the inversion of the next segment. Although different strategies [21, 22] have been proposed to sample the relaxation process more efficiently, it still takes minutes to achieve multislice T1 maps. Therefore, a spiral acquisition with through-time generalized autocalibrating partial parallel acquisitions  and the use of spatiotemporal sparsity constraints  have been proposed to accelerate volumetric T1 mapping for abdominal applications to about 20 s breath holds. More recently, radial spoke-sequential FLASH with NLINV reconstruction resulted in 5 T1 maps of the human brain within a single IR . However, because all these methods rely on reconstructions in image space, they have to provide a sufficient number of good-quality images during IR in order to allow for subsequent T1 fitting. This condition not only enforces restrictions on the choice of k-space trajectories sampling the entire IR but also renders image-space methods unlikely to fully exploit the inherent redundancy in the data (i.e., about 1000 radial spokes for a 4 s acquisition). This becomes obvious when considering that image-space methods have to reconstruct up to 80 images for a single IR, while all what is needed are only three quantitative maps (i.e., equilibrium magnetization, effective flip angle, and T1).
This study therefore aimed at developing an efficient model-based reconstruction technique which takes advantage of the inherent data redundancy by allowing for simultaneous (i.e., spoke-interleaved) multislice T1 mapping at high spatial resolution. The approach combines a single-shot IR multislice radial FLASH acquisition and a model-based nonlinear inverse reconstruction with joint sparsity constraints. Apart from validations using numerical simulations and experimental phantoms, the performance of the proposed method is demonstrated for studies of the human brain and liver.
2.1. Single-Shot Multislice Data Acquisition
The data acquisition scheme is demonstrated in Figure 1. Multislice radial k-space data are continuously acquired after a single nonselective inversion pulse using a small-angle golden-angle sampling scheme [25, 26] in a spoke-interleaved manner. The interslice distance is chosen to be at least the slice thickness to avoid cross talk between slices. The relaxation process for each individual section follows the three-parameter model according to the Bloch equations as follows [2, 3]:where represents the steady-state magnetization, is the equilibrium magnetization, is the position in image space, and is the effective relaxation rate given bywith the repetition time between successive RF pulses as shown in Figure 1, the flip angle, and the total number of slices. The sought quantity T1 is then calculated according to the following exact formula after estimation of , , and :
2.2. Model-Based Reconstruction
As described in , the signal received from multiple receiver coils can be written aswith the jth coil sensitivity map, the chosen k-space trajectory, the prescribed model in (1), and the acquired data. Both the parameter maps and coil sensitivity maps are directly estimated from the acquired k-space data using a single-step model-based reconstruction, that is,where F is the nonlinear forward model mapping all unknowns to the measured data y. and represent the parameter maps and the coil sensitivity maps, respectively. is a L1 regularization to exploit joint sparsity in the parameter maps following the ideas of compressed sensing, while is the Sobolev norm on the coil sensitivities to enforce its intrinsic smoothness. and are the corresponding regularization parameters. D is a convex set representing the acceptable domain for the unknowns. Similar to , D is set to be as the relaxation rate should be nonnegative. The nonlinear inverse problem in (5) is solved by the iteratively regularized Gauss–Newton method (IRGNM)  where the nonlinear problem is linearized in each Gauss–Newton step and solved by the fast iterative shrinkage-thresholding algorithm .
2.3. Numerical Simulations
The numerical phantom designed in  was used for validations. T1 relaxation times of the phantom are 300 ms, 800 ms, 1500 ms (circular objects), and 2000 ms (background). The k-space data were derived from the analytical Fourier representation of an ellipse assuming an array of four circular receiver coils surrounding the phantom without overlap . The simulations employed the spoke-interleaved multislice IR FLASH sequence with the acquisition parameters as listed in Table 1. Complex white Gaussian noise with a standard deviation of 0.1 was added to the simulated k-space data.
All MRI measurements were performed on a human 3T MRI system (Magnetom Prisma fit, Siemens Healthineers, Erlangen, Germany). Phantom and brain studies were conducted with a standard 64-channel head coil, while abdominal scans were performed with an 18-element thorax coil in conjunction with 18 elements of the 32-element spine coil. During technical developments, 6 subjects without known illness were recruited among the students of the local university. Written informed consent, according to the recommendations of the local ethics committee, was obtained from all subjects prior to MRI.
The proposed method was experimentally validated with a commercial reference phantom (Diagnostic Sonar LTD, Scotland, UK) consisting of 6 compartments with defined T1 values surrounded by water. Reference T1 maps are taken from the single-slice acquisition whose accuracy and precision have been validated relative to T1 maps calculated from fully sampled data sets . All the acquisition parameters are summarized in Table 1. For IR FLASH, B0 shimming and center-frequency adjustment were performed using the standard procedures provided by the manufacturer. Any residual transverse magnetization was spoiled by random RF phases . Abdominal studies were performed during a brief breath hold.
In the spoke-interleaved multislice acquisition scheme, the RF repetition time increases for any given slice compared to a single-slice measurement (2). Accordingly, the flip angle can also be increased within a certain range to optimize the overall signal-to-noise ratio (SNR). However, with increasing flip angles, T1 accuracy may suffer as the applied correction formula loses validity . Therefore, we varied the flip angles systematically based on experimental phantom studies. An optimal flip angle is then chosen that balances T1 accuracy and SNR for each acquisition. The validations showed that flip angles of 10°, 12°, and 14° turned out to be optimal for 3-slice, 5-slice, and 7-slice acquisitions, respectively. These values were used for all subsequent studies.
Model-based reconstruction techniques in general do not attempt to reconstruct images, and consequently, data binning prior to reconstruction is not necessary. However, a certain degree of temporal discretization may effectively reduce the computational demand without posing too rigid restrictions on probing of the IR. Here, the number of binned spokes was chosen such that the temporal bin size does not exceed 80 ms and T1 accuracy is not compromised . Small golden angles  were then adjusted for each data acquisition to ensure a homogeneous coverage of each k-space.
At this stage, all data processing was done offline. Multicoil raw data for each slice were first corrected for gradient delays  and then compressed to 8 virtual channels using a principle component analysis. A convolution-based gridding  without density compensation was used to interpolate the radial samples onto a Cartesian grid on which all successive computations were performed. Gradient delay correction, channel compression, and gridding were done in MATLAB (MathWorks, Natik, MA), while the iterative optimization was implemented in C/CUDA using GeForce GTX TITAN (NVIDIA, Santa Clara, CA).
Regularization parameters α and are initially set to 1 and subsequently reduced by a factor of 3 in each Gauss–Newton step. A minimum value of α was introduced to control the noise in higher Gauss–Newton steps. The chosen value of = 0.0015 for applications to the brain was defined by optimizing SNR without compromising quantitative accuracy or delineation of structural details. With similar settings, = 0.001 was chosen for abdominal studies. Constants in the Sobolev norm were the same as in . As was insensitive to the final results, no minimum value was set for . Similar to , 10 Gauss–Newton steps were used for IRGNM to ensure convergence.
Figure 2 shows T1 maps of a numerical phantom obtained by model-based reconstruction from single-slice, 3-slice, 5-slice, and 7-slice acquisitions. The total number of spokes per slice decreased from 1064 to 364, 225, and 156, respectively. Visual inspection of the T1 and difference maps reveals no detectable differences between single-slice and multislice acquisitions. Quantitative ROI analyses in Table 2 confirm the T1 accuracy of all multislice acquisitions, while T1 precision slightly reduces (i.e., standard deviations slightly increase) for increasing number of slices and mostly for shorter T1 values.
Figure 3 and Table 3 offer a similar comparison of single-slice and multislice results for an experimental T1 phantom comprising six tubes with different T1 values. All measurements represent single-shot IR acquisitions with a duration of 4 s. Again, there is no visual difference across T1 and difference maps, while Table 3 confirms excellent T1 accuracy and precision (i.e., small standard deviations) for all measurements and the entire range of T1 values.
Figure 4 compares high-resolution T1 maps, enlarged views, and difference maps of the human brain (same section) obtained by single-slice and multislice acquisitions in 4 s. In line with phantom validations, multislice T1 maps are of comparable quality to the single-slice map. This is best demonstrated by the enlarged T1 maps (Figure 4, middle row) and the quantitative analyses in Table 4. Figure 5 summarizes all 7 maps of a 7-slice T1 mapping experiment using the proposed model-based reconstruction. Similarly, Figure 6 presents the results of a 3-slice T1 mapping study of the liver obtained within a single 4 s breath hold. The magnified views clearly demonstrate good spatial definition of structural details as well as adequate SNR, that is, good precision. The quantitative results yield the regional liver T1 values of 775 ± 34 ms, 783 ± 26 ms, and 766 ± 29 ms (analyzed on regions-of-interests (ROIs) of 144 pixels per T1 map) for the three sections in agreement with literature data of 767 to 812 ms [38, 39].
This work presents a novel multislice T1 mapping technique which combines a single-shot IR spoke-interleaved multislice FLASH acquisition with a sparsity-constrained model-based nonlinear inverse reconstruction. In contrast to mapping techniques based on image space with subsequent fitting, the combined model-based approach effectively exploits the data redundancy in single-slice IR acquisitions. In the present study, this refers to the acquisition of about 1000 radial spokes per single-shot 4-second IR which are available for the reconstruction of only three-parametric maps rather than for up to 80 images (both cases also include all coil sensitivity maps). The present results demonstrate that the radial information available in a single-shot IR acquisition may effectively be used for multislice T1 mapping without compromising resolution, accuracy, and precision as validated by numerical simulations and experimental phantoms.
While 7 simultaneous T1 maps were achievable for the brain, only 3 T1 maps of the liver were presented. This is mainly because of the more pronounced presence of off-resonance contributions (e.g., air-tissue susceptibility differences) within the abdomen which, at higher degrees of undersampling for 5-slice or 7-slice acquisitions, eventually cause residual streaking artifacts. Another influential factor is the need for a higher number of spokes to avoid streaks caused by objects outside the chosen field-of-view, for example, arms lying next to the body . Both effects are also taken care off by sampling the IR with more spokes in a 3-slice acquisition.
Simultaneous multislicing (SMS) is another valuable technique to accelerate multislice imaging in general  as well as multislice parametric mapping [42, 43]. For example, a combination of SMS and magnetic resonance fingerprinting achieved three simultaneous T1 maps within 12 s . SMS has also been shown to accomplish 3-slice myocardial T1 mapping . In this respect, the combination of the present model-based reconstruction and SMS would be an interesting next step. It might further increase the number of T1 maps or offer access to even higher spatial resolution.
At this stage, a major limitation of the proposed method is the long offline computational time. This seems to be a general problem for model-based reconstructions (at least for T1 mapping) as all the data have to be kept in memory during iterative optimization. However, several acceleration techniques are under development, including the extension of the current single-GPU version to a multi-GPU implementation or integration of model-based reconstructions into Berkeley Advanced Reconstruction Toolbox .
The proposed method offers rapid and robust single-shot multislice T1 mapping with high accuracy and precision by combining the advantages of model-based nonlinear inverse reconstruction, radial sampling, parallel imaging, and joint sparsity constraints. This novel approach effectively exploits the inherent redundancy of a single-shot spoke-interleaved IR LL experiment. It now warrants clinical evaluations as well as extensions to other applications such as single-shot multislice myocardial T1 mapping.
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
Drs. Frahm and Uecker are coinventors of a patent describing the real-time MRI technique used here. The other authors declare no conflicts of interest.
All the authors contributed to the development of the method and the design of the study. Xiaoqing Wang and Volkert Roeloffs performed the experiments and data analysis. Dirk Voit performed the sequence implementation. All the authors contributed to the preparation of the manuscript.
X. Wang, A. A. Joseph, O. Kalentev et al., “High-resolution myocardial T1 mapping using single-shot inversion-recovery fast low-angle shot MRI with radial undersampling and iterative reconstruction,” British Journal of Radiology, vol. 89, no. 1068, Article ID 20160255, 2016.View at: Publisher Site | Google Scholar
J. C. Moon, D. R. Messroghli, P. Kellman et al., “Myocardial T1 mapping and extracellular volume quantification: a Society for Cardiovascular Magnetic Resonance (SCMR) and CMR Working Group of the European Society of Cardiology consensus statement,” Journal of Cardiovascular Magnetic Resonance, vol. 15, no. 1, p. 92, 2013.View at: Publisher Site | Google Scholar
R. Treier, A. Steingoetter, M. Fried, W. Schwizer, and P. Boesiger, “Optimized and combined T1 and B1 mapping technique for fast and accurate T1 quantification in contrast-enhanced abdominal MRI,” Magnetic Resonance in Medicine, vol. 57, no. 3, pp. 568–576, 2007.View at: Publisher Site | Google Scholar
X. Wang, V. Roeloffs, K. D. Merboldt, D. Voit, S. Schätz, and J. Frahm, “Singleshot multi-slice T1 mapping at high spatial resolution–inversion recovery FLASH with radial undersampling and iterative reconstruction,” Open Medical Imaging Journal, vol. 9, no. 1, pp. 1–8, 2015.View at: Publisher Site | Google Scholar
A. B. Bakushinsky and M. Y. Kokurin, Iterative Methods for Approximate Solution of Inverse Problems, Springer, Dordrecht, Netherlands, 2004.
K. T. Block and M. Uecker, “Simple method for adaptive gradient-delay compensation in radial MRI,” in Proceedings of 19th Annual Meeting of ISMRM, p. 2816, Montreal, QC, Canada, May 2011.View at: Google Scholar
F. Wajer and K. P. Pruessmann, “Major speedup of reconstruction for sensitivity encoding with arbitrary trajectories,” in Proceedings of 9th Annual Meeting ISMRM, p. 767, Glasgow, UK, April 2001.View at: Google Scholar
N. Gelman, J. R. Ewing, J. M. Gorell, E. M. Spickler, and E. G. Solomon, “Interregional variation of longitudinal relaxation rates in human brain at 3.0 T: relation to estimated iron and water contents,” Magnetic Resonance in Medicine, vol. 45, no. 1, pp. 71–79, 2001.View at: Publisher Site | Google Scholar
M. Uecker, F. Ong, J. Tamir et al., “Berkeley advanced reconstruction toolbox,” in Proceedings of 23rd Annual Meeting of ISMRM, p. 2486, Toronto, ON, Canada, May-June 2015.View at: Google Scholar