Whole-heart coronary MR angiography using a 3D cones phyllotaxis trajectory

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.

Automated detection of cardiac rest period for trigger delay calculation for image-based navigator coronary magnetic resonance angiography

BACKGROUND

Coronary magnetic resonance angiography (coronary MRA) is increasingly being considered as a clinically viable method to investigate coronary artery disease (CAD). Accurate determination of the trigger delay to place the acquisition window within the quiescent part of the cardiac cycle is critical for coronary MRA in order to reduce cardiac motion. This is currently reliant on operator-led decision making, which can negatively affect consistency of scan acquisition. Recently developed deep learning (DL) derived software may overcome these issues by automation of cardiac rest period detection.

METHODS

Thirty individuals (female, n = 10) were investigated using a 0.9 mm isotropic image-navigator (iNAV)-based motion-corrected coronary MRA sequence. Each individual was scanned three times utilising different strategies for determination of the optimal trigger delay: (1) the DL software, (2) an experienced operator decision, and (3) a previously utilised formula for determining the trigger delay. Methodologies were compared using custom-made analysis software to assess visible coronary vessel length and coronary vessel sharpness for the entire vessel length and the first 4 cm of each vessel.

RESULTS

There was no difference in image quality between any of the methodologies for determination of the optimal trigger delay, as assessed by visible coronary vessel length, coronary vessel sharpness for each entire vessel and vessel sharpness for the first 4 cm of the left mainstem, left anterior descending or right coronary arteries. However, vessel length of the left circumflex was slightly greater using the formula method. The time taken to calculate the trigger delay was significantly lower for the DL-method as compared to the operator-led approach (106 ± 38.0 s vs 168 ± 39.2 s, p < 0.01, 95% CI of difference 25.5-98.1 s).

CONCLUSIONS

Deep learning-derived automated software can effectively and efficiently determine the optimal trigger delay for acquisition of coronary MRA and thus may simplify workflow and improve reproducibility.

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.

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.

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.

Three-dimensional Yarnball k-space acquisition for accelerated MRI

PURPOSE

To introduce an efficient sampling technique named Yarnball, which may serve as a direct alternative to 3D Cones.

METHODS

Yarnball evolves through 3D k-space with increasing loop size, and the differential equations defining this flexible trajectory are presented in detail. The sampling efficiencies of Yarnball and 3D Cones were compared through point spread function analysis and simulated imaging (which highlights undersampling in the absence of other scanning effects). The feasibility of Yarnball implementation was demonstrated for fully sampled T₁ -weighted images of the human head at 3 T.

RESULTS

The mostly large 3D loops of the Yarnball trajectory facilitate rapid sampling under peripheral nerve stimulation constraint, an advantage that increases with readout duration (T_RO ). Point spread function analysis yielded 89% (T_RO = 2 ms) and 77% (T_RO = 10 ms) of Yarnball voxels with magnitude less than 0.01% of the point spread function peak. For 3D Cones, these values were only 52% and 29%. The 3D-Cones technique required 1.4 times (T_RO = 2 ms) and 1.8 times (T_RO = 10 ms) more trajectories than Yarnball to produce simulated images of a sphere free from undersampling artifact. For a prolate spheroidal (head-like) object, 1.75 times and 2.6 times more trajectories were required for 3D Cones. Yarnball produced 0.72 mm (1/2k_max ) isotropic T₁ -weighted human brain images free from undersampling artifact in only 98 seconds at 3 T.

CONCLUSION

Yarnball demonstrated greater k-space sampling efficiency than directly comparable 3D Cones, and may have value wherever 3D Cones has been considered. Yarnball may also have value in the context of rapid T₁ -weighted brain imaging.

Combined T2 -preparation and multidimensional outer volume suppression for coronary artery imaging with 3D cones trajectories

PURPOSE

To develop a modular magnetization preparation sequence for combined T2 -preparation and multidimensional outer volume suppression (OVS) for coronary artery imaging.

METHODS

A combined T2 -prepared 1D OVS sequence with fat saturation was defined to contain a 90°-60 180°60 composite nonselective tip-down pulse, two 180°Y hard pulses for refocusing, and a -90° spectral-spatial sinc tip-up pulse. For 2D OVS, 2 modules were concatenated, selective in X and then Y. Bloch simulations predicted robustness of the sequence to B0 and B1 inhomogeneities. The proposed sequence was compared with a T2 -prepared 2D OVS sequence proposed by Luo et al, which uses a spatially selective 2D spiral tip-up. The 2 sequences were compared in phantom studies and in vivo coronary artery imaging studies with a 3D cones trajectory.

RESULTS

Phantom results demonstrated superior OVS for the proposed sequence compared with the Luo sequence. In studies on 15 healthy volunteers, the proposed sequence had superior image edge profile acutance values compared with the Luo sequence for the right (P < .05) and left (P < .05) coronary arteries, suggesting superior vessel sharpness. The proposed sequence also had superior signal-to-noise ratio (P < .05) and passband-to-stopband ratio (P < .05). Reader scores and reader preference indicated superior coronary image quality of the proposed sequence for both the right (P < .05) and left (P < .05) coronary arteries.

CONCLUSION

The proposed sequence with concatenated 1D spatially selective tip-ups and integrated fat saturation has superior image quality and suppression compared with the Luo sequence with 2D spatially selective tip-up.

Reconstruction of undersampled 3D non-Cartesian image-based navigators for coronary MRA using an unrolled deep learning model

PURPOSE

To rapidly reconstruct undersampled 3D non-Cartesian image-based navigators (iNAVs) using an unrolled deep learning (DL) model, enabling nonrigid motion correction in coronary magnetic resonance angiography (CMRA).

METHODS

An end-to-end unrolled network is trained to reconstruct beat-to-beat 3D iNAVs acquired during a CMRA sequence. The unrolled model incorporates a nonuniform FFT operator in TensorFlow to perform the data-consistency operation, and the regularization term is learned by a convolutional neural network (CNN) based on the proximal gradient descent algorithm. The training set includes 6,000 3D iNAVs acquired from 7 different subjects and 11 scans using a variable-density (VD) cones trajectory. For testing, 3D iNAVs from 4 additional subjects are reconstructed using the unrolled model. To validate reconstruction accuracy, global and localized motion estimates from DL model-based 3D iNAVs are compared with those extracted from 3D iNAVs reconstructed with l 1 -ESPIRiT. Then, the high-resolution coronary MRA images motion corrected with autofocusing using the l 1 -ESPIRiT and DL model-based 3D iNAVs are assessed for differences.

RESULTS

3D iNAVs reconstructed using the DL model-based approach and conventional l 1 -ESPIRiT generate similar global and localized motion estimates and provide equivalent coronary image quality. Reconstruction with the unrolled network completes in a fraction of the time compared to CPU and GPU implementations of l 1 -ESPIRiT (20× and 3× speed increases, respectively).

CONCLUSIONS

We have developed a deep neural network architecture to reconstruct undersampled 3D non-Cartesian VD cones iNAVs. Our approach decreases reconstruction time for 3D iNAVs, while preserving the accuracy of nonrigid motion information offered by them for correction.