3632

Deep learning for undersampled spiral DENSE reconstruction
Samuel W Fielden1,2, Eric D Carruth1, Christopher D Nevius1, Christopher M Haggerty1, and Brandon K Fornwalt1,3
1Imaging Sciences & Innovation, Geisinger, Danville, PA, United States, 2Medical & Health Physics, Geisinger, Danville, PA, United States, 3Radiology, Geisinger, Danville, PA, United States

Synopsis

Displacement Encoding with Stimulated Echoes (DENSE) is a powerful technique that has found great utility in accurately measuring cardiac tissue displacement. However, DENSE remains time-consuming to acquire, particularly for 3-dimensionally encoded or higher resolution schemes, and so methods to accelerate image acquisition are needed. Here, we apply the Deep Cascade of Convolutional Neural Networks (DCCNN) to the complex-valued, non-Cartesian data of DENSE to show that accelerated imaging via k-space undersampling is feasible using a deep learning-based reconstruction.

Introduction

Cine Displacement Encoding with Stimulated Echoes (DENSE) is a powerful but time-consuming technique used to obtain accurate and reproducible tissue displacement measurements1. Efforts to speed up DENSE have included an early switch to spiral trajectories, as well as parallel-imaging and compressed-sensing based approaches to reconstruction of undersampled data2. To reduce the number of phase cycles needed to suppress artifact-generating echoes, deep learning (DL) has been applied to DENSE3, but to date, these approaches have not been applied for undersampled image reconstruction. The Deep Cascade of Convolutional Neural Networks (DCCNN) technique makes use of concatenated CNNs, with maintenance of data fidelity through regular updates with the acquired k-space data4. It is extensible to non-Cartesian data in a straightforward manner and has the further advantage of maintaining the complex nature of the MRI data throughout the reconstruction pipeline, as well as handling multiple coils/contrasts. As DENSE is a spiral phase-contrast method, we sought to extend the DCCNN method to simulated undersampled data to explore the feasibility of accelerated DENSE imaging.

Methods

DENSE data from 5 volunteers obtained previously5 was used to generate training, validation, and test image sets. Between 10 and 14 slices from each volunteer were available, encoded in 3 dimensions over 31-39 cardiac phases. Two randomly selected slices (all cardiac phases) from each volunteer were reserved to serve as the validation and test data, respectively, resulting in total dataset sizes of 1678 images used for training, 171 for validation, and 171 for testing.

The 2D DCCNN method was extended for non-Cartesian trajectories by replacing the Fourier transform operator at all applicable points with the nuFFT6 while maintaining appropriate scaling and mask operations (Fig. 1). The number of convolution layers (nd) in each network as well as the depth of the network (nc) were each set to 5. Pseudo raw data was generated by inverse gridding complex image data using a 12-interleaf, 6000-sample dual-density spiral trajectory. Undersampled data was generated by deleting every other interleaf. Trajectories were designed such that the center 33% of k-space remained fully-sampled after undersampling. X-, Y-, and Z- encoded images were treated as separate channels, and the resulting images were analyzed for general similarity to the reference using PSNR. Magnitude and phase SNR of the myocardium was recorded, and differences assessed by t-test. DENSEanalysis7 was used to interrogate the magnitude and variance of in-plane myocardial strain measurements generated by the reconstruction in comparison to the fully-sampled data.

Results

The myocardial magnitude SNR fell by 62% when undersampled, with nearly full recovery after reconstruction. Phase SNR was less affected; it was reduced by 10% when undersampled, but also recovered by the reconstruction (Fig. 2). Magnitude images from the sixth cardiac phase of each subject are shown in Figure 3. Fully-sampled, undersampled, and reconstructed images are shown left-to-right for each subject. The denoising performance varied with subject, with Subjects 1-3 experiencing a higher amount of denoising than 4 and 5, as measured by the PSNR of each image, compared to the fully-sampled reference image (Fig. 4). No differences in measured circumferential or radial strains were observed.

Discussion

The reconstruction improved magnitude image appearance in every case; however, the relative gain varied across subjects. The relative improvement for Subject 4, in particular, was lower than in the others. It is notable that the orientation of the heart for Subject 4 led to a more axially-oriented slice than the others. This may highlight the importance of training on a wide variety of body types.

The dual-density trajectories used for this study provide fully-sampled, low-resolution image data, contaminated with high frequency artifacts and noise. As the reconstruction can be viewed as a denoising problem, this type of trajectory proved sufficient to successfully reconstruct complex DENSE images from noisy undersampled data. The reconstruction failed when presented with undersampled constant-density spiral data (not shown) due to the resulting highly coherent artifacts. No difference in measured strains was observed at the current conservative acceleration factor. Future work might include optimizing the trajectory to minimize structured unsampling artifacts as higher acceleration rates are investigated.

Conclusion

Here, we have shown a DL approach with promise to reconstruct undersampled complex-valued spiral data, such as that used in DENSE. More investigation as to the type and number of anatomies needed for training, as well as the optimal pairing of trajectory design and DL algorithm, is needed.

Acknowledgements

No acknowledgement found.

References

1. Kim D, Gilson WD, Kramer CM, Epstein FH. Myocardial tissue tracking with two-dimensional cine displacement-encoded MR imaging: development and initial evaluation. Radiology. 2004 Mar;230(3):862-71.

2. Chen X, Yang Y, Cai X, Auger DA, Meyer CH, Salerno M, Epstein FH. Accelerated two-dimensional cine DENSE cardiovascular magnetic resonance using compressed sensing and parallel imaging. J Magn Reson Imaging. 2017 Sep;46(3):887-96.

3. Abdishektaei M, Feng X, Meyer CH, Epstein FH. DAS-Net: A generative adversarial net to suppress artifact-generating echoes in DENSE MRI. In ISMRM, Montreal, Canada, 2019, p1089.

4. Schlemper J, Caballero J, Hajnal JV, Price A, Rueckert D. A deep cascade of convolutional neural networks for MR image reconstruction. IEEE Trans Med Imaging. 2017;37(2):491-503.

5. Suever JD, Wehner GJ, Jing L, Powell DK, Hamlet SM, Grabau JD, Mojsejenko D, Andres KN, Haggerty CM, Fornwalt BK. Right ventricular strain, torsion, and dyssynchony in healthy subjects using 3D spiral cine DENSE magnetic resonance imaging. IEEE Trans Med Imaging. 2017 May;36(5):1076-85.

6. Lin J-M. Python non-uniform fast Fourier transform (PyNUFFT): An accelerated non-Cartesian MRI package on a heterogeneous platform (CPU/GPU)." J Imaging 2018;4(3):51.

7. Spottiswoode BS, Zhong X, Hess AT, Kramer CM, Meintjes EM, Mayosi BM, Epstein FH. Tracking myocardial motion from cine DENSE images using spatiotemporal phase unwrapping and temporal fitting. IEEE Trans Med imaging. 2007 Jan;26(1):15-30. Code available at: https://github.com/denseanalysis/denseanalysis

Figures

A deep cascade of convolutional neural networks (DCCNN) for 2D non-Cartesian undersampled reconstruction. x represents data in the image domain, y represents acquired k-space samples, tk is the k-space trajectory used with the nuFFT operator (Fnu) to enforce data consistency (DC) with the acquired spiral k-space samples through the function, fdc. The number of layers within each CNN is controlled through the variable nd, while the depth of the overall network is set by nc. Adapted from Schlemper, et al4.

Myocardial SNR normalized to fully-sampled reference, averaged across subjects. Magnitude SNR was significantly impacted by undersampling and recovered after reconstructing with the DCCNN. Phase SNR was less affected but still showed a small decrease and recovery. All differences were significant (p < 0.001).

Test slice from each of 5 volunteers, scaled such that myocardial signal is approximately equal in each image.

Mean magnitude PSNR for each subject across cardiac phase. The PSNR of all images for all subjects was significantly improved after reconstruction (p < 0.001).

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
3632