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.

Application of 1,000 fps High-Speed Angiography to In-Vitro Hemodynamic Evaluation of Left Ventricular Assist Device Outflow Graft Configurations

Left ventricular assist device (LVAD)-induced hemodynamics are characterized by fast-moving flow with large variations in velocity, making quantitative assessments difficult with existing imaging methods. This study demonstrates the ability of 1,000 fps high-speed angiography (HSA) to quantify the effect of the surgical implantation angle of a LVAD outflow graft on the hemodynamics within the ascending aorta in vitro . High-speed angiography was performed on patient-derived, three-dimensional-printed optically opaque aortic models using a nonsoluble contrast media, ethiodol, as a flow tracer. Outflow graft configuration angles of 45° and 90° with respect to the central aortic axis were considered. Projected velocity distributions were calculated from the high-speed experimental sequences using two methods: a physics-based optical flow algorithm and tracking of radio-opaque particles. Particle trajectories were also used to evaluate accumulated shear stress. Results were then compared with computational fluid dynamics (CFD) simulations to confirm the results of the high-speed imaging method. Flow patterns derived from HSA coincided with the impingement regions and recirculation zones formed in the aortic root as seen in the CFD for both graft configurations. Compared with the 45° graft, the 90° configuration resulted in 81% higher two-dimensional-projected velocities (over 100 cm/s) along the contralateral wall of the aorta. Both graft configurations suggest elevated accumulated shear stresses along individual trajectories. Compared with CFD simulations, HSA successfully characterized the fast-moving flow and hemodynamics in each LVAD graft configuration in vitro , demonstrating the potential utility of this technology as a quantitative imaging modality.

A Prototype Software System for Intra-procedural Staff Dose Monitoring and Virtual Reality Training for Fluoroscopically Guided Interventional Procedures

Staff dose management in fluoroscopically guided interventional procedures is a continuing problem. The scattered radiation display system (SDS), which our group has developed, provides in-room visual feedback of scatter dose to staff members during fluoroscopically guided interventional (FGI) procedures as well as extra-procedure staff and resident training. There have been a number of virtual safety training systems developed that provide detailed feedback for staff, but utilize expensive graphics processing units (GPUs) and dosimeter systems, or interaction with the x-ray system in a manner which entails additional radiation exposure and is not compatible with the As Low as Reasonably Achievable paradigm. The SDS, on the other hand, incorporates a library of look-up-table (LUT) room scatter distributions determined using the EGSnrc Monte Carlo software, which facilitates accurate and rapid system update without the need for GPUs. Real-time display of these distributions is provided for feedback to staff during a procedure. After a procedure is completed, machine parameter and staff position log files are stored, retaining all of the exposure and geometric information for future review. A graphic user interface (GUI) in Unity3D enables procedure playback and interactive virtual-reality (VR) staff and resident training with virtual control of exposure conditions using an Oculus headset and controller. Improved staff and resident awareness using this system should lead to increased safety and reduced occupational dose.

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

Purpose

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

Approach

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

Results

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

Conclusions

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

Determining 3D Distributions of Pulsatile Blood Flow Using Orthogonal Simultaneous Biplane High-Speed Angiography (SB-HSA) with 1000 fps CdTe Photon Counting Detectors for 3D X-ray Particle Image Velocimetry (3D-XPIV) compared to Results Using Computational Fluid Dynamics (CFD)

3D hemodynamic distributions are useful for the diagnosis and treatment of aneurysms. Detailed blood-flow patterns and derived velocity maps can be obtained using 1000 fps High Speed Angiography (HSA). The novel orthogonal Simultaneous Biplane High-Speed Angiography (SB-HSA) system enables flow information to be quantified in multiple planes, and with additional components of flow at depth, accurate 3D flow distributions are available. Computational Fluid Dynamics (CFD) is the current standard for derivation of volumetric flow distributions, but obtaining solution convergence is computationally expensive and time intensive. More importantly, matching in-vivo boundary conditions is non-trivial. Therefore, an experimentally derived 3D flow distribution method could offer realistic results with less computation time. Using SB-HSA image sequences, 3D X-Ray Particle Image Velocimetry (3D-XPIV) was explored as a new method for assessing 3D flow. 3D-XPIV was demonstrated using an in-vitro setup, where a patient-specific internal carotid artery aneurysm model was attached to a flow loop, and an automated injection of iodinated microspheres was used as a flow tracer. Two 1000 fps photon-counting detectors were placed orthogonally with the aneurysm model in the FOV of both planes. Frame-synchronization of the two detectors made correlation of single-particle velocity components at a given timepoint possible. With frame-rates of 1000 fps, small particle displacements between frames resolved realistic time varying flow, where accurate velocity distributions depended on near-instantaneous velocities. 3D-XPIV velocity distributions were compared to CFD velocity distributions, where the simulation boundary conditions matched the in-vitro setup. Results showed similar velocity distributions between CFD and 3D-XPIV.

Multi-angled simultaneous biplane High-Speed Angiography (HSA) of patient-specific 3D-printed aneurysm phantoms using 1000 fps CdTe Photon-Counting Detectors (PCD's)

1000 fps HSA enables visualization of flow details, which may be important in accurately guiding interventional procedures; however, single-plane imaging may lack clear visualization of vessel geometry and flow detail. The previously presented high-speed orthogonal biplane imaging may overcome these limitations but may still result in foreshortening of vessel morphology. In certain morphologies, acquiring two non-orthogonal biplane projections at multiple angles can provide better flow detail rather than a standard orthogonal biplane acquisition. Flow studies of aneurysm models were performed, where simultaneous biplane acquisitions at various angles separating the two detector views allowed for better evaluation of morphology and flow. 3D-printed, patient-specific internal carotid artery aneurysm models were imaged with various non-orthogonal angles between the two high-speed photon-counting detectors (7.5 cm x 5 cm FOV) to provide frame-correlated simultaneous 1000-fps image sequences. Fluid dynamics were visualized in multi-angled planes of each model using automated injections of iodine contrast media. The resulting dual simultaneous frame-correlated 1000-fps acquisitions from multiple planes of each aneurysm model provided improved visualization of complex aneurysm geometries and flow streamlines. Multi-angled biplane acquisitions with frame correlation allows for further understanding of aneurysm morphology and flow details: additionally, the ability to recover fluid dynamics at depth enables accurate analysis of 3D flow streamlines, and it is expected that multiple-planar views will enable better volumetric flow visualization and quantification. Such better visualization has the potential to improve interventional procedures.

Use of high-speed angiography HSA-derived boundary conditions and Physics Informed Neural Networks (PINNs) for comprehensive estimation of neurovascular hemodynamics

Purpose

Physics-informed neural networks (PINNs) and computational fluid dynamics (CFD) have both demonstrated an ability to derive accurate hemodynamics if boundary conditions (BCs) are known. Unfortunately, patient-specific BCs are often unknown, and assumptions based upon previous investigations are used instead. High speed angiography (HSA) may allow extraction of these BCs due to the high temporal fidelity of the modality. We propose to investigate whether PINNs using convection and Navier-Stokes equations with BCs derived from HSA data may allow for extraction of accurate hemodynamics in the vasculature.

Materials and Methods

Imaging data generated from in vitro 1000 fps HSA, as well as simulated 1000 fps angiograms generated using CFD were utilized for this study. Calculations were performed on a 3D lattice comprised of 2D projections temporally stacked over the angiographic sequence. A PINN based on an objective function comprised of the Navier-Stokes equation, the convection equation, and angiography-based BCs was used for estimation of velocity, pressure and contrast flow at every point in the lattice.

Results

Imaging-based PINNs show an ability to capture such hemodynamic phenomena as vortices in aneurysms and regions of rapid transience, such as outlet vessel blood flow within a carotid artery bifurcation phantom. These networks work best with small solution spaces and high temporal resolution of the input angiographic data, meaning HSA image sequences represent an ideal medium for such solution spaces.

Conclusions

The study shows the feasibility of obtaining patient-specific velocity and pressure fields using an assumption-free data driven approach based purely on governing physical equations and imaging data.

2D versus 3D comparison of angiographic imaging biomarkers using computational fluid dynamics simulations of contrast injections

Quantitative angiography (QAngio) may provide hemodynamic information during neurointerventional procedures through imaging biomarkers related to contrast flow. The standard clinical implementation of QAngio is limited by projection imaging: analysis of contrast motion within complex 3D geometries is restricted to 1-2 projection views, truncating the potential wealth of imaging biomarkers related to disease progression or efficacy of treatment. To understand the limitations of 2D biomarkers, we propose the use of in-silico contrast distributions to investigate the potential benefits of 3D-QAngio within the context of neurovascular hemodynamics. Ground-truth in-silico contrast distributions were generated in two patient-specific intracranial aneurysm models, accounting for the physical interactions of contrast media and blood. A short bolus of contrast was utilized to obtain full a wash-in/ wash-out cycle within the aneurysm ROI. Simulated angiograms mimicking clinical cone-beam CT (CBCT) acquisitions were then generated, and volumetric contrast distributions were reconstructed to analyze bulk contrast flow. The ground-truth 3D-CFD, reconstructed 3D-CBCT-DSA, and 2D-DSA projections were used to extract QAngio parameters related to contrast time dilution curves, such as area under the curve (AUC), peak height (PH), mean-transit-time (MTT), time-to-peak (TTP), and time to arrival (TTA). An initial comparison of quantitative flow parameters in both 2D and 3D, in a smaller and larger aneurysm, indicated that 3D-QAngio can provide a good description of bulk flow characteristics (TTA, TTP, MTT), but recovery of integral parameters (PH, AUC) aneurysms is limited. Nonetheless, incorporation of 3D-QAngio methods may provide additional insight into our understanding of abnormal vascular flow patterns.

Eye-Lens Dose Reduction using Region of Interest (ROI) Attenuators in Neuroimaging

Lens dose can be high during neuro-interventional procedures, increasing the risk of cataractogenesis. Although beam collimation can be effective in reducing lens dose, it also restricts the FOV. ROI imaging with a reduced-dose peripheral field permits full-field information with reduced lens dose. This work investigates the magnitude of lens-dose reduction possible with ROI imaging. EGSnrc Monte-Carlo calculations of lens dose were made for the Zubal head phantom as a function of gantry angulation and head shift from isocenter for both large and small FOV's. The lens dose for ROI attenuators of varying transmission was simulated as the weighted sum of the lens dose from the small ROI FOV and that from the attenuated larger FOV. Image intensity and quantum mottle differences between ROI and periphery can be equalized by image processing. The lens dose varies considerably with beam angle, head shift, and field size. For both eyes, the lens-dose reduction with an ROI attenuator increases with LAO angulation, being highest for lateral projections and lowest for PA. For an attenuator with small ROI field (5 × 5 cm) and 20% transmission, the lens dose for lateral projections is reduced by about 75% compared to a full dose 10 ×10 cm FOV, while the reduction ranges between 30 and 40% for PA projections. Use of ROI attenuators can substantially reduce the dose to the lens of the eye for all gantry angles and head shifts, while allowing peripheral information to be seen in a larger FOV.

Comparison of pulsatile flow dynamics before and after endovascular intervention in 3D-printed patient-specific internal carotid artery aneurysm models using 1000 fps photon-counting detectors for Simultaneous Biplane High Speed Angiography (SB-HSA)

A significant challenge regarding the treatment of aneurysms is the variability in morphology and analysis of abnormal flow. With conventional DSA, low frame rates limit the flow information available to clinicians at the time of the vascular intervention. With 1000 fps High-Speed Angiography (HSA), high frame rates enable flow details to be better resolved for endovascular interventional guidance. The purpose of this work is to demonstrate how 1000 fps biplane-HSA can be used to differentiate flow features, such as vortex formation and endoleaks, amongst patient-specific internal carotid artery aneurysm phantoms pre- and post-endovascular intervention using an in-vitro flow setup. The aneurysm phantoms were attached to a flow loop configured to a carotid waveform, with automated injections of contrast media. Simultaneous Biplane High-Speed Angiographic (SB- HSA) acquisitions were obtained at 1000 fps using two photon-counting detectors with the respective aneurysm and inflow/ outflow vasculature in the FOV. After x-rays were turned on, the detector acquisitions occurred simultaneously, during which iodine contrast was injected at a continuous rate. A pipeline stent was then deployed to divert flow from the aneurysm, and image sequences were once again acquired using the same parameters. Optical Flow, an algorithm that calculates velocity based on spatial-temporal intensity changes between pixels, was used to derive velocity distributions from HSA image sequences. Both the image sequences and velocity distributions indicate detailed changes in flow features amongst the aneurysms before and after deployment of the interventional device. SB-HSA can provide detailed flow analysis, including streamline and velocity changes, which may be beneficial for interventional guidance.

Investigating Angiographic Injection Parameters for Cerebral Aneurysm Hemodynamic Characterization Using Patient-Specific Simulated Angiograms

Cerebral aneurysm (CA) rupture is one of the major causes of hemorrhagic stroke. During endovascular therapy (ET), neurointerventionalists rely on qualitative image sequences and do not have access to crucial quantitative hemodynamic information. Quantifying angiographic image sequences can provide vital information, but it is not possible to perform this in a controlled manner in vivo. Computational fluid dynamics (CFD) is a valuable tool capable of providing high fidelity quantitative data by replicating the blood flow physics within the cerebrovasculature. In this work, we use simulated angiograms (SA) to quantify the hemodynamic interaction with a clinically utilized contrast agent. SA enables extraction of time density curves (TDC) within the desired region of interest to analyze hemodynamic parameters such as time to peak (TTP) and mean transit time (MTT) within the aneurysm. We present on the quantification of several hemodynamic parameters of interest for multiple, clinically-relevant scenarios such as variable contrast injection duration and bolus volumes for 7 patient-specific CA geometries. Results indicate that utilizing these analyses provides valuable hemodynamic information relating vascular and aneurysm morphology, contrast flow conditions and injection variability. The injected contrast circulates for multiple cardiac cycles within the aneurysmal region, especially for larger aneurysms and tortuous vasculature. The SA approach enables determination of angiographic parameters for each scenario. Together, these have the potential to overcome the existing barriers in quantifying angiographic procedures in vitro or in vivo, and can provide clinically valuable hemodynamic insights for CA treatment.

Geometrically independent contrast dilution gradient (CDG) velocimetry using photon-counting 1000 fps High Speed Angiography (HSA) for 2D velocity distribution estimation

Purpose

Previous studies have demonstrated the efficacy of contrast dilution gradient (CDG) analysis in determining large vessel velocity distributions from 1000 fps high-speed angiography (HSA). However, the method required vessel centerline extraction, which made it applicable only to non-tortuous geometries using a highly specific contrast injection technique. This study seeks to remove the need for a priori knowledge regarding the direction of flow and modify the vessel sampling method to make the algorithm more robust to non-linear geometries.

Materials and Methods

1000 fps HSA acquisitions were obtained in vitro with a benchtop flow loop using the XC-Actaeon (Varex Inc.) photon-counting detector, and in silico using a passive-scalar transport model within a computational fluid dynamics (CFD) simulation. CDG analyses were obtained using gridline sampling across the vessel, and subsequent 1D velocity measurement in both the x- and y-directions. The velocity magnitudes derived from the component CDG velocity vectors were aligned with CFD results via co-registration of the resulting velocity maps and compared using mean absolute percent error (MAPE) between pixels values in each method after temporal averaging of the 1-ms velocity distributions.

Results

Regions well-saturated with contrast throughout the acquisition showed agreement when compared to CFD (MAPE of 18% for the carotid bifurcation inlet and MAPE of 27% for the internal carotid aneurysm), with respective completion times of 137 seconds and 5.8 seconds.

Conclusions

CDG may be used to obtain velocity distributions in and surrounding vascular pathologies provided the contrast injection is sufficient to provide a gradient, and diffusion of contrast through the system is negligible.

Semi-automatic Co-Registration of 3D CFD Vascular Geometry to 1000 FPS High-Speed Angiographic (HSA) Projection Images for Flow Determination Comparisons

Image co-registration is an important tool that is commonly used to quantitatively or qualitatively compare information from images or data sets that vary in time, origin, etc. This research proposes a method for the semi-automatic co-registration of the 3D vascular geometry of an intracranial aneurysm to novel high-speed angiographic (HSA) 1000 fps projection images. Using the software Tecplot 360, 3D velocimetry data generated from computational fluid dynamics (CFD) for patient-specific vasculature models can be extracted and uploaded into Python. Dilation, translation, and angular rotation of the 3D velocimetry data can then be performed in order to co-register its geometry to corresponding 2D HSA projection images of the 3D printed vascular model. Once the 3D CFD velocimetry data is geometrically aligned, a 2D velocimetry plot can be generated and the Sørensen-Dice coefficient can be calculated in order to determine the success of the co-registration process. The co-registration process was performed ten times for two different vascular models and had an average Sørensen-Dice coefficient of 0.84 ± 0.02. The method presented in this research allows for a direct comparison between 3D CFD velocimetry data and in-vitro 2D velocimetry methods. From the 3D CFD, we can compare various flow characteristics in addition to velocimetry data with HSA-derived flow metrics. The method is robust to other vascular geometries as well.

Use of 1000fps High Speed X-ray Angiography (HSAngio) to quantify differences in flow diversion effects of three stents with different coverage densities in a cerebral aneurysm invitro model

High temporal resolution images acquired using 1000fps HSAngio can be used to visualize blood flow patterns and derive flow velocities during neurointerventional procedures. In this work we use this technology to quantify the changes in the blood flow velocities inside a cerebral aneurysm after treatment with three different stents with varying degrees of metal coverage density; stent A : <2%, stent B: 23% and stent C: 40%. A 3D printed in-vitro model of internal carotid artery aneurysm was connected to a flow loop (60% water, 40% glycerol solution used as circulation fluid, circulation flow rate 8 L/s). An automatic programmable injector (KD Scientific Legato 110) was used to inject iodine contrast agent at a rate of 88 mL/min in 3secs. 1000 fps HSAngio sequences of the contrast injection were acquired using an Aries single photon counting detector (Direct Conversion Inc., Stockholm). From these images blood flow velocities were calculated using an optical flow algorithm. As expected the biggest reduction in blood flow velocity inside the aneurysm was 32.4% after deployment of stent C. However, the velocity profile distribution indicated there was still a significant inflow jet into the aneurysm which could be caused by a endoluminal leak between the stent and the vessel wall. The average reduction was only 14% after placement of stent B and 3% after placement of stent A. Blood velocity distribution maps derived using 1000fps HSAngiography technology can be used to evaluate the quality of flow diversion within the aneurysm after placement of stent. Critical information such as endo luminal leakage which can cause treatment failure can also be detected.

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

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

Initial evaluation of 2D and 3D simulated high-speed 1000 fps vascular contrast-flow image sequences using computational fluid dynamics (CFD)

Digital subtraction angiography (DSA) remains the clinical standard for detailed visualization of the neurovasculature due to its high-spatial resolution; however, detailed blood-flow quantification is impaired by its low-temporal resolution. 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. We have implemented a physics-based optical flow method to extract such information from HSA, but validation of the angiography-derived velocity distributions is not straightforward. Computational fluid dynamics (CFD) is widely regarded as the benchmark for hemodynamic analysis, as it provides a multitude of quantitative flow parameters throughout the volume of interest. However, there are several limitations with this method related to over-simplification of boundary conditions and suboptimal meshing (spatial resolution), that make CFD simulation results an inexact criterion for validation. To overcome this issue for HSA validation, CFD was used to generate both simulated high-speed angiograms and the corresponding ground-truth 3D flow fields to better understand the relationship between the 3D volumetric-flow distribution and the 2D projected-flow distribution as is obtained with angiography, and the subsequent 2D approximation of flow velocity. Several geometries were investigated, ranging from simple pipe models to complex patient-specific aneurysms. Simulated datasets were analyzed with the optical flow algorithm, and the effects of flow divergence, quantum mottle, and intensity gradient on the calculation were evaluated. From these simulations, we can evaluate whether flow fields reconstructed from HSA are representative of significant flow patterns in the 3D vasculature.

Initial investigations of scatter cross-talk in simultaneous biplane high-speed 1000 frames per second neuro-angiography using Monte Carlo simulations

The Aries photon counting detector (PCD) by Direct Conversion Inc. can image up to 1000 frames per second and is used to track contrast bolus in neuro-vasculature for hemodynamic calculations. For 3D tracking, synchronized biplane imaging with 1 ms acquisition times is used such that both imaging planes are exposed simultaneously. This leads to cross-scattered radiation being detected and a degradation of image quality compared to single-plane imaging. In this study, we utilize Monte Carlo (MC) methods to quantify the increase in scatter due to cross-talk without the use of a radiographic grid. EGSnrc biplane simulations were performed with the Zubal anthropomorphic head phantom. The total scatter plus primary and cross-scatter was calculated in the imaging planes for two orthogonal AP and lateral beams with a field size consistent with the 7.5×5 cm Aries detector, while the primary was determined with a 1×1 mm beam. The forward scatter was then determined from the difference between total and primary. The scatter is seen to increase by 4%-56% for AP projections and 48%-71% for lateral projections depending on detector orientation during simultaneous exposure. Scatter degradation from cross-talk can be reduced using an anti-scatter grid as well as the energy thresholding capabilities of the Aries PCD.

Comparison of quantitative imaging characteristics between a new, larger-FOV 1000 fps high-speed angiographic (HSA) photon-counting detector (Aries) with a smaller HSA detector (Actaeon)

High Speed Angiography (HSA) requires imaging detectors with both high-temporal and high-spatial resolution. Both the Aries and Acteon detectors by Direct Conversion (Stockholm, Sweden) are CdTe direct photon-counting detectors (PCD) that have acquisition frame rates of up to 1000-fps and a 100-micrometer pixel pitch; however, the new Aries detector offers a larger field of view (512 × 768 pixels) compared to the smaller Actaeon detector (256 × 256 pixels). An expanded field of view is required for imaging of larger vasculature, thus the Aries offers this advantage. Evaluations were performed of both detectors under Anti-Coincidence Circuitry (ACC-ON) mode, which corrects for charge sharing between pixels. Initial evaluations of instrumentation noise and detector energy-threshold calibration using Am-241 gamma spectroscopy were performed for the new Aries detector. Linearity was also evaluated for the Aries for each of the 12 individual modules that compose the detector field to check for homogeneity in response to exposure throughout the detector. Finally, Normalized Noise Power Spectrum (NNPS), Modulation Transfer Function (MTF) and Detective Quantum Efficiency (DQE) were then compared between the Aries and Actaeon detectors at two different exposures and detector energy thresholds. The detectors are linear up to approximately 1000 μR and have no instrumentation noise above a threshold of 15 keV. As expected, the MTF's and DQE's are similar between the Aries and Actaeon detectors, and there are thus no tradeoff's in image quality to achieve the larger FOV.

Real-time Detection of Patient Head Position and Cephalometric Landmarks from Neuro-Interventional Procedure Images Using Machine Learning for Patient Eye-Lens Dose Prediction

A deep learning (DL) model has been developed to estimate patient-lens dose in real-time for given exposure and geometric conditions during fluoroscopically-guided neuro-interventional procedures. Parameters input into the DL model for dose prediction include the patient head shift from isocenter and cephalometric landmark locations as a surrogate for head size. Machine learning (ML) models were investigated to automatically detect these parameters from the in-procedure fluoroscopic image. Fluoroscopic images of a Kyoto Kagaku anthropomorphic head phantom were taken at various known X (transverse) and Y (longitudinal) shifts, as well as different magnification modes, to create an image database. For each image, anatomical landmark coordinate locations were obtained manually using ImageJ and are used as ground-truth labels for training. This database was then used to train the two separate ML models. One ML model predicts the patient head shift in both the X and Y directions and the other model predicts the coordinates of the anatomical landmarks. From the coordinates, the distance between these anatomical landmarks is calculated, and input into the DL dose-prediction model. Model performance was evaluated using mean absolute error (MAE) and mean absolute percentage error (MAPE) for the head-shift and landmark-coordinate models, respectively. The goal is to implement these two separate models into the Dose Tracking System (DTS) developed by our group. This would allow the DTS to automatically detect the patient head size and position for eye-lens dose prediction and eliminate the need for manual input by the clinical staff.

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

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

Leveraging Patient-Specific Simulated Angiograms to Characterize Cerebral Aneurysm Hemodynamics using Computational Fluid Dynamics

Cerebral aneurysms (CA) affect nearly 6% of the US population and its rupture is one of the major causes of hemorrhagic stroke. Neurointerventionalists performing endovascular therapy (ET) to treat CA rely on qualitative image sequences obtained under fluoroscopy guidance alone, and do not have access to crucial quantitative information regarding blood flow before, during and after treatment - partially contributing to a failure rate of up to 30%. Computational fluid dynamics (CFD) is a powerful tool that can provide a wealth of quantitative data; however, CFD has found limited utility in the clinic due to the challenges in obtaining hemodynamic boundary conditions for each patient. In this work, we present a novel CFD-based simulated angiogram approach (SAA) that resolves the blood flow physics and interaction between blood and injected contrast agent to extract quantitative hemodynamic parameters which can be used to design real-time parametric imaging analysis. The SAA enables correlating contrast agent transport to the underlying hemodynamic conditions via time-density curves (TDC) obtained at several points in the region of interest. The ability of the TDC and the SAA to provide critical hemodynamic parameters in and around CA anatomies, such as washout and local flow changes is explored and presented. This provides invaluable quantitative data to the clinician at the time of intervention, since it incorporates the physics of blood flow and correlates the contrast transport to hemodynamic parameters quantitatively - thereby enabling the clinician to take informed decisions that improve treatment outcomes.

Estimation of patient skin dose in fluoroscopy: summary of a joint report by AAPM TG357 and EFOMP

BACKGROUND

Physicians use fixed C-arm fluoroscopy equipment with many interventional radiological and cardiological procedures. The associated effective dose to a patient is generally considered low risk, as the benefit-risk ratio is almost certainly highly favorable. However, X-ray-induced skin injuries may occur due to high absorbed patient skin doses from complex fluoroscopically guided interventions (FGI). Suitable action levels for patient-specific follow-up could improve the clinical practice. There is a need for a refined metric regarding follow-up of X-ray-induced patient injuries and the knowledge gap regarding skin dose-related patient information from fluoroscopy devices must be filled. The most useful metric to indicate a risk of erythema, epilation or greater skin injury that also includes actionable information is the peak skin dose, that is, the largest dose to a region of skin.

MATERIALS AND METHODS

The report is based on a comprehensive review of best practices and methods to estimate peak skin dose found in the scientific literature and situates the importance of the Digital Imaging and Communication in Medicine (DICOM) standard detailing pertinent information contained in the Radiation Dose Structured Report (RDSR) and DICOM image headers for FGI devices. Furthermore, the expertise of the task group members and consultants have been used to bridge and discuss different methods and associated available DICOM information for peak skin dose estimation.

RESULTS

The report contributes an extensive summary and discussion of the current state of the art in estimating peak skin dose with FGI procedures with regard to methodology and DICOM information. Improvements in skin dose estimation efforts with more refined DICOM information are suggested and discussed.

CONCLUSIONS

The endeavor of skin dose estimation is greatly aided by the continuing efforts of the scientific medical physics community, the numerous technology enhancements, the dose-controlling features provided by the FGI device manufacturers, and the emergence and greater availability of the DICOM RDSR. Refined and new dosimetry systems continue to evolve and form the infrastructure for further improvements in accuracy. Dose-related content and information systems capable of handling big data are emerging for patient dose monitoring and quality assurance tools for large-scale multihospital enterprises.

A parametric fitting technique for rapid determination of a skin-dose correction factor for angle of beam incidence during image-guided endovascular procedures

Skin dose is dependent on the incident beam angle and corrections are needed for accurate estimation of the risk of deterministic effects of the skin. To obtain the angular correction factors (ACF's), EGSnrc Monte Carlo (MC) software was used to calculate the skin dose as a function of incident x-ray beam angle at the center of the field for beam energies from 60 to 120 kVp, field sizes from 5 to 15 cm, and thicknesses of Cu beam filters from 0.2 to 0.5 mm. All MC simulations used 3×10¹⁰ incident photons. The dose was averaged over a 1 mm depth on the entrance surface of a 40×40 cm by 20 cm thick water phantom and was then normalized to the incident primary dose which was calculated using NIST mass energy absorption coefficients and by integrating over the beam energy spectrum. The Matlab tool, 'cftool', was used to fit these normalized dose values to power law equations as a function of incident beam angle, with coefficients that were fit to polynomials as a function of kVp. Separate fitting was done for different beam sizes and beam filters. The skin dose values calculated using the ACF determined from the fitted functional formulas agreed with that calculated by MC with a mean absolute percentage error (MAPE) less than 3% over the entire range of incident angles and kVp values. This fitting technique allows an ACF to be quickly determined for accurate skin dose calculation.

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.

Characterization of velocity patterns produced by pulsatile and constant flows using 1000 fps high-speed angiography (HSA)

In order to accurately quantify rapidly changing blood flow velocities, as typically seen in the neurovasculature, high temporal resolution is necessary. Current methods to extract velocity data from angiographic image sequences are generally limited to 30 fps or less. High-speed angiography (HSA) with a maximal frame rate of 1000 fps can be used to evaluate time-dependent flow details normally averaged out with lower frame rates. For new HSA image sequences, two different quantitative methods were utilized to extract high-temporal resolution velocity changes: X-Ray Particle Image Velocimetry (X-PIV) and optical flow (OF). A variety of flow conditions were examined in a range of patient-specific 3D-printed phantoms. Both pulsatile and constant flow settings were investigated. X-PIV was performed using radiopaque sub-millimeter microspheres, which were tracked throughout the image sequence to provide accurate, but limited sampling of the velocity field within the 3D-printed models. Also, an open source optical flow algorithm, OpenOpticalFlow, was used to perform velocity estimation based on the spatio-temporal intensity changes of iodinated contrast wavefronts. Periodic changes in velocity within each phantom ROI can be illustrated throughout the pulsatile cycle capture by the high-speed detector. In the constant flow sequences, changes in velocity across the phantom geometry can be seen. The ability to accurately measure detailed velocity distributions and velocity changes throughout various flow conditions at high temporal resolution enables further insight into the evaluation and treatment of neurovascular disease states.

The effect of underlying bone on the beam angular correction in calculating the skin dose of the head in neuro-interventional imaging

Skin dose depends on the surface shape, underlying tissue, beam energy, field size, and incident beam angle. These dependencies were determined in order to apply corrections in the skin-dose-tracking system (DTS) for accurate estimation of the risk of deterministic skin effects during fluoroscopically-guided neuro-interventional procedures. The primary-plus-scatter dose was calculated averaged over the skin thickness with underlying subcutaneous fat, and various thicknesses of skull bone on the surface of a cylindrical water phantom to simulate the head. The skin dose was calculated using EGSnrc Monte-Carlo (MC) software with 2×10¹⁰ incident photons and was normalized to the incident primary dose. Simulations were done for beam incident angles from 90 to 10 degrees with the skin surface, field sizes from 5 to 15 cm, bone thicknesses of 0, 1, 5, and 9 mm, and beam energies from 60 to 120 kVp. The results show the scatter-plus-primary to incident-primary dose ratio decreases with decreasing incident angle to the skin and with increasing thickness of underlying bone, while it increases with increasing field size and with increasing beam energy. The correction factor reduces the skin dose for angled rays and the reduction can be substantial for small angles of incidence, especially for angles below 50 degrees. For neuro-interventional procedures, the skin dose-area product (SDAP) with angular and bone correction is shown to be less than that without correction. The results of this study can be used to increase the accuracy of patient-skin-dose estimation for the head during fluoroscopic procedures.

Estimating Compton scatter distributions with a regressional neural network for use in a real-time staff dose management system for fluoroscopic procedures

Staff-dose management in fluoroscopic procedures is a continuing concern due to insufficient awareness of radiation dose levels. To maintain dose as low as reasonably achievable (ALARA), we have developed a software system capable of monitoring the procedure room scattered radiation and the dose to staff members in real-time during fluoroscopic procedures. The scattered-radiation display system (SDS) acquires imaging-system signal inputs to update technique and geometric parameters used to provide a color-coded mapping of room scatter. We have calculated a discrete look-up-table (LUT) of scatter distributions using Monte-Carlo (MC) software and developed an interpolation technique for the multiple parameters known to alter the spatial shape of the distribution. However, the file size for the LUT's can be large (~2GB), leading to long SDS installation times in the clinic. Instead, this work investigated the speed and accuracy of a regressional neural network (RNN) that we developed for predicting the scatter distribution from imaging-system inputs without the need for the LUT and interpolation. This method greatly reduces installation time while maintaining real-time performance. Results using error maps derived from the structural similarity index indicate high visual accuracy of predicted matrices when compared to the MC-calculated distributions. Dose error is also acceptable with a matrix element-averaged percent error of 31%. This dose-monitoring system for staff members can lead to improved radiation safety due to immediate visual feedback of high-dose regions in the room during the procedure as well as enhanced reporting of individual doses post-procedure.

Estimation of Patient Eye-Lens Dose During Neuro-Interventional Procedures using a Dense Neural Network (DNN)

The patient's eye-lens dose changes for each projection view during fluoroscopically-guided neuro-interventional procedures. Monte-Carlo (MC) simulation can be done to estimate lens dose but MC cannot be done in real-time to give feedback to the interventionalist. Deep learning (DL) models were investigated to estimate patient-lens dose for given exposure conditions to give real-time updates. MC simulations were done using a Zubal computational phantom to create a dataset of eye-lens dose values for training the DL models. Six geometric parameters (entrance-field size, LAO gantry angulation, patient x, y, z head position relative to the beam isocenter, and whether patient's right or left eye) were varied for the simulations. The dose for each combination of parameters was expressed as lens dose per entrance air kerma (mGy/Gy). Geometric parameter combinations associated with high-dose values were sampled more finely to generate more high-dose values for training purposes. Additionally, dose at intermediate parameter values was calculated by MC in order to validate the interpolation capabilities of DL. Data was split into training, validation and testing sets. Stacked models and median algorithms were implemented to create more robust models. Model performance was evaluated using mean absolute percentage error (MAPE). The goal for this DL model is that it be implemented into the Dose Tracking System (DTS) developed by our group. This would allow the DTS to infer the patient's eye-lens dose for real-time feedback and eliminate the need for a large database of pre-calculated values with interpolation capabilities.

Using a convolutional neural network for human recognition in a staff dose management software for fluoroscopic interventional procedures

Staff dose management is a continuing concern in fluoroscopically-guided interventional (FGI) procedures. Being unaware of radiation scatter levels can lead to unnecessarily high stochastic and deterministic risks due to the effects of absorbed dose by staff members. Our group has developed a scattered-radiation display system (SDS) capable of monitoring system parameters in real-time using a controller-area network (CAN) bus interface and displaying a color-coded mapping of the Compton-scatter distribution. This system additionally uses a time-of-flight depth sensing camera to track staff member positional information for dose rate updates. The current work capitalizes on our body tracking methodology to facilitate individualized dose recording via human recognition using 16-bit grayscale depth maps acquired using a Microsoft Kinect V2. Background features are removed from the images using a depth threshold technique and connected component analysis, which results in a body silhouette binary mask. The masks are then fed into a convolutional neural network (CNN) for identification of unique body shape features. The CNN was trained using 144 binary masks for each of four individuals (total of 576 images). Initial results indicate high-fidelity prediction (97.3% testing accuracy) from the CNN irrespective of obstructing objects (face masks and lead aprons). Body tracking is still maintained when protective attire is introduced, albeit with a slight increase in positional data error. Dose reports are then able to be produced which contain cumulative dose to each staff member at the eye lens level, waist level, and collar level. Individualized cumulative dose reporting through the use of a CNN in addition to real-time feedback in the clinic will lead to improved radiation dose management.

Comparison of skin dose calculated by the dose tracking system (DTS) with a beam angular correction factor and that calculated by Monte-Carlo

Skin dose is dependent on the incident beam angle and corrections are needed for accurate estimation of the risk of deterministic effects of the skin. Angular-correction factors (ACF) were calculated and incorporated into our skin-dose-tracking system (DTS) and the results compared to Monte-Carlo simulations for a neuro-interventional procedure. To obtain the ACF's, EGSnrc Monte-Carlo (MC) software was used to calculate the dose averaged over 0.5, 1, 2, 3, 4 and 5 mm depth into the entrance surface of a water phantom at the center of the field as a function of incident beam to skin angle from 90-10 degrees for beam field sizes from 5-15 cm and for beam energies from 60-120 kVp. These values were normalized to the incident primary dose to obtain the ACF. The angle of incidence at each mesh vertex in the beam on the surface of the DTS patient graphic was calculated as the complement of the angle between the normal vector and the vector of the intersecting ray from the tube focal spot; skin dose at that vertex was calculated using the corresponding ACF. The skin-dose values with angular correction were compared to those calculated using MC with a matching voxelized phantom. The results show the ACF decreases with decreasing incident angle and skin thickness, and increases with increasing field size and kVp. Good agreement was obtained between the skin dose calculated by the angular-corrected DTS and MC, while use of the ACF allows the real-time performance of the DTS to be maintained.

Assessment of Eye Lens Dose Reduction When Using Lateral Lead Shields on the Patient's Head during Neurointerventional Fluoroscopic Procedures and Cone-beam Computed Tomography (CBCT) Scans

The purpose of this study was to evaluate the effect of placing small lead shields on the temple region of the skull to reduce radiation dose to the lens of the eye during interventional fluoroscopically-guided procedures and cone-beam computed tomography (CBCT) scans of the head. EGSnrc Monte-Carlo code was used to determine the eye lens dose reduction when using lateral lead shields for single x-ray projections, CBCT scans with different protocols, and interventional neuroradiology procedures with the Zubal computational head phantom. A clinical C-Arm system was used to take radiographic projections and CBCT scans of anthropomorphic head phantoms without and with lead patches, and the images were compared to assess the effect of the shields. For single lateral projections, a 0.1 (0.3)-mm-thick lead patch reduced the dose to the left-eye lens by 40% to 60% (55% to 80%) from 45° to 90° RAO and to the right-eye lens by around 30% (55%) from 70° to 90° RAO. For different CBCT protocols, the reduction of lens dose with a 0.3-mm-thick lead patch ranged from 20% to 53% at 110 kVp. For CBCT scans of the anthropomorphic phantom, the lead patch introduced streak artifacts that were mainly in the orbital regions but were insignificant in the brain region where most neurointerventional activity occurs. The dose to the patient's eye lens can be reduced considerably by placing small lead shields over the temple region of the head without substantially compromising image quality in neuro-imaging procedures.

Single-center experience of using high definition (Hi-Def) imaging during neurointervention treatment of intracranial aneurysms using flow diverters

BACKGROUND

A new dual resolution imaging x-ray detector system (Canon Medical Systems Corporation, Tochigi, Japan) has a standard resolution 194 µm pixel conventional flat-panel detector (FPD) mode and a high-resolution 76 µm high-definition (Hi-Def) mode in a single unit. The Hi-Def mode enhances the visualization of the intravascular devices.

OBJECTIVE

We report the clinical experience and physician evaluation of this new detector system with Hi-Def mode for the treatment of intracranial aneurysms using a Pipeline embolization device (PED).

METHODS

During intervention at our institute, under large field of view (FOV) regular resolution FPD mode imaging, the catheter systems and devices were first guided to the proximity of the treatment area. Final placement and deployment of the PED was performed under Hi-Def mode guidance. A post-procedure 9-question physician survey was conducted to qualitatively assess the impact of Hi-Def mode visualization on physicians' intraoperative decision-making. One-sample t-test was performed on the responses from the survey. Dose values reported by the x-ray unit were also recorded.

RESULTS

Twenty-five cases were included in our study. The survey results indicated that, for each of the nine questions, the physicians in all cases indicated that the Hi-Def mode improved visualization compared with the FPD mode. For the 25 cases, the mean cumulative entrance air kerma was 2.35 Gy, the mean dose area product (DAP) was 173.71 Gy.cm², and the mean x-ray exposure time was 39.30 min.

CONCLUSIONS

The Hi-Def mode improves visualization of flow diverters and may help in achieving more accurate placement and deployment of devices.

High-Definition Zoom Mode: A High Resolution X-ray Microscope for Neurointerventional Treatment Procedures

BACKGROUND AND PURPOSE

Visualization of structural details of treatment devices during neurointerventional procedures can be challenging. A new true two-resolution imaging X-ray detector system features a 194 µm pixel conventional flat-panel detector (FPD) mode and a 76 µm pixel high-resolution high-definition (Hi-Def) zoom mode in one detector panel. The Hi-Def zoom mode was developed for use in interventional procedures requiring superior image quality over a small field of view (FOV). We report successful use of this imaging system during intracranial aneurysm treatment in 1 patient with a Pipeline-embolization device and 1 patient with a low-profile visualized intramural support (LVIS Blue) device plus adjunctive coiling.

METHODS

A guide catheter was advanced from the femoral artery insertion site to the proximity of each lesion using standard FPD mode. Under magnified small FOV Hi-Def imaging mode, an intermediate catheter and microcatheters were guided to the treatment site, and the PED and LVIS Blue plus coils were deployed. Radiation doses were tracked intraprocedurally.

RESULTS

Critical details, including structural changes in the PED and LVIS Blue and position and movement of the microcatheter tip within the coil mass, were more readily apparent in Hi-Def mode. Skin-dose mapping indicated that Hi-Def mode limited radiation exposure to the smaller FOV of the treatment area.

CONCLUSIONS

Visualization of device structures was much improved in the high-resolution Hi-Def mode, leading to easier, more controlled deployment of stents and coils than conventional FPD mode.

Flow-Pattern Details in an Aneurysm Model Using High-Speed 1000-Frames-per-Second Angiography

Traditional digital subtraction angiography provides rather limited evaluation of contrast flow dynamics when studying and treating intracranial brain aneurysms. A 1000-frames-per-second photon-counting x-ray detector was used to image detailed iodine-contrast flow patterns in an internal carotid artery aneurysm of a 3D-printed vascular phantom. High-speed imaging revealed differences in vortex and inflow patterns with and without a Pipeline Embolization Device flow diverter in more detail and clarity than could be seen in standard pulsed angiography. Improved temporal imaging has the potential to impact the outcomes of endovascular interventions by allowing clinicians to better understand and act on flow dynamics in real-time.

High-Definition Zoom Mode, a High-Resolution X-Ray Microscope for Neurointerventional Treatment Procedures: A Blinded-Rater Clinical-Utility Study

BACKGROUND AND PURPOSE

Quality of visualization of treatment devices during critical stages of endovascular interventions, can directly impact their safety and efficacy. Our aim was to compare the visualization of neurointerventional procedures and treatment devices using a 194-μm pixel flat panel detector mode and a 76-μm pixel complementary metal oxide semiconductor detector mode (high definition) of a new-generation x-ray detector system using a blinded-rater study.

MATERIALS AND METHODS

Deployment of flow-diversion devices for the treatment of internal carotid artery aneurysms was performed under flat panel detector and high-definition-mode image guidance in a neurointerventional phantom simulating patient cranium and tissue attenuation, embedded with 3D-printed intracranial vascular models, each with an aneurysm in the ICA segment. Image-sequence pairs of device deployments for each detector mode, under similar exposure and FOV conditions, were evaluated by 2 blinded experienced neurointerventionalists who independently selected their preferred image on the basis of visualization of anatomic features, image noise, and treatment device. They rated their selection as either similar, better, much better, or substantially better than the other choice. Inter- and intrarater agreement was calculated and categorized as poor, moderate, and good.

RESULTS

Both raters demonstrating good inter- and intrarater agreement selected high-definition-mode images with a frequency of at least 95% each and, on average, rated the high-definition images as much better than flat panel detector images with a frequency of 73% from a total of 60 image pairs.

CONCLUSIONS

Due to their higher resolution, high-definition-mode images are sharper and visually preferred compared with the flat panel detector images. The improved imaging provided by the high-definition mode can potentially provide an advantage during neurointerventional procedures.

Evaluation of Methods of Displaying the Real-Time Scattered Radiation Distribution during Fluoroscopically-Guided Interventions for Staff Dose Reduction

2D and 3D scatter dose display options are evaluated for usefulness and ease of interpretation for real-time feedback to staff to facilitate changes in individual positioning for dose reduction purposes, as well as improving staff consciousness of radiation presence. Room-sized scatter dose 3D matrices are obtained utilizing Monte Carlo simulations in EGSnrc. These distributions are superimposed on either a ceiling-view 2D graphic of the patient and table for reference or a 3D augmented reality (AR) display featuring a real-time video feed of the interventional room. A slice of the scatter dose matrix, at a selectable distance above the floor, is color-coded and superimposed on the computer graphic or AR display. The 3D display obtains depth information from a ceiling mounted Microsoft Kinect camera, which is equipped with a 1080p visual camera, as well as an IR emitter/receiver to generate a depth map of the interventional suite and persons within it. The 3D depth information allows parts of objects above the 2D dose map to pass through the map without being colorized by it so the height perspective of the dose map can be maintained. The 2D and 3D displays incorporate network information from the imaging system to scale the scatter dose with exposure factors and adjust rotation of the distribution to match the gantry. Demonstration images were displayed to neurosurgery interventional staff and survey responses were collected. Results from the survey indicated that scatter distribution displays would be desirable and helpful in managing staff dose.

Calculation of Forward Scatter Dose Distribution at the skin entrance from the patient table for fluoroscopically guided interventions using a pencil beam convolution kernel

The forward-scatter dose distribution generated by the patient table during fluoroscopic interventions and its contribution to the skin dose is studied. The forward-scatter dose distribution to skin generated by a water table-equivalent phantom and the patient table are calculated using EGSnrc Monte-Carlo and Gafchromic film as a function of x-ray field size and beam penetrability. Forward scatter point spread function's (PSFn) were generated with EGSnrc from a 1×1 mm simulated primary pencil beam incident on the water model and patient table. The forward-scatter point spread function normalized to the primary is convolved over the primary-dose distribution to generate scatter-dose distributions. The utility of PSFn to calculate the entrance skin dose distribution using DTS (dose tracking system) software is investigated. The forward-scatter distribution calculations were performed for 2.32 mm, 3.10 mm, 3.84 mm and 4.24 mm Al HVL x-ray beams for 5×5 cm, 9×9 cm, 13.5×13.5 cm sized x-ray fields for water and 3.1 mm Al HVL x-ray beam for 16.5×16.5 cm field for the patient table. The skin dose is determined with DTS by convolution of the scatter dose PSFn's and with Gafchromic film under PMMA "patient-simulating" blocks for uniform and for shaped x-ray fields. The normalized forward-scatter distribution determined using the convolution method for water table-equivalent phantom agreed with that calculated for the full field using EGSnrc within ±6%. The normalized forward-scatter dose distribution calculated for the patient table for a 16.5×16.5 cm FOV, agreed with that determined using film within ±2.4%. For the homogenous PMMA phantom, the skin dose using DTS was calculated within ±2 % of that measured with the film for both uniform and non-uniform x-ray fields. The convolution method provides improved accuracy over using a single forward-scatter value over the entire field and is a faster alternative to performing full-field Monte-Carlo calculations.

Initial investigations of a special high-definition (Hi-Def) zoom capability in a new detector system for neuro-interventional procedures

Real-time visualization of fine details ranging to 100 um or less in neuro-vascular imaging guided interventions is important. A separate high-resolution detector mounted on a standard flat panel detector (FPD) was previously reported. This device had to be rotated mechanically into position over the FPD for high resolution imaging. Now, the new detector reported here has a high definition (Hi-Def) zoom capability along with the FPD built into one unified housing. The new detector enables rapid switching, by the operator between Hi-Def and FPD modes. Standard physical metrics comparing the new Hi-Def modes with those of the FPD are reported, demonstrating improved imaging resolution and noise capability at patient doses similar to those used for the FPD. Semi-quantitative subjective studies involving qualitative clinician feedback on images of interventional devices such as a Pipeline Embolization Device (PED) acquired in both Hi-Def and FPD modes are presented. The PED is deployed in a patient specific 3D printed neuro-vascular phantom embedded inside realistic bone and with tissue attenuating material. Field-of-view (FOV), exposure and magnification were kept constant for FPD and Hi-Def modes. Static image comparisons of the same view of the PED within the phantom were rated by expert interventionalists who chose from the following ratings: Similar, Better, or Superior. Generally, the Hi-Def zoomed images were much preferred over the FPD, indicating the potential to improve endovascular procedures and hence outcomes using such a Hi-Def feature.

Spectroscopy with a CdTe-based photon-counting imaging detector (PCD) having charge sharing correction capability

The spectroscopic capabilities of a newly upgraded version of a prototype imaging photon counting detector (PCD) was investigated. The XCounter Actaeon has four acquisition modes in which signal processing is varied including one mode having a charge sharing correction so that neighboring pixels that share a detected event will not be erroneously counted twice, hence it is designated the Anti-Coincidence Circuit On or ACC On mode. Since this CdTe-based direct conversion PCD has 100 μm pixels, such charge sharing may frequently occur for typical medical x-ray energies. Each pixel of this PCD has two scalers and two energy discriminators that enable counting without instrumentation noise of events above each threshold level; hence, a spectrum can be obtained by sequentially moving the thresholds of both discriminators. It became evident from the spectra for the various acquisition modes that only those obtained with the charge sharing correction enabled, compared favorably with theoretically predicted spectra. After verifying the energy calibration using the mono-energetic emissions from an Am-241 source, spectra at various kVps from a standard medical x-ray generator were obtained. The spectra generated by ACC On mode for 70 kVp and 110 kVp were the closest match to the theoretical spectra generated by SpekCal. For dual energy applications, ACC On mode with charge sharing correction circuitry would be the best choice among various acquisition modes. Also investigated was the dual energy imaging capability of the Actaeon PCD with ACC On mode to separate Aluminum and Iodine while imaging an artery stenosis phantom.

A simulation platform using 3D printed neurovascular phantoms for clinical utility evaluation of new imaging technologies

Modern 3D printing technology allows rapid prototyping of vascular phantoms based on an actual human patient with a high degree of precision. Using this technology, we present a platform to accurately simulate clinical views of neuro-endovascular interventions and devices. The neuro-endovascular interventional phantom has a 3D printed cerebrovasculature model derived from a patient CT angiogram and embedded inside a human skull providing bone attenuation. Acrylic layers were placed underneath and on top of the skull, simulating entrance and exit tissue attenuation and also simulating forward scatter. The 3D model was connected to a pulsatile flow loop for simulating interventions using clinical devices such as catheters and stents. To validate the x-ray attenuation and establish clinical accuracy, the automatic exposure selection by a clinical c-arm system for the phantom was compared with that for a commercial anthropomorphic head phantom (SK-150, Phantom Labs). The percentage difference between automatic exposure selection for the neuro-intervention phantom and the SK-150 phantom was under 10%. By changing 3D printed models, various patient diseased anatomies can be simulated accurately with the necessary x-ray attenuation. Using this platform various interventional procedures were performed using new imaging technologies such as a high-resolution x-ray fluoroscope and a dose-reduced region-of-interest attenuator and differential temporally filtered display for enhanced interventional imaging. Simulated clinical views from such phantom-based procedures were used to evaluate the potential clinical performance of such new technologies.

Developing a database of 3-D scattered radiation distributions for a c-arm fluoroscope as a function of exposure parameters and phantom

The purpose of this work is to develop a database of 3D scattered radiation dose-rate distributions to estimate the staff dose by location around a C-Arm fluoroscopic system in an interventional procedure room. The primary x-ray beam of a Toshiba Infinix fluoroscopy machine was modeled using EGSnrc Monte Carlo code and the scattered radiation distributions were calculated using 5 × 10⁹ photons per simulation. These 3D distributions were determined over the volume of the room as a function of various parameters such as the beam kVp and beam filter, the size and shape of the field, the angulation of the C-arm, and the phantom size and shape. Two phantom shapes were used in this study: cylindrical and super-ellipses. The results show that shape of the phantom will affect the dose-rate distribution at distances less than 100 cm, with a higher intensity for the super-ellipse. The scatter intensity per entrance air kerma is seen to be approximately proportional to field area and to increase with increasing kVp. The scatter changes proportionally with increases in primary entrance air kerma for factors such as pulse rate, mA and pulse width. This database will allow estimation of the scatter distribution in the procedure room and, when displayed to the staff during a procedure, may facilitate a reduction of occupational dose.

A Patient Dose-Reduction Technique for Neuroendovascular Image-Guided Interventions: Image-Quality Comparison Study

BACKGROUND AND PURPOSE

The ROI-dose-reduced intervention technique represents an extension of ROI fluoroscopy combining x-ray entrance skin dose reduction with spatially different recursive temporal filtering to reduce excessive image noise in the dose-reduced periphery in real-time. The aim of our study was to compare the image quality of simulated neurointerventions with regular and reduced radiation doses using a standard flat panel detector system.

MATERIALS AND METHODS

Ten 3D-printed intracranial aneurysm models were generated on the basis of a single patient vasculature derived from intracranial DSA and CTA. The incident dose to each model was reduced using a 0.7-mm-thick copper attenuator with a circular ROI hole (10-mm diameter) in the middle mounted inside the Infinix C-arm. Each model was treated twice with a primary coiling intervention using ROI-dose-reduced intervention and regular-dose intervention protocols. Eighty images acquired at various intervention stages were shown twice to 2 neurointerventionalists who independently scored imaging qualities (visibility of aneurysm-parent vessel morphology, associated vessels, and/or devices used). Dose-reduction measurements were performed using an ionization chamber.

RESULTS

A total integral dose reduction of 62% per frame was achieved. The mean scores for regular-dose intervention and ROI dose-reduced intervention images did not differ significantly, suggesting similar image quality. Overall intrarater agreement for all scored criteria was substantial (Kendall τ = 0.62887; P < .001). Overall interrater agreement for all criteria was fair (κ = 0.2816; 95% CI, 0.2060-0.3571).

CONCLUSIONS

Substantial dose reduction (62%) with a live peripheral image was achieved without compromising feature visibility during neuroendovascular interventions.

Gantry rotational motion-induced blur in cone-beam computed tomography

As neuro-endovascular image-guided interventions (EIGIs) make use of higher resolution detectors, gantry rotational motion-induced blur becomes more noticeable in acquired projections as well as reconstructed images by reducing the visibility of vascular and device features whose visualization could be critical in the treatment of vascular pathology. Motion-induced blur in projections views is a function of an object's position in the field-of-view (FOV), gantry rotational speed, and frame capture or exposure time. In this work different frame rates were used to investigate the effects of blurring from individual projections on the reconstructed image. To test the effects of these parameters on reconstructed images, a regular pattern phantom of small objects was simulated and projection views were generated at various different frame rates for a given simulated rotational velocity. The reconstruction was made using a linear interpolation of filtered backprojections. Images reconstructed from lower frame rates showed significant blurring in the azimuthal direction, increasingly worse towards the periphery of the image. However, those reconstructed from higher frame rates showed significantly less blur throughout the entire FOV. While lower frame rates could be used with slower gantry speeds this would increase the risk of voluntary or involuntary patient motion contributing to blur over the entire FOV. A high frame rate used with high gantry speeds could reliable provide images without gantry-motion blur while reducing the risk of patient-motion blur. Frame rates exceeding 2000 fps available with photon counting detectors such as the X-counter Actaeon¹ are available.

Anti-scatter grid artifact elimination for high-resolution x-ray imaging detectors without a prior scatter distribution profile

Using anti-scatter grids with high-resolution imaging detectors could result in grid-line artifacts, with increasing severity as detector resolution improves. Grid-line mask subtraction can result in residual artifacts that are due to residual scatter penetrating the grid and not subtracted. By subtracting this residual scatter, the grid artifacts can be minimized. In the previous works, an initial residual-scatter estimate was derived by placing lead markers on a test object; however, any change in the object geometry requires a new scatter estimate. Such a method is impractical to implement during a clinical procedure. In this work, we present a new method to derive the initial scatter estimate to eliminate grid-line artifacts during a procedure. A standard stationary Smit-Roentgen x-ray grid (line density - 70 lines/cm, grid ratio - 13:1) was used with a high-resolution CMOS detector (Dexela Model 1207, pixel size - 75 μm) to image an anthropomorphic head phantom. The initial scatter estimate was derived from the image itself and the grid artifacts were eliminated using recursive correction estimation; this result was compared to that with the estimate derived from placing lead markers on the phantom. In both cases, the contrast-to-noise ratio (CNR) was improved compared to the original image with grid artifacts. Percentage differences in CNR's for three regions between the images corrected with the two estimates were less than 5%. With the new method no a priori scatter distribution profiles are needed, eliminating the need to have libraries of pre-calculated scatter profiles and making the implementation more clinically practical.

Investigation of organ dose variation with adult head size and pediatric age for neuro-interventional projections

The purpose of this study was to evaluate the effect of patient head size on radiation dose to radiosensitive organs, such as the eye lens, brain and spinal cord in fluoroscopically guided neuro-interventional procedures and CBCT scans of the head. The Toshiba Infinix C-Arm System was modeled in BEAMnrc/EGSnrc Monte-Carlo code and patient organ and effective doses were calculated in DOSxynrc/EGSnrc for CBCT and interventional procedures. X-ray projections from different angles, CBCT scans, and neuro-interventional procedures were simulated on a computational head phantom for the range of head sizes in the adult population and for different pediatric ages. The difference of left-eye lens dose between the mean head size and the mean ± 1 standard deviation (SD) ranges from 20% to 300% for projection angles of 0° to 90° RAO. The differences for other organs do not vary as much and is only about 10% for the brain. For a LCI-High CBCT protocol, the difference between mean and mean ± 1 SD head size is about 100% for lens dose and only 10% for mean and peak brain dose; the difference between 20 and 3 year-old mean head size is an increase of about 200% for the eye lens dose and only 30% for mean and peak brain dose. Dose for all organs increases with decreasing head size for the same reference point air kerma. These results will allow size-specific dose estimates to be made using software such as our dose tracking system (DTS).

Use of material decomposition in the context of neurovascular intervention using standard flat panel and a high-resolution CMOS detector

The imaging of endovascular devices during neurovascular procedures such as the coiling of aneurysms guided with CBCT imaging may be challenging due to the presence of highly attenuating materials such as platinum in the coil and stent marker, nickel-titanium in the stent, iodine in the contrast agent, and tantalum in the embolization agent. The use of dual-energy imaging followed by a basis material decomposition image processing-scheme may improve the feature separation and recognition. Two sets of testing were performed to validate this concept. The first trial was the acquisition of dual-energy micro-CBCT data of a 3D-printed simple aneurysm model using a 49.5 μm pixel size CMOS detector (Teledyne DALSA, Waterloo, ON.). Two sets of projection data were acquired using beam energies of 35 kVp and 70 kVp. Axial slices were reconstructed and used to carry out the material decomposition processing. The second trial was the acquisition of dual-energy CBCT images of a RS-240T angiographic head phantom (Radiology Support Devices Inc., CA.) with an iodine vascular insert using a Toshiba Infinix BiPlane C-arm system coupled to a flat panel detector. Two sets of image data were acquired using beam energies of 80 kVp and 120 kVp. Following image reconstruction, slices of the phantom were decomposed using the same processing as previously. The resulting image data over both trials indicate that the decomposition process was successful in separating the kinds of materials commonly used during a neurovascular intervention, such as platinum, cobalt-chromium, and iodine. The normalized root mean square error metric was used to quantitatively assess this. This indicates a basis for future more clinically relevant testing of our methods.

Evaluation of a new photon-counting imaging detector (PCD) with various acquisition modes

The prospect of improved low noise, high speed, and dual-energy imaging that may be associated with the use of photon-counting imaging detectors (PCD) has motivated this evaluation of a newly upgraded version of a prototype PCD. The XCounter Actaeon was evaluated in its four acquisition modes each based upon varying signal processing firmware including a mode with charge sharing correction that enables neighboring pixels that share the energy from one incident x-ray photon detection to be counted only once at the proper summed energy in the pixel with the largest charge deposition. Since this PCD is a CdTe-based direct detector with 100 μm pixels, such charge sharing for typical medical x-ray energy photons may occur frequently and must be corrected to achieve more accurate counts. This charge sharing correction is achieved with an Anti-Coincidence Circuit (ACC) which prevents double pixel counting from one event as well as prevents counting from either event if they are below a preset threshold. Various physical parameters of the PCD were evaluated including linearity, sensitivity, pulse pile-up effects, dark noise, spatial resolution, noise power spectrum, and detective quantum efficiency.

Assessment of organ and effective dose when using region-of-interest attenuators in cone-beam CT and interventional fluoroscopy

In some medical-imaging procedures using cone-beam CT (CBCT) and fluoroscopy, only the center of the field of view (FOV) may be needed to be visualized with optimal image quality. To reduce the dose to the patient while maintaining visualization of the entire FOV, a Cu attenuator with a circular aperture for the region of interest (ROI) is used. The potential organ and effective dose reductions of ROI imaging when applied to CBCT and interventional fluoroscopic procedures were determined using EGSnrc Monte Carlo code. The Monte Carlo model was first validated by comparing the surface dose distribution in a solid-water block phantom with measurement by Gafchromic film. The dependence of dose reduction on the ROI attenuator thickness, the opening size of the ROI, the axial beam position, and the location of the different organs for both neuro and thoracic imaging was evaluated. The results showed a reduction in most organ doses of 45% to 70% and in effective dose of 46% to 66% compared to the dose in a CBCT scan and in an interventional procedure without the ROI attenuator. This work provides evidence of a substantial reduction of organ and effective doses when using an ROI attenuator during CBCT and fluoroscopic procedures.

Calculation of the entrance skin dose distribution for fluoroscopically guided interventions using a pencil beam backscatter model

Radiation backscattered from the patient can contribute substantially to skin dose in fluoroscopically guided interventions (FGIs). The distribution of backscatter is not spatially uniform, and use of a single backscatter factor cannot provide an accurate determination of skin dose. This study evaluates a method to determine the backscatter spatial distribution through convolution of a backscatter-to-primary (BP) point spread function (PSFn). The PSFn is derived for a pencil beam using EGSnrc Monte Carlo software and is convolved with primary distributions using a dose-tracking system. The backscatter distribution calculated using the convolution method is validated with Monte Carlo-derived distributions for three different size "uniform" fields and with XR-QA2 Gafchromic film for nonuniform x-ray fields obtained using region-of-interest (ROI) attenuators and compensation filters, both with homogenous poly-methyl methacrylate and nonhomogenous head phantoms. The BP ratios inside uniform fields were calculated within [Formula: see text] of that determined using EGSnrc. For shaped fields, the BP ratio in the unattenuated ROI was calculated within [Formula: see text] of that measured with film; in the beam-attenuated periphery, agreement was within [Formula: see text], due to the larger uncertainty of the dose-response curve of the film in the low-dose region. This backscatter PSFn convolution method is much faster than performing full-field Monte Carlo calculations and provides improved accuracy in skin dose distribution determination for FGI procedures.