Highly Efficient Nonrigid Motion Corrected 3D Whole-Heart Coronary Vessel Lumen and Wall Imaging
Gastao Cruz1, David Atkinson2, Markus Henningsson1, René Botnar1, and Claudia Prieto1

1Division of Imaging Sciences & Biomedical Engineering, King's College London, London, United Kingdom, 2University College London, London, United Kingdom

Synopsis

Non-invasive visualization of both coronary lumen and vessel wall is desired for assessment of coronary atherosclerosis. An interleaved acquisition was recently proposed to obtain both 3D images with MRI. However, this approach is susceptible to motion artifacts and dual respiratory gating results in long and unpredictable scan times. Here, we propose a ~100% scan efficiency, two-step motion correction method using translational and nonrigid correction to produce co-registered coronary lumen and vessel wall images. The proposed method shows significant improvements over translational correction and similar lumen quality to a reference navigator-gated (6mm) scan, despite a scan time reduction of ~1.8x.

INTRODUCTION

Coronary atherosclerosis is not necessarily stenotic due to outward remodelling of the vessel wall. Coronary plaque burden has been shown to correlate with risk of future coronary events1 and thus direct, non-invasive visualization of both coronary lumen and vessel wall is desired. A 3D flow independent approach for vessel wall imaging was proposed recently2, based on an interleaved acquisition and subtraction of data with (T2prep(+), lumen) and without (T2prep(-)) T2-preparation prepulse. This approach requires dual gating of both datasets to compensate for respiratory motion, resulting in long and unpredictable scan times. Here, we propose a novel nonrigid motion correction framework for free-breathing 3D whole-heart imaging using an interleaved scanning framework to produce co-registered coronary lumen and vessel wall images at near 100% scan efficiency. The proposed approach was compared with a 2D translational correction (TC) and no motion correction (NMC) in 9 subjects. Additionally, the proposed method was compared against a navigator gated (6mm) and tracked coronary lumen scan.

METHODS

Data was acquired using an interleaved scanning framework3 with three scans acquired in an alternated fashion: a set of 3D segmented whole-heart with and without T2 preparation and a 2D single shot golden-radial (GR) coronal image navigator (2D iNAV) (Fig. 1a). Motion correction was performed in two steps: beat-to-beat translational correction and bin-to-bin nonrigid correction (Fig. 1b). iNAVs were reconstructed and rigidly registered to estimate superior-inferior (SI) and right-left (RL) translational motion, allowing beat-to-beat translational correction. A respiratory signal was derived from the SI component and data was binned according to respiratory position. Bins were reconstructed with soft-gated4 iterative SENSE5 and registered to estimate 3D nonrigid motion. These motion fields were used in a General Matrix Description6 (GMD) reconstruction by solving: $$$\widehat{I} = arg min_I\left\{||\sum_b \bf A_b \bf F \bf S \bf U_b \bf I - \bf K||_2^2\right\}$$$, where $$$\widehat{I}$$$ is the motion corrected volume, $$$K$$$ the translational corrected k-space, $$$A_b$$$ the sampling matrix for bin b, $$$F$$$ the Fourier transform, $$$S$$$ the coil sensitivities and $$$U_b$$$ the nonrigid motion fields. Outliers due to deep breaths were removed prior to reconstruction. Both T2prep(+) and T2prep(-) were reconstructed with the framework described and nonrigidly registered to guarantee alignment. Vessel wall images were obtained via image subtraction as described by Andia et al2.

EXPERIMENTS

Nine healthy subjects were scanned free-breathing on a 1.5T Philips scanner using a 32-channel coil. T2prep(+) and T2prep(-) data were acquired with an ECG-triggered 3D Cartesian b-SSFP (1x1 mm in-plane resolution, 2 mm slice thickness, 300x300x90 mm FOV, TR/TE = 5.3/2.6 ms, flip angle = 70°) interleaved with a 2D golden radial iNAV spoiled gradient echo sequence (4x4 mm in-plane resolution, 25mm slice thickness, 300x300 FOV, TR/TE = 2.4/1.07 ms, flip angle = 5°). Non-motion correction (NMC), 2D translational correction (TC) and the proposed translation plus nonrigid (TC+GMD) reconstructions were obtained from these data. Additionally a 6mm diaphragmatic respiratory gated and tracked lumen scan (T2prep(+)) was performed for comparison. Coronary lumen was evaluated with vessel length and vessel sharpness (normalized to Gated), coronary vessel wall was evaluated by vessel wall thickness and vessel wall sharpness (normalized to TC+GMD); both images were evaluated by two experts on a scale of 0 (extreme blurring) to 4 (no blurring).

RESULTS

Reformatted lumen images for Gated, TC+GMD, TC and NMC for 2 subjects are shown in Fig. 2. Motion artifacts in NMC are reduced with TC and further reduced with TC+GMD. The corresponding set of vessel wall images is shown in Fig. 3, along with cross sectional views. Vessel wall appears blurred in NMC; delineation is significantly improved by TC and sharpness further improved by TC+GMD. Fig. 4 shows metrics for coronary lumen evaluation. Lumen vessel length sharpness is reduced for NMC and increases with TC, TC+GMD and Gated. Similar results were obtained from visual scores. Analogous results were found for the vessel wall in Fig. 5: blurring increases the measured wall thickness for NMC and is reduced with TC and TC+GMD; correspondingly, wall sharpness is significantly improved with TC and TC+GMD from NMC. Visual score of the vessel wall was in agreement with the remaining metrics. TC+GMD and Gated show similar image quality and non-significant differences in any metric studied. The proposed TC+GMD improved scan time by 1.8x (on average) compared to Gated.

CONCLUSION

A novel framework for simultaneous nonrigid correction for 3D coronary lumen and vessel wall was introduced. The proposed approach allows near 100% scan efficiency, while maintaining similar image quality to a 6mm gated and tracked lumen acquisition. Significant improvements were found with the proposed approach compared to translational correction alone.

Acknowledgements

This work was supported by the Medical Research Council (MRC), grant MR/L009676/1.

References

1. Kubo T, Imanishi T, Takarada S, Kuroi A, Ueno S, Yamano T, Tanimoto T, Matsuo Y, Masho T, Kitabata H, Tsuda K, Tomobuchi Y, Akasaka T. Assessment of culprit lesion morphology in acute myocardial infarction: ability of optical coherence tomography compared with intravascular ultrasound and coronary angioscopy. J Am Coll Cardiol 2007;50:933–939.

2. Andia ME, Henningsson M, Hussain T, Phinikaridou A, Protti A, Greil G, Botnar RM. Flow-independent 3D whole-heart vessel wall imaging using an interleaved T2-preparation acquisition. Magn Reson Med 2013; 69:150–157.

3. Henningsson M, Mens G, Koken P, Smink J, Botnar RM. A new framework for interleaved scanning in cardiovascular MR: application to image-based respiratory motion correction in coronary MR angiography. Magn Reson Med 2014. doi: 10.1002/mrm.25149.

4. Johnson KM, Block WF, Reeder SB, Samsonov A. Improved least squares MR image reconstruction using estimates of k-space data consistency. Magn Reson Med 2012;67:1600–1608.

5. Pruessmann KP, Weiger M, Börnert P, Boesiger P. Advances in sensitivity encoding with arbitrary k-space trajectories. Magn Reson Med 2001;46:638–651.

6. Batchelor PG, Atkinson D, Irarrazaval P, Hill DLG, Hajnal J, Larkman D. Matrix description of general motion correction applied to multishot images. Magn Reson Med 2005;54:1273–1280.

Figures

Fig1. Proposed approach. a) Acquisition: Data is acquired using interleaved scanning, allowing for datasets with/without T2 preparation (T2prep(+)/T2prep(-), respectively) to be acquired simultaneously with a 2D image navigator (2D iNAV). b) Reconstruction: 2D iNAV-based beat-to-beat translational motion correction followed by 3D bin-to-bin nonrigid motion correction incorporated in the reconstruction (GMD).

Fig2. Reformatted coronary lumen images (T2prep(+)) for Gated, proposed nonrigid approach (TC+GMD), translation correction (TC) and no motion correction (NMC) for two subjects. Blurring in NMC images is reduced with TC and sharpness further increased with TC+GMD (magnified boxes). TC and TC+GMD have similar image quality to the Gated.

Fig3. Vessel wall images for proposed nonrigid correction (TC+GMD), translation correction (TC) and no motion correction (NMC) for two subjects. Cross vessel planes (dotted markings) are shown in magnified boxes. The vessel wall is obscured in NMC. A significant improvement is obtained with TC and sharpness further improved with TC+GMD.

Fig 4. Image metrics for coronary lumen for Gated, TC+GMD, TC and NMC. Statistically significant differences on a P-value > 0.01 are marked with (*). Vessel length for LCA (a) and RCA (b) for 9 subjects (coloured). Vessel sharpness for LCA (c) and RCA (d). Visual score (e).

Fig 5. Image metrics for coronary wall for Gated, TC+GMD, TC and NMC. Statistically significant differences on a P-value > 0.01 are marked with (*).Vessel wall thickness for LCA (a) and RCA (b) for 9 subjects (coloured). Vessel wall sharpness for LCA (c) and RCA (d). Visual score (e).



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