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

In vivo validation of 4D flow MRI for assessing the hemodynamics of portal hypertension

PURPOSE

To implement and validate in vivo radial 4D flow MRI for quantification of blood flow in the hepatic arterial, portal venous, and splanchnic vasculature of healthy volunteers and patients with portal hypertension.

MATERIALS AND METHODS

Seventeen patients with portal hypertension and seven subjects with no liver disease were included in this Health Insurance Portability and Accountability Act (HIPAA)-compliant and Institutional Review Board (IRB)-approved study. Exams were conducted at 3T using a 32-channel body coil with large volumetric coverage and 1.4 mm isotropic true spatial resolution. Using postprocessing software, cut-planes orthogonal to vessels were used to quantify flow (L/min) in the hepatic and splanchnic vasculature.

RESULTS

Flow quantification was successful in all cases. Portal vein and supraceliac aorta flow demonstrated high variability among patients. Measurements were validated indirectly using internal consistency at three different locations within the portal vein (error = 4.2 ± 3.9%) and conservation of mass at the portal confluence (error = 5.9 ± 2.5%) and portal bifurcation (error = 5.8 ± 3.1%).

CONCLUSION

This work demonstrates the feasibility of radial 4D flow MRI to quantify flow in the hepatic and splanchnic vasculature. Flow results agreed well with data reported in the literature, and conservation of mass provided indirect validation of flow quantification. Flow in patients with portal hypertensions demonstrated high variability, with patterns and magnitude consistent with the hyperdynamic state that commonly occurs in portal hypertension.

Intracranial artery velocity measurement using 4D PC MRI at 3 T: comparison with transcranial ultrasound techniques and 2D PC MRI

INTRODUCTION

4D phase contrast MR imaging (4D PC MRI) has been introduced for spatiotemporal evaluation of intracranial hemodynamics in various cerebrovascular diseases. However, it still lacks validation with standards of reference. Our goal was to compare blood flow quantification derived from 4D PC MRI with transcranial ultrasound and 2D PC MRI.

METHODS

Velocity measurements within large intracranial arteries [internal carotid artery (ICA), basilar artery (BA), and middle cerebral artery (MCA)] were obtained in 20 young healthy volunteers with 4D and 2D PC MRI, transcranial Doppler sonography (TCD), and transcranial color-coded duplex sonography (TCCD). Maximum velocities at peak systole (PSV) and end diastole (EDV) were compared using regression analysis and Bland-Altman plots.

RESULTS

Correlation of 4D PC MRI measured velocities was higher in comparison with TCD (r = 0.49-0.66) than with TCCD (0.35-0.44) and 2D PC MRI (0.52-0.60). In mid-BA and ICA C7 segment, a significant correlation was found with TCD (0.68-0.81 and 0.65-0.71, respectively). No significant correlation was found in carotid siphon. On average over all volunteers, PSVs and EDVs in MCA were minimally underestimated compared with TCD/TCCD. Minimal overestimation of velocities was found compared to TCD in mid-BA and ICA C7 segment.

CONCLUSION

4D PC MRI appears as valid alternative for intracranial velocity measurement consistent with previous reference standards, foremost with TCD. Spatiotemporal averaging effects might contribute to vessel size-dependent mild underestimation of velocities in smaller (MCA), and overestimation in larger-sized (BA and ICA) arteries, respectively. Complete spatiotemporal flow analysis may be advantageous in anatomically complex regions (e.g. carotid siphon) relative to restrictions of ultrasound techniques.

Dynamic quantitative phase imaging for biological objects using a pixelated phase mask

This paper describes research in developing a dynamic quantitative phase imaging microscope providing instantaneous measurements of dynamic motions within and among live cells without labels or contrast agents. It utilizes a pixelated phase mask enabling simultaneous measurement of multiple interference patterns derived using the polarization properties of light to track dynamic motions and morphological changes. Optical path difference (OPD) and optical thickness (OT) data are obtained from phase images. Two different processing routines are presented to remove background surface shape to enable quantification of changes in cell position and volume over time. Data from a number of different moving biological organisms and cell cultures are presented.

Volumetric motion quantification by 3D tissue phase mapped CMR

BACKGROUND

The objective of this study was the quantification of myocardial motion from 3D tissue phase mapped (TPM) CMR. Recent work on myocardial motion quantification by TPM has been focussed on multi-slice 2D acquisitions thus excluding motion information from large regions of the left ventricle. Volumetric motion assessment appears an important next step towards the understanding of the volumetric myocardial motion and hence may further improve diagnosis and treatments in patients with myocardial motion abnormalities.

METHODS

Volumetric motion quantification of the complete left ventricle was performed in 12 healthy volunteers and two patients applying a black-blood 3D TPM sequence. The resulting motion field was analysed regarding motion pattern differences between apical and basal locations as well as for asynchronous motion pattern between different myocardial segments in one or more slices. Motion quantification included velocity, torsion, rotation angle and strain derived parameters.

RESULTS

All investigated motion quantification parameters could be calculated from the 3D-TPM data. Parameters quantifying hypokinetic or asynchronous motion demonstrated differences between motion impaired and healthy myocardium.

CONCLUSIONS

3D-TPM enables the gapless volumetric quantification of motion abnormalities of the left ventricle, which can be applied in future application as additional information to provide a more detailed analysis of the left ventricular function.

A multi-center inter-manufacturer study of the temporal stability of phase-contrast velocity mapping background offset errors

BACKGROUND

Phase-contrast velocity images often contain a background or baseline offset error, which adds an unknown offset to the measured velocities. For accurate flow measurements, this offset must be shown negligible or corrected. Some correction techniques depend on replicating the clinical flow acquisition using a uniform stationary phantom, in order to measure the baseline offset at the region of interest and subtract it from the clinical study. Such techniques assume that the background offset is stable over the time of a patient scan, or even longer if the phantom scans are acquired later, or derived from pre-stored background correction images. There is no published evidence regarding temporal stability of the background offset.

METHODS

This study assessed the temporal stability of the background offset on 3 different manufacturers' scanners over 8 weeks, using a retrospectively-gated phase-contrast cine acquisition with fixed parameters and at a fixed location, repeated 5 times in rapid succession each week. A significant offset was defined as 0.6 cm/s within 50 mm of isocenter, based upon an accuracy of 10% in a typical cardiac shunt measurement.

RESULTS

Over the 5 repeated cine acquisitions, temporal drift in the baseline offset was insignificant on two machines (0.3 cm/s, 0.2 cm/s), and marginally insignificant on the third machine (0.5 cm/s) due to an apparent heating effect. Over a longer timescale of 8 weeks, insignificant drift (0.4 cm/s) occurred on one, with larger drifts (0.9 cm/s, 0.6 cm/s) on the other machines.

CONCLUSIONS

During a typical patient study, background drift was insignificant. Extended high gradient power scanning with work requires care to avoid drift on some machines. Over the longer term of 8 weeks, significant drift is likely, preventing accurate correction by delayed phantom corrections or derivation from pre-stored background offset data.

Age, gender, blood pressure, and ventricular geometry influence normal 3D blood flow characteristics in the left heart

AIMS

The aim of this study was to assess the effect of age, gender, physiological, and global cardiac function parameters on differences in normal 3D blood flow in the left ventricle (LV) and atrium (LA) using 4D flow magnetic resonance imaging (MRI).

METHODS AND RESULTS

Four-dimensional flow MRI was acquired in healthy volunteers of two age and gender groups: <30 years (6 women, n = 12) and >50 years (6 women, n = 12). Systolic and early to mid-diastolic vortex flow (number of vortices, duration, area, peak velocity inside the vortex) in the LA and LV was assessed using intra-cardiac flow visualization based on 3D particle traces and velocity vector fields. A larger number of vortices in the LA were found in young compared with older individuals (number of diastolic vortices: 1.6 ± 0.8 vs. 0.7 ± 0.7, P = 0.01) with higher velocities (54 ± 12 cm/s vs. 41 ± 11 cm/s in systole, 47 ± 13 vs. 31 ± 8 cm/s in diastole, P < 0.05). Vortices in the LV base were smaller in women compared with men (369 ± 133 vs. 543 ± 176 mm(2), P = 0.009), while vortex size was increased in mid-ventricular locations (maximum area: 546 ± 321 vs. 293 ± 174 mm(2), P < 0.05). Correlation analysis revealed significant relationships (P = 0.005-0.048, correlation coefficients = 0.44-0.84) between LA and LV vortex characteristics (number, size, vortex velocities) and blood pressure as well as end-diastolic volume, LV length, and ejection fraction.

CONCLUSIONS

Flow patterns in the left heart demonstrated differences related to age, gender, blood pressure, and ventricular geometry. The findings constitute a prerequisite for the understanding of the impact of cardiac disease on intra-cardiac haemodynamics.

Dynamic 4-dimensional microscope system with automated background leveling

This paper describes recent advances in developing an automatic background leveling algorithm for a new, novel interference microscope system and presents images and data of live biological samples. The specially designed optical system enables instantaneous 4-dimensional video measurements of dynamic motions within and among live cells without the need for contrast agents. "Label-free" measurements of biological objects in reflection using harmless light levels are possible without the need for scanning and vibration isolation. This instrument utilizes a pixelated phase mask enabling simultaneous measurement of multiple interference patterns taking advantage of the polarization properties of light enabling phase image movies in real time at video rates to track dynamic motions and volumetric changes. Optical thickness data are derived from phase images. This data is processed with an automatic background leveling routine which separates the objects from the background by thresholding the calculated gradient magnitude of the optical thickness data. Low-order Zernike surfaces are fit to the unmasked background pixels and the resulting background shape is removed. This method effectively eliminates background shape for datasets containing both large and small objects. By applying this method to many sequential frames, it results in all the frames having the same mean background value across all frames which is essential for quantitatively montoring time-dependent processes.

Evaluation of valvular insufficiency and shunts with parallel-imaging compressed-sensing 4D phase-contrast MR imaging with stereoscopic 3D velocity-fusion volume-rendered visualization

PURPOSE

To assess the potential of compressed-sensing parallel-imaging four-dimensional (4D) phase-contrast magnetic resonance (MR) imaging and specialized imaging software in the evaluation of valvular insufficiency and intracardiac shunts in patients with congenital heart disease.

MATERIALS AND METHODS

Institutional review board approval was obtained for this HIPAA-compliant study. Thirty-four consecutive retrospectively identified patients in whom a compressed-sensing parallel-imaging 4D phase-contrast sequence was performed as part of routine clinical cardiac MR imaging between March 2010 and August 2011 and who had undergone echocardiography were included. Multiplanar, volume-rendered, and stereoscopic three-dimensional velocity-fusion visualization algorithms were developed and implemented in Java and OpenGL. Two radiologists independently reviewed 4D phase-contrast studies for each of 34 patients (mean age, 6 years; age range, 10 months to 21 years) and tabulated visible shunts and valvular regurgitation. These results were compared with color Doppler echocardiographic and cardiac MR imaging reports, which were generated without 4D phase-contrast visualization. Cohen κ statistics were computed to assess interobserver agreement and agreement with echocardiographic results.

RESULTS

The 4D phase-contrast acquisitions were performed, on average, in less than 10 minutes. Among 123 valves seen in 34 4D phase-contrast studies, 29 regurgitant valves were identified, with good agreement between observers (k=0.85). There was also good agreement with the presence of at least mild regurgitation at echocardiography (observer 1, κ=0.76; observer 2, κ=0.77) with high sensitivity (observer 1, 75%; observer 2, 82%) and specificity (observer 1, 97%; observer 2, 95%) relative to the reference standard. Eight intracardiac shunts were identified, four of which were not visible with conventional cardiac MR imaging but were detected with echocardiography. No intracardiac shunts were found with echocardiography alone.

CONCLUSION

With velocity-fusion visualization, the compressed-sensing parallel-imaging 4D phase-contrast sequence can augment conventional cardiac MR imaging by improving sensitivity for and depiction of hemodynamically significant shunts and valvular regurgitation.

Sparsity transform k-t principal component analysis for accelerating cine three-dimensional flow measurements

Time-resolved three-dimensional flow measurements are limited by long acquisition times. Among the various acceleration techniques available, k-t methods have shown potential as they permit significant scan time reduction even with a single receive coil by exploiting spatiotemporal correlations. In this work, an extension of k-t principal component analysis is proposed utilizing signal differences between the velocity encodings of three-directional flow measurements to further compact the signal representation and hence improve reconstruction accuracy. The effect of sparsity transform in k-t principal component analysis is demonstrated using simulated and measured data of the carotid bifurcation. Deploying sparsity transform for 8-fold undersampled simulated data, velocity root-mean-square errors were found to decrease by 52 ± 14%, 59 ± 11%, and 16 ± 32% in the common, external, and internal carotid artery, respectively. In vivo, errors were reduced by 15 ± 17% in the common carotid artery with sparsity transform. Based on these findings, spatial resolution of three-dimensional flow measurements was increased to 0.8 mm isotropic resolution with prospective 8-fold undersampling and sparsity transform k-t principal component analysis reconstruction. Volumetric data were acquired in 6 min. Pathline visualization revealed details of helical flow patterns partially hidden at lower spatial resolution.

Combination of tagging and tissue phase mapping to accelerate myocardial motion measurements in three directions

OBJECT

Until now, a three-directional velocity field has mostly been obtained by velocity encoding in three directions, which is very time-consuming and hence not usually used in clinical routine. We show the feasibility of combining in-plane tagging with through-plane tissue phase mapping (TPM) to encode a three-directional velocity field at 3 T with reduced overall acquisition time.

MATERIALS AND METHODS

Assessment of a three-directional velocity field was performed for 10 healthy volunteers. The motion patterns obtained by use of five different sequences including three-directional TPM, TPM in the through-plane direction, TPM in the through-plane direction with horizontal or vertical tagging lines, and TPM in the through-plane direction combined with a tagging grid were evaluated and compared.

RESULTS

A three-dimensional velocity field can be obtained in approximately half the acquisition time by combining through-plane TPM with in-plane tagging. Although the velocity information is derived by different means, differences between the information obtained by three-directional TPM encoding and the suggested technique are only minor.

CONCLUSION

The combination of tagging and TPM enables assessment of the three-directional velocity field in nearly half the time taken when the conventional three-directional TPM sequence is used.

Rapid pediatric cardiac assessment of flow and ventricular volume with compressed sensing parallel imaging volumetric cine phase-contrast MRI

OBJECTIVE

The quantification of cardiac flow and ventricular volumes is an essential goal of many congenital heart MRI examinations, often requiring acquisition of multiple 2D phase-contrast and bright-blood cine steady-state free precession (SSFP) planes. Scan acquisition, however, is lengthy and highly reliant on an imager who is well-versed in structural heart disease. Although it can also be lengthy, 3D time-resolved (4D) phase-contrast MRI yields global flow patterns and is simpler to perform. We therefore sought to accelerate 4D phase contrast and to determine whether equivalent flow and volume measurements could be extracted.

MATERIALS AND METHODS

Four-dimensional phase contrast was modified for higher acceleration with compressed sensing. Custom software was developed to process 4D phase-contrast images. We studied 29 patients referred for congenital cardiac MRI who underwent a routine clinical protocol, including cine short-axis stack SSFP and 2D phase contrast, followed by contrast-enhanced 4D phase contrast. To compare quantitative measurements, Bland-Altman analysis, paired Student t tests, and F tests were used.

RESULTS

Ventricular end-diastolic, end-systolic, and stroke volumes obtained from 4D phase contrast and SSFP were well correlated (ρ = 0.91-0.95; r(2) = 0.83-0.90), with no statistically significant difference. Ejection fractions were well correlated in a subpopulation that underwent higher-resolution compressed-sensing 4D phase contrast (ρ = 0.88; r(2) = 0.77). Four-dimensional phase contrast and 2D phase contrast flow rates were also well correlated (ρ = 0.90; r(2) = 0.82). Excluding ventricles with valvular insufficiency, cardiac outputs derived from outlet valve flow and stroke volumes were more consistent by 4D phase contrast than by 2D phase contrast and SSFP.

CONCLUSION

Combined parallel imaging and compressed sensing can be applied to 4D phase contrast. With custom software, flow and ventricular volumes may be extracted with comparable accuracy to SSFP and 2D phase contrast. Furthermore, cardiac outputs were more consistent by 4D phase contrast.

Normal and altered three-dimensional portal venous hemodynamics in patients with liver cirrhosis

PURPOSE

To compare time-resolved three-dimensional (3D) phase-contrast magnetic resonance (MR) imaging with three-directional velocity encoding (flow-sensitive four-dimensional [4D] MR imaging), with Doppler ultrasonography (US) as standard of reference, for investigating alterations in 3D portal venous hemodynamics in patients with liver cirrhosis compared with healthy age-matched control subjects and healthy young volunteers.

MATERIAL & METHODS

This prospective study was approved by the local ethics committee, and written informed consent was obtained from all participants. Three-dimensional portal venous hemodynamics was assessed, employing flow-sensitive 4D MR imaging with a 3-T MR system (spatial resolution, approximately 2 mm(3); temporal resolution, approximately 45 msec) in 20 patients with hepatic cirrhosis, 20 healthy age-matched control subjects, and 21 healthy young volunteers. Flow characteristics were analyzed by using 3D streamlines and time-resolved particle traces. Quantitative analyses were performed by retrospectively evaluating regional peak and mean velocities, flow volume, and vessel area. Doppler US was used as standard of reference. Independent-sample t tests or Wilcoxon-Mann-Whitney tests were applied for comparing each subject group. Paired-sample t tests or Wilcoxon tests were applied when comparing MR imaging and US.

RESULTS

Three-dimensional visualization of portal venous hemodynamics was successful, with complete visualization of the vessels in 18 patients and 35 volunteers, with limitations in the left intrahepatic branches (87%, reader A; 89%, reader B). A moderate but significant correlation was observed between 4D MR imaging and Doppler US in nearly all maximum and mean velocities, flow volumes, and vessel areas (r = 0.24-0.64, P = .001-.044). With MR imaging, significant underestimation was observed of intrahepatic flow velocities and flow volumes, except vessel area, which Doppler US represented as even lower (P < .001 to P = .045). Six patients had collateralization with reopened umbilical vein, while one had flow reversal in the superior mesenteric vein visible at MR imaging only.

CONCLUSION

Flow-sensitive 4D MR imaging may constitute a promising, alternative technique to Doppler US for evaluating hemodynamics in the portal venous system of patients with liver cirrhosis and may be a means of assessing pathologic changes in flow characteristics.

Improved MR phase-contrast velocimetry using a novel nine-point balanced motion-encoding scheme with increased robustness to eddy current effects

Phase-contrast MRI (PC-MRI) velocimetry is a noninvasive, high-resolution motion assessment tool. However, high motion sensitivity requires strong motion-encoding magnetic gradients, making phase-contrast-MRI prone to baseline shift artifacts due to the generation of eddy currents. In this study, we propose a novel nine-point balanced velocity-encoding strategy, designed to be more accurate in the presence of strong and rapidly changing gradients. The proposed method was validated using a rotating phantom, and its robustness and precision were explored and compared with established approaches through computer simulations and in vivo experiments. Computer simulations yielded a 39-57% improvement in velocity-noise ratio (corresponding to a 27-33% reduction in measurement error), depending on which method was used for comparison. Moreover, in vivo experiments confirmed this by demonstrating a 26-53% reduction in accumulated velocity error over the R-R interval. The nine-point balanced phase-contrast-MRI-encoding strategy is likely useful for settings where high spatial and temporal resolution and/or high motion sensitivity is required, such as in high-resolution rodent myocardial tissue phase mapping.

Phase-contrast MRI and applications in congenital heart disease

A review of phase-contrast magnetic resonance imaging techniques, with specific application to congenital heart disease, is presented. Theory, pitfalls, advantages, and specific examples of multiple, well-described congenital heart disease presentations are discussed.

Improved method for quantification of regional cardiac function in mice using phase-contrast MRI

Phase-contrast magnetic resonance imaging is a technique that allows for characterization of regional cardiac function and for measuring transmural myocardial velocities in human hearts with high temporal and spatial resolution. The application of this technique (also known as tissue phase mapping) to murine hearts has been very limited so far. The aim of our study was to implement and to optimize tissue phase mapping for a comprehensive assessment of murine transmural wall motion. Baseline values for regional motion patterns in mouse hearts, based on the clinically used American Heart Association's 17-segment model, were established, and a detailed motion analysis of mouse heart for the entire cardiac cycle (including epicardial and endocardial motion patterns) is provided. Black-blood contrast was found to be essential to obtain reproducible velocity encoding. Tissue phase mapping of the mouse heart permits the detailed assessment of regional myocardial velocities. While a proof-of-principle application in a murine ischemia-reperfusion model was performed, future studies are warranted to assess its potential for the investigation of systolic and diastolic functions in genetically and surgically manipulated mouse models of human heart disease.

Four-dimensional phase contrast MRI with accelerated dual velocity encoding

PURPOSE

To validate a novel approach for accelerated four-dimensional phase contrast MR imaging (4D PC-MRI) with an extended range of velocity sensitivity.

MATERIALS AND METHODS

4D PC-MRI data were acquired with a radially undersampled trajectory (PC-VIPR). A dual V(enc) (dV(enc) ) processing algorithm was implemented to investigate the potential for scan time savings while providing an improved velocity-to-noise ratio. Flow and velocity measurements were compared with a flow pump, conventional 2D PC MR, and single V(enc) 4D PC-MRI in the chest of 10 volunteers.

RESULTS

Phantom measurements showed excellent agreement between accelerated dV(enc) 4D PC-MRI and the pump flow rate (R(2) ≥ 0.97) with a three-fold increase in measured velocity-to-noise ratio (VNR) and a 5% increase in scan time. In volunteers, reasonable agreement was found when combining 100% of data acquired with V(enc) = 80 cm/s and 25% of the high V(enc) data, providing the VNR of a 80 cm/s acquisition with a wider velocity range of 160 cm/s at the expense of a 25% longer scan.

CONCLUSION

Accelerated dual V(enc) 4D PC-MRI was demonstrated in vitro and in vivo. This acquisition scheme is well suited for vascular territories with wide ranges of flow velocities such as congenital heart disease, the hepatic vasculature, and others.

Analysis of pulse wave velocity in the thoracic aorta by flow-sensitive four-dimensional MRI: reproducibility and correlation with characteristics in patients with aortic atherosclerosis

PURPOSE

To measure aortic pulse wave velocity (PWV) using flow-sensitive four-dimensional (4D) MRI and to evaluate test-retest reliability, inter- and intra-observer variability in volunteers and correlation with characteristics in patients with aortic atherosclerosis.

MATERIALS AND METHODS

Flow-sensitive 4D MRI was performed in 12 volunteers (24 ± 3 years) and 86 acute stroke patients (68 ± 9 years) with aortic atherosclerosis. Retrospectively positioned 28 ± 4 analysis planes along the entire aorta (inter-slice-distance = 10 mm) and frame wise lumen segmentation yielded flow-time-curves for each plane. Global aortic PWV was calculated from time-shifts and distances between the upslope portions of all available flow-time curves.

RESULTS

Inter- and intra-observer variability of PWV measurements in volunteers (7% and 8%) was low while test-retest reliability (22%) was moderate. PWV in patients was significantly higher compared with volunteers (5.8 ± 2.9 versus 3.8 ± 0.8 m/s; P = 0.02). Among 17 patient characteristics considered, statistical analysis revealed significant (P < 0.05) but low correlation of PWV with age (r = 0.25), aortic valve insufficiency (r = 0.29), and pulse pressure (r = 0.28). Multivariate modeling indicated that aortic valve insufficiency and elevated pulse pressure were significantly associated with higher PWV (adjusted R(2) = 0.13).

CONCLUSION

Flow-sensitive 4D MRI allows for estimating aortic PWV with low observer dependence and moderate test-retest reliability. PWV in patients correlated with age, aortic valve insufficiency, and pulse pressure.

Interdependencies of aortic arch secondary flow patterns, geometry, and age analysed by 4-dimensional phase contrast magnetic resonance imaging at 3 Tesla

OBJECTIVE

It was the aim to analyse the impact of age, aortic arch geometry, and size on secondary flow patterns such as helix and vortex flow derived from flow-sensitive magnetic resonance imaging (4D PC-MRI).

METHODS

62 subjects (age range = 20-80 years) without circumscribed pathologies of the thoracic aorta (ascending aortic (AAo) diameter: 3.2 ± 0.6 cm [range 2.2-5.1]) were examined by 4D PC-MRI after IRB-approval and written informed consent. Blood flow visualisation based on streamlines and time-resolved 3D particle traces was performed. Aortic diameter, shape (gothic, crook-shaped, cubic), angle, and age were correlated with existence and extent of secondary flow patterns (helicity, vortices); statistical modelling was performed.

RESULTS

Helical flow was the typical pattern in standard crook-shaped aortic arches. With altered shapes and increasing age, helicity was less common. AAo diameter and age had the highest correlation (r = 0.69 and 0.68, respectively) with number of detected vortices. None of the other arch geometric or demographic variables (for all, P ≥ 0.177) improved statistical modelling.

CONCLUSION

Substantially different secondary flow patterns can be observed in the normal thoracic aorta. Age and the AAo diameter were the parameters correlating best with presence and amount of vortices. Findings underline the importance of age- and geometry-matched control groups for haemodynamic studies.

KEY POINTS

• Secondary blood flow patterns (helices, vortices) are commonly observed in the aorta • Secondary flow patterns predominantly depend on patient age and aortic diameter • Geometric factors show a lesser impact on blood flow patterns than age and diameter • Future analyses of flow patterns should incorporate age- and diameter dependencies.

Four-dimensional flow MRI using spiral acquisition

Time-resolved three-dimensional phase-contrast MRI is an important tool for physiological as well as clinical studies of blood flow in the heart and vessels. The application of the technique is, however, limited by the long scan times required. In this work, we investigate the feasibility of using spiral readouts to reduce the scan time of four-dimensional flow MRI without sacrificing quality. Three spiral approaches are presented and evaluated in vivo and in vitro against a conventional Cartesian acquisition. In vivo, the performance of each method was assessed in the thoracic aorta in 10 volunteers using pathline-based analysis and cardiac output analysis. Signal-to-noise ratio and background phase errors were investigated in vitro. Using spiral readouts, the scan times of a four-dimensional flow acquisition of the thoracic aorta could be reduced 2-3-fold, with no statistically significant difference in pathline validity or cardiac output. The shortened scan time improves the applicability of four-dimensional flow MRI, which may allow the technique to become a part of a clinical workflow for cardiovascular functional imaging.

Acceleration of tissue phase mapping with sensitivity encoding at 3T

BACKGROUND

The objective of this study was to investigate the impact of sensitivity encoding on the quantitative assessment of cardiac motion in black blood cine tissue phase mapping (TPM) sequences. Up to now whole volume coverage of the heart is still limited by the long acquisition times. Therefore, a significant increase in imaging speed without deterioration of quantitative motion information is indispensable.

METHODS

20 volunteers were enrolled in this study. Each volunteer underwent myocardial short-axis TPM scans with different SENSE acceleration factors. The influence of SENSE acceleration on the measured motion curves was investigated.

RESULTS

It is demonstrated that all TPM sequences with SENSE acceleration have only minimum influence on the motion curves. Even with a SENSE factor of four, the decrease in the amplitude of the motion curve was less than 3%. No significant difference was observed for the global correlation coefficient and deviation between the motion curves obtained by the reproducibility and the SENSE accelerated measurements.

CONCLUSIONS

It is feasible to accelerate myocardial TPM measurements with SENSE factors up to 4 without losing substantial information of the motion pattern.

Comparison of four-dimensional flow parameters for quantification of flow eccentricity in the ascending aorta

PURPOSE

To compare quantitative parameters for assessing the degree of eccentric systolic blood flow in the ascending thoracic aorta (AsAo).

MATERIALS AND METHODS

Forty-one patients were studied with three-dimensional (3D), cine phase-contract MRI (4D Flow). Analysis was performed at peak systole for a cross-sectional plane in the AsAo just distal to the sinotubular junction. AsAo flow was graded as normal, mildly, or markedly eccentric based on qualitative visual assessment. For quantitative analysis, flow jet angle and normalized flow displacement from the vessel center were calculated.

RESULTS

Patients with normal AsAo systolic flow (n = 25) had an average flow jet angle of 13.7 degrees and flow displacement 0.04. These parameters were significantly elevated for patients with mild eccentric systolic flow (n = 6): 24.6 degrees (P = 0.012) and 0.12 (P = 0.001), respectively. However, for patients with marked eccentric flow (n = 10), only flow displacement was significantly elevated compared with the mild eccentric group (0.18; P = 0.04); flow angle was 25.7 degrees.

CONCLUSION

Flow displacement is a more reliable quantitative parameter for measuring eccentric AsAo systolic flow than flow jet angle, and should be evaluated in studies investigating the role of eccentric flow in the promotion of aortic pathology.

Reproducibility of flow and wall shear stress analysis using flow-sensitive four-dimensional MRI

PURPOSE

To systematically investigate the scan-rescan reproducibility and observer variability of flow-sensitive four-dimensional (4D) MRI in the aorta for the assessment of blood flow and global and segmental wall shear stress.

MATERIALS AND METHODS

ECG and respiration-synchronized flow-sensitive 4D MRI data (spatio-temporal resolution = 1.7 × 2.0 × 2.2 mm(3) /40.8 ms) were acquired in 12 healthy volunteers. To analyze scan-rescan variability, flow-sensitive 4D MRI was repeated in 10 volunteers during a second visit. Data analysis included calculation of time-resolved and total flow, peak systolic velocity, and regional and global wall shear stress (WSS) in up to 24 analysis planes distributed along the aorta.

RESULTS

Scan-rescan, inter-observer, and intra-observer agreement was excellent for the calculation of total flow and peak systolic velocity (mean differences <5% of the average flow parameter). Global WSS demonstrated moderate agreement and increased variability regarding wall shear stress (scan-rescan, inter-observer, and intra-observer agreement; mean differences <10% of the average WSS parameters). The segmental distribution of wall shear stress in the thoracic aorta could reliably be reproduced (r > 0.87; P < 0.001) for different observers and examinations.

CONCLUSION

Flow-sensitive 4D MRI-based analysis of aortic blood flow can be performed with good reproducibility. Robustness of global and regional WSS quantification was limited, but spatio-temporal WSS distributions could reliably be replicated.

Analysis and correction of background velocity offsets in phase-contrast flow measurements using magnetic field monitoring

The value of phase-contrast magnetic resonance imaging for quantifying tissue motion and blood flow has been long recognized. However, the sensitivity of the method to system imperfections can lead to inaccuracies limiting its clinical acceptance. A key source of error relates to eddy current-induced phase fluctuations, which can offset the measured object velocity significantly. A higher-order dynamic field camera was used to study the spatiotemporal evolution of background phases in cine phase-contrast measurements. It is demonstrated that eddy current-induced offsets in phase-difference data are present up to the second spatial order. Oscillatory temporal behaviors of offsets in the kHz range suggest mechanical resonances of the MR system to be non-negligible in phase-contrast imaging. By careful selection of the echo time, their impact can be significantly reduced. When applying field monitoring data for correcting eddy current and mechanically induced velocity offsets, errors decrease to less than 0.5% of the maximum velocity for various sequence settings proving the robustness of the correction approach. In vivo feasibility is demonstrated for aortic and pulmonary flow measurements in five healthy subjects. Using field monitoring data, mean error in stroke volume was reduced from 10% to below 3%.

Dampening of blood-flow pulsatility along the carotid siphon: does form follow function?

BACKGROUND AND PURPOSE

The tortuous distal part of the ICA may have an attenuating effect on pulsatile arterial flow. We investigated local arterial blood flow patterns in the ICA proximal and distal to the carotid siphon to detect quantitative waveform changes.

MATERIALS AND METHODS

Arterial flow patterns were analyzed by using flow-sensitized 4D PC MR imaging (time-resolved 3D PCMR) at 3T in 17 healthy volunteers. Time-resolved blood flow velocities were extracted from the source data at the C4 and C7 segments of the ICA. PI, RI, and PA were calculated by using time-resolved flow volume. A linear mixed-effects model was applied to compare values at C4 and C7. Furthermore, 3D blood flow visualization was performed for all 34 ICAs.

RESULTS

PI, RI, and PA were significantly lower at the distal C7 segment compared with the proximal C4 segment of the ICA (P < .0001). Helical flow patterns were observed in 5 ICAs of 4 subjects.

CONCLUSIONS

Arterial flow patterns showed a significant reduction in PI, RI, and PA when compared distal to proximal to the carotid siphon. The observed attenuation of flow pulsatility is most likely related to the contorted shape of the distal ICA and may bear a protective effect for downstream cerebral vasculature.

Improved cardiovascular flow quantification with time-resolved volumetric phase-contrast MRI

BACKGROUND

Cardiovascular flow is commonly assessed with two-dimensional, phase-contrast MRI (2-D PC-MRI). However, scan prescription and acquisition over multiple planes is lengthy, often requires direct physician oversight and has inconsistent results. Time-resolved volumetric PC-MRI (4-D flow) may address these limitations.

OBJECTIVE

We assess the degree of agreement and internal consistency between 2-D and 4-D flow quantification in our clinical population.

MATERIALS AND METHODS

Software enabling flow calculation from 4-D flow was developed in Java. With IRB approval and HIPAA compliance, 18 consecutive patients without shunts were identified who underwent both (1) conventional 2-D PC-MRI of the aorta and main pulmonary artery and (2) 4-D flow imaging. Aortic and pulmonary flow rates were assessed with both techniques.

RESULTS

Both methods showed general agreement in flow rates (ρ: 0.87-0.90). Systemic and pulmonary arterial flow rates were well-correlated (ρ: 4-D 0.98-0.99, 2-D 0.93), but more closely matched with 4-D (P < 0.05, Brown-Forsythe). Pulmonary flow rates were lower than systemic rates for 2-D (P < 0.05, two-sample t-test). In a sub-analysis of patients without pulmonary or aortic regurgitation, 2-D showed improved correlation of flow rates while 4-D phase-contrast remained tightly correlated (ρ: 4-D 0.99-1.00, 2-D 0.99).

CONCLUSION

4-D PC-MRI demonstrates greater consistency than conventional 2-D PC-MRI for flow quantification.

Sequence optimization to reduce velocity offsets in cardiovascular magnetic resonance volume flow quantification--a multi-vendor study

PURPOSE

Eddy current induced velocity offsets are of concern for accuracy in cardiovascular magnetic resonance (CMR) volume flow quantification. However, currently known theoretical aspects of eddy current behavior have not led to effective guidelines for the optimization of flow quantification sequences. This study is aimed at identifying correlations between protocol parameters and the resulting velocity error in clinical CMR flow measurements in a multi-vendor study.

METHODS

Nine 1.5T scanners of three different types/vendors were studied. Measurements were performed on a large stationary phantom. Starting from a clinical breath-hold flow protocol, several protocol parameters were varied. Acquisitions were made in three clinically relevant orientations. Additionally, a time delay between the bipolar gradient and read-out, asymmetric versus symmetric velocity encoding, and gradient amplitude and slew rate were studied in adapted sequences as exploratory measurements beyond the protocol. Image analysis determined the worst-case offset for a typical great-vessel flow measurement.

RESULTS

The results showed a great variation in offset behavior among scanners (standard deviation among samples of 0.3, 0.4, and 0.9 cm/s for the three different scanner types), even for small changes in the protocol. Considering the absolute values, none of the tested protocol settings consistently reduced the velocity offsets below the critical level of 0.6 cm/s neither for all three orientations nor for all three scanner types. Using multilevel linear model analysis, oblique aortic and pulmonary slices showed systematic higher offsets than the transverse aortic slices (oblique aortic 0.6 cm/s, and pulmonary 1.8 cm/s higher than transverse aortic). The exploratory measurements beyond the protocol yielded some new leads for further sequence development towards reduction of velocity offsets; however those protocols were not always compatible with the time-constraints of breath-hold imaging and flow-related artefacts.

CONCLUSIONS

This study showed that with current systems there was no generic protocol which resulted into acceptable flow offset values. Protocol optimization would have to be performed on a per scanner and per protocol basis. Proper optimization might make accurate (transverse) aortic flow quantification possible for most scanners. Pulmonary flow quantification would still need further (offline) correction.

SAR reduced black-blood cine TPM for increased temporal resolution at 3T

OBJECT

The objective was to improve the temporal resolution in black-blood CINE tissue phase mapping sequences at high field MR systems. The temporal resolution is limited due to SAR constraints causing idle times into the sequence. The aim was to avoid these idle times and therefore providing an increased number of heart phases.

MATERIALS AND METHODS

Thirteen volunteers were enrolled in this study. Each volunteer underwent different myocardial short-axis scans comprising scans with application of both presaturation pulses, with alternating application of presaturation pulses and with an attenuation of the excitation angle. The last two approaches enable a SAR reduction or increased temporal resolution. The contrast to noise ratio (CNR) between myocardium and blood and the influence on the measured tissue motion were investigated.

RESULTS

High CNR between myocardium and blood could be obtained with the application of alternating presaturation-pulses. Reduction of the flip angle of the presaturation-pulses provided reduced CNR relative to both the original and the alternated presaturation-pulses approach. More details of the myocardial motion were observed with increased temporal resolution.

CONCLUSION

It is feasible to increase the temporal resolution at high field strength by reducing the SAR with either alternating presaturation-pulses or decreased flip angle of these pulses.

Comprehensive 4D velocity mapping of the heart and great vessels by cardiovascular magnetic resonance

BACKGROUND

Phase contrast cardiovascular magnetic resonance (CMR) is able to measure all three directional components of the velocities of blood flow relative to the three spatial dimensions and the time course of the heart cycle. In this article, methods used for the acquisition, visualization, and quantification of such datasets are reviewed and illustrated.

METHODS

Currently, the acquisition of 3D cine (4D) phase contrast velocity data, synchronized relative to both cardiac and respiratory movements takes about ten minutes or more, even when using parallel imaging and optimized pulse sequence design. The large resulting datasets need appropriate post processing for the visualization of multidirectional flow, for example as vector fields, pathlines or streamlines, or for retrospective volumetric quantification.

APPLICATIONS

Multidirectional velocity acquisitions have provided 3D visualization of large scale flow features of the healthy heart and great vessels, and have shown altered patterns of flow in abnormal chambers and vessels. Clinically relevant examples include retrograde streams in atheromatous descending aortas as potential thrombo-embolic pathways in patients with cryptogenic stroke and marked variations of flow visualized in common aortic pathologies. Compared to standard clinical tools, 4D velocity mapping offers the potential for retrospective quantification of flow and other hemodynamic parameters.

CONCLUSIONS

Multidirectional, 3D cine velocity acquisitions are contributing to the understanding of normal and pathologically altered blood flow features. Although more rapid and user-friendly strategies for acquisition and analysis may be needed before 4D velocity acquisitions come to be adopted in routine clinical CMR, their capacity to measure multidirectional flows throughout a study volume has contributed novel insights into cardiovascular fluid dynamics in health and disease.

Acceleration of tissue phase mapping by k-t BLAST: a detailed analysis of the influence of k-t-BLAST for the quantification of myocardial motion at 3T

BACKGROUND

The assessment of myocardial motion with tissue phase mapping (TPM) provides high spatiotemporal resolution and quantitative motion information in three directions. Today, whole volume coverage of the heart by TPM encoding at high spatial and temporal resolution is limited by long data acquisition times. Therefore, a significant increase in imaging speed without deterioration of the quantitative motion information is required. For this purpose, the k-t BLAST acceleration technique was combined with TPM black-blood functional imaging of the heart. Different k-t factors were evaluated with respect to their impact on the quantitative assessment of cardiac motion.

RESULTS

It is demonstrated that a k-t BLAST factor of two can be used with a marginal, but statistically significant deterioration of the quantitative motion data. Further increasing the k-t acceleration causes substantial alteration of the peak velocities and the motion pattern, but the temporal behavior of the contraction is well maintained up to an acceleration factor of six.

CONCLUSIONS

The application of k-t BLAST for the acceleration of TPM appears feasible. A reduction of the acquisition time of almost 45% could be achieved without substantial loss of quantitative motion information.

MR-based visualization and quantification of three-dimensional flow characteristics in the portal venous system

PURPOSE

To evaluate the feasibility of time-resolved flow-sensitive MRI for the three-dimensional (3D) visualization and quantification of normal and pathological portal venous (PV) hemodynamics.

MATERIALS AND METHODS

Portal venous hemodynamics were evaluated in 18 healthy volunteers and 5 patients with liver cirrhosis. ECG- and adaptive respiratory navigator gated flow-sensitive 4D MRI (time-resolved 3D MRI with three-directional velocity encoding) was performed on a 3 Tesla MR system (TRIO, Siemens, Germany). Qualitative flow analysis was achieved using 3D streamlines and time-resolved particle traces originating from seven emitter planes precisely placed at anatomical landmarks in the PV system. Quantitative analysis included retrospective extraction of regional peak and mean velocities and vessel area. Results were compared with standard 2D flow-sensitive MRI and to the reference standard Doppler ultrasound.

RESULTS

Qualitative flow analysis was successfully used in the entire PV system. Venous hemodynamics in all major branches in 17 of 18 volunteers and 3 of 5 patients were reliably depicted with good interobserver agreement (kappa = 0.62). Quantitative analysis revealed no significant differences and moderate agreement for peak velocities between 3D MR and 2D MRI (r = 0.46) and Doppler ultrasound (US) (r = 0.35) and for mean velocities between 3D and 2D MRI (r = 0.41). The PV area was significantly (P < 0.01) higher in 3D and 2D MRI compared with US.

CONCLUSION

We successfully applied 3D MR velocity mapping in the PV system, providing a detailed qualitative and quantitative analysis of normal and pathological hemodynamics.

Clinical evaluation of aortic coarctation with 4D flow MR imaging

PURPOSE

To show that 4D Flow is a clinically viable tool for evaluation of collateral blood flow and demonstration of distorted blood flow patterns in patients with treated and untreated aortic coarctation.

MATERIALS AND METHODS

Time-resolved, 3D phase contrast magnetic resonance imaging (MRI) (4D Flow) was used to assess blood flow in the thoracic aorta of 34 individuals: 26 patients with coarctation (22 after surgery or stent placement) and eight healthy volunteers.

RESULTS

Direct comparison of blood flow calculated with 2D and 4D phase contrast data at standard levels for analysis in coarctation patients showed good correlation and agreement (correlation coefficient r = 0.99, limits of agreement = -20% to 20% for collateral blood flow calculations). Abnormal blood flow patterns were demonstrated at peak systole with 4D Flow visualization techniques in the descending thoracic aorta of patients but not volunteers. Marked helical flow was seen in 9 of 13 patients with angulated aortic arch geometries after coarctation repair. Vortical flow was seen in regions of poststenotic dilation.

CONCLUSION

4D Flow is a fast and reliable means of evaluating collateral blood flow in patients with aortic coarctation in order to establish hemodynamic significance. It also can detect distorted blood flow patterns in the descending aorta after coarctation repair.

Improved accuracy in flow mapping of congenital heart disease using stationary phantom technique

BACKGROUND

Flow mapping by cardiovascular magnetic resonance has become the gold standard for non-invasively defining cardiac output (CO), shunt flow and regurgitation. Previous reports have highlighted the presence of inherent errors in flow mapping that are improved with the use of a stationary phantom control. To our knowledge, these studies have only been performed in healthy volunteers.

RESULTS

We analyzed the variation in flow measurements made with and without stationary phantom correction in 31 patients with congenital heart disease. Variation in stroke volume (SV) measurements was seen in all vessels across all patient groups. The variation was largest when analyzing the right ventricular outflow tract (RVOT), with a range of absolute differences in SV from 0.2 to 70 ml and in CO from 0.02 to 4.8 L/min. In patients with repaired Tetrology of Fallot (ToF), the average ratio of pulmonary to systemic blood flow (Qp:Qs) was 1.18 without and 1.02 with phantom correction. Without performing phantom correction, 23% of the repaired ToF patients were classified as having a residual shunt as compared to 0% when flow mapping was performed with phantom correction. Similarly, in patients with known atrial level shunting (ASD/PAPVR) 20% of patients had no shunt when flow mapping was performed without phantom correction as compared to 0% with phantom correction. In patients with bicuspid aortic valves (BAV), the differences in the regurgitant fraction between measuring flow with and without phantom correction ranged from 0 to 30%, while the regurgitant fraction in the RVOT of ToF patients varied by as much as 31%.

CONCLUSION

The impact of inherent errors in CMR flow mapping should not be underestimated. While the variation across a population may not display a significant trend, for any individual patient it can be quite large. Failure to correct for such variation can lead to clinically significant misinterpretation of flow data. The use of the stationary phantom correction technique appears to improve accuracy both in normal patients as well as those with congenital heart disease.

[Analysis of flow in artificial stenosis models of mid-sized arteries using 3D PC-MRI]

Phase contrast MRI allows access to tri-directional encoded velocity information and therefore, measurement of flow in the human hemodynamic system. The aim of this work was to investigate whether this technology could be applied to support the grading of stenosis in mid-size arteries. Using a specially constructed flow phantom and a stenosis model with tube diameter of 5mm and 8mm and a stenosis of 50%, experiments at different flow rates (180-640 ml/min), slice thickness (1-4 mm), field strength (1.5 and 3.0 T), and multi-slice as well as 3D volume acquisition were performed. The observations were assessed visually and evaluated by signal-to-noise (SNR) ratios in regions before and after the stenosis. The obtained results show that examinations should be performed at high field (3.0 T) and at flow rates up to 500 ml/min without hampering the measurements by areas of signal loss. In comparison, no detectable differences in the flow patterns of the two acquisition schemes could be observed. However, the SNR was higher using the 3D volume acquisition and thick slices. In summary, 3D PC-MRI of mid-size vessels with stenosis is feasible for certain flow rates. The presented results could be seen as guidance for in vivo situations to assess if an examination of a patient is reasonable in terms of outcome.

Regional pulmonary blood flow: Comparison of dynamic contrast-enhanced MR perfusion and phase-contrast MR

Regional pulmonary blood flow can be assessed using both dynamic contrast-enhanced (DCE) MR and phase-contrast (PC) MR. These methods provide somewhat complementary information: DCE MR can assess flow over the entire lung while PC MR can detect rapid changes in flow to a targeted region. Although both methods are considered accurate, one may be more feasible than the other depending on pathology, patient condition, and availability of an intravenous route. The objective of this study was to establish a consensus between the two methods by comparing paired DCE MR and PC MR measurements of relative blood flow in Yorkshire piglets (N = 5, age = 7 days, weight = 3.3 +/- 0.6 kg) under various physiological states including regional lung collapse. A strong correlation (R(2) = 0.71, P < 0.01) was observed between the methods. In conclusion, DCE MR and PC MR provide a consistent measure of changes in regional pulmonary blood flow.

Evaluation of intracranial stenoses and aneurysms with accelerated 4D flow

The aim of this study was to evaluate intracranial arterial stenoses and aneurysms with accelerated time-resolved three-dimensional (3D) phase-contrast MRI or 4D flow. The 4D flow technique was utilized to image four normal volunteers, two patients with intracranial stenoses and two patients with intracranial aneurysms. In order to reduce scan time, parallel imaging was combined with an acquisition strategy that eliminates the corners of k-space. In the two patients with intracranial stenoses, 4D flow velocity measurements showed that one patient had normal velocity profiles in agreement with a previous magnetic resonance angiogram (MRA), while the second showed increased velocities that indicated a less significant narrowing than suspected on a previous MRA, as confirmed by catheter angiography. This result may have prevented an invasive angiogram. In the two patients with 4-mm intracranial aneurysm, one had a stable helical flow pattern with a large jet, while the other had a temporally unstable flow pattern with a more focal jet possibly indicating that the second aneurysm may have a higher likelihood of rupture. Accelerated 4D flow provides time-resolved 3D velocity data in an 8- to 10-min scan. In the stenosis patients, the addition of 4D flow to a traditional MRA adds the velocity data provided from transcranial Doppler ultrasound (TCD) possibly allowing for more accurate grading of stenoses. In the aneurysm patients, visualization of flow patterns may help to provide prognostic information about future risk of rupture.

3D blood flow characteristics in the carotid artery bifurcation assessed by flow-sensitive 4D MRI at 3T

To determine three-dimensional (3D) blood flow patterns in the carotid bifurcation, 10 healthy volunteers and nine patients with internal carotid artery (ICA) stenosis > or =50% were examined by flow-sensitive 4D MRI at 3T. Absolute and mean blood velocities, pulsatility index (PI), and resistance index (RI) were measured in the common carotid arteries (CCAs) by duplex sonography (DS) and compared with flow-sensitive 4D MRI. Furthermore, 3D MRI blood flow patterns in the carotid bifurcation of volunteers and patients before and after recanalization were graded by two independent readers. Blood flow velocities measured by MRI were 31-39% lower than in DS. However, PI and RI differed by only 13-16%. Rating of 3D flow characteristics in the ICA revealed consistent patterns for filling and helical flow in volunteers. In patients with ICA stenosis, 3D blood flow visualization was successfully employed to detect markedly altered filling and helical flow patterns (forward-moving spiral flow) in the ICA bulb and to evaluate the effect of revascularization, which restored filling and helical flow. Our results demonstrate the feasibility of flow-sensitive 4D MRI for the quantification and 3D visualization of physiological and pathological flow patterns in the carotid artery bifurcation.

Optimum fuzzy filters for phase-contrast magnetic resonance imaging segmentation

PURPOSE

To develop and validate a multidimensional segmentation and filtering methodology for accurate blood flow velocity field reconstruction from phase-contrast magnetic resonance imaging (PC MRI).

MATERIALS AND METHODS

The proposed technique consists of two steps: (1) the boundary of the vessel is automatically segmented using the active contour approach; and (2) the noise embedded within the segmented vector field is selectively removed using a novel fuzzy adaptive vector median filtering (FAVMF) technique. This two-step segmentation process was tested and validated on 111 synthetically generated PC MRI slices and on 10 patients with congenital heart disease.

RESULTS

The active contour technique was effective for segmenting blood vessels having a sensitivity and specificity of 93.1% and 92.1% using manual segmentation as a reference standard. FAVMF was the superior technique in filtering out noise vectors, when compared with other commonly used filters in PC MRI (P < 0.05). The peak wall shear rate calculated from the PC MRI data (248 +/- 39 sec(-1)), was significantly decreased to (146 +/- 26 sec(-1)) after the filtering process.

CONCLUSION

The proposed two-step segmentation and filtering methodology is more accurate compared to a single-step segmentation process for post-processing of PC MRI data.

Quantitative 2D and 3D phase contrast MRI: optimized analysis of blood flow and vessel wall parameters

Quantification of CINE phase contrast (PC)-MRI data is a challenging task because of the limited spatiotemporal resolution and signal-to-noise ratio (SNR). The method presented in this work combines B-spline interpolation and Green's theorem to provide optimized quantification of blood flow and vessel wall parameters. The B-spline model provided optimal derivatives of the measured three-directional blood velocities onto the vessel contour, as required for vectorial wall shear stress (WSS) computation. Eight planes distributed along the entire thoracic aorta were evaluated in a 19-volunteer study using both high-spatiotemporal-resolution planar two-dimensional (2D)-CINE-PC ( approximately 1.4 x 1.4 mm(2)/24.4 ms) and lower-resolution 3D-CINE-PC ( approximately 2.8 x 1.6 x 3 mm(3)/48.6 ms) with three-directional velocity encoding. Synthetic data, error propagation, and interindividual, intermodality, and interobserver variability were used to evaluate the reliability and reproducibility of the method. While the impact of MR measurement noise was only minor, the limited resolution of PC-MRI introduced systematic WSS underestimations. In vivo data demonstrated close agreement for flow and WSS between 2D- and 3D-CINE-PC as well as observers, and confirmed the reliability of the method. WSS analysis along the aorta revealed the presence of a circumferential WSS component accounting for 10-20%. Initial results in a patient with atherosclerosis suggest the potential of the method for understanding the formation and progression of cardiovascular diseases.

Complete intracranial arterial and venous blood flow evaluation with 4D flow MR imaging

SUMMARY

Time-resolved, 3D velocity-encoded MR imaging (4D Flow) allows for the acquisition of dynamic, multidirectional data on blood flow and has recently been used for the evaluation of intracranial arterial flow. Using a 3T system with optimization of both temporal resolution and k-space subsampling with a combination of parallel imaging and cut-corner acquisition, we present the clinical assessment of a patient with an arteriovenous malformation by providing complete intracranial arterial and venous coverage in a reasonable scan time.

Improved estimation and visualization of two-dimensional myocardial strain rate using MR velocity mapping

PURPOSE

To estimate regional myocardial strain rate, with reduced sensitivity to noise and velocities outside the region of interest, and provide a visualization of the spatial variation of the obtained tensor field within the myocardium.

MATERIALS AND METHODS

Myocardial velocities were measured using two-dimensional phase contrast velocity mapping. Velocity gradients were estimated using normalized convolution and the calculated 2D strain rate tensor field was visualized using a glyph representation. Validation utilized a numerical phantom with known strain rate distribution. Strain rate glyph visualizations were created for normal myocardium in both systole and diastole and compared to a patient with an anteroseptal infarction.

RESULTS

In the phantom study the strain rate calculated with normalized convolution showed a very good agreement with the analytic solution, while traditional methods for gradient estimation were shown to be sensitive to both noise and surrounding velocity data. Normal myocardium showed a homogenous strain rate distribution, while a heterogeneous strain rate can be clearly seen in the patient data.

CONCLUSION

The proposed approach for quantification and visualization of the regional myocardial strain rate can provide an objective measure of regional myocardial contraction and relaxation that may be valuable for the assessment of myocardial heart disease.