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 parameters
1. 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 sequences
2 (T
1 range:
200-1500 ms, ΔT
1=10 ms; T
2 range: 50-400 ms, ΔT
2=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 ΔT
2=0 ms
and ΔT
1=10 ms calculated at fifty equidistant locations in the
dictionary. For a high encoding capability to T
1, this average dot
product should be as low as possible. (ii) The same as method (i) using two
consecutive entries with ΔT
1=0 ms and ΔT
2=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 (ΔT
1=10 ms
and/or ΔT
2=5 ms) of the examined entries. These local dictionaries
were expanded in T
1 and/or T
2 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 T
1/T
2
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 T
1 and T
2 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 T
1
and T
2 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 T
1 and T
2, 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)