Quantification of internal carotid artery flow with digital subtraction angiography: validation of an optical flow approach with Doppler ultrasound

Simulation-informed learning for time-resolved angiographic contrast agent concentration reconstruction

Three-dimensional Digital Subtraction Angiography (3D-DSA) is a well-established X-ray-based technique for visualizing vascular anatomy. Recently, four-dimensional DSA (4D-DSA) reconstruction algorithms have been developed to enable the visualization of volumetric contrast flow dynamics through time-series of volumes. This reconstruction problem is ill-posed mainly due to vessel overlap in the projection direction and geometric vessel foreshortening, which leads to information loss in the recorded projection images. However, knowledge about the underlying fluid dynamics can be leveraged to constrain the solution space. In our work, we implicitly include this information in a neural network-based model that is trained on a dataset of image-based blood flow simulations. The model predicts the spatially averaged contrast agent concentration for each centerline point of the vasculature over time, lowering the overall computational demand. The trained network enables the reconstruction of relative contrast agent concentrations with a mean absolute error of 0.02±0.02 and a mean absolute percentage error of 5.31±9.25 %. Moreover, the network is robust to varying degrees of vessel overlap and vessel foreshortening. Our approach demonstrates the potential of the integration of machine learning and blood flow simulations in time-resolved angiographic contrast agent concentration reconstruction.

Cerebral Sinus Hemodynamics in Adults Revealed by 4D Flow MRI

BACKGROUND

The hemodynamics of the cerebral sinuses play a vital role in understanding blood flow-related diseases, yet the hemodynamics of the cerebral sinuses in normal adults remains an unresolved issue.

PURPOSE

To evaluate hemodynamics in the cerebral sinus of adults using 4-dimensional flow MRI (4D Flow MRI).

STUDY TYPE

Cross-sectional.

POPULATION

Ninety-nine healthy volunteers (mean age, 42.88 ± 13.16 years old; females/males, 55/44).

FIELD STRENGTH/SEQUENCE

3 T/4D Flow MRI.

ASSESSMENT

The blood flow velocity, average blood flow rate (Q), and vortexes at the superior sagittal sinus (SSS), straight sinus (STS), transverse sinus, sigmoid sinus, and jugular bulb of each volunteer were evaluated by two independent neuroradiologists. The relationship between the total cerebral Q and sex and age was also assessed. Twelve volunteers underwent two scans within a month.

STATISTICAL TESTS

The intraclass correlation coefficient (ICC) evaluated the inter-observer agreement. Blood flow parameters among volunteers were compared by the independent-sample t-test or Mann-Whitney U test. The multiple linear regression equation was used to evaluate the relationship between total cerebral Q and age and sex. P < 0.05 indicated statistical significance.

RESULTS

The test-retest and interobserver reliability of average velocity and Q were moderate to high (ICC: 0.54-0.99). Cerebral sinus velocity varied by segment and cardiac cycle. The SSS's velocity and Q increased downstream and Q near torcular herophili was 3.5 times that through the STS. The total cerebral Q decreased by 0.06 mL/s per year (β = -0.06 ± 0.013) and was sex-independent within the group. Vortexes were found in 12.12%, 8.9%, and 59.8% of torcular herophili, transverse-sigmoid junction, and jugular bulb, respectively, and were related to higher upstream flow.

DATA CONCLUSION

Cerebral sinuses could be measured visually and quantitatively in vivo by 4D Flow MRI, providing a basis for future research on pulsating tinnitus, multiple sclerosis, and other related diseases.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Application of Digital Subtraction Angiography in Predicting the Outcomes of Intraoperative Hemorrhage of Juvenile Nasopharyngeal Angiofibroma

OBJECTIVE

Juvenile nasopharyngeal angiofibroma (JNA) is a very rare hemorrhagic vascular tumor that predominantly affects adolescent boys. The tumor is relatively large when detected, and the risk of intraoperative bleeding is high. We aimed to examine factors associated with intraoperative blood loss in JNA surgery.

METHODS

Thirteen patients with JNA who underwent surgery at the Jikei University Hospital between 2009 and 2020 were retrospectively reviewed, and factors associated with blood loss were examined by single regression analysis.

RESULTS

The mean age was 20.8 ± 7.7 years. Preoperative angiographic images were evaluated in 9 of the 13 cases. The 6 patients with the largest bleeding volumes, all had residual nutrient vessels from the internal carotid artery (ICA), with an average number of 2.5 vessels. The mean blood loss of patients with residual nutrient vessels from the ICA was 3037 ± 2568 mL. Single regression analysis of bleeding volume against the number of remaining nutrient vessels from the ICA and the total peak contrast density of nutrient vessels (Cmax) standardized by region of interest showed that the coefficient was positive (P < 0.05 for both), confirming a significant correlation between the 2, respectively.

CONCLUSIONS

The amount of bleeding significantly correlated with the number of remaining nutrient vessels from the ICA after preoperative embolization and with the total Cmax/region of interest. The ability to predict the amount of preoperative blood loss using this study will facilitate proposals for external incisions in patients with JNA.

Angiographic velocimetry analysis using contrast dilution gradient method with a 1000 frames per second photon-counting detector

Purpose

Contrast dilution gradient (CDG) analysis is a quantitative method allowing blood velocity estimation using angiographic acquisitions. Currently, CDG is restricted to peripheral vasculature due to the suboptimal temporal resolution of current imaging systems. We investigate extension of CDG methods to the flow conditions of proximal vasculature using 1000 frames per second (fps) high-speed angiographic (HSA) imaging.

Approach

We performed in-vitro HSA acquisitions using the XC-Actaeon detector and 3D-printed patient-specific phantoms. The CDG approach was used for blood velocity estimation expressed as the ratio of temporal and spatial contrast gradients. The gradients were extracted from 2D contrast intensity maps synthesized by plotting intensity profiles along the arterial centerline at each frame. In-vitro results obtained at various frame rates via temporal binning of 1000 fps data were retrospectively compared to computational fluid dynamics (CFD) velocimetry. Full-vessel velocity distributions were estimated at 1000 fps via parallel line expansion of the arterial centerline analysis.

Results

Using HSA, the CDG method displayed agreement with CFD at or above 250 fps [mean-absolute error (MAE): 2.6 ± 6.3    cm / s , p = 0.05 ]. Relative velocity distributions correlated well with CFD at 1000 fps with universal underapproximation due to effects of pulsatile contrast injection (MAE: 4.3 cm/s).

Conclusions

Using 1000 fps HSA, CDG-based extraction of velocities across large arteries is possible. The method is sensitive to noise; however, image processing techniques and a contrast injection, which adequately fills the vessel assist algorithm accuracy. The CDG method provides high resolution quantitative information for rapidly transient flow patterns observed in arterial circulation.

A novel angiographic method to estimate arterial blood flow rates using contrast reflux: Effect of injection parameters

BACKGROUND

Contrast reflux, which is the retrograde movement of contrast against flow direction, is commonly observed during angiography. Despite a vast body of literature on angiography, the hemodynamic factors affecting contrast reflux have not been studied. Numerous methods have been developed to extract flow from angiography, but the reliability of these methods is not yet sufficient to be of routine clinical use.

PURPOSE

To evaluate the effect of baseline blood flow rates and injection conditions on the extent of contrast reflux. To estimate arterial flow rates based on measurement of contrast reflux length.

MATERIALS AND METHODS

Iodinated contrast was injected into an idealized tube as well as a physiologically accurate model of the cervico-cerebral vasculature. A total of 194 high-speed angiograms were acquired under varying "blood" flow rates and injection conditions (catheter size, injection rate, and injection time). The length of contrast reflux was compared to the input variables and to dimensionless fluid dynamics parameters at the catheter-tip. Arterial blood flow rates were estimated using contrast reflux length as well as a traditional transit-time method and compared to measured flow rates.

RESULTS

Contrast reflux lengths were significantly affected by contrast injection rate (p < 0.0001), baseline blood flow rate (p = 0.0004), and catheter size (p = 0.04), but not by contrast injection time (p = 0.4). Reflux lengths were found to be correlated to dimensionless fluid dynamics parameters by an exponential function (R²  = 0.6-0.99). When considering the entire dataset in unison, flow estimation errors with the reflux-length method (39% ± 33%) were significantly higher (p = 0.003) than the transit-time method (33% ± 36%). However, when subgrouped by catheter, the error with the reflux-length method was substantially reduced and was significantly lower (14% ± 14%, p < 0.0001) than the transit-time method.

CONCLUSION

Results show correlations between contrast reflux length and baseline hemodynamic parameters that have not been reported previously. Clinically relevant blood flow rate estimation is feasible by simple measurement of reflux length. In vivo and clinical studies are required to confirm these correlations and to refine the methodology of estimating blood flow by reflux.

Extraction of patient-specific boundary conditions from 4D-DSA and their influence on CFD simulations of cerebral aneurysms

We developed a new technique for extracting patient-specific inflow conditions, such as the pulse cycle duration and blood flow velocity, from four-dimensional digital subtraction angiography images and experimentally examined its validity. The maximum error between the values extracted by the technique and measured values was 14.3%. We performed blood flow simulations and calculated representative haemodynamic parameters. The maximum differences between the parameters obtained using general and patient-specific inflow conditions were approximately 400%, 150%, and 50% for the velocity, normalised wall shear stress, and pressure loss coefficient, respectively. These results indicate that patient-specific conditions are critical for accurately reproducing aneurysmal haemodynamics.

2D vessel contrast dilution gradient (CDG) analysis using 1000 fps high speed angiography (HSA) for velocity distribution estimation

PURPOSE

Contrast dilution gradient (CDG) analysis is a technique used to extract velocimetric 2D information from digitally subtracted angiographic (DSA) acquisitions. This information may then be used by clinicians to quantitatively assess the effects of endovascular treatment on flow conditions surrounding pathologies of interest. The method assumes negligible diffusion conditions, making 1000 fps high speed angiography (HSA), in which diffusion between 1 ms frames may be neglected, a strong candidate for velocimetric analysis using CDG. Previous studies have demonstrated the success of CDG analysis in obtaining velocimetric one-dimensional data at the arterial centerline of simple vasculature. This study seeks to resolve velocity distributions across the entire vessel using 2D-CDG analysis with HSA acquisitions.

MATERIALS AND METHODS

HSA acquisitions for this study were obtained in vitro with a benchtop flow loop at 1000 fps using the XC-Actaeon (Direct Conversion Inc.) photon counting detector. 2D-CDG analyses were compared with computational fluid dynamics (CFD) via automatic co-registration of the results from each velocimetry method. This comparison was performed using mean absolute error between pixel values in each method (after temporal averaging).

RESULTS

CDG velocity magnitudes were slightly under approximated relative to CFD results (mean velocity: 27 cm/s, mean absolute error: 4.3 cm/s) as a result of incomplete contrast filling. Relative 2D spatial velocity distributions in CDG analysis agreed well with CFD distributions qualitatively.

CONCLUSIONS

CDG may be used to obtain velocity distributions in and surrounding vascular pathologies provided diffusion is negligible relative to convection in the flow, given a continuous gradient of contrast.

Comparison of algorithms for estimating blood flow velocities in cerebral arteries based on the transport information of contrast agent: An in silico study

While many algorithms have been proposed to estimate blood flow velocities based on the transport information of contrast agent acquired by digital subtraction angiography (DSA), most relevant studies focused on a single vessel, leaving a question open as to whether the algorithms would be suitable for estimating blood flow velocities in arterial systems with complex topological structures. In this study, a one-dimensional (1-D) modeling method was developed to simulate the transport of contrast agent in cerebral arterial networks with various anatomical variations or having occlusive disease, thereby generating an in silico database for examining the accuracies of some typical algorithms (i.e., time-of-center of gravity (TCG), shifted least-squares (SLS), and cross correlation (CC) algorithms) that estimate blood flow velocity based on the concentration-time curves (CTCs) of contrast agent. The results showed that the TCG algorithm had the best performance in estimating blood flow velocities in most cerebral arteries, with the accuracy being only mildly affected by anatomical variations of the cerebral arterial network. Nevertheless, the presence of a stenosis of moderate to high severity in the internal carotid artery could considerably impair the accuracy of the TCG algorithm in estimating blood flow velocities in some cerebral arteries where the transport of contrast agent was disturbed by strong collateral flows. In summary, the study suggests that the TCG algorithm may offer a promising means for estimating blood flow velocities based on CTCs of contrast agent monitored in cerebral arteries, provided that the shapes of CTCs are not highly distorted by collateral flows.

Mathematical Models for Blood Flow Quantification in Dialysis Access Using Angiography: A Comparative Study

Blood flow rate in dialysis (vascular) access is the key parameter to examine patency and to evaluate the outcomes of various endovascular interve7ntions. While angiography is extensively used for dialysis access–salvage procedures, to date, there is no image-based blood flow measurement application commercially available in the angiography suite. We aim to calculate the blood flow rate in the dialysis access based on cine-angiographic and fluoroscopic image sequences. In this study, we discuss image-based methods to quantify access blood flow in a flow phantom model. Digital subtraction angiography (DSA) and fluoroscopy were used to acquire images at various sampling rates (DSA—3 and 6 frames/s, fluoroscopy—4 and 10 pulses/s). Flow rates were computed based on two bolus tracking algorithms, peak-to-peak and cross-correlation, and modeled with three curve-fitting functions, gamma variate, lagged normal, and polynomial, to correct errors with transit time measurement. Dye propagation distance and the cross-sectional area were calculated by analyzing the contrast enhancement in the vessel. The calculated flow rates were correlated versus an in-line flow sensor measurement. The cross-correlation algorithm with gamma-variate curve fitting had the best accuracy and least variability in both imaging modes. The absolute percent error (mean ± SEM) of flow quantification in the DSA mode at 6 frames/s was 21.4 ± 1.9%, and in the fluoroscopic mode at 10 pulses/s was 37.4 ± 3.6%. The radiation dose varied linearly with the sampling rate in both imaging modes and was substantially low to invoke any tissue reactions or stochastic effects. The cross-correlation algorithm and gamma-variate curve fitting for DSA acquisition at 6 frames/s had the best correlation with the flow sensor measurements. These findings will be helpful to develop a software-based vascular access flow measurement tool for the angiography suite and to optimize the imaging protocol amenable for computational flow applications.

Endovascular Treatment Versus Best Medical Therapy in Acute Ischemic Stroke Patients with Mild Symptoms

BACKGROUND

The benefit of endovascular treatment (EVT) for acute ischemic stroke patients with mild deficits is unknown. We sought to evaluate the natural history of patients with a low National Institute of Health Stroke Score (NIHSS) and an intracranial occlusion.

METHODS

We included patients with a computed tomography angiogram-proven intracranial arterial occlusion who presented within 24 hours of symptom onset with an NIHSS of ≤6. We compared outcomes of patients who were treated with EVT and those who were not by performing propensity score-matched analysis. Primary outcome was modified Rankin score (mRS) at 90 days.

RESULTS

A total of 66 patients were included: 38 were men (57.6%) with a median age of 69 (interquartile range [IQR], 57-79.5) years. Median NIHSS was 3 (IQR, 2-5). Median time from symptom onset to presentation was 2.87 hours (IQR, 1.3-5.9). Forty of the total cohort underwent best medical therapy alone (60.6%), whereas 26 underwent EVT (39.4%). Nineteen of the 26 patients who underwent EVT had a good clinical outcome (mRS ≤2) (73.1%), compared with 29 of 40 best medical therapy patients (72.5%) (odds ratio, 0.833 with 95% confidence interval, 0.263-2.631; P = 0.755). Following propensity score adjustment there was a tendency toward lower mRS following EVT (P = 0.051).

CONCLUSIONS

Despite the higher number of proximal occlusions in the EVT group, overall outcomes were similar, with >70% of patients in each cohort having a good outcome at 90 days.

Blood flow quantification in dialysis access using digital subtraction angiography: A retrospective study

BACKGROUND AND OBJECTIVE

Vascular access is the "lifeline" of end-stage renal disease patients, which is surgically constructed to remove blood-waste and return artificially filtered blood into circulation. The arteriovenous shunting causes an abrupt change in blood flow and results in increased fluidic stress, which predisposes to access stenosis and thrombosis. While access flow is crucial to evaluate interventional endpoint, application to measure flow using digital angiogram is not yet available. The goal of this study was to determine the feasibility of flow quantification in dialysis access using a software tool and to guide the design of an imaging protocol.

METHODS

173 digital subtraction angiographic (DSA) images were retrospectively analyzed to evaluate access flow in a custom-programming environment. Four bolus transit time algorithms and a distance calculation method were assessed for flow computation. Gamma variate function was applied to remove secondary flow and intensity outliers in the bolus time-intensity curves and evaluated for enhancement in computational accuracy. The percent deviations of flow rates computed from dilution of iodinated radio-contrast material were compared with in situ catheter-based flow measurement.

RESULTS

Among the implemented bolus transit time algorithms, quantification error (mean ± standard error) of cross-correlation algorithm without and with gamma variate curve fitting was 35 ± 1% and 22 ± 1%, respectively. All other algorithms had quantification error >27%. The bias and limits of agreement of the cross-correlation algorithm with gamma variate curve fit was -94 ml/min and [-353, 165] mL/min, respectively.

CONCLUSIONS

The cross-correlation algorithm with gamma variate curve fit had the best accuracy and reproducibility for image-based blood flow computation. To further enhance accuracy, images may need to be acquired with a dedicated injection protocol with predetermined parameters such as the duration, rate and mode of bolus injection, and the acquisition frame rate.

Optimizing the Quality of 4D-DSA Temporal Information

BACKGROUND AND PURPOSE

Quantification of blood flow using a 4D-DSA would be useful in the diagnosis and treatment of cerebrovascular diseases. A protocol optimizing identification of density variations in the time-density curves of a 4D-DSA has not been defined. Our purpose was to determine the contrast injection protocol most likely to result in the optimal pulsatility signal strength.

MATERIALS AND METHODS

Two 3D-printed patient-specific models were used and connected to a pulsatile pump and flow system, which delivered 250-260 mL/min to the model. Contrast medium (Isovue, 370 mg I/mL, 75% dilution) was injected through a 6F catheter positioned upstream from the inlet of the model. 4D-DSA acquisitions were performed for the following injection rates: 1.5, 2.0, 2.5, 3.0 and 3.5 mL/s for 8 seconds. To determine pulsatility, we analyzed the time-density curve at the inlets using the oscillation amplitude and a previously described numeric metric, the sideband ratio. Vascular geometry from 4D-DSA reconstructions was compared with ground truth and micro-CT measurements of the model. Dimensionless numbers that characterize hemodynamics, Reynolds and Craya-Curtet, were calculated for each injection rate.

RESULTS

The strongest pulsatility signal occurred with the 2.5 mL/s injections. The largest oscillation amplitudes were found with 2.0- and 2.5-mL/s injections. Geometric accuracy was best preserved with injection rates of >1.5 mL/s.

CONCLUSIONS

An injection rate of 2.5 mL/s provided the strongest pulsatility signal in the 4D-DSA time-density curve. Geometric accuracy was best preserved with injection rates above 1.5 mL/s. These results may be useful in future in vivo studies of blood flow quantification.

Fluoroscopic angiography quantifies delay in cerebral circulation time and requires less radiation in carotid stenosis patients: A pilot study

BACKGROUND

Quantitative digital subtraction angiography (DSA) facilitates in-room assessment of flow changes in various cerebrovascular diseases and improves patient safety. The purpose of this study was to compare the diagnostic accuracy of quantitative fluoroscopic angiography (FA) and DSA.

METHODS

Twenty-two patients with >70% carotid stenosis according to NASCET criteria were prospectively included in the study. All patients received DSA and FA (ArtisZee, Siemens Healthcare, Forchheim, Germany) before and after carotid stenting in the same angiosuite. The regions of interest (ROIs) included the extracranial internal carotid artery (eICA), first segment of the middle cerebral artery (MCA1), and sigmoid sinus in the anterior-posterior view; cavernous portion of the ICA (cICA), parietal vein, and jugular vein in the lateral views. The time-to-peak (TTP) for all ROIs and cerebral circulation time (CCT) were measured from FA and DSA scans. TTP, CCT, and radiation doses from DSA were compared with those from FA.

RESULTS

The mean age of the patients were 69 ± 9.5 years old. The average stenosis was 89.7% ± 7.8% before stenting and 31% ± 3.6% after stenting. No patient suffered from periprocedural stroke. The intermethod correlation for TTP for all ROIs except the eICA and cICA ranged from 0.46 to 0.65 before stenting and 0.57 to 0.73 after stenting, and that for CCT was 0.65 before stenting and 0.57 after stenting. The radiation doses were significantly lower for FA than for DSA regardless of views or periprocedural timing (p < 0.001).

CONCLUSION

Stenosis facilitated the creation of a bolus by manual injection and therefore increased the accuracy of cerebral flow quantification in FA. Cerebral hemodynamic assessment by FA is quicker and associated with less radiation.

Flow diversion for treatment of intracranial aneurysms: Mechanism and implications

Flow diverters are new generation stents that have recently garnered a large amount of interest for use in treatment of intracranial aneurysms. Flow diverters reduce blood flow into the aneurysm, with redirection along the path of the parent vessel. Flow stagnation into the aneurysm and neck coverage with subsequent endothelialization are the important synergistic mechanisms by which the therapy acts. Several studies have examined the mechanisms by which flow diverters subsequently lead to aneurysm occlusion. This review aims to provide a general overview of the flow diverters and their mechanism of action and potential implications. ANN NEUROL 2019;85:793-800.

Velocimetry based on dye visualization for a pulsatile tubing flow measurement

A dye visualization experiment was designed and operated as an analogy of x-ray angiographic imaging. The velocimetry technique based on dye visualization for a pulsatile tubing flow measurement is introduced. The optical flow method was utilized to recover the velocity field from the visualization images. For the same flow condition, a digital particle image velocimetry (PIV) system was also employed to measure the same flow field in the middle plane. The purpose of this study was to determine the accuracy of the velocity field estimation from the transmittance-based two-dimensional projection image of the three-dimensional volumetric flow field by dye visualization in comparison with the PIV measurement results. Compared to the PIV results in the middle plane, the averaged velocity magnitude from the dye visualization measurement was underestimated by about 16%-24% in the central region and by about 29%-43% in the outer region across the tube at two time instants of the cyclic pulsatile flow.

A multiscale computational modeling for cerebral blood flow with aneurysms and/or stenoses

A 1-dimensional (1D)-3-dimensional (3D) multiscale model for the human vascular network was proposed by combining a low-fidelity 1D modeling of blood circulation to account for the global hemodynamics with a detailed 3D simulation of a zonal vascular segment. The coupling approach involves a direct exchange of flow and pressure information at interfaces between the 1D and 3D models and thus enables patient-specific morphological models to be inserted into flow network with minimum computational efforts. The proposed method was validated with good agreements against 3 simplified test cases where experimental data and/or full 3D numerical solution were available. The application of the method in aneurysm and stenosis studies indicated that the deformation of the geometry caused by the diseases may change local pressure loss and as a consequence lead to an alteration of flow rate to the vessel segment.

Quantification of Blood Velocity with 4D Digital Subtraction Angiography Using the Shifted Least-Squares Method

BACKGROUND AND PURPOSE

4D-DSA provides time-resolved 3D-DSA volumes with high temporal and spatial resolutions. The purpose of this study is to investigate a shifted least squares method to estimate the blood velocity from the 4D DSA images. Quantitative validation was performed using a flow phantom with an ultrasonic flow probe as ground truth. Quantification of blood velocity in human internal carotid arteries was compared with measurements generated from 3D phase-contrast MR imaging.

MATERIALS AND METHODS

The centerlines of selected vascular segments and the time concentration curves of each voxel along the centerlines were determined from the 4D-DSA dataset. The temporal shift required to achieve a minimum difference between any point and other points along the centerline of a segment was calculated. The temporal shift as a function of centerline point position was fit to a straight line to generate the velocity. The proposed shifted least-squares method was first validated using a flow phantom study. Blood velocities were also estimated in the 14 ICAs of human subjects who had both 4D-DSA and phase-contrast MR imaging studies. Linear regression and correlation analysis were performed on both the phantom study and clinical study, respectively.

RESULTS

Mean velocities of the flow phantom calculated from 4D-DSA matched very well with ultrasonic flow probe measurements with 11% relative root mean square error. Mean blood velocities of ICAs calculated from 4D-DSA correlated well with phase-contrast MR imaging measurements with Pearson correlation coefficient r = 0.835.

CONCLUSIONS

The availability of 4D-DSA provides the opportunity to use the shifted least-squares method to estimate velocity in vessels within a 3D volume.

Measuring blood velocity using 4D-DSA: A feasibility study

PURPOSE

Four-dimensional (4D) DSA reconstruction provides three-dimensional (3D) time-resolved visualization of contrast bolus passage through arterial vasculature in the interventional setting. The purpose of this study was to evaluate the feasibility of using these data in measuring blood velocity and flow.

METHODS

The pulsatile signals in the time concentration curves (TCCs) measured at different points along a vessel are markers of the movement of a contrast bolus and thus of blood flow. When combined with the spatial content, that is, geometry of the vasculature, this information then provides the data required to determine blood velocity. A Fourier-based algorithm was used to identify and follow the pulsatility signal. A Side Band Ratio (SBR) metric was used to reduce uncertainty in identifying the pulsatility in regions where the signal was weak. We tested this method using 4D-DSA reconstructions from vascular phantoms as well as from human studies.

RESULTS

In five studies using 3D printed patient-specific cerebrovascular phantoms, velocities calculated from the 4D-DSAs were found to be within 10% of velocities measured with a flow meter. Calculated velocity and flow values from three human studies were within the range of those reported in the literature.

CONCLUSIONS

4D-DSA provides temporal and spatial information about blood flow and vascular geometry. This information is obtained using conventional rotational angiographic systems. In this small feasibility study, these data allowed calculations of velocity values that correlated well with measured values. The availability of velocity and blood flow information in the interventional setting would support a more quantitative approach to diagnosis, treatment planning and post-treatment evaluations of a variety of cerebrovascular diseases.

Automatic flow analysis of digital subtraction angiography using independent component analysis in patients with carotid stenosis

PURPOSE

Current time-density curve analysis of digital subtraction angiography (DSA) provides intravascular flow information but requires manual vasculature selection. We developed an angiographic marker that represents cerebral perfusion by using automatic independent component analysis.

MATERIALS AND METHODS

We retrospectively analyzed the data of 44 patients with unilateral carotid stenosis higher than 70% according to North American Symptomatic Carotid Endarterectomy Trial criteria. For all patients, magnetic resonance perfusion (MRP) was performed one day before DSA. Fixed contrast injection protocols and DSA acquisition parameters were used before stenting. The cerebral circulation time (CCT) was defined as the difference in the time to peak between the parietal vein and cavernous internal carotid artery in a lateral angiogram. Both anterior-posterior and lateral DSA views were processed using independent component analysis, and the capillary angiogram was extracted automatically. The full width at half maximum of the time-density curve in the capillary phase in the anterior-posterior and lateral DSA views was defined as the angiographic mean transient time (aMTT; i.e., aMTTAP and aMTTLat). The correlations between the degree of stenosis, CCT, aMTTAP and aMTTLat, and MRP parameters were evaluated.

RESULTS

The degree of stenosis showed no correlation with CCT, aMTTAP, aMTTLat, or any MRP parameter. CCT showed a strong correlation with aMTTAP (r = 0.67) and aMTTLat (r = 0.72). Among the MRP parameters, CCT showed only a moderate correlation with MTT (r = 0.67) and Tmax (r = 0.40). aMTTAP showed a moderate correlation with Tmax (r = 0.42) and a strong correlation with MTT (r = 0.77). aMTTLat also showed similar correlations with Tmax (r = 0.59) and MTT (r = 0.73).

CONCLUSION

Apart from vascular anatomy, aMTT estimates brain parenchyma hemodynamics from DSA and is concordant with MRP. This process is completely automatic and provides immediate measurement of quantitative peritherapeutic brain parenchyma changes during stenting.

In vitro digital subtraction angiographic evaluation of flow diverters in a patient-specific aneurysm

Background The importance of both porosity and pore density of a flow diverter is well recognized in treatment of intracranial aneurysms; however, understanding of the effect of individual wire (wire number and size) is critical in improving device design and use. Methods A total of 10 multi-layered flow diverters with different wire numbers (32, 48, 56, and 72) and sizes (30, 35, and 40 µm) were implanted into identical patient-specific middle cerebral artery aneurysm models. Digital subtraction angiography was acquired at 30 f/s and X-ray signals at three selected regions of interest were compared to determine the amount of intra-aneurysmal flow. Results Flow reduction ranged from 19% for a high porosity (82%) and low pore density (5 pores/mm²) to nearly 80% for a low porosity (49%) and high pore density (36 pores/mm²). An increase in the wire number from 32 to 72 lowers intra-aneurysmal flow and redirects the flow jet; however, the effect of wire size is not observed. Conclusions In our in vitro angiographic study, flow jet is influenced by the wire number in a device qualitatively; quantitatively, intra-aneurysmal flow is affected by both the porosity and pore density. A 2.5 mm device performs better in flow diversion of a middle cerebral artery (MCA) aneurysm than a 3 mm device with the same wire size and wire number, but thicker wires do not lead to better flow diversion.

Divergence Compensatory Optical Flow Method for Blood Velocimetry

Detailed blood velocity map in the vascular system can be obtained by applying the optical flow method (OFM) in processing fluoroscopic digital subtracted catheter angiographic images; however, there are still challenges with the accuracy of this method. In the present study, a divergence compensatory optical flow method (DC-OFM), in which a nonzero divergence of velocity is assumed due to the finite resolution of the image, was explored and applied to the digital subtraction angiography (DSA) images of blood flow. The objective of this study is to examine the applicability and evaluate the accuracy of DC-OFM in assessing the blood flow velocity in vessels. First, an Oseen vortex flow was simulated on the standard particle image to generate an image pair. Then, the DC-OFM was applied on the particle image pair to recover the velocity field for validation. Second, DSA images of intracranial arteries were used to examine the accuracy of the current method. For each set of images, the first image is the in vivo DSA image, and the second image is generated by superimposing a given flow field. The recovered velocity map by DC-OFM agrees well with the exact velocity for both the particle images and the angiographic images. In comparison with the traditional OFM, the present method can provide more accurate velocity estimation. The accuracy of the velocity estimation can also be improved by implementing preprocess techniques including image intensification, Gaussian filtering, and "image-shift."

4D Flow MRI in Neuroradiology: Techniques and Applications

Assessment of the intracranial flow is important for the understanding and management of cerebral vascular diseases. From brain aneurysms and arteriovenous malformations lesions to intracranial and cervical stenosis, the appraisal of the blood flow can be crucial and influence positively on patients' management. The determination of the intracranial hemodynamics and the collateral pattern seems to play to a major role in the management of these lesions. 4D flow magnetic resonance imaging is a noninvasive phase contrast derived method that has been developed and applied in neurovascular diseases. It has a great potential if followed by further technical improvements and comprehensive and systematic clinical studies.

Does Arterial Flow Rate Affect the Assessment of Flow-Diverter Stent Performance?

BACKGROUND AND PURPOSE

Our aim was to assess the performance of flow-diverter stents. The pre- and end-of-treatment angiographies are commonly compared. However, the arterial flow rate may change between acquisitions; therefore, a better understanding of its influence on the local intra-aneurysmal hemodynamics before and after flow-diverter stent use is required.

MATERIALS AND METHODS

Twenty-five image-based aneurysm models extracted from 3D rotational angiograms were conditioned for computational fluid dynamics simulations. Pulsatile simulations were performed at different arterial flow rates, covering a wide possible range of physiologic flows among 1-5 mL/s. The effect of flow-diverter stents on intra-aneurysmal hemodynamics was numerically simulated with a porous medium model. Spatiotemporal-averaged intra-aneurysmal flow velocity and flow rate were calculated for each case to quantify the hemodynamics after treatment. The short-term flow-diverter stent performance was characterized by the relative velocity reduction inside the aneurysm.

RESULTS

Spatiotemporal-averaged intra-aneurysmal flow velocity before and after flow-diverter stent use is linearly proportional to the mean arterial flow rate (minimum R² > 0.983 of the linear regression models for untreated and stented models). Relative velocity reduction asymptotically decreases with increasing mean arterial flow rate. When the most probable range of arterial flow rate was considered (3-5 mL/s), instead of the wide possible flow range, the mean SD of relative velocity reduction was reduced from 3.6% to 0.48%.

CONCLUSIONS

Both intra-aneurysmal aneurysm velocity and flow-diverter stent performance depend on the arterial flow rate. The performance could be considered independent of the arterial flow rates within the most probable range of physiologic flows.

Normal ranges and test-retest reproducibility of flow and velocity parameters in intracranial arteries measured with phase-contrast magnetic resonance imaging

INTRODUCTION

The purpose of the present study was to investigate normal ranges and test-retest reproducibility of phase-contrast MRI (PC-MRI)-measured flow and velocity parameters in intracranial arteries.

METHODS

Highest flow (HF), lowest flow (LF), peak systolic velocity (PSV), and end diastolic velocity (EDV) were measured at two dates in the anterior (ACA), middle (MCA), and posterior (PCA) cerebral arteries of 30 healthy volunteers using two-dimensional PC-MRI at 3 T. Least detectable difference (LDD) was calculated.

RESULTS

In the left ACA, HF was (mean (range, LDD)) 126 ml/min (36-312, 59 %), LF 61 ml/min (0-156, 101 %), PSV 64 cm/s (32-141, 67 %), and EDV 35 cm/s (18-55, 42 %); in the right ACA, HF was 154 ml/min (42-246, 49 %), LF 77 ml/min (0-156, 131 %), PSV 75 cm/s (26-161, 82 %), and EDV 39 cm/s (7-59, 67 %). In the left MCA, HF was 235 ml/min (126-372, 35 %), LF 116 ml/min (42-186, 48 %), PSV 90 cm/s (55-183, 39 %), and EDV 46 cm/s (20-66, 28 %); in the right MCA, HF was 238 ml/min (162-342, 44 %), LF 120 ml/min (72-216, 48 %), PSV 88 cm/s (55-141, 35 %), and EDV 45 cm/s (26-67, 23 %). In the left PCA, HF was 108 ml/min (42-168, 54 %), LF 53 ml/min (18-108, 64 %), PSV 50 cm/s (24-77, 63 %), and EDV 28 cm/s (14-40, 45 %); in the right PCA, HF was 98 ml/min (30-162, 49 %), LF 49 ml/min (12-84, 55 %), PSV 47 cm/s (27-88, 59 %), and EDV 27 cm/s (16-41, 45 %).

CONCLUSION

PC-MRI-measured flow and velocity parameters in the main intracranial arteries have large normal ranges. Reproducibility is highest in MCA.

Intra-aneurysmal hemodynamics: evaluation of pCONus and pCANvas bifurcation aneurysm devices using DSA optical flow imaging

Background

Implantation of self-expanding stents from the parent artery into the sac of a bifurcation aneurysm is regularly used to facilitate endovascular coil occlusion with the so-called waffle cone technique (WCT). Self-expanding aneurysm bridging stents like Solitaire AB, can be used; however, bifurcation devices like pCONus and pCANvas are especially designed for WCT. These devices provide additional support for coil implantation owing to intraluminal nylon fibers (pCONus) or membranes (pCANvas) covering the intracranial aneurysm neck.

Objective

Assessment of the intra-aneurysmal hemodynamic impact of these three devices: a regular intracranial stent (Solitaire AB) and two bifurcation devices (pCONus and pCANvas).

Material and methods

An in vitro experiment was set up using a silicone model of a basilar tip aneurysm filled with blood mimicking fluid under a pulsatile circulation. Solitaire AB, pCONus, and pCANvas were successively implanted in the model for hemodynamic evaluation. High frame rate DSA series were acquired under various conditions. Intra-aneurysmal flow changes, including mean aneurysm flow amplitude ratio (R), were subsequently assessed by the optical flow method, measuring the detector velocity field before and after device implantations.

Results

pCONus and Solitaire minimally reduced the intra-aneurysmal flow (R=0.96, p=0.17 and R=0.91, p=0.01, respectively), whereas pCANvas strongly diminished the intra-aneurysmal flow (R=0.41, p=5×10⁻¹²).

Conclusions

Waffle cone deployment of stents and technique-specific devices had no undesirable effect on the intra-aneurysmal flow. In particular, no increased flow was redirected into the aneurysm sac. The intraluminal membrane of the pCANvas strongly reduced the intra-aneurysmal flow, potentially preventing recanalization problems.

Contrast-enhanced angiographic cone-beam computed tomography without pre-diluted contrast medium

INTRODUCTION

Contrast-enhanced cone-beam computed tomography (CBCT) has been introduced and accepted as a useful technique to evaluate delicate vascular anatomy and neurovascular stents. Current protocol for CBCT requires quantitative dilution of contrast medium to obtain adequate quality images. Here, we introduce simple methods to obtain contrast-enhanced CBCT without quantitative contrast dilution.

METHODS

A simple experiment was performed to estimate the change in flow rate in the internal carotid artery during the procedure. Transcranial doppler (TCD) was used to evaluate the velocity change before and after catheterization and fluid infusion. In addition, 0.3 cm(3)/s (n = 3) and 0.2 cm(3)/s (n = 7) contrast infusions were injected and followed by saline flushes using a 300 mmHg pressure bag to evaluate neurovascular stent and host arteries.

RESULTS

Flow velocities changed -15 ± 6.8 % and +17 ± 5.5 % from baseline during catheterization and guiding catheter flushing with a 300 mmHg pressure bag, respectively. Evaluation of the stents and vascular structure was feasible using this technique in all patients. Quality assessment showed that the 0.2 cm(3)/s contrast infusion protocol was better for evaluating the stent and host artery.

CONCLUSION

Contrast-enhanced CBCT can be performed without quantitative contrast dilution. Adequate contrast dilution can be achieved with a small saline flush and normal blood flow.

Changes of time-attenuation curve blood flow parameters in patients with and without carotid stenosis

BACKGROUND AND PURPOSE

From the time-attenuation curves of DSA flow parameters, maximal intensity, maximal slope, and full width at half maximum of selected vascular points are defined. The study explores the reliability of defining the flow parameters by the time-attenuation curves of DSA.

MATERIALS AND METHODS

Seventy patients with unilateral carotid artery stenosis (group A) and 56 healthy controls (group B) were retrospectively enrolled. Fixed contrast injection protocols and DSA acquisition parameters were used with all patients. The M1, sigmoid sinus, and internal jugular vein on anteroposterior view DSA and the M2, parietal vein, and superior sagittal sinus on lateral view DSA were chosen as ROI targets for measuring flow parameters. The difference of time of maximal intensity between 2 target points was defined as the circulation time between the target points.

RESULTS

The maximal intensity difference of 2 selected points from the ICA to the M1, sigmoid sinus, internal jugular vein, M2, parietal vein, and superior sagittal sinus was significantly longer in group A than in group B. The maximum slope of M1, M2, and the superior sagittal sinus was significantly lower in group A than in group B. The full width at half maximum of M1 and M2 was significantly larger in group A than in group B. The maximal slope of M1 demonstrated the best diagnostic performance.

CONCLUSIONS

The maximal intensity difference of 2 selected points derived from DSA can be used as a definitive alternative flow parameter for intracranial circulation time measurement. Maximal slope and full width at half maximum complement the maximal intensity difference of 2 selected points in defining flow characteristics of healthy subjects and patients with carotid stenosis.

Intra-aneurysmal flow patterns: illustrative comparison among digital subtraction angiography, optical flow, and computational fluid dynamics

BACKGROUND AND PURPOSE

Digital subtraction angiography is the gold standard vascular imaging and it is used for all endovascular treatment of intracranial anerysms. Optical flow imaging has been described as a potential method to evaluate cerebral hemodynamics through DSA. In this study, we aimed to compare the flow patterns measured during angiography, by using an optical flow method, with those measured by using computational fluid dynamics in intracranial aneurysms.

MATERIALS AND METHODS

A consecutive series of 21 patients harboring unruptured saccular intracranial aneurysms who underwent diagnostic angiography before treatment was considered. High-frame-rate digital subtraction angiography was performed to obtain an intra-aneurysmal velocity field by following the cardiac-modulated contrast wave through the vascular structures by using optical flow principles. Additionally, computational fluid dynamics modeling was performed for every case by using patient-specific inlet-boundary conditions measured with the optical flow method from both DSA and 3D rotational angiography datasets. Three independent observers compared qualitatively both the inflow direction and the apparent recirculation in regular DSA, optical flow images, and computational fluid dynamics flow patterns for each patient; κ statistics were estimated.

RESULTS

We included 21 patients. In 14 of these 21, the flow patterns were conclusive and matching between the optical flow images and computational fluid dynamics within the same projection view (κ = .91). However, in only 8 of these 14 patients the optical flow images were conclusive and matching regular DSA images (observer κ = 0.87). In 7 of the 21 patients, the flow patterns in the optical flow images were inconclusive, possibly due to improper projection angles.

CONCLUSIONS

The DSA-based optical flow technique was considered qualitatively consistent with computational fluid dynamics outcomes in evaluating intra-aneurysmal inflow direction and apparent recirculation. Moreover, the optical flow technique may provide the premises for new solutions for improving the visibility of flow patterns when contrast motion in DSA is not apparent. This technique is a diagnostic method to evaluate intra-aneurysmal flow patterns and could be used in the future for validation and patient evaluation.