4D Golden-Angle Radial MRI at Subsecond Temporal Resolution

Technical Note: Swing golden angle - A navigator-interleaved golden angle trajectory with eddy current suppression - Application in free-running cardiac MRI

BACKGROUND

Golden angle (GA) radial trajectory is advantageous for dynamic magnetic resonance imaging (MRI). Recently, several advanced algorithms have been developed based on navigator-interleaved GA trajectory to realize free-running cardiac MRI. However, navigator-interleaved GA trajectory suffers from the eddy-current effect, which reduces the image quality.

PURPOSE

This work aims to integrate the navigator-interleaved GA trajectory with clinical cardiac MRI acquisition, with the minimum eddy-current artifacts. The ultimate goal is to realize a high-quality free-running cardiac imaging technique.

METHODS

In this paper, we propose a new "swing golden angle" (swingGA) radial profile order. SwingGA samples the k-space by rotating back and forth at the generalized golden ratio interval, with smoothly interleaved navigator readouts. The sampling efficiency and angle increment distributions were investigated by numerical simulations. Static phantom imaging experiments were conducted to evaluate the eddy current effect, compared with cartesian, golden angle radial (GA), and tiny golden angle (tGA) trajectories. Furthermore, 12 heart-healthy subjects (aged 21-25 years) were recruited for free-running cardiac imaging with different sampling trajectories. Dynamic images were reconstructed by a low-rank subspace-constrained algorithm. The image quality was evaluated by signal-to-noise-ratio and spectrum analysis in the heart region, and compared with traditional clinical cardiac MRI images.

RESULTS

SwingGA pattern achieves the highest sampling efficiency (mSE > 0.925) and the minimum azimuthal angle increment (mAD < 1.05). SwingGA can effectively suppress eddy currents in static phantom images, with the lowest normalized root mean square error (nRMSE) values among radial trajectories. For the in-vivo cardiac images, swingGA enjoys the highest SNR both in the blood pool and myocardium, and contains the minimum level of high-frequency artifacts. The free-running cardiac images have good consistency with traditional clinical cardiac MRI, and the swingGA sampling pattern achieves the best image quality among all sampling patterns.

CONCLUSIONS

The proposed swingGA sampling pattern can effectively improve the sampling efficiency and reduce the eddy currents for the navigator-interleaved GA sequence. SwingGA is a promising sampling pattern for free-running cardiac MRI.

Movienet: Deep space-time-coil reconstruction network without k-space data consistency for fast motion-resolved 4D MRI

PURPOSE

To develop a novel deep learning approach for 4D-MRI reconstruction, named Movienet, which exploits space-time-coil correlations and motion preservation instead of k-space data consistency, to accelerate the acquisition of golden-angle radial data and enable subsecond reconstruction times in dynamic MRI.

METHODS

Movienet uses a U-net architecture with modified residual learning blocks that operate entirely in the image domain to remove aliasing artifacts and reconstruct an unaliased motion-resolved 4D image. Motion preservation is enforced by sorting the input image and reference for training in a linear motion order from expiration to inspiration. The input image was collected with a lower scan time than the reference XD-GRASP image used for training. Movienet is demonstrated for motion-resolved 4D MRI and motion-resistant 3D MRI of abdominal tumors on a therapeutic 1.5T MR-Linac (1.5-fold acquisition acceleration) and diagnostic 3T MRI scanners (2-fold and 2.25-fold acquisition acceleration for 4D and 3D, respectively). Image quality was evaluated quantitatively and qualitatively by expert clinical readers.

RESULTS

The reconstruction time of Movienet was 0.69 s (4 motion states) and 0.75 s (10 motion states), which is substantially lower than iterative XD-GRASP and unrolled reconstruction networks. Movienet enables faster acquisition than XD-GRASP with similar overall image quality and improved suppression of streaking artifacts.

CONCLUSION

Movienet accelerates data acquisition with respect to compressed sensing and reconstructs 4D images in less than 1 s, which would enable an efficient implementation of 4D MRI in a clinical setting for fast motion-resistant 3D anatomical imaging or motion-resolved 4D imaging.

Radial magnetic resonance image reconstruction with a deep unrolled projected fast iterative soft-thresholding network

Radially sampling of magnetic resonance imaging (MRI) is an effective way to accelerate the imaging. How to preserve the image details in reconstruction is always challenging. In this work, a deep unrolled neural network is designed to emulate the iterative sparse image reconstruction process of a projected fast soft-threshold algorithm (pFISTA). The proposed method, an unrolled pFISTA network for Deep Radial MRI (pFISTA-DR), include the preprocessing module to refine coil sensitivity maps and initial reconstructed image, the learnable convolution filters to extract image feature maps, and adaptive threshold to robustly remove image artifacts. Experimental results show that, among the compared methods, pFISTA-DR provides the best reconstruction and achieved the highest PSNR, the highest SSIM and the lowest reconstruction errors.

The CALIPR framework for highly accelerated myelin water imaging with improved precision and sensitivity

Quantitative magnetic resonance imaging (MRI) techniques are powerful tools for the study of human tissue, but, in practice, their utility has been limited by lengthy acquisition times. Here, we introduce the Constrained, Adaptive, Low-dimensional, Intrinsically Precise Reconstruction (CALIPR) framework in the context of myelin water imaging (MWI); a quantitative MRI technique generally regarded as the most rigorous approach for noninvasive, in vivo measurement of myelin content. The CALIPR framework exploits data redundancy to recover high-quality images from a small fraction of an imaging dataset, which allowed MWI to be acquired with a previously unattainable sequence (fully sampled acquisition 2 hours:57 min:20 s) in 7 min:26 s (4.2% of the dataset, acceleration factor 23.9). CALIPR quantitative metrics had excellent precision (myelin water fraction mean coefficient of variation 3.2% for the brain and 3.0% for the spinal cord) and markedly increased sensitivity to demyelinating disease pathology compared to a current, widely used technique. The CALIPR framework facilitates drastically improved MWI and could be similarly transformative for other quantitative MRI applications.

Fast Compressed Sensing of 3D Radial T1 Mapping with Different Sparse and Low-Rank Models

Knowledge of the relative performance of the well-known sparse and low-rank compressed sensing models with 3D radial quantitative magnetic resonance imaging acquisitions is limited. We use 3D radial T1 relaxation time mapping data to compare the total variation, low-rank, and Huber penalty function approaches to regularization to provide insights into the relative performance of these image reconstruction models. Simulation and ex vivo specimen data were used to determine the best compressed sensing model as measured by normalized root mean squared error and structural similarity index. The large-scale compressed sensing models were solved by combining a GPU implementation of a preconditioned primal-dual proximal splitting algorithm to provide high-quality T1 maps within a feasible computation time. The model combining spatial total variation and locally low-rank regularization yielded the best performance, followed closely by the model combining spatial and contrast dimension total variation. Computation times ranged from 2 to 113 min, with the low-rank approaches taking the most time. The differences between the compressed sensing models are not necessarily large, but the overall performance is heavily dependent on the imaged object.

Kz-accelerated variable-density stack-of-stars MRI

This work aimed to develop a modified stack-of-stars golden-angle radial sampling scheme with variable-density acceleration along the slice (kz) dimension (referred to as VD-stack-of-stars) and to test this new sampling trajectory with multi-coil compressed sensing reconstruction for rapid motion-robust 3D liver MRI. VD-stack-of-stars sampling implements additional variable-density undersampling along the kz dimension, so that slice resolution (or volumetric coverage) can be increased without prolonging scan time. The new sampling trajectory (with increased slice resolution) was compared with standard stack-of-stars sampling with fully sampled kz (with standard slice resolution) in both non-contrast-enhanced free-breathing liver MRI and dynamic contrast-enhanced MRI (DCE-MRI) of the liver in volunteers. For both sampling trajectories, respiratory motion was extracted from the acquired radial data, and images were reconstructed using motion-compensated (respiratory-resolved or respiratory-weighted) dynamic radial compressed sensing reconstruction techniques. Qualitative image quality assessment (visual assessment by experienced radiologists) and quantitative analysis (as a metric of image sharpness) were performed to compare images acquired using the new and standard stack-of-stars sampling trajectories. Compared to standard stack-of-stars sampling, both non-contrast-enhanced and DCE liver MR images acquired with VD-stack-of-stars sampling presented improved overall image quality, sharper liver edges and increased hepatic vessel clarity in all image planes. The results have suggested that the proposed VD-stack-of-stars sampling scheme can achieve improved performance (increased slice resolution or volumetric coverage with better image quality) over standard stack-of-stars sampling in free-breathing DCE-MRI without increasing scan time. The reformatted coronal and sagittal images with better slice resolution may provide added clinical value.