Towards Judging the Encoding Capability of Magnetic Resonance Fingerprinting Sequences
Karsten Sommer1, Thomas Amthor1, Peter Koken1, Mariya Doneva1, and Peter Börnert1

1Philips Research Europe, Hamburg, Germany

Synopsis

A key question of magnetic resonance fingerprinting (MRF) is the appropriate choice of sequence parameters to achieve a high sensitivity to the tissue parameters of interest. In this contribution, different candidates for a measure of MRF sequence encoding capability are evaluated. While interpretation of measures that rely on local or global dot products proved difficult, a ‘brute force’ Monte Carlo approach showed good agreement with experimental results. By restricting this Monte Carlo method to small local dictionaries, substantial acceleration could be achieved.

Purpose

Magnetic resonance fingerprinting (MRF) is a new technique promising the efficient quantification of tissue parameters1. A fundamentally unresolved issue of MRF, however, is the appropriate choice of sequence parameters to achieve a high sensitivity to the tissue parameters of interest. In particular, a fast and robust quality measure to determine this encoding capability is desirable to identify the ideal set of sequence parameters for a specific application. In this contribution, we compare and evaluate different candidates for such an MRF sequence performance measure.

Methods

Five different methods were considered to measure the sequence encoding capability. For the evaluation of these approaches, dictionaries of SSFP-like sequences2 (T1 range: 200-1500 ms, ΔT1=10 ms; T2 range: 50-400 ms, ΔT2=5 ms) with an inversion pulse (TI=20 ms) followed by randomly generated flip angle (FA) patterns (range: 0‑70°, 1000 sequence steps, Fig. 1) were calculated. To reduce the dimensionality of the sequence parameter space, a fixed repetition time (TR) of 15 ms was employed. The five methods were: (i) The average ‘local’ dot product of two consecutive dictionary entries with ΔT2=0 ms and ΔT1=10 ms calculated at fifty equidistant locations in the dictionary. For a high encoding capability to T1, this average dot product should be as low as possible. (ii) The same as method (i) using two consecutive entries with ΔT1=0 ms and ΔT2=5 ms. (iii) A matrix containing all ‘global’ dot products of fifty equidistant entries in the dictionary was generated, and the maximum of this matrix was calculated. (iv) ‘Full‘ Monte Carlo (MC) method: Gaussian noise (SNR~10) was added to a dictionary entry, and the resulting signal was matched to the dictionary; this was repeated 2000 times with randomly resampled noise for the same entry. This method was applied for fifty equidistant entries, and the dictionary-wide average error was calculated. (v) Finally, acceleration of the full MC method was realized by using local dictionaries consisting only of the closest neighbors (ΔT1=10 ms and/or ΔT2=5 ms) of the examined entries. These local dictionaries were expanded in T1 and/or T2 direction, until the cumulative matching probability of the outermost entries was below 2%. To evaluate the performance of these measures, an MRF measurement of a gel phantom containing 12 substances with different T1/T2 combinations in the range of the dictionary was performed on a Philips Achieva 1.5 T scanner with an 8-channel head receiver array using the FA pattern in Fig. 1. Afterwards, reduced sequence lengths were simulated by shortening both the acquired data and the dictionary in time.

Results and Discussion

Figures 2a and b depict the mean T1 and T2 values in the 12 tubes of the phantom as a function of sequence length, respectively. The results of methods (i-iii) for the same sequence, as depicted in Fig. 2c, reveal a fundamental problem of the dot product based techniques: due to the normalization of the dictionary entries, absolute signal values of the entries generally decrease with increasing sequence length. Consequently, methods (i) and (iii) predict a decreasing accuracy of the MRF sequence towards higher sequence lengths, which is inconsistent with the experimental results. In contrast, the relative matching errors that were obtained using the full MC method (iv) (Fig. 2d) were consistent with the phantom measurements: the convergence to the final T1 and T2 values (Figs. 2a,b) coincided with sharp decreases of the relative matching errors as indicated by the vertical lines. Figure 3 shows a comparison of the considered methods for the different sequences. While the local dot product based methods (i/ii) yield measures somewhat similar to the full MC method both for T1 and T2, the curve shape for the global dot product method (iii) was markedly different. The accelerated MC method successfully reproduced the matching errors of the full MC method (max. deviation < 0.1 percentage points). On a 6-core Intel Xeon 2.67 GHz system, the evaluation of a single sequence using this method required about 3 to 4 min. Compared with the full MC method, this represented a reduction of the computation time by a factor > 12 (Fig. 3d).

Conclusion

Fast and accurate evaluation of the quality of a given sequence in MRF is generally challenging. While interpretation of the dot product quality measures can be difficult in some cases, the presented Monte Carlo method is able to accurately predict the minimum required duration of the sequence. By restricting this method to small local dictionaries, substantial acceleration could be achieved. This technique may hence be applied to increase the sensitivity of MRF to tissue parameters by optimizing the employed sequence.

Acknowledgements

No acknowledgement found.

References

[1] Ma D et al., Nature 495, 187 (2013) [2] Jiang Y et al., Magn Res Med (2014)

Figures

Fig. 1: The FA pattern that was employed for the phantom measurement.

Fig. 2: Comparison of phantom experiments and encoding capability measures as a function of sequence length. (a) The mean T1 values in the 12 tubes of the phantom as a function of sequence length. (b) The same for T2. (c) The quality measures for methods (i-iii). (d) The relative T1/T2 matching errors for method (iv).

Fig.3: Variation of the quality measures for sequences with different FA patterns. (a) Relative T1 matching errors of all 10 sequences for the full and the accelerated MC as well as the local T1 dot product method. (b) The same for T2. (c) The results of the global dot product method. (d) The speedup in computation time for the accelerated MC method compared to the full MC method.



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