Breath-Hold Real-Time Phase Contrast MRI using Radial k-space Sampling and Compressed Sensing
Hassan Haji-Valizadeh1, Elwin Bassett2, Genesh Adluru3, Edward VR DiBella 3, and Daniel Kim3

1Bioengineering, University of Utah, Salt Lake City, UT, United States, 2Physics, University of Utah, Salt Lake City, UT, United States, 3Radiology,UCAIR, University of Utah, Salt Lake City, UT, United States

Synopsis

Phase contrast (PC) MRI is a useful tool for assessing hemodynamic, but suffers from low data acquisition efficiency. In this study we compared real-time PC MRI between Cartesian and Radial undersampling trajectories. Our results show that both real-time MRI pulse sequences yield velocity measurements that agree well with those produced by reference breath-hold PC MRI pulse sequence. Compared with real-time MRI with Cartesian sampling, Radial sampling produced images with fewer artifacts. This study demonstrates feasibility of real-time PC MRI using radial k-space sampling and constrained reconstruction.

Purpose

Phase contrast (PC) MRI is a useful tool for assessment of hemodynamic consequences caused by a variety of heart diseases. This study seeks to develop and test real-time PC MRI [1, 2] using compressed sensing (CS)[3, 4]. We chose to acquire breath-hold images in order to avoid any artifact due to motion. In this study, we compare two commonly used k-space sampling patterns (Radial and Cartesian). Both trajectories have desirable characteristics which may significantly affect PC reconstruction quality.Cartesian trajectories are less sensitive to gradient imperfection, and are capable of rectangular field of view that may be tailored for the human chest. Radial trajectories are more resistant to motion during acquisition [5], and have better SNR due to higher k-space center sampling. Furthermore, the temporal resolution of radial trajectories can be retrospectively determined when irrational angles such as golden angle ratio are used for acquisition [6]. To the best of our knowledge, no study has compared real time PC MRI between Cartesian and Radial k-space sampling schemes. We sought to compare their performance in human subjects.

Method

(Imaging Acquisition) Based on the work by Joseph et al. [2], we implemented Radial real-time PC MRI pulse sequences using 7 k-space lines per image (i.e., 14 k-space lines per cardiac phase since interleaved reference and velocity-encoded acquisitions [7]). For real-time with Radial sampling, we used golden angle ratio with tiny gold angles [8](Fig. 1). For real-time with Cartesian sampling, we used a sampling pattern with variable density, where two central k-space lines are always sampled [8](Fig. 9). For reference, we performed breath-hold acquisition with GRAPPA (scan time = 9 heart beats). We acquired volunteers data on both 3T (Prisma, Siemens), and 1.5T (Espree, Siemens). The imaging parameters The three acquisitions (reference, radial real-time, Cartesian real-time) used the imaging parameters, which included: FOV between 350 – 400 mm for Cartesian acquisitions with 75% FOV in phase encoding direction, FOV between 275-300 mm for Radial acquisitions, slice thickness = 8 mm, flip angle = 20, real-time scan time = 6 heart beats, matrix size of 128X128 for Radial acquisitions, matrix size of 128X96 for Cartesian acquisitions, and temporal resolution = 74.70 ms for Espree data sets and 72.80 ms for Prisma data sets. We imaged a total of five valve planes. Two aortic valves, and two mitral valves on the Espreee, and 1 aortic valve on the Prisma. Aortic valves were images using venc = 250 cm/s and the mitral valve plane with venc = 150 cm/s. Pulse sequence order (reference, radial real-time, and Cartesian real-time) was randomized.

(Image Reconstruction) Coil sensitivity profiles were obtained by Self-Calibration [10], and Undersampled data sets were reconstructed using CS with temporal Total Variation (TTV) constraint in combination with SENSE parallel imaging (40 iterations). TTV normalized regularization weights were determined empirically using training date sets: Cartesian, 0.03; radial, 0.095 for Espree data sets, and Cartesian 0.04; radial, 0.6 for Prisma data set. Region of interests (ROI) were drawn manually around each valve, and the max velocity was calculated. Pearson’s correlation and Bland Altman analyses were conducted on pooled multi-phase velocity measurements

Results

As shown in figure 2, compared with reference breath-hold acquisition, real-time PC MRI with Cartesian and Radial sampling schemes produced comparable image quality. Compared with radial, Cartesian sampling produced higher artifacts as shown. Compared with reference peak velocity measurements, both real-time PC MRI with Cartesian and radial sampling schemes produced measurements that are strongly correlated and in good agreement (see Fig. 3).

Discussion

This study demonstrates feasibility of real time PC MRI using either Cartesian and Radial k-space sampling pattern. Compared with Cartesian, radial produced less images artifacts. Future studies include assessment of diagnostic performance in patients with limited breath-hold capacity or in patients with arrhythmia (e.g., beat-to-beat variation).

Acknowledgements

This work was supported in part by the following grants:

NIH- 5R01HL116895-02

AHA - 14GRNT18350028

References

[1] Lin HY, Bender JA, Ding Y, Chung YC, Hinton AM, Pennell ML, Simonetti OP. Shared velocity encoding: A method to improve the temporal resolution of phase-contrast velocity measurements. Magn Reson Med 2012, 68(3): 703-710.

[2] Joseph AA, Kowallick JT, Merboldt KD, Voit D, Schaetz S, Zhang S, Frahm J. Real-time flow MRI of the aorta at a resolution of 40 msec. Magn Reson Med 2014; 40(1): 206-213.

[3] Lustig M, Donoho D, Pauly JM. Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med 2007;58(6):1182-1195.

[4] Kim D, Dyvorne HA, Otazo R, Feng L, Sodickson, DK, Lee VS. Accelerated phase-contrast cine MRI using k-t SPARSE-SENSE. Magn Reson Med 2012; 67(4): 1054-1064.

[5] Wild JM, Paley MN, Kasuboski L, Swift A, Fichele S, Woodhouse N, Griffiths PD, van Beek EJ. Dynamic radial projection MRI of inhaled hyperpolarized 3He gas. Magn Reson Med 2003;49(6):991-997.

[6] Winkelmann S, Schaeffter T, Koehler T, Eggers H, Doessel O. An optimal radial profile order based on the Golden Ratio for time-resolved MRI. IEEE transactions on medical imaging 2007;26(1):68-76.

[7] Axel L, Morton D. MR flow imaging by velocity-compensated/uncompensated difference images. J Comput Assist Tomogr 1987;11:31–34.

[8] Wundrak S, Paul J, Ulrici J, Hell E, Geibel MA, Bernhardt P, Rottbauer W, Rasche V. Golden ratio sparse MRI using tiny golden angles. Magn Reson Med 2015.

[9] Feng L, Otazo R, Jung H, Jensen JH, Ye JC, Sodickson DK, Kim D. Accelerated cardiac T2 mapping using breath-hold multiecho fast spin-echo pulse sequence with k-t FOCUSS. Magn Reson Med 2011; 65(6): 1661-1669.

[10] Walsh DO, Gmitro AF, Marcellin MW. Adaptive reconstruction of phased array MR imagery. Magn Reson Med 2000;43:682–690.

Figures

(Table.1) Bland Altman and Pearson’s correlation statistical analyses. Breath-hold PC MRI is used as reference. Both real-time PC MRI with Radial and Cartesian sampling trajectories produced velocity measurements that are strongly correlated and in good agreement with reference measurements

(Figure 1) Undersampled Radial and Cartesian k-space trajectories used for real time PC MRI, where each acquires 7 k-space lines per image. Radial used tiny golden angles to optimize CS performance.

(Figure 2) Comparison of real time Cartesian and Radial sequences to reference images. As shown by the blue arrows, Cartesian real time phase images suffer from aliasing along phase encoding direction which is not seen in real time Radial acquisition.

(Figure 3) Bland Altman and Regression analysis of maximum velocity for Real -time Cartesian and Radial Breath hold acquisitions. When compared to reference VPW, Real-time Radial PC derived VPW was characterized by a 95% confidence interval and correlation of comparable size and value as compared to real-time Cartesian.

(Figure 4) Comparison of maximum Velocity VPW for reference, real-time Radial and real-time Cartesian acquisition schemes. Both Real-time Radial and Cartesian derived VPW exhibited high correlation with respect to reference derived VPW.



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