2890

A fast navigator (fastNAV) for prospective respiratory motion correction in first-pass myocardial perfusion imaging
Ronald Mooiweer1, Radhouene Neji2, Sarah McElroy1, Muhummad Sohaib Nazir1, Reza Razavi1, Amedeo Chiribiri1, and Sébastien Roujol1
1School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom, 2MR Research Collaborations, Siemens Healthcare, Frimley, United Kingdom

Synopsis

A fast respiratory navigator (fastNAV) was developed for dynamic contrast enhanced CMR perfusion imaging by combining spatially non-selective saturation with slice-selective tip-up and slice-selective excitation pulses. A calibration scan was developed to enable the estimation of subject-specific tracking factors. Prospective motion correction using fastNAV was applied to perfusion imaging in 10 patients under free-breathing conditions. Compared to conventional perfusion imaging, fastNAV reduced the effect of respiratory motion while no difference in image quality was observed.

Introduction

First-pass cardiac magnetic resonance (CMR) perfusion imaging is widely used for patients with suspected coronary artery disease. While breath-holding acquisitions are not always effective, prospective respiratory motion correction using a navigator and slice tracking can reduce through-plane motion (1,2). CMR perfusion imaging is based on saturation-recovery, posing a challenge for obtaining a consistent navigator signal while keeping the temporal footprint to a minimum and avoiding a reduction in perfusion image quality (3–7). Furthermore, while a fixed tracking factor is often used to relate right hemi-diaphragm (RHD) navigator motion to cardiac motion (1,2), subject-specific tracking factors can improve motion correction even further (8,9). Here, a fast navigator (fastNAV) was develop for perfusion imaging and tested using subject-specific slice tracking (10).

Methods

Pulse sequence
The prototype navigator signal is obtained by combining spatially non-selective saturation (for perfusion imaging) with slice-selective tip-up and slice-selective excitation pulses (Figure 1A&B). The slice-selective tip-up pulse (-90˚flip-angle, 40mm thickness) was between a BIR4-90 adiabatic saturation pulse (11,12) and a spoiler gradient. The navigator excitation slice (15˚ flip-angle, 20mm thickness) was angulated from the tip-up slice in the transverse plane to generate signal from the overlap of these two pulses on the RHD. Gradient-recalled navigator signal was frequency encoded in the superior-inferior direction (Figure 1C). Additional time per navigator: 15ms, including 5ms for real-time feedback.

Subject-specific tracking factor estimation
FastNAV-based slice tracking in superioinferior direction was applied to each short-axis slice for perfusion imaging. Conventional cross-correlation analysis (13) was used to extract RHD displacement, which was multiplied by a subject-specific tracking factor. A separate sequence was used to calibrate the subject-specific tracking factor. To allow the near-simultaneous measurement of fastNAV signal and left ventricle (LV) displacement for calibration, one coronal image was acquired per heart-beat, followed by one module including saturation and fastNAV acquisition (Figure 1D). This was repeated 30x to capture several breathing cycles. LV motion was determined offline using a rigid motion estimation algorithm (14). The slope between the LV and fastNAV motion traces was found by least-squares-fitting, which was used as the subject-specific tracking factor, and the correlation coefficient (R) was calculated. Calibration was typically done in 1 minute, once the data was transferred off-line.

Experimental validation
10 Patients were included in this study, which was performed at 1.5T (MAGNETOM Aera, Siemens Healthcare, Erlangen, Germany). Two dynamic contrast enhanced rest perfusion protocols were performed in each patient, using the fastNAV-based motion correction and a conventional approach without motion correction, with a 10 minute interval in between them to allow for contrast washout. The order of the two perfusion protocols was randomised across subjects. Both perfusion acquisitions used an ECG-triggered balanced steady state free precession (bSSFP) saturation recovery sequence prescribed in the short axis orientation with the following parameters: TR/TE: 2.38/1.04ms, FA:50˚, FOV:360x270mm2, slices:3, slice thickness:8mm, matrix:192x144, voxel size:1.9x1.9mm2, Partial Fourier factor:3/4, GRAPPA:2 (24 integrated reference lines), bandwidth:1085 Hz/pixel, total acquisition duration per heartbeat: 603.9ms, saturation recovery time:133ms, dynamics:60. Patients were instructed to breathe normally.

In-plane motion of the LV was quantified by measuring the average Dice similarity coefficient (15) over dynamics (avDice) as well as the average displacement of the LV center of mass location (avCoM), following manual delineation in each dynamic. Image quality was assessed for each dataset: 4=excellent, 3=minor artefact but not limiting diagnosis, 2=major artefact but not limiting diagnosis, 1=poor image quality and non-diagnostic.

Results

Figure 2 shows an example of tracking factor calibration. Over all patients, subject-specific tracking factors of 0.46±0.13 (min=0.18, max=0.68) were measured with an R value of 0.94±0.03 (min=0.90, max=0.98), indicative of a robust relationship between measured fastNAV and LV motion. A robust fastNAV signal was acquired during perfusion imaging in all patients (Figure 3). An example of motion-corrected perfusion imaging is shown in Figure 4.

Images with fastNAV display improved motion metrics in most patient (Figure 5), while overall higher avDice similarity (0.94±0.02 vs. 0.91±0.03, p<0.001) and reduced avCoM displacement (4.03±0.84 vs. 5.22±1.22, p<0.001) were reported. Qualitatively, all slices were scored either as “excellent” or “minor artefact but not limiting diagnosis” with no statistically significant difference overall.

Discussion

FastNAV added only 15 ms to each saturation recovery block, shorter than previously reported navigators for CMR perfusion imaging (3,16,17). The average subject-specific tracking factor across subjects was 0.46, whereas the commonly used factor for conventional RHD NAV (18) is 0.6. The discrepancy might be due to the small sample size in this study.

Overall fastNAV-enabled prospective slice tracking resulted in substantial reduction of the global respiratory motion of the heart, but the presence of respiratory-induced local non-rigid deformation as observed in the apical slice of a few patients could not be corrected. The combination of prospective and retrospective motion correction is promising for 3D motion correction.

Future larger studies will be needed to characterise the clinical benefit of fastNAV under stress conditions in patients with coronary artery diseases. The validity of the tracking factor estimated at rest and used during stress conditions remains to be demonstrated.

Conclusion

The fastNAV enabled fast and robust RHD motion tracking in a CMR perfusion sequence. Combined with real-time slice-tracking and subject-specific slice-tracking, fastNAV reduced the effect of respiratory motion.

Acknowledgements

This work was supported by the Engineering and Physical Sciences Research Council (EPSRC) grant (EP/R010935/1), the British Heart Foundation (BHF) (PG/19/11/34243), the Wellcome EPSRC Centre for Medical Engineering at King’s College London (WT 203148/Z/16/Z), the National Institute for Health Research (NIHR) Biomedical Research Centre based at Guy’s and St Thomas’ National Health Service (NHS) Foundation Trust and King’s College London, and Siemens Healthineers. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health.

References

1. Danias PG, McConnell M V., Khasgiwala VC, Chuang ML, Edelman RR, Manning WJ. Prospective navigator correction of image position for coronary MR angiography. Radiology 1997;203:733–736 doi: 10.1148/radiology.203.3.9169696.

2. Wang Y, Riederer SJ, Ehman RL. Respiratory Motion of the Heart: Kinematics and the Implications for the Spatial Resolution in Coronary Imaging. Magn. Reson. Med. 1995;33:713–719 doi: 10.1002/mrm.1910330517.

3. Pedersen H, Kelle S, Ringgaard S, et al. Quantification of myocardial perfusion using free-breathing MRI and prospective slice tracking. Magn. Reson. Med. 2009;61:734–738 doi: 10.1002/mrm.21880.

4. Pedersen H, Ringgaard S, Kim WY. A patient-adapted navigator tracking approach for prospective 3D respiratory motion correction in free-breathing myocardial perfusion imaging. In: International Society for Magnetic Resonance in Medicine (ISMRM). ; 2005. p. 512.

5. Bush M, Chen C, Liu Y, Ahmad R, Jin N, Simonetti OP. Cardiac perfusion imaging with prospective respiratory motion correction. In: Proceedings from the 22nd Annual SCMR Scientific Sessions. Belleveu, WA; 2019. pp. 497–499.

6. Basha TA, Roujol S, Kissinger K V., et al. Free-breathing cardiac MR stress perfusion with real-time slice tracking. Magn. Reson. Med. 2014;72:689–698 doi: 10.1002/mrm.24977.

7. Spuentrup E, Manning WJ, Botnar RM, Kissinger K V., Stuber M. Impact of navigator timing on free-breathing submillimeter 3D coronary magnetic resonance angiography. Magn. Reson. Med. 2002;47:196–201 doi: 10.1002/mrm.10026.

8. Moghari MH, Hu P, Kissinger K V., et al. Subject-specific estimation of respiratory navigator tracking factor for free-breathing cardiovascular MR. Magn. Reson. Med. 2012;67:1665–1672 doi: 10.1002/mrm.23158.

9. Bush MA, Ahmad R, Jin N, Liu Y, Simonetti OP. Patient specific prospective respiratory motion correction for efficient, free-breathing cardiovascular MRI. Magn. Reson. Med. 2019;81:3662–3674 doi: 10.1002/mrm.27681.

10. Mooiweer R, Neji R, Sohaib M, Reza N. A fast navigator ( fastNAV ) for prospective respiratory motion correction in first-pass myocardial perfusion imaging. 2020:1–11 doi: 10.1002/mrm.28617.

11. Staewen R, Johnson A, Ross B, Parrish T, Merkle H, Garwood M. 3-D FLASH Imaging Using a Single Surface Coil and a New Adiabatic Pulse, BIR-4. Investig. Radiol. Radiol. 1990;25:559–567 doi: 10.1097/00004424-199005000-00015.

12. Kim D, Cernicanu A, Axel L. B0 and B1-insensitive uniform T1-weighting for quantitative, first-pass myocardial perfusion magnetic resonance imaging. Magn. Reson. Med. 2005;54:1423–1429 doi: 10.1002/mrm.20704.

13. Ehman RL, Felmlee JP. Adaptive technique for high-definition MR imaging of moving structures. Radiology 1989;173:255–263 doi: 10.1148/radiology.173.1.2781017.

14. Roujol S, Benois-Pineau J, De Senneville BD, Ries M, Quesson B, Moonen CTW. Robust real-time-constrained estimation of respiratory motion for interventional MRI on mobile organs. IEEE Trans. Inf. Technol. Biomed. 2012;16:365–374 doi: 10.1109/TITB.2012.2190366.

15. Dice LR. Measures of the Amount of Ecologic Association Between Species. Ecology 1945;26:297–302 doi: 10.2307/1932409.

16. Akçakaya M, Gulaka P, Basha T a, Ngo LH, Manning WJ, Nezafat R. Free-breathing phase contrast MRI with near 100% respiratory navigator efficiency using k-space-dependent respiratory gating. Magn. Reson. Med. 2013;00:1–8 doi: 10.1002/mrm.24874.

17. Cleppien DEJ, Horstick G, Abegunewardene N, et al. Comparison of the quantitative first pass myocardial perfusion MRI with and without prospective slice tracking: Comparison between breath-hold and free-breathing condition. Magn. Reson. Med. 2010;64:1461–1470 doi: 10.1002/mrm.22513.

18. Kramer CM, Barkhausen J, Bucciarelli-Ducci C, Flamm SD, Kim RJ, Nagel E. Standardized cardiovascular magnetic resonance imaging (CMR) protocols: 2020 update. J. Cardiovasc. Magn. Reson. 2020;22:1–18 doi: 10.1186/s12968-020-00607-1.

Figures

Figure 1. A) Diagram of cardiac triggered perfusion imaging with fastNAV. prep: preparatory pulses for saturation recovery imaging including fastNAV. B) Geometry of slice-selective fastNAV pulses overlaid on an axial localiser image. Orange rectangle: fastNAV tip-up slice, green rectangle: fastNAV excitaiton slice. C) Coronal localiser image showing fastNAV’s frequency encoding direction as blue arrow. D) Pulse sequence diagram of the fastNAV calibration scan, where fastNAV is acquired after a coronal image of the LV per heartbeat.

Figure 2. Tracking factor calibration for subject #3. A&B) Coronal images in inspiration (A) and expiration (B) that are used to detect LV motion. A dotted white line marks the inferior edge of the LV in inspiration to highlight the difference with expiration. C) fastNAV signal during calibration, white lines mark the detected RHD position. D) Superior-inferior fastNAV-based RHD motion and imaged-based LV motion. E) Correlation between the fastNAV-based and image-based motion traces. The slope was found using least-squares fitting and was used as subject-specific tracking factor.

Figure 3. Signal examples of fastNAV acquired during perfusion imaging with fastNAV-based motion correction for all subjects. Consecutive fastNAV signals acquired over the first 19 heartbeats of each dynamic acquisition (3 fastNAV signals per heartbeat, i.e one per slice) are shown for each subject. The fastNAV signal had high SNR, was consistent over time. The visual displacements observed in the fastNAV signal traces were successfully tracked (white horizontal marks).

Figure 4. Video of the entire dynamic perfusion imaging series acquired with the conventional non-motion corrected sequence (top row) and with the fastNAV-enabled prospective motion corrected sequence (bottom row) for subject #3. Base, mid, and apical images are displayed at 10 frames per second, every second frame is shown, and image quality is reduced for file size reduction. Clear vertical movement is seen in the standard sequence of images, while this is largely removed when fastNAV is applied.

Figure 5. Quantitative motion analysis of perfusion images over all subjects using the conventional sequence (“no MoCo”, red circles) and the fastNAV prospective motion corrected sequence (“fastNAV MoCo”, blue squares). A) Average overlap of the left ventricle (LV) across images (avDice). B) Average displacement of the LV center of mass (avCoM), with inverted y-axis to aid the comparison with avDice. All results are shown per subject and per slice (b/m/a: base/mid/apex).

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