Cardiac real-time MRI using a pre-emphasized spiral acquisition based on the gradient system transfer function

Joint image reconstruction and segmentation of real-time cardiovascular magnetic resonance imaging in free-breathing using a model based on disentangled representation learning

BACKGROUND

To investigate image quality and agreement of derived cardiac function parameters in a novel joint image reconstruction and segmentation approach based on disentangled representation learning, enabling real-time cardiac cine imaging during free-breathing.

METHODS

A multi-tasking neural network architecture, incorporating disentangled representation learning, was trained using simulated examinations based on data from a public repository along with cardiovascular magnetic resonance (CMR) scans specifically acquired for model development. An exploratory feasibility study evaluated the method on undersampled real-time acquisitions using an in-house developed spiral balanced steady-state free precession pulse sequence in eight healthy participants and five patients with intermittent atrial fibrillation. Images and predicted left ventricle segmentations were compared to the reference standard of electrocardiography (ECG)-gated segmented Cartesian cine with repeated breath-holds and corresponding manual segmentation.

RESULTS

On a 5-point Likert scale, image quality of the real-time breath-hold approach and Cartesian cine was comparable in healthy participants (RT-BH: 1.99 ± 0.98, Cartesian: 1.94 ± 0.86, p = 0.052), but slightly inferior in free-breathing (RT-FB: 2.40 ± 0.98, p < 0.001). In patients with arrhythmia, both real-time approaches demonstrated favorable image quality (RT-BH: 2.10 ± 1.28, p < 0.001, RT-FB: 2.40 ± 1.13, p < 0.01, Cartesian: 2.68 ± 1.13). Intra-observer reliability was good (intraclass correlation coefficient = 0.77, 95% confidence interval [0.75, 0.79], p < 0.001). In functional analysis, a positive bias was observed for ejection fractions derived from the proposed model compared to the clinical reference standard (RT-BH mean: 58.5 ± 5.6%, bias: +3.47%, 95% confidence interval [-0.86, 7.79%], RT-FB mean: 57.9 ± 10.6%, bias: +1.45%, [-3.02, 5.91%], Cartesian mean: 54.9 ± 6.7%).

CONCLUSION

The introduced real-time CMR imaging technique enables high-quality cardiac cine data acquisitions in 1-2 min, eliminating the need for ECG gating and breath-holds. This approach offers a promising alternative to the current clinical practice of segmented acquisition, with shorter scan times, improved patient comfort, and increased robustness to arrhythmia and patient non-compliance.

Feasibility of undersampled spiral trajectories in MREPT for fast conductivity imaging

PURPOSE

To investigate spiral-based imaging including trajectories with undersampling as a fast and robust alternative for phase-based magnetic resonance electrical properties tomography (MREPT) techniques.

METHODS

Spiral trajectories with various undersampling ratios were prescribed to acquire images from an experimental phantom and a healthy volunteer at 3T. The non-Cartesian acquisitions were reconstructed using SPIRiT, and conductivity maps were derived using phase-based cr-MREPT. The resulting maps were compared between different sampling trajectories. Additionally, a conductivity map was obtained using a Cartesian balanced SSFP acquisition from the volunteer to comparatively demonstrate the robustness of the proposed method.

RESULTS

The phantom and volunteer results illustrate the benefits of the spiral acquisitions. Specifically, undersampled spiral acquisitions display improved robustness against field inhomogeneity artifacts and lowered SD values with shortened readout times. Furthermore, average of conductivity values measured for the cerebrospinal fluid with the spiral acquisitions were 1.703 S/m, indicating a close agreement with the theoretical values of 1.794 S/m.

CONCLUSION

A spiral-based acquisition framework for conductivity imaging with and without undersampling is presented. Overall, spiral-based acquisitions improved robustness against field inhomogeneity artifacts, while achieving whole head coverage with multiple averages in less than a minute.

ECG-free cine MRI with data-driven clustering of cardiac motion for quantification of ventricular function

BACKGROUND

Despite the widespread use of cine MRI for evaluation of cardiac function, existing real-time methods do not easily enable quantification of ventricular function. Moreover, segmented cine MRI assumes periodicity of cardiac motion. We aim to develop a self-gated, cine MRI acquisition scheme with data-driven cluster-based binning of cardiac motion.

METHODS

A Cartesian golden-step balanced steady-state free precession sequence with sorted k-space ordering was designed. Image data were acquired with breath-holding. Principal component analysis and k-means clustering were used for binning of cardiac phases. Cluster compactness in the time dimension was assessed using temporal variability, and dispersion in the spatial dimension was assessed using the Caliński-Harabasz index. The proposed and the reference electrocardiogram (ECG)-gated cine methods were compared using a four-point image quality score, SNR and CNR values, and Bland-Altman analyses of ventricular function.

RESULTS

A total of 10 subjects with sinus rhythm and 8 subjects with arrhythmias underwent cardiac MRI at 3.0 T. The temporal variability was 45.6 ms (cluster) versus 24.6 ms (ECG-based) (p < 0.001), and the Caliński-Harabasz index was 59.1 ± 9.1 (cluster) versus 22.0 ± 7.1 (ECG based) (p < 0.001). In subjects with sinus rhythm, 100% of the end-systolic and end-diastolic images from both the cluster and reference approach received the highest image quality score of 4. Relative to the reference cine images, the cluster-based multiphase (cine) image quality consistently received a one-point lower score (p < 0.05), whereas the SNR and CNR values were not significantly different (p = 0.20). In cases with arrhythmias, 97.9% of the end-systolic and end-diastolic images from the cluster approach received an image quality score of 3 or more. The mean bias values for biventricular ejection fraction and volumes derived from the cluster approach versus reference cine were negligible.

CONCLUSION

ECG-free cine cardiac MRI with data-driven clustering for binning of cardiac motion is feasible and enables quantification of cardiac function.

Beat phenomena of oscillating readouts

PURPOSE

To demonstrate slowly varying, erroneous magnetic field gradients for oscillating readouts due to the mechanically resonant behavior of gradient systems.

METHODS

Projections of a static phantom were acquired using a one-dimensional (1D) EPI sequence with varying EPI frequencies ranging from 1121 to 1580 Hz on clinical 3T systems (30 mT/m, 200 T/m/s). Phase due to static B0 inhomogeneities was eliminated by a complex division of two separate scans with different polarities of the EPI readout. The temporal evolution of phase was evaluated and related to the mechanical resonances of the gradient systems derived from the gradient modulation transfer function. Additionally, the impact of temporally varying mechanical resonance effects on EPI was evaluated using an echo-planar spectroscopic imaging sequence.

RESULTS

A beat phenomenon resulting in a slowly varying phase was observed. Its temporal frequency was given by the difference between the EPI frequency and the mechanical resonance frequency of the activated gradient axis. The maximum erroneous, oscillating phase during phase encoding was ±0.5 rad for an EPI frequency of 1281 Hz. Echo-planar spectroscopic imaging images showed the resulting time-dependent stretching/compression of the FOV.

CONCLUSION

Oscillating readouts such as those used in EPI can result in low-frequency, erroneous phase contributions, which are explained by the beat phenomenon. Therefore, EPI phase-correction approaches may need to include beat effects for accurate image reconstruction.

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.

Do Disinfectant Solutions during Gargling Reach the Pharynx?

INTRODUCTION

In times of COVID-19, gargling disinfectant is commonly used. Disinfectant solutions seem to decrease the infection's symptoms. For disinfection, several techniques are reported. So far, there are no data about the regions in the upper airways achieved by gargled fluid.

METHODS

Ten healthy volunteers without any dysphagia were investigated with a high-sensitivity flexible endoscopic evaluation of swallowing (hsFEES®) during and after gargling colored water. One volunteer repeated the gargling process in fast and real-time MRI.

RESULTS

In all cases, no color accumulation was detected on the posterior pharyngeal wall, epi- or hypopharynx during gargling. The MRI scans confirmed the results.

CONCLUSIONS

hsFEES® and fast MRI provide an insight into the gargling pattern. Data show that during gargling, the fluid covers the soft tissue in the oral cavity and the anterior part of the soft palate, but not the posterior pharyngeal wall nor the epi- and hypopharynx.

Fast measurement of the gradient system transfer function at 7 T

PURPOSE

In this work, a new method to determine the gradient system transfer function (GSTF) with high frequency resolution and high SNR is presented, using fast and simple phantom measurements. The GSTF is an effective instrument for hardware characterization and calibration, which can be used to correct for gradient distortions, or enhance gradient fidelity.

METHODS

The thin-slice approach for phantom-based measurements of the GSTF is expanded by adding excitations that are shifted after the application of the probing gradient, to capture long-lasting field fluctuations with high SNR. A physics-informed regularization procedure is implemented to derive high-quality transfer functions from a small number of measurements. The resulting GSTFs are evaluated by means of gradient time-course estimation and pre-emphasis of a trapezoidal test gradient on a 7T scanner.

RESULTS

The GSTFs determined with the proposed method capture sharp mechanical resonances with a high level of detail. The measured trapezoidal gradient progressions are authentically reproduced by the GSTF estimations on all three axes. The GSTF-based pre-emphasis considerably improves the gradient fidelity in the plateau phase of the test gradient and almost completely eliminates lingering field oscillations.

CONCLUSION

The presented approach allows fast and simple characterization of gradient field fluctuations caused by long-living eddy current and vibration effects, which become more pronounced at ultrahigh field strengths.

Dynamic pressure-volume loop analysis by simultaneous real-time cardiovascular magnetic resonance and left heart catheterization

BACKGROUND

Left ventricular (LV) contractility and compliance are derived from pressure-volume (PV) loops during dynamic preload reduction, but reliable simultaneous measurements of pressure and volume are challenging with current technologies. We have developed a method to quantify contractility and compliance from PV loops during a dynamic preload reduction using simultaneous measurements of volume from real-time cardiovascular magnetic resonance (CMR) and invasive LV pressures with CMR-specific signal conditioning.

METHODS

Dynamic PV loops were derived in 16 swine (n = 7 naïve, n = 6 with aortic banding to increase afterload, n = 3 with ischemic cardiomyopathy) while occluding the inferior vena cava (IVC). Occlusion was performed simultaneously with the acquisition of dynamic LV volume from long-axis real-time CMR at 0.55 T, and recordings of invasive LV and aortic pressures, electrocardiogram, and CMR gradient waveforms. PV loops were derived by synchronizing pressure and volume measurements. Linear regression of end-systolic- and end-diastolic- pressure-volume relationships enabled calculation of contractility. PV loops measurements in the CMR environment were compared to conductance PV loop catheter measurements in 5 animals. Long-axis 2D LV volumes were validated with short-axis-stack images.

RESULTS

Simultaneous PV acquisition during IVC-occlusion was feasible. The cardiomyopathy model measured lower contractility (0.2 ± 0.1 mmHg/ml vs 0.6 ± 0.2 mmHg/ml) and increased compliance (12.0 ± 2.1 ml/mmHg vs 4.9 ± 1.1 ml/mmHg) compared to naïve animals. The pressure gradient across the aortic band was not clinically significant (10 ± 6 mmHg). Correspondingly, no differences were found between the naïve and banded pigs. Long-axis and short-axis LV volumes agreed well (difference 8.2 ± 14.5 ml at end-diastole, -2.8 ± 6.5 ml at end-systole). Agreement in contractility and compliance derived from conductance PV loop catheters and in the CMR environment was modest (intraclass correlation coefficient 0.56 and 0.44, respectively).

CONCLUSIONS

Dynamic PV loops during a real-time CMR-guided preload reduction can be used to derive quantitative metrics of contractility and compliance, and provided more reliable volumetric measurements than conductance PV loop catheters.

Evaluation of combined late gadolinium-enhancement and functional cardiac magnetic resonance imaging using spiral real-time acquisition

The purpose of the current study was to implement and validate joint real-time acquisition of functional and late gadolinium-enhancement (LGE) cardiac magnetic resonance (MR) images during free breathing. Inversion recovery cardiac real-time images with a temporal resolution of 50 ms were acquired using a spiral trajectory (IR-CRISPI) with a pre-emphasis based on the gradient system transfer function during free breathing. Functional and LGE cardiac MR images were reconstructed using a low-rank plus sparse model. Late gadolinium-enhancement appearance, image quality, and functional parameters of IR-CRISPI were compared with clinical standard balanced steady-state free precession breath-hold techniques in 10 patients. The acquisition of IR-CRISPI in free breathing of the entire left ventricle took 97 s on average. Bland-Altman analysis and Wilcoxon tests showed a higher artifact level for the breath-hold technique (p = 0.003), especially for arrhythmic patients or patients with dyspnea, but an increased noise level for IR-CRISPI of the LGE images (p = 0.01). The estimated transmural extent of the enhancement differed by not more than 25% and did not show a significant bias between the techniques (p = 0.50). The ascertained functional parameters were similar for the breath-hold technique and IR-CRISPI, that is, with a minor, nonsignificant (p = 0.16) mean difference of the ejection fraction of 2.3% and a 95% confidence interval from -4.8% to 9.4%. IR-CRISPI enables joint functional and LGE imaging in free breathing with good image quality but distinctly shorter scan times in comparison with breath-hold techniques.

Real-time cardiac MRI using an undersampled spiral k-space trajectory and a reconstruction based on a variational network

PURPOSE

Cardiac MRI represents the gold standard to determine myocardial function. However, the current clinical standard protocol, a segmented Cartesian acquisition, is time-consuming and can lead to compromised image quality in the case of arrhythmia or dyspnea. In this article, a machine learning-based reconstruction of undersampled spiral k-space data is presented to enable free breathing real-time cardiac MRI with good image quality and short reconstruction times.

METHODS

Data were acquired in free breathing with a 2D spiral trajectory corrected by the gradient system transfer function. Undersampled data were reconstructed by a variational network (VN), which was specifically adapted to the non-Cartesian sampling pattern. The network was trained with data from 11 subjects. Subsequently, the imaging technique was validated in 14 subjects by quantifying the difference to a segmented reference acquisition, an expert reader study, and by comparing derived volumes and functional parameters with values obtained using the current clinical gold standard.

RESULTS

The scan time for the entire heart was below 1 min. The VN reconstructed data in about 0.9 s per image, which is considerably shorter than conventional model-based approaches. The VN furthermore performed better than a U-Net and not inferior to a low-rank plus sparse model in terms of achieved image quality. Functional parameters agreed, on average, with reference data.

CONCLUSIONS

The proposed VN method enables real-time cardiac imaging with both high spatial and temporal resolution in free breathing and with short reconstruction time.

Real-Time Spiral CMR Is Superior to Conventional Segmented Cine-Imaging for Left-Ventricular Functional Assessment in Patients with Arrhythmia

(1) Background: Segmented Cartesian Cardiovascular magnetic resonance (CMR) often fails to deliver robust assessment of cardiac function in patients with arrhythmia. We aimed to assess the performance of a tiny golden-angle spiral real-time CMR sequence at 1.5 T for left-ventricular (LV) volumetry in patients with irregular heart rhythm; (2) Methods: We validated the real-time sequence against the standard breath-hold segmented Cartesian sequence in 32 patients, of whom 11 presented with arrhythmia. End-diastolic volume (EDV), end-systolic volume (ESV), stroke volume (SV), and ejection fraction (EF) were assessed. In arrhythmic patients, real-time and standard Cartesian acquisitions were compared against a reference echocardiographic modality; (3) Results: In patients with sinus rhythm, good agreements and correlations were found between the segmented and real-time methods, with only minor, non-significant underestimation of EDV for the real-time sequence (135.95 ± 30 mL vs. 137.15 ± 31, p = 0.164). In patients with arrhythmia, spiral real-time CMR yielded superior image quality to the conventional segmented imaging, allowing for excellent agreement with the reference echocardiographic volumetry. In contrast, in this cohort, standard Cartesian CMR showed significant underestimation of LV-ESV (106.72 ± 63.51 mL vs. 125.47 ± 72.41 mL, p = 0.026) and overestimation of LVEF (42.96 ± 10.81% vs. 39.02 ± 11.72%, p = 0.039); (4) Conclusions: Real-time spiral CMR improves image quality in arrhythmic patients, allowing reliable assessment of LV volumetry.

A data-driven semantic segmentation model for direct cardiac functional analysis based on undersampled radial MR cine series

PURPOSE

Image acquisition and subsequent manual analysis of cardiac cine MRI is time-consuming. The purpose of this study was to train and evaluate a 3D artificial neural network for semantic segmentation of radially undersampled cardiac MRI to accelerate both scan time and postprocessing.

METHODS

A database of Cartesian short-axis MR images of the heart (148,500 images, 484 examinations) was assembled from an openly accessible database and radial undersampling was simulated. A 3D U-Net architecture was pretrained for segmentation of undersampled spatiotemporal cine MRI. Transfer learning was then performed using samples from a second database, comprising 108 non-Cartesian radial cine series of the midventricular myocardium to optimize the performance for authentic data. The performance was evaluated for different levels of undersampling by the Dice similarity coefficient (DSC) with respect to reference labels, as well as by deriving ventricular volumes and myocardial masses.

RESULTS

Without transfer learning, the pretrained model performed moderately on true radial data [maximum number of projections tested, P = 196; DSC = 0.87 (left ventricle), DSC = 0.76 (myocardium), and DSC =0.64 (right ventricle)]. After transfer learning with authentic data, the predictions achieved human level even for high undersampling rates (P = 33, DSC = 0.95, 0.87, and 0.93) without significant difference compared with segmentations derived from fully sampled data.

CONCLUSION

A 3D U-Net architecture can be used for semantic segmentation of radially undersampled cine acquisitions, achieving a performance comparable with human experts in fully sampled data. This approach can jointly accelerate time-consuming cine image acquisition and cumbersome manual image analysis.