X-ray videodensitometric methods for blood flow and velocity measurement: a critical review of literature

Experimental study using phantom models of cerebral aneurysms and 4D-DSA to measure blood flow on 3D-color-coded images

BACKGROUND

The current 3D-iFlow application can only measure the arrival time of contrast media through intensity values. If the flow rate could be estimated by 3D-iFlow, patient-specific hemodynamics could be determined within the scope of normal diagnostic management, eliminating the need for additional resources for blood flow rate estimation.

OBJECTIVE

The aim of this study is to develop and validate a method for measuring the flow rate by data obtained from 3D-iFlow images - a prototype application in Four-dimensional digital subtraction angiography (4D-DSA).

METHODS

Using phantom model and experimental circuit with circulating glycerin solution, an equation for the relationship between contrast media intensity and flow rate was developed. Applying the equation to the aneurysm phantom models, the derived flow rate was evaluated.

RESULTS

The average errors between the derived flow rate and setting flow rate became larger when the glycerin flow and the X-rays from the X-ray tube of the angiography system were parallel to each other or when the measurement point included overlaps with other contrast enhanced areas.

CONCLUSION

Although the error increases dependent on the imaging direction and overlap of contrast enhanced area, the developed equation can estimate the flow rate using the image intensity value measured on 3D-iFlow based on 4D-DSA.

Angiography-based hemodynamic features predict recurrent ischemic events after angioplasty and stenting of intracranial vertebrobasilar atherosclerotic stenosis

OBJECTIVES

To assess the predictive value of hemodynamic features for stroke relapse in patients with intracranial vertebrobasilar atherosclerotic stenosis treated with percutaneous transluminal angioplasty and stenting (PTAS) using quantitative digital subtraction angiography (q-DSA).

METHODS

In this retrospective longitudinal study, patients with intracranial vertebrobasilar atherosclerotic stenosis and who underwent PTAS treatment between January 2012 and May 2020 were enrolled. The q-DSA assessment was performed before and after PTAS. ROIs 1-4 were placed along the vertebral artery, proximal and distal basilar artery, and posterior cerebral artery; ROIs 5-8 were in 5 mm and 10 mm proximal and distal to the lesion, respectively. Relative time to peak (rTTP) was defined as the difference in TTP between ROIs. Cox regression analyses were performed to determine risk factors for recurrent stroke.

RESULTS

A total of 137 patients (mean age, 62 years ± 10 [standard deviation], 83.2% males) were included, and 26 (19.0%) patients had stroke relapse during follow-up (median time of 42.6 months [interquartile range, 19.7-60.7]). Preprocedural rTTP4-1 (adjusted hazard ratio (HR) = 2.270; 95% CI 1.371-3.758; p = 0.001) and preprocedural rTTP8-5 (adjusted HR = 0.240; 95% CI 0.088-0.658; p = 0.006) were independently associated with the recurrent stroke. These hemodynamic parameters provided an incremental prognostic value for stroke relapse (AUC, 0.817 [0.704-0.931]; the net reclassification index, 0.431 [0.057-0.625]; and the integrated discrimination index, 0.140 [0.035-0.292]).

CONCLUSIONS

In patients with intracranial vertebrobasilar atherosclerosis treated with PTAS, preprocedural prolonged TTP of the target vessel and shortened trans-stenotic TTP difference were associated with stroke relapse. Q-DSA-defined hemodynamic parameters provided incremental predictive value over conventional parameters for stroke recurrence.

CLINICAL RELEVANCE STATEMENT

Quantitative DSA analysis enables intuitive observation and semi-quantitative evaluation of peri-therapeutic cerebral blood flow. More importantly, quantitative DSA-defined hemodynamic parameters have the potential for risk stratification of patients with intracranial atherosclerotic stenosis.

KEY POINTS

Semi-quantitative angiography-based parameters can reflect pre- and postprocedural subtle changes in blood flow in patients with intracranial atherosclerotic stenosis. Although angioplasty procedures can significantly improve blood flow status, patients with more restricted baseline blood flow still show a higher risk of stroke recurrence. Angiography-based hemodynamic features possess prognostic value and can serve as clinical markers to assess stroke risk of patients with intracranial atherosclerotic stenosis.

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.

Interleaved x-ray imaging: A method for simultaneous acquisition of quantitative and diagnostic digital subtraction angiography

BACKGROUND

Flow altering angiographic procedures suffer from ill-defined, qualitative endpoints. Quantitative digital subtraction angiography (qDSA) is an emerging technology that aims to address this issue by providing intra-procedural blood velocity measurements from time-resolved, 2D angiograms. To date, qDSA has used 30 frame/s DSA imaging, which is associated with high radiation dose rate compared to clinical diagnostic DSA (up to 4 frame/s).

PURPOSE

The purpose of this study is to demonstrate an interleaved x-ray imaging method which decreases the radiation dose rate associated with high frame rate qDSA while simultaneously providing low frame rate diagnostic DSA images, enabling the acquisition of both datasets in a single image sequence with a single injection of contrast agent.

METHODS

Interleaved x-ray imaging combines low radiation dose image frames acquired at a high rate with high radiation dose image frames acquired at a low rate. The feasibility of this approach was evaluated on an x-ray system equipped with research prototype software for x-ray tube control. qDSA blood velocity quantification was evaluated in a flow phantom study for two lower dose interleaving protocols (LD1:3.7 ± 0.02 mGy / s $3.7 \pm 0.02\ {\mathrm{mGy}}/{\mathrm{s}}$ and LD2:1.7 ± 0.04 mGy / s $1.7 \pm 0.04{\mathrm{\ mGy}}/{\mathrm{s}}$ ) and one conventional (full dose) protocol (11.4 ± 0.04 mGy / s ) $11.4 \pm 0.04{\mathrm{\ mGy}}/{\mathrm{s}})$ . Dose was measured at the interventional reference point. Fluid velocities ranging from 24 to 45 cm/s were investigated. Gold standard velocities were measured using an ultrasound flow probe. Linear regression and Bland-Altman analysis were used to compare ultrasound and qDSA.

RESULTS

The LD1 and LD2 interleaved protocols resulted in dose rate reductions of -67.7% and -85.5%, compared to the full dose qDSA scan. For the full dose protocol, the Bland-Altman limits of agreement (LOA) between qDSA and ultrasound velocities were [0.7, 6.7] cm/s with a mean difference of 3.7 cm/s. The LD1 interleaved protocol results were similar (LOA: [0.3, 6.9] cm/s, bias: 3.6 cm/s). The LD2 interleaved protocol resulted in slightly larger LOA: [-2.5, 5.5] cm/s with a decrease in the bias: 1.5 cm/s. Linear regression analysis showed a strong correlation between ultrasound and qDSA derived velocities using the LD1 protocol, with aR 2 ${R}^2$ of0.96 $0.96$ , a slope of1.05 $1.05$ and an offset of1.9 $1.9$  cm/s. Similar values were also found for the LD2 protocol, with aR 2 ${R}^2$ of0.93 $0.93$ , a slope of0.98 $0.98$ and an offset of2.0 $2.0$  cm/s.

CONCLUSIONS

The interleaved method enables simultaneous acquisition of low-dose high-rate images for intra-procedural blood velocity quantification (qDSA) and high-dose low-rate images for vessel morphology evaluation (diagnostic DSA).

Spatiotemporal frequency domain analysis for blood velocity measurement during embolization procedures

BACKGROUND

Currently, determining procedural endpoints and treatment efficacy of vascular interventions is largely qualitative and relies on subjective visual assessment of digital subtraction angiography (DSA) images leading to large interobserver variabilities and poor reproducibility. Quantitative metrics such as the residual blood velocity in embolized vessel branches could help establish objective and reproducible endpoints. Recently, velocity quantification techniques based on a contrast enhanced X-ray sequence such as qDSA and 4D DSA have been proposed. These techniques must be robust, and, to avoid radiation dose concerns, they should be compatible with low dose per frame image acquisition.

PURPOSE

To develop and evaluate a technique for robust blood velocity quantification from low dose contrast enhanced X-ray image sequences that leverages the oscillating signal created by pulsatile blood flow.

METHODS

The proposed spatiotemporal frequency domain (STF) approach quantifies velocities from time attenuation maps (TAMs) representing the oscillating signal over time for all points along a vessel centerline. Due to the time it takes a contrast bolus to travel along the vessel centerline, the resulting TAM resembles a sheared sine wave. The shear angle is related to the velocity and can be determined in the spatiotemporal frequency domain after applying the 2D Fourier transform to the TAM. The approach was evaluated in a straight tube phantom using three different radiation dose levels and compared to ultrasound transit-time-based measurements. The STF velocity results were also compared to previously published approaches for the measurement of blood velocity from contrast enhanced X-ray sequences including shifted least squared (SLS) and phase shift (PHS). Additionally, an in vivo porcine study (n = 8) was performed where increasing amounts of embolic particles were injected into a hepatic or splenic artery with intermittent velocity measurements after each injection to monitor the resulting reduction in velocity.

RESULTS

At the lowest evaluated dose level (average air kerma rate 1.3 mGy/s at the interventional reference point), the Pearson correlation between ultrasound and STF velocity measurements was99 % $99\%$ . This was significantly higher (p < 0.0001 $p < 0.0001$ ) than corresponding correlation results between ultrasound and the previously published SLS and PHS approaches (91 $\hskip.001pt 91$ and93 % $93\%$ , respectively). In the in vivo study, a reduction in velocity was observed in85.7 % $85.7\%$ of cases after injection of 1 mL,96.4 % $96.4\%$ after 3 mL, and100.0 % $100.0\%$ after 4 mL of embolic particles.

CONCLUSIONS

The results show good agreement of the spatiotemporal frequency domain approach with ultrasound even in low dose per frame image sequences. Additionally, the in vivo study demonstrates the ability to monitor the physiological changes due to embolization. This could provide quantitative metrics during vascular procedures to establish objective and reproducible endpoints.

Motion-compensation approach for quantitative digital subtraction angiography and its effect on in-vivo blood velocity measurement

Purpose

Quantitative monitoring of flow-altering interventions has been proposed using algorithms that quantify blood velocity from time-resolved two-dimensional angiograms. These algorithms track the movement of contrast oscillations along a vessel centerline. Vessel motion may occur relative to a statically defined vessel centerline, corrupting the blood velocity measurement. We provide a method for motion-compensated blood velocity quantification.

Approach

The motion-compensation approach utilizes a vessel segmentation algorithm to perform frame-by-frame vessel registration and creates a dynamic vessel centerline that moves with the vasculature. Performance was evaluated in-vivo through comparison with manually annotated centerlines. The method was also compared to a previous uncompensated method using best- and worst-case static centerlines chosen to minimize and maximize centerline placement accuracy. Blood velocities determined through quantitative DSA (qDSA) analysis for each centerline type were compared through linear regression analysis.

Results

Centerline distance errors were 0.3 ± 0.1    mm relative to gold standard manual annotations. For the uncompensated approach, the best- and worst-case static centerlines had distance errors of 1.1 ± 0.6 and 2.9 ± 1.2    mm , respectively. Linear regression analysis found a high R -squared between qDSA-derived blood velocities using gold standard centerlines and motion-compensated centerlines (R 2 = 0.97 ) with a slope of 1.15 and a small offset of - 0.6    cm / s . The use of static centerlines resulted in low coefficients of determination for the best case (R 2 = 0.35 ) and worst-case (R 2 = 0.20 ) scenarios, with slopes close to zero.

Conclusions

In-vivo validation of motion-compensated qDSA analysis demonstrated improved velocity quantification accuracy in vessels with motion, addressing an important clinical limitation of the current qDSA algorithm.

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.

Estimating Pulsatile Flow Velocity using Four-Dimensional Digital Subtraction Angiography. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2023;2023:1-4. [PMID: 38082691 DOI: 10.1109/embc40787.2023.10340843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Digital subtraction angiography (DSA) is a X-ray based imaging modality for interventional procedures, and remains criterion standard for diagnosing vascular diseases. The imaging protocol of DSA involves administration of a foreign contrast medium into the blood vessel that opacifies the vasculature during the imaging. Using two-dimensional (2D) DSA with high temporal resolution, it was recently demonstrated that the pulsatile velocity can be estimated by evaluating the temporal and spatial variations of the contrast medium distributions in the blood vessel. In this paper, we evaluate the feasibility to estimate the pulsatile flow velocity using the four-dimensional (4D) DSA. To overcome the noise and artefacts of 4D-DSA data, a empirical mode decomposition plus autocorrelation based method is proposed to estimate the pulsatile velocities, and the pulsatile velocities estimated using 2D-DSA are used as reference for comparison.Clinical Relevance-4D-DSA encompasses both structural and temporal information; it theoretically reduces the need of multiple scans, hence reducing the radiation doses. The estimated pulsatile flow velocities open up a new parameter for hemodynamic studies and potential for real-time diagnostic and therapeutic monitoring during interventional procedures.

Effects of helical centerline stent vs. straight stent placement on blood flow velocity

As an approach to maintain patency in femoropopliteal stenting, a helical stent configuration was proposed, which showed improved patency in clinical trials. However, the effects of helical stent placement on the flow have not been quantitatively analyzed. The purpose of this study was to estimate flow velocities to quantify the influence of helical stent placement. Helical and straight stents were implanted in three healthy pigs, and the flow velocities were estimated using the time-intensity curve (TIC) in the angiography images. The angiographic images indicated thinning of the leading edge of the contrast medium through the helically deformed artery, which was not observed in the straight stent. The slower rise of the TIC peak in the helical stent indicated faster travel of this thinner edge. Arterial expansion due to stenting was observed in all cases, and the expansion rate varied according to location. All cases of helical stent implantation showed that velocity was maintained (55.0%-71.3% velocity retention), unlike for straight stent implantation (43.0%-68.0% velocity retention); however, no significant difference was observed.

Dispersion of Heterogeneous Medium in Pulsatile Blood Flow and Absolute Pulsatile Flow Velocity Quantification

Heterogeneous medium enhanced angiogr- ams are key diagnostic tools in clinical practice; the associated hemodynamic information is crucial for diagnosing cardiovascular diseases. However, the dynamics of such medium in physiological blood flow are poorly understood. Herein, we report a previously unnoticed dispersion pattern, which is a universal phenomenon, of a medium in pulsatile blood flow. We present a physical theory for studying the dispersion of a steadily injected heterogeneous medium into a thin tubular blood vessel in which the blood flow is pulsatile. In a thin tubular blood vessel, we demonstrate that variations of concentration associated with the heterogeneous medium obey a one-dimensional advection diffusion equation, and the diffusion has limited effect whenever a short vascular segment is considered. A distinct feature of the distribution of the medium in the axial distance-time plane is a "dilation-retraction" pattern. The time evolution signals at different axial positions exhibit distinct concentration waveforms. A numerical scheme is proposed for exploiting this information to estimate the pulsatile velocity. Artificial data are adopted to validate the scheme. Real X-ray angiography is also analyzed to support our theory and method. The theory is applicable whenever imaging protocols involve a heterogeneous medium in pulsatile flow.

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.

Quantitative transport mapping (QTM) for differentiating benign and malignant breast lesion: Comparison with traditional kinetics modeling and semi-quantitative enhancement curve characteristics

PURPOSE

To test the feasibility of using quantitative transport mapping (QTM) method, which is based on the inversion of transport equation using spatial deconvolution without any arterial input function, for automatically postprocessing dynamic contrast enhanced MRI (DCE-MRI) to differentiate malignant and benign breast tumors.

MATERIALS AND METHODS

Breast DCE-MRI data with biopsy confirmed malignant (n = 13) and benign tumors (n = 13) was used to assess QTM velocity (|u|) and diffusion coefficient (D), volume transfer constant (K^trans), volume fraction of extravascular extracellular space (V_e) from kinetics method, and traditional enhancement curve characteristics (ECC: amplitude A, wash-in rate α, wash-out rate β). A Mann-Whitney U test and receiver operating characteristic curve (ROC) analysis were performed to assess the diagnostic performance of these parameters for distinguishing between benign and malignant tumors.

RESULTS

Between malignant and benign tumors, there was a significant difference in |u| and K^trans, (p = 0.0066, 0.0274, respectively), but not in D, V_e, A, α and β (p = 0.1119, 0.2382, 0.4418,0.2592 and 0.9591, respectively). ROC area-under-the-curve was 0.82, 0.75 (95% confidence level 0.60-0.95, 0.51-0.90) for |u| and K^trans, respectively.

CONCLUSION

QTM postprocesses DCE-MRI automatically through deconvolution in space and time to solve the inverse problem of the transport equation. Comparing with traditional kinetics method and ECC, QTM method showed better diagnostic accuracy in differentiating benign from malignant breast tumors in this study.

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.

Flexible Doppler ultrasound device for the monitoring of blood flow velocity

Thrombosis and restenosis after vascular reconstruction procedures may cause complications such as stroke, but a clinical means to continuously monitor vascular conditions is lacking. Conventional ultrasound probes are rigid, particularly for postoperative patients with fragile skin. Techniques based on photoplethysmography or thermal analysis provide only relative changes in flow volume and have a shallow detection depth. Here, we introduce a flexible Doppler ultrasound device for the continuous monitoring of the absolute velocity of blood flow in deeply embedded arteries based on the Doppler effect. The device is thin (1 mm), lightweight (0.75 g), and skin conforming. When the dual-beam Doppler method is used, the influence of the Doppler angle on the velocity measurement is avoided. Experimental studies on ultrasound phantoms and human subjects demonstrate accurate measurement of the flow velocity. The wearable Doppler device has the potential to enhance the quality of care of patients after reconstruction surgery.

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.

Estimation of coronary stenosis severity based on flow distribution ratios

We suggest improving minimally-invasive stenosis severity estimation, using a combination of existing geometry-based methods with Transluminal Attenuation Gradient measurements. Instead of local flow values, the method uses flow distribution ratios along the entire tree. The tree geometry is used to derive a lumped model and predict the 'theoretical' ratios in each bifurcation, while attenuation measurements are used for extracting 'actual' ratios. The discrepancies between the measured and the theoretical values are utilized to assess a functional degree of stenosis. Our experimental and numerical analyses show that the quantitative value of discrepancy is proportional to stenosis severity, regardless of boundary conditions.

Engineering of magnetic nanoparticles as magnetic particle imaging tracers

Magnetic particle imaging (MPI) has recently emerged as a promising non-invasive imaging technique because of its signal linearly propotional to the tracer mass, ability to generate positive contrast, low tissue background, unlimited tissue penetration depth, and lack of ionizing radiation. The sensitivity and resolution of MPI are highly dependent on the properties of magnetic nanoparticles (MNPs), and extensive research efforts have been focused on the design and synthesis of tracers. This review examines parameters that dictate the performance of MNPs, including size, shape, composition, surface property, crystallinity, the surrounding environment, and aggregation state to provide guidance for engineering MPI tracers with better performance. Finally, we discuss applications of MPI imaging and its challenges and perspectives in clinical translation.

A novel diagnostic approach for assessing pulmonary blood flow distribution using conventional X-ray angiography

Objective

Quantitative assessment of pulmonary blood flow distribution is important when determining the clinical indications for treating pulmonary arterial branch stenosis. Lung perfusion scintigraphy is currently the gold standard for quantitative blood flow measurement. However, it is expensive, cannot provide a real-time assessment, requires additional sedation, and exposes the patient to ionizing radiation. The aim of this study was to investigate the feasibility of a novel technology for measuring pulmonary blood flow distribution in each lung by conventional X-ray pulmonary angiography and to compare its performance to that of lung perfusion scintigraphy.

Methods

Contrast-enhanced X-ray pulmonary angiography images were acquired at a frame rate of 30 frames per second. The baseline mask image, obtained before contrast agent injection, was subtracted from subsequent, consecutive images. The time-signal intensity curves of two regions of interest, established at each lung field, were obtained on a frame-to-frame basis. The net increase in signal intensity within each region at the torrent period during the second cardiac cycle before contrast agent enhancement over the total lung field was measured, and the right-to-left ratio of the signal intensity was calculated. The right-to-left ratio obtained with this approach was compared to that obtained with scintigraphy. Agreement of the right-to-left ratio between X-ray angiography and lung scintigraphy measurements was assessed using linear fitting with the Pearson correlation coefficient.

Result

The calculation of the right-to-left ratio of pulmonary blood flow by our kinetic model was feasible for seven children as a pilot study. The right-to-left ratio of pulmonary blood flow distribution calculated from pulmonary angiography was in good agreement with that of lung perfusion scintigraphy, with a Pearson correlation coefficient of 0.91 and a slope of linear fit of 1.2 (p<0.005).

Conclusion

The novel diagnostic technology using X-ray pulmonary angiography from our kinetic model can feasibly quantify the right-to-left ratio of pulmonary blood flow distribution. This technology may serve as a substitute for lung perfusion scintigraphy, which is quite beneficial for small children susceptible to radiation exposure.

High Precision Human Skin Temperature Fluctuations Measuring Instrument

This paper describes the experimental results of testing a prototype of a high precision human skin rapid temperature fluctuations measuring instrument. Based on the author’s work, an original circuit solution on a miniature semiconductor diode sensor has been designed. The proposed circuitry provides operation in the full voltage range with automatic setting and holding the operating point, as well as the necessary slope of the conversion coefficient (up to 2300 mV/°C), which makes it possible to register fast temperature oscillations from the surface of the human body and other biological objects. Simulation results in the Microcap 12 software and laboratory tests have confirmed all declared design specifications: temperature resolution of 0.01 °C, transducer thermal time constant of 0.05 s. An original thermostat and an experimental setup for the simultaneous registration of the electrocardiogram, pulse wave signals from the Biopac polygraph MP36 and a signal of temperature oscillations from the prototype thermometer have been designed for further investigations. The preliminary test results indicates that using the designed measuring instrument gives a possibility to provide an in-depth study of the relationship between micro- and macro-blood circulations manifested in skin temperature fluctuations.

A novel method for estimating the urine drainage time from the renal collecting system

PURPOSE

Partial obstruction of the upper urinary tract is a common urological pathology that leads to progressive atrophy and dysfunction of the kidney. Most methods for evaluating the urine drainage rate, to assess the severity of partial obstruction, involve injection of markers into the blood stream and therefore the filtration rate from the blood effects the drainage rate. This study presents a novel method for assessing the drainage rate from the upper urinary tract by analyzing sequential fluoroscopic images from a routine nephrostogram, in which contrast material is introduced directly into the renal collecting system.

METHODS

Fluoroscopic images from 36 nephrostograms, following percutaneous nephrolithotomy, were retrospectively evaluated, 19 with a dilated renal pelvis. A radiological model for calculating the radiopacity of the renal pelvis, which reflects the amount of contrast material in each sequential image, was developed. Using this model, an algorithm was designed for generating a drainage curve and calculating the "drainage time" t_1/2 in which half of the contrast material has drained from the renal pelvis.

RESULTS

Analysis of images of a step-wedge phantom made of an increasing number of contrast material layers showed that the calculated radiopacity of each step was proportional to the amount of contrast material, independent of the background attenuation. Analysis of the nephrostograms showed that the drainage curves highly fitted an exponential function (R = 0.961), with a significantly higher t_1/2 for dilated cases.

CONCLUSION

The developed method may be used for a quantitative and accurate estimation of the urine drainage rate.

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.

Estimating blood flow velocity using time-resolved 3D angiography and a derived physical law of contrast media

OBJECTIVE

Four-dimensional (4D) digital subtraction angiography (DSA) offers a method for evaluating hemodynamics. It is, however, unclear how the delivered contrast medium interacts with the physiological blood flow, and how hemodynamic information may be inferred from the mixture of the contrast medium and blood. In this study, we present a theoretical explanation of contrast dynamics, and an accompanying algorithm for estimating blood flow velocity.

APPROACH

We retrospectively recruited 23 patients who underwent both 4D DSA and magnetic resonance (MR) phase-contrast imaging. The 4D DSA-reconstructed contrast dynamics were first studied for the internal carotid arteries. Using physical laws governing fluid motion within a curved tube, we showed that the reconstructed contrast dynamics obeyed a simple advection equation. We then proposed an algorithm for estimating the contrast dynamics using angiographic data, and subsequently estimated the axial blood flow velocity using an advection equation.

MAIN RESULTS

The estimated velocities were compared using three techniques: the Fourier technique, Lin's method, and MR phase contrast. Testing with noise-corrupted artificial data showed that the proposed algorithm was noise resistant. The velocities of 23 patients computed by 4D DSA using the proposed algorithm showed a moderate correlation with the MR phase contrast (r = 0.61), and good correlations with the other two techniques (r = 0.75 and r = 0.72).

SIGNIFICANCE

The proposed algorithm and has been applied to blood vessel segments with poor signal-to-noise ratios and axial lengths of less than 3 cm, and has a physical basis for computing axial flow velocities using an advection equation. The results of the proposed algorithm are consistent with existing methods.

A technique for intra-procedural blood velocity quantitation using time-resolved 2D digital subtraction angiography

BACKGROUND

2D digital subtraction angiography (DSA) is utilized qualitatively to assess blood velocity changes that occur during arterial interventions. Quantitative angiographic metrics, such as blood velocity, could be used to standardize endpoints during angiographic interventions.

PURPOSE

To assess the accuracy and precision of a quantitative 2D DSA (qDSA) technique and to determine its feasibility for in vivo measurements of blood velocity.

MATERIALS AND METHODS

A quantitative DSA technique was developed to calculate intra-procedural blood velocity. In vitro validation was performed by comparing velocities from the qDSA method and an ultrasonic flow probe in a bifurcation phantom. Parameters of interest included baseline flow rate, contrast injection rate, projection angle, and magnification. In vivo qDSA analysis was completed in five different branches of the abdominal aorta in two 50 kg swine and compared to 4D Flow MRI. Linear regression, Bland-Altman, Pearson's correlation coefficient and chi squared tests were used to assess the accuracy and precision of the technique.

RESULTS

In vitro validation showed strong correlation between qDSA and flow probe velocities over a range of contrast injection and baseline flow rates (slope = 1.012, 95% CI [0.989,1.035], Pearson's r = 0.996, p < .0001). The application of projection angle and magnification corrections decreased variance to less than 5% the average baseline velocity (p = 0.999 and p = 0.956, respectively). In vivo validation showed strong correlation with a small bias between qDSA and 4D Flow MRI velocities for all five abdominopelvic arterial vessels of interest (slope = 1.01, Pearson's r = 0.880, p = <.01, Bias = 0.117 cm/s).

CONCLUSION

The proposed method allows for accurate and precise calculation of blood velocities, in near real-time, from time resolved 2D DSAs.

In vivo skin thermophysical property testing technology using flexible thermosensor-based 3ω method

Thermophysical properties of human skin surface and subsurface can reflect the hydration state of the skin and the blood flow rate in the near surface microvessels, which reveals important physiological information related to dermatology and overall health status of human body. Although a few techniques have been developed to measure these signs, complicated devices are required and the subjects need to be completely fixed during the test period. Here, a flexible thermosensor-based 3ω technology was used to monitor thermal conductivity of human skins at different states. Through the analysis of these characteristics, the corresponding physiological state can be established, which can provide a new detection method for the evaluation or prediction of human health status.

A Quantitative Digital Subtraction Angiography Technique for Characterizing Reduction in Hepatic Arterial Blood Flow During Transarterial Embolization

OBJECTIVE

There is no standardized and objective method for determining the optimal treatment endpoint (sub-stasis) during transarterial embolization. The objective of this study was to demonstrate the feasibility of using a quantitative digital subtraction angiography (qDSA) technique to characterize intra-procedural changes in hepatic arterial blood flow velocity in response to transarterial embolization in an in vivo porcine model.

MATERIALS AND METHODS

Eight domestic swine underwent bland transarterial embolizations to partial- and sub-stasis angiographic endpoints with intraprocedural DSA acquisitions. Embolized lobes were assessed on histopathology for ischemic damage and tissue embolic particle density. Analysis of target vessels used qDSA and a commercially available color-coded DSA (ccDSA) tool to calculate blood flow velocities and time-to-peak, respectively.

RESULTS

Blood flow velocities calculated using qDSA showed a statistically significant difference (p < 0.01) between partial- and sub-stasis endpoints, whereas time-to-peak calculated using ccDSA did not show a significant difference. During the course of embolizations, the average correlation with volume of particles delivered was larger for qDSA (- 0.86) than ccDSA (0.36). There was a statistically smaller mean squared error (p < 0.01) and larger coefficient of determination (p < 0.01) for qDSA compared to ccDSA. On pathology, the degree of embolization as calculated by qDSA had a moderate, positive correlation (p < 0.01) with the tissue embolic particle density of ischemic regions within the embolized lobe.

CONCLUSIONS

qDSA was able to quantitatively discriminate angiographic embolization endpoints and, compared to a commercially available ccDSA method, improve intra-procedural characterization of blood flow changes. Additionally, the qDSA endpoints correlated with tissue-level changes.

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.

Segmentation-Based Blood Flow Parameter Refinement in Cerebrovascular Structures Using 4-D Arterial Spin Labeling MRA

OBJECTIVE

Cerebrovascular diseases are one of the main global causes of death and disability in the adult population. The preferred imaging modality for the diagnostic routine is digital subtraction angiography, an invasive modality. Time-resolved three-dimensional arterial spin labeling magnetic resonance angiography (4D ASL MRA) is an alternative non-invasive modality, which captures morphological and blood flow data of the cerebrovascular system, with high spatial and temporal resolution. This work proposes advanced medical image processing methods that extract the anatomical and hemodynamic information contained in 4D ASL MRA datasets.

METHODS

A previously published segmentation method, which uses blood flow data to improve its accuracy, is extended to estimate blood flow parameters by fitting a mathematical model to the measured vascular signal. The estimated values are then refined using regression techniques within the cerebrovascular segmentation. The proposed method was evaluated using fifteen 4D ASL MRA phantoms, with ground-truth morphological and hemodynamic data, fifteen 4D ASL MRA datasets acquired from healthy volunteers, and two 4D ASL MRA datasets from patients with a stenosis.

RESULTS

The proposed method reached an average Dice similarity coefficient of 0.957 and 0.938 in the phantom and real dataset segmentation evaluations, respectively. The estimated blood flow parameter values are more similar to the ground-truth values after the refinement step, when using phantoms. A qualitative analysis showed that the refined blood flow estimation is more realistic compared to the raw hemodynamic parameters.

CONCLUSION

The proposed method can provide accurate segmentations and blood flow parameter estimations in the cerebrovascular system using 4D ASL MRA datasets.

SIGNIFICANCE

The information obtained with the proposed method can help clinicians and researchers to study the cerebrovascular system non-invasively.

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.

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.

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.

Noninvasive Assessment of Jugular Venous Pressure via Force-Coupled Single Crystal Ultrasound

OBJECTIVE

we have developed a handheld device for noninvasive quantitative assessment of jugular venous pressure (JVP).

METHODS

we used a single crystal ultrasound coupled to a force-sensing load cell to measure JVP based on the force necessary to collapse the internal jugular vein (IJV) walls. We used a gelatin-based model system of the IJV to test the ability of single crystal ultrasound to identify the IJV and verified the cross-sectional position and diameter of the vessels with conventional imaging ultrasound. We also tested our prototype device on healthy human volunteers.

RESULTS

experiments on model system demonstrated that vessel diameters determined with single crystal ultrasound were in close agreement with the diameters derived from conventional 2-D ultrasound. Proof-of-concept human experiments demonstrate that single crystal ultrasound can detect the IJV in basal and collapsed states, as compared to gold-standard sonography (insert stats). Assessment of JVP in human volunteers was physiologically consistent with and sensitive to postural changes (supine JVP 6.6 ± 2.4 mmHg; standing JVP 4.2 ± 1.9 mmHg (p < 0.0001).

CONCLUSION

noninvasive assessment of JVP could prove valuable in informing rapid clinical decision-making across various pathologies and conditions leading to derangements in intravascular volume status.

Magnetic particle imaging for in vivo blood flow velocity measurements in mice

Magnetic particle imaging (MPI) is a new imaging technology. It is a potential candidate to be used for angiographic purposes, to study perfusion and cell migration. The aim of this work was to measure velocities of the flowing blood in the inferior vena cava of mice, using MPI, and to evaluate it in comparison with magnetic resonance imaging (MRI). A phantom mimicking the flow within the inferior vena cava with velocities of up to 21 cm s^-1 was used for the evaluation of the applied analysis techniques. Time-density and distance-density analyses for bolus tracking were performed to calculate flow velocities. These findings were compared with the calibrated velocities set by a flow pump, and it can be concluded that velocities of up to 21 cm s^-1 can be measured by MPI. A time-density analysis using an arrival time estimation algorithm showed the best agreement with the preset velocities. In vivo measurements were performed in healthy FVB mice (n  =  10). MRI experiments were performed using phase contrast (PC) for velocity mapping. For MPI measurements, a standardized injection of a superparamagnetic iron oxide tracer was applied. In vivo MPI data were evaluated by a time-density analysis and compared to PC MRI. A Bland-Altman analysis revealed good agreement between the in vivo velocities acquired by MRI of 4.0  ±  1.5 cm s^-1 and those measured by MPI of 4.8  ±  1.1 cm s^-1. Magnetic particle imaging is a new tool with which to measure and quantify flow velocities. It is fast, radiation-free, and produces 3D images. It therefore offers the potential for vascular imaging.

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.

In-room assessment of intravascular velocity from time-resolved rotational angiography in patients with arteriovenous malformation: a pilot study

Background

Time-resolved rotational angiography (t-RA) enables interventionists to better comprehend complex arteriovenous malformations (AVMs), thereby facilitating endovascular treatment. However, its use in evaluating hemodynamic changes has rarely been explored.

Objective

This study uses t-RA to estimate intravascular flow in patients with AVM to compare this with flow in the normal population.

Methods

Patients with available t-RA scans were prospectively categorized into one of three groups: hemorrhagic AVM, non-hemorrhagic AVM and control. Pulsatile time–density curves (TDCs) for C1, C6 and VOIMCA were used for amplitude and velocity estimation. C1 was at the cervical internal carotid artery (ICA), 2–3 cm below the carotid canal, C6 was at the paraclinoid segment of the ICA, and VOIMCA was at the junction of the first and second segment of the middle cerebral artery (MCA). A waveform amplitude ratio was defined as (peak − trough)/trough contrast intensity. VICA was defined as the distance between C6 and C1 divided by the time required for the wave to pass, and correspondingly, the average velocity of MCA (VMCA) was defined as the distance between C6 and VOIMCA divided by the duration for the same peak to travel from C6 and VOIMCA, AVM volume was estimated by MR angiography.

Results

Amplitude ratios AC1 and AC6, and average flow velocities VICA and VMCA were significantly larger in the non-hemorrhagic group than in the control group, while the hemorrhagic AVM group was not significantly different from the controls. VICA and VMCA showed moderate to good correlations with AVM volume (r=0.51 and 0.73, respectively). VMCA (33.0±9.1) was significantly lower than VICA (41.3±13.2) in the control group, but not in the two AVM groups.

Conclusion

TDC waveform propagation derived from t-RA can quantify hemodynamic differences between AVM and the control group. t-RA provides both real-time anatomic and hemodynamic evaluation, and can thus potentially improve the interventional workflow.

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.

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

Treatment Evaluation of Flow-Limiting Stenoses of the Superficial Femoral and Popliteal Artery by Parametric Color-Coding Analysis of Digital Subtraction Angiography Series

PURPOSE

To evaluate the hemodynamic effect of percutaneous transluminal intervention (PTI) on stenosis of the superficial femoral (SFA) and popliteal arteries (PA) using time-density curves (TDCs) derived from digital subtraction angiography (DSA) series in correlation with ultrasound peak systolic velocity ratio (PSVR) and ankle brachial index (ABI).

MATERIALS AND METHODS

DSA series of SFA or PA of patients with symptomatic peripheral arterial occlusive disease was obtained with a flat-panel angiography system with intention-to-treat. In DSA series acquired before and after PTI, TDCs were analyzed proximal and distal of each stenosis using parametric color coding (PCC). For correlation, ABI and PSVR measurements pre- and post-PTI were recorded for all patients.

RESULTS

In total, 25 stenoses of the SFA or PA were treated by PTI in 22 patients (17 male, 5 female, mean age 68 years). After treatment, peak-to-peak (PTP) times between TDCs proximal and distal to the treated vessel segment decreased statistically significantly (p = 0.01) on average from PTP = 1.9 ± 1.7 s to mean PTP = 1 ± 1 s. ABI and PSVR also changed statistically significantly after treatment (pretreatment ABI = 0.7 ± 0.2, PSVR = 4.2 ± 1.9; post-ABI = 0.9 ± 0.2, PSVR = 1.3 ± 0.4, both p < 0.05). Correlation parameters did not show a strong correlation between change in TDC and clinical parameters ABI and PSVR.

CONCLUSION

Using PCC for analyzing contrast medium dynamics in DSA series is clinically useful for evaluating stenoses of the SFA and PA and for immediate treatment control after PTA.

LEVEL OF EVIDENCE

Case series, IV.

Vagal denervation inhibits the increase in pulmonary blood flow during partial lung aeration at birth

KEY POINTS

Lung aeration at birth significantly increases pulmonary blood flow, which is unrelated to increased oxygenation or other spatial relationships that match ventilation to perfusion. Using simultaneous X-ray imaging and angiography in near-term rabbits, we investigated the relative contributions of the vagus nerve and oxygenation to the increase in pulmonary blood flow at birth. Vagal denervation inhibited the global increase in pulmonary blood flow induced by partial lung aeration, although high inspired oxygen concentrations can partially mitigate this effect. The results of the present study indicate that a vagal reflex may mediate a rapid global increase in pulmonary blood flow in response to partial lung aeration.

ABSTRACT

Air entry into the lungs at birth triggers major cardiovascular changes, including a large increase in pulmonary blood flow (PBF) that is not spatially related to regional lung aeration. To investigate the possible underlying role of a vagally-mediated stimulus, we used simultaneous phase-contrast X-ray imaging and angiography in near-term (30 days of gestation) vagotomized (n = 15) or sham-operated (n = 15) rabbit kittens. Rabbits were imaged before ventilation, when one lung was ventilated (unilateral) with 100% nitrogen (N₂ ), air or 100% oxygen (O₂ ), and after all kittens were switched to unilateral ventilation in air and then ventilation of both lungs using air. Compared to control kittens, vagotomized kittens had little or no increase in PBF in both lungs following unilateral ventilation when ventilation occurred with 100% N₂ or with air. However, relative PBF did increase in vagotomized animals ventilated with 100% O₂ , indicating the independent stimulatory effects of local oxygen concentration and autonomic innervation on the changes in PBF at birth. These findings demonstrate that vagal denervation inhibits the previously observed increase in PBF with partial lung aeration, although high inspired oxygen concentrations can partially mitigate this effect.

Computational Fluid Dynamics modeling of contrast transport in basilar aneurysms following flow-altering surgeries

In vivo measurement of blood velocity fields and flow descriptors remains challenging due to image artifacts and limited resolution of current imaging methods; however, in vivo imaging data can be used to inform and validate patient-specific computational fluid dynamics (CFD) models. Image-based CFD can be particularly useful for planning surgical interventions in complicated cases such as fusiform aneurysms of the basilar artery, where it is crucial to alter pathological hemodynamics while preserving flow to the distal vasculature. In this study, patient-specific CFD modeling was conducted for two basilar aneurysm patients considered for surgical treatment. In addition to velocity fields, transport of contrast agent was simulated for the preoperative and postoperative conditions using two approaches. The transport of a virtual contrast passively following the flow streamlines was simulated to predict post-surgical flow regions prone to thrombus deposition. In addition, the transport of a mixture of blood with an iodine-based contrast agent was modeled to compare and verify the CFD results with X-ray angiograms. The CFD-predicted patterns of contrast flow were qualitatively compared to in vivo X-ray angiograms acquired before and after the intervention. The results suggest that the mixture modeling approach, accounting for the flow rates and properties of the contrast injection, is in better agreement with the X-ray angiography data. The virtual contrast modeling assessed the residence time based on flow patterns unaffected by the injection procedure, which makes the virtual contrast modeling approach better suited for prediction of thrombus deposition, which is not limited to the peri-procedural state.