The mixability of angiographic contrast with arterial blood

In silico intravascular optical coherence tomography (IVOCT) for quality assured imaging with reduced intervention

In the clinical application of intravascular optical coherence tomography (IVOCT), it is necessary to flush opaque blood during image acquisition. However, there are no specific standards for how to perform low-dose but effective flushing. In this study, computational fluid dynamics (CFD) and optical models were integrated to numerically simulate the complete process of IVOCT, which includes blood flushing with normal saline followed by image acquisition. Moreover, an intermittent injection scheme was proposed, and its advantages over the conventionally adopted scheme of continuous injection were verified. The results show that intermittent injection can significantly reduce the dosage of normal saline (reduced by 44.4%) with only a slight sacrifice of image quality (reduced by 8.7%, but still acceptable). The developed model and key findings in this work can help surgeons practice optimized IVOCT operations and potentially lead to improved designs of the IVOCT equipment.

Enhancing cerebral vasculature analysis with pathlength-corrected 2D angiographic parametric imaging: A feasibility study

BACKGROUND

2D angiographic parametric imaging (API) quantitatively extracts imaging biomarkers related to contrast flow and is conventionally applied to 2D digitally subtracted angiograms (DSA's). In the interventional suite, API is typically performed using 1-2 projection views and is limited by vessel overlap, foreshortening, and depth-integration of contrast motion.

PURPOSE

This work explores the use of a pathlength-correction metric to overcome the limitations of 2D-API: the primary objective was to study the effect of converting 3D contrast flow to projected contrast flow using a simulated angiographic framework created with computational fluid dynamics (CFD) simulations, thereby removing acquisition variability.

METHODS

The pathlength-correction framework was applied to in-silico angiograms, generating a reference (i.e., ground-truth) volumetric contrast distribution in four patient-specific intracranial aneurysm geometries. Biplane projections of contrast flow were created from the reference volumetric contrast distributions, assuming a cone-beam geometry. A Parker-weighted reconstruction was performed to obtain a binary representation of the vessel structure in 3D. Standard ray tracing techniques were then used to track the intersection of a ray from the focal spot with each voxel of the reconstructed vessel wall to a pixel in the detector plane. The lengths of each ray through the 3D vessel lumen were then projected along each ray-path to create a pathlength-correction map, where the pixel intensity in the detector plane corresponds to the vessel width along each source-detector ray. By dividing the projection sequences with this correction map, 2D pathlength-corrected in-silico angiograms were obtained. We then performed voxel-wise (3D) API on the ground-truth contrast distribution and compared it to pixel-wise (2D) API, both with and without pathlength correction for each biplane view. The percentage difference (PD) between the resultant API biomarkers in each dataset were calculated within the aneurysm region of interest (ROI).

RESULTS

Intensity-based API parameters, such as the area under the curve (AUC) and peak height (PH), exhibited notable changes in magnitude and spatial distribution following pathlength correction: these now accurately represent conservation of mass of injected contrast media within each arterial geometry and accurately reflect regions of stagnation and recirculation in each aneurysm ROI. Improved agreement was observed between these biomarkers in the pathlength-corrected biplane maps: the maximum PD within the aneurysm ROI is 3.3% with pathlength correction and 47.7% without pathlength correction. As expected, improved agreement with ROI-averaged ground-truth 3D counterparts was observed for all aneurysm geometries, particularly large aneurysms: the maximum PD for both AUC and PH was 5.8%. Temporal parameters (mean transit time, MTT, time-to-peak, TTP, time-to-arrival, TTA) remained unaffected after pathlength correction.

CONCLUSIONS

This study indicates that the values of intensity-based API parameters obtained with conventional 2D-API, without pathlength correction, are highly dependent on the projection orientation, and uncorrected API should be avoided for hemodynamic analysis. The proposed metric can standardize 2D API-derived biomarkers independent of projection orientation, potentially improving the diagnostic value of all acquired 2D-DSA's. Integration of a pathlength correction map into the imaging process can allow for improved interpretation of biomarkers in 2D space, which may lead to improved diagnostic accuracy during procedures involving the cerebral vasculature.

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.

Anatomically- and physiologically-informed computational model of hepatic contrast perfusion for virtual imaging trials

PURPOSE

Virtual (in silico) imaging trials (VITs), involving computerized phantoms and models of the imaging process, provide a modern alternative to clinical imaging trials. VITs are faster, safer, and enable otherwise-impossible investigations. Current phantoms used in VITs are limited in their ability to model functional behavior such as contrast perfusion which is an important determinant of dose and image quality in CT imaging. In our prior work with the XCAT computational phantoms, we determined and modeled inter-organ (organ to organ) intravenous contrast concentration as a function of time from injection. However, intra-organ concentration, heterogeneous distribution within a given organ, was not pursued. We extend our methods in this work to model intra-organ concentration within the XCAT phantom with a specific focus on the liver.

METHODS

Intra-organ contrast perfusion depends on the organ's vessel network. We modeled the intricate vascular structures of the liver, informed by empirical and theoretical observations of anatomy and physiology. The developed vessel generation algorithm modeled a dual-input-single-output vascular network as a series of bifurcating surfaces to optimally deliver flow within the bounding surface of a given XCAT liver. Using this network, contrast perfusion was simulated within voxelized versions of the phantom by using knowledge of the blood velocities in each vascular structure, vessel diameters and length, and the time since the contrast entered the hepatic artery. The utility of the enhanced phantom was demonstrated through a simulation study with the phantom voxelized prior to CT simulation with the relevant liver vasculature prepared to represent blood and iodinated contrast media. The spatial extent of the blood-contrast mixture was compared to clinical data.

RESULTS

The vascular structures of the liver were generated with size and orientation which resulted in minimal energy expenditure required to maintain blood flow. Intravenous contrast was simulated as having known concentration and known total volume in the liver as calibrated from time-concentration curves (TCC). Measurements of simulated CT ROIs were found to agree with clinically-observed values of early arterial phase contrast enhancement of the parenchyma (∼5 HU). Similarly, early enhancement in the hepatic artery was found to agree with average clinical enhancement (180 HU).

CONCLUSIONS

The computational methods presented here furthered the development of the XCAT phantoms allowing for multi-timepoint contrast perfusion simulations, enabling more anthropomorphic virtual clinical trials intended for optimization of current clinical imaging technologies and applications. This article is protected by copyright. All rights reserved.

Derivation of vascular wall shear stress from 1000 fps high-speed angiography (HSA) velocity distributions

Pathological changes in blood flow lead to altered hemodynamic forces, which are responsible for a number of conditions related to the remodeling and regeneration of the vasculature. More specifically, wall shear stress (WSS) has been shown to be a significant hemodynamic parameter with respect to aneurysm growth and rupture, as well as plaque activation leading to increased risk of stroke. In-vivo measurement of shear stress is difficult due to the stringent requirements on spatial resolution near the wall boundaries, as well as the deviation from the commonly assumed parabolic flow behavior at the wall. In this work, we propose an experimental method of in-vitro WSS calculations from high-temporal resolution velocity distributions, which are derived from 1000 fps high-speed angiography (HSA). The high-spatial and temporal resolution of our HSA detector makes such high-resolution velocity gradient measurements feasible. Presented here is the methodology for calculation of WSS in the imaging plane, as well as initial results for a variety of vascular geometries at physiologically realistic flow rates. Further, the effect of spatial resolution on the gradient calculation is explored using CFD-derived velocity data. Such angiographic-based analysis with HSA has the potential to provide critical hemodynamic feedback in an interventional setting, with the overarching objective of supporting clinical decision-making and improving patient outcomes.

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.

Evaluation of methods to derive blood flow velocity from 1000 fps high-speed angiographic sequences (HSA) using optical flow (OF) and computational fluid dynamics (CFD)

Digital Subtraction Angiography (DSA) is considered the gold standard for imaging and guiding treatment of neurovascular lesions, such as cerebral aneurysms and carotid stenoses. Though DSA can show high-resolution morphology, it remains difficult to extract temporal physiological information, because higher frame-rates are necessary to accurately quantify neurovascular flow details. Recent advances in photon-counting detector technology have led us to develop High-Speed Angiography (HSA), where X-ray images are acquired at 1000 fps for more accurate visualization and quantification of blood flow. Blood flow was imaged using HSA under constant flow conditions within various 3D printed patient-specific phantoms. Blood velocity was quantified using an open source Optical Flow algorithm, OpenOpticalFlow, to perform velocity estimation based on the spatio-temporal intensity changes of iodinated contrast wavefronts. The results of these algorithms are then compared with Computational Fluid Dynamics (CFD) simulations, using the same inlet boundary conditions and model geometries. The performance of these algorithms at lower temporal resolution was then also assessed by simulating lower frame rates from the acquired 1000 fps data. It is important to ascertain the hemodynamic effect of abnormal neurovascular conditions, as well as their effect on treatment of such conditions during the actual clinical interventional procedure. While theoretical CFD results requiring considerable computer capability are delayed for hours or more, it is expected that clinical results from multiple HSA sequences will be available almost immediately while the patient is still under treatment, and even right after flow conditions are changed beneficially by the intervention.

An in vitro study of pressure increases during contrast injections in diagnostic cerebral angiography

BACKGROUND

During diagnostic cerebral angiography, the contrast bolus injected into a vessel can cause substantial changes in baseline pressures and flows. One potential, and serious complication is the re-rupture of aneurysms due to these injections. The goals of this in vitro study were to evaluate the effect of injection conditions on intraneurysmal pressure changes during angiography.

METHODS

A silicone replica of a complete circle of Willis model with ophthalmic, anterior communicating, and basilar tip aneurysms was connected to a physiologically accurate flow pump. Contrast injections were performed under different conditions (carotid or vertebral vessel imaging, catheter diameter, injection rate, injection time, and arterial blood flow rate) and the pressure in each aneurysm was recorded before and during each injection. The effect of injection conditions on percentage increase in aneurysm pressures was statistically assessed. Additionally, the effect of the distance between the aneurysm and the catheter-tip on aneurysmal pressures was assessed.

RESULTS

Mean intraneurysmal pressures during injection (84.5 ± 10.8 mmHg) were significantly higher than pre-injection pressures (80.4 ± 10.6 mmHg, p < 0.0001). Only 3 of the 5 conditions - carotid injections, higher injection rates, and smaller catheter diameters - significantly increased intraneurysmal pressures. The catheter-tip distance showed no correlation to pressure increases.

CONCLUSIONS

Increasing contrast injection rates and decreasing catheter diameters are correlated to intraneurysmal pressure increases during angiography irrespective of the distance to the catheter tip. Future in vivo studies are required to confirm these findings and determine whether the amplitude of pressure increases with commonly used injection rates can be clinically detrimental.

4D-DSA: Development and Current Neurovascular Applications

Originally described by Davis et al in 2013, 4D-Digital Subtraction Angiography (4D-DSA) has developed into a commercially available application of DSA in the angiography suite. 4D-DSA provides the user with 3D time-resolved images, allowing observation of a contrast bolus at any desired viewing angle through the vasculature and at any time point during the acquisition (any view at any time). 4D-DSA mitigates some limitations that are intrinsic to both 2D- and 3D-DSA images. The clinical applications for 4D-DSA include evaluations of AVMs and AVFs, intracranial aneurysms, and atherosclerotic occlusive disease. Recent advances in blood flow quantification using 4D-DSA indicate that these data provide both the velocity and geometric information necessary for the quantification of blood flow. In this review, we will discuss the development, acquisition, reconstruction, and current neurovascular applications of 4D-DSA volumes.

Stagnant venous outflow in ruptured arteriovenous malformations revealed by delayed quantitative digital subtraction angiography

PURPOSE

To investigate the reproducibility of quantitative digital subtraction angiography (QDSA) measurements and their associations with brain arteriovenous malformation (BAVM) hemorrhage.

METHODS

From 2011-2019, 37 patients with BAVMs who had undergone both diagnostic and stereotactic DSA were divided into hemorrhagic and nonhemorrhagic groups. QDSA analysis was performed on the 2 DSA exams. The inter-exam reliabilities of QDSA measurements across the diagnostic and stereotactic DSA were tested using intraclass correlation coefficients (ICCs). Demographics, BAVM characteristics, and QDSA results for the hemorrhagic and nonhemorrhagic groups were compared.

RESULTS

Fifteen of 37 (40.5 %) patients presented with hemorrhage were associated with smaller BAVM volume and the presence of intranidal aneurysm and exclusive deep venous drainage. The median interval between the diagnostic and stereotactic DSA was 49 days and did not differ between the groups. In both groups, the inter-exam QDSA measurements were more reliable for drainage veins and transnidal time (ICCs ranged from 0.38-0.93) than for feeding arteries (ICCs ranged from 0.01-0.74). Among the venous parameters, the hemorrhagic group had lower peak density, area under the curve, inflow gradient, and outflow gradient on both DSA exams and larger full width at half maximum and stasis index on the stereotactic DSA exam than the nonhemorrhagic group.

CONCLUSIONS

In BAVMs, the QDSA measurements for veins are more reliable than those for arteries. QDSA analysis reflecting stagnant venous drainage is associated with BAVM hemorrhage, but may be confounded by the acute hemodynamic change after hemorrhage.

Visualization and Measurements of Blood Cells Flowing in Microfluidic Systems and Blood Rheology: A Personalized Medicine Perspective

Hemorheological alterations in the majority of metabolic diseases are always connected with blood rheology disturbances, such as the increase of blood and plasma viscosity, cell aggregation enhancement, and reduction of the red blood cells (RBCs) deformability. Thus, the visualizations and measurements of blood cells deformability flowing in microfluidic devices (point-of-care devices) can provide vital information to diagnose early symptoms of blood diseases and consequently to be used as a fast clinical tool for early detection of biomarkers. For instance, RBCs rigidity has been correlated with myocardial infarction, diabetes mellitus, hypertension, among other blood diseases. In order to better understand the blood cells behavior in microfluidic devices, rheological properties analysis is gaining interest by the biomedical committee, since it is strongly dependent on the interactions and mechanical cells proprieties. In addition, the development of blood analogue fluids capable of reproducing the rheological properties of blood and mimic the RBCs behavior at in vitro conditions is crucial for the design, performance and optimization of the microfluidic devices frequently used for personalized medicine. By combining the unique features of the hemorheology and microfluidic technology for single-cell analysis, valuable advances in personalized medicine for new treatments and diagnosis approach can be achieved.

Optimization of quantitative time-resolved 3D (4D) digital subtraction angiography in a porcine liver model

BACKGROUND

Time-resolved three-dimensional digital subtraction angiography (4D-DSA) can be used to quantify blood velocity. Contrast pulsatility, a major discriminant on 4D-DSA, is yet to be optimized. We investigated the effects of different imaging and injection parameters on sideband ratio (SBR), a measure of contrast pulsatile strength, within the hepatic vasculature of an in vivo porcine model.

METHODS

Fifty-nine hepatic 4D-DSA procedures were performed in three female domestic swine (mean weight 54 kg). Contrast injections were performed in the common hepatic artery with different combinations of imaging duration (6 s or 12 s), injection rates (from 1.0 to 2.5 mL/s), contrast concentration (50% or 100%), and catheter size (4 Fr or 5 Fr). Reflux was recorded. SBR and vessel cross-sectional areas were calculated in 289 arterial segments. Multiple linear mixed-effects models were estimated to determine the effects of parameters on SBR and cross-sectional vessel area.

RESULTS

Twelve-second acquisitions yielded a SBR higher than 6 s (p < 0.001). No significant differences in SBR were seen between different catheter sizes (p = 0.063) or contrast concentration (p = 0.907). For higher injection rates (2.5 mL/s), SBR was lower (p = 0.007) and cross-sectional area was higher (p < 0.001). Reflux of contrast does not significantly affect SBR (p = 0.087).

CONCLUSIONS

The strength of contrast pulsatility used for flow quantitation with 4D-DSA can be increased by adjusting injection rates and using longer acquisition times. Reduction of contrast concentration to 50% is feasible and reflux of contrast does not significantly hinder contrast pulsatility.

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.

Preliminary in vitro angiographic comparison of the flow diversion behavior of Evolve and Pipeline devices

Background and purpose

Flow diverters are increasingly used to treat a broad category of cerebral aneurysms. We conducted an in vitro study to angiographically compare the flow diversion effect of Surpass Evolve from Stryker Neurovascular with the Pipeline Shield Embolization Device from Medtronic Neurovascular.

Methods

Three copies each of three carotid aneurysm geometries were manufactured from silicone. Evolve and Pipeline flow diverters were deployed in one copy of each geometry; the third copy was used as Control. High-speed angiography was acquired under pulsatile flow in each replica, contrast concentration-time curves within the aneurysms were recorded, and the curves were quantified with six parameters. The parameters were statistically evaluated to compare the flow diversion effect of both devices.

Results

The Evolve showed greater flow diversion trends in almost all intra-geometry comparisons than the Pipeline. When aggregated over the three geometries, the Evolve was statistically significantly better than the Pipeline in four of the six parameters, and about the same or better (not statistically significant) than the Pipeline in the other two parameters.

Conclusions

The Evolve device demonstrated greater in vitro flow diversion effects than Pipeline. Comparative efficacy of the devices will need to be adjudicated based on clinical outcomes.

Angiographic assessment of the efficacy of flow diverter treatment for cerebral aneurysms

BACKGROUND

The recent growth of neuro-endovascular treatment has rekindled interest in the use of angiographic techniques for flow assessment. Aneurysm treatment with flow diverters is particularly amenable to such analysis. We analyze contrast time-density curves - recorded within aneurysms before (pre) and immediately after (post) flow diverter implantation to estimate six-month treatment outcomes.

METHODS

Fifty-six patients with 65 aneurysms were treated with flow diverters at two institutions. A region of interest was drawn around the aneurysm perimeter in image sequences taken both pre and post angiography, and the temporal variation in grayscale intensity within the aneurysm (time-density curve) was recorded. Eleven parameters were quantified from each time-density curve. Aneurysm occlusion status was recorded six months post treatment. The change in parameters from pre to post treatment was statistically evaluated between aneurysm occluded and non-occluded groups.

RESULTS

Of the 11 parameters, eight were significantly different before and immediately after flow diversion. Considering the entire data set, none of the parameters was statistically different between the occluded and non-occluded groups. However, subgroup analyses showed that four variables were significantly different between the aneurysm occluded and non-occluded groups. The sensitivity of these variables to predict aneurysm occlusion at six months ranged from 60% to 89%, while the specificity ranged from 55% to 70%.

CONCLUSIONS

Device-induced intra-aneurysmal flow alterations quantified by simple aneurysmal time-density curves can potentially be used to predict long-term outcomes of flow diversion. Large multi-center studies will be required to confirm these findings. Patient-to-patient variability in coagulation may need to be incorporated for clinically relevant predictive values.

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.

Three-Dimensional Quantitative Color-Coding Analysis of Hepatic Arterial Flow Change during Chemoembolization of Hepatocellular Carcinoma

PURPOSE

To evaluate feasibility of using three-dimensional (3D) quantitative color-coding analysis (QCA) to quantify substasis endpoints after transcatheter arterial chemoembolization of hepatocellular carcinoma (HCC).

MATERIALS AND METHODS

This single-institution prospective study included 20 patients with HCC who had undergone segmental or subsegmental transcatheter arterial chemoembolization between December 2015 and March 2017. The chemoembolization endpoint was a sluggish anterograde tumor-feeding arterial flow without residual tumor stains. Contrast medium bolus arrival time (BAT) was used as an indicator of arterial flow. BAT of the proper hepatic artery was obtained as a reference point. BATs of the proximal right lobar artery, proximal left lobar artery, and segmental artery that received embolization were analyzed before and after chemoembolization. Wilcoxon signed rank test was used to evaluate the difference between BATs before and after chemoembolization.

RESULTS

BATs before and after chemoembolization of the segmental artery that received embolization were 0.47 seconds (interquartile range [IQR], 0.31-0.70 s) and 1.04 seconds (IQR, 0.78-2.01 s; P < .001), respectively. BATs before and after chemoembolization of the proximal left lobar hepatic artery (0.35 s [IQR, 0.11-0.55] and 0.13 s [IQR, 0.05-0.32], P = .025) and right lobar hepatic artery (0.23 s [IQR, 0.13-0.65] and 0.22 s [IQR, 0.08-0.39], P = .027) exhibited no significant change.

CONCLUSIONS

3D QCA is a feasible method for quantifying sluggish segmental arterial flow after transcatheter arterial chemoembolization in patients with HCC.

Reliability and Accuracy of Peri-Interventional Stenosis Grading in Peripheral Artery Disease Using Color-Coded Quantitative Fluoroscopy: A Phantom Study Comparing a Clinical and Scientific Postprocessing Software

Purpose

To assess quantitative stenosis grading by color-coded fluoroscopy using an in vitro pulsatile flow phantom.

Methods

Three different stenotic tubes (80%, 60%, and 40% diameter restriction) and a nonstenotic reference tube were compared regarding their different flow behavior by using contrast-enhanced fluoroscopy with a flat-detector system for visualisation purposes. Time-density curves (TDC), area under the curve (AUC), time-to-peak (TTP), and different ROI sizes were analyzed in three independent measurements using two different postprocessing software solutions. In addition, exemplary TDCs of a patient with a high-grade stenosis before and after stent angioplasty were acquired.

Results

Color-coded fluoroscopy enabled depiction of differences in AUC and TDC between high-grade (80%), middle (60%), low-grade (40%), and nonstenotic tubes. The best correlation between high-, middle-, and low-grade stenosis was appreciated in ROIs behind the stenosis. This effect was enhanced by using longer integration times (5s, 7s) and a maximum frame rate of image acquisition for analysis (correlation coefficient rho=0.9284 at 5s). TTP showed no significant differences between high- and low-grade stenosis.

Conclusions

Various clinical studies in the literature already demonstrated reproducible and reliable stenosis grading by analyzing TDCs acquired with color-coded fluoroscopy. In contrast to TTP, AUC values derived in ROIs behind the stenosis proved to be reliable parameters for stenosis grading. However, our results also demonstrate that several factors are able to significantly impact the evaluation of AUC values. More precisely, accuracy of acquired AUC values can be improved by choosing longer integration times, a large ROI size adapted to the vessel diameter, and a higher frame rate of image acquisition.

Pressure and Flow Rate Changes During Contrast Injections in Cerebral Angiography: Correlation to Reflux Length

Cerebral angiography involves the antegrade injection of contrast media through a catheter into the vasculature to visualize the region of interest under X-ray imaging. Depending on the injection and blood flow parameters, the bolus of contrast can propagate in the upstream direction and proximal to the catheter tip, at which point contrast is said to have refluxed. In this in vitro study, we investigate the relationship of fundamental hemodynamic variables to this phenomenon. Contrast injections were carried out under steady and pulsatile flow using various vessel diameters, catheter sizes, working fluid flow rates, and injection rates. The distance from the catheter tip to the proximal edge of the contrast bolus, called reflux length, was measured on the angiograms; the relation of this reflux length to different hemodynamic parameters was evaluated. Results show that contrast reflux occurs when the pressure distal to the catheter tip increases to be greater than the pressure proximal to the catheter tip. The ratio of this pressure difference to the baseline flow rate, called reflux resistance here, was linearly correlated to the normalized reflux length (reflux length/vessel diameter). Further, the ratio of blood flow to contrast fluid momentums, called the Craya-Curtet number, was correlated to the normalized reflux length via a sigmoid function. A sigmoid function was also found to be representative of the relationship between the ratio of the Reynolds numbers of blood flow to contrast and the normalized reflux length. As described by previous reports, catheter based contrast injections cause substantial increases in local flow and pressure. Contrast reflux should generally be avoided during standard antegrade angiography. Our study shows two specific correlations between contrast reflux length and baseline and intra-injection parameters that have not been published previously. Further studies need to be conducted to fully characterize the phenomena and to extract reliable indicators of clinical utility. Parameters relevant to cerebral angiography are studied here, but the essential principles are applicable to all angiographic procedures involving antegrade catheter injections.

In vitro angiographic comparison of the flow-diversion performance of five neurovascular stents

Background and purpose Data differentiating flow diversion properties of commercially available low- and high-porosity stents are limited. This in vitro study applies angiographic analysis of intra-aneurysmal flow to compare the flow-diversion performance of five neurovascular devices in idealized sidewall and bifurcation aneurysm models. Methods Five commercial devices (Enterprise, Neuroform, LVIS, FRED, and Pipeline) were implanted in silicone sidewall and bifurcation aneurysm models under physiological average flow of blood analog fluid. High-speed angiographic images were acquired pre- and post-device implantation and contrast concentration-time curves within the aneurysm were recorded. The curves were quantified with five parameters to assess changes in contrast transport, and thus aneurysm hemodynamics, due to each device. Results Inter-device flow-diversion performance was more easily distinguished in the sidewall model than the bifurcation model. There were no obvious overall statistical trends in the bifurcation parameters but the Pipeline performed marginally better than the other devices. In the sidewall geometry, overall evidence suggests that the LVIS performed better than the Neuroform and Enterprise. The Pipeline and FRED devices were statistically superior to the three stents and Pipeline was superior to FRED in all sidewall parameters evaluated. Conclusions Based on this specific set of experiments, lower-porosity flow diverters perform significantly better in reducing intra-aneurysmal flow activity than higher-porosity stents in sidewall-type geometries. The LVIS device is potentially a better flow diverter than the Neuroform and Enterprise devices, while the Pipeline is potentially better than the FRED.

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 investigation of contrast flow jet timing in patient-specific intracranial aneurysms

BACKGROUND

The direction and magnitude of intra-aneurysmal flow jet are significant risk factors of subarachnoid hemorrhage, and the change of flow jet during an endovascular procedure has been used for prediction of aneurysm occlusion or whether an additional flow diverter (FD) is warranted. However, evaluation of flow jets is often unreliable due to a large variation of flow jet on the digital subtraction angiograms, and this flow pattern variation may result in incorrect clinical diagnosis Therefore, factors contributing to the variation in flow jet are examined at an in vitro setting, and the findings can help us to understand the nature of flow jet and devise a better plan to quantify the aneurysmal hemodynamics accurately.

METHODS

Intra-aneurysmal flows in three patient-specific aneurysms between 11 and 25 mm were investigated in vitro, and a FD was deployed in each aneurysm model. X-ray imaging of these models were performed at injection rates between 0.2 and 2 mL/s. Pulsatile blood pump and aneurysm model were imaged together to determine the timing of flow jet.

RESULTS

The contrast bolus arrives at the aneurysm early at high contrast injection rates. The flow patterns with slow injection rates exhibit strong inertia that is associated with the systole flow. Flow jets arrive at the aneurysms at the peak systole when the bolus is injected at 0.2 mL/s. The contrast-to-signal ratio is the highest at the injection rate of 0.5 mL/s. Effect of flow diversion can only be assessed at an injection rate greater than 0.5 mL/s.

CONCLUSIONS

Intra-aneurysmal flow jet is highly dependent on the injection rate of the contrast agent. For the internal carotid artery (ICA) aneurysms, the systolic flows can be visualized at slow injection rates (<0.5 mL/s), while the diastolic flow jets are visible at higher injection rates (>1 mL/s). Dependence of flow jet on the contrast injection rate has serious clinical implications and needs to be considered during diagnostic procedures; a protocol with a consistent injection rate is highly recommended.

A Fourier-based approach to the angiographic assessment of flow diverter efficacy in the treatment of cerebral aneurysms

Flow diversion is an emerging endovascular treatment option for cerebral aneurysms. Quantitative assessment of hemodynamic changes induced by flow diversion can aid clinical decision making in the treatment of cerebral aneurysms. In this article, besides summarizing past key research efforts, we propose a novel metric for the angiographic assessment of flow diverter deployments in the treatment of cerebral aneurysms. By analyzing the frequency spectra of signals derived from digital subtraction angiography (DSA) series, the metric aims to quantify the prevalence of frequency components that correspond to the patient-specific heart rate. Indicating the decoupling of aneurysms from healthy blood circulation, our proposed metric could advance clinical guidelines for treatment success prediction. The very promising results of a retrospective feasibility study on 26 DSA series warrant future efforts to study the validity of the proposed metric within a clinical setting.

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

BACKGROUND AND PURPOSE

Digital subtraction angiography is the reference standard technique to evaluate intracranial vascular anatomy and used on the endovascular treatment of vascular diseases. A dedicated optical flow-based algorithm was applied to DSA to measure arterial flow. The first quantification results of internal carotid artery flow validated with Doppler sonography are reported.

MATERIALS AND METHODS

We included 22 consecutive patients who underwent endovascular procedures. To assess the sensitivity of the algorithm to contrast agent-blood mixing dynamics, we acquired high-frame DSA series (60 images/s) with different injection rates: 1.5 mL/s (n = 19), 2.0 mL/s (n = 18), and 3.0 mL/s (n = 13). 3D rotational angiography was used to extract the centerline of the vessel and the arterial section necessary for volume flow calculation. Optical flow was used to measure flow velocities in straight parts of the ICAs; these data were further compared with Doppler sonography data. DSA mean flow rates were linearly regressed on Doppler sonography measurements, and regression slope coefficient bias from value 1 was analyzed within the 95% confidence interval.

RESULTS

DSA mean flow rates measured with the optical flow approach significantly matched Doppler sonography measurements (slope regression coefficient, b = 0.83 ± 0.19, P = .05) for injection rate = 2.0 mL/s and circulating volumetric blood flow <6 mL/s. For injection rate = 1.5 mL/s, volumetric blood flow <3 mL/s correlated well with Doppler sonography (b = 0.67 ± 0.33, P = .05). Injection rate = 3.0 mL/s failed to provide DSA-optical flow measurements correlating with Doppler sonography because of the lack of measurable pulsatility.

CONCLUSIONS

A new model-free optical flow technique was tested reliably on the ICA. DSA-based blood flow velocity measurements were essentially validated with Doppler sonography whenever the conditions of measurable pulsatility were achieved (injection rates = 1.5 and 2.0 mL/s).

Effect of saline injection mixing on accuracy of conductance lumen sizing of peripheral vessels

Transient displacement of blood in vessel lumen with saline injection is necessary in the conductance method for measurement of arterial cross-sectional area (CSA). The displacement of blood is dictated by the interactions between arterial flow hemodynamics and saline injection dynamics. The objective of the present study is to understand how the accuracy of conductance measurements is affected by the saline injection. Computational simulations were performed to assess the error in predictions of arterial CSA using conductance measurements over a range of peripheral artery diameters (i.e., 4, 7, and 10 mm) with an introducing catheter (6 Fr.) for various blood flow and saline injection rates. The simulation results were validated using the conductance measurements of the phantoms with known diameters (i.e., 7 and 10 mm). The results demonstrated that a minimum ratio of saline injection rate to blood flow rate of 3 is needed to fully displace the blood and result in accurate measurement of CSA for the peripheral artery sizes considered. Furthermore, the error was shown to be minimized as the detection electrodes are positioned between the distal to the mixing zone induced by saline injection and far downstream (4–8 cm from the injection catheter tip). The present study shows that even for the large peripheral arteries (7–10 mm) where mixing can occur, an appropriate injection rate and detection position can produce accurate measurement of lumen size.

Comparison of aortic arch and intravenous contrast injection techniques for C-arm cone beam CT: implications for cerebral perfusion imaging in the angiography suite

RATIONALE AND OBJECTIVES

The ability to perform cerebral perfusion imaging (CPI) in the angiography suite has provided a new tool for diagnosis and treatment of neurovascular patients but requires comparable contrast perfusion to each cerebral hemisphere. In the angiography suite, contrast injection may be performed via an intra-arterial or intravenous (IV) route. The purpose of this study was to investigate whether a difference exists between contrast injection in the aortic arch (AA) and a peripheral vein (IV), particularly in the setting of stroke.

MATERIALS AND METHODS

Using three canines, both AA and IV injection protocols compatible with CPI were performed prospectively at three time points after creation of a stroke. The common carotid arteries in the resulting image data sets were segmented and the means and distributions of corresponding pixel intensities analyzed with Student's t-test. Using similar techniques, the internal carotid arteries of three patients (one female, two males, ages 69, 29, and 20) undergoing AA contrast injection with cone beam computed tomography (CBCT) cerebral imaging were analyzed and compared retrospectively with those of three random patients (one female, two males, ages 19, 57, and 35) undergoing standard head CT scans using IV contrast administration. All acquisitions followed institutionally approved protocols and informed consent.

RESULTS

No statistical significance (P < .05) was found when mean values for the right and left carotid artery pixel intensities were compared in the canine model or the clinical studies in which patients underwent imaging after AA or IV contrast administration.

CONCLUSIONS

No statistically significant difference exists between right and left carotid artery filling density using either AA or IV contrast injection methods, making both suitable for CPI in the angiography suite.

A workflow for patient-individualized virtual angiogram generation based on CFD simulation

Increasing interest is drawn on hemodynamic parameters for classifying the risk of rupture as well as treatment planning of cerebral aneurysms. A proposed method to obtain quantities such as wall shear stress, pressure, and blood flow velocity is to numerically simulate the blood flow using computational fluid dynamics (CFD) methods. For the validation of those calculated quantities, virtually generated angiograms, based on the CFD results, are increasingly used for a subsequent comparison with real, acquired angiograms. For the generation of virtual angiograms, several patient-specific parameters have to be incorporated to obtain virtual angiograms which match the acquired angiograms as best as possible. For this purpose, a workflow is presented and demonstrated involving multiple phantom and patient cases.

Dispersive Transport of Angiographic Contrast During Antegrade Arterial Injection

PURPOSE: Angiography is commonly used during endovascular procedures to navigate catheters into a target artery and for evaluation of the arterial luminal geometry. X-ray attenuating contrast material is injected into the arteries and transported into pathologies such as aneurysms or arteriovenous malformations. Images of the transported contrast are used to guide therapeutic decisions. Experience and intuition of the interventionalist are often serving as guide for the injection force, and hence, the speed and volume of the bolus. Forceful injections of small boluses can evoke local turbulence and dispersive mixing in the zone immediately distal to the catheter tip. Turbulence by its nature acts as a strong agitating mechanism such that the bolus of contrast quickly mixes with the flowing blood to occupy the entire lumen so the artery can be visualized. The aims of the present study are (a) to determine the distance from catheter tip beyond which contrast can consider to be fully mixed with the blood during antegrade injection and (b) to determine the thickness of the boundary layer in which contrast concentration is poor, which can contribute to underestimation of vascular diameter using this method. METHODS: We performed in silico experiments to describe blood and angiographic contrast transport in a straight artery model. The conditions investigated are derived from clinical contrast injection rates typically found in cerebral angiography. RESULTS: A recirculation flow exists in the mixing zone distal to the catheter tip issuing the contrast and convective mixing rather than diffusion is dominating the rapid mixing process. In the vicinity of the arterial wall in the mass transfer boundary layer, however, transport is dominated by molecular diffusion. For lower molecule diffusion coefficient, the mass transfer boundary layer contains a lower concentration of contrast than for a higher molecular diffusion coefficient. CONCLUSIONS: These findings imply that contrast visibility near the arterial wall is poor such that arterial dimensions derived from angiograms may be underestimated and consequently sizing of potential implants inaccurate. Outside the mass transfer boundary layer contrast can be considered as fully mixed with the carrying flow in about 10 arterial diameters distal to the injection port.

Comprehensive validation of computational fluid dynamics simulations of in-vivo blood flow in patient-specific cerebral aneurysms

PURPOSE

Recently, image-based computational fluid dynamic (CFD) simulations have been proposed to investigate the local hemodynamics inside human cerebral aneurysms. It was suggested that the knowledge of the computed three-dimensional flow fields can be used to assist clinical risk assessment and treatment decision making. Therefore, it was desired to know the reliability of CFD for cerebral blood flow simulation, and be able to provide clinical feedback. However, the validations were not yet comprehensive as they lack either patient-specific boundary conditions (BCs) required for CFD simulations or quantitative comparison methods.

METHODS

In this study, based on a recently proposed in-vitro quantitative CFD evaluation approach via virtual angiography, the CFD evaluation was extended from phantom to patient studies. In contrast to previous work, patient-specific blood flow rates obtained by transcranial color coded Doppler ultrasound measurements were used to impose CFD BCs. Virtual angiograms (VAs) were constructed which resemble clinically acquired angiograms (AAs). Quantitative measures were defined to thoroughly evaluate the correspondence of the detailed flow features between the AAs and the VAs, and thus, the reliability of CFD simulations.

RESULTS

The proposed simulation pipeline provided a comprehensive validation method of CFD simulation for reproducing cerebral blood flow, with a focus on the aneurysm region. Six patient cases were tested and close similarities were found in terms of spatial and temporal variations of contrast agent (CA) distribution between AAs and VAs. For patient #1 to #5, discrepancies of less than 11% were found for the relative root mean square errors in time intensity curve comparisons from characteristic vasculature positions. For patient #6, where the CA concentration curve at vessel inlet cannot be directly extracted from the AAs and given as a BC, deviations about 20% were found.

CONCLUSIONS

As a conclusion, the reliability of the CFD simulations was well confirmed. Besides, it was shown that the accuracy of CFD simulations was closely related to the input BCs.

Cerebral CT perfusion using an interventional C-arm imaging system: cerebral blood flow measurements

BACKGROUND AND PURPOSE

CTP imaging in the interventional suite could reduce delays to the start of image-guided interventions and help determine the treatment progress and end point. However, C-arms rotate slower than clinical CT scanners, making CTP challenging. We developed a cerebral CTP protocol for C-arm CBCT and evaluated it in an animal study.

MATERIALS AND METHODS

Five anesthetized swine were imaged by using C-arm CBCT and conventional CT. The C-arm rotates in 4.3 seconds plus a 1.25-second turnaround, compared with 0.5 seconds for clinical CT. Each C-arm scan had 6 continuous bidirectional sweeps. Multiple scans each with a different delay to the start of an aortic arch iodinated contrast injection and a novel image reconstruction algorithm were used to increase temporal resolution. Three different scan sets (consisting of 6, 3, or 2 scans) and 3 injection protocols (3-mL/s 100%, 3-mL/s 67%, and 6-mL/s 50% contrast concentration) were studied. CBF maps for each scan set and injection were generated. The concordance and Pearson correlation coefficients (ρ and r) were calculated to determine the injection providing the best match between the following: the left and right hemispheres, and CT and C-arm CBCT.

RESULTS

The highest ρ and r values (both 0.92) for the left and right hemispheres were obtained by using the 6-mL 50% iodinated contrast concentration injection. The same injection gave the best match for CT and C-arm CBCT for the 6-scan set (ρ = 0.77, r = 0.89). Some of the 3-scan and 2-scan protocols provided matches similar to those in CT.

CONCLUSIONS

This study demonstrated that C-arm CBCT can produce CBF maps that correlate well with those from CTP.