Semiautomated method for noise reduction and background phase error correction in MR phase velocity data

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.

On the limitations of partial Fourier acquisition in phase-contrast MRI of turbulent kinetic energy

PURPOSE

To investigate limitations of partial Fourier acquisition in phase-contrast MRI of turbulent kinetic energy (TKE).

METHODS

To assess the validity of partial Fourier reconstruction of TKE and phase images, computational fluid dynamics data of mean and turbulent velocities in a stenotic U-bend phantom was used. Partial Fourier acquisition with 75% k-space coverage was simulated and TKE data were reconstructed using zero-filling, homodyne reconstruction, and the method of projections onto convex sets (POCS). Results were compared to data from fully sampled k-space and 75% symmetric sampling. In addition, compressed sensing (CS) reconstruction was compared for a standard variable density sampling pattern and a variable density sampling pattern combined with 75% partial Fourier. For illustration purposes, in vivo examples of velocity magnitude and TKE maps of aortic flow reconstructed with the different methods are provided.

RESULTS

In accordance with theory, partial Fourier reconstruction of TKE maps from phase-contrast data results in artifacts relative to fully sampled data. It is demonstrated that neither homodyne reconstruction nor POCS can improve reconstruction of TKE data with respect to zero-filling reconstruction when compared to ground-truth (RMS error: 4.70%, 4.34%, and 2.45% for homodyne, POCS, and zero-filling reconstruction of in vivo data, respectively). CS reconstruction from data acquired with partial Fourier did not recover the resolution loss incurred by partial Fourier sampling.

CONCLUSION

Partial Fourier reconstruction of TKE maps from phase-contrast data does not yield a benefit over zero-filling reconstruction. In consequence, symmetric sampling is preferred over partial Fourier acquisition for a given number of phase-encodes in phase-contrast MRI.

Impact of Age and Diastolic Function on Novel, 4D flow CMR Biomarkers of Left Ventricular Blood Flow Kinetic Energy

Two-dimensional (2D) methods of assessing mitral inflow velocities are pre-load dependent, limiting their reliability for evaluating diastolic function. Left ventricular (LV) blood flow kinetic energy (KE) derived from four-dimensional flow cardiovascular magnetic resonance imaging (4D flow CMR) may offer improvements. It remains unclear whether 4D LV blood flow KE parameters are associated with physiological factors, such as age when compared to 2D mitral inflow velocities. Fifty-three healthy volunteers underwent standard CMR, plus 4D flow acquisition. LV blood flow KE parameters demonstrated good reproducibility with mean coefficient of variation of 6 ± 2% and an accuracy of 99% with a precision of 97%. The LV blood flow KEi_EDV E/A ratio demonstrated good association to the 2D mitral inflow E/A ratio (r = 0.77, P < 0.01), with both decreasing progressively with advancing age (P < 0.01). Furthermore, peak E-wave KEi_EDV and A-wave KEi_EDV displayed a stronger association to age than the corresponding 2D metrics, peak E-wave and A-wave velocity (r = −0.51 vs −0.17 and r = 0.65 vs 0.46). Peak E-wave KEi_EDV decreases whilst peak A-wave KEi_EDV increases with advancing age. This study presents values for various LV blood flow KE parameters in health, as well as demonstrating that they show stronger and independent correlations to age than standard diastolic metrics.

Caval to pulmonary 3D flow distribution in patients with Fontan circulation and impact of potential 4D flow MRI error sources

PURPOSE

Uneven flow distribution in patients with Fontan circulation is suspected to lead to complications. 4D flow MRI offers evaluation using time-resolved pathlines; however, the potential error is not well understood. The aim of this study was to systematically assess variability in flow distribution caused by well-known sources of error.

METHODS

4D flow MRI was acquired in 14 patients with Fontan circulation. Flow distribution was quantified by the % of caval venous flow pathlines reaching the left and right pulmonary arteries. Impact of data acquisition and data processing uncertainties were investigated by (1) probabilistic 4D blood flow tracking at varying noise levels, (2) down-sampling to mimic acquisition at different spatial resolutions, (3) pathline calculation with and without eddy current correction, and (4) varied segmentation of the Fontan geometry to mimic analysis errors.

RESULTS

Averaged among the cohort, uncertainties accounted for flow distribution errors from noise ≤3.2%, low spatial resolution ≤2.3% to 3.8%, eddy currents ≤6.4%, and inaccurate segmentation ≤3.9% to 9.1% (dilation and erosion, respectively). In a worst-case scenario (maximum additive errors for all 4 sources), flow distribution errors were as high as 22.5%.

CONCLUSION

Inaccuracies related to postprocessing (segmentation, eddy currents) resulted in the largest potential error (≤15.5% combined) whereas errors related to data acquisition (noise, low spatial resolution) had a lower impact (≤5.5%-7.0% combined). Whereas it is unlikely that these errors will be additive or affect the identification of severe asymmetry, these results illustrate the importance of eddy current correction and accurate segmentation to minimize Fontan flow distribution errors.

In Vitro Validation of 4D Flow MRI for Local Pulse Wave Velocity Estimation

PURPOSE

Arterial stiffness has predictive value for cardiovascular disease (CVD). Local artery stiffness can provide insight on CVD pathology and may be useful for diagnosis and prognosis. However, current methods are invasive, require real-time expertise for measurement, or are limited by arterial region. 4D Flow MRI can non-invasively measure local stiffness by estimating local pulse wave velocity (PWV). This technique can be applied to vascular regions, previously accessible only by invasive stiffness measurement methods. MRI PWV data can also be analyzed post-exam. However, 4D Flow MRI requires validation before it is used in vivo to measure local PWV.

METHODS

PWV, calculated from 4D Flow MRI and a benchtop experiment, was compared with petersons elastic modulus (PEM) of in vitro models. PEM was calculated using high-speed camera images and pressure transducers. Three transit-time algorithms were analyzed for PWV measurement accuracy and precision.

RESULTS

PWV from 4D Flow MRI and reference benchtop experiments show strong correlation with PEM (R² = 0.99). The cross correlation transit-time algorithm showed the lowest percent difference between 4D Flow MRI and benchtop experiments (4-7%), and the point to point of 50% upstroke algorithm had the highest transit-time vs. distance data average R² (0.845).

CONCLUSION

4D Flow MRI is a feasible method for estimating local PWV in simple in vitro models and is a viable tool for clinical analysis. In addition, choice in transit-time algorithm depends on flow waveform shape and arterial region. This study strengthens the validity of 4D Flow MRI local PWV measurement in simple models. However, this technique requires validation in more complex models before it is used in vivo.

Left ventricular blood flow kinetic energy after myocardial infarction - insights from 4D flow cardiovascular magnetic resonance

BACKGROUND

Myocardial infarction (MI) leads to complex changes in left ventricular (LV) haemodynamics that are linked to clinical outcomes. We hypothesize that LV blood flow kinetic energy (KE) is altered in MI and is associated with LV function and infarct characteristics. This study aimed to investigate the intra-cavity LV blood flow KE in controls and MI patients, using cardiovascular magnetic resonance (CMR) four-dimensional (4D) flow assessment.

METHODS

Forty-eight patients with MI (acute-22; chronic-26) and 20 age/gender-matched healthy controls underwent CMR which included cines and whole-heart 4D flow. Patients also received late gadolinium enhancement imaging for infarct assessment. LV blood flow KE parameters were indexed to LV end-diastolic volume and include: averaged LV, minimal, systolic, diastolic, peak E-wave and peak A-wave KEiEDV. In addition, we investigated the in-plane proportion of LV KE (%) and the time difference (TD) to peak E-wave KE propagation from base to mid-ventricle was computed. Association of LV blood flow KE parameters to LV function and infarct size were investigated in all groups.

RESULTS

LV KEiEDV was higher in controls than in MI patients (8.5 ± 3 μJ/ml versus 6.5 ± 3 μJ/ml, P = 0.02). Additionally, systolic, minimal and diastolic peak E-wave KEiEDV were lower in MI (P < 0.05). In logistic-regression analysis, systolic KEiEDV (Beta = - 0.24, P < 0.01) demonstrated the strongest association with the presence of MI. In multiple-regression analysis, infarct size was most strongly associated with in-plane KE (r = 0.5, Beta = 1.1, P < 0.01). In patients with preserved LV ejection fraction (EF), minimal and in-plane KEiEDV were reduced (P < 0.05) and time difference to peak E-wave KE propagation during diastole increased (P < 0.05) when compared to controls with normal EF.

CONCLUSIONS

Reduction in LV systolic function results in reduction in systolic flow KEiEDV. Infarct size is independently associated with the proportion of in-plane LV KE. Degree of LV impairment is associated with TD of peak E-wave KE. In patient with preserved EF post MI, LV blood flow KE mapping demonstrated significant changes in the in-plane KE, the minimal KEiEDV and the TD. These three blood flow KE parameters may offer novel methods to identify and describe this patient population.

Age-related changes of shape and flow dynamics in healthy adult aortas: A 4D flow MRI study

BACKGROUND

Abnormal flow dynamics play an early and causative role in pathologic changes of the ascending aorta.

PURPOSE

To identify: 1) the changes in flow, shape, and size that occur in the ascending aorta with normal human ageing and 2) the influence of these factors on aortic flow dynamics.

STUDY TYPE

Retrospective.

SUBJECTS

In all, 247 subjects (age range 19-86 years, mean 49 ± 17.7, 169 males) free of aortic or aortic valve pathology were included in this study. Subjects were stratified by youngest (18-33 years; n = 64), highest (>60 years, n = 67), and the middle two quartiles (34-60 years, n = 116).

FIELD STRENGTH/SEQUENCE

Subjects underwent a cardiac MRI (3T) exam including 4D-flow MRI of the aorta.

ASSESSMENT

Aortic curvature, arch shape, ascending aortic angle, ascending aortic diameter, and the stroke volume normalized by the aortic volume (nSV) were measured. Velocity, vorticity, and helicity were quantified across the thoracic aorta.

STATISTICAL TESTS

Univariate and multivariate regressions were used to quantify continuous relationships between variables.

RESULTS

Aortic diameter, ascending aortic angle, shape, and curvature all increased across age while nSV decreased (all P < 0.0001). Systolic vorticity in the mid arch decreased by 50% across the age range (P < 0.0001), while peak helicity decreased by 80% (P < 0.0001). Curvature tightly governs optimal flow in the youngest quartile, with an effect size 1.5 to 4 times larger than other parameters in the descending aorta, but had a minimal influence with advancing age. In the upper quartile of age, flow dynamics were almost completely determined by nSV, exerting an effect size on velocity and vorticity >10 times that of diameter and other shape factors.

DATA CONCLUSION

Aortic shape influences flow dynamics in younger subjects. Flow conditions become increasingly disturbed with advancing age, and in these conditions nSV has a more dominant effect on flow patterns than shape factors.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:90-100.

Hemodynamic vascular biomarkers for initiation of paraclinoid internal carotid artery aneurysms using patient-specific computational fluid dynamic simulation based on magnetic resonance imaging

PURPOSE

We performed computational fluid dynamics (CFD) for patients with and without paraclinoid internal carotid artery (ICA) aneurysms to evaluate the distribution of vascular biomarkers at the aneurysm initiation sites of the paraclinoid ICA.

METHODS

This study included 35 patients who were followed up for aneurysms using 3D time of flight (TOF) magnetic resonance angiography (MRA) and 3D cine phase-contrast MR imaging. Fifteen affected ICAs were included in group A with the 15 unaffected contralateral ICAs in group B. Thirty-three out of 40 paraclinoid ICAs free of aneurysms and arteriosclerotic lesions were included in group C. We deleted the aneurysms in group A based on the 3D TOF MRA dataset. We performed CFD based on MR data set and obtained wall shear stress (WSS), its derivatives, and streamlines. We qualitatively evaluated their distributions at and near the intracranial aneurysm initiation site among three groups. We also calculated and compared the normalized highest (nh-) WSS and nh-spatial WSS gradient (SWSSG) around the paraclinoid ICA among three groups.

RESULTS

High WSS and SWSSG distribution were observed at and near the aneurysm initiation site in group A. High WSS and SWSSG were also observed at similar locations in group B and group C. However, nh-WSS and nh-SWSSG were significantly higher in group A than in group C, and nh-SWSSG was significantly higher in group A than in group B.

CONCLUSION

Our findings indicated that nh-WSS and nh-SWSSG were good biomarkers for aneurysm initiation in the paraclinoid ICA.

Segmental biventricular analysis of myocardial function using high temporal and spatial resolution tissue phase mapping

OBJECTIVE

Myocardial dysfunction of the right ventricle (RV) is an important indicator of RV diseases, e.g. RV infarction or pulmonary hypertension. Tissue phase mapping (TPM) has been widely used to determine function of the left ventricle (LV) by analyzing myocardial velocities. The analysis of RV motion is more complicated due to the different geometry and smaller wall thickness. The aim of this work was to adapt and optimize TPM to the demands of the RV.

MATERIALS AND METHODS

TPM measurements were acquired in 25 healthy volunteers using a velocity-encoded phase-contrast sequence and kt-accelerated parallel imaging in combination with optimized navigator strategy and blood saturation. Post processing was extended by a 10-segment RV model and a detailed biventricular analysis of myocardial velocities was performed.

RESULTS

High spatio-temporal resolution (1.0 × 1.0 × 6 mm³, 21.3 ms) and the optimized blood saturation enabled good delineation of the RV and its velocities. Global and segmental velocities, as well as time to peak velocities showed significant differences between the LV and RV. Furthermore, complex timing of the RV could be demonstrated by segmental time to peak analysis.

CONCLUSION

High spatio-temporal resolution TPM enables a detailed biventricular analysis of myocardial motion and might provide a reliable tool for description and detection of diseases affecting left and right ventricular function.

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.

Radiofrequency saturation induced bias in aqueductal cerebrospinal fluid flow quantification obtained using two-dimensional cine phase contrast magnetic resonance imaging

PURPOSE

To explore the extent of bias in cerebrospinal fluid flow estimates due to radiofrequency saturation, and its possible impact on the use of two-dimensional cine phase contrast magnetic resonance imaging in the diagnosis and characterization of normal pressure hydrocephalus in patients.

THEORY AND METHODS

Theoretical signal equations were generated to describe saturation dependence on velocity. An experimental set of phase contrast magnetic resonance imaging scans with two different flip angles was used to show bias in flow estimates in a flow phantom, and in six different healthy volunteers. The cerebral aqueduct was targeted as the flow region of interest.

RESULTS

Data from a constant flow phantom showed a spatial distribution of voxels with significant bias in flow at the periphery of the flow region. The velocity difference (bias) maps of the cerebral aqueduct correlated with the spatial velocity gradients around peak systole and peak diastole, and high correlation with temporal velocity gradients during transition between systole and diastole. The aqueductal stroke volume for θ = 30° were found to be significantly higher than for θ = 10° using a Wilcoxon signed rank test.

CONCLUSION

This work shows the extent of bias in cerebrospinal fluid flow quantification due to radiofrequency saturation effects. This clinical relevance of this error was presented with respect to shunt responsiveness among normal pressure hydrocephalus patients. Magn Reson Med 79:2067-2076, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Comparison of radial 4D Flow-MRI with perivascular ultrasound to quantify blood flow in the abdomen and introduction of a porcine model of pre-hepatic portal hypertension

OBJECTIVES

Objectives of this study were to compare radial time-resolved phase contrast magnetic resonance imaging (4D Flow-MRI) with perivascular ultrasound (pvUS) and to explore a porcine model of acute pre-hepatic portal hypertension (PHTN).

METHODS

Abdominal 4D Flow-MRI and pvUS in portal and splenic vein, hepatic and both renal arteries were performed in 13 pigs of approximately 60 kg. In six pigs, measurements were repeated after partial portal vein (PV) ligature. Inter- and intra-reader comparisons and statistical analysis including Bland-Altman (BA) comparison, paired Student's t tests and linear regression were performed.

RESULTS

PvUS and 4D Flow-MRI measurements agreed well; flow before partial PV ligature was 322 ± 30 ml/min in pvUS and 297 ± 27 ml/min in MRI (p = 0.294), and average BA difference was 25 ml/min [-322; 372]. Inter- and intra-reader results differed very little, revealed excellent correlation (R ² = 0.98 and 0.99, respectively) and resulted in BA differences of -5 ml/min [-161; 150] and -2 ml/min [-28; 25], respectively. After PV ligature, PV flow decreased from 356 ± 50 to 298 ± 61 ml/min (p = 0.02), and hepatic arterial flow increased from 277 ± 36 to 331 ± 65 ml/min (p = n.s.).

CONCLUSION

The successful in vivo comparison of radial 4D Flow-MRI to perivascular ultrasound revealed good agreement of abdominal blood flow although with considerable spread of results. A model of pre-hepatic PHTN was successfully introduced and acute responses monitored.

KEY POINTS

• Radial 4D Flow-MRI in the abdomen was successfully compared to perivascular ultrasound. • Inter- and intra-reader testing demonstrated excellent reproducibility of upper abdominal 4D Flow-MRI. • A porcine model of acute pre-hepatic portal hypertension was successfully introduced. • 4D Flow-MRI successfully monitored acute changes in a model of portal hypertension.

Left Atrial and Left Atrial Appendage 4D Blood Flow Dynamics in Atrial Fibrillation

BACKGROUND

Atrial 4D flow magnetic resonance imaging was used for the characterization of left atrial (LA) and left atrial appendage (LAA) flow dynamics in patients with atrial fibrillation (AF).

METHODS AND RESULTS

4D flow magnetic resonance imaging measured in vivo 3D blood flow velocities in 60 AF patients and 15 controls. Anatomic maps of LA and LAA stasis and velocity were calculated to quantify atrial peak velocity, mean velocity, and stasis (velocities <0.1 m/s). In a substudy with 30 AF patients, 4D flow metrics were compared with Doppler transesophageal echocardiography. For all 15 controls, LAA mean and peak velocities were consistently lower (by 21%/12%; P<0.001) while LAA stasis was higher (by 58%; P<0.001) compared with the LA. In contrast, lower LAA velocity and increased LAA stasis were only found in a fraction (38 of 60) of AF patients. In AF patients, increased CHA2DS2-VASc score was associated with significantly (P<0.043) reduced LA velocities and elevated stasis. There was a heterogeneous expression of atrial flow dynamics, and 25% to 68% of AF patients demonstrated flow in the normal range: 25%/68% for LA/LAA stasis and 38%/60% for LA/LAA peak velocities. Transesophageal echocardiography velocities modestly but significantly (P<0.05) correlated with 4D flow-based LA velocities (r=0.41) and stasis (r=-0.39).

CONCLUSIONS

AF resulted in overall impaired but individually variable flow dynamics in both the LA and LAA. AF patients demonstrated atrial flow in the normal range, despite elevated CHA2DS2-VASc score.

Left ventricular flow propagation velocity measurement: Is it cast in stone?

This study aims to investigate the measurement of left ventricular flow propagation velocity, V _p, using phase contrast magnetic resonance imaging and to assess the discrepancies resulting from inflow jet direction and individual left ventricular size. Three V _p measuring techniques, namely non-adaptive (NA), adaptive positions (AP) and adaptive vectors (AV) method, were suggested and compared. We performed the comparison on nine healthy volunteers and nine post-infarct patients at four measurement positions, respectively, at one-third, one-half, two-thirds and the conventional 4 cm distances from the mitral valve leaflet into the left ventricle. We found that the V _p measurement was affected by both the inflow jet direction and measurement positions. Both NA and AP methods overestimated V _p, especially in dilated left ventricles, while the AV method showed the strongest correlation with the isovolumic relaxation myocardial strain rate (r = 0.53, p < 0.05). Using the AV method, notable difference in mean V _p was also observed between healthy volunteers and post-infarct patients at positions of: one-half (81 ± 31 vs. 58 ± 25 cm/s), two-thirds (89 ± 32 vs. 45 ± 15 cm/s) and 4 cm (98 ± 23 vs. 47 ± 13 cm/s) distances. The use of AV method and measurement position at one-half distance was found to be the most suitable method for assessing diastolic dysfunction given varying left ventricular sizes and inflow jet directions.

Superior Abdominal 4D Flow MRI Data Consistency with Adjusted Preprocessing Workflow and Noncontrast Acquisitions

RATIONALE AND OBJECTIVES

To assess the impact of an alternative preprocessing workflow on noncontrast- and contrast-enhanced abdominal four-dimensional flow magnetic resonance imaging (4D flow MRI) data consistency.

MATERIALS AND METHODS

Twenty patients with cirrhosis and portal hypertension (5 women; 53 ± 10 years old) underwent 4D flow MRI at 3.0T before and after administration of 0.03 mmol/kg of gadofosveset trisodium with velocity sensitivities of 100 and 50 cm/s for arterial and venous flow quantifications, respectively. 4D flow MRI data were preprocessed using the conventional workflow (workflow 1), applying noise filters prior to eddy current correction, and an alternative workflow (workflow 2), first correcting for eddy currents and using noise filtering only if needed for anti-aliasing. Vessel segmentation quality was ranked by independent reviewers and compared via Wilcoxon signed-rank tests. Flow quantification and conservation of mass at two portal and one arterial branch points were compared via paired t tests.

RESULTS

Segmentation quality was significantly higher for workflow 2 (P < 0.05) with excellent interobserver agreement (κ = 0.92). Workflow 2 resulted in larger flow values (P < 0.05) with improved conservation of mass (7.3 ± 6.1% vs. 27.7 ± 25.0%, P < 0.001 [portal]; 10.7 ± 9.0% vs. 21.7 ± 21.6%, P = 0.02 [arterial]). Peak velocities and abdominal aortic flow were similar (P > 0.05). Noncontrast acquisitions yielded significantly smaller portal flow values (P < 0.05) with improved conservation of mass (5.8 ± 4.7% vs. 8.7 ± 6.9%, P = 0.05 [portal]; 6.2 ± 4.5% vs. 13.7 ± 10.2%, P = 0.03 [arterial]).

CONCLUSIONS

Superior abdominal 4D flow MRI data consistency was obtained by applying eddy current correction before any other data manipulation, using noise masking and velocity anti-aliasing cautiously, and using noncontrast acquisitions.

3D phase contrast MRI: Partial volume correction for robust blood flow quantification in small intracranial vessels

PURPOSE

Recent advances in 3D-PCMRI (phase contrast MRI) sequences allow for measuring the complex hemodynamics in cerebral arteries. However, the small size of these vessels vs spatial resolution can lead to non-negligible partial volume artifacts, which must be taken into account when computing blood flow rates. For this purpose, we combined the velocity information provided by 3D-PCMRI with vessel geometry measured with 3DTOF (time of flight MRI) or 3DRA (3D rotational angiography) to correct the partial volume effects in flow rate assessments.

METHODS

The proposed methodology was first tested in vitro on cylindrical and patient specific vessels subject to fully controlled pulsatile flows. Both 2D- and 3D-PCMRI measurements using various spatial resolutions ranging from 20 to 1.3 voxels per vessel diameter were analyzed and compared with flowmeter baseline. Second, 3DTOF, 2D- and 3D-PCMRI measurements were performed in vivo on 35 patients harboring internal carotid artery (ICA) aneurysms indicated for endovascular treatments requiring 3DRA imaging.

RESULTS

The in vitro 2D- and 3D-PCMRI mean flow rates assessed with partial volume correction showed very low sensitivity to the acquisition resolution above ≈2 voxels per vessel diameter while uncorrected flow rates deviated critically when decreasing the spatial resolution. 3D-PCMRI flow rates measured in vivo in ICA agreed very well with 2D-PCMRI data and a good flow conservation was observed at the C7 bifurcation. Globally, partial volume correction led to 10-15% lower flow rates than uncorrected values as those reported in most of the published studies on intracranial flows.

CONCLUSION

Partial volume correction may improve the accuracy of PCMRI flow rate measurements especially in small vessels such as intracranial arteries. Magn Reson Med 79:129-140, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Phase-contrast MR flow imaging: A tool to determine hepatic hemodynamics in rats with a healthy, fibrotic, or cirrhotic liver

PURPOSE

To test a magnetic resonance (MR) scanning protocol as a noninvasive tool to determine hepatic hemodynamics and to assess the degree of liver fibrosis in an animal model of liver fibrosis and cirrhosis.

MATERIALS AND METHODS

Fifty-four male Wistar rats were studied. Thirty-nine received thioacetamide (TAA) in their drinking water for either 12 or 16 weeks. MR measurements were performed using flow-sensitive 2D phase-contrast MRI and a 9.4T preclinical scanner. The following hemodynamic parameters were investigated: portal cross-sectional area, mean portal flow velocity, and portal and aortic flow volume rate. Therefore, rats (n = 46) were divided into three groups: CON (control, n = 13), FIB (fibrosis, n = 25), and CIR (cirrhosis, n = 8). Furthermore, the degree of liver fibrosis was assessed by a self-established MR score and verified by a standardized histological score (n = 48).

RESULTS

Portal and aortic flow parameters could be reliably detected. A significant decrease in portal flow velocity was found in FIB (FIB vs. CON: -21%, P = 0.006 and CIR vs. CON: -17%, P = 0.105) and in portal flow volume rate in FIB and CIR (FIB vs. CON: -20%, P = 0.009 and CIR vs. CON: -25%, P = 0.024). If the histological score is taken as standard, the self-established MR score enabled discrimination between healthy and diseased livers (sensitivity to identify diseased livers: 89% and specificity to identify healthy livers: 100%).

CONCLUSION

This MR scanning protocol presents a noninvasive tool to determine hepatic hemodynamics in healthy and diseased rats.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2017;46:1526-1534.

Flow mixing during peripheral veno-arterial extra corporeal membrane oxygenation - A simulation study

Peripheral veno-arterial extra-corporeal membrane oxygenation (ECMO) is an artificial circulation that supports patients with severe cardiac and respiratory failure. Differential hypoxia during ECMO support has been reported, and it has been suggested that it is due to the mixing of well-perfused retrograde ECMO flow and poorly-perfused antegrade left ventricle (LV) flow in the aorta. This study aims to quantify the relationship between ECMO support level and location of the mixing zone (MZ) of the ECMO and LV flows. Steady-state and transient computational fluid dynamics (CFD) simulations were performed using a patient-specific geometrical model of the aorta. A range of ECMO support levels (from 5% to 95% of total cardiac output) were evaluated. For ECMO support levels above 70%, the MZ was located in the aortic arch, resulting in perfusion of the arch branches with poorly perfused LV flow. The MZ location was stable over the cardiac cycle for high ECMO flows (>70%), but moved 5cm between systole and diastole for ECMO support level of 60%. This CFD approach has potential to improve individual patient care and ECMO design.

Real-time phase-contrast flow cardiovascular magnetic resonance with low-rank modeling and parallel imaging

BACKGROUND

Conventional phase-contrast cardiovascular magnetic resonance (PC-CMR) employs cine-based acquisitions to assess blood flow condition, in which electro-cardiogram (ECG) gating and respiration control are generally required. This often results in lower acquisition efficiency, and limited utility in the presence of cardiovascular pathology (e.g., cardiac arrhythmia). Real-time PC-CMR, without ECG gating and respiration control, is a promising alternative that could overcome limitations of the conventional approach. But real-time PC-CMR involves image reconstruction from highly undersampled (k, t)-space data, which is very challenging. In this study, we present a novel model-based imaging method to enable high-resolution real-time PC-CMR with sparse sampling.

METHODS

The proposed method captures spatiotemporal correlation among flow-compensated and flow-encoded image sequences with a novel low-rank model. The image reconstruction problem is then formulated as a low-rank matrix recovery problem. With proper temporal subspace modeling, it results in a convex optimization formulation. We further integrate this formulation with the SENSE-based parallel imaging model to handle multichannel acquisitions. The performance of the proposed method was systematically evaluated in 2D real-time PC-CMR with flow phantom experiments and in vivo experiments (with healthy subjects). Additionally, we performed a feasibility study of the proposed method on patients with cardiac arrhythmia.

RESULTS

The proposed method achieves a spatial resolution of 1.8 mm and a temporal resolution of 18 ms for 2D real-time PC-CMR with one directional flow encoding. For the flow phantom experiments, both regular and irregular flow patterns were accurately captured. For the in vivo experiments with healthy subjects, flow dynamics obtained from the proposed method correlated well with those from the cine-based acquisitions. For the experiments with the arrhythmic patients, the proposed method demonstrated excellent capability of resolving the beat-by-beat flow variations, which cannot be obtained from the conventional cine-based method.

CONCLUSION

The proposed method enables high-resolution real-time PC-CMR at 2D without ECG gating and respiration control. It accurately resolves beat-by-beat flow variations, which holds great promise for studying patients with irregular heartbeats.

Accelerated dual-venc 4D flow MRI for neurovascular applications

PURPOSE

To improve velocity-to-noise ratio (VNR) and dynamic velocity range of 4D flow magnetic resonance imaging (MRI) by using dual-velocity encoding (dual-venc) with k-t generalized autocalibrating partially parallel acquisition (GRAPPA) acceleration.

MATERIALS AND METHODS

A dual-venc 4D flow MRI sequence with k-t GRAPPA acceleration was developed using a shared reference scan followed by three-directional low- and high-venc scans (repetition time / echo time / flip angle = 6.1 msec / 3.4 msec / 15°, temporal/spatial resolution = 43.0 msec/1.2 × 1.2 × 1.2 mm³ ). The high-venc data were used to correct for aliasing in the low-venc data, resulting in a single dataset with the favorable VNR of the low-venc but without velocity aliasing. The sequence was validated with a 3T MRI scanner in phantom experiments and applied in 16 volunteers to investigate its feasibility for assessing intracranial hemodynamics (net flow and peak velocity) at the major intracranial vessels. In addition, image quality and image noise were assessed in the in vivo acquisitions.

RESULTS

All 4D flow MRI scans were acquired successfully with an acquisition time of 20 ± 4 minutes. The shared reference scan reduced the total acquisition time by 12.5% compared to two separate scans. Phantom experiments showed 51.4% reduced noise for dual-venc compared to high-venc and an excellent agreement of velocities (ρ = 0.8, P < 0.001). The volunteer data showed decreased noise in dual-venc data (54.6% lower) compared to high-venc, and improved image quality, as graded by two observers: fewer artifacts (P < 0.0001), improved vessel conspicuity (P < 0.0001), and reduced noise (P < 0.0001).

CONCLUSION

Dual-venc 4D flow MRI exhibits the superior VNR of the low-venc acquisition and reliably incorporates low- and high-velocity fields simultaneously. In vitro and in vivo data demonstrate improved flow visualization, image quality, and image noise.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:102-114.

Influence of Spatial Resolution in Three-dimensional Cine Phase Contrast Magnetic Resonance Imaging on the Accuracy of Hemodynamic Analysis

Introduction:

We aim to elucidate the effect of spatial resolution of three-dimensional cine phase contrast magnetic resonance (3D cine PC MR) imaging on the accuracy of the blood flow analysis, and examine the optimal setting for spatial resolution using flow phantoms.

Materials and Methods:

The flow phantom has five types of acrylic pipes that represent human blood vessels (inner diameters: 15, 12, 9, 6, and 3 mm). The pipes were fixed with 1% agarose containing 0.025 mol/L gadolinium contrast agent. A blood-mimicking fluid with human blood property values was circulated through the pipes at a steady flow. Magnetic resonance (MR) images (three-directional phase images with speed information and magnitude images for information of shape) were acquired using the 3-Tesla MR system and receiving coil. Temporal changes in spatially-averaged velocity and maximum velocity were calculated using hemodynamic analysis software. We calculated the error rates of the flow velocities based on the volume flow rates measured with a flowmeter and examined measurement accuracy.

Results:

When the acrylic pipe was the size of the thoracicoabdominal or cervical artery and the ratio of pixel size for the pipe was set at 30% or lower, spatially-averaged velocity measurements were highly accurate. When the pixel size ratio was set at 10% or lower, maximum velocity could be measured with high accuracy. It was difficult to accurately measure maximum velocity of the 3-mm pipe, which was the size of an intracranial major artery, but the error for spatially-averaged velocity was 20% or less.

Conclusions:

Flow velocity measurement accuracy of 3D cine PC MR imaging for pipes with inner sizes equivalent to vessels in the cervical and thoracicoabdominal arteries is good. The flow velocity accuracy for the pipe with a 3-mm-diameter that is equivalent to major intracranial arteries is poor for maximum velocity, but it is relatively good for spatially-averaged velocity.

Surgical planning for living donor liver transplant using 4D flow MRI, computational fluid dynamics and in vitro experiments

This study used magnetic resonance imaging (MRI), computational fluid dynamics (CFD) modeling, and in vitro experiments to predict patient-specific alterations in hepatic hemodynamics in response to partial hepatectomy in living liver donors. 4D Flow MRI was performed on three donors before and after hepatectomy and models of the portal venous system were created. Virtual surgery was performed to simulate (1) surgical resection and (2) post-surgery vessel dilation. CFD simulations were conducted using in vivo flow data for boundary conditions. CFD results showed good agreement with in vivo data, and in vitro experimental values agreed well with imaging and simulation results. The post-surgery models predicted an increase in all measured hemodynamic parameters, and the dilated virtual surgery model predicted post-surgery conditions better than the model that only simulated resection. The methods used in this study have potential significant value for the surgical planning process for the liver and other vascular territories.

Qualitative grading of aortic regurgitation: a pilot study comparing CMR 4D flow and echocardiography

Over the past 10 years there has been intense research in the development of volumetric visualization of intracardiac flow by cardiac magnetic resonance (CMR). This volumetric time resolved technique called CMR 4D flow imaging has several advantages over standard CMR. It offers anatomical, functional and flow information in a single free-breathing, ten-minute acquisition. However, the data obtained is large and its processing requires dedicated software. We evaluated a cloud-based application package that combines volumetric data correction and visualization of CMR 4D flow data, and assessed its accuracy for the detection and grading of aortic valve regurgitation using transthoracic echocardiography as reference. Between June 2014 and January 2015, patients planned for clinical CMR were consecutively approached to undergo the supplementary CMR 4D flow acquisition. Fifty four patients (median age 39 years, 32 males) were included. Detection and grading of the aortic valve regurgitation using CMR 4D flow imaging were evaluated against transthoracic echocardiography. The agreement between 4D flow CMR and transthoracic echocardiography for grading of aortic valve regurgitation was good (κ = 0.73). To identify relevant, more than mild aortic valve regurgitation, CMR 4D flow imaging had a sensitivity of 100 % and specificity of 98 %. Aortic regurgitation can be well visualized, in a similar manner as transthoracic echocardiography, when using CMR 4D flow imaging.

Unravelling cardiovascular disease using four dimensional flow cardiovascular magnetic resonance

Knowledge of normal and abnormal flow patterns in the human cardiovascular system increases our understanding of normal physiology and may help unravel the complex pathophysiological mechanisms leading to cardiovascular disease. Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has emerged as a suitable technique that enables visualization of in vivo blood flow patterns and quantification of parameters that could potentially be of prognostic value in the disease process. In this review, current image processing tools that are used for comprehensive visualization and quantification of blood flow and energy distribution in the heart and great vessels will be discussed. Also, imaging biomarkers extracted from 4D flow CMR will be reviewed that have been shown to distinguish between normal and abnormal flow patterns. Furthermore, current applications of 4D flow CMR in the heart and great vessels will be discussed, showing its potential as an additional diagnostic modality which could aid in disease management and timing of surgical intervention.

Accuracy of 4D Flow Measurement of Cerebrospinal Fluid Dynamics in the Cervical Spine: An In Vitro Verification Against Numerical Simulation

Abnormal alterations in cerebrospinal fluid (CSF) flow are thought to play an important role in pathophysiology of various craniospinal disorders such as hydrocephalus and Chiari malformation. Three directional phase contrast MRI (4D Flow) has been proposed as one method for quantification of the CSF dynamics in healthy and disease states, but prior to further implementation of this technique, its accuracy in measuring CSF velocity magnitude and distribution must be evaluated. In this study, an MR-compatible experimental platform was developed based on an anatomically detailed 3D printed model of the cervical subarachnoid space and subject specific flow boundary conditions. Accuracy of 4D Flow measurements was assessed by comparison of CSF velocities obtained within the in vitro model with the numerically predicted velocities calculated from a spatially averaged computational fluid dynamics (CFD) model based on the same geometry and flow boundary conditions. Good agreement was observed between CFD and 4D Flow in terms of spatial distribution and peak magnitude of through-plane velocities with an average difference of 7.5 and 10.6% for peak systolic and diastolic velocities, respectively. Regression analysis showed lower accuracy of 4D Flow measurement at the timeframes corresponding to low CSF flow rate and poor correlation between CFD and 4D Flow in-plane velocities.

Left Atrial 4-Dimensional Flow Magnetic Resonance Imaging: Stasis and Velocity Mapping in Patients With Atrial Fibrillation

OBJECTIVES

Left atrial (LA) 4-dimensional flow magnetic resonance imaging (MRI) was used to derive anatomic maps of LA stasis, peak velocity, and time-to-peak (TTP) velocity in patients with atrial fibrillation (AF) and to identify relationships between LA flow with LA volume and patient characteristics.

MATERIALS AND METHODS

Four-dimensional flow MRI for the in vivo assessment of time-resolved 3-dimensional LA blood flow velocities was performed in 111 subjects: 42 patients with a history of AF and in sinus rhythm (AF-sinus), 39 patients with persistent AF (AF-afib), 10 young healthy volunteers (HVs), and 20 age-appropriate controls (CTRL). Data analysis included the 3-dimensional segmentation of the LA and the calculation of LA stasis, peak velocity, and TTP maps. Regional LA flow dynamics were quantified by calculating mean stasis, peak velocity, and TTP in the LA center region and the region adjacent to the LA wall.

RESULTS

A sensitivity analysis identified thresholds for global LA stasis (<0.1 m/s) and peak velocity (top 5% LA velocities), which detected significant differences between AF patients and controls for global LA stasis (HV, 25% ± 5%; CTRL, 29% ± 10%; AF-sinus, 41% ± 13%; AF-afib, 52% ± 17%) and peak velocity (HV, 0.43 ± 0.02 m/s; CTRL, 0.37 ± 0.04 m/s; AF-sinus, 0.33 ± 0.05 m/s; AF-afib, 0.30 ± 0.05 m/s). Regional analysis revealed significantly increased stasis at both LA center and wall for AF patients compared with age-appropriate controls (29%-84% difference, P < 0.006) and for AF-afib versus AF-sinus patients (22%-30% difference, P < 0.004). In addition, stasis close to the LA wall was significantly elevated (P < 0.001) compared with the LA center for all subject groups. Multiple regressions revealed significant (RAdj = 0.45-0.50, P < 0.001) relationships between impaired global LA flow (reduced velocity and increased stasis) with age (|β| = 0.27-0.50, P < 0.002) and LA volume (|β| = 0.26-0.50, P < 0.003).

CONCLUSIONS

Atrial 4-dimensional flow MRI detected changes in global and regional LA flow dynamics associated with AF, age, and LA volume. Longitudinal studies are needed to test the diagnostic value of LA flow metrics as potential risk factors for thromboembolic events.

Evaluation of blood flow distribution asymmetry and vascular geometry in patients with Fontan circulation using 4-D flow MRI

BACKGROUND

Asymmetrical caval to pulmonary blood flow is suspected to cause complications in patients with Fontan circulation. The aim of this study was to test the feasibility of 4-D flow MRI for characterizing the relationship between 3-D blood flow distribution and vascular geometry.

OBJECTIVE

We hypothesized that both flow distribution and geometry can be calculated with low interobserver variability and will detect a direct relationship between flow distribution and Fontan geometry.

MATERIALS AND METHODS

Four-dimensional flow MRI was acquired in 10 Fontan patients (age: 16 ± 4 years [mean ± standard deviation], range: 9-21 years). The Fontan connection was isolated by 3-D segmentation to evaluate flow distribution from the inferior vena cava (IVC) and superior vena cava (SVC) to the left and right pulmonary arteries (LPA, RPA) and to characterize geometry (cross-sectional area, caval offset, vessel angle).

RESULTS

Flow distribution results indicated SVC flow tended toward the RPA while IVC flow was more evenly distributed (SVC to RPA: 78% ± 28 [9-100], IVC to LPA: 54% ± 28 [4-98]). There was a significant relationship between pulmonary artery cross-sectional area and flow distribution (IVC to RPA: R(2)=0.50, P=0.02; SVC to LPA: R(2)=0.81, P=0.0004). Good agreement was found between observers and for flow distribution when compared to net flow values.

CONCLUSION

Four-dimensional flow MRI was able to detect relationships between flow distribution and vessel geometry. Future studies are warranted to investigate the potential of patient specific hemodynamic analysis to improve diagnostic capability.

Hemodynamic evaluation in patients with transposition of the great arteries after the arterial switch operation: 4D flow and 2D phase contrast cardiovascular magnetic resonance compared with Doppler echocardiography

BACKGROUND

Peak velocity measurements are used to evaluate the significance of stenosis in patients with transposition of the great arteries after the arterial switch operation (TGA after ASO). 4D flow cardiovascular magnetic resonance (CMR) provides 3-directional velocity encoding and full volumetric coverage of the great arteries and may thus improve the hemodynamic evaluation in these patients. The aim of this study was to compare peak velocities measured by 4D flow CMR with 2D phase contrast (PC) CMR and the gold standard Doppler echocardiography (echo) in patients with TGA after ASO.

METHODS

Nineteen patients (mean age 13 ± 9 years, range 1-25 years) with TGA after ASO who underwent 2D PC CMR and 4D flow CMR were included in this study. Peak velocities were measured with 4D flow CMR in the aorta and pulmonary arteries and compared to peak velocities measured with 2D PC CMR and Doppler echo. 2D PC CMR data were available in the ascending aorta, main, right and left pulmonary arteries (AAO/MPA/RPA/LPA) for 19/18/17/17 scans, respectively, and Doppler echo data were available for 13/9/6/6 scans, respectively. Peak velocities were measured with: 1) a single cross section for 2D PC CMR, 2) velocity maximum intensity projections (MIPs) for 4D flow CMR and 3) Doppler echo.

RESULTS

Significantly higher peak velocities were found with 4D flow CMR than 2D PC CMR in the AAO (p = 0.003), MPA (p = 0.002) and RPA (p = 0.005) but not in the LPA (p = 0.200). No difference in peak velocity was found between 4D flow CMR and Doppler echo (p > 0.46) or 2D PC CMR and echo (p > 0.11) for all analyzed vessel segments.

CONCLUSIONS

4D flow CMR evaluation of patients with TGA after ASO detected higher peak velocities than 2D PC CMR, indicating the potential of 4D flow CMR to provide improved stenosis assessment in these patients.

Altered aortic shape in bicuspid aortic valve relatives influences blood flow patterns

AIMS

Bicuspid aortic valve (BAV) is known to exhibit familial inheritance and is associated with aortopathy and altered aortic haemodynamics. However, it remains unclear whether BAV-related aortopathy can be inherited independently of valve morphology.

METHODS AND RESULTS

Four-dimensional flow magnetic resonance imaging for the in vivo assessment of thoracic aortic 3D blood flow was performed in 24 BAV relatives with trileaflet aortic valves (age = 40 ± 14 years) and 15 healthy controls (age = 37 ± 10 years). Data analysis included aortic dimensions, shape (round/gothic/cubic), and 3D blood flow characteristics (semi-quantitative vortex/helix grading and peak velocities). Cubic and gothic aortic shapes were markedly more prevalent in BAV relatives compared with controls (38 vs. 7%). Ascending aorta (AAo) vortex flow in BAV relatives was significantly increased compared with controls (grading = 1.5 ± 1.0 vs. 0.6 ± 0.9, P = 0.015). Aortic haemodynamics were influenced by aortic shape: peak velocities were reduced for gothic aortas vs. round aortas (P = 0.003); vortex flow was increased for cubic aortas in the AAo (P < 0.001) and aortic arch (P = 0.004); vortex and helix flows were elevated for gothic aortas in the AAo and descending aorta (P = 0.003, P = 0.029). Logistic regression demonstrated significant associations of shape with severity of vortex flow in AAo (P < 0.001) and aortic arch (P = 0.016) in BAV relatives.

CONCLUSION

BAV relatives expressed altered aortic shape and increased vortex flow despite the absence of valvular disease or aortic dilatation. These data suggest a heritable component of BAV-related aortopathy affecting aortic shape and aberrant blood flow, independent of valve morphology.

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).

Hemodynamic Measurement Using Four-Dimensional Phase-Contrast MRI: Quantification of Hemodynamic Parameters and Clinical Applications

Recent improvements have been made to the use of time-resolved, three-dimensional phase-contrast (PC) magnetic resonance imaging (MRI), which is also named four-dimensional (4D) PC-MRI or 4D flow MRI, in the investigation of spatial and temporal variations in hemodynamic features in cardiovascular blood flow. The present article reviews the principle and analytical procedures of 4D PC-MRI. Various fluid dynamic biomarkers for possible clinical usage are also described, including wall shear stress, turbulent kinetic energy, and relative pressure. Lastly, this article provides an overview of the clinical applications of 4D PC-MRI in various cardiovascular regions.

Advanced Analysis Techniques for Intra-cardiac Flow Evaluation from 4D Flow MRI

PURPOSE OF THE REVIEW

Time-resolved 3D velocity-encoded MR imaging with velocity encoding in three directions (4D Flow) has emerged as a novel MR acquisition technique providing detailed information on flow in the cardiovascular system. In contrast to other clinically available imaging techniques such as echo-Doppler, 4D Flow MRI provides the 3D Flow velocity field within a volumetric region of interest over the cardiac cycle. This work reviews the most recent advances in the development and application of dedicated image analysis techniques for the assessment of intra-cardiac flow features from 4D Flow MRI.

RECENT FINDINGS

Novel image analysis techniques have been developed for extraction of relevant intra-cardiac flow features from 4D Flow MRI, which have been successfully applied in various patient cohorts and volunteer studies. Disturbed flow patterns have been linked with valvular abnormalities and ventricular dysfunction. Recent technical advances have resulted in reduced scan times and improvements in image quality, increasing the potential clinical applicability of 4D Flow MRI.

SUMMARY

4D Flow MRI provides unique capabilities for 3D visualization and quantification of intra-cardiac blood flow. Contemporary knowledge on 4D Flow MRI shows promise for further exploration of the potential use of the technique in research and clinical applications.

4D flow MR imaging of the portal venous system: a feasibility study in children

OBJECTIVES

To determine the feasibility of 4D flow MRI for visualization and quantification of the portal venous haemodynamics in children and young adults.

METHODS

4D flow was performed in 28 paediatric patients (median age, 8.5 years; interquartile range, 5.2-16.5), 15 with non-operated native portal system and 13 with surgically created portal shunt. Image quality assessment for 3D flow visualization and flow pattern analyses was performed. Regional 4D flow peak velocity and net flow were compared with 2D-cine phase contrast MRI (2D-PC MR) in the post-surgical patients.

RESULTS

Mean 3D flow visualization quality score was excellent (mean ± SD, 4.2 ± 0.9) with good inter-rater agreement (κ,0.67). Image quality in children aged >10 years was better than children ≤10 years (p < 0.05). Flow pattern was defined for portal, superior mesenteric, splenic veins and splenic artery in all patients. 4D flow and 2D-PC MR peak velocity and net flow were similar with good correlation (peak velocity: 4D flow 22.2 ± 9.1 cm/s and 2D-PC MR 25.2 ± 11.2 cm/s, p = 0.46; r = 0.92, p < 0.0001; net flow: 4D flow 9.5 ± 7.4 ml/s and 2D-PC MR 10.1 ± 7.3 ml/s, p = 0.65; r = 0.81, p = 0.0007).

CONCLUSIONS

4D flow MRI is feasible and holds promise for the comprehensive 3D visualization and quantification of portal venous flow dynamics in children and young adults.

KEY POINTS

• 4D flow MRI is feasible in children and young adults. • 4D flow MRI has the ability to non-invasively characterize portal haemodynamics. • Image quality of 4D flow MRI is better is older children. • 4D flow MRI can accurately quantify portal flow compared to 2D-cine PC MRI.

Emerging Applications of Abdominal 4D Flow MRI

OBJECTIVE

Comprehensive assessment of abdominal hemodynamics is crucial for many clinical diagnoses but is challenged by a tremendous complexity of anatomy, normal physiology, and a wide variety of pathologic abnormalities. This article introduces 4D flow MRI as a powerful technique for noninvasive assessment of the hemodynamics of abdominal vascular territories.

CONCLUSION

Four-dimensional flow MRI provides clinicians with a more extensive and straightforward approach to evaluate disorders that affect blood flow in the abdomen. This review presents a series of clinical cases to illustrate the utility of 4D flow MRI in the comprehensive assessment of the abdominal circulation.

Improved Semiautomated 4D Flow MRI Analysis in the Aorta in Patients With Congenital Aortic Valve Anomalies Versus Tricuspid Aortic Valves

OBJECTIVE

The aim of this study was to systematically investigate a newly developed semiautomated workflow for the analysis of aortic 4-dimensional flow MRI and its ability to detect hemodynamic differences in patients with congenitally altered aortic valve (bicuspid or quadricuspid valves) compared with tricuspid aortic valves.

METHODS

Four-dimensional flow MRI data were acquired in 20 patients with aortic dilatation (9 tricuspid aortic valves, 11 congenitally altered aortic valves). A semiautomated workflow was evaluated regarding interobserver variability, accuracy of net flow, regurgitant fraction and peak systolic velocity, and the ability to detect differences between cohorts. Results were compared with manual segmentation of vessel contours.

RESULTS

Despite the significantly reduced analysis time, a good interobserver agreement was found for net flow and peak systolic velocity, and a moderate agreement was found for regurgitation. Significant differences in peak velocities in the descending aorta (P = 0.014) could be detected.

CONCLUSIONS

Four-dimensional flow MRI-based semiautomated analysis of aortic hemodynamics can be performed with good reproducibility and accuracy.

Longitudinal Monitoring of Hepatic Blood Flow before and after TIPS by Using 4D-Flow MR Imaging

Purpose To demonstrate the feasibility of four-dimensional (4D)-flow magnetic resonance (MR) imaging for noninvasive longitudinal hemodynamic monitoring of hepatic blood flow before and after transjugular intrahepatic portosystemic shunt (TIPS) placement. Materials and Methods The institutional review board approved this prospective Health Insurance Portability and Accountability Act compliant study with written informed consent. Four-dimensional-flow MR imaging was performed in seven patients with portal hypertension and refractory ascites before and 2 and 12 weeks after TIPS placement by using a time-resolved three-dimensional radial phase-contrast acquisition. Flow and peak velocity measurements were obtained in the superior mesenteric vein (SMV), splenic vein (SV), portal vein (PV), and the TIPS. Flow volumes and peak velocities in each vessel, as well as the ratio of in-stent to PV flow, were compared before and after TIPS placement by using analysis of variance. Results Flow volumes significantly increased in the SMV (0.24 L/min; 95% confidence interval [CI]: 0.07, 0.41), SV (0.31 L/min; 95% CI: 0.07, 0.54), and PV (0.88 L/min; 95% CI: 0.06, 1.70) after TIPS placement (all P < .05), with no significant difference between the first and second post-TIPS placement acquisitions (all P > .11). Ascites resolved in six of seven patients. In those with resolved ascites, the TIPS-to-PV flow ratio was 0.8 ± 0.2 and 0.9 ± 0.2 at the two post-TIPS time points, respectively, while the observed ratios were 4.6 and 4.3 in the patient with refractory ascites at the two post-TIPS time points, respectively. In this patient, 4D-flow MR imaging demonstrated arterio-portal-venous shunting, with draining into the TIPS. Conclusion Four-dimensional-flow MR imaging is feasible for noninvasive longitudinal hemodynamic monitoring of hepatic blood flow before and after TIPS placement. ^© RSNA, 2016 Online supplemental material is available for this article.

On the influence of respiratory motion in radial tissue phase mapping cardiac MRI

PURPOSE

To investigate the impact of respiratory motion on radial tissue phase mapping (TPM) measurements, and to improve image quality and scan efficiency without compromising velocity fidelity by increasing the respiratory acceptance window with and without motion correction.

MATERIALS AND METHODS

A radial golden angle TPM sequence was measured in 10 healthy volunteers in three short axis slices at 3T. Ungated ( CFREE), self-gated with a single acceptance window ( CREF), motion-corrected averaging using all ( CMCall), or selected ( CMC) data reconstructions were compared by means of various image quality measures and resulting velocities.

RESULTS

Using all data ( CFREE) resulted in significantly higher perceived signal-to-noise ratio (SNR) (P < 0.001), but significantly reduced sharpness (P < 0.001) and contrast (P = 0.02), when compared to CREF. Coefficient of variation (CV) and perceived sharpness were not significantly different (P > 0.05). With motion-correction, perceived sharpness could be significantly improved ( CMC: P = 0.002; CMCall: P = 0.002) in comparison to CFREE. Velocity peaks of CFREE were significantly reduced compared to CREF (all peaks: P < 0.001; except the longitudinal "E" peak: P = 0.03). The peak velocities in CMC and CMCall were not significantly different from CREF (all peaks: P > 0.08; except longitudinal "E"/"A" peaks: P > 0.01).

CONCLUSION

Free-breathing reconstruction results in good perceived image sharpness and velocity information with slightly, but significantly, reduced peak velocities. For achieving velocities and image quality comparable to data from a single acceptance window, but higher gating efficiency, selected motion-corrected TPM (CMC) can be applied. J. Magn. Reson. Imaging 2016;44:1218-1228.

Accelerated phase contrast flow imaging with direct complex difference reconstruction

PURPOSE

To propose and evaluate a new model-based reconstruction method for highly accelerated phase-contrast magnetic resonance imaging (PC-MRI) with sparse sampling.

THEORY AND METHODS

This work presents a new constrained reconstruction method based on low-rank and sparsity constraints to accelerate PC-MRI. More specifically, we formulate the image reconstruction problem into separate reconstructions of flow-reference image sequence and complex differences. We then utilize the joint partial separability and sparsity constraints to enable high quality reconstruction from highly undersampled (k,t)-space data. We further integrate the proposed method with ESPIRiT based parallel imaging model to effectively handle multichannel acquisition.

RESULTS

The proposed method was evaluated with in vivo data acquired from both 2D and 3D PC flow imaging experiments, and compared with several state-of-the-art methods. Experimental results demonstrate that the proposed method leads to more accurate velocity reconstruction from highly undersampled (k,t)-space data, and particularly superior capability of capturing the peak velocity of blood flow. In terms of flow visualization, blood flow patterns obtained from the proposed reconstruction also exhibit better agreement with those obtained from the fully sampled reference.

CONCLUSION

The proposed method achieves improved accuracy over several state-of-the-art methods for velocity reconstruction with highly accelerated (k,t)-space data. Magn Reson Med 77:1036-1048, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Advanced flow MRI: emerging techniques and applications

Magnetic resonance imaging (MRI) techniques provide non-invasive and non-ionising methods for the highly accurate anatomical depiction of the heart and vessels throughout the cardiac cycle. In addition, the intrinsic sensitivity of MRI to motion offers the unique ability to acquire spatially registered blood flow simultaneously with the morphological data, within a single measurement. In clinical routine, flow MRI is typically accomplished using methods that resolve two spatial dimensions in individual planes and encode the time-resolved velocity in one principal direction, typically oriented perpendicular to the two-dimensional (2D) section. This review describes recently developed advanced MRI flow techniques, which allow for more comprehensive evaluation of blood flow characteristics, such as real-time flow imaging, 2D multiple-venc phase contrast MRI, four-dimensional (4D) flow MRI, quantification of complex haemodynamic properties, and highly accelerated flow imaging. Emerging techniques and novel applications are explored. In addition, applications of these new techniques for the improved evaluation of cardiovascular (aorta, pulmonary arteries, congenital heart disease, atrial fibrillation, coronary arteries) as well as cerebrovascular disease (intra-cranial arteries and veins) are presented.

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.

Does altered aortic flow in marfan syndrome relate to aortic root dilatation?

Purpose

To examine possible hemodynamic alterations in adolescent to adult Marfan syndrome (MFS) patients with aortic root dilatation.

Materials and Methods

Four‐dimensional flow MRI was performed in 20 MFS patients and 12 age‐matched normal subjects with a 3T system. The cross‐sectional areas of 10 planes along the aorta were segmented for calculating the axial and circumferential wall shear stress (WSS_axial, WSS_circ), oscillatory shear index (OSI_axial, OSI_circ), and the nonroundness (NR), presenting the asymmetry of segmental WSS. Pearson's correlation analysis was performed to present the correlations between the quantified indices and the body surface area (BSA), aortic root diameter (ARD), and Z score of the ARD. P < 0.05 indicated statistical significance.

Results

Patients exhibited lower WSS_axial in the aortic root and the WSS_circ in the arch (P < 0.05–0.001). MFS patients exhibited higher OSI_axial and OSI_circ in the sinotubular junction and arch, but lower OSI_circ in the descending aorta (all P < 0.05). The NR values were lower in patients (P < 0.05). The WSS_axial or WSS_circ exhibited moderate to strong correlations with BSA, ARD, or Z score (R² = 0.50–0.72) in MFS patients.

Conclusion

The significant differences in the quantified indices, which were associated with BSA, ARD, or Z score, in MFS were opposite to previous reports for younger MFS patients, indicating that altered flows in MFS patients may depend on the disease progress. The possible time dependency of hemodynamic alterations in MFS patients strongly suggests that longitudinal follow‐up of 4D Flow is needed to comprehend disease progress. J. Magn. Reson. Imaging 2016;44:500–508.