Unwrapping phase contrast MRI by iterative graph cuts

PURPOSE

To develop and evaluate a phase unwrapping method for cine phase contrast MRI based on graph cuts.

METHODS

A proposed Iterative Graph Cuts method was evaluated in 10 cardiac patients with two-dimensional flow quantification which was repeated at low venc settings to provoke wrapping. The images were also unwrapped by a path-following method (ROMEO), and a Laplacian-based method (LP). Net flow was quantified using semi-automatic vessel segmentation. High venc images were also wrapped retrospectively to asses the residual amount of wrapped voxels.

RESULTS

The absolute net flow error after unwrapping at venc = 100 cm/s was 1.8 mL, which was 0.83 mL smaller than for LP. The repeatability error at high venc without unwrapping was 2.5 mL. The error at venc = 50 cm/s was 7.5 mL, which was 8.2 mL smaller than for ROMEO and 5.7 mL smaller than for LP. For retrospectively wrapped images with synthetic venc of 100/50/25 cm/s, the residual amount of wrapped voxels was 0.00/0.12/0.79%, which was 0.09/0.26/8.0 percentage points smaller than for LP. With synthetic venc of 25 cm/s, omitting magnitude information resulted in 3.2 percentage points more wrapped voxels, and only spatial/temporal unwrapping resulted in 4.6/21 percentage points more wrapped voxels compared to spatiotemporal unwrapping.

CONCLUSION

Iterative Graph Cuts enables unwrapping of cine phase contrast MRI with very small errors, except for at extreme blood velocities, with equal or better performance compared to ROMEO and LP. The use of magnitude information and spatiotemporal unwrapping is recommended.

Unraveling diurnal and technical variability in cerebral hemodynamics from neurovascular 4D-Flow MRI

Neurovascular 4D-Flow MRI enables non-invasive evaluation of cerebral hemodynamics including measures of cerebral blood flow (CBF), vessel pulsatility index (PI), and cerebral pulse wave velocity (PWV). 4D-Flow measures have been linked to various neurovascular disorders including small vessel disease and Alzheimer's disease; however, physiological and technical sources of variability are not well established. Here, we characterized sources of diurnal physiological and technical variability in cerebral hemodynamics using 4D-Flow in a retrospective study of cognitively unimpaired older adults (N = 750) and a prospective study of younger adults (N = 10). Younger participants underwent repeated MRI sessions at 7am, 4 pm, and 10 pm. In the older cohort, having an MRI earlier on the day was significantly associated with higher CBF and lower PI. In prospective experiments, time of day significantly explained variability in CBF and PI; however, not in PWV. Test-retest experiments showed high CBF intra-session repeatability (repeatability coefficient (RPC) =7.2%), compared to lower diurnal repeatability (RPC = 40%). PI and PWV displayed similar intra-session and diurnal variability (PI intra-session RPC = 22%, RPC = 24% 7am vs 4 pm; PWV intra-session RPC = 17%, RPC = 21% 7am vs 4 pm). Overall, CBF measures showed low technical variability, supporting diurnal variability is from physiology. PI and PWV showed higher technical variability but less diurnal variability.

Four-dimensional flow MRI for quantitative assessment of cerebrospinal fluid dynamics: Status and opportunities

Neurological disorders can manifest with altered neurofluid dynamics in different compartments of the central nervous system. These include alterations in cerebral blood flow, cerebrospinal fluid (CSF) flow, and tissue biomechanics. Noninvasive quantitative assessment of neurofluid flow and tissue motion is feasible with phase contrast magnetic resonance imaging (PC MRI). While two-dimensional (2D) PC MRI is routinely utilized in research and clinical settings to assess flow dynamics through a single imaging slice, comprehensive neurofluid dynamic assessment can be limited or impractical. Recently, four-dimensional (4D) flow MRI (or time-resolved three-dimensional PC with three-directional velocity encoding) has emerged as a powerful extension of 2D PC, allowing for large volumetric coverage of fluid velocities at high spatiotemporal resolution within clinically reasonable scan times. Yet, most 4D flow studies have focused on blood flow imaging. Characterizing CSF flow dynamics with 4D flow (i.e., 4D CSF flow) is of high interest to understand normal brain and spine physiology, but also to study neurological disorders such as dysfunctional brain metabolite waste clearance, where CSF dynamics appear to play an important role. However, 4D CSF flow imaging is challenged by the long T1 time of CSF and slower velocities compared with blood flow, which can result in longer scan times from low flip angles and extended motion-sensitive gradients, hindering clinical adoption. In this work, we review the state of 4D CSF flow MRI including challenges, novel solutions from current research and ongoing needs, examples of clinical and research applications, and discuss an outlook on the future of 4D CSF flow.

Super-resolution Left Ventricular Flow and Pressure Mapping by Navier-Stokes-Informed Neural Networks

Intraventricular vector flow mapping (VFM) is a growingly adopted echocardiographic modality that derives time-resolved two-dimensional flow maps in the left ventricle (LV) from color-Doppler sequences. Current VFM models rely on kinematic constraints arising from planar flow incompressibility. However, these models are not informed by crucial information about flow physics; most notably the pressure and shear forces within the fluid and the resulting accelerations. This limitation has rendered VFM unable to combine information from different time frames in an acquisition sequence or derive fluctuating pressure maps. In this study, we leveraged recent advances in artificial intelligence (AI) to develop AI-VFM, a vector flow mapping modality that uses physics-informed neural networks (PINNs) encoding mass conservation and momentum balance inside the LV, and no-slip boundary conditions at the LV endocardium. AI-VFM recovers the flow and pressure fields in the LV from standard echocardiographic scans. It performs phase unwrapping and recovers flow data in areas without input color-Doppler data. AI-VFM also recovers complete flow maps at time points without color-Doppler input data, producing super-resolution flow maps. We show that informing the PINNs with momentum balance is essential to achieving temporal super-resolution and significantly increases the accuracy of AI-VFM compared to informing the PINNs only with mass conservation. AI-VFM is solely informed by each patient's flow physics; it does not utilize explicit smoothness constraints or incorporate data from other patients or flow models. AI-VFM takes 15 minutes to run in off-the-shelf graphics processing units and its underlying PINN framework could be extended to map other flow-associated metrics like blood residence time or the concentration of coagulation species.

A comparison of phase unwrapping methods in velocity-encoded MRI for aortic flows

PURPOSE

The phase of a MRI signal is used to encode the velocity of blood flow. Phase unwrapping artifacts may appear when aiming to improve the velocity-to-noise ratio (VNR) of the measured velocity field. This study aims to compare various unwrapping algorithms on ground-truth synthetic data generated using computational fluid dynamics (CFD) simulations.

METHODS

We compare four different phase unwrapping algorithms on two different synthetic datasets of four-dimensional flow MRI and 26 datasets of 2D PC-MRI acquisitions including the ascending and descending aorta. The synthetic datasets are constructed using CFD simulations of an aorta with a coarctation, with different levels of spatiotemporal resolutions and noise. The error of the unwrapped images was assessed by comparison against the ground truth velocity field in the synthetic data and dual-VENC reconstructions in the in vivo data.

RESULTS

Using the unwrapping algorithms, we were able to remove aliased voxels in the data almost entirely, reducing the L2-error compared to the ground truth by 50%-80%. Results indicated that the best choice of algorithm depend on the spatiotemporal resolution and noise level of the dataset. Temporal unwrapping is most successful with a high temporal and low spatial resolution (δ t = 30 $$ \delta t=30 $$  ms,h = 2 . 5 $$ h=2.5 $$  mm), reducing the L2-error by 70%-85%, while Laplacian unwrapping performs better with a lower temporal or better spatial resolution (δ t = 60 $$ \delta t=60 $$ ms,h = 1 . 5 $$ h=1.5 $$ mm), especially for signal-to-noise ratio (SNR) 12 as opposed to SNR 15, with an error reduction of 55%-85% compared to the 50%-75% achieved by the Temporal method. The differences in performance between the methods are statistically significant.

CONCLUSIONS

The temporal method and spatiotemporal Laplacian method provide the best results, with the spatiotemporal Laplacian being more robust. However, single-V enc $$ {V}_{\mathrm{enc}} $$ methods only situationally and not generally reach the performance of dual-V enc $$ {V}_{\mathrm{enc}} $$ unwrapping methods.

Haemodynamic effects of hypertension and type 2 diabetes: Insights from a 4D flow MRI-based personalized cardiovascular mathematical model

Type 2 diabetes (T2D) and hypertension increase the risk of cardiovascular diseases mediated by whole-body changes to metabolism, cardiovascular structure and haemodynamics. The haemodynamic changes related to hypertension and T2D are complex and subject-specific, however, and not fully understood. We aimed to investigate the haemodynamic mechanisms in T2D and hypertension by comparing the haemodynamics between healthy controls and subjects with T2D, hypertension, or both. For all subjects, we combined 4D flow magnetic resonance imaging data, brachial blood pressure and a cardiovascular mathematical model to create a comprehensive subject-specific analysis of central haemodynamics. When comparing the subject-specific haemodynamic parameters between the four groups, the predominant haemodynamic difference is impaired left ventricular relaxation in subjects with both T2D and hypertension compared to subjects with only T2D, only hypertension and controls. The impaired relaxation indicates that, in this cohort, the long-term changes in haemodynamic load of co-existing T2D and hypertension cause diastolic dysfunction demonstrable at rest, whereas either disease on its own does not. However, through subject-specific predictions of impaired relaxation, we show that altered relaxation alone is not enough to explain the subject-specific and group-related differences; instead, a combination of parameters is affected in T2D and hypertension. These results confirm previous studies that reported more adverse effects from the combination of T2D and hypertension compared to either disease on its own. Furthermore, this shows the potential of personalized cardiovascular models in providing haemodynamic mechanistic insights and subject-specific predictions that could aid in the understanding and treatment planning of patients with T2D and hypertension. KEY POINTS: The combination of 4D flow magnetic resonance imaging data and a cardiovascular mathematical model allows for a comprehensive analysis of subject-specific haemodynamic parameters that otherwise cannot be derived non-invasively. Using this combination, we show that diastolic dysfunction in subjects with both type 2 diabetes (T2D) and hypertension is the main group-level difference between controls, subjects with T2D, subjects with hypertension, and subjects with both T2D and hypertension. These results suggest that, in this relatively healthy population, the additional load of both hypertension and T2D affects the haemodynamic function of the left ventricle, whereas each disease on its own is not enough to cause significant effects under resting conditions. Finally, using the subject-specific model, we show that the haemodynamic effects of diastolic dysfunction alone are not sufficient to explain all the observed haemodynamic differences. Instead, additional subject-specific variations in cardiac and vascular function combine to explain the complex haemodynamics of subjects affected by hypertension and/or T2D.

Hemodynamic forces from 4D flow magnetic resonance imaging predict left ventricular remodeling following cardiac resynchronization therapy

BACKGROUND

Patients with heart failure and left bundle branch block (LBBB) may receive cardiac resynchronization therapy (CRT), but current selection criteria are imprecise, and many patients have limited treatment response. Hemodynamic forces (HDF) have been suggested as a marker for CRT response. The aim of this study was therefore to investigate left ventricular (LV) HDF as a predictive marker for LV remodeling after CRT.

METHODS

Patients with heart failure, EF < 35% and LBBB (n = 22) underwent CMR with 4D flow prior to CRT. LV HDF were computed in three directions using the Navier-Stokes equations, reported in median N [interquartile range], and the ratio of transverse/longitudinal HDF was calculated for systole and diastole. Transthoracic echocardiography was performed before and 6 months after CRT. Patients with end-systolic volume reduction ≥ 15% were defined as responders.

RESULTS

Non-responders had smaller HDF than responders in the inferior-anterior direction in systole (0.06 [0.03] vs. 0.07 [0.03], p = 0.04), and in the apex-base direction in diastole (0.09 [0.02] vs. 0.1 [0.05], p = 0.047). Non-responders had larger diastolic HDF ratio compared to responders (0.89 vs. 0.67, p = 0.004). ROC analysis of diastolic HDF ratio for identifying CRT non-responders had AUC of 0.88 (p = 0.005) with sensitivity 57% and specificity 100% for ratio > 0.87. Intragroup comparison found higher HDF ratio in systole compared to diastole for responders (p = 0.003), but not for non-responders (p = 0.8).

CONCLUSION

Hemodynamic force ratio is a potential marker for identifying patients with heart failure and LBBB who are unlikely to benefit from CRT. Larger-scale studies are required before implementation of HDF analysis into clinical practice.

Four-Dimensional Flow MR Imaging: Technique and Advances

4D Flow MRI is an advanced imaging technique for comprehensive non-invasive assessment of the cardiovascular system. The capture of the blood velocity vector field throughout the cardiac cycle enables measures of flow, pulse wave velocity, kinetic energy, wall shear stress, and more. Advances in hardware, MRI data acquisition and reconstruction methodology allow for clinically feasible scan times. The availability of 4D Flow analysis packages allows for more widespread use in research and the clinic and will facilitate much needed multi-center, multi-vendor studies in order to establish consistency across scanner platforms and to enable larger scale studies to demonstrate clinical value.

4D Flow cardiovascular magnetic resonance consensus statement: 2023 update

Hemodynamic assessment is an integral part of the diagnosis and management of cardiovascular disease. Four-dimensional cardiovascular magnetic resonance flow imaging (4D Flow CMR) allows comprehensive and accurate assessment of flow in a single acquisition. This consensus paper is an update from the 2015 '4D Flow CMR Consensus Statement'. We elaborate on 4D Flow CMR sequence options and imaging considerations. The document aims to assist centers starting out with 4D Flow CMR of the heart and great vessels with advice on acquisition parameters, post-processing workflows and integration into clinical practice. Furthermore, we define minimum quality assurance and validation standards for clinical centers. We also address the challenges faced in quality assurance and validation in the research setting. We also include a checklist for recommended publication standards, specifically for 4D Flow CMR. Finally, we discuss the current limitations and the future of 4D Flow CMR. This updated consensus paper will further facilitate widespread adoption of 4D Flow CMR in the clinical workflow across the globe and aid consistently high-quality publication standards.

A clinician's guide to understanding aortic 4D flow MRI

Four-dimensional flow magnetic resonance imaging is an emerging technique which may play a role in diagnosis and risk-stratification of aortic disease. Some knowledge of flow dynamics and related parameters is necessary to understand and apply this technique in clinical workflows. The purpose of the current review is to provide a guide for clinicians to the basics of flow imaging, frequently used flow-related parameters, and their relevance in the context of aortic disease.Clinical relevance statement Understanding normal and abnormal aortic flow could improve clinical care in patients with aortic disease.

Fetal cardiovascular blood flow MRI: techniques and applications

Fetal cardiac MRI is challenging due to fetal and maternal movements as well as the need for a reliable cardiac gating signal and high spatiotemporal resolution. Ongoing research and recent technical developments to address these challenges show the potential of MRI as an adjunct to ultrasound for the assessment of the fetal heart and great vessels. MRI measurements of blood flow have enabled the assessment of normal fetal circulation as well as conditions with disrupted circulations, such as congenital heart disease, along with associated organ underdevelopment and hemodynamic instability. This review provides details of the techniques used in fetal cardiovascular blood flow MRI, including single slice and volumetric imaging sequences, post-processing and analysis, along with a summary of applications in human studies and animal models.

Cerebral hemodynamics comparison using transcranial doppler ultrasound and 4D flow MRI

Introduction: Age-related changes in cerebral hemodynamics are controversial and discrepancies may be due to experimental techniques. As such, the purpose of this study was to compare cerebral hemodynamics measurements of the middle cerebral artery (MCA) between transcranial Doppler ultrasound (TCD) and four-dimensional flow MRI (4D flow MRI). Methods: Twenty young (25 ± 3 years) and 19 older (62 ± 6 years) participants underwent two randomized study visits to evaluate hemodynamics at baseline (normocapnia) and in response to stepped hypercapnia (4% CO₂, and 6% CO₂) using TCD and 4D flow MRI. Cerebral hemodynamic measures included MCA velocity, MCA flow, cerebral pulsatility index (PI) and cerebrovascular reactivity to hypercapnia. MCA flow was only assessed using 4D flow MRI. Results: MCA velocity between the TCD and 4D flow MRI methods was positively correlated across the normocapnia and hypercapnia conditions (r = 0.262; p = 0.004). Additionally, cerebral PI was significantly correlated between TCD and 4D flow MRI across the conditions (r = 0.236; p = 0.010). However, there was no significant association between MCA velocity using TCD and MCA flow using 4D flow MRI across the conditions (r = 0.079; p = 0.397). When age-associated differences in cerebrovascular reactivity using conductance were compared using both methodologies, cerebrovascular reactivity was greater in young adults compared to older adults when using 4D flow MRI (2.11 ± 1.68 mL/min/mmHg/mmHg vs. 0.78 ± 1.68 mL/min/mmHg/mmHg; p = 0.019), but not with TCD (0.88 ± 1.01 cm/s/mmHg/mmHg vs. 0.68 ± 0.94 cm/s/mmHg/mmHg; p = 0.513). Conclusion: Our results demonstrated good agreement between the methods at measuring MCA velocity during normocapnia and in response to hypercapnia, but MCA velocity and MCA flow were not related. In addition, measurements using 4D flow MRI revealed effects of aging on cerebral hemodynamics that were not apparent using TCD.

Normative Cerebral Hemodynamics in Middle-aged and Older Adults Using 4D Flow MRI: Initial Analysis of Vascular Aging

Background Characterizing cerebrovascular hemodynamics in older adults is important for identifying disease and understanding normal neurovascular aging. Four-dimensional (4D) flow MRI allows for a comprehensive assessment of cerebral hemodynamics in a single acquisition. Purpose To establish reference intracranial blood flow and pulsatility index values in a large cross-sectional sample of middle-aged (45-65 years) and older (>65 years) adults and characterize the effect of age and sex on blood flow and pulsatility. Materials and Methods In this retrospective study, patients aged 45-93 years (cognitively unimpaired) underwent cranial 4D flow MRI between March 2010 and March 2020. Blood flow rates and pulsatility indexes from 13 major arteries and four venous sinuses and total cerebral blood flow were collected. Intraobserver and interobserver reproducibility of flow and pulsatility measures was assessed in 30 patients. Descriptive statistics (mean ± SD) of blood flow and pulsatility were tabulated for the entire group and by age and sex. Multiple linear regression and linear mixed-effects models were used to assess the effect of age and sex on total cerebral blood flow and vessel-specific flow and pulsatility, respectively. Results There were 759 patients (mean age, 65 years ± 8 [SD]; 506 female patients) analyzed. For intra- and interobserver reproducibility, median intraclass correlation coefficients were greater than 0.90 for flow and pulsatility measures across all vessels. Regression coefficients β ± standard error from multiple linear regression showed a 4 mL/min decrease in total cerebral blood flow each year (age β = -3.94 mL/min per year ± 0.44; P < .001). Mixed effects showed a 1 mL/min average annual decrease in blood flow (age β = -0.95 mL/min per year ± 0.16; P < .001) and 0.01 arbitrary unit (au) average annual increase in pulsatility over all vessels (age β = 0.011 au per year ± 0.001; P < .001). No evidence of sex differences was observed for flow (β = -1.60 mL/min per male patient ± 1.77; P = .37), but pulsatility was higher in female patients (sex β = -0.018 au per male patient ± 0.008; P = .02). Conclusion Normal reference values for blood flow and pulsatility obtained using four-dimensional flow MRI showed correlations with age. © RSNA, 2023 Supplemental material is available for this article. See also the editorial by Steinman in this issue.

Automated hemodynamic assessment for cranial 4D flow MRI

Cranial 4D flow MRI post-processing typically involves manual user interaction which is time-consuming and associated with poor repeatability. The primary goal of this study is to develop a robust quantitative velocity tool (QVT) that utilizes threshold-based segmentation techniques to improve segmentation quality over prior approaches based on centerline processing schemes (CPS) that utilize k-means clustering segmentation. This tool also includes an interactive 3D display designed for simplified vessel selection and automated hemodynamic visualization and quantification. The performances of QVT and CPS were compared in vitro in a flow phantom and in vivo in 10 healthy participants. Vessel segmentations were compared with ground-truth computed tomography in vitro (29 locations) and manual segmentation in vivo (13 locations) using linear regression. Additionally, QVT and CPS MRI flow rates were compared to perivascular ultrasound flow in vitro using linear regression. To assess internal consistency of flow measures in vivo, conservation of flow was assessed at vessel junctions using linear regression and consistency of flow along vessel segments was analyzed by fitting a Gaussian distribution to a histogram of normalized flow values. Post-processing times were compared between the QVT and CPS using paired t-tests. Vessel areas segmented in vitro (CPS: slope = 0.71, r = 0.95 and QVT: slope = 1.03, r = 0.95) and in vivo (CPS: slope = 0.61, r = 0.96 and QVT: slope = 0.93, r = 0.96) were strongly correlated with ground-truth area measurements. However, CPS (using k-means segmentation) consistently underestimated vessel areas. Strong correlations were observed between QVT and ultrasound flow (slope = 0.98, r = 0.96) as well as flow at junctions (slope = 1.05, r = 0.98). Mean and standard deviation of flow along vessel segments was 9.33e-16 ± 3.05%. Additionally, the QVT demonstrated excellent interobserver agreement and significantly reduced post-processing by nearly 10 min (p < 0.001). By completely automating post-processing and providing an easy-to-use 3D visualization interface for interactive vessel selection and hemodynamic quantification, the QVT offers an efficient, robust, and repeatable means to analyze cranial 4D flow MRI. This software is freely available at: https://github.com/uwmri/QVT.

Evaluation of aortic stenosis: From Bernoulli and Doppler to Navier-Stokes

Uni-dimensional Doppler echocardiography data provide the mainstay of quantative assessment of aortic stenosis, with the transvalvular pressure drop a key indicator of haemodynamic burden. Sophisticated methods of obtaining velocity data, combined with improved computational analysis, are facilitating increasingly robust and reproducible measurement. Imaging modalities which permit acquisition of three-dimensional blood velocity vector fields enable angle-independent valve interrogation and calculation of enhanced measures of the transvalvular pressure drop. This manuscript clarifies the fundamental principles of physics that underpin the evaluation of aortic stenosis and explores modern techniques that may provide more accurate means to grade aortic stenosis and inform appropriate management.

Advances in machine learning applications for cardiovascular 4D flow MRI

Four-dimensional flow magnetic resonance imaging (MRI) has evolved as a non-invasive imaging technique to visualize and quantify blood flow in the heart and vessels. Hemodynamic parameters derived from 4D flow MRI, such as net flow and peak velocities, but also kinetic energy, turbulent kinetic energy, viscous energy loss, and wall shear stress have shown to be of diagnostic relevance for cardiovascular diseases. 4D flow MRI, however, has several limitations. Its long acquisition times and its limited spatio-temporal resolutions lead to inaccuracies in velocity measurements in small and low-flow vessels and near the vessel wall. Additionally, 4D flow MRI requires long post-processing times, since inaccuracies due to the measurement process need to be corrected for and parameter quantification requires 2D and 3D contour drawing. Several machine learning (ML) techniques have been proposed to overcome these limitations. Existing scan acceleration methods have been extended using ML for image reconstruction and ML based super-resolution methods have been used to assimilate high-resolution computational fluid dynamic simulations and 4D flow MRI, which leads to more realistic velocity results. ML efforts have also focused on the automation of other post-processing steps, by learning phase corrections and anti-aliasing. To automate contour drawing and 3D segmentation, networks such as the U-Net have been widely applied. This review summarizes the latest ML advances in 4D flow MRI with a focus on technical aspects and applications. It is divided into the current status of fast and accurate 4D flow MRI data generation, ML based post-processing tools for phase correction and vessel delineation and the statistical evaluation of blood flow.

Motion-corrected 4D-Flow MRI for neurovascular applications

Neurovascular 4D-Flow MRI has emerged as a powerful tool for comprehensive cerebrovascular hemodynamic characterization. Clinical studies in at risk populations such as aging adults indicate hemodynamic markers can be confounded by motion-induced bias. This study develops and characterizes a high fidelity 3D self-navigation approach for retrospective rigid motion correction of neurovascular 4D-Flow data. A 3D radial trajectory with pseudorandom ordering was combined with a multi-resolution low rank regularization approach to enable high spatiotemporal resolution self-navigators from extremely undersampled data. Phantom and volunteer experiments were performed at 3.0T to evaluate the ability to correct for different amounts of induced motions. In addition, the approach was applied to clinical-research exams from ongoing aging studies to characterize performance in the clinical setting. Simulations, phantom and volunteer experiments with motion correction produced images with increased vessel conspicuity, reduced image blurring, and decreased variability in quantitative measures. Clinical exams revealed significant changes in hemodynamic parameters including blood flow rates, flow pulsatility index, and lumen areas after motion correction in probed cerebral arteries (Flow: P<0.001 Lt ICA, P=0.002 Rt ICA, P=0.004 Lt MCA, P=0.004 Rt MCA; Area: P<0.001 Lt ICA, P<0.001 Rt ICA, P=0.004 Lt MCA, P=0.004 Rt MCA; flow pulsatility index: P=0.042 Rt ICA, P=0.002 Lt MCA). Motion induced bias can lead to significant overestimation of hemodynamic markers in cerebral arteries. The proposed method reduces measurement bias from rigid motion in neurovascular 4D-Flow MRI in challenging populations such as aging adults.

Simple phase unwrapping method with continuous convex minimization

Phase unwrapping is a problem to reconstruct true phase values from modulo 2π phase values measured using various phase imaging techniques. This procedure is essentially formulated as a discrete optimization problem. However, most energy minimization methods using continuous optimization techniques have ignored the discrete nature and solved it as a continuous minimization problem directly, leading to losing exactness of the algorithms. We propose a new minimum norm method that can yield the optimal solution of the discrete problem by minimizing a continuous energy function. In contrast to the graph-cuts method, which is state of the art in this field, the proposed method requires much less memory space and a very simple implementation. Therefore, it can be simply extended to 3D or 4D phase unwrapping problems.

Kinetic energy of left ventricular blood flow across heart failure phenotypes and in subclinical diastolic dysfunction

BACKGROUND

Kinetic energy (KE) of intracardiac blood flow reflects myocardial work spent on accelerating blood and provides a mechanistic window into diastolic filling dynamics. Diastolic dysfunction may represent an early stage in the development of heart failure (HF). Here we evaluated the hemodynamic effects of impaired diastolic function in subjects with and without HF, testing the hypothesis that left ventricular KE differs between controls, subjects with subclinical diastolic dysfunction (SDD), and HF patients.

METHODS

We studied 77 subjects (16 controls, 20 subjects with SDD, 16 HFpEF, 9 HFmrEF, and 16 HFrEF patients, age- and sex-matched at the group level). Cardiac magnetic resonance at 1.5T included intracardiac 4D flow and cine imaging. Left ventricular KE was calculated as 0.5*m*v².

RESULTS

Systolic KE was similar between groups (p>0.4), also after indexing to stroke volume (p=0.25), and was primarily driven by ventricular emptying rate (p<0.0001, R²=0.52). Diastolic KE was higher in heart failure patients than controls (p<0.05) but similar between SDD and HFpEF (p>0.18), correlating with inflow conditions (E-wave velocity, p<0.0001, R²=0.24) and end-diastolic volume (p=0.0003, R²=0.17) but not with average e' (p=0.07).

CONCLUSIONS

Diastolic KE differs between controls and heart failure, suggesting more work is spent filling the failing ventricle, while systolic KE does not differentiate between well-matched groups with normal ejection fraction even in the presence of relaxation abnormalities and heart failure. Mechanistically, KE reflects the acceleration imparted on the blood and is driven by variations in ventricular emptying and filling rates, volumes, and heart rate, regardless of underlying pathology.

Accelerated dual-venc 4D flow MRI with variable high-venc spatial resolution for neurovascular applications

Purpose

Dual‐velocity encoded (dual‐venc or DV) 4D flow MRI achieves wide velocity dynamic range and velocity‐to‐noise ratio (VNR), enabling accurate neurovascular flow characterization. To reduce scan time, we present interleaved dual‐venc 4D Flow with independently prescribed, prospectively undersampled spatial resolution of the high‐venc (HV) acquisition: Variable Spatial Resolution Dual Venc (VSRDV).

Methods

A prototype VSRDV sequence was developed based on a Cartesian acquisition with eight‐point phase encoding, combining PEAK‐GRAPPA acceleration with zero‐filling in phase and partition directions for HV. The VSRDV approach was optimized by varying z, the zero‐filling fraction of HV relative to low‐venc, between 0%–80% in vitro (realistic neurovascular model with pulsatile flow) and in vivo (n = 10 volunteers). Antialiasing precision, mean and peak velocity quantification accuracy, and test–retest reproducibility were assessed relative to reference images with equal‐resolution HV and low venc (z = 0%).

Results

In vitro results for all z demonstrated an antialiasing true positive rate at least 95% for RPEAK−GRAPPA = 2 and 5, with no linear relationship to z (p = 0.62 and 0.13, respectively). Bland–Altman analysis for z = 20%, 40%, 60%, or 80% versus z = 0% in vitro and in vivo demonstrated no bias >1% of venc in mean or peak velocity values at any RZF. In vitro mean and peak velocity, and in vivo peak velocity, had limits of agreement within 15%.

Conclusion

VSRDV allows up to 34.8% scan time reduction compared to PEAK‐GRAPPA accelerated DV 4D Flow MRI, enabling large spatial coverage and dynamic range while maintaining VNR and velocity measurement accuracy.

Click here for author‐reader discussions

Turbulent Intensity of Blood Flow in the Healthy Aorta Increases With Dobutamine Stress and is Related to Cardiac Output

Introduction: The blood flow in the normal cardiovascular system is predominately laminar but operates close to the threshold to turbulence. Morphological distortions such as vascular and valvular stenosis can cause transition into turbulent blood flow, which in turn may cause damage to tissues in the cardiovascular system. A growing number of studies have used magnetic resonance imaging (MRI) to estimate the extent and degree of turbulent flow in different cardiovascular diseases. However, the way in which heart rate and inotropy affect turbulent flow has not been investigated. In this study we hypothesized that dobutamine stress would result in higher turbulence intensity in the healthy thoracic aorta.

Method: 4D flow MRI data were acquired in twelve healthy subjects at rest and with dobutamine, which was infused until the heart rate increased by 60% when compared to rest. A semi-automatic segmentation method was used to segment the thoracic aorta in the 4D flow MR images. Subsequently, flow velocity and several turbulent kinetic energy (TKE) parameters were calculated in the ascending aorta, aortic arch, descending aorta and whole thoracic aorta.

Results: With dobutamine infusion there was an increase in heart rate (66 ± 9 vs. 108 ± 13 bpm, p < 0.001) and stroke volume (88 ± 13 vs. 102 ± 25 ml, p < 0.01). Additionally, there was an increase in Peak Average velocity (0.7 ± 0.1 vs. 1.2 ± 0.2 m/s, p < 0.001, Peak Max velocity (1.3 ± 0.1 vs. 2.0 ± 0.2 m/s, p < 0.001), Peak Total TKE (2.9 ± 0.7 vs. 8.0 ± 2.2 mJ, p < 0.001), Peak Median TKE (36 ± 7 vs. 93 ± 24 J/m3, p = 0.002) and Peak Max TKE (176 ± 33 vs. 334 ± 69 J/m3, p < 0.001). The relation between cardiac output and Peak Total TKE in the whole thoracic aorta was very strong (R² = 0.90, p < 0.001).

Conclusion: TKE of blood flow in the healthy thoracic aorta increases with dobutamine stress and is strongly related to cardiac output. Quantification of such turbulence intensity parameters with cardiac stress may serve as a risk assessment of aortic disease development.

Deep learning-based velocity antialiasing of 4D-flow MRI

Purpose

To develop a convolutional neural network (CNN) for the robust and fast correction of velocity aliasing in 4D‐flow MRI.

Methods

This study included 667 adult subjects with aortic 4D‐flow MRI data with existing velocity aliasing (n = 362) and no velocity aliasing (n = 305). Additionally, 10 controls received back‐to‐back 4D‐flow scans with systemically varied velocity‐encoding sensitivity (vencs) at 60, 100, and 175 cm/s. The no‐aliasing data sets were used to simulate velocity aliasing by reducing the venc to 40%–70% of the original, alongside a ground truth locating all aliased voxels (153 training, 152 testing). The 152 simulated and 362 existing aliasing data sets were used for testing and compared with a conventional velocity antialiasing algorithm. Dice scores were calculated to quantify CNN performance. For controls, the venc 175‐cm/s scans were used as the ground truth and compared with the CNN‐corrected venc 60 and 100 cm/s data sets

Results

The CNN required 176 ± 30 s to perform compared with 162 ± 14 s for the conventional algorithm. The CNN showed excellent performance for the simulated data compared with the conventional algorithm (median range of Dice scores CNN: [0.89–0.99], conventional algorithm: [0.84–0.94], p < 0.001, across all simulated vencs) and detected more aliased voxels in existing velocity aliasing data sets (median detected CNN: 159 voxels [31–605], conventional algorithm: 65 [7–417], p < 0.001). For controls, the CNN showed Dice scores of 0.98 [0.95–0.99] and 0.96 [0.87–0.99] for venc = 60 cm/s and 100 cm/s, respectively, while flow comparisons showed moderate‐excellent agreement.

Conclusion

Deep learning enabled fast and robust velocity anti‐aliasing in 4D‐flow MRI.

Clinical Feasibility of Structural and Functional MRI in Free-Breathing Neonates and Infants

BACKGROUND

Evaluation of structural lung abnormalities with magnetic resonance imaging (MRI) has previously been shown to be predictive of clinical neonatal outcomes in preterm birth. MRI during free-breathing with phase-resolved functional lung (PREFUL) may allow for complimentary functional information without exogenous contrast.

PURPOSE

To investigate the feasibility of structural and functional pulmonary MRI in a cohort of neonates and infants with no cardiorespiratory disease. Macrovascular pulmonary blood flows were also evaluated.

STUDY TYPE

Prospective.

POPULATION

Ten term infants with no clinically defined cardiorespiratory disease were imaged. Infants recruited from the general population and neonatal intensive care unit (NICU) were studied.

FIELD STRENGTH/SEQUENCE

T₁ -weighted VIBE, T₂ -weighted BLADE uncorrected for motion. Ultrashort echo time (UTE) and 3D-flow data were acquired during free-breathing with self-navigation and retrospective reconstruction. Single slice 2D-gradient echo (GRE) images were acquired during free-breathing for PREFUL analysis. Imaging was performed at 3 T.

ASSESSMENT

T₁ , T₂ , and UTE images were scored according to the modified Ochiai scheme by three pediatric body radiologists. Ventilation/perfusion-weighted maps were extracted from free-breathing GRE images using PREFUL analysis. Ventilation and perfusion defect percent (VDP, QDP) were calculated from the segmented ventilation and perfusion-weighted maps. Time-averaged cardiac blood velocities from three-dimensional-flow were evaluated in major pulmonary arteries and veins.

STATISTICAL TEST

Intraclass correlation coefficient (ICC).

RESULTS

The ICC of replicate structural scores was 0.81 (95% CI: 0.45-0.95) across three observers. Elevated Ochiai scores, VDP, and QDP were observed in two NICU participants. Excluding these participants, mean ± standard deviation structural scores were 1.2 ± 0.8, while VDP and QDP were 1.0% ± 1.1% and 0.4% ± 0.5%, respectively. Main pulmonary arterial blood flows normalized to body surface area were 3.15 ± 0.78 L/min/m² .

DATA CONCLUSION

Structural and functional pulmonary imaging is feasible using standard clinical MRI hardware (commercial whole-body 3 T scanner, table spine array, and flexible thoracic array) in free-breathing infants.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Multiple motion encoding in phase-contrast MRI: A general theory and application to elastography imaging

While MRI allows to encode the motion of tissue in the magnetization's phase, it remains yet a challenge to obtain high fidelity motion images due to wraps in the phase for high encoding efficiencies. Therefore, we propose an optimal multiple motion encoding method (OMME) and exemplify it in Magnetic Resonance Elastography (MRE) data. OMME is formulated as a non-convex least-squares problem for the motion using an arbitrary number of phase-contrast measurements with different motion encoding gradients (MEGs). The mathematical properties of OMME are proved in terms of standard deviation and dynamic range of the motion's estimate for arbitrary MEGs combination which are confirmed using synthetically generated data. OMME's performance is assessed on MRE data from in vivo human brain experiments and compared to dual encoding strategies. The unwrapped images are further used to reconstruct stiffness maps and compared to the ones obtained using conventional unwrapping methods. OMME allowed to successfully combine several MRE phase images with different MEGs, outperforming dual encoding strategies in either motion-to-noise ratio (MNR) or number of successfully reconstructed voxels with good noise stability. This lead to stiffness maps with greater resolution of details than obtained with conventional unwrapping methods. The proposed OMME method allows for a flexible and noise robust increase in the dynamic range and thus provides wrap-free phase images with high MNR. In MRE, the method may be especially suitable when high resolution images with high MNR are needed.

Hemodynamic force analysis is not ready for clinical trials on HFpEF

Hemodynamic force analysis has been proposed as a novel tool for early detection of subclinical systolic dysfunction in heart failure with preserved ejection fraction (HFpEF). Here we investigated the ability of hemodynamic forces to discriminate between healthy subjects and heart failure patients with varying degrees of systolic dysfunction. We studied 34 controls, 16 HFpEF patients, and 25 heart failure patients with mid-range (HFmrEF) or reduced ejection fraction (HFrEF) using cardiac magnetic resonance with acquisition of cine images and 4D flow at 1.5 T. The Navier-Stokes equation was used to compute global left ventricular hemodynamic forces over the entire cardiac cycle. Forces were analyzed for systole, diastole, and the entire heartbeat, with and without normalization to left ventricular volume. Volume-normalized hemodynamic forces demonstrated significant positive correlation with EF (r² = 0.47, p < 0.0001) and were found significantly lower in heart failure with reduced ejection fraction compared to controls (p < 0.0001 for systole and diastole). No difference was seen between controls and HFpEF (p > 0.34). Non-normalized forces displayed no differences between controls and HFpEF (p > 0.24 for all analyses) and did not correlate with EF (p = 0.36). Left ventricular hemodynamic force analysis, whether indexed to LV volumes or not, is not ready for clinical trials on HFpEF assessment.

Clinical Applications of 4D Flow MRI in the Portal Venous System

Evaluation of the hemodynamics in the portal venous system plays an essential role in many hepatic pathologies. Changes in portal flow and vessel morphology are often indicative of disease.Routinely used imaging modalities, such as CT, ultrasound, invasive angiography, and MRI, often focus on either hemodynamics or anatomical imaging. In contrast, 4D flow MRI facilitiates a more comprehensive understanding of pathophysiological mechanisms by simultaneously and noninvasively acquiring time-resolved flow and anatomical information in a 3D imaging volume.Though promising, 4D flow MRI in the portal venous system is especially challenging due to small vessel calibers, slow flow velocities, and breathing motion. In this review article, we will discuss how to account for these challenges when planning and conducting 4D flow MRI acquisitions in the upper abdomen. We will address patient preparation, sequence acquisition, postprocessing, quality control, and analysis of 4D flow data.In the second part of this article, we will review potential clinical applications of 4D flow MRI in the portal venous system. The most promising area for clinical utilization is the diagnosis and grading of liver cirrhosis and its complications. Relevant parameters acquired by 4D flow MRI include the detection of reduced or reversed flow in the portal venous system, characterization of portosystemic collaterals, and impaired response to a meal challenge. In patients with cirrhosis, 4D flow MRI has the potential to address the major unmet need of noninvasive detection of gastroesophageal varices at high risk for bleeding. This could replace many unnecessary, purely diagnostic, and invasive esophagogastroduodenoscopy procedures, thereby improving patient compliance with follow-up. Moreover, 4D flow MRI offers unique insights and added value for surgical planning and follow-up of multiple hepatic interventions, including transjugular intrahepatic portosystemic shunts, liver transplantation, and hepatic disease in children. Lastly, we will discuss the path to clinical implementation and remaining challenges.

The Impact of Aging on the Association Between Aortic Stiffness and Cerebral Pulsatility Index

The central arteries dampen the pulsatile forces from myocardial contraction, limiting the pulsatility that reaches the cerebral vasculature, although there are limited data on this relationship with aging in humans. The purpose of this study was to determine the association between aortic stiffness and cerebral artery pulsatility index in young and older adults. We hypothesized that cerebral pulsatility index would be associated with aortic stiffness in older adults, but not in young adults. We also hypothesized that both age and aortic stiffness would be significant predictors for cerebral pulsatility index. This study included 23 healthy older adults (aged 62 ± 6 years) and 33 healthy young adults (aged 25 ± 4 years). Aortic stiffness was measured using carotid-femoral pulse wave velocity (cfPWV), while cerebral artery pulsatility index in the internal carotid arteries (ICAs), middle cerebral arteries (MCAs), and basilar artery were assessed using 4D Flow MRI. Cerebral pulsatility index was calculated as (maximum flow - minimum flow) / mean flow. In the combined age group, there was a positive association between cfPWV and cerebral pulsatility index in the ICAs (r = 0.487; p < 0.001), MCAs (r = 0.393; p = 0.003), and basilar artery (r = 0.576; p < 0.001). In young adults, there were no associations between cfPWV and cerebral pulsatility index in any of the arteries of interest (ICAs: r = 0.253; p = 0.156, MCAs: r = -0.059; p = 0.743, basilar artery r = 0.171; p = 0.344). In contrast, in older adults there was a positive association between cfPWV and cerebral pulsatility index in the MCAs (r = 0.437; p = 0.037) and basilar artery (r = 0.500; p = 0.015). However, the relationship between cfPWV and cerebral pulsatility index in the ICAs of the older adults did not reach the threshold for significance (r = 0.375; p = 0.078). In conclusion, age and aortic stiffness are significant predictors of cerebral artery pulsatility index in healthy adults. This study highlights the importance of targeting aortic stiffness in our increasingly aging population to reduce the burden of age-related changes in cerebral hemodynamics.

Cerebrovascular stiffness and flow dynamics in the presence of amyloid and tau biomarkers

INTRODUCTION

This work investigated the relationship between cerebrovascular disease (CVD) markers and Alzheimer's disease (AD) biomarkers of amyloid beta deposition, and neurofibrillary tau tangles in subjects spanning the AD clinical spectrum.

METHODS

A total of 136 subjects participated in this study. Four groups were established based on AD biomarker positivity from positron emission tomography (amyloid [A] and tau [T]) and clinical diagnosis (cognitively normal [CN] and impaired [IM]). CVD markers were derived from structural and quantitative magnetic resonance imaging data.

RESULTS

Transcapillary pulse wave delay was significantly longer in controls compared to AT biomarker-confirmed groups (A+/T-/CN P < .001, A+/T+/CN P < .001, A+/T+/IM P = .003). Intracranial low-frequency oscillations were diminished in AT biomarker-confirmed groups both CN and impaired (A+/T-/CN P = .039, A+/T+/CN P = .007, A+/T+/IM P = .011). A significantly higher presence of microhemorrhages was measured in A+/T+/CN compared to controls (P = .006).

DISCUSSION

Cerebrovascular markers indicate increased vessel stiffness and reduced vasomotion in AT biomarker-positive subjects during preclinical AD.

Virtual injections using 4D flow MRI with displacement corrections and constrained probabilistic streamlines

PURPOSE

Streamlines from 4D-flow MRI have been used clinically for intracranial blood-flow tracking. However, deterministic and stochastic errors degrade streamline quality. The purpose of this study is to integrate displacement corrections, probabilistic streamlines, and novel fluid constraints to improve selective blood-flow tracking and emulate "virtual bolus injections."

METHODS

Both displacement artifacts (deterministic) and velocity noise (stochastic) inherently occur during phase-contrast MRI acquisitions. Here, two displacement correction methods, single-step and iterative, were tested in silico with simulated displacements and were compared with ground-truth velocity fields. Next, the effects of combining displacement corrections and constrained probabilistic streamlines were performed in 10 healthy volunteers using time-averaged 4D-flow data. Measures of streamline length and depth into vasculature were then compared with streamlines generated with no corrections and displacement correction alone using one-way repeated-measures analysis of variance and Friedman's tests. Finally, virtual injections with improved streamlines were generated for three intracranial pathology cases.

RESULTS

Iterative displacement correction outperformed the single-step method in silico. In volunteers, the combination of displacement corrections and constrained probabilistic streamlines allowed for significant improvements in streamline length and increased the number of streamlines entering the circle of Willis relative to streamlines with no corrections and displacement correction alone. In the pathology cases, virtual injections with improved streamlines were qualitatively similar to dynamic arterial spin labeling images and allowed for forward/reverse selective flow tracking to characterize cerebrovascular malformations.

CONCLUSION

Virtual injections with improved streamlines from 4D-flow MRI allow for flexible, robust, intracranial flow tracking.

Divergence-Free Constrained Phase Unwrapping and Denoising for 4D Flow MRI Using Weighted Least-Squares

A novel divergence-free constrained phase unwrapping method was proposed and evaluated for 4D flow MRI. The unwrapped phase field was obtained by integrating the phase variations estimated from the wrapped phase data using weighted least-squares. The divergence-free constraint for incompressible blood flow was incorporated to regulate and denoise the resulting phase field. The proposed method was tested on synthetic phase data of left ventricular flow and in vitro 4D flow measurement of Poiseuille flow. The method was additionally applied to in vivo 4D flow measurements in the thoracic aorta from 30 human subjects. The performance of the proposed method was compared to the state-of-the-art 4D single-step Laplacian algorithm. The synthetic phase data were completely unwrapped by the proposed method for all the cases with velocity encoding (venc) as low as 20% of the maximum velocity and signal-to-noise ratio as low as 5. The in vitro Poiseuille flow data were completely unwrapped with a 60% increase in the velocity-to-noise ratio. For the in-vivo aortic datasets with venc ratio less than 0.4, the proposed method significantly improved the success rate by as much as 40% and reduced the velocity error levels by a factor of 10 compared to the state-of-the-art method. The divergence-free constrained method exhibits reliability and robustness on phase unwrapping and shows improved accuracy of velocity and hemodynamic quantities by unwrapping the low-venc 4D flow MRI data.

Velocity vector reconstruction for real-time phase-contrast MRI with radial Maxwell correction

PURPOSE

To develop an auto-calibrated image reconstruction for highly accelerated multi-directional phase-contrast (PC) MRI that compensates for (1) reconstruction instabilities occurring for phase differences near ± π and (2) phase errors by concomitant magnetic fields that differ for individual radial spokes.

THEORY AND METHODS

A model-based image reconstruction for real-time PC MRI based on nonlinear inversion is extended to multi-directional flow by exploiting multiple flow-encodings for the estimation of velocity vectors. An initial smoothing constraint during iterative optimization is introduced to resolve the ambiguity of the solution space by penalizing phase wraps. Maxwell terms are considered as part of the signal model on a line-by-line basis to address phase errors by concomitant magnetic fields. The reconstruction methods are evaluated using simulated data and cross-sectional imaging of a rotating-disc, as well as in vivo for the aortic arch and cervical spinal canal at 3T.

RESULTS

Real-time three-directional velocity mapping in the aortic arch is achieved at 1.8 × 1.8 × 6 mm³ spatial and 60 ms temporal resolution. Artificial phase wraps are avoided in all cases using the smoothness constraint. Inter-spoke differences of concomitant magnetic fields are effectively compensated for by the model-based image reconstruction with integrated radial Maxwell correction.

CONCLUSION

Velocity vector reconstructions based on nonlinear inversion allow for high degrees of radial data undersampling paving the way for multi-directional PC MRI in real time. Whether a spoke-wise treatment of Maxwell terms is required or a computationally cheaper frame-wise approach depends on the individual application.

Whole-Heart 4D Flow MRI for Evaluation of Normal and Regurgitant Valvular Flow: A Quantitative Comparison Between Pseudo-Spiral Sampling and EPI Readout

Background

Pseudo‐spiral Cartesian sampling with compressed sensing reconstruction has facilitated highly accelerated 4D flow magnetic resonance imaging (MRI) in various cardiovascular structures. However, unlike echo planar imaging (EPI)‐accelerated 4D flow MRI, it has not been validated in whole‐heart applications.

Hypothesis

Pseudo‐spiral 4D flow MRI (PROUD [PROspective Undersampling in multiple Dimensions]) is comparable to EPI in robustness of valvular flow measurements and remains comparable as the undersampling factor is increased and scan time reduced.

Study Type

Prospective.

Population

Twelve healthy subjects and eight patients with valvular regurgitation.

Field Strength/Sequence

3.0 T; PROUD and EPI 4D flow sequences, 2D flow and balanced steady‐state free precession sequences.

Assessment

Valvular blood flow was quantified using valve tracking. PROUD‐ and EPI‐based measurements of aortic (AV) and pulmonary (PV) flow volumes and left and right ventricular stroke volumes were tested for agreement with 2D MRI‐based measurements. PROUD reconstructions with undersampling factors (R) of 9, 14, 28, and 56 were tested for intervalve consistency (per valve, compared to the other valves) and preservation of peak velocities and E/A ratios.

Statistical Tests

We used repeated measures ANOVA, Bland‐Altman, Wilcoxon signed rank, and intraclass correlation coefficients. P < 0.05 was considered statistically significant.

Results

PROUD and EPI intervalve consistencies were not significantly different both in healthy subjects (valve‐averaged mean difference [limits of agreement width]: 3.2 ± 0.8 [8.7 ± 1.1] mL/beat for PROUD, 5.5 ± 2.9 [13.7 ± 2.3] mL/beat for EPI, P = 0.07) and in patients with valvular regurgitation (2.3 ± 1.2 [15.3 ± 5.9] mL/beat for PROUD, 0.6 ± 0.6 [19.3 ± 2.9] mL/beat for EPI, P = 0.47). Agreement between EPI and PROUD was higher than between 4D flow (EPI or PROUD) and 2D MRI for forward flow, stroke volumes, and regurgitant volumes. Up to R = 28 in healthy subjects and R = 14 in patients with valvular regurgitation, PROUD intervalve consistency remained comparable to that of EPI. Peak velocities and E/A ratios were preserved up to R = 9.

Conclusion

PROUD is comparable to EPI in terms of intervalve consistency and may be used with higher undersampling factors to shorten scan times further.

Level of Evidence

1

Technical Efficacy Stage

2

Open or closed: Changes in ductus arteriosus flow patterns at birth using 4D flow MRI in newborn piglets

The ductus arteriosus (DA) functionally closes during the transition from fetal to postnatal life because of lung aeration and corresponding cardiovascular changes. The thorough and explicit measurement and visualization of the right and left heart output during this transition has not been previously accomplished. We combined 4D flow MRI (dynamic volumetric blood flow measurements) and T2 relaxometry (oxygen delivery quantification) in surgically instrumented newborn piglets to assess the DA. This was performed in Large White-Landrace-Duroc piglets (n = 34) spanning four age groups: term-9 days, term-3, term+1, and term+5. Subject's DA status was classified using 4D flow: closed DA, forward DA flow, reversed DA flow, and bidirectional DA flow. Over all states, vessel diameters and flows normalized to body weight increased with age (for example in the ascending aorta flow-term-9: 126.6 ± 45.4; term+5: 260.2 ± 80.0 ml/min per kg; p = 0.0005; ascending aorta diameter-term-9: 5.2 ± 0.8; term+5: 7.7 ± 0.4 mm; p = 0.0004). In subjects with reversed DA blood flow there was lower common carotid artery blood flow (forward: 37.5 ± 9.0; reversed: 20.0 ± 7.4 ml/min per kg; p = 0.032). Linear regression revealed that as net DA flow decreases, common carotid artery flow decreases (R2  = 0.32, p = 0.004), and left (R2  = 0.33, p = 0.003) and right (R2  = 0.34, p = 0.003) pulmonary artery flow increases. Bidirectional DA blood flow changed oxygen saturation as determined by MRI between the ascending and descending aorta (ascending aorta: 90.1% ± 8.4%; descending aorta: 75.6% ± 14.2%; p < 0.05). Expanded use of these techniques in preterm animal models will aid in providing new understandings of normal versus abnormal DA transition, as well as to test the effectiveness of related clinical interventions.

Structural Heart 4D Flow MRI for Hemodynamic Assessment: How We Do It

MRI is an essential diagnostic tool in the anatomic and functional evaluation of cardiovascular disease. In many practices, 2D phase-contrast (2D-PC) has been used for blood flow quantification. 4D Flow MRI is a time-resolved volumetric acquisition that captures the vector field of blood flow along with anatomic images. 4D Flow MRI provides a simpler acquisition compared to 2D-PC and facilitates a more accurate and comprehensive hemodynamic assessment. Advancements in accelerated imaging have significantly shortened scan times of 4D Flow MRI while preserving image quality, enabling this technology to transition from the research arena to routine clinical practice. In this article, we review technical optimization based on our clinical experience of over 10 years with 4D Flow MRI. We also present pearls and pitfalls in the practical application of 4D Flow MRI, including how to quantify cardiovascular shunts, valvular or vascular stenosis, and valvular regurgitation. As experience increases, and as 4D Flow sequences and post-processing software become more broadly available, 4D Flow MRI will likely become an essential component of cardiac imaging for practices involved in the management of congenital and acquired structural heart disease.

Left atrial 4D flow cardiovascular magnetic resonance: a reproducibility study in sinus rhythm and atrial fibrillation

BACKGROUND

Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) allows sophisticated quantification of left atrial (LA) blood flow, and could yield novel biomarkers of propensity for intra-cardiac thrombus formation and embolic stroke. As reproducibility is critically important to diagnostic performance, we systematically investigated technical and temporal variation of LA 4D flow in atrial fibrillation (AF) and sinus rhythm (SR).

METHODS

Eighty-six subjects (SR, n = 64; AF, n = 22) with wide-ranging stroke risk (CHA2DS2VASc 0-6) underwent LA 4D flow assessment of peak and mean velocity, vorticity, vortex volume, and stasis. Eighty-five (99%) underwent a second acquisition within the same session, and 74 (86%) also returned at 30 (27-35) days for an interval scan. We assessed variability attributable to manual contouring (intra- and inter-observer), and subject repositioning and reacquisition of data, both within the same session (same-day scan-rescan), and over time (interval scan). Within-subject coefficients of variation (CV) and bootstrapped 95% CIs were calculated and compared.

RESULTS

Same-day scan-rescan CVs were 6% for peak velocity, 5% for mean velocity, 7% for vorticity, 9% for vortex volume, and 10% for stasis, and were similar between SR and AF subjects (all p > 0.05). Interval-scan variability was similar to same-day scan-rescan variability for peak velocity, vorticity, and vortex volume (all p > 0.05), and higher for stasis and mean velocity (interval scan CVs of 14% and 8%, respectively, both p < 0.05). Longitudinal changes in heart rate and blood pressure at the interval scan in the same subjects were associated with significantly higher variability for LA stasis (p = 0.024), but not for the remaining flow parameters (all p > 0.05). SR subjects showed significantly greater interval-scan variability than AF patients for mean velocity, vortex volume, and stasis (all p < 0.05), but not peak velocity or vorticity (both p > 0.05).

CONCLUSIONS

LA peak velocity and vorticity are the most reproducible and temporally stable novel LA 4D flow biomarkers, and are robust to changes in heart rate, blood pressure, and differences in heart rhythm.

Retrospective Camera-Based Respiratory Gating in Clinical Whole-Heart 4D Flow MRI

Background

Respiratory gating is generally recommended in 4D flow MRI of the heart to avoid blurring and motion artifacts. Recently, a novel automated contact‐less camera‐based respiratory motion sensor has been introduced.

Purpose

To compare camera‐based respiratory gating (CAM) with liver‐lung‐navigator‐based gating (NAV) and no gating (NO) for whole‐heart 4D flow MRI.

Study Type

Retrospective.

Subjects

Thirty two patients with a spectrum of cardiovascular diseases.

Field Strength/Sequence

A 3T, 3D‐cine spoiled‐gradient‐echo‐T1‐weighted‐sequence with flow‐encoding in three spatial directions.

Assessment

Respiratory phases were derived and compared against each other by cross‐correlation. Three radiologists/cardiologist scored images reconstructed with camera‐based, navigator‐based, and no respiratory gating with a 4‐point Likert scale (qualitative analysis). Quantitative image quality analysis, in form of signal‐to‐noise ratio (SNR) and liver‐lung‐edge (LLE) for sharpness and quantitative flow analysis of the valves were performed semi‐automatically.

Statistical Tests

One‐way repeated measured analysis of variance (ANOVA) with Wilks's lambda testing and follow‐up pairwise comparisons. Significance level of P ≤ 0.05. Krippendorff's‐alpha‐test for inter‐rater reliability.

Results

The respiratory signal analysis revealed that CAM and NAV phases were highly correlated (C = 0.93 ± 0.09, P < 0.01). Image scoring showed poor inter‐rater reliability and no significant differences were observed (P ≥ 0.16). The image quality comparison showed that NAV and CAM were superior to NO with higher SNR (P = 0.02) and smaller LLE (P < 0.01). The quantitative flow analysis showed significant differences between the three respiratory‐gated reconstructions in the tricuspid and pulmonary valves (P ≤ 0.05), but not in the mitral and aortic valves (P > 0.05). Pairwise comparisons showed that reconstructions without respiratory gating were different in flow measurements to either CAM or NAV or both, but no differences were found between CAM and NAV reconstructions.

Data Conclusion

Camera‐based respiratory gating performed as well as conventional liver‐lung‐navigator‐based respiratory gating. Quantitative image quality analysis showed that both techniques were equivalent and superior to no‐gating‐reconstructions. Quantitative flow analysis revealed local flow differences (tricuspid/pulmonary valves) in images of no‐gating‐reconstructions, but no differences were found between images reconstructed with camera‐based and navigator‐based respiratory gating.

Level of Evidence

3

Technical Efficacy

Stage 2

Vertebral artery hypoplasia influences age-related differences in blood flow of the large intracranial arteries

Our purpose was to compare cerebral blood flow in the large intracranial vessels between healthy adults with (VAH+) and without (No VAH) vertebral artery hypoplasia. We also evaluated age-related differences in regional blood flow through the large cerebral arteries. Healthy young (n = 20; age = 25 ± 3 years) and older adults (n = 19; age = 61 ± 5 years) underwent 4D flow MRI scans to evaluate blood flow in the internal carotid arteries (ICA) and basilar artery (BA). VAH was determined retrospectively from 4D flow MRI using both structural (vessel diameter ≤ 2 mm) and flow criteria (flow ≤ 50 mL/min). We identified 5 young and 5 older adults with unilateral VAH (prevalence = 26%). ICA flow was lower in the VAH+ group compared with the No VAH group (367 ± 75 mL/min vs. 432 ± 92 mL/min, respectively; p < 0.05). There was no difference in BA flow between VAH+ and No VAH (110 ± 20 mL/min vs. 126 ± 40 mL/min, respectively; p = 0.24). When comparing age-related differences in blood flow in the No VAH group, older adults demonstrated lower BA flow compared with young adults (111 ± 38 mL/min vs. 140 ± 38 mL/min, respectively; p < 0.05) but not ICA flow (428 ± 89 mL/min vs. 436 ± 98 mL/min, respectively; p = 0.82). In contrast, in the VAH+ group, older adults had lower ICA flow compared with young adults (312 ± 65 mL/min vs. 421 ± 35 mL/min, respectively; p < 0.01), but not BA flow (104 ± 16 mL/min vs. 117 ± 23 mL/min, respectively; p = 0.32). Our results suggest that the presence of VAH is associated with lower ICA blood flow. Furthermore, VAH may contribute to the variability in the age-related differences in cerebral blood flow in healthy adults.

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.

A Dual-VENC Four-Dimensional Flow MRI Framework for Analysis of Subject-Specific Heterogeneous Nonlinear Vessel Deformation

Advancement of subject-specific in silico medicine requires new imaging protocols tailored to specific anatomical features, paired with new constitutive model development based on structure/function relationships. In this study, we develop a new dual-velocity encoding coefficient (VENC) 4D flow MRI protocol that provides unprecedented spatial and temporal resolution of in vivo aortic deformation. All previous dual-VENC 4D flow MRI studies in the literature focus on an isolated segment of the aorta, which fail to capture the full spectrum of aortic heterogeneity that exists along the vessel length. The imaging protocol developed provides high sensitivity to all blood flow velocities throughout the entire cardiac cycle, overcoming the challenge of accurately measuring the highly unsteady nonuniform flow field in the aorta. Cross-sectional area change, volumetric flow rate, and compliance are observed to decrease with distance from the heart, while pulse wave velocity (PWV) is observed to increase. A nonlinear aortic lumen pressure-area relationship is observed throughout the aorta such that a high vessel compliance occurs during diastole, and a low vessel compliance occurs during systole. This suggests that a single value of compliance may not accurately represent vessel behavior during a cardiac cycle in vivo. This high-resolution MRI data provide key information on the spatial variation in nonlinear aortic compliance, which can significantly advance the state-of-the-art of in-silico diagnostic techniques for the human aorta.

Technique for comprehensive fetal hepatic blood flow assessment in sheep using 4D flow MRI

KEY POINTS

The comprehensive visualization and quantification of in vivo fetal hepatic haemodynamics, particularly the shunting of ductus venosus blood, has been elusive and is not yet fully understood. We introduce the combination of chronically instrumented fetal sheep and 4D flow MRI of the whole fetal liver, which allows retrospective blood flow measurement in all visible vessels as well as qualitative assessment. The applicability and usefulness of this technique is exhibited in normally grown fetal Merino sheep in mid- and late-gestation with detailed dynamic distribution of hepatic blood flow presented. The feasibility of this approach in clinical pathology is demonstrated in two growth-restricted fetuses at mid-gestation. Further exemplification of blood flow quantification is performed over major hepatic vessels.

ABSTRACT

Although the fetal vasculature has been demarcated and well understood for several decades, the corresponding haemodynamics permitting oxygen- and nutrient-rich blood delivery to the fetal organs has been comparatively difficult to study. We married two well-established methods: 4D flow MRI, a volumetric and dynamic blood-flow measurement technique, and chronically instrumented sheep to broadly assess fetal hepatic circulation. We performed this technique in mid- and late-gestation fetal Merino sheep under normoxemic conditions and major hepatic vasculature was segmented to quantify blood flow and related parameters. Dynamic blood flow was visualized, exhibiting an acceleration of umbilical vein blood through the ductus venosus as well as spiralling into the inferior vena cava where its stream remained separate from that of the hepatic veins and lower body. Ductus venosus changes from mid- to late-gestation included larger diameter (mid: 5.8 ± 0.9 vs. late: 7.1 ± 1.1 mm; P = 0.003) and cross-sectional area (mid: 27.1 ± 8.6 vs. late: 40.4 ± 11.8 mm2 ; P = 0.003), and lower velocity averaged over the cardiac cycle (mid: 15.7 ± 5.4 vs. late: 9.8 ± 7.0 cm s-1 ; P = 0.020). This resulted in higher magnitude blood flow (indexed to umbilical vein input) at mid-gestation in the ductus venosus (mid: 0.73 ± 0.21; late: 0.46 ± 0.21; P = 0.008). The visualization and quantification results support the further use of this technique to better understand regional blood flow changes during normal or abnormal fetal growth, as well as to observe acute haemodynamic responses to physiological challenges or drug interventions.

Phase unwrapping based on adaptive image in-painting of fringe patterns in measuring gear tooth flanks by laser interferometry

Phase unwrapping in regions of abnormal fringes remains an unresolved issue. In this paper, we present an approach that combines an image-inpainting strategy based on an adaptive window which is obtained according to the density and orientation of fringe patterns and a quality-guided algorithm for phase unwrapping. First, a threshold is set to a quality map to detect the target region. Second, the target region is filled with new phase values by the adaptive image-inpainting method. Then, a quality-guided phase unwrapping algorithm is applied to this newly generated wrapped phase map. Finally, postprocessing of the unwrapped result is performed. The method is validated through several simulation and experiments. The results demonstrate that the proposed algorithm is effective in the presence of abnormal fringes.

Simultaneous Multi-VENC and Simultaneous Multi-Slice Phase Contrast Magnetic Resonance Imaging

This work develops a novel, simultaneous multi-VENC and simultaneous multi-slice (SMV+SMS) imaging in a single acquisition for robust phase contrast (PC) MRI. To this end, the pulse sequence was designed to permit concurrent acquisition of multiple VENCs as well as multiple slices on a shared frequency encoding gradient, in which each effective echo time for multiple VENCs was controlled by adjusting net gradient area while multiple slices were simultaneously excited by employing multiband resonance frequency (RF) pulses. For VENC and slice separation, RF phase cycling and gradient blip were applied to create both inter-VENC and inter-slice shifts along phase encoding direction, respectively. With an alternating RF phase cycling that generates oscillating steady-state with low and high signal amplitude, the acquired multi-VENC k-space was reformulated into 3D undersampled k-space by generating a virtual dimension along VENC direction for modulation induced artifact reduction. In vivo studies were conducted to validate the feasibility of the proposed method in comparison with conventional PC MRI. The proposed method shows comparable performance to the conventional method in delineating both low and high flow velocities across cardiac phases with high spatial coverage without apparent artifacts. In the presence of high flow velocity that is above the VENC value, the proposed method exhibits clear depiction of flow signals over conventional method, thereby leading to high VNR image with improved velocity dynamic range.

Age-Related Reductions in Cerebrovascular Reactivity Using 4D Flow MRI

Cerebrovascular reactivity (CVR), is important for determining future risk of cerebrovascular disease. It is unclear if primary aging is associated with reductions in CVR because previous studies often include participants with vascular risk factors. Additionally, the inconsistency in the literature may be due to the inherent difficulty in quantifying intracranial cerebral blood flow and CVR. To address these limitations, we determined the effect of age on CVR in the large intracranial vessels in adults with low vascular risk using state-of-the-art MRI techniques. We also determined if the effect of age on CVR was sex-specific. Young (n = 20; 25 ± 3 years) and older (n = 19; 61 ± 5 years) healthy, physically active adults participated in the study. CVR was measured in response to hypercapnia using 4D flow MRI, which allows for simultaneous angiographic and quantitative blood flow measurements in the intracranial arteries. Older adults had lower global CVR and CVR in multiple intracranial arteries [right and left internal carotid arteries (ICA), right and left middle cerebral arteries (MCA), and basilar artery (BA)] compared with young adults (p < 0.05 for all). In addition, the MCA dilated significantly in response to hypercapnia in young (p < 0.05), but not older adults. Young men demonstrated higher global CVR and CVR in multiple intracranial arteries (ICAs, MCAs, and BA) compared with young women and older men (p < 0.05 for both); however, CVR did not differ between young women and older women. Our results demonstrate that, using 4D flow MRI, primary aging is associated with lower CVR in adults with low vascular risk. In addition, the effect of age on CVR may be driven by men. The 4D flow MRI technique may provide a promising new alternative to measure cerebrovascular physiology without the limitations of commonly used techniques. Future studies could utilize this MRI technique to examine interventions to maintain CVR with advancing age. This study was registered under clinicaltrials.gov # NCT02840851.

Motion robust respiratory-resolved 3D radial flow MRI and its application in neonatal congenital heart disease

PURPOSE

To test and implement a motion-robust and respiratory-resolved 3D Radial Flow framework that addresses the need for rapid, high resolution imaging in neonatal patients with congenital heart disease.

METHODS

A 4-point velocity encoding and 3D radial trajectory with double-golden angle ordering was combined with bulk motion correction (from projection center of mass) and respiration phase detection (from principal component analysis of heartbeat-averaged data) to create motion-robust 3D velocity cardiac time-averaged data. This framework was tested in a whole-chest digital phantom with simulated bulk and realistic physiological motion. In vivo imaging was performed in 20 congenital heart disease infants under feed-and-sleep with submillimeter isotropic resolution in ~3 min. Flows were validated against clinical 2D PCMRI and whole-heart visualizations of blood flow were performed.

RESULTS

The proposed framework resolved all simulated digital phantom motion states (mean ± standard error: rotation - azimuthal = 0.29 ± 0.02°; translation - T_y = 1.29 ± 0.12 mm, T_z = -0.27 ± 0.13 mm; rotation+translation - polar = 0.49 ± 0.16°, T_x = -2.47 ± 0.51 mm, T_z = 5.78 ± 1.33 mm). Measured timing errors of peak expiration across all signal-to-noise ratio values were 22% of the true respiratory period (range = [404-489 ± 298-334] ms). For in vivo imaging, motion correction improved 3D Radial Flow measurements (no correction: R² = 0.62, root mean square error = 0.80 L/min/m² , Bland-Altman bias [limits of agreement] = -0.21 [-1.40, 0.94] L/min/m² ; motion corrected, expiration: R² = 0.90, root mean square error = 0.46 L/min/m² , bias [limits of agreement] = 0.06 [-0.49, 0.62] L/min/m² ). Respiratory-resolved 3D velocity visualizations were achieved in various neonatal pathologies pre- and postsurgical correction.

CONCLUSION

3D cardiac flow may be visualized and accurately quantified in neonatal subjects using the proposed framework. This technique may enable more comprehensive hemodynamic studies in small infants.

Efficient triple-VENC phase-contrast MRI for improved velocity dynamic range

PURPOSE

To evaluate the utility of an efficient triple velocity-encoding (VENC) 4D flow MRI implementation to improve velocity unwrapping of 4D flow MRI data with the same scan time as an interleaved dual-VENC acquisition.

METHODS

A balanced 7-point acquisition was used to derive 3 sets of 4D flow images corresponding to 3 different VENCs. These 3 datasets were then used to unwrap the aliased lowest VENC into a minimally aliased, triple-VENC dataset. Triple-VENC MRI was evaluated and compared with dual-VENC MRI over 3 different VENC ranges (50-150, 60-150, and 60-180 cm/s) in vitro in a steadily rotating phantom as well as in a pulsatile flow phantom. In vivo, triple-VENC data of the thoracic aorta were also evaluated in 3 healthy volunteers (2 males, 26-44 years old) with VENC = 50/75/150 cm/s. Two triple-VENC (triconditional and biconditional) and 1 dual-VENC unwrapping algorithms were quantitatively assessed through comparison to a reference, unaliased, single-VENC scan.

RESULTS

Triple-VENC 4D flow constant rotation phantom results showed high correlation with the analytical solution (intraclass correlation coefficient = 0.984-0.995, P < .001) and up to a 61% reduction in velocity noise compared with the corresponding single-VENC scans (VENC = 150, 180 cm/s). Pulsatile flow phantom experiments demonstrated good agreement between triple-VENC and single-VENC acquisitions (peak flow < 0.8% difference; peak velocity < 11.7% difference). Triconditional triple-VENC unwrapping consistently outperformed dual-VENC unwrapping, correctly unwrapping more than 83% and 46%-66% more voxels in vitro and in vivo, respectively.

CONCLUSION

Triple-VENC 4D flow MRI adds no additional scan time to dual-VENC MRI and has the potential for improved unwrapping to extend the velocity dynamic range beyond dual-VENC methods.

Optimal Dual-VENC Unwrapping in Phase-Contrast MRI

Dual-VENC strategies have been proposed to improve the velocity-to-noise ratio in phase-contrast MRI. However, they are based on aliasing-free high-VENC data. The aim of this paper is to propose a dual-VENC velocity estimation method allowing high-VENC aliased data. For this purpose, we reformulate the phase-contrast velocity as a least squares estimator, providing a natural framework for including multiple encoding gradient measurements. By analyzing the mathematical properties of both single- and dual-VENC problems, we can justify theoretically high/low-VENC ratios such that the aliasing velocity can be minimized. The resulting reconstruction algorithm was assessed using three types of data: numerical, experimental, and volunteers. In clinical practice, this method would allow shorter examination times by avoiding tedious adaptation of VENC values by repeated scans.

Computationally effective 2D and 3D fast phase unwrapping algorithms and their applications to Doppler optical coherence tomography

We propose a simplification for a robust and easy to implement fast phase unwrapping (FPU) algorithm that is used to solve the phase wrapping problem encountered in various fields of optical imaging and metrology. We show that the number of necessary computations using the algorithm can be reduced compared to its original version. FPU can be easily extended from two to three spatial dimensions. We demonstrate the applicability of the two- and three-dimensional FPU algorithm for Doppler optical coherence tomography (DOCT) in numerical simulations, and in the imaging of a flow phantom and blood flow in the human retina in vivo. We introduce an FPU applicability plot for use as a guide in the selection of the most suitable version of the algorithm depending on the phase noise in the acquired data. This plot uses the circular standard deviation of the wrapped phase distribution as a measure of noise and relates it to the root-mean-square error of the recovered, unwrapped phase.

Validation and reproducibility of cardiovascular 4D-flow MRI from two vendors using 2 × 2 parallel imaging acceleration in pulsatile flow phantom and in vivo with and without respiratory gating

Background

4D-flow magnetic resonance imaging (MRI) is increasingly used.

Purpose

To validate 4D-flow sequences in phantom and in vivo, comparing volume flow and kinetic energy (KE) head-to-head, with and without respiratory gating.

Material and Methods

Achieva dStream (Philips Healthcare) and MAGNETOM Aera (Siemens Healthcare) 1.5-T scanners were used. Phantom validation measured pulsatile, three-dimensional flow with 4D-flow MRI and laser particle imaging velocimetry (PIV) as reference standard. Ten healthy participants underwent three cardiac MRI examinations each, consisting of cine-imaging, 2D-flow (aorta, pulmonary artery), and 2 × 2 accelerated 4D-flow with (Resp+) and without (Resp−) respiratory gating. Examinations were acquired consecutively on both scanners and one examination repeated within two weeks. Volume flow in the great vessels was compared between 2D- and 4D-flow. KE were calculated for all time phases and voxels in the left ventricle.

Results

Phantom results showed high accuracy and precision for both scanners. In vivo, higher accuracy and precision (P < 0.001) was found for volume flow for the Aera prototype with Resp+ (–3.7 ± 10.4 mL, r = 0.89) compared to the Achieva product sequence (–17.8 ± 18.6 mL, r = 0.56). 4D-flow Resp− on Aera had somewhat larger bias (–9.3 ± 9.6 mL, r = 0.90) compared to Resp+ (P = 0.005). KE measurements showed larger differences between scanners on the same day compared to the same scanner at different days.

Conclusion

Sequence-specific in vivo validation of 4D-flow is needed before clinical use. 4D-flow with the Aera prototype sequence with a clinically acceptable acquisition time (<10 min) showed acceptable bias in healthy controls to be considered for clinical use. Intra-individual KE comparisons should use the same sequence.

Fetal hemodynamics and cardiac streaming assessed by 4D flow cardiovascular magnetic resonance in fetal sheep

BACKGROUND

To date it has not been possible to obtain a comprehensive 3D assessment of fetal hemodynamics because of the technical challenges inherent in imaging small cardiac structures, movement of the fetus during data acquisition, and the difficulty of fusing data from multiple cardiac cycles when a cardiac gating signal is absent. Here we propose the combination of volumetric velocity-sensitive cardiovascular magnetic resonance imaging ("4D flow" CMR) and a specialized animal preparation (catheters to monitor fetal heart rate, anesthesia to immobilize mother and fetus) to examine fetal sheep cardiac hemodynamics in utero.

METHODS

Ten pregnant Merino sheep underwent surgery to implant arterial catheters in the target fetuses. Anesthetized ewes underwent 4D flow CMR with acquisition at 3 T for fetal whole-heart coverage with 1.2-1.5 mm spatial resolution and 45-62 ms temporal resolution. Flow was measured in the heart and major vessels, and particle traces were used to visualize circulatory patterns in fetal cardiovascular shunts. Conservation of mass was used to test internal 4D flow consistency, and comparison to standard 2D phase contrast (PC) CMR was performed for validation.

RESULTS

Streaming of blood from the ductus venosus through the foramen ovale was visualized. Flow waveforms in the major thoracic vessels and shunts displayed normal arterial and venous patterns. Combined ventricular output (CVO) was 546 mL/min per kg, and the distribution of flows (%CVO) were comparable to values obtained using other methods. Internal 4D flow consistency across 23 measurement locations was established with differences of 14.2 ± 12.1%. Compared with 2D PC CMR, 4D flow showed a strong correlation (R2 = 0.85) but underestimated flow (bias = - 21.88 mL/min per kg, p < 0.05).

CONCLUSIONS

The combination of fetal surgical preparation and 4D flow CMR enables characterization and quantification of complex flow patterns in utero. Visualized streaming of blood through normal physiological shunts confirms the complex mechanism of substrate delivery to the fetal heart and brain. Besides offering insight into normal physiology, this technology has the potential to qualitatively characterize complex flow patterns in congenital heart disease phenotypes in a large animal model, which can support the development of new interventions to improve outcomes in this population.