Retrospective reconstruction of cardiac cine images from golden-ratio radial MRI using one-dimensional navigators

High-resolution CINE imaging of active guided knee motion using continuously acquired golden-angle radial MRI and rotary sensor information

To explore and extend on dynamic imaging of joint motion, an MRI-safe device guiding knee motion with an attached rotary encoder was used in MRI measurements of multiple knee flexion-extension cycles using radial gradient echo imaging with the golden-angle as azimuthal angle increment. Reproducibility of knee motion was investigated. Real-time and CINE mode anatomical images were reconstructed for different knee flexion angles by synchronizing the encoder information with the MRI data, and performing flexion angle selective gating across multiple motion cycles. When investigating the influence of the rotation angle window width on reconstructed CINE images, it was found that angle windows between 0.5° and 3° exhibited acceptable image sharpness without suffering from significant motion-induced blurring. Furthermore, due to flexible retrospective image reconstruction afforded by the radial golden-angle imaging, the number of motion cycles included in the reconstruction could be retrospectively reduced to investigate the corresponding influence of acquisition time on image quality. Finally, motion reproducibility between motion cycles and accuracy of the flexion angle selective gating were sufficient to acquire whole-knee 3D dynamic imaging with a retrospectively gated 3D cone UTE sequence.

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.

A novel multipurpose device for guided knee motion and loading during dynamic magnetic resonance imaging

INTRODUCTION

This work aimed to develop a novel multipurpose device for guided knee flexion-extension, both passively using a motorized pneumatic system and actively (muscle-driven) with the joint unloaded or loaded during dynamic MRI. Secondary objectives were to characterize the participant experience during device use, and present preliminary dynamic MRI data to demonstrate the different device capabilities.

MATERIAL AND METHODS

Self-reported outcomes were used to characterize the pain, physical exertion and discomfort levels experienced by 10 healthy male participants during four different active knee motion and loading protocols using the novel device. Knee angular data were recorded during the protocols to determine the maximum knee range of motion achievable. Dynamic MRI was acquired for three healthy volunteers during passive, unloaded knee motion using 2D Cartesian TSE, 2D radial GRE and 3D UTE sequences; and during active, unloaded and loaded knee motion using 2D radial GRE imaging. Because of the different MRI sequences used, spatial resolution was inherently lower for active knee motion than for passive motion acquisitions.

RESULTS

Depending on the protocol, some participants reported slight pain, mild discomfort and varying levels of physical exertion. On average, participants achieved ∼40° of knee flexion; loaded conditions can create knee moments up to 27Nm. High quality imaging data were obtained during different motion and loading conditions. Dynamic 3D data allowed to retrospectively extract arbitrarily oriented slices.

CONCLUSION

A novel multipurpose device for guided, physiologically relevant knee motion and loading during dynamic MRI was developed. Device use was well tolerated and suitable for acquiring high quality images during different motion and loading conditions. Different bone positions between loaded and unloaded conditions were likely due to out-of-plane motion, particularly because image registration was not performed. Ultimately, this device could be used to advance our understanding of physiological and pathological joint mechanics.

Phase contrast coronary blood velocity mapping with both high temporal and spatial resolution using triggered Golden Angle rotated Spiral k-t Sparse Parallel imaging (GASSP) with shifted binning

PURPOSE

High temporal and spatial resolutions are required for coronary blood flow measures. Current spiral breath-hold phase contrast (PC) MRI at 3T focus on either high spatial or high temporal resolution. We propose a golden angle (GA) rotated Spiral k-t Sparse Parallel imaging (GASSP) sequence for both high spatial (0.8 mm) and high temporal (<21 ms) resolutions.

METHODS

GASSP PC data are acquired in left anterior descending and right coronary arteries of eight healthy subjects. Binning of GA rotated spiral data into cardiac frames may lead to large k-space gaps. To reduce those gaps, the binning window is shifted and a triggered GA scheme that resets the rotation angle every heartbeat is proposed. The gap reductions are evaluated in simulations and all subjects. Peak systolic velocity (PSV), peak diastolic velocity (PDV), coronary blood flow rate, and vessel area are validated against two reference scans, and repeatability/reproducibility are determined.

RESULTS

Shifted binning reduced the mean k-space gaps of the triggered GA scheme by 14°-22° in simulations and about 20° in vivo. The k-space gap across three cardiac frames was reduced with the triggered GA scheme compared to the standard GA scheme (35.3°± 3.6° vs. 43°± 13.7°, t-test P = .04). PSV, PDV, flow rate, and area had high intra-scan repeatability (0.92 ≤ intraclass correlation coefficient [ICC] ≤ 0.99), and inter-scan (0.78 ≤ ICC ≤ 0.91) and intra-observer (0.91 ≤ ICC ≤ 0.98) reproducibility.

CONCLUSION

GASSP enables single breath-hold coronary PC MRI with high temporal and spatial resolutions. Shifted binning and a triggered GA scheme reduce k-space gaps. Quantitative coronary flow metrics are highly reproducible, especially within the same scanning session.

Elastic AlignedSENSE for Dynamic MR Reconstruction: A Proof of Concept in Cardiac Cine

Numerous methods in the extensive literature on magnetic resonance imaging (MRI) reconstruction exploit temporal redundancy to accelerate cardiac cine. Some of them include motion compensation, which involves high computational costs and long runtimes. In this work, we proposed a method—elastic alignedSENSE (EAS)—for the direct reconstruction of a motion-free image plus a set of nonrigid deformations to reconstruct a 2D cardiac sequence. The feasibility of the proposed approach was tested in 2D Cartesian and golden radial multi-coil breath-hold cardiac cine acquisitions. The proposed approach was compared against parallel imaging compressed sense (sPICS) and group-wise motion corrected compressed sense (GWCS) reconstructions. EAS provides better results on objective measures with considerable less runtime when an acceleration factor is higher than 10×. Subjective assessment of an expert, however, invited proposing the combination of EAS and GWCS as a preferable alternative to GWCS or EAS in isolation.

Self-calibrated interpolation of non-Cartesian data with GRAPPA in parallel imaging

PURPOSE

To develop a non-Cartesian k-space reconstruction method using self-calibrated region-specific interpolation kernels for highly accelerated acquisitions.

METHODS

In conventional non-Cartesian GRAPPA with through-time GRAPPA (TT-GRAPPA), the use of region-specific interpolation kernels has demonstrated improved reconstruction quality in dynamic imaging for highly accelerated acquisitions. However, TT-GRAPPA requires the acquisition of a large number of separate calibration scans. To reduce the overall imaging time, we propose Self-calibrated Interpolation of Non-Cartesian data with GRAPPA (SING) to self-calibrate region-specific interpolation kernels from dynamic undersampled measurements. The SING method synthesizes calibration data to adapt to the distinct shape of each region-specific interpolation kernel geometry, and uses a novel local k-space regularization through an extension of TT-GRAPPA. This calibration approach is used to reconstruct non-Cartesian images at high acceleration rates while mitigating noise amplification. The reconstruction quality of SING is compared with conjugate-gradient SENSE and TT-GRAPPA in numerical phantoms and in vivo cine data sets.

RESULTS

In both numerical phantom and in vivo cine data sets, SING offers visually and quantitatively similar reconstruction quality to TT-GRAPPA, and provides improved reconstruction quality over conjugate-gradient SENSE. Furthermore, temporal fidelity in SING and TT-GRAPPA is similar for the same acceleration rates. G-factor evaluation over the heart shows that SING and TT-GRAPPA provide similar noise amplification at moderate and high rates.

CONCLUSION

The proposed SING reconstruction enables significant improvement of acquisition efficiency for calibration data, while matching the reconstruction performance of TT-GRAPPA.

Reconstruction techniques for cardiac cine MRI

The present survey describes the state-of-the-art techniques for dynamic cardiac magnetic resonance image reconstruction. Additionally, clinical relevance, main challenges, and future trends of this image modality are outlined. Thus, this paper aims to provide a general vision about cine MRI as the standard procedure in functional evaluation of the heart, focusing on technical methodologies.

Noncontrast MR angiography: An update

Both computed tomography (CT) angiography (CTA) and contrast-enhanced MR angiography (CEMRA) have proven to be useful and accurate cross-sectional imaging modalities over a wide range of vascular territories and vascular disorders. A key advantage of MRA is that, unlike CTA, it can be performed without the administration of a contrast agent. In this review article we consider the motivations for using noncontrast MRA, potential contrast mechanisms, imaging techniques, advantages, and drawbacks with respect to CTA and CEMRA, and the level of evidence for using the various MRA techniques. In addition, we explore new developments that promise to expand the reliability and range of clinical applications for noncontrast MRA, along with functional MRA capabilities not available with CTA or CEMRA. Level of Evidence: 1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:355-373.

Improved respiratory self-navigation for 3D radial acquisitions through the use of a pencil-beam 2D-T2 -prep for free-breathing, whole-heart coronary MRA

PURPOSE

In respiratory self-navigation (SN), signal from static structures, such as the chest wall, may complicate motion detection or introduce post-correction artefacts. Suppressing signal from superfluous tissues may therefore improve image quality. We thus test the hypothesis that SN whole-heart coronary magnetic resonance angiography (MRA) will benefit from an outer-volume suppressing 2D-T2 -Prep and present both phantom and in vivo results.

METHODS

A 2D-T2 -Prep and a conventional T2 -Prep were used prior to a free-breathing 3D-radial SN sequence. Both techniques were compared by imaging a home-built moving cardiac phantom and by performing coronary MRA in nine healthy volunteers. Reconstructions were performed using both a reference-based and a reference-independent approach to motion tracking, along with several coil combinations. Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were compared, along with vessel sharpness (VS).

RESULTS

In phantoms, using the 2D-T2 -Prep increased SNR by 16% to 53% and mean VS by 8%; improved motion tracking precision was also achieved. In volunteers, SNR increased by an average of 29% to 33% in the blood pool and by 15% to 25% in the myocardium, depending on the choice of reconstruction coils and algorithm, and VS increased by 34%.

CONCLUSION

A 2D-T2 -Prep significantly improves image quality in both phantoms and volunteers when performing SN coronary MRA. Magn Reson Med 79:1293-1303, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Cardiac 4D phase-contrast CMR at 9.4 T using self-gated ultra-short echo time (UTE) imaging

BACKGROUND

Time resolved 4D phase contrast (PC) cardiovascular magnetic resonance (CMR) in mice is challenging due to long scan times, small animal ECG-gating and the rapid blood flow and cardiac motion of small rodents. To overcome several of these technical challenges we implemented a retrospectively self-gated 4D PC radial ultra-short echo-time (UTE) acquisition scheme and assessed its performance in healthy mice by comparing the results with those obtained with an ECG-triggered 4D PC fast low angle shot (FLASH) sequence.

METHODS

Cardiac 4D PC CMR images were acquired at 9.4 T in healthy mice using the proposed self-gated radial center-out UTE acquisition scheme (TE/TR of 0.5 ms/3.1 ms) and a standard Cartesian 4D PC imaging sequence (TE/TR of 2.1 ms/5.0 ms) with a four-point Hadamard flow encoding scheme. To validate the proposed UTE flow imaging technique, experiments on a flow phantom with variable pump rates were performed.

RESULTS

The anatomical images and flow velocity maps of the proposed 4D PC UTE technique showed reduced artifacts and an improved SNR (left ventricular cavity (LV): 8.9 ± 2.5, myocardium (MC): 15.7 ± 1.9) compared to those obtained using a typical Cartesian FLASH sequence (LV: 5.6 ± 1.2, MC: 10.1 ± 1.4) that was used as a reference. With both sequences comparable flow velocities were obtained in the flow phantom as well as in the ascending aorta (UTE: 132.8 ± 18.3 cm/s, FLASH: 134.7 ± 13.4 cm/s) and pulmonary artery (UTE: 78.5 ± 15.4 cm/s, FLASH: 86.6 ± 6.2 cm/s) of the animals. Self-gated navigator signals derived from information of the oversampled k-space center were successfully extracted for all animals with a higher gating efficiency of time spent on acquiring gated data versus total measurement time (UTE: 61.8 ± 11.5%, FLASH: 48.5 ± 4.9%).

CONCLUSIONS

The proposed self-gated 4D PC UTE sequence enables robust and accurate flow velocity mapping of the mouse heart in vivo at high magnetic fields. At the same time SNR, gating efficiency, flow artifacts and image quality all improved compared to the images obtained using the well-established, ECG-triggered, 4D PC FLASH sequence.

Retrospective 3D motion correction using spherical navigator echoes

PURPOSE

To develop and evaluate a rapid spherical navigator echo (SNAV) motion correction technique, then apply it for retrospective correction of brain images.

METHODS

The pre-rotated, template matching SNAV method (preRot-SNAV) was developed in combination with a novel hybrid baseline strategy, which includes acquired and interpolated templates. Specifically, the SNAV templates are only rotated around X- and Y-axis; for each rotated SNAV, simulated baseline templates that mimic object rotation about the Z-axis were interpolated. The new method was first evaluated with phantom experiments. Then, a customized SNAV-interleaved gradient echo sequence was used to image three volunteers performing directed head motion. The SNAV motion measurements were used to retrospectively correct the brain images. Experiments were performed using a 3.0T whole-body MRI scanner and both single and 8-channel head coils.

RESULTS

Phantom rotations and translations measured using the hybrid baselines agreed to within 0.9° and 1mm compared to those measured with the original preRot-SNAV method. Retrospective motion correction of in vivo images using the hybrid preRot-SNAV effectively corrected for head rotation up to 4° and 4mm.

CONCLUSIONS

The presented hybrid approach enables the acquisition of pre-rotated baseline templates in as little as 2.5s, and results in accurate measurement of rotations and translations. Retrospective 3D motion correction successfully reduced motion artifacts in vivo.

USPIO-enhanced 3D-cine self-gated cardiac MRI based on a stack-of-stars golden angle short echo time sequence: Application on mice with acute myocardial infarction

PURPOSE

To develop and assess a 3D-cine self-gated method for cardiac imaging of murine models.

MATERIALS AND METHODS

A 3D stack-of-stars (SOS) short echo time (STE) sequence with a navigator echo was performed at 7T on healthy mice (n = 4) and mice with acute myocardial infarction (MI) (n = 4) injected with ultrasmall superparamagnetic iron oxide (USPIO) nanoparticles. In all, 402 spokes were acquired per stack with the incremental or the golden angle method using an angle increment of (360/402)° or 222.48°, respectively. A cylindrical k-space was filled and repeated with a maximum number of repetitions (NR) of 10. 3D cine cardiac images at 156 μm resolution were reconstructed retrospectively and compared for the two methods in terms of contrast-to-noise ratio (CNR). The golden angle images were also reconstructed with NR = 10, 6, and 3, to assess cardiac functional parameters (ejection fraction, EF) on both animal models.

RESULTS

The combination of 3D SOS-STE and USPIO injection allowed us to optimize the identification of cardiac peaks on navigator signal and generate high CNR between blood and myocardium (15.3 ± 1.0). The golden angle method resulted in a more homogeneous distribution of the spokes inside a stack (P < 0.05), enabling reducing the acquisition time to 15 minutes. EF was significantly different between healthy and MI mice (P < 0.05).

CONCLUSION

The method proposed here showed that 3D-cine images could be obtained without electrocardiogram or respiratory gating in mice. It allows precise measurement of cardiac functional parameters even on MI mice. J. Magn. Reson. Imaging 2016;44:355-365.

Segmented golden ratio radial reordering with variable temporal resolution for dynamic cardiac MRI

PURPOSE

Golden ratio (GR) radial reordering allows for retrospective choice of temporal resolution by providing a near-uniform k-space sampling within any reconstruction window. However, when applying GR to electrocardiogram (ECG)-gated cardiac imaging, the k-space coverage may not be as uniform because a single reconstruction window is broken into several temporally isolated ones. The goal of this study was to investigate the image artifacts caused by applying GR to ECG-gated cardiac imaging and to propose a segmented GR method to address this issue.

METHODS

Computer simulation and phantom experiments were used to evaluate the image artifacts resulting from three k-space sampling patterns (ie, uniform radial, conventional GR, and segmented GR). Two- and three-dimensional cardiac cine images were acquired in seven healthy subjects. Imaging artifacts due to k-space sampling nonuniformity were graded on a 5-point scale by an experienced cardiac imaging reader.

RESULTS

Segmented GR provides more uniform k-space sampling that is independent of heart-rate variation than conventional GR. Cardiac cine images using segmented GR have significantly higher and more reliable image quality than conventional GR.

CONCLUSION

Segmented GR successfully addresses the nonuniform sampling that occurs with combining conventional GR with ECG gating. This technique can potentially be applied to any ECG-gated cardiac imaging application to allow for retrospective selection of a reconstruction window. Magn Reson Med 76:94-103, 2016. © 2015 Wiley Periodicals, Inc.

Self-gated cardiac Cine MRI of the rat on a clinical 3 T MRI system

The ability to perform small animal functional cardiac imaging on clinical MRI scanners may be of particular value in cases in which the availability of a dedicated high field animal MRI scanner is limited. Here, we propose radial MR cardiac imaging in the rat on a whole-body clinical 3 T scanner in combination with interspersed projection navigators for self-gating without any additional external triggering requirements for electrocardiogram (ECG) and respiration. Single navigator readouts were interspersed using the same TR and a high navigator frequency of 54 Hz into a radial golden-angle acquisition. The extracted navigator function was thresholded to exclude data for reconstruction from inhalation phases during the breathing cycle, enabling free breathing acquisition. To minimize flow artifacts in the dynamic cine images a center-out half echo radial acquisition scheme with ramp sampling was used. Navigator functions were derived from the corresponding projection navigator data from which both respiration and cardiac cycles were extracted. Self-gated cine acquisition resulted in high-quality cardiac images which were free of major artifacts with spatial resolution of up to 0.21 × 0.21 × 1.00 mm(3) and a contrast-to-noise ratio (CNR) of 21 ± 3 between the myocardium and left ventricle. Self-gated golden ratio based radial acquisition successfully acquired cine images of the rat heart on a clinical MRI system without the need for dedicated animal ECG equipment.

Intrinsic correction of system delays for radial magnetic resonance imaging

INTRODUCTION

When using radial MR image acquisition techniques gradient or sampling delays due to hardware imperfections can cause mismatch between the expected and the actual k-space trajectory along the readout direction. To provide a robust and simple correction of such system delays we developed a new calibration method which is independent of using any reference data or applying sequence modifications.

MATERIAL AND METHODS

Radial data obtained with 180°, 360° and golden-angle radial ordering schemes were deliberately shifted along the readout direction for a discrete range of gradient delays. Following 2D regridding, images were reconstructed and analyzed in image space for all applied shifts to estimate the optimal system delay. Phantom and in vivo measurements were performed to test the robustness of the algorithm.

RESULTS

Using the 360° and golden-angle radial ordering schemes system delays in the range of 3.3μs to 6.3μs were estimated and corrected for several imaging applications and different conditions, including cardiac and real-time MRI as well as multiple acquisitions using different imaging parameters and slice orientations. When using the standard 180° radial acquisition scheme no automated correction was possible. With a mean computation time of 23.2±14.0s for the delay estimation computational demands were moderate allowing implementation of the algorithm on the image reconstruction system of any modern MR system.

CONCLUSION

We have demonstrated that radial data acquired with a 360° or golden-angle ordering scheme can be used for reliable intrinsic correction of system delays. The proposed technique enables a per-scan correction of system delays without the need for additional calibration data or modifications of the radial imaging sequence.

Self-gated fat-suppressed cardiac cine MRI

PURPOSE

To develop a self-gated alternating repetition time balanced steady-state free precession (ATR-SSFP) pulse sequence for fat-suppressed cardiac cine imaging.

METHODS

Cardiac gating is computed retrospectively using acquired magnetic resonance self-gating data, enabling cine imaging without the need for electrocardiogram (ECG) gating. Modification of the slice-select rephasing gradients of an ATR-SSFP sequence enables the acquisition of a one-dimensional self-gating readout during the unused short repetition time (TR). Self-gating readouts are acquired during every TR of segmented, breath-held cardiac scans. A template-matching algorithm is designed to compute cardiac trigger points from the self-gating signals, and these trigger points are used for retrospective cine reconstruction. The proposed approach is compared with ECG-gated ATR-SSFP and balanced steady-state free precession in 10 volunteers and five patients.

RESULTS

The difference of ECG and self-gating trigger times has a variability of 13 ± 11 ms (mean ± SD). Qualitative reviewer scoring and ranking indicate no statistically significant differences (P > 0.05) between self-gated and ECG-gated ATR-SSFP images. Quantitative blood-myocardial border sharpness is not significantly different among self-gated ATR-SSFP ( 0.61±0.15 mm -1), ECG-gated ATR-SSFP ( 0.61±0.15 mm -1), or conventional ECG-gated balanced steady-state free precession cine MRI ( 0.59±0.15 mm -1).

CONCLUSION

The proposed self-gated ATR-SSFP sequence enables fat-suppressed cardiac cine imaging at 1.5 T without the need for ECG gating and without decreasing the imaging efficiency of ATR-SSFP.