Correction of phase offset errors in main pulmonary artery flow quantification

Pulmonary 4D-flow MRI imaging in landrace pigs under rest and stress

4D-flow MRI is a promising technique for assessing vessel hemodynamics. However, its utilization is currently limited by the lack of reference values, particularly for pulmonary vessels. In this work, we have analysed flow and velocity in the pulmonary trunk (PT), left and right pulmonary arteries (LPA and RPA, respectively) in Landrace pigs at both rest and stress through the software MEVISFlow. Nine healthy Landrace pigs were acutely instrumented closed-chest and transported to the CMR facility for evaluation. After rest measurements, dobutamine was administered to achieve a 25% increase in heart rate compared to rest. 4D-flow MRI images have been analysed through MEVISFlow by two independent observers. Inter- and intra-observer reproducibility was quantified using intraclass correlation coefficient. A significant difference between rest and stress regarding flow and velocity in all the pulmonary vessels was observed. Mean flow increased 55% in PT, 75% in LPA and 40% in RPA. Mean peak velocity increased 55% in PT, 75% in LPA and 66% in RPA. A good-to-excellent reproducibility was observed in rest and stress for flow measurements in all three arteries. An excellent reproducibility for velocity was found in PT at rest and stress, a good one for LPA and RPA at rest, while poor reproducibility was found at stress. The current study showed that pulmonary flow and velocity assessed through 4D-flow MRI follow the physiological alterations during cardiac cycle and after stress induced by dobutamine. A clinical translation to assess pulmonary diseases with 4D-flow MRI under stress conditions needs investigation.

Interrelation Between Cerebrospinal Fluid Pressure, Intracranial Morphology and Venous Hemodynamics Studied by 4D Flow MRI

PURPOSE

To quantify the effects of CSF pressure alterations on intracranial venous morphology and hemodynamics in idiopathic intracranial hypertension (IIH) and spontaneous intracranial hypotension (SIH) and assess reversibility when the underlying cause is resolved.

METHODS

We prospectively examined venous volume, intracranial venous blood flow and velocity, including optic nerve sheath diameter (ONSD) as a noninvasive surrogate of CSF pressure changes in 11 patients with IIH, 11 age-matched and sex-matched healthy controls and 9 SIH patients, before and after neurosurgical closure of spinal dural leaks. We applied multiparametric MRI including 4D flow MRI, time-of-flight (TOF) and T2-weighted half-Fourier acquisition single-shot turbo-spin echo (HASTE).

RESULTS

Sinus volume overlapped between groups at baseline but decreased after treatment of intracranial hypotension (p = 0.067) along with a significant increase of ONSD (p = 0.003). Blood flow in the middle and dorsal superior sagittal sinus was remarkably lower in patients with higher CSF pressure (i.e., IIH versus controls and SIH after CSF leak closure) but blood flow velocity was comparable cross-sectionally between groups and longitudinally in SIH.

CONCLUSION

We were able to demonstrate the interaction of CSF pressure, venous volumetry, venous hemodynamics and ONSD using multiparametric brain MRI. Closure of CSF leaks in SIH patients resulted in symptoms suggestive of increased intracranial pressure and caused a subsequent decrease of intracranial venous volume and of blood flow within the superior sagittal sinus while ONSD increased. In contrast, blood flow parameters from 4D flow MRI did not discriminate IIH, SIH and controls as hemodynamics at baseline overlapped at most vessel cross-sections.

Correction of phase offset errors and quantification of background noise, signal-to-noise ratio, and encoded-displacement uncertainty on DENSE MRI for kinematics of the descending thoracic and abdominal aorta

Displacement encoding with stimulated echoes (DENSE) MRI is a phase contrast technique that allows the encoding of tissue displacement into the phase of the magnetic resonance signal. Recent developments in this technique allow the imaging of relatively thin structures such as the aortic wall. Quantifying background noise associated to DENSE MRI is required to assess the uncertainty of derived displacement measurements and for the design and implementation of adequate noise-reduction techniques. Although noise and error management of cardiac DENSE MRI has been previously studied, developments for aortic applications are scarce. Herein, we evaluate the noise and uncertainty of DENSE MRI scans at three different locations along the descending aorta: the distal aortic arch (DAA), the descending thoracic aorta (DTA), and infrarenal abdominal aorta (IAA). Additionally, we analyze three datasets from in vitro validation experiments with polyvinyl alcohol phantoms. We implement and evaluate the effectiveness of an offset-error correction algorithm and noise filtering techniques on DENSE MRI for aortic motion applications. Our results show that the phase signal of pixels composing the static background was normally distributed, centered on average at 0.003 ± 0.02 rad and - 0.02 ± 0.024 rad for each phase directions, suggesting that background noise is random, isotropic, and DENSE MRI has little offset errors. However, background signal noise significantly increased with elapsed time of the cardiac cycle; and was spatially heterogeneous consistently increased towards the anterior space. Background noise showed no significant differences between the 3 aortic locations and the in vitro experiments. However, SNR depended on the displacement of the region of interest, in consequence it was found significantly larger at DAA (16.7 ± 8.5, p = 0.003) and DTA (15.4 ± 7.6, p = 0.008) than at the IAA (8.0 ± 4.1), but not significantly different than the SNR of in vitro experiments (8.0 ± 3.7), and had an overall average of 13 ± 7. The applied methods significantly reduced the offset error and effect of noise on the estimation of encoded displacements. Finally, this analysis suggests that the implemented DENSE MRI protocol is adequate to assess the motion of healthy human aortas. However, the relative effect of noise increased considerably on the analysis of an ageing or diseased aortas with impaired mobility, calling for further analyses on pathologically stiffened aortas.

Flow evaluation software for four-dimensional flow MRI: a reliability and validation study

PURPOSE

Four-dimensional time-resolved phase-contrast cardiovascular magnetic resonance imaging (4D flow MRI) enables blood flow quantification in multiple vessels, which is crucial for patients with congenital heart disease (CHD). We investigated net flow volumes in the ascending aorta and pulmonary arteries by four different postprocessing software packages for 4D flow MRI in comparison with 2D cine phase-contrast measurements (2D PC).

MATERIAL AND METHODS

4D flow and 2D PC datasets of 47 patients with biventricular CHD (median age 16, range 0.6-52 years) were acquired at 1.5 T. Net flow volumes in the ascending aorta, the main, right, and left pulmonary arteries were measured using four different postprocessing software applications and compared to offset-corrected 2D PC data. Reliability of 4D flow postprocessing software was assessed by Bland-Altman analysis and intraclass correlation coefficient (ICC). Linear regression of internal flow controls was calculated. Interobserver reproducibility was evaluated in 25 patients.

RESULTS

Correlation and agreement of flow volumes were very good for all software compared to 2D PC (ICC ≥ 0.94; bias ≤ 5%). Internal controls were excellent for 2D PC (r ≥ 0.95, p < 0.001) and 4D flow (r ≥ 0.94, p < 0.001) without significant difference of correlation coefficients between methods. Interobserver reliability was good for all vendors (ICC ≥ 0.94, agreement bias < 8%).

CONCLUSION

Haemodynamic information from 4D flow in the large thoracic arteries assessed by four commercially available postprocessing applications matches routinely performed 2D PC values. Therefore, we consider 4D flow MRI-derived data ready for clinical use in patients with CHD.

Temporal evolution of mechanical stimuli from vascular remodeling in response to the severity and duration of aortic coarctation in a preclinical model

Coarctation of the aorta (CoA) is one of the most common congenital cardiovascular diseases. CoA patients frequently undergo surgical repair, but hypertension (HTN) is still common. The current treatment guideline has revealed irreversible changes in structure and function, yet revised severity guidelines have not been proposed. Our objective was to quantify temporal alterations in mechanical stimuli and changes in arterial geometry in response to the range of CoA severities and durations (i.e. age of treatment) seen clinically. Rabbits were exposed to CoA resulting in peak-to-peak blood pressure gradient (BPGpp) severities of ≤ 10, 10-20, and ≥ 20 mmHg for a duration of ~ 1, 3, or 20 weeks using permanent, dissolvable, and rapidly dissolvable sutures. Elastic moduli and thickness were estimated from imaging and longitudinal fluid-structure interaction (FSI) simulations were conducted at different ages using geometries and boundary conditions from experimentally measured data. Mechanical stimuli were characterized including blood flow velocity patterns, wall tension, and radial strain. Experimental results show vascular alternations including thickening and stiffening proximal to the coarctation with increasing severity and/or duration of CoA. FSI simulations indicate wall tension in the proximal region increases markedly with coarctation severity. Importantly, even mild CoA induced stimuli for remodeling that exceeds values seen in adulthood if not treated early and using a BPGpp lower than the current clinical threshold. The findings are aligned with observations from other species and provide some guidance for the values of mechanical stimuli that could be used to predict the likelihood of HTN in human patients with CoA.

A Fluid-Structure Interaction Analysis of Blood Clot Motion in a Branch of Pulmonary Arteries

INTRODUCTION

Pulmonary embolism (PE) is one of the most prevalent diseases amid hospitalized patients taking many people's lives annually. This phenomenon, however, has not been investigated via numerical simulations.

METHODS

In this study, an image-based model of pulmonary arteries has been constructed from a 44-year-old man's computed tomography images. The fluid-structure interaction method was used to simulate the motion of the blood clot. In this regard, Navier-Stokes equations, as the governing equations, have been solved in an arbitrary Lagrangian-Eulerian (ALE) formulation.

RESULTS

According to our results, the velocity of visco-hyperelastic model of the emboli was relatively higher than the emboli with hyperelastic model, despite their similar behavioral pattern. The stresses on the clot were also investigated and showed that the blood clot continuously sustained stresses greater than 165 Pa over an about 0.01 s period, which can cause platelets to leak and make the clot grow or tear apart.

CONCLUSIONS

It could be concluded that in silico analysis of the cardiovascular diseases initiated from clot motion in blood flow is a valuable tool for a better understanding of these phenomena.

Accelerated two-dimensional phase-contrast for cardiovascular MRI using deep learning-based reconstruction with complex difference estimation

PURPOSE

To develop and validate a deep learning-based reconstruction framework for highly accelerated two-dimensional (2D) phase contrast (PC-MRI) data with accurate and precise quantitative measurements.

METHODS

We propose a modified DL-ESPIRiT reconstruction framework for 2D PC-MRI, comprised of an unrolled neural network architecture with a Complex Difference estimation (CD-DL). CD-DL was trained on 155 fully sampled 2D PC-MRI pediatric clinical datasets. The fully sampled data ( n = 29 $$ n=29 $$ ) was retrospectively undersampled (6-11 × $$ \times $$ ) and reconstructed using CD-DL and a parallel imaging and compressed sensing method (PICS). Measurements of peak velocity and total flow were compared to determine the highest acceleration rate that provided accuracy and precision within ± 5 % $$ \pm 5\% $$ . Feasibility of CD-DL was demonstrated on prospectively undersampled datasets acquired in pediatric clinical patients ( n = 5 $$ n=5 $$ ) and compared to traditional parallel imaging (PI) and PICS.

RESULTS

The retrospective evaluation showed that 9 × $$ \times $$ accelerated 2D PC-MRI images reconstructed with CD-DL provided accuracy and precision (bias, [95 % $$ \% $$ confidence intervals]) within ± 5 % $$ \pm 5\% $$ . CD-DL showed higher accuracy and precision compared to PICS for measurements of peak velocity (2.8 % $$ \% $$ [ - 2 . 9 $$ -2.9 $$ , 4.5] vs. 3.9 % $$ \% $$ [ - 11 . 0 $$ -11.0 $$ , 4.9]) and total flow (1.8 % $$ \% $$ [ - 3 . 9 $$ -3.9 $$ , 3.4] vs. 2.9 % $$ \% $$ [ - 7 . 1 $$ -7.1 $$ , 6.9]). The prospective feasibility study showed that CD-DL provided higher accuracy and precision than PICS for measurements of peak velocity and total flow.

CONCLUSION

In a retrospective evaluation, CD-DL produced quantitative measurements of 2D PC-MRI peak velocity and total flow with ≤ 5 % $$ \le 5\% $$ error in both accuracy and precision for up to 9 × $$ \times $$ acceleration. Clinical feasibility was demonstrated using a prospective clinical deployment of our 8 × $$ \times $$ undersampled acquisition and CD-DL reconstruction in a cohort of pediatric patients.

Fully automated background phase correction using M-estimate SAmple consensus (MSAC)-Application to 2D and 4D flow

PURPOSE

Flow quantification by phase-contrast MRI is hampered by spatially varying background phase offsets. Correction performance by polynomial regression on stationary tissue may be affected by outliers such as wrap-around or constant flow. Therefore, we propose an alternative, M-estimate SAmple Consensus (MSAC) to reject outliers, and improve and fully automate background phase correction.

METHODS

The MSAC technique fits polynomials to randomly drawn small samples from the image. Over several trials, it aims to find the best consensus set of valid pixels by rejecting outliers to the fit and minimizing the residuals of the remaining pixels. The robustness of MSAC to its few parameters was investigated and verified using third-order polynomial correction fits on a total of 118 2D flow (97 with wrap-around) and 18 4D flow data sets (14 with wrap-around), acquired at 1.5 T and 3 T. Background phase was compared with standard stationary correction and phantom correction. Pulmonary/systemic flow ratios in 2D flow were derived, and exemplary 4D flow analysis was performed.

RESULTS

The MSAC technique is robust over a range of parameter choices, and a unique set of parameters is suitable for both 2D and 4D flow. In 2D flow, phase errors were significantly reduced by MSAC compared with stationary correction (p = 0.005), and stationary correction shows larger errors in pulmonary/systemic flow ratios compared with MSAC. In 4D flow, MSAC shows similar performance as stationary correction.

CONCLUSIONS

The MSAC method provides fully automated background phase correction to 2D and 4D flow data and shows improved robustness over stationary correction, especially with outliers present.

Whole-brain mapping of mouse CSF flow via HEAP-METRIC phase-contrast MRI

Purpose

CSF plays important roles in clearing brain waste and homeostasis. However, mapping whole‐brain CSF flow in the rodents is difficult, primarily due to its assumed very low velocity. Therefore, we aimed to develop a novel phase‐contrast MRI method to map whole‐brain CSF flow in the mouse brain.

Methods

A novel generalized Hadamard encoding–based multi‐band scheme (dubbed HEAP‐METRIC, Hadamard Encoding APproach of Multi‐band Excitation for short TR Imaging aCcelerating) using complex Hadamard matrix was developed and incorporated into conventional phase contrast (PC)‐MRI to significantly increase SNR.

Results

Slow flow phantom imaging validated HEAP‐METRIC PC‐MRI’s ability to achieve fast and accurate mapping of slow flow velocities (~10² µm/s). With the SNR gain afforded by HEAP‐METRIC scheme, high‐resolution (0.08 × 0.08 mm in‐plane resolution and 36 0.4 mm slices) PC‐MRI was completed in 21 min for whole‐brain CSF flow mapping in the mouse. Using this novel method, we provide the first report of whole‐brain CSF flow in the awake mouse brain with an average flow velocity of ~200 µm/s. Furthermore, HEAP‐METRIC PC‐MRI revealed CSF flow was reduced by isoflurane anesthesia, accompanied by reduction of glymphatic function as measured by dynamic contrast‐enhanced MRI.

Conclusion

We developed and validated a generalized HEAP‐METRIC PC‐MRI for mapping low velocity flow. With this method, we have achieved the first whole‐brain mapping of awake mouse CSF flow and have further revealed that anesthesia reduces CSF flow velocity.

A subject-specific assessment of measurement errors and their correction in cerebrospinal fluid velocity maps using 4D flow MRI

PURPOSE

Phase-contrast MRI (PC-MRI) of cerebrospinal fluid (CSF) velocity is used to evaluate the characteristics of intracranial diseases, such as normal-pressure hydrocephalus (NPH). Nevertheless, PC-MRI has several potential error sources, with eddy-current-based phase offset error being non-negligible in CSF measurement. In this study, we assess the measurement error of CSF velocity maps obtained using 4D flow MRI and evaluate correction methods.

METHODS

CSF velocity maps of 10 patients with NPH were acquired using 4D flow MRI (velocity-encoding = 5 cm/s). Distributed phase offset error was estimated for a whole 3D background field by polynomial fitting using robust regression analysis. This estimated phase offset error was then used to correct the CSF velocity maps. The estimated error profiles were compared with those obtained using an existing 2D correction approach involving local background information near the region of interest.

RESULTS

The residual standard error of the polynomial fitting against the phase offset error extracted from the measured velocities was within 0.2 cm/s. The spatial dependencies of the phase offset errors showed similar tendencies in all cases, but sufficient differences in these values were found to indicate requirement of velocity correction. Differences of the estimated errors among other correction approaches were in the order of 10^-2 cm/s, and the estimated errors were in good agreement with those obtained using existing approaches.

CONCLUSION

Our method is capable of estimating the measurement error of CSF velocity maps obtained from 4D flow MRI and provides quantitatively reasonable characteristics for the main CSF profile in the cerebral aqueduct in patients with NPH.

Importance of through-plane heart motion correction for the assessment of aortic regurgitation severity using phase contrast magnetic resonance imaging

PURPOSE

To elucidate the influence of through-plane heart motion on the assessment of aortic regurgitation (AR) severity using phase contrast magnetic resonance imaging (PC-MRI).

APPROACH

A patient cohort with chronic AR (n = 34) was examined with PC-MRI. The regurgitant volume (RVol) and fraction (RFrac) were extracted from the PC-MRI data before and after through-plane heart motion correction and was then used for assessment of AR severity.

RESULTS

The flow volume errors were strongly correlated to aortic diameter (R = 0.80, p < 0.001) with median (IQR 25%;75%): 16 (14; 17) ml for diameter>40mm, compared with 9 (7; 10) ml for normal aortic size (p < 0.001). RVol and RFrac were underestimated (uncorrected:64 ± 37 ml and 39 ± 17%; corrected:76 ± 37 ml and 44 ± 15%; p < 0.001) and ~ 20% of the patients received lower severity grade without correction.

CONCLUSION

Through-plane heart motion introduces relevant flow volume errors, especially in patients with aortic dilatation that may result in underestimation of the severity grade in patients with chronic AR.

Luminal enhancement in intracranial aneurysms: fact or feature?-A quantitative multimodal flow analysis

Purpose

Intracranial aneurysm (IA) wall enhancement on post-contrast vessel wall magnetic resonance imaging (VW-MRI) is assumed to be a biomarker for vessel wall inflammation and aneurysm instability. However, the exact factors contributing to enhancement are not yet clarified. This study investigates the relationship between luminal enhancement and intra-aneurysmal flow behaviour to assess the suitability of VW-MRI as a surrogate method for determining quantitative and qualitative flow behaviour in the aneurysm sac.

Methods

VW-MRI signal is measured in the lumen of three patient-specific IA flow models and compared with the intra-aneurysmal flow fields obtained using phase-contrast magnetic resonance imaging (PC-MRI) and computational fluid dynamics (CFD). The IA flow models were supplied with two different time-varying flow regimes.

Results

Overall, the velocity fields acquired using PC-MRI or CFD were in good agreement with the VW-MRI enhancement patterns. Generally, the regions with slow-flowing blood show higher VW-MRI signal intensities, whereas high flow leads to a suppression of the signal. For all aneurysm models, a signal value above three was associated with velocity values below three cm/s.

Conclusion

Regions with lower enhancements have been correlated with the slow and high flow at the same time. Thus, further factors like flow complexity and stability can contribute to flow suppression in addition to the flow magnitude. Nevertheless, VW-MRI can qualitatively assess intra-aneurysmal flow phenomena and estimate the velocity range present in the corresponding region.

Blood Flow Velocity Pulsatility and Arterial Diameter Pulsatility Measurements of the Intracranial Arteries Using 4D PC-MRI

4D phase contrast magnetic resonance imaging (PC-MRI) allows for the visualization and quantification of the cerebral blood flow. A drawback of software that is used to quantify the cerebral blood flow is that it oftentimes assumes a static arterial luminal area over the cardiac cycle. Quantifying the lumen area pulsatility index (aPI), i.e. the change in lumen area due to an increase in distending pressure over the cardiac cycle, can provide insight in the stiffness of the arteries. Arterial stiffness has received increased attention as a predictor in the development of cerebrovascular disease. In this study, we introduce software that allows for measurement of the aPI as well as the blood flow velocity pulsatility index (vPI) from 4D PC-MRI. The internal carotid arteries of seven volunteers were imaged using 7 T MRI. The aPI and vPI measurements from 4D PC-MRI were validated against measurements from 2D PC-MRI at two levels of the internal carotid arteries (C3 and C7). The aPI and vPI computed from 4D PC-MRI were comparable to those measured from 2D PC-MRI (aPI: mean difference: 0.03 (limits of agreement: -0.14 - 0.23); vPI: 0.03 (-0.17-0.23)). The measured blood flow rate for the C3 and C7 segments was similar, indicating that our proposed software correctly captures the variation in arterial lumen area and blood flow velocity that exists along the distal end of the carotid artery. Our software may potentially aid in identifying changes in arterial stiffness of the intracranial arteries caused by pathological changes to the vessel wall.

Velocity quantification in 44 healthy volunteers using accelerated multi-VENC 4D flow CMR

BACKGROUND

To evaluate the feasibility of a k-t accelerated multi-VENC 4D phase contrast flow MRI acquisition of the main heart-surrounding vessels, its benefits over a traditional single-VENC acquisition and to present reference flow and velocity values in a large cohort of volunteers.

METHODS

44 healthy volunteers were examined on a 3 T MRI scanner (Ingenia, Philips, Best, The Netherlands). 4D flow measurements were obtained with a FOV including the aorta and the pulmonary arteries. VENC values were set to 40, 100 and 200 cm/s and unfolded based on an MRI signal model. Unfolded multi-VENC data was compared to the single-VENC with VENC 200 cm/s. Flow and velocity quantification was performed in several regions of interest (ROI) placed in the ascending aorta and in the main pulmonary artery. Conservation of mass analysis was performed for single- and multi-VENC datasets. Values for mean and maximal flow velocity and stroke volume were calculated and compared to the literature.

RESULTS

Mean scan time was 13.8 ± 4 min. Differences between stroke volumes between the ascending aorta and the main pulmonary artery were significantly lower in multi-VENC datasets compared to single-VENC datasets (9.6 ± 7.8 mL vs. 25.4 ± 26.4 mL, p < 0.001). This was also true for differences in stroke volume between up- and downstream ROIs in the ascending aorta and pulmonary artery. Values for mean and maximal velocities and stroke volume were in-line with previous studies. To highlight potential clinical applications two exemplary 4D flow measurements in patients with different pathologies are shown and compared to single-VENC datasets.

CONCLUSIONS

k-t accelerated multi-VENC 4D phase contrast flow MRI acquisition of the great vessels is feasible in a clinically acceptable scan duration. It offers improvements over traditional single-VENC 4D flow, expectedly being valuable when vessels with different flow velocities or complex flow phenomena are evaluated.

Cardiac magnetic resonance systematically overestimates mitral regurgitations by the indirect method

Objective

Cardiac MRI is quickly emerging as the gold standard for assessment of mitral regurgitation, most commonly with the indirect method subtracting forward flow in aorta from volumetric segmentation of the left ventricle. We aimed to investigate how aortic flow measurements with increasing distance from the aortic valve affect calculated mitral regurgitations and whether measurements were influenced by breath-hold regimen.

Methods

Free-breathing and breath-hold phase contrast flows were measured in aorta at valve level, sinotubular (ST) junction, mid-ascending aorta and in the pulmonary trunk. Flow measurements were pairwise compared, and subsequently, after exclusion of patients with visible mitral and tricuspid regurgitations for left-sided and right-sided comparisons, respectively, flow-measured stroke volumes were compared with ventricular volumetric segmentations.

Results

Thirty-nine participants without arrhythmias or structural abnormalities of the large vessels were included. Stroke volumes measured with free-breathing and breath-hold flow decreased equally with increasing distance to the aortic valves (breath-hold flow: aortic valve 105.6±20.8 mL, ST junction 101.5±20.7 mL, mid-ascending aorta 98.1±21.5 mL). After exclusion of atrioventricular regurgitations, stroke volumes determined by volumetric measurements were higher compared with values determined by flow measurements, corresponding to ‘false’ atrioventricular regurgitations of 8.0%±5.8% with flow measured at valve level, 11.6%±5.2% at the ST junction and 15.3%±5.0% at the mid-ascending aorta.

Conclusions

Stroke volumes determined by flow decrease throughout the proximal aorta and are systematically lower than volumetrically measured stroke volumes. The indirect method systematically overestimates mitral regurgitations, especially with increasing distance from the aortic valves.

Reliable phase-contrast flow volume magnetic resonance measurements are feasible without adjustment of the velocity encoding parameter

Purpose: To show that adjustment of velocity encoding (VENC) for phase-contrast (PC) flow volume measurements is not necessary in modern MR scanners with effective background velocity offset corrections. Approach: The independence on VENC was demonstrated theoretically, but also experimentally on dedicated phantoms and on patients with chronic aortic regurgitation ( n = 17 ) and one healthy volunteer. All PC measurements were performed using a modern MR scanner, where the pre-emphasis circuit but also a subsequent post-processing filter were used for effective correction of background velocity offset errors. Results: The VENC level strongly affected the velocity noise level in the PC images and, hence, the estimated peak flow velocity. However, neither the regurgitant blood flow volume nor the mean flow velocity displayed any clinically relevant dependency on the VENC level. Also, the background velocity offset was shown to be close to zero ( < 0.6    cm / s ) for a VENC range of 150 to 500    cm / s , adding no significant errors to the PC flow volume measurement. Conclusions: Our study shows that reliable PC flow volume measurements are feasible without adjustment of the VENC parameter. Without the need for VENC adjustments, the scan time can be reduced for the benefit of the patient.

Super-resolution and denoising of 4D-Flow MRI using physics-Informed deep neural nets

BACKGROUND AND OBJECTIVE

Time resolved three-dimensional phase contrast magnetic resonance imaging (4D-Flow MRI) has been used to non-invasively measure blood velocities in the human vascular system. However, issues such as low spatio-temporal resolution, acquisition noise, velocity aliasing, and phase-offset artifacts have hampered its clinical application. In this research, we developed a purely data-driven method for super-resolution and denoising of 4D-Flow MRI.

METHODS

The flow velocities, pressure, and the MRI image magnitude are modeled as a patient-specific deep neural net (DNN). For training, 4D-Flow MRI images in the complex Cartesian space are used to impose data-fidelity. Physics of fluid flow is imposed through regularization. Creative loss function terms have been introduced to handle noise and super-resolution. The trained patient-specific DNN can be sampled to generate noise-free high-resolution flow images. The proposed method has been implemented using the TensorFlow DNN library and tested on numerical phantoms and validated in-vitro using high-resolution particle image velocitmetry (PIV) and 4D-Flow MRI experiments on transparent models subjected to pulsatile flow conditions.

RESULTS

In case of numerical phantoms, we were able to increase spatial resolution by a factor of 100 and temporal resolution by a factor of 5 compared to the simulated 4D-Flow MRI. There is an order of magnitude reduction of velocity normalized root mean square error (vNRMSE). In case of the in-vitro validation tests with PIV as reference, there is similar improvement in spatio-temporal resolution. Although the vNRMSE is reduced by 50%, the method is unable to negate a systematic bias with respect to the reference PIV that is introduced by the 4D-Flow MRI measurement.

CONCLUSIONS

This work has demonstrated the feasibility of using the readily available machinery of deep learning to enhance 4D-Flow MRI using a purely data-driven method. Unlike current state-of-the-art methods, the proposed method is agnostic to geometry and boundary conditions and therefore eliminates the need for tedious tasks such as accurate image segmentation for geometry, image registration, and estimation of boundary flow conditions. Arbitrary regions of interest can be selected for processing. This work will lead to user-friendly analysis tools that will enable quantitative hemodynamic analysis of vascular diseases in a clinical setting.

The clinical impact of phase offset errors and different correction methods in cardiovascular magnetic resonance phase contrast imaging: a multi-scanner study

BACKGROUND

Cardiovascular magnetic resonance (CMR) phase contrast (PC) flow measurements suffer from phase offset errors. Background subtraction based on stationary phantom measurements can most reliably be used to overcome this inaccuracy. Stationary tissue correction is an alternative and does not require additional phantom scanning. The aim of this study was 1) to compare measurements with and without stationary tissue correction to phantom corrected measurements on different GE Healthcare CMR scanners using different software packages and 2) to evaluate the clinical implications of these methods.

METHODS

CMR PC imaging of both the aortic and pulmonary artery flow was performed in patients on three different 1.5 T CMR scanners (GE Healthcare) using identical scan parameters. Uncorrected, first, second and third order stationary tissue corrected flow measurement were compared to phantom corrected flow measurements, our reference method, using Medis QFlow, Circle cvi42 and MASS software. The optimal (optimized) stationary tissue order was determined per scanner and software program. Velocity offsets, net flow, clinically significant difference (deviation > 10% net flow), and regurgitation severity were assessed.

RESULTS

Data from 175 patients (28 (17-38) years) were included, of which 84% had congenital heart disease. First, second and third order and optimized stationary tissue correction did not improve the velocity offsets and net flow measurements. Uncorrected measurements resulted in the least clinically significant differences in net flow compared to phantom corrected data. Optimized stationary tissue correction per scanner and software program resulted in net flow differences (> 10%) in 19% (MASS) and 30% (Circle cvi42) of all measurements compared to 18% (MASS) and 23% (Circle cvi42) with no correction. Compared to phantom correction, regurgitation reclassification was the least common using uncorrected data. One CMR scanner performed worse and significant net flow differences of > 10% were present both with and without stationary tissue correction in more than 30% of all measurements.

CONCLUSION

Phase offset errors had a significant impact on net flow quantification, regurgitation assessment and varied greatly between CMR scanners. Background phase correction using stationary tissue correction worsened accuracy compared to no correction on three GE Healthcare CMR scanners. Therefore, careful assessment of phase offset errors at each individual scanner is essential to determine whether routine use of phantom correction is necessary.

TRIAL REGISTRATION

Observational Study.

Mitral Valve Flow and Myocardial Motion Assessed with Dual-Echo Dual-Velocity Cardiac MRI

Purpose

To develop a dual-echo phase-contrast (DEPC) MRI approach with which each echo is acquired by using a different velocity sensitivity within one repetition time (TR) and demonstrate the feasibility of this approach to measure transmitral blood flow (E) and myocardial tissue (E _m) velocities.

Materials and Methods

The flow across tubes of known diameter was measured by using the proposed DEPC method and compared with flowmeter measurements and theoretic predictions. Then, with both the DEPC MRI sequence and the conventional single-echo phase-contrast (SEPC) MRI sequence, E, E _m, and E/E _m were measured in six healthy volunteers (mean age, 49 years ± 13 [standard deviation]) and eight patients (mean age, 54 years ± 15) being evaluated for cardiac disease. Differences between the DEPC and conventional SEPC MRI methods were assessed by percent error, Pearson correlation, and Bland-Altman analyses.

Results

Velocities measured in vitro and in vivo by using the SEPC and DEPC MRI approaches were well correlated (r ² > 0.97), with negligible bias (<0.5 cm/sec) and comparable velocity-to-noise ratios. Imaging times were approximately 19% shorter with the DEPC method (TR, 5.7 msec) than with the SEPC method (TR, 2.8 msec ± 4.2) (P < .05).

Conclusion

The proposed DEPC method was sensitive to two velocity regimes within a single TR, resulting in a shorter imaging time compared with the imaging time in conventional SEPC MRI. Preliminary human study results suggest the feasibility of using this approach to estimate E/E _m.Supplemental material is available for this article.© RSNA, 2020.

Evaluation of self-calibrated non-linear phase-contrast correction in pediatric and congenital cardiovascular magnetic resonance imaging

BACKGROUND

The need for background error correction in phase-contrast flow analysis has historically posed a challenge in cardiac magnetic resonance (MR) imaging. While previous studies have shown that phantom correction improves flow measurements, it impedes scanner workflow.

OBJECTIVE

To evaluate the efficacy of self-calibrated non-linear phase-contrast correction on flows in pediatric and congenital cardiac MR compared to phantom correction as the standard.

MATERIALS AND METHODS

We retrospectively identified children who had great-vessel phase-contrast and static phantom sequences acquired between January 2015 and June 2015. We applied a novel correction method to each phase-contrast sequence post hoc. Uncorrected, non-linear, and phantom-corrected flows were compared using intraclass correlation. We used paired t-tests to compare how closely non-linear and uncorrected flows approximated phantom-corrected flows. In children without intra- or extracardiac shunts or significant semilunar valvular regurgitation, we used paired t-tests to compare how closely the uncorrected pulmonary-to-systemic flow ratio (Qp:Qs) and non-linear Qp:Qs approximated phantom-corrected Qp:Qs.

RESULTS

We included 211 diagnostic-quality phase-contrast sequences (93 aorta, 74 main pulmonary artery [MPA], 21 left pulmonary artery [LPA], 23 right pulmonary artery [RPA]) from 108 children (median age 15 years, interquartile range 11-18 years). Intraclass correlation showed strong agreement between non-linear and phantom-corrected flow measurements but also between uncorrected and phantom-corrected flow measurements. Non-linear flow measurements did not more closely approximate phantom-corrected measurements than did uncorrected measurements for any vessel. In 39 children without significant shunting or regurgitation, mean non-linear Qp:Qs (1.07; 95% confidence interval [CI] = 1.01, 1.13) was no closer than mean uncorrected Qp:Qs (1.06; 95% CI = 1.00, 1.13) to mean phantom-corrected Qp:Qs (1.02; 95% CI = 0.98, 1.06).

CONCLUSION

Despite strong agreement between self-calibrated non-linear and phantom correction, cardiac flows and shunt calculations with non-linear correction were no closer to phantom-corrected measurements than those without background correction. However, phantom-corrected flows also demonstrated minimal differences from uncorrected flows. These findings suggest that in the current era, more accurate phase-contrast flow measurements might limit the need for background correction. Further investigation of the clinical impact and optimal methods of background correction in the pediatric and congenital cardiac population is needed.

Stationary tissue background correction increases the precision of clinical evaluation of intra-cardiac shunts by cardiovascular magnetic resonance

We aimed to evaluate the clinical utility of stationary tissue background phase correction for affecting precision in the measurement of Qp/Qs by cardiovascular magnetic resonance (CMR). We enrolled consecutive patients (n = 91) referred for CMR at 1.5T without suspicion of cardiac shunt, and patients (n = 10) with verified cardiac shunts in this retrospective study. All patients underwent phase contrast flow quantification in the ascending aorta and pulmonary trunk. Flow was quantified using two semi-automatic software platforms (SyngoVia VA30, Vendor 1; Segment 2.0R4534, Vendor 2). Measurements were performed both uncorrected and corrected for linear (Vendor 1 and Vendor 2) or quadratic (Vendor 2) background phase. The proportion of patients outside the normal range of Qp/Qs was compared using the McNemar’s test. Compared to uncorrected measurements, there were fewer patients with a Qp/Qs outside the normal range following linear correction using Vendor 1 (10% vs 18%, p < 0.001), and Vendor 2 (10% vs 18%, p < 0.001), and following quadratic correction using Vendor 2 (7% vs 18%, p < 0.001). No patient with known shunt was reclassified as normal following stationary background correction. Therefore, we conclude that stationary tissue background correction reduces the number of patients with a Qp/Qs ratio outside the normal range in a consecutive clinical population, while simultaneously not reclassifying any patient with known cardiac shunts as having a normal Qp/Qs. Stationary tissue background correction may be used in clinical patients to increase diagnostic precision.

Assessment of mitral valve regurgitation by cardiovascular magnetic resonance imaging

Mitral regurgitation (MR) is a common valvular heart disease and is the second most frequent indication for heart valve surgery in Western countries. Echocardiography is the recommended first-line test for the assessment of valvular heart disease, but cardiovascular magnetic resonance imaging (CMR) provides complementary information, especially for assessing MR severity and to plan the timing of intervention. As new CMR techniques for the assessment of MR have arisen, standardizing CMR protocols for research and clinical studies has become important in order to optimize diagnostic utility and support the wider use of CMR for the clinical assessment of MR. In this Consensus Statement, we provide a detailed description of the current evidence on the use of CMR for MR assessment, highlight its current clinical utility, and recommend a standardized CMR protocol and report for MR assessment.

In this Consensus Statement, Garg and colleagues describe the current evidence on the use of cardiovascular magnetic resonance imaging for the assessment of mitral regurgitation, highlight its current clinical utility, and recommend a standardized imaging protocol and report.

Comparison of 4D Flow MRI to 2D Flow MRI in the pulmonary arteries in healthy volunteers and patients with pulmonary hypertension

Purpose

4D and 2D phase-contrast MRI (2D Flow MRI, 4D Flow MRI, respectively) are increasingly being used to noninvasively assess pulmonary hypertension (PH). The goals of this study were i) to evaluate whether established quantitative parameters in 2D Flow MRI associated with pulmonary hypertension can be assessed using 4D Flow MRI; ii) to compare results from 4D Flow MRI on a digital broadband 3T MR system with data from clinically established MRI-techniques as well as conservation of mass analysis and phantom correction and iii) to elaborate on the added value of secondary flow patterns in detecting PH.

Methods

11 patients with PH (4f, 63 ± 16y), 15 age-matched healthy volunteers (9f, 56 ± 11y), and 20 young healthy volunteers (13f, 23 ± 2y) were scanned on a 3T MR scanner (Philips Ingenia). Subjects were examined with a 4D Flow, a 2D Flow and a bSSFP sequence. For extrinsic comparison, quantitative parameters measured with 4D Flow MRI were compared to i) a static phantom, ii) 2D Flow acquisitions and iii) stroke volume derived from a bSSFP sequence. For intrinsic comparison conservation of mass-analysis was employed. Dedicated software was used to extract various flow, velocity, and anatomical parameters. Visualization of blood flow was performed to detect secondary flow patterns.

Results

Overall, there was good agreement between all techniques, 4D Flow results revealed a considerable spread. Data improved after phantom correction. Both 4D and 2D Flow MRI revealed concordant results to differentiate patients from healthy individuals, especially based on values derived from anatomical parameters. The visualization of a vortex, indicating the presence of PH was achieved in 9 /11 patients and 2/35 volunteers.

Discussion

This study confirms that quantitative parameters used for characterizing pulmonary hypertension can be gathered using 4D Flow MRI within clinically reasonable limits of agreement. Despite its unfavorable spatial and lesser temporal resolution and a non-neglible spread of results, the identification of diseased study participants was possible.

Data Quality and Optimal Background Correction Order of Respiratory-Gated k-Space Segmented Spoiled Gradient Echo (SGRE) and Echo Planar Imaging (EPI)-Based 4D Flow MRI

Background

A reduction in scan time of 4D Flow MRI would facilitate clinical application. A recent study indicates that echo‐planar imaging (EPI) 4D Flow MRI allows for a reduction in scan time and better data quality than the recommended k‐space segmented spoiled gradient echo (SGRE) sequence. It was argued that the poor data quality of SGRE was related to the nonrecommended absence of respiratory motion compensation. However, data quality can also be affected by the background offset compensation.

Purpose

To compare the data quality of respiratory motion‐compensated SGRE and EPI 4D Flow MRI and their dependence on background correction (BC) order.

Study Type

Retrospective.

Subjects

Eighteen healthy subjects (eight female, mean age 32 ± 5 years).

Field Strength and Sequence

1.5 T. [Correction added on July 26, 2019, after first online publication: The preceding field strength was corrected.] SGRE and EPI‐based 4D Flow MRI.

Assessment

Data quality was investigated visually and by comparing flows through the cardiac valves and aorta. Measurements were obtained from transvalvular flow and pathline analysis.

Statistical Tests

Linear regression and Bland–Altman analysis were used. Wilcoxon test was used for comparison of visual scoring. Student's t‐test was used for comparison of flow volumes.

Results

No significant difference was found by visual inspection (P = 0.08). Left ventricular (LV) flows were strongly and very strongly associated with SGRE and EPI, respectively (R² = 0.86–0.94 SGRE; 0.71–0.79 EPI, BC0–4). LV and right ventricular (RV) outflows and LV pathline flows were very strongly associated (R² = 0.93–0.95 SGRE; 0.88–0.91 EPI, R² = 0.91–0.95 SGRE; 0.91–0.93 EPI, BC1–4). EPI LV outflow was lower than the short‐axis‐based stroke volume. EPI RV outflow and proximal descending aortic flow were lower than SGREs.

Data Conclusion

Both sequences yielded good internal data consistency when an adequate background correction was applied. Second and first BC order were considered sufficient for transvalvular flow analysis in SGRE and EPI, respectively. Higher BC orders were preferred for particle tracing.

Level of Evidence 4

Technical Efficacy Stage 1

J. Magn. Reson. Imaging 2020;51:885–896.

In-vivo validation of interpolation-based phase offset correction in cardiovascular magnetic resonance flow quantification: a multi-vendor, multi-center study

BACKGROUND

A velocity offset error in phase contrast cardiovascular magnetic resonance (CMR) imaging is a known problem in clinical assessment of flow volumes in vessels around the heart. Earlier studies have shown that this offset error is clinically relevant over different systems, and cannot be removed by protocol optimization. Correction methods using phantom measurements are time consuming, and assume reproducibility of the offsets which is not the case for all systems. An alternative previously published solution is to correct the in-vivo data in post-processing, interpolating the velocity offset from stationary tissue within the field-of-view. This study aims to validate this interpolation-based offset correction in-vivo in a multi-vendor, multi-center setup.

METHODS

Data from six 1.5 T CMR systems were evaluated, with two systems from each of the three main vendors. At each system aortic and main pulmonary artery 2D flow studies were acquired during routine clinical or research examinations, with an additional phantom measurement using identical acquisition parameters. To verify the phantom acquisition, a region-of-interest (ROI) at stationary tissue in the thorax wall was placed and compared between in-vivo and phantom measurements. Interpolation-based offset correction was performed on the in-vivo data, after manually excluding regions of spatial wraparound. Correction performance of different spatial orders of interpolation planes was evaluated.

RESULTS

A total of 126 flow measurements in 82 subjects were included. At the thorax wall the agreement between in-vivo and phantom was - 0.2 ± 0.6 cm/s. Twenty-eight studies were excluded because of a difference at the thorax wall exceeding 0.6 cm/s from the phantom scan, leaving 98. Before correction, the offset at the vessel as assessed in the phantom was - 0.4 ± 1.5 cm/s, which resulted in a - 5 ± 16% error in cardiac output. The optimal order of the interpolation correction plane was 1st order, except for one system at which a 2nd order plane was required. Application of the interpolation-based correction revealed a remaining offset velocity of 0.1 ± 0.5 cm/s and 0 ± 5% error in cardiac output.

CONCLUSIONS

This study shows that interpolation-based offset correction reduces the offset with comparable efficacy as phantom measurement phase offset correction, without the time penalty imposed by phantom scans.

TRIAL REGISTRATION

The study was registered in The Netherlands National Trial Register (NTR) under TC 4865 . Registered 19 September 2014. Retrospectively registered.

Left ventricular blood flow patterns at rest and under dobutamine stress in healthy pigs

Intracardiac blood flow patterns are affected by the morphology of cardiac structures and are set up to support the heart's pump function. Exercise affects contractility and chamber size as well as pre- and afterload. The aim of this study was to test the feasibility of four-dimensional phase contrast cardiovascular MRI under pharmacological stress and to study left ventricular blood flow under stress. 4D flow data were successfully acquired and analysed in 12 animals. During dobutamine infusion, heart rate and ejection fraction increased (82 ± 5 bpm versus 124 ± 3 bpm/46 ± 9% versus 65 ± 7%; both p < 0.05). A decrease in left ventricular end-diastolic volume (72 ± 14 mL versus 55 ± 8 mL; p < 0.05) and end-systolic volume (40 ± 15 mL versus 19 ± 6 mL; p < 0.05) but no change in stroke volume were observed. Trans-mitral diastolic inflow velocity increased under dobutamine and the trajectory of inflowing blood was directed towards the anterior septum with increased inflow angle (26 ± 5°) when compared with controls (15 ± 2°). In 5/6 animals undergoing stress diastolic vortices developed later, and in 3/6 animals vortices collapsed earlier with significantly smaller cross-sectional area during diastole. The vorticity index was not affected. Under the stress condition direct flow (% ejection within the next heart beat) increased from 43 ± 6% to 53 ± 8%. 4D MRI blood flow acquisition and analysis are feasible in pig hearts under dobutamine-induced stress. Flow patterns characterized by high blood velocity and antero-septally oriented diastolic inflow as well as decreased ventricular volumes are unfavourable conditions for diastolic vortex development under pharmacological stress, and cardiac output is increased by a rise in heart rate and directly ejected left ventricular blood volume.

A method to correct background phase offset for phase-contrast MRI in the presence of steady flow and spatial wrap-around artifact

PURPOSE

Background phase offsets in phase-contrast MRI are often corrected using polynomial regression; however, correction performance degrades when temporally invariant outliers such as steady flow or spatial wrap-around artifact are present. We describe and validate an iterative method called automatic rejection of temporally invariant outliers (ARTO), which excludes these outliers from the fitting process.

METHODS

The ARTO method iteratively removes pixels with large polynomial regression errors analyzed by a Gaussian mixture model fitting of the residual distribution. A total of 150 trials of a simulated phantom (75 with wrap-around artifact) and 125 phase-contrast MRI cines from 22 healthy subjects (48 with wrap-around artifact) were used for validation. Background phase offsets were corrected using second-order weighted regularized least squares (WRLS) with and without ARTO. Flow volumes after WRLS and WRLS+ARTO corrections were compared with the known truth (phantom) and stationary phantom reference (in vivo) using Bland-Altman analysis. The ratio between the pulmonary flow and the systemic flow was also computed in a subset of 6 subjects.

RESULTS

In the simulated phantom, compared with WRLS and no correction, correction with WRLS+ARTO produced superior agreement in volumetric flow quantification with the known truth. In vivo, WRLS+ARTO also produced superior agreement with stationary phantom-corrected volumetric flow compared with WRLS and no correction. In data sets with wrap-around artifact, WRLS produced significantly larger variance in the pulmonary flow and systemic flow ratio than stationary phantom correction (P = .0008).

CONCLUSION

The proposed method provides automatic exclusion of temporally invariant outliers and produces flow quantification results comparable to stationary phantom correction.

Visual Analysis of Aneurysm Data using Statistical Graphics

This paper presents a framework to explore multi-field data of aneurysms occurring at intracranial and cardiac arteries by using statistical graphics. The rupture of an aneurysm is often a fatal scenario, whereas during treatment serious complications for the patient can occur. Whether an aneurysm ruptures or whether a treatment is successful depends on the interaction of different morphological such as wall deformation and thickness, and hemodynamic attributes like wall shear stress and pressure. Therefore, medical researchers are very interested in better understanding these relationships. However, the required analysis is a time-consuming process, where suspicious wall regions are difficult to detect due to the time-dependent behavior of the data. Our proposed visualization framework enables medical researchers to efficiently assess aneurysm risk and treatment options. This comprises a powerful set of views including 2D and 3D depictions of the aneurysm morphology as well as statistical plots of different scalar fields. Brushing and linking aids the user to identify interesting wall regions and to understand the influence of different attributes on the aneurysm's state. Moreover, a visual comparison of pre- and post-treatment as well as different treatment options is provided. Our analysis techniques are designed in collaboration with domain experts, e.g., physicians, and we provide details about the evaluation.

Lower cerebral blood flow in subjects with Alzheimer's dementia, mild cognitive impairment, and subjective cognitive decline using two-dimensional phase-contrast magnetic resonance imaging

Introduction

In this cross-sectional study, we aimed to detect differences in cerebral blood flow (CBF) between subjects with Alzheimer's disease (AD), mild cognitive impairment (MCI), and subjective cognitive decline (SCD), using two-dimensional phase-contrast magnetic resonance imaging.

Methods

We included 74 AD patients (67 years, 51% female), 36 MCI patients (66 years, 33% female), and 62 patients with SCD (60 years, 32% female) from the Amsterdam Dementia Cohort. Patients with SCD are those who visited the memory clinic with subjective cognitive complaints without objective cognitive impairment. Whole-brain CBF (mL/100 g/min) was calculated using total volume flow measured with two-dimensional phase-contrast magnetic resonance imaging and normalized for brain volume.

Results

Mean CBF values (SD) were lower in AD compared to SCD (age and sex adjusted 70 ± 26 vs. 82 ± 24 mL/100 g/min, P < .05). Mean CBF values of MCI were comparable to AD. Across clinical groups, lower CBF was associated with lower scores on the Mini–Mental State Examination (age and sex adjusted stβ = 0.19 per mL/100 g/min; P = .02).

Discussion

Lower whole-brain CBF is seen in AD patients compared to SCD patients and is associated with worse cognitive function.

•

The study consisted of a large sample of patients with AD, MCI, and controls.

•

CBF measured with 2D PC MRI differed between AD patients and controls.

•

Lower CBF was associated with worse cognitive function measured with MMSE.

•

2D PC MRI may be used as a marker for disease severity in a memory clinic.

Voxel-by-voxel 4D flow MRI-based assessment of regional reverse flow in the aorta

BACKGROUND

Complex and reverse flow in the aorta has been implicated in aneurysm development and stroke via retrograde embolization.

PURPOSE

To evaluate global and regional differences between standard 2D plane-based and volumetric voxel-based quantification of regional forward/reverse flow, and reverse flow fraction (RFF) in the aorta.

STUDY TYPE

Retrospective.

SUBJECTS

In all, 35 subjects: 10 healthy controls (age: 57 ± 7 years, nine male), nine patients without aortic valve regurgitation (AR) (age: 63 ± 10 years, seven male), six patients with mild AR (age: 66 ± 6 years, five male), and 10 with moderate or severe AR (age: 60 ± 16 years, eight male).

FIELD STRENGTH/SEQUENCE

4D flow MRI (3T and 1.5T) was employed to acquire 3D blood flow velocities with entire thoracic aorta in all subjects.

ASSESSMENT

Data analysis included standard 2D plane-based quantification of forward/reverse flow, and RFF-plane. In addition, a new semiautomatic workflow based on 3D segmentation and extraction of an aorta centerline was developed for voxel-by-voxel visualization (forward/reverse flow and RFF-voxel maps) and quantification of regional voxel-by-voxel forward/reverse flow in the entire thoracic aorta.

STATISTICAL TESTS

Kruskal-Wallis tests were performed to test for differences between groups. A two-sample t-test or Wilcoxon rank sum test was used to compare voxel-based and plane-based results.

RESULTS

Semiautomatic plane-based analysis showed excellent agreement with standard manual plane-based analysis for net flow and RFF-plane (RFF-plane: y = 0.99x-0.0, net flow: y = 1.00x-0.21, R > 0.99, P < 0.0001). Voxel-by-voxel maps demonstrated marked regional flow reversal in the ascending aorta in all patients and RFF-voxel was significantly increased (P < 0.001) compared to RFF-plane for all four groups, with the most pronounced differences for mild AR (18.0 ± 15.2% vs. 4.7 ± 5.4%). Voxel-based flow and RFF-voxel along the aorta showed areas with marked regional flow reversal (eg, vortex flow) compared to plane-based analysis.

DATA CONCLUSION

Voxel-based analysis demonstrated regional flow reversal that was not detected by plane-based analysis.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;47:1276-1286.

Pixel-wise derivation of pulmonary regurgitation index could influence clinical decision: A phase-contrast MR imaging study on patients with repaired tetralogy of Fallot

PURPOSE

Regurgitant fraction (RF) measured from 2D phase-contrast MRI has been used as a standard to quantitate the degree of pulmonary regurgitation and serves as a determinant indicator of prognosis for tetralogy of Fallot after surgical repair. This study demonstrated the potential underestimate of RF using the conventional definition and its impact on clinical decision when backward flow occurred during systolic periods.

METHODS

Quantitative flow parameters, including forward flow volume (FFV), backward flow volume (BFV), and RF were estimated by two approaches: One derived from conventional ROI-averaged curve (bulk quantity) and the other in a pixel-wise manner to spatially reflect inhomogeneous flow profile (pixel-wise quantity). Eccentricity at systolic peak (Ecc_sys) was adopted as an index reflecting spatial flow inhomogeneity.

RESULTS

Flow parameters derived from ROI-averaged curves on main pulmonary artery were significantly smaller than that of pixel-wise measurement (P<0.001). Difference between RF_bulk and RF_px for the group of Ecc_sys > 0.3 appears greater compared to the group with Ecc_sys < 0.3. Thirteen out of 68 RF values (19%) were underrated while using bulk analysis.

CONCLUSIONS

The spatial-related flow parameters showed more consistency with the qualitative flow profile pattern for pulmonary arteries, and could be a decisive complement for diagnostic classification.

Comparison of 4D flow and 2D velocity-encoded phase contrast MRI sequences for the evaluation of aortic hemodynamics

The purpose of this study was to compare aortic flow and velocity quantification using 4D flow MRI and 2D CINE phase-contrast (PC)-MRI with either one-directional (2D-1dir) or three-directional (2D-3dir) velocity encoding. 15 healthy volunteers (51 ± 19 years) underwent MRI including (1) breath-holding 2D-1dir and (2) free breathing 2D-3dir PC-MRI in planes orthogonal to the ascending (AA) and descending (DA) aorta, as well as (3) free breathing 4D flow MRI with full thoracic aorta coverage. Flow quantification included the co-registration of the 2D PC acquisition planes with 4D flow MRI data, AA and DA segmentation, and calculation of AA and DA peak systolic velocity, peak flow and net flow volume for all sequences. Additionally, the 2D-3dir velocity taking into account the through-plane component only was used to obtain results analogous to a free breathing 2D-1dir acquisition. Good agreement was found between 4D flow and 2D-3dir peak velocity (differences = -3 to 6 %), peak flow (-7 %) and net volume (-14 to -9 %). In contrast, breath-holding 2D-1dir measurements exhibited indices significantly lower than free breathing 2D-3dir and 2D-1dir (differences = -35 to -7 %, p < 0.05). Finally, high correlations (r ≥ 0.97) were obtained for indices estimated with or without eddy current correction, with the lowest correlation observed for net volume. 4D flow and 2D-3dir aortic hemodynamic indices were in concordance. However, differences between respiration state and 2D-1dir and 2D-3dir measurements indicate that reference values should be established according to the PC-MRI sequence, especially for the widely used net flow (e.g. stroke volume in the AA).

MR phase-contrast imaging in pulmonary hypertension

Pulmonary hypertension (PH) is a life-threatening, multifactorial pathophysiological haemodynamic condition, diagnosed when the mean pulmonary arterial pressure equals or exceeds 25 mmHg at rest during right heart catheterization. Cardiac MRI, in general, and MR phase-contrast (PC) imaging, in particular, have emerged as potential techniques for the standardized assessment of cardiovascular function, morphology and haemodynamics in PH. Allowing the quantification and characterization of macroscopic cardiovascular blood flow, MR PC imaging offers non-invasive evaluation of haemodynamic alterations associated with PH. Techniques used to study the PH include both the routine two-dimensional (2D) approach measuring predominant velocities through an acquisition plane and the rapidly evolving four-dimensional (4D) PC imaging, which enables the assessment of the complete time-resolved, three-directional blood-flow velocity field in a volume. Numerous parameters such as pulmonary arterial mean velocity, vessel distensibility, flow acceleration time and volume and tricuspid regurgitation peak velocity, as well as the duration and onset of vortical blood flow in the main pulmonary artery, have been explored to either diagnose PH or find non-invasive correlates to right heart catheter parameters. Furthermore, PC imaging-based analysis of pulmonary arterial pulse-wave velocities, wall shear stress and kinetic energy losses grants novel insights into cardiopulmonary remodelling in PH. This review aimed to outline the current applications of 2D and 4D PC imaging in PH and show why this technique has the potential to contribute significantly to early diagnosis and characterization of PH.

Phase Error Correction in Time-Averaged 3D Phase Contrast Magnetic Resonance Imaging of the Cerebral Vasculature

Purpose

Volume flow rate (VFR) measurements based on phase contrast (PC)-magnetic resonance (MR) imaging datasets have spatially varying bias due to eddy current induced phase errors. The purpose of this study was to assess the impact of phase errors in time averaged PC-MR imaging of the cerebral vasculature and explore the effects of three common correction schemes (local bias correction (LBC), local polynomial correction (LPC), and whole brain polynomial correction (WBPC)).

Methods

Measurements of the eddy current induced phase error from a static phantom were first obtained. In thirty healthy human subjects, the methods were then assessed in background tissue to determine if local phase offsets could be removed. Finally, the techniques were used to correct VFR measurements in cerebral vessels and compared statistically.

Results

In the phantom, phase error was measured to be <2.1 ml/s per pixel and the bias was reduced with the correction schemes. In background tissue, the bias was significantly reduced, by 65.6% (LBC), 58.4% (LPC) and 47.7% (WBPC) (p < 0.001 across all schemes). Correction did not lead to significantly different VFR measurements in the vessels (p = 0.997). In the vessel measurements, the three correction schemes led to flow measurement differences of -0.04 ± 0.05 ml/s, 0.09 ± 0.16 ml/s, and -0.02 ± 0.06 ml/s. Although there was an improvement in background measurements with correction, there was no statistical difference between the three correction schemes (p = 0.242 in background and p = 0.738 in vessels).

Conclusions

While eddy current induced phase errors can vary between hardware and sequence configurations, our results showed that the impact is small in a typical brain PC-MR protocol and does not have a significant effect on VFR measurements in cerebral vessels.

Effects of contour propagation and background corrections in different MRI flow software packages

Background

Velocity-encoded magnetic resonance imaging (VENC-MRI) is a commonly used technique in cardiac examinations. This technique utilizes the phase shift properties of protons moving along a magnetic field gradient. VENC-MRI offers a unique way of measuring the severity of valve regurgitation by directly quantifying the regurgitation flow volume.

Purpose

To compare flow analysis results of different software programs and to assess the effect of background correction in sample patient cases.

Material and Methods

A phantom was built out of Polymethyl methacrylate (PMMA) which provides tubes of different diameters. These tubes can be connected to an external water circuit to generate a water flow inside the tubes. Expected absolute flow quantities inside the tubes were determined from preset tube- and flow-parameters. Different flow conditions were measured with a VENC-MRI sequence and the images evaluated using different software packages. In a second step six randomly selected patients showing different degrees of aortic insufficiency were evaluated in clinical terms.

Results

The contour propagation algorithms used in the software packages performed differently even on static phantom geometry. In terms of clinical evaluation the software packages performed similarly. Enabling background correction or leaving out manual correction of propagated contours changed results for severity of aortic insufficiency.

Conclusion

Turning on background correction and manual correction of propagated contours in MRI flow volume measurements is strongly recommended.

Cardiovascular magnetic resonance phase contrast imaging

Cardiovascular magnetic resonance (CMR) phase contrast imaging has undergone a wide range of changes with the development and availability of improved calibration procedures, visualization tools, and analysis methods. This article provides a comprehensive review of the current state-of-the-art in CMR phase contrast imaging methodology, clinical applications including summaries of past clinical performance, and emerging research and clinical applications that utilize today's latest technology.

Pulmonary regurgitant volume is superior to fraction using background-corrected phase contrast MRI in determining the severity of regurgitation in repaired tetralogy of Fallot

In the assessment of pulmonary regurgitation (PR) using phase contrast MRI, phase offset errors affect the accuracy of flow. This study evaluated the use of automated background correction for phase offset in the quantification of PR fraction and volume in patients with repaired tetralogy of Fallot (TOF), and to assess its clinical impact. We retrospectively analyzed 203 cardiac MRI studies, performed on 1.5-T scanner. Pulmonary flow (Q_P) and systemic flow (Q_S) was assessed both with and without background correction. Non-corrected and corrected Q_P was correlated with Q_S. PR was correlated with (1) indexed right ventricular end-diastolic volume (RVEDVi) and (2) with differential right and left ventricular stroke volumes (PR_SV). Both PR fraction and volume showed major change after correction (−43 to +36 % and −13 to +13 ml/m²). Corrected Q_P and Q_S were stronger correlated with each other than non-corrected Q_P and Q_S [r = 0.78 vs. 0.73 (p < 0.001)]. Both PR fraction and volume were stronger correlated with RVEDVi, compared to their non-corrected counterparts (p < 0.001). PR volume was stronger correlated with RVEDVi, compared to PR fraction [r = 0.74 vs. 0.69 (p < 0.001)]. When patients were divided according to PR severity, 12 % of patients reclassified after correction. Background correction for phase offset significantly changed the quantification of PR. Non-corrected assessment of PR may result in the misclassification of patients. Our data suggest that the use of PR volume is favourable in the follow-up of patients with repaired TOF.

In vivo analysis of physiological 3D blood flow of cerebral veins

OBJECTIVES

To visualize and quantify physiological blood flow of intracranial veins in vivo using time-resolved, 3D phase-contrast MRI (4D flow MRI), and to test measurement accuracy.

METHODS

Fifteen healthy volunteers underwent repeated ECG-triggered 4D flow MRI (3 Tesla, 32-channel head coil). Intracranial venous blood flow was analysed using dedicated software allowing for blood flow visualization and quantification in analysis planes at the superior sagittal, straight, and transverse sinuses. MRI was evaluated for intra- and inter-observer agreement and scan-rescan reproducibility. Measurements of the transverse sinuses were compared with transcranial two-dimensional duplex ultrasound.

RESULTS

Visualization of 3D blood flow within cerebral sinuses was feasible in 100 % and within at least one deep cerebral vein in 87 % of the volunteers. Blood flow velocity/volume increased along the superior sagittal sinus and was lower in the left compared to the right transverse sinus. Intra- and inter-observer reliability and reproducibility of blood flow velocity (mean difference 0.01/0.02/0.02 m/s) and volume (mean difference 0.0002/-0.0003/0.00003 l/s) were good to excellent. High/low velocities were more pronounced (8 % overestimation/9 % underestimation) in MRI compared to ultrasound.

CONCLUSIONS

Four-dimensional flow MRI reliably visualizes and quantifies three-dimensional cerebral venous blood flow in vivo and is promising for studies in patients with sinus thrombosis and related diseases.

KEY POINTS

• 4D flow MRI can be used to visualize and quantify physiological cerebral venous haemodynamics • Flow quantification within cerebral sinuses reveals high reliability and accuracy of 4D flow MRI • Blood flow volume and velocity increase along the superior sagittal sinus • Limited spatial resolution currently precludes flow quantification in small cerebral veins.

Baseline correction does not improve flow quantification in phase-contrast velocity measurement for routine clinical practice

INTRODUCTION

Velocity offset errors may influence flow measurement in phase-contrast cardiovascular magnetic resonance (CMR). By using a stationary gel phantom, offset errors probably may be corrected. We tested its impact on flow measurement and, in particular, on shunt calculation in patients proven not to have any shunt.

METHODS

Flow measurements were carried out in 24 patients with congenital heart disease. Baseline correction was performed by using a stationary gel phantom.

RESULTS

Significantly more patients without shunts incorrectly showed a calculated shunt after baseline correction.

CONCLUSIONS

Baseline correction did not improve flow measurement and was clinically not relevant for routine CMR.

Analysis of temperature dependence of background phase errors in phase-contrast cardiovascular magnetic resonance

BACKGROUND

The accuracy of phase-contrast cardiovascular magnetic resonance (PC-CMR) can be compromised by background phase errors. It is the objective of the present work to provide an analysis of the temperature dependence of background phase errors in PC-CMR by means of gradient mount temperature sensing and magnetic field monitoring.

METHODS

Background phase errors were measured for various temperatures of the gradient mount using magnetic field monitoring and validated in a static phantom. The effect of thermal changes during k-space acquisition was simulated and confirmed with measurements in a stationary phantom.

RESULTS

The temperature of the gradient mount was found to increase by 20-30 K during PC-CMR measurements of 6-12 min duration. Associated changes in background phase errors of up to 11% or 0.35 radian were measured at 10 cm from the magnet's iso-center as a result of first order offsets. Zeroth order phase errors exhibited little thermal dependence.

CONCLUSIONS

It is concluded that changes in gradient mount temperature significantly modify background phase errors during PC-CMR with high gradient duty cycle. Since temperature increases significantly during the first minutes of scanning the results presented are also of relevance for single-slice or multi-slice PC-CMR scans. The findings prompt for further studies to investigate advanced correction methods taking into account gradient temperature and/or the use of concurrent field-monitoring to map gradient-induced fields throughout the scan.

Quantifying errors in flow measurement using phase contrast magnetic resonance imaging: comparison of several boundary detection methods

Quantifying flow from phase-contrast MRI (PC-MRI) data requires that the vessels of interest be segmented. The estimate of the vessel area will dictate the type and magnitude of the error sources that affect the flow measurement. These sources of errors are well understood, and mathematical expressions have been derived for them in previous work. However, these expressions contain many parameters that render them difficult to use for making practical error estimates. In this work, some realistic assumptions were made that allow for the simplification of such expressions in order to make them more useful. These simplified expressions were then used to numerically simulate the effect of segmentation accuracy and provide some criteria that if met, would keep errors in flow quantification below 10% or 5%. Four different segmentation methods were used on simulated and phantom MRA data to verify the theoretical results. Numerical simulations showed that including partial volumed edge pixels in vessel segmentation provides less error than missing them. This was verified with MRA simulations, as the best performing segmentation method generally included such pixels. Further, it was found that to obtain a flow error of less than 10% (5%), the vessel should be at least 4 (5) pixels in diameter, have an SNR of at least 10:1 and have a peak velocity to saturation cut-off velocity ratio of at least 5:3.

Method for calculating confidence intervals for phase contrast flow measurements

BACKGROUND

Phase contrast (PC) measurements play an important role in several cardiovascular magnetic resonance (CMR) protocols but considerable variation is observed in such measurements. Part of this variation stems from the propagation of thermal noise from the measurement data through the image reconstruction to the region of interest analysis used in flow measurement, which limits the precision. The purpose of this study was to develop a method for direct estimation of the variation caused by thermal noise and to validate this method in phantom and in vivo data.

METHODS

The estimation of confidence intervals in flow measurements is complicated by noise correlation among the image pixels and cardiac phases. This correlation is caused by sequence and reconstruction parameters. A method for the calculation of the standard deviation of region of interest measurements was adapted and expanded to accommodate typical clinical PC measurements and the region-of-interest analysis used for such measurements. This included the dependency between cardiac phases that arises due to retrospective cardiac gating used in such studies. The proposed method enables calculation of standard deviations of flow measurements without the need for repeated experiments or repeated reconstructions. The method was compared to repeated trials in phantom measurements and pseudo replica reconstructions of in vivo data. Three different flow protocols (free breathing and breath hold with various accelerations) were compared in terms of the confidence interval ranges caused by thermal noise in the measurement data.

RESULTS

Using the proposed method it was possible to accurately predict confidence intervals for flow measurements. The method was in good agreement with repeated measurements in phantom experiments and there was also good agreement with confidence intervals predicted by pseudo replica reconstructions in both phantom and in vivo data. The proposed method was used to demonstrate that the variation in cardiac output caused by thermal noise is on the order of 1% in clinically used free breathing protocols, and on the order of 3-5% in breath-hold protocols with higher parallel imaging factors.

CONCLUSIONS

It is possible to calculate confidence intervals for Cartesian PC contrast flow measurements directly without the need for time-consuming pseudo replica reconstructions.

Towards highly accelerated Cartesian time-resolved 3D flow cardiovascular magnetic resonance in the clinical setting

BACKGROUND

The clinical applicability of time-resolved 3D flow cardiovascular magnetic resonance (CMR) remains compromised by the long scan times associated with phase-contrast imaging. The present work demonstrates the applicability of 8-fold acceleration of Cartesian time-resolved 3D flow CMR in 10 volunteers and in 9 patients with different congenital heart diseases (CHD). It is demonstrated that accelerated 3D flow CMR data acquisition and image reconstruction using k-t PCA (principal component analysis) can be implemented into clinical workflow and results are sufficiently accurate relative to conventional 2D flow CMR to permit for comprehensive flow quantification in CHD patients.

METHODS

The fidelity of k-t PCA was first investigated on retrospectively undersampled data for different acceleration factors and compared to k-t SENSE and fully sampled reference data. Subsequently, k-t PCA with 8-fold nominal undersampling was applied on 10 healthy volunteers and 9 CHD patients on a clinical 1.5 T MR scanner. Quantitative flow validation was performed in vessels of interest on the 3D flow datasets and compared to 2D through-plane flow acquisitions. Particle trace analysis was used to qualitatively visualise flow patterns in patients.

RESULTS

Accelerated time-resolved 3D flow data were successfully acquired in all subjects with 8-fold nominal scan acceleration. Nominal scan times excluding navigator efficiency were on the order of 6 min and 7 min in patients and volunteers. Mean differences in stroke volume in selected vessels of interest were 2.5 ± 8.4 ml and 1.63 ± 4.8 ml in volunteers and patients, respectively. Qualitative flow pattern analysis in the time-resolved 3D dataset revealed valuable insights into hemodynamics including circular and helical patterns as well as flow distributions and origin in the Fontan circulation.

CONCLUSION

Highly accelerated time-resolved 3D flow using k-t PCA is readily applicable in clinical routine protocols of CHD patients. Nominal scan times of 6 min are well tolerated and allow for quantitative and qualitative flow assessment in all great vessels.

Volumetric velocity measurements in restricted geometries using spiral sampling: a phantom study

OBJECT

The aim of this study was to evaluate the accuracy of maximum velocity measurements using volumetric phase-contrast imaging with spiral readouts in a stenotic flow phantom.

MATERIALS AND METHODS

In a phantom model, maximum velocity, flow, pressure gradient, and streamline visualizations were evaluated using volumetric phase-contrast magnetic resonance imaging (MRI) with velocity encoding in one (extending on current clinical practice) and three directions (for characterization of the flow field) using spiral readouts. Results of maximum velocity and pressure drop were compared to computational fluid dynamics (CFD) simulations, as well as corresponding low-echo-time (TE) Cartesian data. Flow was compared to 2D through-plane phase contrast (PC) upstream from the restriction.

RESULTS

Results obtained with 3D through-plane PC as well as 4D PC at shortest TE using a spiral readout showed excellent agreements with the maximum velocity values obtained with CFD (<1 % for both methods), while larger deviations were seen using Cartesian readouts (-2.3 and 13 %, respectively). Peak pressure drop calculations from 3D through-plane PC and 4D PC spiral sequences were respectively 14 and 13 % overestimated compared to CFD.

CONCLUSION

Identification of the maximum velocity location, as well as the accurate velocity quantification can be obtained in stenotic regions using short-TE spiral volumetric PC imaging.

Analysis of an automated background correction method for cardiovascular MR phase contrast imaging in children and young adults

BACKGROUND

Phase contrast magnetic resonance imaging (MRI) is a powerful tool for evaluating vessel blood flow. Inherent errors in acquisition, such as phase offset, eddy currents and gradient field effects, can cause significant inaccuracies in flow parameters. These errors can be rectified with the use of background correction software.

OBJECTIVE

To evaluate the performance of an automated phase contrast MRI background phase correction method in children and young adults undergoing cardiac MR imaging.

MATERIALS AND METHODS

We conducted a retrospective review of patients undergoing routine clinical cardiac MRI including phase contrast MRI for flow quantification in the aorta (Ao) and main pulmonary artery (MPA). When phase contrast MRI of the right and left pulmonary arteries was also performed, these data were included. We excluded patients with known shunts and metallic implants causing visible MRI artifact and those with more than mild to moderate aortic or pulmonary stenosis. Phase contrast MRI of the Ao, mid MPA, proximal right pulmonary artery (RPA) and left pulmonary artery (LPA) using 2-D gradient echo Fast Low Angle SHot (FLASH) imaging was acquired during normal respiration with retrospective cardiac gating. Standard phase image reconstruction and the automatic spatially dependent background-phase-corrected reconstruction were performed on each phase contrast MRI dataset. Non-background-corrected and background-phase-corrected net flow, forward flow, regurgitant volume, regurgitant fraction, and vessel cardiac output were recorded for each vessel. We compared standard non-background-corrected and background-phase-corrected mean flow values for the Ao and MPA. The ratio of pulmonary to systemic blood flow (Qp:Qs) was calculated for the standard non-background and background-phase-corrected data and these values were compared to each other and for proximity to 1. In a subset of patients who also underwent phase contrast MRI of the MPA, RPA, and LPA a comparison was made between standard non-background-corrected and background-phase-corrected mean combined flow in the branch pulmonary arteries and MPA flow. All comparisons were performed using the Wilcoxon sign rank test (α = 0.05).

RESULTS

Eighty-five children and young adults (mean age 14 years; range 10 days to 32 years) met the criteria for inclusion. Background-phase-corrected mean flow values for the Ao and MPA were significantly lower than those for non-background-corrected standard Ao (P = 0.0004) and MPA flow values (P < 0.0001), respectively. However, no significant difference was seen between the standard non-background (P = 0.295) or background-phase-corrected (P = 0.0653) mean Ao and MPA flow values. Neither the mean standard non-background-corrected (P = 0.408) nor the background-phase-corrected (P = 0.0684) Qp:Qs was significantly different from 1. However in the 27 patients with standard non-background-corrected data, the difference between the Ao and MPA flow values was greater than 10%. There were 19 patients with background-phase-corrected data in which the difference between the Ao and MPA flow values was greater than 10%. In the subset of 43 patients who underwent MPA and branch pulmonary artery phase contrast MRI, the sum of the standard non-background-corrected mean RPA and LPA flow values was significantly different from the standard non-background-corrected mean MPA flow (P = 0.0337). The sum of the background-phase-corrected mean RPA and LPA flow values was not significantly different from the background-phase-corrected mean MPA flow value (P = 0.1328), suggesting improvement in pulmonary artery flow calculations using background-phase-correction.

CONCLUSION

Our data suggest that background phase correction of phase contrast MRI data does not significantly change Qp:Qs quantification, and there are residual errors in expected Qp:Qs quantification despite background phase correction. However the use of background phase correction does improve quantification of MPA flow relative to combined RPA and LPA flow. Further work is needed to validate these findings in other patient populations, using other MRI units, and across vendors.