Optimized Inversion-Time Schedules For High-Resolution Multi-Inversion EPI Quantitative Measurements of T1
Ouri Cohen1,2, Ville Renvall3, and Jonathan Polimeni1,2

1Athinoula A. Martinos Center, Charlestown, MA, United States, 2Radiology, Massachusetts General Hospital, Boston, MA, United States, 3Department of Neuroscience and Biomedical Engineering, Aalto University School of Science, Espoo, Finland

Synopsis

A novel optimized method for high-resolution quantitative EPI measurements of T1 is introduced and validated on a 3T clinical scanner in a phantom and a healthy volunteer. The method offers a 5-fold acceleration in scan time over previous techniques allowing fully quantitative 1.2 mm3 isotropic T1 maps in less than 30 seconds.

Introduction

We have previously demonstrated a fast quantitative, T1 mapping method based on inversion recovery echo planar imaging (IR-EPI) [1]. The method relies on the acquisition of multiple EPI slices following a broadly slab-selective inversion. The slice acquisition order is then permuted by a fixed ‘skip-factor’ for each of the subsequent inversion recovery periods. By rapidly acquiring multiple inversion times (TI) for each slice, accurate T1 maps can be obtained. Nevertheless, the reordering scheme used is empirically chosen and may not be optimal for the desired range of tissue parameters. Rather than a fixed skip-factor, our goal was to find an optimal reordering scheme, i.e. one that maximizes the discrimination between different tissues. In this work we demonstrate that optimizing the reordering scheme permits a 5-fold reduction in the number of inversions required resulting in significant scan time savings without adversely affecting the quantitative maps obtained.

Methods

To find the optimal reordering we simulated the magnetization for a given slice ordering of all possible tissue T1 values. The dot product matrix was then computed and a similarity metric was defined as the sum of all off-diagonal elements, similarly to [2]. A genetic optimization algorithm [3] was then used to find an ordering that minimized the chosen similarity metric. The original skip-factor ordering with 21 inversions and a skip-factor of 7 was used as a baseline. To minimize computation time, the optimization was done assuming 4 inversions and was only allowed to run for 30 generations (iterations). All experiments were conducted on a 3T Siemens Trio whole body scanner (Siemens Healthcare, Erlangen, Germany) equipped with a 32-channel phased-array head coil. A phantom composed of vials with different concentrations of doped water was used to mimic the variety of T1 values in the brain. The phantom was scanned at a resolution of 1.2 mm3 using both the unoptimized (21 inversions) and the optimized (4 inversions) ordering with the following acquisition parameters: TR/TE/flip angle/BW/matrix size/slices/R = 6690 ms/26.9 ms/ 90°/ 1240 Hz/pixel / 192×192/100/4 with FLEET autocalibration [3] (with α=10° and 5 preparatory pulses) and online GRAPPA reconstruction. Total scan time for the IR-EPI data was 140/27 seconds for the unoptimized/optimized sequence, plus a fixed 30 s for GRAPPA autocalibration. The acquired data was processed offline in Matlab. The signal in each voxel was compared to a pre-computed dictionary of signal magnetizations and the best match selected [3]. The mean of the T1 in each vial obtained using both orderings was used to calculate a best-fit curve whose R2 value was used to measure the quantitative similarity between the two orderings. A healthy 24 y.o. female volunteer was recruited and provided written informed consent. The subject was scanned with both slice orderings.

Results

A visual comparison between the skip-factor 7 unoptimized and optimized slice orderings is shown in Figure 1. A comparison between the T1 values obtained with each ordering is shown in Figure 2 along with a best-fit curve whose R2 value was 0.9994 indicating the closeness of the quantitative values obtained with the two orderings. The standard deviation for the optimized ordering, however, is substantially smaller. Representative images from the phantom for both acquisitions are shown in Figure 3, showing that the T1 values match known values for the solutions. Maps from the human subject for both slice ordering are shown in Figure 4, where the T1 values approximate reported values for brain tissue at 3T [5] and are similar to those of the unoptimized ordering.

Discussion

Despite its clinical importance, quantitative tissue mapping is hindered by the long scan time required. Our optimized IR-EPI method requires less than 30 seconds for fully quantitative, 1.2 mm3 isotropic T1 maps representing a 5-fold saving over previous techniques. Similar protocols have been successfully used at higher (7T) fields making this method easily portable to higher fields. Current limitations include the long processing time required for the optimization but may be resolved by use of improved optimization algorithms and parallel processing.

Conclusion

An optimization method to reduce scan time in IR-EPI was developed and tested on a clinical 3-T scanner. Potential applications for the increased acceleration include clinical scanning and imaging of dynamic processes such as dynamic contrast enhancement [6].

Acknowledgements

Supported by NIH NIBIB K01-EB011498, P41-EB015896, and R01-EB019437 and the Athinoula A. Martinos Center for Biomedical Imaging, and made possible by NIH NCRR Shared Instrumentation Grants S10-RR023401 and S10-RR020948.

References

[1] Renvall et al, ISMRM 2014 #4282 [2] Cohen et al, ISMRM 2014 #0027, [3]Polimeni et al, MRM, 2015; PMID 25809559 [4] Haupt, R.L. and Haupt, S.E. Practical Genetic Algorithms , 2003, Wiley & Sons [5] Ma et al, Nature 2013; 495:187-192 [5] http://www.ncbi.nlm.nih.gov/pubmed/18259791 [6] Kalpathy-Cramer et al, ISMRM 2015 #4392.

Figures

Figure 1: The unoptimized (a) and optimized (b) slice ordering used. The unoptimized ordering used a skip-factor 7 and 21 inversions. Note that, in contrast, the optimized ordering required only 4 inversions.

Figure 2: Mean and standard deviation of the vial phantom T1 values obtained using each ordering method.The best-fit curve plotted had an R2 of 0.9994, indicative of the closeness between the values obtained by the two orderings. Note that despite the 5-fold smaller number of inversions, the optimized ordering yielded standard deviations that were substantially smaller.

Figure 3: Representative sample images acquired using the optimized and unoptimized orderings. Note the similarity between the maps obtained.

Figure 4: Sample in vivo slices from the optimized and unoptimized acquisitions. The quantitative T1 values obtained approximated reported values for brain tissue at 3T and are similar to those of the unoptimized ordering.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0334