5220

Initial implementation of magnetic resonance fingerprinting on a preclinical 14.1 T scanner
Yasuhiko Terada1

1University of Tsukuba, Tsukuba, Japan

Synopsis

Magnetic resonance fingerprinting (MRF) is a technique that enables simultaneous quantification of tissue parameters in a single scan. So far, MRF has been realized mostly for clinical 1.5 T and 3 T scanners, and a preclinical 7 T scanner. Application to a higher field scanner is limited and challenging because of the increased sensitivity to system hardware imperfections, such as B0 and B1 inhomogeneities and slice-profile imperfection. In this study, we assessed the feasibility of the MRF approach in a preclinical, wide-bore 14.1 T system.

Introduction

Magnetic resonance fingerprinting (MRF)1 is a promising technique that enables multiparameter mapping in a single scan. So far, MRF has been realized mostly for clinical 1.5 T and 3 T scanners, and a preclinical 7 T scanner2. Application to a higher field scanner is limited and challenging because of the increased sensitivity to system hardware imperfections, such as B0 and B1 inhomogeneities and slice-profile imperfection. In this study, we examined and confirmed the feasibility of the MRF approach for a preclinical, wide-bore 14.1 T scanner.

Methods

The MRF was implemented on the 14.1 T MRI system (Bruker AVANCE III 600 MHz Wide Bore NMR/MRI Spectrometer, Bruker Inc) with a 40 mm inner diameter gradient insert with up to 1.47 T/m gradient strengths (Bruker Micro2.5 gradient system) and a 30 mm inner diameter quadrature coil (MicWB40, Bruker Inc.). The samples were CuSO4–doped water phantoms with varied concentrations and a tomato fruit.

The MRF acquisition was developed from a fast imaging with steady-state free precession (FISP) 2,3. The MRF acquisition was initiated with an adiabatic inversion pulse (4.84 ms duration), which was followed by 300 successive FISP acquisition periods with varying flip angel (FA) and repetition time (TR). Echo time (TE) was held constant (4 ms). FA and TR were determined by pseudo-randomly according to the method proposed by Gao et al.3. FA was varied sinusoidally (5 to 90o) and TR was varied using a Perlin noise pattern (20 to 30 ms). For the acquisition of the MRF signal evolutions, the same line of k space for each of 300 sequential images during each MRF scan repetition period (whole TR) was measured.

The MRF dictionary was created using a home-written extended phase graph simulator, which was implemented in a C/C++ programming environment. The MRF dictionary consisted of 137280 profiles generated from 24 T1 values (20-600 ms), 44 T2 values (20-600 ms: increment, 20 ms) and 25 B1 values for the CuSO4-doped water measurement and 218736 for the tomato measurement. Unrealistic profiles with T1<T2 were excluded. The dictionary was matched to the acquired MRF signal evolution profile using vector-based inner product comparisons.

Results and discussion

The MRF-FISP results varied, depending on the slice thickness (ST) and spoiler gradient moment (Fig. 1). Although the MRF-T1 estimates were not significantly different, the MRF-T2 estimates were different from each other. The MRF-T2 map suffered from the banding artifacts for ST = 0.4 mm and the smaller moment. These artifacts, which also appeared in the trueFISP image, were caused by the field inhomogeneity. In general, the banding artifacts tend to be remarkable in the high-field-strength and wide field-of-view scanners. Under the current experimental setup, the field homogeneity needed to be no lower than ~0.03 ppm over a 30 mm diameter spherical volume to avoid banding artifacts. For the system we used, this could not be satisfied even with the automated, localized third-order shimming.

As shown in Fig. 1, the banding artifacts could be alleviated by increasing the spoiler gradient moment, but it also caused a decrease in the T2 value possibly because of the diffusion weighting. Therefore, there is a tradeoff between the banding-artifact alleviation and T2 accuracy. For ST = 1.2 mm, the MRF-T2 map showed another artifacts, possibly because of the large B0 inhomogeneity in a voxel. The MRF-T2 estimate also depended on the slice profile (Fig. 2). The banding artifacts less appeared for the profiles with the smaller time-bandwidth products (TBPs), possibly because spins in one voxel were less coherent and they were effectively dephased.

Figure 3 shows the MRF results obtained under an optimum condition (ST = 0.8 mm; moment = 8π/0.8 mm; SLR3 pulse). Figure 4 shows the MRF-T1 and T2 maps of a tomato fruit, which was obtained under the above optimum condition. There were less banding artifacts in the images and the MRF-T2 estimates corresponded to the spin-echo based estimates. The MRF-T1 estimates were smaller than the inversion-recovery based estimates, possibly because the whole TR was short and T1 saturation was not fully recovered. This could be corrected by selecting a longer whole TR.

Conclusion

We examined the feasibility of the MRF-FISP approach for the wide-bore 14.1 T system. We tested different slice thicknesses, spoiler gradient moments, and slice profiles, and concluded that the MRF-FISP is feasible under certain conditions.

Acknowledgements

The author would like to thank Prof. William S. Price for his kind permission to use the MRI system and Dr. Allan Torres and Dr. Mikhail Zubkov for their help to use the software.

References

1. Ma D, et al. Magnetic resonance fingerprinting, Nature 2013;495:187-193.

2. Gao Y, et al. Preclinical MR fingerprinting (MRF) at 7 T: effective quantitative imaging for rodent disease models, NMR Biomed. 2015;28:384-394.

3. Jiang Y, et al. MR Fingerprinting Using Fast Imaging with Steady State Precession (FISP) with Spiral Readout, Magn. Reson. Med. 2015;74(6):1621-31.

Figures

Fig. 1: MRF-T1 and T2 maps for different slice thicknesses (STs). ST = 0.4, 0.8, and 1.2 mm. FOV = 30 mm x 30 mm. Matrix size = 128 x 128. Whole TR for MRF = 10 s. The trueFISP image was measured with TE/TR = 12.5 ms/25 ms. For comparison, T1 and T2 maps obtained with standard methods (SE-T1 and SE-T2) were shown. The SE-T1 map was generated from inversion recovery spin echo images. The SE-T2 map was generated from spin echo images.

Fig. 2: MRF-T2 maps for different slice profiles. SLR: Shinnar-Le Roux pulse; PW: RF pulse width; TBP: time bandwidth product.

Fig. 3: MRF-based T1, T2, proton density (PD), and B1 maps obtained under the optimum condition for the CuSO4 phantom. Slice thickness = 0.8 mm. FOV = 30 mm x 30 mm. Matrix size = 128 x 128. Whole TR for MRF = 10 s. The SE-T1 and SE-T2 maps were the same as shown in Fig. 1.

Fig. 4: T1 and T2 maps generated from MRF and SE for a tomato fruit. FOV = 30 x 30 mm. Matrix size = 128 x 128. Whole TR = 20 s. Inversion recovery spin echo images (inversion time = 500, 1000, 2000, and 3000 ms; TE/TR = 10 ms/10 s) were used to generate SE-T1 map. Spin echo images (TE = 10, 50, 100, and 200 ms; TR = 10 s) were used to generate SE-T2 map.

Proc. Intl. Soc. Mag. Reson. Med. 25 (2017)
5220