Coronary MRA performed using respiratory self-navigation techniques that apply retrospective motion-correction suffer from artifacts originating from static anatomical structures surrounding the heart. This becomes even more complicated at 3T due to insufficient fat suppression. Here we implemented free-breathing self-navigated 3D coronary MRA at 3T. Instead of using a motion model, we resolved the motion using 4D k-t sparse SENSE. Resulting images were then quantitatively compared with those reconstructed from the same datasets but using a 1D motion correction. We demonstrate that non-contrast whole-heart coronary MRA can be performed at 3T and that 4D motion-resolved reconstruction effectively minimizes adverse effects of respiratory-motion.
Non-contrast enhanced coronary MRA was performed in (n=7) healthy
adult volunteers on a 3T clinical scanner (Magnetom PRISMA, Siemens
Healthcare). Data were acquired using a prototype ECG-triggered respiratory self-navigated
free-breathing 3D radial GRE imaging sequence preceded by an adiabatic T2
preparation module (T2-Prep) to improve blood tissue contrast [6]. Off-resonant
water excitation (WE) as described previously [5] was used as a lipid-nulling strategy
to enhance contrast at the level of the coronary arteries and to minimize image
artifacts originating from background tissue. 3D volumes with an isotropic voxel
size of 1.1 mm3 were acquired with the following imaging
parameters: field of view (220 mm)3, matrix size 1923, TET2-Prep
= 40 ms, RF excitation angle 18°, TE/TR = 2.5 ms/5.6 ms, with 24 radial readouts
per segment for a total of ~15k k-space lines. All datasets were 1D respiratory
motion-corrected using a superior-inferior (SI) projection acquired every 24
k-space lines as previously described [2]. The same data also underwent motion-resolved
4D (x-y-z-respiratory) k-t
sparse SENSE reconstruction [4]. For this, respiratory
signals were extracted from the filtered k-space center modulations and used to
sort all k-space readouts in 4 different respiratory bins (Fig. 1), ranging from end-expiration to end-inspiration. Respiratory motion-resolved
images were then obtained by solving:
$$\underset{m}{\arg\min}\parallel F\cdot C \cdot m - s \parallel ^2 _2 + \lambda _1 \parallel D_1m \parallel _1 \quad , $$
where F represents the non-uniform fast Fourier transform (NUFFT)
operator, C the coil sensitivity
maps, m the 4D image data set to
reconstruct (where the 4th dimension is the respiratory dimension), s the radial k-space data, D1 the finite difference
operators applied along the respiratory dimension, and λ1 = 0.05 a regularization parameter that was empirically
determined. Both the motion-corrected and motion-resolved
images were qualitatively compared by visual inspection. For quantitative
comparison, vessel
sharpness of the first 4 cm of both the right coronary artery (RCA) and left
anterior descending (LAD), and maximum visible length of both coronaries were
calculated using SoapBubble [7]. Statistics were computed via two-tailed
student’s t-test for paired data.
[1] Stehning et al., Magn Reson Med 54, 476-480 (2005)
[2] Piccini et al., Magn Reson Med 68, 571-579 (2012)
[3] Feng et al., Magn Reson Med, doi: 10.1002/mrm.25665 (2015)
[4] Piccini et al., Magn Reson Med, In press (2016)
[5] Bastiaansen et al., Proc Intl Soc Mag Reson Med 24, 1072 (2016)
[6] Nezafat et al., Magn Reson Med 55, 858-864 (2006)
[7] Etienne et al., Magn Reson Med 48, 658-666 (2002)