Intra-vascular blood velocity and volumetric flow rate calculated from dynamic 4D CT angiography using a time of flight technique

Reconstruction of blood flow velocity with deep learning information fusion from spectral ct projections and vessel geometry

In this work, we investigate a new deep learning reconstruction method of blood flow velocity within deformed vessels from contrast enhanced X-ray projections and vessel geometry. The principle of the method is to perform linear or nonlinear dimension reductions on the Radon projections and on the mesh of the vessel. These low dimensional projections are then fused to obtain the velocity field in the vessel. The accuracy of the reconstruction method is proved using various neural network architectures with realistic unsteady blood flows. The approach leverages the vessel geometry information and outperforms the simple PCA-net.

Deep learning methods for blood flow reconstruction in a vessel with contrast enhanced x-ray computed tomography

The reconstruction of blood velocity in a vessel from contrast enhanced x-ray computed tomography projections is a complex inverse problem. It can be formulated as reconstruction problem with a partial differential equation constraint. A solution can be estimated with the a variational adjoint method and proper orthogonal decomposition (POD) basis. In this work, we investigate new inversion approaches based on PODs coupled with deep learning methods. The effectiveness of the reconstruction methods is shown with simulated realistic stationary blood flows in a vessel. The methods outperform the reduced adjoint method and show large speed-up at the online stage.

Quantitative hemodynamic assessment of stenotic below-the-knee arteries using spatio-temporal bolus tracking on 4D-CT angiography

BACKGROUND

Peripheral arterial disease (PAD) is a chronic occlusive disease that restricts blood flow in the lower limbs, causing partial or complete blockages of the blood flow. While digital subtraction angiography (DSA) has traditionally been the preferred method for assessing blood flow in the lower limbs, advancements in wide beam Computed Tomography (CT), allowing successive acquisition at high frame rate, might enable hemodynamic measurements.

PURPOSE

To quantify the arterial blood flow in stenotic below-the-knee (BTK) arteries. To this end, we propose a novel method for contrast bolus tracking and assessment of quantitative hemodynamic parameters in stenotic arteries using 4D-CT.

METHODS

Fifty patients with suspected PAD underwent 4D-CT angiography in addition to the clinical run-off computed tomography angiography (CTA). From these dynamic acquisitions, the BTK arteries were segmented and the region of maximum blood flow was extracted. Time attenuation curves (TAC) were estimated using 2D spatio-temporal B-spline regression, enforcing both spatial and temporal smoothness. From these curves, quantitative hemodynamic parameters, describing the shape of the propagating contrast bolus were automatically extracted. We evaluated the robustness of the proposed TAC fitting method with respect to interphase delay and imaging noise and compared it to commonly used approaches. Finally, to illustrate the potential value of 4D-CT, we assessed the correlation between the obtained hemodynamic parameters and the presence of PAD.

RESULTS

280 out of 292 arteries were successfully segmented, with failures mainly due to a delayed contrast arrival. The proposed method led to physiologically plausible hemodynamic parameters and was significantly more robust compared to 1D temporal regression. A significant correlation between the presence of proximal stenoses and several hemodynamic parameters was found.

CONCLUSIONS

The proposed method based on spatio-temporal bolus tracking was shown to lead to stable and physiologically plausible estimation of quantitative hemodynamic parameters, even in the case of stenotic arteries. These parameters may provide valuable information in the evaluation of PAD and contribute to its diagnosis.

A Noninvasive Assessment of Flow Based on Contrast Dispersion in Computed Tomography Angiography: A Computational and Experimental Phantom Study

Transluminal attenuation gradient (TAG), defined as the gradient of the contrast agent attenuation drop along the vessel, is an imaging biomarker that indicates stenosis in the coronary arteries. The transluminal attenuation flow encoding (TAFE) equation is a theoretical platform that quantifies blood flow in each coronary artery based on computed tomography angiography (CTA) imaging. This formulation couples TAG (i.e., contrast dispersion along the vessel) with fluid dynamics. However, this theoretical concept has never been validated experimentally. The aim of this proof-of-principle phantom study is to validate TAFE based on CTA imaging. Dynamic CTA images were acquired every 0.5 s. The average TAFE estimated flow rates were compared against four predefined pump values in a straight (20, 25, 30, 35, and 40 ml/min) and a tapered phantom (25, 35, 45, and 55 ml/min). Using the TAFE formulation with no correction, the flow rates were underestimated by 33% and 81% in the straight and tapered phantoms, respectively. The TAFE formulation was corrected for imaging artifacts focusing on partial volume averaging and radial variation of contrast enhancement. After corrections, the flow rates estimated in the straight and tapered phantoms had an excellent Pearson correlation of r = 0.99 and 0.87 (p < 0.001), respectively, with only a 0.6%±0.2 mL/min difference in estimation of the flow rate. In this proof-of-concept phantom study, we corrected the TAFE formulation and showed a good agreement with the actual pump values. Future clinical validations are needed for feasibility of TAFE in clinical use.

Angiographic Pulse Wave Coherence in the Human Brain

A stroke volume of arterial blood that arrives to the brain housed in the rigid cranium must be matched over the cardiac cycle by an equivalent volume of ejected venous blood. We hypothesize that the brain maintains this equilibrium by organizing coherent arterial and venous pulse waves. To test this hypothesis, we applied wavelet computational methods to diagnostic cerebral angiograms in four human patients, permitting the capture and analysis of cardiac frequency phenomena from fluoroscopic images acquired at faster than cardiac rate. We found that the cardiac frequency reciprocal phase of a small region of interest (ROI) in a named artery predicts venous anatomy pixel-wise and that the predicted pixels reconstitute venous bolus passage timing. Likewise, a small ROI in a named vein predicts arterial anatomy and arterial bolus passage timing. The predicted arterial and venous pixel groups maintain phase complementarity across the bolus travel. We thus establish a novel computational method to analyze vascular pulse waves from minimally invasive cerebral angiograms and provide the first direct evidence of arteriovenous coupling in the intact human brain. This phenomenon of arteriovenous coupling may be a physiologic mechanism for how the brain precisely maintains mechanical equilibrium against volume displacement and kinetic energy transfer resulting from cyclical deformations with each heartbeat. The study also paves the way to study deranged arteriovenous coupling as an underappreciated pathophysiologic disturbance in a myriad of neurological pathologies linked by mechanical disequilibrium.

Reconstruction of vascular blood flow in a vessel from tomographic projections

In this work, we study the measurement of blood velocity with contrast enhanced computed tomography. The reconstruction is based on CT projections perpendicular to the main axis of the vessel and on a partial differential equation describing the propagation of the contrast agent. The inverse problem is formulated as an optimal control problem with the transport equation as constraint. The velocity field is obtained with stationary and unstationary Navier-Stokes equations and it is reconstructed with the adjoint method. The velocity and the density of the contrast agent are well reconstructed. The reconstruction results obtained are better for the axial component of the velocity than for transverse components.

Dynamic modes of inflow jet in brain aneurysms

The transition of the inflow jet to turbulence is crucial in understanding the pathology of brain aneurysms. Previous works Le et al. (2010, 2013) have shown evidence for a highly dynamic inflow jet in the ostium of brain aneurysms. While it is highly desired to investigate this inflow jet dynamics in clinical practice, the constraints on spatial and temporal resolutions of in vivo data do not allow a detailed analysis of this transition. In this work, Dynamic Mode Decomposition (DMD) is used to identify the most energetic modes of the inflow jet in patient-specific models of internal carotid aneurysms via the utilization of high-resolution simulation data. It is hypothesized that dynamic modes are not solely controlled by the blood flow waveform at the parent artery. They are also dependent on jet-wall interaction phenomena. DMD analysis shows that the spatial extent of low- frequency modes corresponds well to the most energetic areas of the inflow jet. The high-frequency modes are short-lived and correspond to the flow separation at the proximal neck and the jet's impingement onto the aneurysmal wall. Low-frequency modes can be reconstructed at relatively low spatial and temporal resolutions comparable to ones of in vivo data. The current results suggest that DMD can be practically useful in analyzing blood flow patterns of brain aneurysms with in vivo data.

Four-dimensional computed tomography angiography analysis of internal carotid arteries opacification at the skull base to detect delayed cerebral ischemia: a feasibility study

PURPOSE

Delayed cerebral ischemia represents a significant cause of poor functional outcome for patients with vasospasm after subarachnoid hemorrhage. We investigated whether delayed cerebral ischemia could be detected by the arterial opacification of internal carotid artery at the level of the skull base.

METHODS

In this exploratory, nested retrospective cohort diagnostic accuracy study, patients with clinical and/or transcranial Doppler suspicion of vasospasm who underwent four-dimensional computed tomography angiography were included. They were split into two groups for the main endpoint analysis, according to the actually adopted morphological (cerebral infarction) and clinical criteria (neurologic deterioration) of delayed cerebral ischemia. Opacification with a temporal resolution of 0.15 s of both internal carotid arteries at the skull base level was obtained through a semi-automated segmentation method based on skeletonization, and analyzed by a wavelet transform (rbio2.2, level 1). The results obtained by k-means clustering were analyzed with regard to the state of delayed cerebral infarction.

RESULTS

Over ten patients included and analyzed, five patients presented a delayed cerebral ischemia, two of them in both side. The semi-automated processing and analysis clustered two different types of opacification curves. The obtaining of a nonlinear opacification pattern was associated (p < 0.001) with delayed cerebral ischemia.

CONCLUSIONS

The analysis of arterial opacification of internal carotid arteries at skull base by the proposed processing is feasible and leads to cluster two types of opacification that may help to early detect and prevent delayed cerebral ischemia, in particularly when examinations are artifacted by aneurysm treatment materials.

Determination of renal function and injury using near-infrared fluorimetry in experimental cardiorenal syndrome

Cardiorenal syndrome type 1 causes acute kidney injury but is poorly understood; animal models and diagnostic aids are lacking. Robust noninvasive measurements of glomerular filtration rate are required for injury models and clinical use. Several have been described but are untested in translational models and suffer from biologic interference. We developed a mouse model of cardiorenal syndrome and tested the novel near-infrared fluorophore ZW800-1 to assess renal and cardiac function. We performed murine cardiac arrest and cardiopulmonary resuscitation followed by transthoracic echocardiography, 2 and 24 h later. Transcutaneous fluorescence of ZW800-1 bolus dispersion and clearance was assessed with whole animal imaging and compared with glomerular filtration rate (GFR; inulin clearance), tubular cell death (using unbiased stereology), and serum creatinine. Correlation, Bland-Altman, and polar analyses were used to compare GFR with ZW800-1 clearance. Cardiac arrest and cardiopulmonary resuscitation caused reversible cardiac failure, halving fractional shortening of the left ventricle (n = 12, P = 0.03). Acute kidney injury resulted with near-zero GFR and sixfold increase in serum creatinine 24 h later (n = 16, P < 0.01). ZW800-1 biodistribution and clearance were exclusively renal. ZW800-1 t_1/2 and clearance correlated with GFR (r = 0.92, n = 31, P < 0.0001). ZW800-1 fluorescence was reduced in cardiac arrest, and cardiopulmonary resuscitation-treated mice compared with sham animals 810 s after injection (P < 0.01) and bolus time-dispersion curves demonstrated that ZW800-1 fluorescence dispersion correlated with left ventricular function (r = 0.74, P < 0.01). Cardiac arrest and cardiopulmonary resuscitation lead to experimental cardiorenal syndrome type 1. ZW800-1, a small near-infrared fluorophore being developed for clinical intraoperative imaging, is favorable for evaluating cardiac and renal function noninvasively.

Contrast Gradient-Based Blood Velocimetry With Computed Tomography: Theory, Simulations, and Proof of Principle in a Dynamic Flow Phantom

OBJECTIVES

The aim of this study was to introduce a new theoretical framework describing the relationship between the blood velocity, computed tomography (CT) acquisition velocity, and iodine contrast enhancement in CT images, and give a proof of principle of contrast gradient-based blood velocimetry with CT.

MATERIALS AND METHODS

The time-averaged blood velocity (v(blood)) inside an artery along the axis of rotation (z axis) is described as the mathematical division of a temporal (Hounsfield unit/second) and spatial (Hounsfield unit/centimeter) iodine contrast gradient. From this new theoretical framework, multiple strategies for calculating the time-averaged blood velocity from existing clinical CT scan protocols are derived, and contrast gradient-based blood velocimetry was introduced as a new method that can calculate v(blood) directly from contrast agent gradients and the changes therein. Exemplarily, the behavior of this new method was simulated for image acquisition with an adaptive 4-dimensional spiral mode consisting of repeated spiral acquisitions with alternating scan direction. In a dynamic flow phantom with flow velocities between 5.1 and 21.2 cm/s, the same acquisition mode was used to validate the simulations and give a proof of principle of contrast gradient-based blood velocimetry in a straight cylinder of 2.5 cm diameter, representing the aorta.

RESULTS

In general, scanning with the direction of blood flow results in decreased and scanning against the flow in increased temporal contrast agent gradients. Velocity quantification becomes better for low blood and high acquisition speeds because the deviation of the measured contrast agent gradient from the temporal gradient will increase. In the dynamic flow phantom, a modulation of the enhancement curve, and thus alternation of the contrast agent gradients, can be observed for the adaptive 4-dimensional spiral mode and is in agreement with the simulations. The measured flow velocities in the downslopes of the enhancement curves were in good agreement with the expected values, although the accuracy and precision worsened with increasing flow velocities.

CONCLUSIONS

The new theoretical framework increases the understanding of the relationship between the blood velocity, CT acquisition velocity, and iodine contrast enhancement in CT images, and it interconnects existing blood velocimetry methods with research on transluminary attenuation gradients. With these new insights, novel strategies for CT blood velocimetry, such as the contrast gradient-based method presented in this article, may be developed.