Optimization and validation of accelerated golden-angle radial sparse MRI reconstruction with self-calibrating GRAPPA operator gridding

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.

High-resolution free-breathing hepatobiliary phase MRI of the liver using XD-GRASP

PURPOSE

To evaluate the performance of high-resolution free-breathing (FB) hepatobiliary phase imaging of the liver using the eXtra-Dimension Golden-angle RAdial Sparse Parallel (XD-GRASP) MRI technique.

METHODS

Fifty-eight clinical patients (41 males, mean age = 52.9 ± 12.9) with liver lesions who underwent dynamic contrast-enhanced MRI with a liver-specific contrast agent were prospectively recruited for this study. Both breath-hold volumetric interpolated examination (BH-VIBE) imaging and FB imaging were performed during the hepatobiliary phase. FB images were acquired using a stack-of-stars golden-angle radial sequence and were reconstructed using the XD-GRASP method. Two experienced radiologists blinded to acquisition schemes independently scored the overall image quality, liver edge sharpness, hepatic vessel clarity, conspicuity of lesion, and overall artifact level of each image. The non-parametric paired two-tailed Wilcoxon signed-rank test was used for statistical analysis.

RESULTS

Compared to BH-VIBE images, XD-GRASP images received significantly higher scores (P < 0.05) for the liver edge sharpness (4.83 ± 0.45 vs 4.29 ± 0.46), the hepatic vessel clarity (4.64 ± 0.67 vs 4.15 ± 0.56) and the conspicuity of lesion (4.75 ± 0.53 vs 4.31 ± 0.50). There were no significant differences (P > 0.05) between BH-VIBE and XD-GRASP images for the overall image quality (4.61 ± 0.50 vs 4.74 ± 0.47) and the overall artifact level (4.13 ± 0.44 vs 4.05 ± 0.61).

CONCLUSION

Compared to conventional BH-VIBE MRI, FB radial acquisition combined with XD-GRASP reconstruction facilitates higher spatial resolution imaging of the liver during the hepatobiliary phase. This enhancement can significantly improve the visualization and evaluation of the liver.

T1 mapping for Head and Neck Cancer Patients undergoing Chemoradiotherapy: Feasibility of 3D Stack of Star Imaging

BACKGROUND

Measuring tissue oxygen concentration is crucial in understanding the pathophysiological process of hypoxia in head and neck cancer (HNC) and its significant role in cancer biology. This study aimed to determine the feasibility of T1 mapping using a variable flip angle (VFA) technique with stack of stars (SOS) trajectory sampling in HNC patients undergoing chemoradiotherapy (CRT).

METHODS

To evaluate the ability of SOS acquisition to detect T1, a phantom study was conducted and compared to conventional Cartesian acquisition (CART). Additionally, four newly diagnosed patients were recruited and underwent two scans each at baseline and inter-treatment. The repeatability of SOS and CART acquisitions was assessed by comparing the T1 measurements of CSF from the baseline and intra-treatment MRI studies. The changes in ∆T1 of the tumors during air and oxygen inhalation between baseline and inter-treatment scans were also evaluated.

RESULTS

Our study found that the 3D VFA SOS sequence was effective in reducing motion artifacts compared to the conventional VFA sequence with CART sampling and the same scan time, as demonstrated by the results from the phantom and patient studies. In terms of repeatability, no significant correlation was observed between the variability in ΔT1 measurements of CSF obtained from SOS T1 maps. The SOS ΔT1 measurements showed higher consistency, as evidenced by the ICC values ranging from 0.52 to 0.92. The ∆T1 measurements on the primary tumors increased after the first CRT (p<0.05) for all patients who showed a positive treatment response, except for one patient (0.05

CONCLUSION

The 3D VFA SOS sequence is a feasible and reliable method for T1 mapping in HNC patients undergoing CRT. The use of this technique could potentially aid in the assessment of treatment response and contribute to improving patient outcomes.

Collapse

Key Words

• 3D Stack of Star
• Head and Neck Cancer
• Oxygen-Enhancement
• T1 measurements
• Tumor Hypoxia
• Variable Flip Angle Method

Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Seong-Eun Kim Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, UT, USA.
John A Roberts Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, UT, USA
Eugene G Kholmovski Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, UT, USA; Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
Ying Hitchcock Department of Radiation Oncology, University of Utah, Salt Lake City, UT, USA
Yoshimi Anzai Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, UT, USA

Collapse

Assessment of Beat-To-Beat Variability in Left Atrial Hemodynamics Using Real Time Phase Contrast MRI in Patients With Atrial Fibrillation

BACKGROUND

Hemodynamic assessment of left atrial (LA) flow using phase contrast MRI provides insight into thromboembolic risk in atrial fibrillation (AF). However, conventional flow imaging techniques are averaged over many heartbeats.

PURPOSE

To evaluate beat-to-beat variability and LA hemodynamics in patients with AF using real time phase contrast (RTPC) MRI.

STUDY TYPE

Prospective.

SUBJECTS

Thirty-five patients with history of AF (68 ± 10 years, nine female), 10 healthy controls (57 ± 19 years, four female).

FIELD STRENGTH/SEQUENCE

5T, 2D RTPC with through-plane velocity-encoded gradient echo sequence and 4D flow MRI with three-directional velocity-encoded gradient echo sequence.

ASSESSMENT

RTPC was continuously acquired for a mid-LA slice in all subjects. 4D flow data were interpolated at the RTPC location and normally projected for comparison with RTPC. RR intervals extracted from RTPC were used to calculate heart rate variability (HRV = interquartile range over median × 100%). Patients were classified into low (<9.7%) and high (>9.7%) HRV groups. LA peak/mean velocity and stasis (%velocities < 5.8 cm/sec) were calculated from segmented 2D images. Variability in RTPC flow metrics was quantified by coefficient of variation (CV) over all cycles.

STATISTICAL TESTS

Pearson's correlation coefficient (r), Bland-Altman analysis, Kruskal-Wallis test. A P value < 0.05 was considered statistically significant.

RESULTS

RTPC and 4D flow measurements were strongly/significantly correlated for all hemodynamic parameters (R2  = 0.75-0.83) in controls. Twenty-four patients had low HRV (mean = 4 ± 2%) and 11 patients had high HRV (27 ± 9%). In patients, increased HRV was significantly correlated with CV of peak velocity (r = 0.67), mean velocity (r = 0.51), and stasis (r = 0.41). A stepwise decrease in peak/mean velocity and increase in stasis was observed when comparing controls vs. low HRV vs. high HRV groups. Mean velocity and stasis differences were significant for control vs. high HRV groups.

CONCLUSIONS

RTPC may be suitable for assessing the impact of HRV on hemodynamics and provide insight for AF management in highly arrhythmic patients.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 2.

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.

Deep Learning-Based Reconstruction for Cardiac MRI: A Review

Cardiac magnetic resonance (CMR) is an essential clinical tool for the assessment of cardiovascular disease. Deep learning (DL) has recently revolutionized the field through image reconstruction techniques that allow unprecedented data undersampling rates. These fast acquisitions have the potential to considerably impact the diagnosis and treatment of cardiovascular disease. Herein, we provide a comprehensive review of DL-based reconstruction methods for CMR. We place special emphasis on state-of-the-art unrolled networks, which are heavily based on a conventional image reconstruction framework. We review the main DL-based methods and connect them to the relevant conventional reconstruction theory. Next, we review several methods developed to tackle specific challenges that arise from the characteristics of CMR data. Then, we focus on DL-based methods developed for specific CMR applications, including flow imaging, late gadolinium enhancement, and quantitative tissue characterization. Finally, we discuss the pitfalls and future outlook of DL-based reconstructions in CMR, focusing on the robustness, interpretability, clinical deployment, and potential for new methods.

4D Golden-Angle Radial MRI at Subsecond Temporal Resolution

Intraframe motion blurring, as a major challenge in free-breathing dynamic MRI, can be reduced if high temporal resolution can be achieved. To address this challenge, this work proposes a highly accelerated 4D (3D + time) dynamic MRI framework with subsecond temporal resolution that does not require explicit motion compensation. The method combines standard stack-of-stars golden-angle radial sampling and tailored GRASP-Pro (Golden-angle RAdial Sparse Parallel imaging with imProved performance) reconstruction. Specifically, 4D dynamic MRI acquisition is performed continuously without motion gating or sorting. The k-space centers in stack-of-stars radial data are organized to guide estimation of a temporal basis, with which GRASP-Pro reconstruction is employed to enforce joint low-rank subspace and sparsity constraints. This new basis estimation strategy is the new feature proposed for subspace-based reconstruction in this work to achieve high temporal resolution (e.g., subsecond/3D volume). It does not require sequence modification to acquire additional navigation data, it is compatible with commercially available stack-of-stars sequences, and it does not need an intermediate reconstruction step. The proposed 4D dynamic MRI approach was tested in abdominal motion phantom, free-breathing abdominal MRI, and dynamic contrast-enhanced MRI (DCE-MRI). Our results have shown that GRASP-Pro reconstruction with the new basis estimation strategy enables highly-accelerated 4D dynamic imaging at subsecond temporal resolution (with five spokes or less for each dynamic frame per image slice) for both free-breathing non-DCE-MRI and DCE-MRI. In the abdominal phantom, better image quality with lower root mean square error and higher structural similarity index was achieved using GRASP-Pro compared with standard GRASP. With the ability to acquire each 3D image in less than 1 s, intraframe respiratory blurring can be intrinsically reduced for body applications with our approach, which eliminates the need for explicit motion detection and motion compensation.

Repeatability and robustness of MP-GRASP T1 mapping

PURPOSE

To demonstrate the repeatability of fast 3D T1 mapping using Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI and its robustness to variation of imaging parameters including flip angle and spatial resolution in phantoms and the brain.

THEORY AND METHODS

Multiple imaging experiments were performed to (1) assess the robustness of MP-GRASP T1 mapping to B1 inhomogeneity using a single tube phantom filled with uniform MnCl2 liquid; (2) compare the repeatability of T1 mapping between MP-GRASP and inversion recovery-based spin-echo (IR-SE; over 12 scans), using a commercial T1MES phantom; (3) evaluate the longitudinal variation of T1 estimation using MP-GRASP with varying imaging parameters, including spatial resolution, flip angle, TR/TE, and acceleration rate, using the T1MES phantom (106 scans performed over a period of 12 months); and (4) evaluate the variation of T1 estimation using MP-GRASP with varying imaging parameters in the brain (24 scans in a single visit). In addition, the accuracy of MP-GRASP T1 mapping was also validated against IR-SE by performing linear correlation and calculating the Lin's concordance correlation coefficient (CCC).

RESULTS

MP-GRASP demonstrates good robustness to B1 inhomogeneity, with intra-slice variability below 1% in the single tube phantom experiment. The longitudinal variability is good both in the phantom (below 2.5%) and in the brain (below 2%) with varying imaging parameters. The T1 values estimated from MP-GRASP are accurate compared to that from the IR-SE imaging (R2 = 0.997, Lin's CCC = 0.996).

CONCLUSION

MP-GRASP shows excellent repeatability of T1 estimation over time, and it is also robust to variation of different imaging parameters evaluated in this study.

Golden-Angle Radial MRI: Basics, Advances, and Applications

In recent years, golden‐angle radial sampling has received substantial attention and interest in the magnetic resonance imaging (MRI) community, and it has become a popular sampling trajectory for both research and clinical use. However, although the number of relevant techniques and publications has grown rapidly, there is still a lack of a review paper that provides a comprehensive overview and summary of the basics of golden‐angle rotation, the advantages and challenges/limitations of golden‐angle radial sampling, and recommendations in using different types of golden‐angle radial trajectories for MRI applications. Such a review paper is expected to be helpful both for clinicians who are interested in learning the potential benefits of golden‐angle radial sampling and for MRI physicists who are interested in exploring this research direction. The main purpose of this review paper is thus to present an overview and summary about golden‐angle radial MRI sampling. The review consists of three sections. The first section aims to answer basic questions such as: what is a golden angle; how is the golden angle calculated; why is golden‐angle radial sampling useful, and what are its limitations. The second section aims to review more advanced trajectories of golden‐angle radial sampling, including tiny golden‐angle rotation, stack‐of‐stars golden‐angle radial sampling, and three‐dimensional (3D) kooshball golden‐angle radial sampling. Their respective advantages and limitations and potential solutions to address these limitations are also discussed. Finally, the third section reviews MRI applications that can benefit from golden‐angle radial sampling and provides recommendations to readers who are interested in implementing golden‐angle radial trajectories in their MRI studies.

Deep learning for radial SMS myocardial perfusion reconstruction using the 3D residual booster U-net

PURPOSE

To develop an end-to-end deep learning solution for quickly reconstructing radial simultaneous multi-slice (SMS) myocardial perfusion datasets with comparable quality to the pixel tracking spatiotemporal constrained reconstruction (PT-STCR) method.

METHODS

Dynamic contrast enhanced (DCE) radial SMS myocardial perfusion data were obtained from 20 subjects who were scanned at rest and/or stress with or without ECG gating using a saturation recovery radial CAIPI turboFLASH sequence. Input to the networks consisted of complex coil combined images reconstructed using the inverse Fourier transform of undersampled radial SMS k-space data. Ground truth images were reconstructed using the PT-STCR pipeline. The performance of the residual booster 3D U-Net was tested by comparing it to state-of-the-art network architectures including MoDL, CRNN-MRI, and other U-Net variants.

RESULTS

Results demonstrate significant improvements in speed requiring approximately 8 seconds to reconstruct one radial SMS dataset which is approximately 200 times faster than the PT-STCR method. Images reconstructed with the residual booster 3D U-Net retain quality of ground truth PT-STCR images (0.963 SSIM/40.238 PSNR/0.147 NRMSE). The residual booster 3D U-Net has superior performance compared to existing network architectures in terms of image quality, temporal dynamics, and reconstruction time.

CONCLUSION

Residual and booster learning combined with the 3D U-Net architecture was shown to be an effective network for reconstructing high-quality images from undersampled radial SMS datasets while bypassing the reconstruction time of the PT-STCR method.

Magnetization-prepared GRASP MRI for rapid 3D T1 mapping and fat/water-separated T1 mapping

PURPOSE

This study aimed to (i) develop Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI using a stack-of-stars trajectory for rapid free-breathing T1 mapping and (ii) extend MP-GRASP to multi-echo acquisition (MP-Dixon-GRASP) for fat/water-separated (water-specific) T1 mapping.

METHODS

An adiabatic non-selective 180° inversion-recovery pulse was added to a gradient-echo-based golden-angle stack-of-stars sequence for magnetization-prepared 3D single-echo or 3D multi-echo acquisition. In combination with subspace-based GRASP-Pro reconstruction, the sequence allows for standard T1 mapping (MP-GRASP) or fat/water-separated T1 mapping (MP-Dixon-GRASP), respectively. The accuracy of T1 mapping using MP-GRASP was evaluated in a phantom and volunteers (brain and liver) against clinically accepted reference methods. The repeatability of T1 estimation was also assessed in the phantom and volunteers. The performance of MP-Dixon-GRASP for water-specific T1 mapping was evaluated in a fat/water phantom and volunteers (brain and liver).

RESULTS

ROI-based mean T1 values are correlated between the references and MP-GRASP in the phantom (R2 = 1.0), brain (R2 = 0.96), and liver (R2 = 0.73). MP-GRASP achieved good repeatability of T1 estimation in the phantom (R2 = 1.0), brain (R2 = 0.99), and liver (R2 = 0.82). Water-specific T1 is different from in-phase and out-of-phase composite T1 (composite T1 when fat and water signal are mixed in phase or out of phase) both in the phantom and volunteers.

CONCLUSION

This work demonstrated the initial performance of MP-GRASP and MP-Dixon-GRASP MRI for rapid 3D T1 mapping and 3D fat/water-separated T1 mapping in the brain (without motion) and in the liver (during free breathing). With fat/water-separated T1 estimation, MP-Dixon-GRASP could be potentially useful for imaging patients with fatty-liver diseases.

Intraindividual Comparison of Compressed Sensing-Accelerated Cartesian and Radial Arterial Phase Imaging of the Liver in an Experimental Tumor Model

OBJECTIVES

The aim of this study was to intraindividually compare the performance of 2 compressed sensing (CS)-accelerated magnetic resonance imaging (MRI) sequences, 1 featuring Cartesian (compressed sensing volumetric interpolated breath-hold examination [CS-VIBE]) and the other radial (golden-angle radial sparse parallel [GRASP]) k-space sampling in continuous dynamic imaging during hepatic vascular phases, using extracellular and hepatocyte-specific contrast agents.

MATERIALS AND METHODS

Seven New Zealand white rabbits, with induced VX2 liver tumors (median number of lesions, 2 ± 0.83; range, 1-3), received 2 continuously acquired T1-weighted prototype CS-accelerated MRI sequences (CS-VIBE and GRASP) with high spatial (0.8 × 0.8 × 1.5 mm) and temporal resolution (3.5 seconds) in randomized order on 2 separate days using a 1.5-T scanner. In all animals, imaging was performed using first gadobutrol at a dose of 0.1 mmol/kg and, then 45 minutes later, gadoxetic acid at a dose of 0.025 mmol/kg.The following qualitative parameters were assessed using 3- and 5-point Likert scales (3 and 5 being the highest scores respectively): image quality (IQ), arterial and venous vessel delineation, tumor enhancement, motion artifacts, and sequence-specific artifacts. Furthermore, the following quantitative parameters were obtained: relative peak signal enhancement, time to peak, mean transit time, and plasma flow ratios. Paired sampled t tests and Wilcoxon signed rank tests were used for intraindividual comparison. Image analysis was performed by 2 radiologists.

RESULTS

Six of 7 animals underwent the full imaging protocol and obtained data were analyzed statistically. Overall IQ was rated moderate to excellent, not differing significantly between the 2 sequences.Gadobutrol-enhanced CS-VIBE examinations revealed the highest mean Likert scale values in terms of vessel delineation and tumor enhancement (arterial 4.4 [4-5], venous 4.3 [3-5], and tumor 2.9 [2-3]). Significantly, more sequence-specific artifacts were seen in GRASP examinations (P = 0.008-0.031). However, these artifacts did not impair IQ. Excellent Likert scale ratings were found for motion artifacts in both sequences. In both sequences, a maximum of 4 hepatic arterial dominant phases were obtained. Regarding the relative peak signal enhancement, CS-VIBE and GRASP showed similar results. The relative peak signal enhancement values did not differ significantly between the 2 sequences in the aorta, the hepatic artery, or the inferior vena cava (P = 0.063-0.536). However, significantly higher values were noted for CS-VIBE in gadoxetic acid-enhanced examinations in the portal vein (P = 0.031) and regarding the tumor enhancement (P = 0.005). Time to peak and mean transit time or plasma flow ratios did not differ significantly between the sequences.

CONCLUSIONS

Both CS-VIBE and GRASP provide excellent results in dynamic liver MRI using extracellular and hepatocyte-specific contrast agents, in terms of IQ, peak signal intensity, and presence of artifacts.

Artifact reduction in free-breathing, free-running myocardial perfusion imaging with interleaved non-selective RF excitations

PURPOSE

To reduce inflow and motion artifacts in free-breathing, free-running, steady-state spoiled gradient echo T₁ -weighted (SPGR) myocardial perfusion imaging.

METHOD

Unsaturated spins from inflowing blood or out-of-plane motion cause flashing artifacts in free-running SPGR myocardial perfusion. During free-running SPGR, 1 non-selective RF excitation was added after every 3 slice-selective RF excitations to suppress inflow artifacts by forcing magnetization in neighboring regions to steady-state. Bloch simulations and phantom experiments were performed to evaluate the impact of the flip angle and non-selective RF frequency on inflowing spins and tissue contrast. Free-running perfusion with (n = 11) interleaved non-selective RF or without (n = 11) were studied in 22 subjects (age = 60.2 ± 14.3 years, 11 male). Perfusion images were graded on a 5-point Likert scale for conspicuity of wall enhancement, inflow/motion artifact, and streaking artifact and compared using Wilcoxon sum-rank testing.

RESULT

Numeric simulation showed that 1 non-selective RF excitation applied after every 3 slice-selective RF excitations produced superior out-of-plane signal suppression compared to 1 non-selective RF excitation applied after every 6 or 9 slice-selective RF excitations. In vitro experiments showed that a 30° flip angle produced near-optimal myocardial contrast. In vivo experiments demonstrated that the addition of interleaved non-selective RF significantly (P < .01) improved conspicuity of wall enhancement (mean score = 4.4 vs. 3.2) and reduced inflow/motion (mean score = 4.5 vs. 2.5) and streaking (mean score = 3.9 vs. 2.4) artifacts.

CONCLUSION

Non-selective RF excitations interleaved between slice-selective excitations can reduce image artifacts in free-breathing, ungated perfusion images. Further studies are warranted to assess the diagnostic accuracy of the proposed solution for evaluating myocardial ischemia.

Highly accelerated free-breathing real-time phase contrast cardiovascular MRI via complex-difference deep learning

PURPOSE

To develop and evaluate a real-time phase contrast (PC) MRI protocol via complex-difference deep learning (DL) framework.

METHODS

DL used two 3D U-nets to separately filter aliasing artifact from radial real-time velocity-compensated and complex-difference images. U-nets were trained with synthetic real-time PC generated from electrocardiograph (ECG) -gated, breath-hold, segmented PC (ECG-gated segmented PC) acquired at the ascending aorta of 510 patients. In 21 patients, free-breathing, ungated real-time (acceleration rate = 28.8) and ECG-gated segmented (acceleration rate = 2) PC were prospectively acquired at the ascending aorta. Hemodynamic parameters (cardiac output [CO], stroke volume [SV], and mean velocity at peak systole [peak mean velocity]) were measured for ECG-gated segmented and DL-filtered synthetic real-time PC and compared using Bland-Altman and linear regression analyses. Additionally, hemodynamic parameters were quantified from DL-filtered, compressed-sensing (CS) -reconstructed, and gridding reconstructed prospective real-time PC and compared to ECG-gated segmented PC.

RESULTS

Synthetic real-time PC with DL showed strong correlation (R > 0.98) and good agreement with ECG-gated segmented PC for quantified hemodynamic parameters (mean-difference: CO = -0.3 L/min, SV = -4.3 mL, peak mean velocity = -2.3 cm/s). On average, DL required 0.39 s/frame to filter prospective real-time PC, which was 4.6-fold faster than CS. Compared to CS, DL showed superior correlation, tighter limits of agreement (LOAs), better bias for peak mean velocity, and worse bias for CO and SV. Compared to gridding, DL showed similar correlation, tighter LOAs for CO and SV, similar bias for CO, and worse bias for SV and peak mean velocity.

CONCLUSION

The complex-difference DL framework accelerated real-time PC-MRI by nearly 28-fold, enabling rapid free-running real-time assessment of flow hemodynamics.

Whole-heart, ungated, free-breathing, cardiac-phase-resolved myocardial perfusion MRI by using Continuous Radial Interleaved simultaneous Multi-slice acquisitions at sPoiled steady-state (CRIMP)

PURPOSE

To develop a whole-heart, free-breathing, non-electrocardiograph (ECG)-gated, cardiac-phase-resolved myocardial perfusion MRI framework (CRIMP; Continuous Radial Interleaved simultaneous Multi-slice acquisitions at sPoiled steady-state) and test its quantification feasibility.

METHODS

CRIMP used interleaved radial simultaneous multi-slice (SMS) slice groups to cover the whole heart in 9 or 12 short-axis slices. The sequence continuously acquired data without magnetization preparation, ECG gating or breath-holding, and captured multiple cardiac phases. Images were reconstructed by a motion-compensated patch-based locally low-rank reconstruction. Bloch simulations were performed to study the signal-to-noise ratio/contrast-to-noise ratio (SNR/CNR) for CRIMP and to study the steady-state signal under motion. Seven patients were scanned with CRIMP at stress and rest to develop the sequence. One human and two dogs were scanned at rest with a dual-bolus method to test the quantification feasibility of CRIMP. The dual-bolus scans were performed using both CRIMP and an ungated radial SMS saturation recovery (SMS-SR) sequence with injection dose = 0.075 mmol/kg to compare the sequences in terms of SNR, cardiac phase resolution and quantitative myocardial blood flow (MBF).

RESULTS

Perfusion images with multiple cardiac phases in all image slices with a temporal resolution of 72 ms/frame were obtained. Simulations and in-vivo acquisitions showed CRIMP kept the inner slices in steady-state regardless of motion. CRIMP outperformed SMS-SR in slice coverage (9 over 6), SNR (mean 20% improvement), and provided cardiac phase resolution. CRIMP and SMS-SR sequences provided comparable MBF values (rest systolic CRIMP = 0.58 ± 0.07, SMS-SR = 0.61 ± 0.16).

CONCLUSION

CRIMP allows for whole-heart, cardiac-phase-resolved myocardial perfusion images without ECG-gating or breath-holding. The sequence can provide MBF if an accurate arterial input function is obtained separately.

Quantitative 3D myocardial perfusion with an efficient arterial input function

PURPOSE

The purpose of this study was to further develop and combine several innovative sequence designs to achieve quantitative 3D myocardial perfusion. These developments include an optimized 3D stack-of-stars readout (150 ms per beat), efficient acquisition of a 2D arterial input function, tailored saturation pulse design, and potential whole heart coverage during quantitative stress perfusion.

THEORY AND METHODS

All studies were performed free-breathing on a Prisma 3T MRI scanner. Phantom validation was used to verify sequence accuracy. A total of 21 subjects (3 patients with known disease) were scanned, 12 with a rest only protocol and 9 with both stress (regadenoson) and rest protocols. First pass quantitative perfusion was performed with gadoteridol (0.075 mmol/kg).

RESULTS

Implementation and quantitative perfusion results are shown for healthy subjects and subjects with known coronary disease. Average rest perfusion for the 15 included healthy subjects was 0.79 ± 0.19 mL/g/min, the average stress perfusion for 6 healthy subject studies was 2.44 ± 0.61 mL/g/min, and the average global myocardial perfusion reserve ratio for 6 healthy subjects was 3.10 ± 0.24. Perfusion deficits for 3 patients with ischemia are shown. Average resting heart rate was 59 ± 7 bpm and the average stress heart rate was 81 ± 10 bpm.

CONCLUSION

This work demonstrates that a quantitative 3D myocardial perfusion sequence with the acquisition of a 2D arterial input function is feasible at high stress heart rates such as during stress. T1 values and gadolinium concentrations of the sequence match the reference standard well in a phantom, and myocardial rest and stress perfusion and myocardial perfusion reserve values are consistent with those published in literature.

Highly accelerated, real-time phase-contrast MRI using radial k-space sampling and GROG-GRASP reconstruction: a feasibility study in pediatric patients with congenital heart disease

Retrospective electrocardiogram-gated, 2D phase-contrast (PC) flow MRI is routinely used in clinical evaluation of valvular/vascular disease in pediatric patients with congenital heart disease (CHD). In patients not requiring general anesthesia, clinical standard PC is conducted with free breathing for several minutes per slice with averaging. In younger patients under general anesthesia, clinical standard PC is conducted with breath-holding. One approach to overcome this limitation is using either navigator gating or self-navigation of respiratory motion, at the expense of lengthening scan times. An alternative approach is using highly accelerated, free-breathing, real-time PC (rt-PC) MRI, which to date has not been evaluated in CHD patients. The purpose of this study was to develop a 38.4-fold accelerated 2D rt-PC pulse sequence using radial k-space sampling and compressed sensing with 1.5 × 1.5 × 6.0 mm³ nominal spatial resolution and 40 ms nominal temporal resolution, and evaluate whether it is capable of accurately measuring flow in 17 pediatric patients (aortic valve, pulmonary valve, right and left pulmonary arteries) compared with clinical standard 2D PC (either breath-hold or free breathing). For clinical translation, we implemented an integrated reconstruction pipeline capable of producing DICOMs of the order of 2 min per time series (46 frames). In terms of association, forward volume, backward volume, regurgitant fraction, and peak velocity at peak systole measured with standard PC and rt-PC were strongly correlated (R² > 0.76; P < 0.001). Compared with clinical standard PC, in terms of agreement, forward volume (mean difference = 1.4% (3.0% of mean)) and regurgitant fraction (mean difference = -2.5%) were in good agreement, whereas backward volume (mean difference = -1.1 mL (28.2% of mean)) and peak-velocity at peak systole (mean difference = -21.3 cm/s (17.2% of mean)) were underestimated by rt-PC. This study demonstrates that the proposed rt-PC with the said spatial resolution and temporal resolution produces relatively accurate forward volumes and regurgitant fractions but underestimates backward volumes and peak velocities at peak systole in pediatric patients with CHD.

Isotropic 3D Cartesian single breath-hold CINE MRI with multi-bin patch-based low-rank reconstruction

PURPOSE

To develop a novel acquisition and reconstruction framework for isotropic 3D Cartesian cardiac CINE within a single breath-hold for left ventricle (LV) and whole-heart coverage.

METHODS

A variable-density Cartesian acquisition with spiral profile ordering, out-inward sampling, and acquisition-adaptive alternating tiny golden/golden angle increment between spiral arms is proposed to provide incoherent and nonredundant sampling within and among cardiac phases. A novel multi-bin patch-based low-rank reconstruction, named MB-PROST, is proposed to exploit redundant information on a local (within a patch), nonlocal (similar patches within a spatial neighborhood), and temporal (among all cardiac phases) scale with an implicit motion alignment among patches. The proposed multi-bin patch-based low-rank reconstruction reconstruction is compared against compressed sensing reconstruction, whereas LV function parameters derived from the proposed 3D CINE framework are compared against those estimated from conventional multislice 2D CINE imaging in 10 healthy subjects and 15 patients.

RESULTS

The proposed framework provides 3D cardiac CINE images with high spatial (1.9 mm³ ) and temporal resolution (˜50 ms) in a single breath-hold of ˜20 s for LV and ˜26 s for whole-heart coverage in healthy subjects. Shorter breath-hold durations of ˜13 to 15 s are feasible for LV coverage with slightly anisotropic resolution (1.9 × 1.9 × 2.5 mm) in patients. LV function parameters derived from 3D CINE were in good agreement with 2D CINE, with a bias of -0.1 mL/0.1 mL, -0.9 mL/-1.0 mL, -0.1%/-0.8%; and confidence intervals of ±1.7 mL/±3.7 mL, ±1.2 mL/±2.6 mL, and ±1.2%/±3.6% (10 healthy subjects/15 patients) for end-systolic volume, end-diastolic volume, and ejection fraction, respectively.

CONCLUSION

The proposed framework enables 3D isotropic cardiac CINE in a single breath-hold scan of ˜20 s/˜26 s for LV/whole-heart coverage, showing good agreement with clinical 2D CINE scans in terms of LV functional assessment.

MRSIGMA: Magnetic Resonance SIGnature MAtching for real-time volumetric imaging

PURPOSE

To propose a real-time 3D MRI technique called MR SIGnature MAtching (MRSIGMA) for high-resolution volumetric imaging and motion tracking with very low imaging latency.

METHODS

MRSIGMA consists of two steps: (1) offline learning of a database of possible 3D motion states and corresponding motion signature ranges and (2) online matching of new motion signatures acquired in real time with prelearned motion states. Specifically, the offline learning step (non-real-time) reconstructs motion-resolved 4D images representing different motion states and assigns a unique motion range to each state. The online matching step (real-time) acquires motion signatures only and selects one of the prelearned 3D motion states for each newly acquired signature, which generates 3D images efficiently in real time. The MRSIGMA technique was evaluated on 15 golden-angle stack-of-stars liver data sets, and the performance of respiratory motion tracking with the online-generated real-time 3D MRI was compared with the corresponding 2D projections acquired in real time.

RESULTS

The total latency of generating each 3D image during online matching was about 300 ms, including acquisition of the motion signature data (~138 ms) and corresponding matching process (~150 ms). Linear correlation assessment suggested excellent correlation (R²  = 0.948) between motion displacement measured from the online-generated real-time 3D images and the 2D real-time projections.

CONCLUSION

This proof-of-concept study demonstrates the feasibility of MRSIGMA for high-resolution real-time volumetric imaging, which shifts the acquisition and reconstruction burden to an offline learning step and leaves fast online matching for online imaging with very low imaging latency. The MRSIGMA technique can potentially be used for real-time motion tracking in MRI-guided radiation therapy.

Improving dynamic contrast-enhanced MRI of the lung using motion-weighted sparse reconstruction: Initial experiences in patients

PURPOSE

The purpose of this study was to evaluate the performance of motion-weighted Golden-angle RAdial Sparse Parallel MRI (motion-weighted GRASP) for free-breathing dynamic contrast-enhanced MRI (DCE-MRI) of the lung.

METHODS

Motion-weighted GRASP incorporates a soft-gating motion compensation algorithm into standard GRASP reconstruction, so that motion-corrupted motion k-space (e.g., k-space acquired in inspiratory phases) contributes less to the final reconstructed images. Lung MR data from 20 patients (mean age = 57.9 ± 13.5) with known pulmonary lesions were retrospectively collected for this study. Each subject underwent a free-breathing DCE-MR scan using a fat-statured T₁-weighted stack-of-stars golden-angle radial sequence and a post-contrast breath-hold MR scan using a Cartesian volumetric-interpolated imaging sequence (BH-VIBE). Each radial dataset was reconstructed using GRASP without motion compensation and motion-weighted GRASP. All MR images were visually evaluated by two experienced radiologists blinded to reconstruction and acquisition schemes independently. In addition, the influence of motion-weighted reconstruction on dynamic contrast-enhancement patterns was also investigated.

RESULTS

For image quality assessment, motion-weighted GRASP received significantly higher visual scores than GRASP (P < 0.05) for overall image quality (3.68 vs. 3.39), lesion conspicuity (3.54 vs. 3.18) and overall artifact level (3.53 vs. 3.15). There was no significant difference (P > 0.05) between the breath-hold BH-VIBE and motion-weighted GRASP images. For assessment of temporal fidelity, motion-weighted GRASP maintained a good agreement with respect to GRASP.

CONCLUSION

Motion-weighted GRASP achieved better reconstruction performance in free-breathing DCE-MRI of the lung compared to standard GRASP, and it may enable improved assessment of pulmonary lesions.

SANTIS: Sampling-Augmented Neural neTwork with Incoherent Structure for MR image reconstruction

PURPOSE

To develop and evaluate a novel deep learning-based reconstruction framework called SANTIS (Sampling-Augmented Neural neTwork with Incoherent Structure) for efficient MR image reconstruction with improved robustness against sampling pattern discrepancy.

METHODS

With a combination of data cycle-consistent adversarial network, end-to-end convolutional neural network mapping, and data fidelity enforcement for reconstructing undersampled MR data, SANTIS additionally utilizes a sampling-augmented training strategy by extensively varying undersampling patterns during training, so that the network is capable of learning various aliasing structures and thereby removing undersampling artifacts more effectively and robustly. The performance of SANTIS was demonstrated for accelerated knee imaging and liver imaging using a Cartesian trajectory and a golden-angle radial trajectory, respectively. Quantitative metrics were used to assess its performance against different references. The feasibility of SANTIS in reconstructing dynamic contrast-enhanced images was also demonstrated using transfer learning.

RESULTS

Compared to conventional reconstruction that exploits image sparsity, SANTIS achieved consistently improved reconstruction performance (lower errors and greater image sharpness). Compared to standard learning-based methods without sampling augmentation (e.g., training with a fixed undersampling pattern), SANTIS provides comparable reconstruction performance, but significantly improved robustness, against sampling pattern discrepancy. SANTIS also achieved encouraging results for reconstructing liver images acquired at different contrast phases.

CONCLUSION

By extensively varying undersampling patterns, the sampling-augmented training strategy in SANTIS can remove undersampling artifacts more robustly. The novel concept behind SANTIS can particularly be useful for improving the robustness of deep learning-based image reconstruction against discrepancy between training and inference, an important, but currently less explored, topic.

GRASP-Pro: imProving GRASP DCE-MRI through self-calibrating subspace-modeling and contrast phase automation

PURPOSE

To propose a highly accelerated, high-resolution dynamic contrast-enhanced MRI (DCE-MRI) technique called GRASP-Pro (golden-angle radial sparse parallel imaging with imProved performance) through a joint sparsity and self-calibrating subspace constraint with automated selection of contrast phases.

METHODS

GRASP-Pro reconstruction enforces a combination of an explicit low-rank subspace-constraint and a temporal sparsity constraint. The temporal basis used to construct the subspace is learned from an intermediate reconstruction step using the low-resolution portion of radial k-space, which eliminates the need for generating the basis using auxiliary data or a physical signal model. A convolutional neural network was trained to generate the contrast enhancement curve in the artery, from which clinically relevant contrast phases are automatically selected for evaluation. The performance of GRASP-Pro was demonstrated for high spatiotemporal resolution DCE-MRI of the prostate and was compared against standard GRASP in terms of overall image quality, image sharpness, and residual streaks and/or noise level.

RESULTS

Compared to GRASP, GRASP-Pro reconstructed dynamic images with enhanced sharpness, less residual streaks and/or noise, and finer delineation of the prostate without prolonging reconstruction time. The image quality improvement reached statistical significance (P < 0.05) in all the assessment categories. The neural network successfully generated contrast enhancement curves in the artery, and corresponding peak enhancement indexes correlated well with that from the manual selection.

CONCLUSION

GRASP-Pro is a promising method for rapid and continuous DCE-MRI. It enables superior reconstruction performance over standard GRASP and allows reliable generation of artery enhancement curve to guide the selection of desired contrast phases for improving the efficiency of GRASP MRI workflow.

Feasibility of multiple-view myocardial perfusion MRI using radial simultaneous multi-slice acquisitions

Purpose

Dynamic contrast enhanced MRI of the heart typically acquires 2–4 short-axis (SA) slices to detect and characterize coronary artery disease. This acquisition scheme is limited by incomplete coverage of the left ventricle. We studied the feasibility of using radial simultaneous multi-slice (SMS) technique to achieve SA, 2-chamber and/or 4-chamber long-axis (2CH LA and/or 4CH LA) coverage with and without electrocardiography (ECG) gating using a motion-robust reconstruction framework.

Methods

12 subjects were scanned at rest and/or stress, free breathing, with or without ECG gating. Multiple sets of radial SMS k-space were acquired within each cardiac cycle, and each SMS set sampled 3 parallel slices that were either SA, 2CH LA, or 4CH LA slices. The radial data was interpolated onto Cartesian space using an SMS GRAPPA operator gridding method. Self-gating and respiratory states binning of the data were done. The binning information as well as a pixel tracking spatiotemporal constrained reconstruction method were applied to obtain motion-robust image reconstructions. Reconstructions with and without the pixel tracking method were compared for signal-to-noise ratio and contrast-to-noise ratio.

Results

Full coverage of the heart (at least 3 SA and 3 LA slices) during the first pass of contrast at every heartbeat was achieved by using the radial SMS acquisition. The proposed pixel tracking reconstruction improves the average SNR and CNR by 21% and 30% respectively, and reduces temporal blurring for both gated and ungated acquisitions.

Conclusion

Acquiring simultaneous multi-slice SA, 2CH LA and/or 4CH LA myocardial perfusion images in every heartbeat is feasible in both gated and ungated acquisitions. This can add confidence when detecting and characterizing coronary artery disease by revealing ischemia in different views, and by providing apical coverage that is improved relative to SA slices alone. The proposed pixel tracking framework improves the reconstruction while adding little computational cost.

Accelerated, free-breathing, noncontrast, electrocardiograph-triggered, thoracic MR angiography with stack-of-stars k-space sampling and GRASP reconstruction

PURPOSE

To develop an accelerated, free-breathing, noncontrast, electrocardiograph-triggered, thoracic MR angiography (NC-MRA) pulse sequence capable of achieving high spatial resolution at clinically acceptable scan time and test whether it produces clinically acceptable image quality in patients with suspected aortic disease.

METHODS

We modified a "coronary" MRA pulse sequence to use a stack-of-stars k-space sampling pattern and combined it with golden-angle radial sparse parallel (GRASP reconstruction to enable self-navigation of respiratory motion and high data acceleration. The performance of the proposed NC-MRA was evaluated in 13 patients, where clinical standard contrast-enhanced MRA (CE-MRA) was used as control. For visual analysis, two readers graded the conspicuity of vessel lumen, artifacts, and noise level on a 5-point scale (overall score index = sum of three scores). The aortic diameters were measured at seven standardized locations. The mean visual scores, inter-observer variability, and vessel diameters were compared using appropriate statistical tests.

RESULTS

The overall mean visual score index (12.1 ± 1.7 for CE-MRA versus 12.1 ± 1.0 for NC-MRA) scores were not significantly different (P > 0.16). The two readers' scores were significantly different for CE-MRA (P = 0.01) but not for NC-MRA (P = 0.21). The mean vessel diameters were not significantly different, except at the proximal aortic arch (P < 0.03). The mean diameters were strongly correlated (R² ≥ 0.96) and in good agreement (absolute mean difference ≤ 0.01 cm and 95% confidence interval ≤ 0.62 cm).

CONCLUSION

This study shows that the proposed NC-MRA produces clinically acceptable image quality in patients at high spatial resolution (1.5 mm × 1.5 mm × 1.5 mm) and clinically acceptable scan time (~6 min).