0692

Fast high-resolution isotropic whole-heart T2 mapping using focused navigation
Simone Rumac1, Christopher W. Roy2, Jérôme Yerly2,3, John Heerfordt2,4, Davide Piccini2,4, Matthias Stuber3,5, and Ruud B. van Heeswijk2
1Department of Radiology, Department of Radiology, Lausanne University Hospital (CHUV) and University of Lausanne (UNIL), Laus, Lausanne, Switzerland, 2Department of Radiology, Department of Radiology, Lausanne University Hospital (CHUV) and University of Lausanne (UNIL), Lausanne, Switzerland, Lausanne, Switzerland, 3CIBM Center for BioMedical Imaging, Lausanne, Switzerland, Lausanne, Switzerland, 4Advanced Clinical Imaging Technology, Siemens Healthcare AG, Lausanne, Switzerland, Lausanne, Switzerland, 5Department of Radiology, Department of Radiology, Lausanne University Hospital (CHUV) and University of Lausanne, Lausanne, Switzerland

Synopsis

Cardiac parametric mapping techniques are gaining traction for the clinical routine assessment of various pathologies. Despite the complex 3D patterns of many myocardial conditions, most current techniques are breath-held single-slice 2D acquisitions. We propose a free-breathing high-resolution isotropic 3D T2 mapping technique for the heart where breathing motion is corrected in k-space before image reconstruction. In 4 healthy volunteers and one patient with myocardial infarction, we found that our technique produced sharp and accurate T2 maps but had slightly lower precision than routine techniques.

Introduction

T2 mapping of the myocardium is gaining acceptance as a routine technique for the assessment of acute edema.1 Even though most myocardial pathologies have unpredictable 3D patterns throughout the myocardium, most current clinical techniques are breath-held single-slice 2D acquisitions. To characterize these complex disease patterns, several 3D T2 mapping techniques have been proposed.2–5 These techniques usually come with a sacrifice in either scan duration, precision, or motion robustness. In order to avoid these compromises, we here propose an isotropic 3D T2 mapping that 1) robustly corrects for respiratory motion through autofocusing, 2) produces precise maps through compressed sensing (CS) with combined local low-rank (LLR) and total variation (TV) regularization and 3) has high accuracy through dictionary fitting of the exact time course of the pulse sequence.

Methods

The acquisition (Figure 1) consisted of a prototype ECG-triggered free-breathing 3D radial bSSFP sequence with a phyllotaxis trajectory,6 field of view=(220mm)3, isotropic spatial resolution=(1.5mm)3, TR/TE=3.4/1.7ms, α=40°, acquisition window=150ms, five T2-preparation times TE-T2prep=0/25/40/55/70ms, 3420 readouts/T2prep, and a total acquisition time of 5×76 heartbeats=~6min. All data were acquired on a 1.5T clinical scanner (MAGNETOM Sola, Siemens Healthcare, Erlangen, Germany). Respiratory motion was corrected using a recently developed autofocusing technique that derives 3D respiratory motion from a periodically acquired radial readout called focused navigation (fNAV).7 Source images were reconstructed with LLR regularization along the T2prep dimension (λ=0.015) and TV denoising (λ=0.05). For each scan, Bloch equation simulations were used to generate a dictionary with a range of T2 values. A pixel-wise fit of the source images with this dictionary then resulted in the map. The accuracy of the technique was established through comparison to spin-echo T2 mapping in the ISMRM-NIST phantom (QalibreMD). The lowest considered T2 values in the phantom was 15ms, and the highest was 430ms. The average relaxation times in each vial were compared to the reference by linear regression and Bland-Altman analysis.
Next, the technique was applied to the heart of 4 healthy volunteers (age=29±7y, 2 female) with IRB approval and written informed consent, and a Student’s t-test was used to compare the average values in the corresponding myocardial area to routine breath-held 2D T2 mapping.8 Finally, we applied this technique in one patient (age=87y, female) with a right coronary artery (RCA) myocardial infarction.

Results

The phantom T2 maps demonstrated good agreement with the gold standard over the relevant T2 range (y=0.98x+3.28, R2=0.99, Fig.2A-B). The fNAV correction resulted in well-aligned source images with no visible motion artifacts. The CS-based reconstruction provided sharp and denoised source images in the healthy volunteers (Fig.3A-E), which translated into sharp and precise T2 maps (Fig.3F-H). The 3D myocardial T2 values in vivo were in agreement with those of routine 2D maps (T2=49.5±4.6ms vs. 46.6±2.8ms, respectively, p=0.57), but had a lower precision: the average standard deviation in the segmented myocardium was 7.8ms for the proposed sequence, and 4.1ms for routine 2D maps. In the patient, our technique matched the T2 values obtained with routine breath-held 2D T2 mapping, despite the small T2 elevation in the infarcted region (Fig.4).

Discussion

The proposed technique showed good correlation to the reference T2 relaxation times in the NIST phantom, as well as a low bias. The in-vivo T2 values agreed well with those of the routine 2D technique, both in the healthy volunteers and the patient. The lower precision that was found is mostly likely caused by the 11-fold smaller pixel size.
In conclusion, we demonstrated the high accuracy of a novel fast high-resolution 3D T2 mapping technique for the heart. The benefit of the motion robustness of the technique due to the fNAV and radial acquisition needs to be demonstrated in a large patient study. Furthermore, the image and map reconstruction can be further optimized to increase the mapping precision.

Acknowledgements

No acknowledgement found.

References

1. Messroghli, D. R. et al. Clinical recommendations for cardiovascular magnetic resonance mapping of T1, T2, T2* and extracellular volume: A consensus statement by the Society for Cardiovascular Magnetic Resonance (SCMR) endorsed by the European Association for Cardiovascular Imaging (EACVI). J Cardiovasc Magn Reson 19, 75 (2017).

2. van Heeswijk, R. B. et al. Self-navigated isotropic three-dimensional cardiac T2 mapping. Magn. Reson. Med. 73, 1549–54 (2015).

3. Ding, H. et al. Three-dimensional whole-heart T2 mapping at 3T. Magn. Reson. Med. 74, 803–16 (2015).

4. Milotta, G. et al. 3D Whole-heart free-breathing qBOOST-T2 mapping. Magn. Reson. Med. 83, 1673–1687 (2020).

5. Bustin, A. et al. Accelerated free‐breathing whole‐heart 3D T 2 mapping with high isotropic resolution. Magn. Reson. Med. 83, 988–1002 (2020).

6. Piccini, D., Littmann, A., Nielles-Vallespin, S. & Zenge, M. O. Spiral phyllotaxis: the natural way to construct a 3D radial trajectory in MRI. Magn. Reson. Med. 66, 1049–1056 (2011).

7. Roy, C. W. et al. Motion Compensated Whole-Heart Coronary Magnetic Resonance Angiography using Focused Navigation (fNAV). arXiv:2010.14206 (2020).

8. Giri, S. et al. T2 quantification for improved detection of myocardial edema. J. Cardiovasc. Magn. Reson. 11, 56 (2009). 9. Hamilton, J. I. et al. MR Fingerprinting for Quantification of Myocardial T1, T2, and M0. in (2015).

Figures

Schematic overview of the proposed T2 mapping technique. Five 3D bSSFP acquisitions with increasing T2prep times were acquired. The respiratory motion was then 3D translationally corrected through focused navigation (fNAV). The resulting k-spaces were transformed into source images through compressed sensing with LLR and TV regularization. Finally, the T2 map was calculated by fitting a dictionary that is generated from the timings of the individual raw data.

Accurate T2 mapping in the NIST-ISMRM phantom. A) Linear regression plot of the T2 values obtained with the proposed technique. The T2 relaxation times match the spin-echo reference over the entire range of interest (15ms-430ms). B) The Bland-Altman analysis demonstrates a negligible bias. The dotted lines represent the 95% limits of agreement.

Sharp source images and precise T2 maps in healthy volunteers. A-E) The source images contain sharp features (papillary muscles, edge definition) and are well aligned to one another. F-H) The resulting short-axis T2 map and its two perpendicular views are similarly sharp.

T2 mapping in a patient with myocardial infarction. A) LGE can be observed in the septum. B) The corresponding T2 map obtained with 2D T2-prepared bSSFP shows a small T2 elevation in the infarcted area. C) The proposed technique accurately confirms the T2 elevation.

Proc. Intl. Soc. Mag. Reson. Med. 29 (2021)
0692