5D whole-heart sparse MRI

Reduction of ringing artifacts induced by diaphragm drifting in free-breathing dynamic pulmonary MRI using 3D koosh-ball acquisition

PURPOSE

To reduce the ringing artifacts of the motion-resolved images in free-breathing dynamic pulmonary MRI.

METHODS

A golden-step based interleaving (GSI) technique was proposed to reduce ringing artifacts induced by diaphragm drifting. The pulmonary MRI data were acquired using a superior-inferior navigated 3D radial UTE sequence in an interleaved manner during free breathing. Successive interleaves were acquired in an incoherent fashion along the polar direction. Four-dimensional images were reconstructed from the motion-resolved k-space data obtained by retrospectively binning. The reconstruction algorithms included standard nonuniform fast Fourier transform (NUFFT), Voronoi-density-compensated NUFFT, extra-dimensional UTE, and motion-state weighted motion-compensation reconstruction. The proposed interleaving technique was compared with a conventional sequential interleaving (SeqI) technique on a phantom and eight subjects.

RESULTS

The quantified ringing artifacts level in the motion-resolved image is positively correlated with the quantified nonuniformity level of the corresponding k-space. The nonuniformity levels of the end-expiratory and end-inspiratory k-space binned from GSI data (0.34 ± 0.07, 0.33 ± 0.05) are significantly lower with statistical significance (p < 0.05) than that binned from SeqI data (0.44 ± 0.11, 0.42 ± 0.12). Ringing artifacts are substantially reduced in the dynamic images of eight subjects acquired using the proposed technique in comparison with that acquired using the conventional SeqI technique.

CONCLUSION

Ringing artifacts in the motion-resolved images induced by diaphragm drifting can be reduced using the proposed GSI technique for free-breathing dynamic pulmonary MRI. This technique has the potential to reduce ringing artifacts in free-breathing liver and kidney MRI based on full-echo interleaved 3D radial acquisition.

Quiescent frame, contrast-enhanced coronary magnetic resonance angiography reconstructed using limited number of physiologic frames from 5D free-running acquisitions

BACKGROUND

5D, free-running imaging resolves sets of 3D whole-heart images in both cardiac and respiratory dimensions. In an application such as coronary imaging when a single, static image is of interest, computationally expensive offline iterative reconstruction is still needed to compute the multiple 3D datasets.

PURPOSE

Evaluate how the number of physiologic bins included in the reconstruction affects the computational cost and resulting image quality of a single, static volume reconstruction.

STUDY TYPE

Retrospective.

SUBJECTS

15 pediatric patients following Ferumoxytol infusion (4 mg/kg).

FIELD STRENGTH/SEQUENCE

1.5 T/Ungated 5D free-running GRE sequence.

ASSESSMENT

The raw data of each subject were binned and reconstructed into a 5D (x-y-z-cardiac-respiratory) images. 1, 3, 5, 7, and 9 bins adjacent to both sides of the retrospectively determined cardiac resting phase and 1, 3 bins adjacent to the end-expiration phase are used for limited frame reconstructions. The static volume within each limited reconstruction was compared with the corresponding full 5D reconstruction using the structural similarity index measure (SSIM). A non-linear regression model was used to fit SSIM with the percentage of data used compared to full reconstruction (% data). A linear regression model was used to fit computation time with % raw data used. Coronary artery sharpness is measured on each limited reconstructed images to determine the minimal number of cardiac and respiratory bins needed to preserve image quality.

STATISTICAL TESTS

The coefficient of determination (R2) is computed for each regression model.

RESULTS

The % of data used in the reconstruction was linearly related to the computational time (R2 = 0.99). The SSIM of the static image from the limited reconstructions is non-linearly related with the % of data used (R2 = 0.80). Over the 15 patients, the model showed SSIM of 0.9 with 18% of data, and SSIM of 0.96 with 30% of data. The coronary artery sharpness of images reconstructed using no less than 5 cardiac and all respiratory phases is not significantly different from the full reconstructed images using all cardiac and respiratory bins.

DATA CONCLUSION

Reconstruction using only a limited number of acquired physiological states can linearly reduce the computational cost while preserving similarity to the full reconstruction image. It is suggested to use no less than 5 cardiac and all respiratory phases in the limited reconstruction to best preserve the original quality seen on the full reconstructed images.

Free-Running Cardiac and Respiratory Motion-Resolved Imaging: A Paradigm Shift for Managing Motion in Cardiac MRI?

Cardiac magnetic resonance imaging (MRI) is widely used for non-invasive assessment of cardiac morphology, function, and tissue characteristics due to its exquisite soft-tissue contrast. However, it remains time-consuming and requires proficiency, making it costly and limiting its widespread use. Traditional cardiac MRI is inefficient as signal acquisition is often limited to specific cardiac phases and requires complex view planning, parameter adjustments, and management of both respiratory and cardiac motion. Recent efforts have aimed to make cardiac MRI more efficient and accessible. Among these innovations, the free-running framework enables 5D whole-heart imaging without the need for an electrocardiogram signal, respiratory breath-holding, or complex planning. It uses a fully self-gated approach to extract cardiac and respiratory signals directly from the acquired image data, allowing for more efficient coverage in time and space without the need for electrocardiogram gating, triggering, navigators, or breath-holds. This review provides a comprehensive overview of the free-running framework, detailing its history, concepts, recent improvements, and clinical applications.

Motion-robust free-running volumetric cardiovascular MRI

PURPOSE

To present and assess an outlier mitigation method that makes free-running volumetric cardiovascular MRI (CMR) more robust to motion.

METHODS

The proposed method, called compressive recovery with outlier rejection (CORe), models outliers in the measured data as an additive auxiliary variable. We enforce MR physics-guided group sparsity on the auxiliary variable, and jointly estimate it along with the image using an iterative algorithm. For evaluation, CORe is first compared to traditional compressed sensing (CS), robust regression (RR), and an existing outlier rejection method using two simulation studies. Then, CORe is compared to CS using seven three-dimensional (3D) cine, 12 rest four-dimensional (4D) flow, and eight stress 4D flow imaging datasets.

RESULTS

Our simulation studies show that CORe outperforms CS, RR, and the existing outlier rejection method in terms of normalized mean square error and structural similarity index across 55 different realizations. The expert reader evaluation of 3D cine images demonstrates that CORe is more effective in suppressing artifacts while maintaining or improving image sharpness. Finally, 4D flow images show that CORe yields more reliable and consistent flow measurements, especially in the presence of involuntary subject motion or exercise stress.

CONCLUSION

An outlier rejection method is presented and tested using simulated and measured data. This method can help suppress motion artifacts in a wide range of free-running CMR applications.

Assessing Aortic Motion with Automated 3D Cine Balanced SSFP MRI Segmentation

PURPOSE

To apply free-running three-dimensional (3D) cine balanced steady state free precession (bSSFP) CMR framework in combination with AI segmentations to quantify time-resolved aortic displacement, diameter and diameter change.

METHODS

In this prospective study, we implemented a free-running 3D cine bSSFP sequence with scan time of about 4minutes facilitated by pseudo-spiral Cartesian undersampling and compressed-sensing reconstruction. Automated segmentation of all cardiac timeframes was applied through the use of nnU-Net. Dynamic 3D motion maps were created for three repeated scans per volunteer, leading to the detailed quantification of motion, as well as the measurement and change in diameter of the ascending aorta.

RESULTS

A total of 14 adult healthy volunteers (median age, 28 years (IQR: 26.0-31.3), 6 female) were included. Automated segmentation compared to manual segmentation of the aorta test set showed a Dice score of 0.93 ± 0.02. The median (interquartile range) over all volunteers for the largest maximum and mean ascending aorta (AAo) displacement in the first scan was 13.0 (4.4) mm and 5.6 (2.4) mm, respectively. Peak mean diameter in the AAo was 25.9 (2.2) mm and peak mean diameter change was 1.4 (0.5) mm. The maximum individual variability over the three repeated scans of maximum and mean AAo displacement was 3.9 (1.6) mm and 2.2 (0.8) mm, respectively. The maximum individual variability of mean diameter and diameter change were 1.2 (0.5) mm and 0.9 (0.4) mm.

CONCLUSION

A free-running 3D cine bSSFP CMR scan with a scan time of four minutes combined with an automated nnU-net segmentation consistently captured the aorta's cardiac motion-related 4D displacement, diameter, and diameter change.

Free-breathing time-resolved 4D MRI with improved T1-weighting contrast

This work proposes MP-Grasp4D (magnetization-prepared golden-angle radial sparse parallel 4D) MRI, a free-breathing, inversion recovery (IR)-prepared, time-resolved 4D MRI technique with improved T1-weighted contrast. MP-Grasp4D MRI acquisition incorporates IR preparation into a radial gradient echo sequence. MP-Grasp4D employs a golden-angle navi-stack-of-stars sampling scheme, where imaging data of rotating radial stacks and navigator stacks (acquired at a consistent rotation angle) are alternately acquired. The navigator stacks are used to estimate a temporal basis for low-rank subspace-constrained reconstruction. This allows for the simultaneous capture of both IR-induced contrast changes and respiratory motion. One temporal frame of the imaging volume in MP-Grasp4D MRI is reconstructed from a single stack and an adjacent navigator stack on average, resulting in a nominal temporal resolution of 0.16 seconds per volume. Images corresponding to the optimal inversion time (TI) can be retrospectively selected for providing the best image contrast. Reader studies were conducted to assess the performance of MP-Grasp4D MRI in liver imaging across 30 subjects in comparison with standard Grasp4D MRI without IR preparation. MP-Grasp4D MRI received significantly higher scores (P < 0.05) than Grasp4D in all assessment categories. There was a moderate to almost perfect agreement (kappa coefficient from 0.42 to 0.9) between the two readers for image quality assessment. When the scan time is reduced, MP-Grasp4D MRI preserves image contrast and quality, demonstrating additional acceleration capability. MP-Grasp4D MRI improves T1-weighted contrast for free-breathing time-resolved 4D MRI and eliminates the need for explicit motion compensation. This method is expected to be valuable in different MRI applications such as MR-guided radiotherapy.

Intrapericardial Administration to Achieve Localized and Targeted Treatment for Cardiac Disease

Intrapericardial administration has been proposed as an alternative delivery route of pharmacological agents via the bilaminar sac of pericardium surrounding the heart. To date, intrapericardial administration has entailed the localized administration of a broad spectrum of therapeutic agents. These agents include stem cells, extracellular matrix, growth factor, drugs, bioactive materials, and genetic materials, to the heart and coronary arteries. The route not only overcomes the limitations associated with traditional systemic administration methods, but also presents multiple intrinsic advantages over the other approaches, allowing greater therapeutic actions. Intrapericardial administration exhibits versatility in addressing certain cardiac conditions and ongoing research in this field certainly holds promise for further innovations and advancements to improve cardiac treatment. Thus, this review discusses the anatomy and physiology of the pericardium, the intrapericardial administration access routes, the recent application of intrapericardial delivery in the context of cardiac repair as well as the challenges associated with the approach.

Respiratory motion-corrected T1 mapping of the abdomen

OBJECTIVE

The purpose of this study was to investigate an approach for motion-corrected T1 mapping of the abdomen that allows for free breathing data acquisition with 100% scan efficiency.

MATERIALS AND METHODS

Data were acquired using a continuous golden radial trajectory and multiple inversion pulses. For the correction of respiratory motion, motion estimation based on a surrogate was performed from the same data used for T1 mapping. Image-based self-navigation allowed for binning and reconstruction of respiratory-resolved images, which were used for the estimation of respiratory motion fields. Finally, motion-corrected T1 maps were calculated from the data applying the estimated motion fields. The method was evaluated in five healthy volunteers. For the assessment of the image-based navigator, we compared it to a simultaneously acquired ultrawide band radar signal. Motion-corrected T1 maps were evaluated qualitatively and quantitatively for different scan times.

RESULTS

For all volunteers, the motion-corrected T1 maps showed fewer motion artifacts in the liver as well as sharper kidney structures and blood vessels compared to uncorrected T1 maps. Moreover, the relative error to the reference breathhold T1 maps could be reduced from up to 25% for the uncorrected T1 maps to below 10% for the motion-corrected maps for the average value of a region of interest, while the scan time could be reduced to 6-8 s.

DISCUSSION

The proposed approach allows for respiratory motion-corrected T1 mapping in the abdomen and ensures accurate T1 maps without the need for any breathholds.

Large-scale 3D non-Cartesian coronary MRI reconstruction using distributed memory-efficient physics-guided deep learning with limited training data

OBJECT

To enable high-quality physics-guided deep learning (PG-DL) reconstruction of large-scale 3D non-Cartesian coronary MRI by overcoming challenges of hardware limitations and limited training data availability.

MATERIALS AND METHODS

While PG-DL has emerged as a powerful image reconstruction method, its application to large-scale 3D non-Cartesian MRI is hindered by hardware limitations and limited availability of training data. We combine several recent advances in deep learning and MRI reconstruction to tackle the former challenge, and we further propose a 2.5D reconstruction using 2D convolutional neural networks, which treat 3D volumes as batches of 2D images to train the network with a limited amount of training data. Both 3D and 2.5D variants of the PG-DL networks were compared to conventional methods for high-resolution 3D kooshball coronary MRI.

RESULTS

Proposed PG-DL reconstructions of 3D non-Cartesian coronary MRI with 3D and 2.5D processing outperformed all conventional methods both quantitatively and qualitatively in terms of image assessment by an experienced cardiologist. The 2.5D variant further improved vessel sharpness compared to 3D processing, and scored higher in terms of qualitative image quality.

DISCUSSION

PG-DL reconstruction of large-scale 3D non-Cartesian MRI without compromising image size or network complexity is achieved, and the proposed 2.5D processing enables high-quality reconstruction with limited training data.

Motion robust coronary MR angiography using zigzag centric ky-kz trajectory and high-resolution deep learning reconstruction

PURPOSE

To develop a new MR coronary angiography (MRCA) technique by employing a zigzag fan-shaped centric ky-kz k-space trajectory combined with high-resolution deep learning reconstruction (HR-DLR).

METHODS

All imaging data were acquired from 12 healthy subjects and 2 patients using two clinical 3-T MR imagers, with institutional review board approval. Ten healthy subjects underwent both standard 3D fast gradient echo (sFGE) and centric ky-kz k-space trajectory FGE (cFGE) acquisitions to compare the scan time and image quality. Quantitative measures were also performed for signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) as well as sharpness of the vessel. Furthermore, the feasibility of the proposed cFGE sequence was assessed in two patients. For assessing the feasibility of the centric ky-kz trajectory, the navigator-echo window of a 30-mm threshold was applied in cFGE, whereas sFGE was applied using a standard 5-mm threshold. Image quality of MRCA using cFGE with HR-DLR and sFGE without HR-DLR was scored in a 5-point scale (non-diagnostic = 1, fair = 2, moderate = 3, good = 4, and excellent = 5). Image evaluation of cFGE, applying HR-DLR, was compared with sFGE without HR-DLR. Friedman test, Wilcoxon signed-rank test, or paired t tests were performed for the comparison of related variables.

RESULTS

The actual MRCA scan time of cFGE with a 30-mm threshold was acquired in less than 5 min, achieving nearly 100% efficiency, showcasing its expeditious and robustness. In contrast, sFGE was acquired with a 5-mm threshold and had an average scan time of approximately 15 min. Overall image quality for MRCA was scored 3.3 for sFGE and 2.7 for cFGE without HR-DLR but increased to 3.6 for cFGE with HR-DLR and (p < 0.05). The clinical result of patients obtained within 5 min showed good quality images in both patients, even with a stent, without artifacts. Quantitative measures of SNR, CNR, and sharpness of vessel presented higher in cFGE with HR-DLR.

CONCLUSION

Our findings demonstrate a robust, time-efficient solution for high-quality MRCA, enhancing patient comfort and increasing clinical throughput.

Similarity-driven motion-resolved reconstruction for ferumoxytol-enhanced whole-heart MRI in congenital heart disease

A similarity-driven multi-dimensional binning algorithm (SIMBA) reconstruction of free-running cardiac magnetic resonance imaging data was previously proposed. While very efficient and fast, the original SIMBA focused only on the reconstruction of a single motion-consistent cluster, discarding the remaining data acquired. However, the redundant data clustered by similarity may be exploited to further improve image quality. In this work, we propose a novel compressed sensing (CS) reconstruction that performs an effective regularization over the clustering dimension, thanks to the integration of inter-cluster motion compensation (XD-MC-SIMBA). This reconstruction was applied to free-running ferumoxytol-enhanced datasets from 24 patients with congenital heart disease, and compared to the original SIMBA, the same XD-MC-SIMBA reconstruction but without motion compensation (XD-SIMBA), and a 5D motion-resolved CS reconstruction using the free-running framework (FRF). The resulting images were compared in terms of lung-liver and blood-myocardium sharpness, blood-myocardium contrast ratio, and visible length and sharpness of the coronary arteries. Moreover, an automated image quality score (IQS) was assigned using a pretrained deep neural network. The lung-liver sharpness and blood-myocardium sharpness were significantly higher in XD-MC-SIMBA and FRF. Consistent with these findings, the IQS analysis revealed that image quality for XD-MC-SIMBA was improved in 18 of 24 cases, compared to SIMBA. We successfully tested the hypothesis that multiple motion-consistent SIMBA clusters can be exploited to improve the quality of ferumoxytol-enhanced cardiac MRI when inter-cluster motion-compensation is integrated as part of a CS reconstruction.

Trajectory correction enables free-running chemical shift encoded imaging for accurate cardiac proton-density fat fraction quantification at 3T

BACKGROUND

Metabolic diseases can negatively alter epicardial fat accumulation and composition, which can be probed using quantitative cardiac chemical shift encoded (CSE) cardiovascular magnetic resonance (CMR) by mapping proton-density fat fraction (PDFF). To obtain motion-resolved high-resolution PDFF maps, we proposed a free-running cardiac CSE-CMR framework at 3T. To employ faster bipolar readout gradients, a correction for gradient imperfections was added using the gradient impulse response function (GIRF) and evaluated on intermediate images and PDFF quantification.

METHODS

Ten minutes free-running cardiac 3D radial CSE-CMR acquisitions were compared in vitro and in vivo at 3T. Monopolar and bipolar readout gradient schemes provided 8 echoes (TE1/ΔTE = 1.16/1.96 ms) and 13 echoes (TE1/ΔTE = 1.12/1.07 ms), respectively. Bipolar-gradient free-running cardiac fat and water images and PDFF maps were reconstructed with or without GIRF correction. PDFF values were evaluated in silico, in vitro on a fat/water phantom, and in vivo in 10 healthy volunteers and 3 diabetic patients.

RESULTS

In monopolar mode, fat-water swaps were demonstrated in silico and confirmed in vitro. Using bipolar readout gradients, PDFF quantification was reliable and accurate with GIRF correction with a mean bias of 0.03% in silico and 0.36% in vitro while it suffered from artifacts without correction, leading to a PDFF bias of 4.9% in vitro and swaps in vivo. Using bipolar readout gradients, in vivo PDFF of epicardial adipose tissue was significantly lower compared to subcutaneous fat (80.4 ± 7.1% vs 92.5 ± 4.3%, P < 0.0001).

CONCLUSIONS

Aiming for an accurate PDFF quantification, high-resolution free-running cardiac CSE-MRI imaging proved to benefit from bipolar echoes with k-space trajectory correction at 3T. This free-breathing acquisition framework enables to investigate epicardial adipose tissue PDFF in metabolic diseases.

Optimal Protocol for Contrast-enhanced Free-running 5D Whole-heart Coronary MR Angiography at 3T

Free-running 5D whole-heart coronary MR angiography (MRA) is gaining in popularity because it reduces scanning complexity by removing the need for specific slice orientations, respiratory gating, or cardiac triggering. At 3T, a gradient echo (GRE) sequence is preferred in combination with contrast injection. However, neither the injection scheme of the gadolinium (Gd) contrast medium, the choice of the RF excitation angle, nor the dedicated image reconstruction parameters have been established for 3T GRE free-running 5D whole-heart coronary MRA. In this study, a Gd injection scheme, RF excitation angles of lipid-insensitive binominal off-resonance RF excitation (LIBRE) pulse for valid fat suppression and continuous data acquisition, and compressed-sensing reconstruction regularization parameters were optimized for contrast-enhanced free-running 5D whole-heart coronary MRA using a GRE sequence at 3T. Using this optimized protocol, contrast-enhanced free-running 5D whole-heart coronary MRA using a GRE sequence is feasible with good image quality at 3T.

Intra-bin correction and inter-bin compensation of respiratory motion in free-running five-dimensional whole-heart magnetic resonance imaging

BACKGROUND

Free-running cardiac and respiratory motion-resolved whole-heart five-dimensional (5D) cardiovascular magnetic resonance (CMR) can reduce scan planning and provide a means of evaluating respiratory-driven changes in clinical parameters of interest. However, respiratory-resolved imaging can be limited by user-defined parameters which create trade-offs between residual artifact and motion blur. In this work, we develop and validate strategies for both correction of intra-bin and compensation of inter-bin respiratory motion to improve the quality of 5D CMR.

METHODS

Each component of the reconstruction framework was systematically validated and compared to the previously established 5D approach using simulated free-running data (N = 50) and a cohort of 32 patients with congenital heart disease. The impact of intra-bin respiratory motion correction was evaluated in terms of image sharpness while inter-bin respiratory motion compensation was evaluated in terms of reconstruction error, compression of respiratory motion, and image sharpness. The full reconstruction framework (intra-acquisition correction and inter-acquisition compensation of respiratory motion [IIMC] 5D) was evaluated in terms of image sharpness and scoring of image quality by expert reviewers.

RESULTS

Intra-bin motion correction provides significantly (p < 0.001) sharper images for both simulated and patient data. Inter-bin motion compensation results in significant (p < 0.001) lower reconstruction error, lower motion compression, and higher sharpness in both simulated (10/11) and patient (9/11) data. The combined framework resulted in significantly (p < 0.001) sharper IIMC 5D reconstructions (End-expiration (End-Exp): 0.45 ± 0.09, End-inspiration (End-Ins): 0.46 ± 0.10) relative to the previously established 5D implementation (End-Exp: 0.43 ± 0.08, End-Ins: 0.39 ± 0.09). Similarly, image scoring by three expert reviewers was significantly (p < 0.001) higher using IIMC 5D (End-Exp: 3.39 ± 0.44, End-Ins: 3.32 ± 0.45) relative to 5D images (End-Exp: 3.02 ± 0.54, End-Ins: 2.45 ± 0.52).

CONCLUSION

The proposed IIMC reconstruction significantly improves the quality of 5D whole-heart MRI. This may be exploited for higher resolution or abbreviated scanning. Further investigation of the diagnostic impact of this framework and comparison to gold standards is needed to understand its full clinical utility, including exploration of respiratory-driven changes in physiological measurements of interest.

Multi-Device Parallel MRI Reconstruction: Efficient Partitioning for Undersampled 5D Cardiac CINE

Cardiac CINE, a form of dynamic cardiac MRI, is indispensable in the diagnosis and treatment of heart conditions, offering detailed visualization essential for the early detection of cardiac diseases. As the demand for higher-resolution images increases, so does the volume of data requiring processing, presenting significant computational challenges that can impede the efficiency of diagnostic imaging. Our research presents an approach that takes advantage of the computational power of multiple Graphics Processing Units (GPUs) to address these challenges. GPUs are devices capable of performing large volumes of computations in a short period, and have significantly improved the cardiac MRI reconstruction process, allowing images to be produced faster. The innovation of our work resides in utilizing a multi-device system capable of processing the substantial data volumes demanded by high-resolution, five-dimensional cardiac MRI. This system surpasses the memory capacity limitations of single GPUs by partitioning large datasets into smaller, manageable segments for parallel processing, thereby preserving image integrity and accelerating reconstruction times. Utilizing OpenCL technology, our system offers adaptability and cross-platform functionality, ensuring wider applicability. The proposed multi-device approach offers an advancement in medical imaging, accelerating the reconstruction process and facilitating faster and more effective cardiac health assessment.

Combined free-running four-dimensional anatomical and flow magnetic resonance imaging with native contrast using Synchronization of Neighboring Acquisitions by Physiological Signals

BACKGROUND

Four-dimensional (4D) flow magnetic resonance imaging (MRI) often relies on the injection of gadolinium- or iron-oxide-based contrast agents to improve vessel delineation. In this work, a novel technique is developed to acquire and reconstruct 4D flow data with excellent dynamic visualization of blood vessels but without the need for contrast injection. Synchronization of Neighboring Acquisitions by Physiological Signals (SyNAPS) uses pilot tone (PT) navigation to retrospectively synchronize the reconstruction of two free-running three-dimensional radial acquisitions, to create co-registered anatomy and flow images.

METHODS

Thirteen volunteers and two Marfan syndrome patients were scanned without contrast agent using one free-running fast interrupted steady-state (FISS) sequence and one free-running phase-contrast MRI (PC-MRI) sequence. PT signals spanning the two sequences were recorded for retrospective respiratory motion correction and cardiac binning. The magnitude and phase images reconstructed, respectively, from FISS and PC-MRI, were synchronized to create SyNAPS 4D flow datasets. Conventional two-dimensional (2D) flow data were acquired for reference in ascending (AAo) and descending aorta (DAo). The blood-to-myocardium contrast ratio, dynamic vessel area, net volume, and peak flow were used to compare SyNAPS 4D flow with Native 4D flow (without FISS information) and 2D flow. A score of 0-4 was given to each dataset by two blinded experts regarding the feasibility of performing vessel delineation.

RESULTS

Blood-to-myocardium contrast ratio for SyNAPS 4D flow magnitude images (1.5 ± 0.3) was significantly higher than for Native 4D flow (0.7 ± 0.1, p < 0.01) and was comparable to 2D flow (2.3 ± 0.9, p = 0.02). Image quality scores of SyNAPS 4D flow from the experts (M.P.: 1.9 ± 0.3, E.T.: 2.5 ± 0.5) were overall significantly higher than the scores from Native 4D flow (M.P.: 1.6 ± 0.6, p = 0.03, E.T.: 0.8 ± 0.4, p < 0.01) but still significantly lower than the scores from the reference 2D flow datasets (M.P.: 2.8 ± 0.4, p < 0.01, E.T.: 3.5 ± 0.7, p < 0.01). The Pearson correlation coefficient between the dynamic vessel area measured on SyNAPS 4D flow and that from 2D flow was 0.69 ± 0.24 for the AAo and 0.83 ± 0.10 for the DAo, whereas the Pearson correlation between Native 4D flow and 2D flow measurements was 0.12 ± 0.48 for the AAo and 0.08 ± 0.39 for the DAo. Linear correlations between SyNAPS 4D flow and 2D flow measurements of net volume (r2 = 0.83) and peak flow (r2 = 0.87) were larger than the correlations between Native 4D flow and 2D flow measurements of net volume (r2 = 0.79) and peak flow (r2 = 0.76).

CONCLUSION

The feasibility and utility of SyNAPS were demonstrated for joint whole-heart anatomical and flow MRI without requiring electrocardiography gating, respiratory navigators, or contrast agents. Using SyNAPS, a high-contrast anatomical imaging sequence can be used to improve 4D flow measurements that often suffer from poor delineation of vessel boundaries in the absence of contrast agents.

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.

Reproducibility of 3D thoracic aortic displacement from 3D cine balanced SSFP at 3 T without contrast enhancement

PURPOSE

Aortic motion has direct impact on the mechanical stresses acting on the aorta. In aortic disease, increased stiffness of the aorta may lead to decreased aortic motion over time, which could be a predictor for aortic dissection or rupture. This study investigates the reproducibility of obtaining 3D displacement and diameter maps quantified using accelerated 3D cine MRI at 3 T.

METHODS

A noncontrast-enhanced, free-breathing 3D cine sequence based on balanced SSFP and pseudo-spiral undersampling with high spatial isotropic resolution was developed (spatial/temporal resolution [1.6 mm]3 /67 ms). The thoracic aorta of 14 healthy volunteers was prospectively scanned three times at 3 T: twice on the same day and a third time 2 weeks later. Aortic displacement was calculated using iterative closest point nonrigid registration of manual segmentations of the 3D aorta at end-systole and mid-diastole. Interexamination and interobserver regional analysis of mean displacement for five regions of interest was performed using Bland-Altman analysis. Additionally, a complementary voxel-by-voxel analysis was done, allowing a more local inspection of the method.

RESULTS

No significant differences were found in mean and maximum displacement for any of the regions of interest for the interexamination and interobserver analysis. The maximum displacement measured in the lower half of the ascending aorta was 11.0 ± 3.4 mm (range: 3.0-17.5 mm) for the first scan. The smallest detectable change in mean displacement in the lower half of the ascending aorta was 3 mm.

CONCLUSION

Detailed 3D cine balanced SSFP at 3 T allows for reproducible quantification of systolic-diastolic mean aortic displacement within acceptable limits.

On the Determination of Lagrange Multipliers for a Weighted LASSO Problem Using Geometric and Convex Analysis Techniques

Compressed Sensing (CS) encompasses a broad array of theoretical and applied techniques for recovering signals, given partial knowledge of their coefficients, cf. Candés (C. R. Acad. Sci. Paris, Ser. I 346, 589-592 (2008)), Candés et al. (IEEE Trans. Inf. Theo (2006)), Donoho (IEEE Trans. Inf. Theo. 52(4), (2006)), Donoho et al. (IEEE Trans. Inf. Theo. 52(1), (2006)). Its applications span various fields, including mathematics, physics, engineering, and several medical sciences, cf. Adcock and Hansen (Compressive Imaging: Structure, Sampling, Learning, p. 2021), Berk et al. (2019 13th International conference on Sampling Theory and Applications (SampTA) pp. 1-5. IEEE (2019)), Brady et al. (Opt. Express 17(15), 13040-13049 (2009)), Chan (Terahertz imaging with compressive sensing. Rice University, USA (2010)), Correa et al. (2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP) pp. 7789-7793 (2014, May) IEEE), Gao et al. (Nature 516(7529), 74-77 (2014)), Liu and Kang (Opt. Express 18(21), 22010-22019 (2010)), McEwen and Wiaux (Mon. Notices Royal Astron. Soc. 413(2), 1318-1332 (2011)), Marim et al. (Opt. Lett. 35(6), 871-873 (2010)), Yu and Wang (Phys. Med. Biol. 54(9), 2791 (2009)), Yu and Wang (Phys. Med. Biol. 54(9), 2791 (2009)). Motivated by our interest in the mathematics behind Magnetic Resonance Imaging (MRI) and CS, we employ convex analysis techniques to analytically determine equivalents of Lagrange multipliers for optimization problems with inequality constraints, specifically a weighted LASSO with voxel-wise weighting. We investigate this problem under assumptions on the fidelity term A x - b 2 2 , either concerning the sign of its gradient or orthogonality-like conditions of its matrix. To be more precise, we either require the sign of each coordinate of 2 ( A x - b ) T A to be fixed within a rectangular neighborhood of the origin, with the side lengths of the rectangle dependent on the constraints, or we assume A T A to be diagonal. The objective of this work is to explore the relationship between Lagrange multipliers and the constraints of a weighted variant of LASSO, specifically in the mentioned cases where this relationship can be computed explicitly. As they scale the regularization terms of the weighted LASSO, Lagrange multipliers serve as tuning parameters for the weighted LASSO, prompting the question of their potential effective use as tuning parameters in applications like MR image reconstruction and denoising. This work represents an initial step in this direction.

A characterization of cardiac-induced noise in R2 * maps of the brain

PURPOSE

Cardiac pulsation increases the noise level in brain maps of the transverse relaxation rate R2 *. Cardiac-induced noise is challenging to mitigate during the acquisition of R2 * mapping data because its characteristics are unknown. In this work, we aim to characterize cardiac-induced noise in brain maps of the MRI parameter R2 *.

METHODS

We designed a sampling strategy to acquire multi-echo 3D data in 12 intervals of the cardiac cycle, monitored with a fingertip pulse-oximeter. We measured the amplitude of cardiac-induced noise in this data and assessed the effect of cardiac pulsation on R2 * maps computed across echoes. The area of k-space that contains most of the cardiac-induced noise in R2 * maps was then identified. Based on these characteristics, we introduced a tentative sampling strategy that aims to mitigate cardiac-induced noise in R2 * maps of the brain.

RESULTS

In inferior brain regions, cardiac pulsation accounts for R2 * variations of up to 3 s-1 across the cardiac cycle (i.e., ∼35% of the overall variability). Cardiac-induced fluctuations occur throughout the cardiac cycle, with a reduced intensity during the first quarter of the cycle. A total of 50% to 60% of the overall cardiac-induced noise is localized near the k-space center (k < 0.074 mm-1 ). The tentative cardiac noise mitigation strategy reduced the variability of R2 * maps across repetitions by 11% in the brainstem and 6% across the whole brain.

CONCLUSION

We provide a characterization of cardiac-induced noise in brain R2 * maps that can be used as a basis for the design of mitigation strategies during data acquisition.

Dynamic cardiac MRI with high spatiotemporal resolution using accelerated spiral-out and spiral-in/out bSSFP pulse sequences at 1.5 T

OBJECTIVE

To develop two spiral-based bSSFP pulse sequences combined with L + S reconstruction for accelerated ungated, free-breathing dynamic cardiac imaging at 1.5 T.

MATERIALS AND METHODS

Tiny golden angle rotated spiral-out and spiral-in/out bSSFP sequences combined with view-sharing (VS), compressed sensing (CS), and low-rank plus sparse (L + S) reconstruction were evaluated and compared via simulation and in vivo dynamic cardiac imaging studies. The proposed methods were then validated against the standard cine, in terms of quantitative image assessment and qualitative quality rating.

RESULTS

The L + S method yielded the least residual artifacts and the best image sharpness among the three methods. Both spiral cine techniques showed clinically diagnostic images (score > 3). Compared to standard cine, there were significant differences in global image quality and edge sharpness for spiral cine techniques, while there was significant difference in image contrast for the spiral-out cine but no significant difference for the spiral-in/out cine. There was good agreement in left ventricular ejection fraction for both the spiral-out cine (- 1.6 [Formula: see text] 3.1%) and spiral-in/out cine (- 1.5 [Formula: see text] 2.8%) against standard cine.

DISCUSSION

Compared to the time-consuming standard cine (~ 5 min) which requires ECG-gating and breath-holds, the proposed spiral bSSFP sequences achieved ungated, free-breathing cardiac movies at a similar spatial (1.5 × 1.5 × 8 mm3) and temporal resolution (36 ms) per slice for whole heart coverage (10-15 slices) within 45 s, suggesting the clinical potential for improved patient comfort or for imaging patients with arrhythmias or who cannot hold their breath.

Motion Compensated Unsupervised Deep Learning for 5D MRI

We propose an unsupervised deep learning algorithm for the motion-compensated reconstruction of 5D cardiac MRI data from 3D radial acquisitions. Ungated free-breathing 5D MRI simplifies the scan planning, improves patient comfort, and offers several clinical benefits over breath-held 2D exams, including isotropic spatial resolution and the ability to reslice the data to arbitrary views. However, the current reconstruction algorithms for 5D MRI take very long computational time, and their outcome is greatly dependent on the uniformity of the binning of the acquired data into different physiological phases. The proposed algorithm is a more data-efficient alternative to current motion-resolved reconstructions. This motion-compensated approach models the data in each cardiac/respiratory bin as Fourier samples of the deformed version of a 3D image template. The deformation maps are modeled by a convolutional neural network driven by the physiological phase information. The deformation maps and the template are then jointly estimated from the measured data. The cardiac and respiratory phases are estimated from 1D navigators using an auto-encoder. The proposed algorithm is validated on 5D bSSFP datasets acquired from two subjects.

Three-dimensional sector-wise golden angle-improved k-space uniformity after electrocardiogram binning

PURPOSE

To develop and evaluate a 3D sector-wise golden-angle (3D-SWIG) profile ordering scheme for cardiovascular MR cine imaging that maintains high k-space uniformity after electrocardiogram (ECG) binning.

METHOD

Cardiovascular MR (CMR) was performed at 1.5 T. A balanced SSFP pulse sequence was implemented with a novel 3D-SWIG radial ordering, where k-space was divided into wedges, and each wedge was acquired in a separate heartbeat. The high uniformity of k-space coverage after physiological binning can be used to perform functional imaging using a very short acquisition. The 3D-SWIG was compared with two commonly used 3D radial trajectories for CMR (i.e., double golden angle and spiral phyllotaxis) in numerical simulations. Free-breathing 3D-SWIG and conventional breath-held 2D cine were compared in patients (n = 17) referred clinically for CMR. Quantitative comparison was performed based on left ventricular segmentation.

RESULTS

Numerical simulations showed that 3D-SWIG both required smaller steps between successive readouts and achieved better k-space sampling uniformity after binning than either the double golden angle or spiral phyllotaxis trajectories. In vivo evaluation showed that measurements of left ventricular ejection fraction calculated from a 48 heart-beat free-breathing 3D-SWIG acquisition were highly reproducible and agreed with breath-held 2D-Cartesian cine (mean ± SD difference of -3.1 ± 3.5% points).

CONCLUSIONS

The 3D-SWIG acquisition offers a simple solution for highly improved k-space uniformity after physiological binning. The feasibility of the 3D-SWIG method is demonstrated in this study through whole-heart cine imaging during free breathing with an acquisition time of less than 1 min.

Retrospective motion correction through multi-average k-space data elimination (REMAKE) for free-breathing cardiac cine imaging

PURPOSE

To develop a motion-robust reconstruction technique for free-breathing cine imaging with multiple averages.

METHOD

Retrospective motion correction through multiple average k-space data elimination (REMAKE) was developed using iterative removal of k-space segments (from individual k-space samples) that contribute most to motion corruption while combining any remaining segments across multiple signal averages. A variant of REMAKE, termed REMAKE+, was developed to address any losses in SNR due to k-space information removal. With REMAKE+, multiple reconstructions using different initial conditions were performed, co-registered, and averaged. Both techniques were validated against clinical "standard" signal averaging reconstruction in a static phantom (with simulated motion) and 15 patients undergoing free-breathing cine imaging with multiple averages. Quantitative analysis of myocardial sharpness, blood/myocardial SNR, myocardial-blood contrast-to-noise ratio (CNR), as well as subjective assessment of image quality and rate of diagnostic quality images were performed.

RESULTS

In phantom, motion artifacts using "standard" (RMS error [RMSE]: 2.2 ± 0.5) were substantially reduced using REMAKE/REMAKE+ (RMSE: 1.5 ± 0.4/1.0 ± 0.4, p < 0.01). In patients, REMAKE/REMAKE+ led to higher myocardial sharpness (0.79 ± 0.09/0.79 ± 0.1 vs. 0.74 ± 0.12 for "standard", p = 0.004/0.04), higher image quality (1.8 ± 0.2/1.9 ± 0.2 vs. 1.6 ± 0.4 for "standard", p = 0.02/0.008), and a higher rate of diagnostic quality images (99%/100% vs. 94% for "standard"). Blood/myocardial SNR for "standard" (94 ± 30/33 ± 10) was higher vs. REMAKE (80 ± 25/28 ± 8, p = 0.002/0.005) and tended to be lower vs. REMAKE+ (105 ± 33/36 ± 12, p = 0.02/0.06). Myocardial-blood CNR for "standard" (61 ± 22) was higher vs. REMAKE (53 ± 19, p = 0.003) and lower vs. REMAKE+ (69 ± 24, p = 0.007).

CONCLUSIONS

Compared to "standard" signal averaging reconstruction, REMAKE and REMAKE+ provide improved myocardial sharpness, image quality, and rate of diagnostic quality images.

High-resolution, respiratory-resolved coronary MRA using a Phyllotaxis-reordered variable-density 3D cones trajectory

PURPOSE

To develop a respiratory-resolved motion-compensation method for free-breathing, high-resolution coronary magnetic resonance angiography (CMRA) using a 3D cones trajectory.

METHODS

To achieve respiratory-resolved 0.98 mm resolution images in a clinically relevant scan time, we undersample the imaging data with a variable-density 3D cones trajectory. For retrospective motion compensation, translational estimates from 3D image-based navigators (3D iNAVs) are used to bin the imaging data into four phases from end-expiration to end-inspiration. To ensure pseudo-random undersampling within each respiratory phase, we devise a phyllotaxis readout ordering scheme mindful of eddy current artifacts in steady state free precession imaging. Following binning, residual 3D translational motion within each phase is computed using the 3D iNAVs and corrected for in the imaging data. The noise-like aliasing characteristic of the combined phyllotaxis and cones sampling pattern is leveraged in a compressed sensing reconstruction with spatial and temporal regularization to reduce aliasing in each of the respiratory phases.

RESULTS

In initial studies of six subjects, respiratory motion compensation using the proposed method yields improved image quality compared to non-respiratory-resolved approaches with no motion correction and with 3D translational correction. Qualitative assessment by two cardiologists and quantitative evaluation with the image edge profile acutance metric indicate the superior sharpness of coronary segments reconstructed with the proposed method (P < 0.01).

CONCLUSION

We have demonstrated a new method for free-breathing, high-resolution CMRA based on a variable-density 3D cones trajectory with modified phyllotaxis ordering and respiratory-resolved motion compensation with 3D iNAVs.

Free-breathing half-radial dual-echo balanced steady-state free precession thoracic imaging with wobbling Archimedean spiral pole trajectories

PURPOSE

To demonstrate free-breathing thoracic MRI with a minimal-TR balanced steady-state free precession (bSSFP) technique using wobbling Archimedean spiral pole (WASP) trajectories.

METHODS

Phantom and free-breathing in vivo chest imaging in healthy volunteers was performed at 1.5T with a half-radial, dual-echo, bSSFP sequence, termed bSTAR. For maximum sampling efficiency, a single analog-to-digital converter window along the full bipolar readout was used. To ensure a homogeneous coverage of the k-space over multiple breathing cycles, radial k-space sampling followed short-duration Archimedean spiral interleaves that were randomly titled by a small polar angle and rotated by a golden angle about the polar axis; depticting a wobbling Archimedean spiral pole (WASP) trajectory. In phantom and in vivo experiments, WASP trajectories were compared to spiral phyllotaxis sampling in terms of eddy currents and were used to generate in vivo thorax images at different respiratory phases.

RESULTS

WASP trajectories provided artifact-free bSTAR imaging in both phantom and in vivo and respiratory self-gated reconstruction was successfully performed in all subjects. The amount of the acquired data allowed the reconstruction of 10 volumes at different respiratory levels with isotropic resolution of 1.77mm from a scan of 5.5minutes (using a TR of 1.32ms), and one high-resolution 1.16mm end-expiratory volume from a scan of 4.7minutes (using a TR of 1.42ms). The very short TR of bSTAR mitigated off-resonance artifacts despite the large field-of-view.

CONCLUSION

We have demonstrated the feasibility of high-resolution free-breathing thoracic imaging with bSTAR using the wobbling Archimedean spiral pole in healthy subjects at 1.5T.

MR Multitasking-based multi-dimensional assessment of cardiovascular system (MT-MACS) with extended spatial coverage and water-fat separation

PURPOSE

To extend the MR MultiTasking-based Multidimensional Assessment of Cardiovascular System (MT-MACS) technique with larger spatial coverage and water-fat separation for comprehensive aortocardiac assessment.

METHODS

MT-MACS adopts a low-rank tensor image model for 7D imaging, with three spatial dimensions for volumetric imaging, one cardiac motion dimension for cine imaging, one respiratory motion dimension for free-breathing imaging, one T2-prepared inversion recovery time dimension for multi-contrast assessment, and one T2*-decay time dimension for water-fat separation. Nine healthy subjects were recruited for the 3T study. Overall image quality was scored on bright-blood (BB), dark-blood (DB), and gray-blood (GB) contrasts using a 4-point scale (0-poor to 3-excellent) by two independent readers, and their interreader agreement was evaluated. Myocardial wall thickness and left ventricular ejection fraction (LVEF) were quantified on DB and BB contrasts, respectively. The agreement in these metrics between MT-MACS and conventional breath-held, electrocardiography-triggered 2D sequences were evaluated.

RESULTS

MT-MACS provides both water-only and fat-only images with excellent image quality (average score = 3.725/3.780/3.835/3.890 for BB/DB/GB/fat-only images) and moderate to high interreader agreement (weighted Cohen's kappa value = 0.727/0.668/1.000/1.000 for BB/DB/GB/fat-only images). There were good to excellent agreements in myocardial wall thickness measurements (intraclass correlation coefficients [ICC] = 0.781/0.929/0.680/0.878 for left atria/left ventricle/right atria/right ventricle) and LVEF quantification (ICC = 0.716) between MT-MACS and 2D references. All measurements were within the literature range of healthy subjects.

CONCLUSION

The refined MT-MACS technique provides multi-contrast, phase-resolved, and water-fat imaging of the aortocardiac systems and allows evaluation of anatomy and function. Clinical validation is warranted.

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.

Feasibility of 3.0 T balanced fast field echo non-contrast-enhanced whole-heart coronary magnetic resonance angiography

Background

Coronary artery disease (CAD) is one of the most common diseases seriously harmful to human health caused by atherosclerosis. Besides coronary computed tomography angiography (CCTA) and invasive coronary angiography (ICA), coronary magnetic resonance angiography (CMRA) has become an alternative examination. The purpose of this study was to prospectively evaluate the feasibility of 3.0 T free-breathing whole-heart non-contrast-enhanced coronary magnetic resonance angiography (NCE-CMRA).

Methods

After Institutional Review Board approval, the NCE-CMRA data sets of 29 patients acquired successfully at 3.0 T were evaluated independently by two blinded readers for visualization and image quality of coronary arteries using the subjective quality grade. The acquisition times were recorded in the meantime. A part of the patients had undergone CCTA, we represented stenosis by scores and used the Kappa to evaluate the consistency between CCTA and NCE-CMRA.

Results

Six patients did not get diagnostic image quality because of severe artifacts. The image quality score assessed by both radiologists is 3.2±0.7, which means the NCE-CMRA can show the coronary arteries excellently. The main vessels of the coronary artery on NCE-CMRA images are considered reliably assessable. The acquisition time of NCE-CMRA, is 8.8±1.2 min. The Kappa of CCTA and NCE-CMRA on detecting stenosis is 0.842 (P<0.001).

Conclusions

The NCE-CMRA results in reliable image quality and visualization parameters of coronary arteries within a short scan time. The NCE-CMRA and CCTA have a good agreement for detecting stenosis.

Latest Advances in Image Acceleration: All Dimensions are Fair Game

Magnetic resonance imaging (MRI) is a versatile modality that can generate high-resolution images with a variety of tissue contrasts. However, MRI is a slow technique and requires long acquisition times, which increase with higher temporal and spatial resolution and/or when multiple contrasts and large volumetric coverage is required. In order to speedup MR data acquisition, several approaches have been introduced in the literature. Most of these techniques acquire less data than required and exploit intrinsic redundancies in the MR images to recover the information that was not sampled. This article presents a review of MR acquisition and reconstruction methods that have exploited redundancies in the temporal, spatial, and contrast/parametric dimensions to accelerate image data acquisition, focusing on cardiac and abdominal MR imaging applications. The review describes how each of these dimensions has been separately exploited for speeding up MR acquisition to then discuss more advanced techniques where multiple dimensions are exploited together for further reducing scan times. Finally, future directions for multidimensional image acceleration and remaining technical challenges are discussed. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: 1.

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.

Fetal 3D cardiovascular cine image acquisition using radial sampling and compressed sensing

PURPOSE

To explore a fetal 3D cardiovascular cine acquisition using a radial image acquisition and compressed-sensing reconstruction and compare image quality and scan time with conventional multislice 2D imaging.

METHODS

Volumetric fetal cardiac data were acquired in 26 volunteers using a radial 3D balanced SSFP pulse sequence. Cardiac gating was performed using a Doppler ultrasound device. Images were reconstructed using a parallel-imaging and compressed-sensing algorithm. Multiplanar reformatting to standard cardiac views was performed before image analysis. Clinical 2D images were used for comparison. Qualitative and quantitative image evaluation were performed by two experienced observers (scale: 1-4). Volumes, mass, and function were assessed.

RESULTS

Average scan time for the 3D imaging was 6 min, including one localizer. A 2D imaging stack covering the entire heart including localizer sequences took at least 6.5 min, depending on planning complexity. The 3D acquisition was successful in 7 of 26 subjects (27%). Overall image contrast and perceived resolution were lower in the 3D images. Nonetheless, the 3D images had, on average, a moderate cardiac diagnostic quality (median [range]: 3 [1-4]). Standard clinical 2D acquisitions had a high cardiac diagnostic quality (median [range]: 4 [3, 4]). Cardiac measurements were not different between 2D and 3D images (all p > 0.16).

CONCLUSION

The presented free-breathing whole-heart fetal 3D radial cine MRI acquisition and reconstruction method enables retrospective visualization of all cardiac views while keeping examination times short. This proof-of-concept work produced images with diagnostic quality, while at the same time reducing the planning complexity to a single localizer.

The mechanisms of arterial signal intensity profile in non-contrast coronary MRA (NC-MRCA): a 3D printed phantom investigation and clinical translations

Signal intensity (SI) drop has been proposed as an indirect stenosis assessment in non-contrast coronary MRA (NC-MRCA) but it uses unproven assumptions. We aimed to clarify the mechanisms that govern the SI in vitro and develop a stenosis detection method in vivo. Flow phantom tubes with/without stenosis were scanned under two spatial resolutions (0.5/1.0 mm³) on a 3.0 T MRI. Thirty-two coronary arteries from 11 volunteers were prospectively scanned with an EKG- and respiratory-gated 3D NC-MRCA with a resolution of 1.0 mm³, with coronary computed tomography angiography (CTA) as reference. The normalized SI along the centerline of the tubes or the coronary arteries was assessed against the distance from the orifice using a linear regression model. Its coefficient (SI decay slope) and goodness-of-fit (R2) were extracted to assess the effect of flow velocity and stenosis on the SI profile curve. The R2 was utilized for the stenosis detection. Phantom study: A slow flow velocity caused a steep SI decay slope. The SI drop revealed only at the inlet and outlet of stenosis due to the flow turbulence/vortex and yielded low R2, in which shape changed by the resolution. Clinical study: The R2 cutoff to detect ≥ 50% stenosis for the left and right coronary arteries were 0.64 and 0.20 with a sensitivity/specificity of 71.5/71.5 and 66.7/100 (%), respectively. The SI drop did not reflect the actual stenosis position and not suitable for the stenosis localization. The R2 cutoff represents an alternative method to detect stenoses on NC-MRCA at vessel level.Trial registration: ClinicalTrials.gov; NCT03768999, registered on December 7, 2018.

Simultaneous multi-parametric acquisition and reconstruction techniques in cardiac magnetic resonance imaging: Basic concepts and status of clinical development

Simultaneous multi-parametric acquisition and reconstruction techniques (SMART) are gaining attention for their potential to overcome some of cardiovascular magnetic resonance imaging's (CMR) clinical limitations. The major advantages of SMART lie within their ability to simultaneously capture multiple "features" such as cardiac motion, respiratory motion, T1/T2 relaxation. This review aims to summarize the overarching theory of SMART, describing key concepts that many of these techniques share to produce co-registered, high quality CMR images in less time and with less requirements for specialized personnel. Further, this review provides an overview of the recent developments in the field of SMART by describing how they work, the parameters they can acquire, their status of clinical testing and validation, and by providing examples for how their use can improve the current state of clinical CMR workflows. Many of the SMART are in early phases of development and testing, thus larger scale, controlled trials are needed to evaluate their use in clinical setting and with different cardiac pathologies.

Non-contrast free-breathing whole-heart 3D cine cardiovascular magnetic resonance with a novel 3D radial leaf trajectory

PURPOSE

To develop and validate a non-contrast free-breathing whole-heart 3D cine steady-state free precession (SSFP) sequence with a novel 3D radial leaf trajectory.

METHODS

We used a respiratory navigator to trigger acquisition of 3D cine data at end-expiration to minimize respiratory motion in our 3D cine SSFP sequence. We developed a novel 3D radial leaf trajectory to reduce gradient jumps and associated eddy-current artifacts. We then reconstructed the 3D cine images with a resolution of 2.0mm³ using an iterative nonlinear optimization algorithm. Prospective validation was performed by comparing ventricular volumetric measurements from a conventional breath-hold 2D cine ventricular short-axis stack against the non-contrast free-breathing whole-heart 3D cine dataset in each patient (n = 13).

RESULTS

All 3D cine SSFP acquisitions were successful and mean scan time was 07:09 ± 01:31 min. End-diastolic ventricular volumes for left ventricle (LV) and right ventricle (RV) measured from the 3D datasets were smaller than those from 2D (LV: 159.99 ± 42.99 vs. 173.16 ± 47.42; RV: 180.35 ± 46.08 vs. 193.13 ± 49.38; p-value≤0.044; bias<8%), whereas ventricular end-systolic volumes were more comparable (LV: 79.12 ± 26.78 vs. 78.46 ± 25.35; RV: 97.18 ± 32.35 vs. 102.42 ± 32.53; p-value≥0.190, bias<6%). The 3D cine data had a lower subjective image quality score.

CONCLUSION

Our non-contrast free-breathing whole-heart 3D cine sequence with novel leaf trajectory was robust and yielded smaller ventricular end-diastolic volumes compared to 2D cine imaging. It has the potential to make examinations easier and more comfortable for patients.

Ferumoxytol-enhanced 4D multiphase, steady-state imaging with magnetic resonance in congenital heart disease: ventricular volume and function across 2D and 3D software platforms

Background

Quantitative ventricular volumetry and function are important in the management of congenital heart disease (CHD). Ferumoxytol-enhanced (FE) 4D multiphase, steady state imaging with contrast enhancement (MUSIC) enables high-resolution, 3D cardiac phase-resolved magnetic resonance imaging (MRI) of the beating heart and extracardiac vessels in a single acquisition and without concerns about renal impairment. We aim to evaluate the semi-automatic quantification of ventricular volumetry and function of 4D MUSIC MRI using 2D and 3D software platforms.

Methods

This HIPAA-compliant and IRB-approved study prospectively recruited 50 children with CHD (3 days to 18 years) who underwent 4D MUSIC MRI at 3.0T between 2013-2017 for clinical indications. Each patient was either intubated in the neonatal intensive care unit (NICU) or underwent general anesthesia at MRI suite. For 2D analysis, we reformatted MUSIC images in Digital Imaging and Communications in Medicine (DICOM) format into ventricular short-axis slices with zero interslice gap. For 3D analysis, we imported DICOMs into a commercially available 3D software platform. Using semi-automatic thresholding, we quantified biventricular volume and ejection fraction (EF). We assessed the bias between MUSIC-derived 2D vs. 3D measurements and correlation between MUSIC vs. conventional 2D balanced steady-state free precession (bSSFP) cine images. We evaluated intra- and inter-observer agreement.

Results

There was a high degree of correlation between MUSIC-derived volumetric and functional measurements using 2D vs. 3D software (r=0.99, P<0.001). Volumes derived using 3D software platforms were larger than 2D by 0.2 to 2.0 mL/m2 whereas EF measurements were higher by 1.2-3.0%. MUSIC volumetric and functional measures derived from 2D and 3D software platforms corresponded highly with those derived from multi-slice SSFP cine images (r=0.99, P<0.001). The mean difference in volume for reformatted 4D MUSIC relative to bSSFP cine was 1.5 to 3.9 mL/m2. Intra- and inter-observer reliability was excellent.

Conclusions

Accurate and reliable ventricular volumetry and function can be derived from FE 4D MUSIC MRI studies using commercially available 2D and 3D software platforms. If fully validated in multicenter studies, the FE 4D-MUSIC pulse sequence may supercede conventional multislice 2D cine cardiovascular MRI acquisition protocols for functional evaluation of children with complex CHD.

Visualization and Analysis of Multidimensional Cardiovascular Magnetic Resonance Imaging: Challenges and Opportunities

Recent advances in magnetic resonance imaging are enabling the efficient creation of high-dimensional, multiparametric images, containing a wealth of potential information about the structure and function of many organs, including the cardiovascular system. However, the sizes of these rich data sets are so large that they are outstripping our ability to adequately visualize and analyze them, thus limiting their clinical impact. While there are some intrinsic limitations of human perception and of conventional display devices which hamper our ability to effectively use these data, newer computational methods for handling the data may aid our ability to extract and visualize the salient components of these high-dimensional data sets.

Cardio-respiratory motion-corrected 3D cardiac water-fat MRI using model-based image reconstruction

PURPOSE

Myocardial fat infiltrations are associated with a range of cardiomyopathies. The purpose of this study was to perform cardio-respiratory motion-correction for model-based water-fat separation to image fatty infiltrations of the heart in a free-breathing, non-cardiac-triggered high-resolution 3D MRI acquisition.

METHODS

Data were acquired in nine patients using a free-breathing, non-cardiac-triggered high-resolution 3D Dixon gradient-echo sequence and radial phase encoding trajectory. Motion correction was combined with a model-based water-fat reconstruction approach. Respiratory and cardiac motion models were estimated using a dual-mode registration algorithm incorporating both motion-resolved water and fat information. Qualitative comparisons of fat structures were made between 2D clinical routine reference scans and reformatted 3D motion-corrected images. To evaluate the effect of motion correction the local sharpness of epicardial fat structures was analyzed for motion-averaged and motion-corrected fat images.

RESULTS

The reformatted 3D motion-corrected reconstructions yielded qualitatively comparable fat structures and fat structure sharpness in the heart as the standard 2D breath-hold. Respiratory motion correction improved the local sharpness on average by 32% ± 24% with maximum improvements of 81% and cardiac motion correction increased the sharpness further by another 15% ± 11% with maximum increases of 31%. One patient showed a fat infiltration in the myocardium and cardio-respiratory motion correction was able to improve its visualization in 3D.

CONCLUSION

The 3D water-fat separated cardiac images were acquired during free-breathing and in a clinically feasible and predictable scan time. Compared to a motion-averaged reconstruction an increase in sharpness of fat structures by 51% ± 27% using the presented motion correction approach was observed for nine patients.

Free-running cardiac and respiratory motion-resolved 5D whole-heart coronary cardiovascular magnetic resonance angiography in pediatric cardiac patients using ferumoxytol

BACKGROUND

Coronary cardiovascular magnetic resonance angiography (CCMRA) of congenital heart disease (CHD) in pediatric patients requires accurate planning, adequate sequence parameter adjustments, lengthy scanning sessions, and significant involvement from highly trained personnel. Anesthesia and intubation are commonplace to minimize movements and control respiration in younger subjects. To address the above concerns and provide a single-click imaging solution, we applied our free-running framework for fully self-gated (SG) free-breathing 5D whole-heart CCMRA to CHD patients after ferumoxytol injection. We tested the hypothesis that spatial and motion resolution suffice to visualize coronary artery ostia in a cohort of CHD subjects, both for intubated and free-breathing acquisitions.

METHODS

In 18 pediatric CHD patients, non-electrocardiogram (ECG) triggered 5D free-running gradient echo CCMRA with whole-heart 1 mm3 isotropic spatial resolution was performed in seven minutes on a 1.5T CMR scanner. Eleven patients were anesthetized and intubated, while seven were breathing freely without anesthesia. All patients were slowly injected with ferumoxytol (4 mg/kg) over 15 minutes. Cardiac and respiratory motion-resolved 5D images were reconstructed with a fully SG approach. To evaluate the performance of motion resolution, visibility of coronary artery origins was assessed. Intubated and free-breathing patient sub-groups were compared for image quality using coronary artery length and conspicuity as well as lung-liver interface sharpness.

RESULTS

Data collection using the free-running framework was successful in all patients in less than 8 min; scan planning was very simple without the need for parameter adjustments, while no ECG lead placement and triggering was required. From the resulting SG 5D motion-resolved reconstructed images, coronary artery origins could be retrospectively extracted in 90% of the cases. These general findings applied to both intubated and free-breathing pediatric patients (no difference in terms of lung-liver interface sharpness), while image quality and coronary conspicuity between both cohorts was very similar.

CONCLUSIONS

A simple-to-use push-button framework for 5D whole-heart CCMRA was successfully employed in pediatric CHD patients with ferumoxytol injection. This approach, working without any external gating and for a wide range of heart rates and body sizes provided excellent definition of cardiac anatomy for both intubated and free-breathing patients.

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.

Cardiac MR: From Theory to Practice

Cardiovascular disease (CVD) is the leading single cause of morbidity and mortality, causing over 17. 9 million deaths worldwide per year with associated costs of over $800 billion. Improving prevention, diagnosis, and treatment of CVD is therefore a global priority. Cardiovascular magnetic resonance (CMR) has emerged as a clinically important technique for the assessment of cardiovascular anatomy, function, perfusion, and viability. However, diversity and complexity of imaging, reconstruction and analysis methods pose some limitations to the widespread use of CMR. Especially in view of recent developments in the field of machine learning that provide novel solutions to address existing problems, it is necessary to bridge the gap between the clinical and scientific communities. This review covers five essential aspects of CMR to provide a comprehensive overview ranging from CVDs to CMR pulse sequence design, acquisition protocols, motion handling, image reconstruction and quantitative analysis of the obtained data. (1) The basic MR physics of CMR is introduced. Basic pulse sequence building blocks that are commonly used in CMR imaging are presented. Sequences containing these building blocks are formed for parametric mapping and functional imaging techniques. Commonly perceived artifacts and potential countermeasures are discussed for these methods. (2) CMR methods for identifying CVDs are illustrated. Basic anatomy and functional processes are described to understand the cardiac pathologies and how they can be captured by CMR imaging. (3) The planning and conduct of a complete CMR exam which is targeted for the respective pathology is shown. Building blocks are illustrated to create an efficient and patient-centered workflow. Further strategies to cope with challenging patients are discussed. (4) Imaging acceleration and reconstruction techniques are presented that enable acquisition of spatial, temporal, and parametric dynamics of the cardiac cycle. The handling of respiratory and cardiac motion strategies as well as their integration into the reconstruction processes is showcased. (5) Recent advances on deep learning-based reconstructions for this purpose are summarized. Furthermore, an overview of novel deep learning image segmentation and analysis methods is provided with a focus on automatic, fast and reliable extraction of biomarkers and parameters of clinical relevance.

Motion-compensated fat-water imaging for 3D cardiac MRI at ultra-high fields

PURPOSE

Respiratory motion-compensated (MC) 3D cardiac fat-water imaging at 7T.

METHODS

Free-breathing bipolar 3D triple-echo gradient-recalled-echo (GRE) data with radial phase-encoding (RPE) trajectory were acquired in 11 healthy volunteers (7M\4F, 21-35 years, mean: 30 years) with a wide range of body mass index (BMI; 19.9-34.0 kg/m² ) and volunteer tailored B 1 + shimming. The bipolar-corrected triple-echo GRE-RPE data were binned into different respiratory phases (self-navigation) and were used for the estimation of non-rigid motion vector fields (MF) and respiratory resolved (RR) maps of the main magnetic field deviations (ΔB₀ ). RR ΔB₀ maps and MC ΔB₀ maps were compared to a reference respiratory phase to assess respiration-induced changes. Subsequently, cardiac binned fat-water images were obtained using a model-based, respiratory motion-corrected image reconstruction.

RESULTS

The 3D cardiac fat-water imaging at 7T was successfully demonstrated. Local respiration-induced frequency shifts in MC ΔB₀ maps are small compared to the chemical shifts used in the multi-peak model. Compared to the reference exhale ΔB₀ map these changes are in the order of 10 Hz on average. Cardiac binned MC fat-water reconstruction reduced respiration induced blurring in the fat-water images, and flow artifacts are reduced in the end-diastolic fat-water separated images.

CONCLUSION

This work demonstrates the feasibility of 3D fat-water imaging at UHF for the entire human heart despite spatial and temporal B 1 + and B₀ variations, as well as respiratory and cardiac motion.

Low-Rank and Framelet Based Sparsity Decomposition for Interventional MRI Reconstruction

OBJECTIVE

Interventional MRI (i-MRI) is crucial for MR image-guided therapy. Current image reconstruction methods for dynamic MR imaging are mostly retrospective that may not be suitable for i-MRI in real-time. Therefore, an algorithm to reconstruct images without a temporal pattern as in dynamic imaging is needed for i-MRI.

METHODS

We proposed a low-rank and sparsity (LS) decomposition algorithm with framelet transform to reconstruct the interventional feature with a high temporal resolution. Different from the existing LS based algorithm, the spatial sparsity of both the low-rank and sparsity components was used. We also used a primal dual fixed point (PDFP) method for optimization of the objective function to avoid solving sub-problems. Intervention experiments with gelatin and brain phantoms were carried out for validation.

RESULTS

The LS decomposition with framelet transform and PDFP could provide the best reconstruction performance compared with those without. Satisfying reconstruction results were obtained with only 10 radial spokes for a temporal resolution of 60 ms.

CONCLUSION AND SIGNIFICANCE

The proposed method has the potential for i-MRI in many different application scenarios.

Coronary Flow Assessment Using Accelerated 4D Flow MRI With Respiratory Motion Correction

Magnetic resonance imaging (MRI) can potentially be used for non-invasive screening of patients with stable angina pectoris to identify probable obstructive coronary artery disease. MRI-based coronary blood flow quantification has to date only been performed in a 2D fashion, limiting its clinical applicability. In this study, we propose a framework for coronary blood flow quantification using accelerated 4D flow MRI with respiratory motion correction and compressed sensing image reconstruction. We investigate its feasibility and repeatability in healthy subjects at rest. Fourteen healthy subjects received 8 times-accelerated 4D flow MRI covering the left coronary artery (LCA) with an isotropic spatial resolution of 1.0 mm³. Respiratory motion correction was performed based on 1) lung-liver navigator signal, 2) real-time monitoring of foot-head motion of the liver and LCA by a separate acquisition, and 3) rigid image registration to correct for anterior-posterior motion. Time-averaged diastolic LCA flow was determined, as well as time-averaged diastolic maximal velocity (V_MAX) and diastolic peak velocity (V_PEAK). 2D flow MRI scans of the LCA were acquired for reference. Scan-rescan repeatability and agreement between 4D flow MRI and 2D flow MRI were assessed in terms of concordance correlation coefficient (CCC) and coefficient of variation (CV). The protocol resulted in good visibility of the LCA in 11 out of 14 subjects (six female, five male, aged 28 ± 4 years). The other 3 subjects were excluded from analysis. Time-averaged diastolic LCA flow measured by 4D flow MRI was 1.30 ± 0.39 ml/s and demonstrated good scan-rescan repeatability (CCC/CV = 0.79/20.4%). Time-averaged diastolic V_MAX (17.2 ± 3.0 cm/s) and diastolic V_PEAK (24.4 ± 6.5 cm/s) demonstrated moderate repeatability (CCC/CV = 0.52/19.0% and 0.68/23.0%, respectively). 4D flow- and 2D flow-based diastolic LCA flow agreed well (CCC/CV = 0.75/20.1%). Agreement between 4D flow MRI and 2D flow MRI was moderate for both diastolic V_MAX and V_PEAK (CCC/CV = 0.68/20.3% and 0.53/27.0%, respectively). In conclusion, the proposed framework of accelerated 4D flow MRI equipped with respiratory motion correction and compressed sensing image reconstruction enables repeatable diastolic LCA flow quantification that agrees well with 2D flow MRI.

Phase2Phase: Respiratory Motion-Resolved Reconstruction of Free-Breathing Magnetic Resonance Imaging Using Deep Learning Without a Ground Truth for Improved Liver Imaging

OBJECTIVES

Respiratory binning of free-breathing magnetic resonance imaging data reduces motion blurring; however, it exacerbates noise and introduces severe artifacts due to undersampling. Deep neural networks can remove artifacts and noise but usually require high-quality ground truth images for training. This study aimed to develop a network that can be trained without this requirement.

MATERIALS AND METHODS

This retrospective study was conducted on 33 participants enrolled between November 2016 and June 2019. Free-breathing magnetic resonance imaging was performed using a radial acquisition. Self-navigation was used to bin the k-space data into 10 respiratory phases. To simulate short acquisitions, subsets of radial spokes were used in reconstructing images with multicoil nonuniform fast Fourier transform (MCNUFFT), compressed sensing (CS), and 2 deep learning methods: UNet3DPhase and Phase2Phase (P2P). UNet3DPhase was trained using a high-quality ground truth, whereas P2P was trained using noisy images with streaking artifacts. Two radiologists blinded to the reconstruction methods independently reviewed the sharpness, contrast, and artifact-freeness of the end-expiration images reconstructed from data collected at 16% of the Nyquist sampling rate. The generalized estimating equation method was used for statistical comparison. Motion vector fields were derived to examine the respiratory motion range of 4-dimensional images reconstructed using different methods.

RESULTS

A total of 15 healthy participants and 18 patients with hepatic malignancy (50 ± 15 years, 6 women) were enrolled. Both reviewers found that the UNet3DPhase and P2P images had higher contrast (P < 0.01) and fewer artifacts (P < 0.01) than the CS images. The UNet3DPhase and P2P images were reported to be sharper than the CS images by 1 reviewer (P < 0.01) but not by the other reviewer (P = 0.22, P = 0.18). UNet3DPhase and P2P were similar in sharpness and contrast, whereas UNet3DPhase had fewer artifacts (P < 0.01). The motion vector lengths for the MCNUFFT800 and P2P800 images were comparable (10.5 ± 4.2 mm and 9.9 ± 4.0 mm, respectively), whereas both were significantly larger than CS2000 (7.0 ± 3.9 mm; P < 0.0001) and UNnet3DPhase800 (6.9 ± 3.2; P < 0.0001) images.

CONCLUSIONS

Without a ground truth, P2P can reconstruct sharp, artifact-free, and high-contrast respiratory motion-resolved images from highly undersampled data. Unlike the CS and UNet3DPhase methods, P2P did not artificially reduce the respiratory motion range.

Volumetric coronary endothelial function assessment: a feasibility study exploiting stack-of-stars 3D cine MRI and image-based respiratory self-gating

Abnormal coronary endothelial function (CEF), manifesting as depressed vasoreactive responses to endothelial-specific stressors, occurs early in atherosclerosis, independently predicts cardiovascular events, and responds to cardioprotective interventions. CEF is spatially heterogeneous along a coronary artery in patients with atherosclerosis, and thus recently developed and tested non-invasive 2D MRI techniques to measure CEF may not capture the extent of changes in CEF in a given coronary artery. The purpose of this study was to develop and test the first volumetric coronary 3D MRI cine method for assessing CEF along the proximal and mid-coronary arteries with isotropic spatial resolution and in free-breathing. This approach, called 3D-Stars, combines a 6 min continuous, untriggered golden-angle stack-of-stars acquisition with a novel image-based respiratory self-gating method and cardiac and respiratory motion-resolved reconstruction. The proposed respiratory self-gating method agreed well with respiratory bellows and center-of-k-space methods. In healthy subjects, 3D-Stars vessel sharpness was non-significantly different from that by conventional 2D radial in proximal segments, albeit lower in mid-portions. Importantly, 3D-Stars detected normal vasodilatation of the right coronary artery in response to endothelial-dependent isometric handgrip stress in healthy subjects. Coronary artery cross-sectional areas measured using 3D-Stars were similar to those from 2D radial MRI when similar thresholding was used. In conclusion, 3D-Stars offers good image quality and shows feasibility for non-invasively studying vasoreactivity-related lumen area changes along the proximal coronary artery in 3D during free-breathing.

Respiratory Motion-Registered Isotropic Whole-Heart T2 Mapping in Patients With Acute Non-ischemic Myocardial Injury

Background: T₂ mapping is a magnetic resonance imaging technique that can be used to detect myocardial edema and inflammation. However, the focal nature of myocardial inflammation may render conventional 2D approaches suboptimal and make whole-heart isotropic 3D mapping desirable. While self-navigated 3D radial T₂ mapping has been demonstrated to work well at a magnetic field strength of 3T, it results in too noisy maps at 1.5T. We therefore implemented a novel respiratory motion-resolved compressed-sensing reconstruction in order to improve the 3D T₂ mapping precision and accuracy at 1.5T, and tested this in a heterogeneous patient cohort. Materials and Methods: Nine healthy volunteers and 25 consecutive patients with suspected acute non-ischemic myocardial injury (sarcoidosis, n = 19; systemic sclerosis, n = 2; acute graft rejection, n = 2, and myocarditis, n = 2) were included. The free-breathing T₂ maps were acquired as three ECG-triggered T₂-prepared 3D radial volumes. A respiratory motion-resolved reconstruction was followed by image registration of the respiratory states and pixel-wise T₂ mapping. The resulting 3D maps were compared to routine 2D T₂ maps. The T₂ values of segments with and without late gadolinium enhancement (LGE) were compared in patients. Results: In the healthy volunteers, the myocardial T₂ values obtained with the 2D and 3D techniques were similar (45.8 ± 1.8 vs. 46.8 ± 2.9 ms, respectively; P = 0.33). Conversely, in patients, T₂ values did differ between 2D (46.7 ± 3.6 ms) and 3D techniques (50.1 ± 4.2 ms, P = 0.004). Moreover, with the 2D technique, T₂ values of the LGE-positive segments were similar to those of the LGE-negative segments (T_2LGE-= 46.2 ± 3.7 vs. T_2LGE+ = 47.6 ± 4.1 ms; P = 0.49), whereas the 3D technique did show a significant difference (T_2LGE- = 49.3 ± 6.7 vs. T_2LGE+ = 52.6 ± 8.7 ms, P = 0.006). Conclusion: Respiratory motion-registered 3D radial imaging at 1.5T led to accurate isotropic 3D whole-heart T₂ maps, both in the healthy volunteers and in a small patient cohort with suspected non-ischemic myocardial injury. Significantly higher T₂ values were found in patients as compared to controls in 3D but not in 2D, suggestive of the technique's potential to increase the sensitivity of CMR at earlier stages of disease. Further study will be needed to demonstrate its accuracy.

Pilot tone navigation for respiratory and cardiac motion-resolved free-running 5D flow MRI

Purpose

In this work, we integrated the pilot tone (PT) navigation system into a reconstruction framework for respiratory and cardiac motion‐resolved 5D flow. We tested the hypotheses that PT would provide equivalent respiratory curves, cardiac triggers, and corresponding flow measurements to a previously established self‐gating (SG) technique while being independent from changes to the acquisition parameters.

Methods

Fifteen volunteers and 9 patients were scanned with a free‐running 5D flow sequence, with PT integrated. Respiratory curves and cardiac triggers from PT and SG were compared across all subjects. Flow measurements from 5D flow reconstructions using both PT and SG were compared to each other and to a reference electrocardiogram‐gated and respiratory triggered 4D flow acquisition. Radial trajectories with variable readouts per interleave were also tested in 1 subject to compare cardiac trigger quality between PT and SG.

Results

The correlation between PT and SG respiratory curves were 0.95 ± 0.06 for volunteers and 0.95 ± 0.04 for patients. Heartbeat duration measurements in volunteers and patients showed a bias to electrocardiogram measurements of, respectively, 0.16 ± 64.94 ms and 0.01 ± 39.29 ms for PT versus electrocardiogram and of 0.24 ± 63.68 ms and 0.09 ± 32.79 ms for SG versus electrocardiogram. No significant differences were reported for the flow measurements between 5D flow PT and from 5D flow SG. A decrease in the cardiac triggering quality of SG was observed for increasing readouts per interleave, whereas PT quality remained constant.

Conclusion

PT has been successfully integrated in 5D flow MRI and has shown equivalent results to the previously described 5D flow SG technique, while being completely acquisition‐independent.

Coronary Magnetic Resonance Angiography in Chronic Coronary Syndromes

Cardiovascular disease is the leading cause of mortality worldwide, with atherosclerotic coronary artery disease (CAD) accounting for the majority of cases. X-ray coronary angiography and computed tomography coronary angiography (CCTA) are the imaging modalities of choice for the assessment of CAD. However, the use of ionising radiation and iodinated contrast agents remain drawbacks. There is therefore a clinical need for an alternative modality for the early identification and longitudinal monitoring of CAD without these associated drawbacks. Coronary magnetic resonance angiography (CMRA) could be a potential alternative for the detection and monitoring of coronary arterial stenosis, without exposing patients to ionising radiation or iodinated contrast agents. Further advantages include its versatility, excellent soft tissue characterisation and suitability for repeat imaging. Despite the early promise of CMRA, widespread clinical utilisation remains limited due to long and unpredictable scan times, onerous scan planning, lower spatial resolution, as well as motion related image quality degradation. The past decade has brought about a resurgence in CMRA technology, with significant leaps in image acceleration, respiratory and cardiac motion estimation and advanced motion corrected or motion-resolved image reconstruction. With the advent of artificial intelligence, great advances are also seen in deep learning-based motion estimation, undersampled and super-resolution reconstruction promising further improvements of CMRA. This has enabled high spatial resolution (1 mm isotropic), 3D whole heart CMRA in a clinically feasible and reliable acquisition time of under 10 min. Furthermore, latest super-resolution image reconstruction approaches which are currently under evaluation promise acquisitions as short as 1 min. In this review, we will explore the recent technological advances that are designed to bring CMRA closer to clinical reality.