Search

OBJECTIVES: 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. METHODS: 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. RESULTS: 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 neuro-interventional activity occurs. CONCLUSION: 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.

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.

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.

The imaging parameters used in neurointerventional procedures were evaluated to better understand the exposure techniques used clinically and their impact on patient dose. All parameters are available on the imaging system’s network bus in real time for each exposure pulse during a procedure. The Canon Dose Tracking System (DTS), which we developed, records the parameters of each exposure event in a raw data log file of controller area network (CAN) packets. We have collected such log files for 120 neurointerventional cases. Parameters are extracted by converting the raw data log file to a reformatted text file using a MATLAB script. The text file is input into a Microsoft Visual Studio project which outputs a new text file, which, with a reference table, allows the parameters to be identified. A written Python script is used to extract the specific parameters that were to be evaluated and output a .csv file. These were then input into MATLAB for analysis. The parameters extracted were the kVp, beam filter type, mAs, and the cranial/caudal angle as well as the RAO/LAO angle for the frontal and lateral gantries for DA and pulsed fluoroscopy (PF) modes. The gantry angles ranged from 34⁰ CRA to 42⁰ CAU and from 114⁰ RAO to 91⁰ LAO for DA and PF, respectively. The median kVp was 84 and 73 and the average per frame mAs was about 11 and 1.8 for DA and PF, respectively. This analysis should allow a better understanding of clinical practice in order to relate technique to patient and staff dose.

Angiographic parametric imaging (API) is a quantitative imaging method which uses digital subtraction angiography (DSA) to calculate biomarkers related to hemodynamics. This method has been used for neurovascular disease diagnosis and treatment outcome predictions in clinical settings but results are regarded with caution since derived biomarkers are strongly correlated with contrast injection parameters. This study aimed to assess utilization of truncated singular value decomposition (TSVD) in correcting API maps across various injection rates. Digital angiography data was collected using two neurovascular phantoms embedded in a simulated flow loop. Contrast volumes of 5 and 10 mL along with injection rates of 5, 10, 15, and 20 mL/sec were utilized during testing. API maps were generated with baseline and stenosis models using gamma variate fitting along with TSVD of the arterial input of the phantom. Surrogate regional blood flow (sRBF) and regional blood volume (sRBV) maps indicate consistent values across varying injection rates along with decreases in flow and volumes following introduction of a stenosis (Baseline: sRBF=53.3±4.20 arbitrary volume units (AVU)/min, sRBV=2.66±0.14 AVU, Stenosis: sRBF=28.6±3.78 AVU/min, sRBV=1.75±0.45 AVU). Mean transit time (MTT) and time of maximal residue function (Tmax) maps indicate consistent and decreasing parameter values respectively as injection rates increase along with increases in each parameter in the presence of a stenosis (Baseline: MTT=0.72±0.14 sec, Tmax=1.36±0.14 sec, Stenosis: MTT=0.94±0.13 sec, Tmax=1.71±0.18 sec). This study indicates TSVD has the potential to normalize API parameter maps across various injection rates potentially allowing for the implementation of API in ischemic stroke diagnosis and treatment.

It is important to determine the patient’s skin dose accurately for fluoroscopic interventional procedures in order to estimate the risk of deterministic injury. The purpose of this study is to investigate how the patient’s skin dose changes as a function of x-ray beam incident angle for flat and curved surfaces. The primary and scatter dose was calculated averaged over a 2.0 mm depth at the surfaces of both cubic and cylindrical water phantoms to simulate different patient curvature. The total skin dose was calculated using EGSnrc Monte-Carlo (MC) software with 1010 photons incident and the primary dose was calculated at the central axis using the mass energy absorption coefficients published by NIST and integrated over the beam-energy spectrum. Simulations were done for incident angles from 90 to 10 degrees, beam field sizes from 5 to 15 cm, cylinder diameters from 20 to 30 cm, 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 due to increased primary attenuation and decreases from cubic to cylindrical phantom and with decreasing cylinder diameter at all angles due to reduced backscatter. These results can be used to determine angular correction factors needed to accurately estimate patient skin dose when the beam is not normal to the entrance surface during fluoroscopic procedures.

Purpose: 3D printed Patient-Specific Neurovascular Phantoms (3DP-PSNP) containing significant portions of the neurovasculature can be used to develop and test new diagnostic tools. The purpose of this research is to assess the use of 3DP-PSNP to study the correlation between Angiographic Parametric Imaging (API) features and the severity of carotid artery disease. Materials and Methods: We developed 3DP-PSNP for twenty patients with carotid atherosclerosis and performed two studies. In the first study, we used three phantoms with complete Circle of Willis (COW) with none, moderate and severe stenosis respectively. In the second experiment, all phantoms regardless of the COW structure were used. 3DP-PSNPs were connected in a simulated physiological pulsatile flow loop and Digital Subtractive Angiography (DSA) was performed by injecting 10 ml of contrast at 10 mL/s. An API software calculated imaging biomarkers: time-to-peak TTP, mean transit time MTT, time to arrival TTA, peak height PH, and area under the curve AUC for both carotids. Results: For none, moderate, and severe stenosis respectively, absolute mean percent differences between diseased vs. contralateral carotids were: 1.9%, 8.1%, 32.3% TTP; 6.8%, 7.5%, 41.0% TTA; 13.2%, 12.4%, 6.4% MTT; 8.7%, 11.3%, 97.6% PH; and 10.3%, 22.6%, 100% AUC. Injection parameters did not cause significant change to the difference between diseased and contralateral carotids. The second experiment showed strong correlation between the TTA and location of the stenosis, regardless of COW configuration. Conclusions: Overall, API was assessed in 3DP-PSNPs and shown to have increasing differences between diseased and contralateral carotid arteries with increasing severity.

Purpose: 3D-printing of patient-specific phantoms such as the mitral valve (MV) is challenging due to inability of current imaging systems to reconstruct fine moving features and 3D printing constraints. We investigated methods to 3D-print MV structures using ex-vivo micro-CT. Materials and Methods: A dissected porcine MV was imaged using micro-CT in diastole, using a special fixation holder. The holder design was based on a patient ECG gated cardiac CT scan using as reference points the papillary muscles and annulus. Next the micro-CT volume was segmented and 3D-printed in various elastic materials. We tested different postprocessing techniques for support material removal and surface coatings to preserve the MV integrity. To test the error a Cloud Comparison of the porcine valve-mesh file and the valve-mesh file from the patient ECG gated cardiac CT scan was performed. Results: Best results for the 3D-printed models were achieved using TangoPlus poly-jet material with a Objet Eden printer. The error computation yielded a 2.6mm deviation-distance between the two aligned valves indicating adequate alignment. The post-processing methods for support removal were challenging and required 24+ hours sample-emersion in slow agitating sodium hydroxide baths. Conclusions: The most challenging part for MV manufacturing is 3D volume acquisition and the post-printing methods during support cleaning. We developed methods to circumvent both, the imaging and the 3D-printing challenges and to ensure that the final phantom includes the fine chordae and valve geometry. Using these solutions, we were able to create complete MV structures which could benefit medical research and device testing.

Vascular pathologies such as stenosis cause changes in blood flow patterns and rates in patient vasculature. These changes are difficult to observe using conventional x-ray imaging techniques that typically operate at no more than 30 fps. However, 1000 fps High-Speed Angiography (HSA) provides sufficient temporal resolution to observe flow patterns in patient vasculature. The Actaeon detector from XCounter provides 1000 fps HSA with high spatial resolution (100 μm pixel pitch) enabling observation of detailed blood flow in 3D printed phantoms. To observe the effects of stenosis on flow, three carotid artery phantoms beginning 7.6 cm proximal to the bifurcation of the common carotid and extending 6.9 cm into the internal and external carotids with different degrees of stenosis (0%, 33%, and 66% stenosis) in the internal carotid were printed. Iodine contrast was injected into the phantom proximal to the bifurcation and images were recorded at 1000 fps. Time-density curves were generated for ROIs distal and proximal to the location of the stenosis for both branches of the carotid artery. The curves for the control phantom (no stenosis) showed little difference between the branches of the carotid while the stenotic phantoms showed an increased velocity in the flow of contrast and a lower maximum signal intensity for the stenosed internal carotid. 1000 fps HSA could be used during treatment of vascular pathologies to observe changes in flow patterns following an endovascular treatment.

While angiography may be considered the gold standard for evaluating diseases of the human vasculature, vascular flow details are unavailable due to the low temporal resolution of flat panel detectors (FPDs) which operate at a maximum of 15 – 30 fps. Higher frame rates are necessary to extract any meaningful flow detail, which may act as additional information that can be used to characterize flow-dependent disease states. These higher rates have become available with recent advances in photon-counting detector (PCD) technology. The XCounter Actaeon was used to perform high frame rate imaging at 1000 fps. The Actaeon also provides superior spatial resolution due to its 100 um pixel size and electronic charge sharing correction, making it a good candidate for small ROI imaging. With this detector, “High Speed Angiography” (HSA) was performed on a variety of 3D printed patient-specific vasculature and interventional devices, using a simulated flow loop and iodinated contrast media. The images, which illustrate pathology-dependent flow detail, were recorded in a sequence of 1 ms frames. In addition, the energy discrimination capabilities of the Actaeon were used such that with a lower energy threshold, instrumentation noise was virtually negligible. The per frame noise quality and overall patient dose were acceptable as compared to standard angiography dose rates using FPDs. The previously unseen flow detail may give new insight into the diagnosis, progression, and treatment of neuro and cardiovascular pathologies.

BackgroundA 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.ObjectiveWe 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).MethodsDuring 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.ResultsTwenty-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.cm2, and the mean x-ray exposure time was 39.30 min.ConclusionsThe Hi-Def mode improves visualization of flow diverters and may help in achieving more accurate placement and deployment of devices.

BACKGROUND AND PURPOSE: Visualization of structural details of treatment devices during neurointerventional procedures can be challenging. A new true 2-resolution imaging x-ray detector system (Canon Medical Systems Corporation, Tochigi, Japan) 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 1 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 (PED; Medtronic, Dublin, Ireland) and 1 patient with a low-profile visualized intramural support (LVIS Blue) device (MicroVention-Terumo, Somerset, New Jersey) 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.

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.

Purpose: Various CT-FFR methods are being proposed as a non-invasive method to estimate cardiac disease severity. 3D printed patient specific cardiovascular models with high geometric accuracy can be used to simulate blood flow conditions and perform precise and repeatable benchtop flow experiments for validation of such methods. Materials and Methods: Twelve patient-specific 3D printed cardiac phantoms were created from CT Angiography (CTA) scans using a compliant 3D printing material. Pressure sensors were connected to the aortic root and distal ends of the three main coronary arteries to measure benchtop pressure gradients for each stenosed vessel. The patient geometries were used in Canon Medical Systems 1D fluid dynamics algorithm to calculate the CT- FFR. In addition, a 3D computational fluid dynamics simulation was done using ANSYS to estimate pressure gradients across the coronary arteries. Experimental data and 1D and 3D flow simulations were compared to the standard catheter lab FFR measurement (Invasive-FFR). Results: The average percent difference in Benchtop FFR/Invasive FFR, CT-FFR/Invasive FFR, and CFDFFR/Invasive FFR was 0.05, 0.06, and 0.07 respectively. The average time it took for the CT-FFR simulation was ~35 minutes and it took ~15 hours for the CFD-FFR simulation but can vary based on the number of iterations the user defines the software to run. Conclusions: Benchtop FFR proved to be highly accurate when compared to both 1D and 3D CFD software and therefore, 3D printing of patient specific coronary phantoms is a quality tool for CT-FFR software validation.

Region-of-Interest (ROI) fluoroscopy uses a differential x-ray attenuator to reduce the dose to the patient in the periphery region while maintaining regular dose within the ROI treatment area, resulting in an image with differential brightness regions. The brightness difference can be corrected by subtracting a mask of the ROI attenuator from the dose-reduced image. The purpose of this work is to implement a Convolutional-Neural-Network (CNN) capable of deriving a mask of the ROI attenuator from the dose-reduced images, which can be used to equalize the brightness in the dose-reduced images without a pre-acquired mask. A data set of 10,000 ROI dose reduced images of various objects including anthropomorphic head and chest phantoms were generated with different ROI positions and sizes. A 22 layer CNN was developed to derive a mask of the ROI attenuator from the dose reduced image. The network was trained on 30% and tested on 15% of the images from the ROI image data set. The trained CNN was used on the remaining 55% of the data set to generate the ROI mask, and the average computation time for each image was calculated to be 70 ms. The Mean-Square-Error (MSE) for the testing data set was calculated to be 2.03e-05. For the remaining data set of 5500 images the average MSE between the network output and the corresponding expected ROI mask was calculated to be 2.21e-05. The masks generated by the CNN were successfully used to restore and equalize the brightness of the dose reduced images.

We have been working on the development of a Scatter Display System (SDS) for monitoring and displaying scatterrelated dose to staff members during fluoroscopic interventional procedures. We have considered various methods for such a display using augmented reality (AR) and computer-generated virtual reality (VR). The current work focuses on development of the VR SDS display, which shows the color-coded scattered dose distribution in a horizontal plane at a selected height above the floor in a top-down view of the interventional suite. Reported is the first development of the methodology for real-time functionality of this software via integration of controller area network (CAN) bus digital signals from the Canon C-Arm Biplane System. Importing the CAN bus information allows immediate selection of the appropriate pre-calculated scatter dose distribution consistent with the x-ray beam orientation and characteristics as well as selection of the proper gantry and table graphic for the display. The Python CAN interface module was used for streamlining the integration of the CAN bus interface. Development of real-time functionality for the SDS allows it to provide feedback to staff during clinical procedures for informed dose management; the SDS can work alongside the patient skin dose tracking system (DTS) for complete clinical monitoring of staff and patient dose.

Purpose: The purpose of this study is to apply targeted Parametric Imaging on aneurysms to quantitatively investigate contrast flow changes at pre-, post-treatment and follow-up with outcome scoring. Methods: The angiograms for 50 patients were acquired, 25 treated with coil embolization and 25 treated using a flow diverter. API was performed by synthesizing the time density curve (TDC) at every pixel. Based on the TDCs, we calculated various parameters for the quantitative characterization of contrast flow through the vascular network and aneurysms and displayed them using color encoded maps. The parameters included were : Time to Peak (TTP), Mean Transit Time (MTT), Time of Arrival (TTA), Peak Height (PH) and Area Under the Curve (AUC). Two Regions of Interest (ROI) were manually marked over the aneurysm dome and the main artery. Average aneurysm parameter values were normalized to those values recorded in the main artery and recorded pre-/post-treatment and follow-up and compared to Raymond Roy scores and flow diverter stent scoring. Results: The normalized mean values were as follows (pre and post treatment): TTP (1.09+/-0.14, 1.55+/-1.36), MTT (1.07+/-0.23, 1.27+/-0.42), TTA (0.14+/-0.15, 0.26+/-0.23), PH (1.2+/-0.54, 0.95+/-0.83) and AUC (1.29+/-0.69, 1.44+/- 1.92). The neural network gave a validation accuracy of 0.8036 with a loss of 0.0927. A receiver operating characteristic curve with an AUC of 0.866 was obtained. Conclusions: API can quantitatively describe the flow in the aneurysm for initial investigation of the radiomics of intracranial aneurysms. It also shows a clear demarcation between pre and post treatment. Statistical modelling and a machine learning network is used to prove the success of our model.

Purpose: Neurovascular interventional procedures using biplane fluoroscopic imaging systems can lead to increased risk of radiation-induced skin injuries. The authors developed a biplane dose tracking system (Biplane-DTS) to calculate the cumulative skin dose distribution from the frontal and lateral x-ray tubes and display it in real-time as a color-coded map on a 3D graphic of the patient for immediate feedback to the physician. The agreement of the calculated values with the dose measured on phantoms was evaluated. Methods: The Biplane-DTS consists of multiple components including 3D graphic models of the imaging system and patient, an interactive graphical user interface, a data acquisition module to collect geometry and exposure parameters, the computer graphics processing unit, and functions for determining which parts of the patient graphic skin surface are within the beam and for calculating dose. The dose is calculated to individual points on the patient graphic using premeasured calibration files of entrance skin dose per mAs including backscatter; corrections are applied for field area, distance from the focal spot and patient table and pad attenuation when appropriate. The agreement of the calculated patient skin dose and its spatial distribution with measured values was evaluated in 2D and 3D for simulated procedure conditions using a PMMA block phantom and an SK-150 head phantom, respectively. Dose values calculated by the Biplane-DTS were compared to the measurements made on the phantom surface with radiochromic film and a calibrated ionization chamber, which was also used to calibrate the DTS. The agreement with measurements was specifically evaluated with variation in kVp, gantry angle, and field size. Results: The dose tracking system that was developed is able to acquire data from the two x-ray gantries on a biplane imaging system and calculate the skin dose for each exposure pulse to those vertices of a patient graphic that are determined to be in the beam. The calculations are done in real-time with a typical graphic update time of 30 ms and an average vertex separation of 3 mm. With appropriate corrections applied, the Biplane-DTS was able to determine the entrance dose within 6% and the spatial distribution of the dose within 4% compared to the measurements with the ionization chamber and film for the SK150 head phantom. The cumulative dose for overlapping fields from both gantries showed similar agreement. Conclusions: The Biplane-DTS can provide a good estimate of the peak skin dose and cumulative skin dose distribution during biplane neurointerventional procedures. Real-time display of this information should help the physician manage patient dose to reduce the risk of radiation-induced skin injuries.

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

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 × 109 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.

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.

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.

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 10p. 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.

Purpose: 3D printed patient specific vascular models provide the ability to perform precise and repeatable benchtop experiments with simulated physiological blood flow conditions. This approach can be applied to CT-derived patient geometries to determine coronary flow related parameters such as Fractional Flow Reserve (FFR). To demonstrate the utility of this approach we compared bench-top results with non-invasive CT-derived FFR software based on a computational fluid dynamics algorithm and catheter based FFR measurements. Materials and Methods: Twelve patients for whom catheter angiography was clinically indicated signed written informed consent to CT Angiography (CTA) before their standard care that included coronary angiography (ICA) and conventional FFR (Angio-FFR). The research CTA was used first to determine CT-derived FFR (Vital Images) and second to generate patient specific 3D printed models of the aortic root and three main coronary arteries that were connected to a programmable pulsatile pump. Benchtop FFR was derived from pressures measured proximal and distal to coronary stenosis using pressure transducers. Results: All 12 patients completed the clinical study without any complication, and the three FFR techniques (Angio-FFR, CT-FFR, and Benchtop FFR) are reported for one or two main coronary arteries. The Pearson correlation among Benchtop FFR/ Angio-FFR, CT-FFR/ Benchtop FFR, and CT-FFR/ Angio-FFR are 0.871, 0.877, and 0.927 respectively. Conclusions: 3D printed patient specific cardiovascular models successfully simulated hyperemic blood flow conditions, matching invasive Angio-FFR measurements. This benchtop flow system could be used to validate CTderived FFR diagnostic software, alleviating both cost and risk during invasive procedures.

3D printing has been used to create complex arterial phantoms to advance device testing and physiological condition evaluation. Stereolithographic (STL) files of patient-specific cardiovascular anatomy are acquired to build cardiac vasculature through advanced mesh-manipulation techniques. Management of distal branches in the arterial tree is important to make such phantoms practicable. We investigated methods to manage the distal arterial flow resistance and pressure thus creating physiologically and geometrically accurate phantoms that can be used for simulations of image-guided interventional procedures with new devices. Patient specific CT data were imported into a Vital Imaging workstation, segmented, and exported as STL files. Using a mesh-manipulation program (Meshmixer) we created flow models of the coronary tree. Distal arteries were connected to a compliance chamber. The phantom was then printed using a Stratasys Connex3 multimaterial printer: the vessel in TangoPlus and the fluid flow simulation chamber in Vero. The model was connected to a programmable pump and pressure sensors measured flow characteristics through the phantoms. Physiological flow simulations for patient-specific vasculature were done for six cardiac models (three different vasculatures comparing two new designs). For the coronary phantom we obtained physiologically relevant waves which oscillated between 80 and 120 mmHg and a flow rate of ~125 ml/min, within the literature reported values. The pressure wave was similar with those acquired in human patients. Thus we demonstrated that 3D printed phantoms can be used not only to reproduce the correct patient anatomy for device testing in image-guided interventions, but also for physiological simulations. This has great potential to advance treatment assessment and diagnosis.

Purpose: Accurate patient-specific phantoms for device testing or endovascular treatment planning can be 3D printed. We expand the applicability of this approach for cardiovascular disease, in particular, for CT-geometry derived benchtop measurements of Fractional Flow Reserve, the reference standard for determination of significant individual coronary artery atherosclerotic lesions. Materials and Methods: Coronary CT Angiography (CTA) images during a single heartbeat were acquired with a 320x0.5mm detector row scanner (Toshiba Aquilion ONE). These coronary CTA images were used to create 4 patientspecific cardiovascular models with various grades of stenosis: severe

Following new trends in precision medicine, Juxatarenal Abdominal Aortic Aneurysm (JAAA) treatment has been enabled by using patient-specific fenestrated endovascular grafts. The X-ray guided procedure requires precise orientation of multiple modular endografts within the arteries confirmed via radiopaque markers. Patient-specific 3D printed phantoms could familiarize physicians with complex procedures and new devices in a risk-free simulation environment to avoid periprocedural complications and improve training. Using the Vascular Modeling Toolkit (VMTK), 3D Data from a CTA imaging of a patient scheduled for Fenestrated EndoVascular Aortic Repair (FEVAR) was segmented to isolate the aortic lumen, thrombus, and calcifications. A stereolithographic mesh (STL) was generated and then modified in Autodesk MeshMixer for fabrication via a Stratasys Eden 260 printer in a flexible photopolymer to simulate arterial compliance. Fluoroscopic guided simulation of the patient-specific FEVAR procedure was performed by interventionists using all demonstration endografts and accessory devices. Analysis compared treatment strategy between the planned procedure, the simulation procedure, and the patient procedure using a derived scoring scheme. Results: With training on the patient-specific 3D printed AAA phantom, the clinical team optimized their procedural strategy. Anatomical landmarks and all devices were visible under x-ray during the simulation mimicking the clinical environment. The actual patient procedure went without complications. Conclusions: With advances in 3D printing, fabrication of patient specific AAA phantoms is possible. Simulation with 3D printed phantoms shows potential to inform clinical interventional procedures in addition to CTA diagnostic imaging.

X-ray image-guided interventions are increasingly used to treat surgically inaccessible anomalies. Many of these interventions can be prolonged and result in a significant exposure to the patient’s skin.1 The biological damage to the tissue is associated with threshold dose levels for various deterministic effects. A substantial component of dose to the skin is due to backscattered radiation from the patient. Backscatter factors have been measured and calculated by a number of investigators; their results have shown the backscatter factor to be dependent on the size of the field, the thickness of the patient, and the beam energy.2–11 These backscatter factors can be used to convert the incident air kerma to the entrance surface dose at the center of a uniform beam of a given field size since the backscatter factors are determined under these conditions. However, the backscatter contribution to the skin dose will not be uniform across the field from the center to the edge but will decrease toward the periphery of the field and extend beyond the field boundary.12 Furthermore, when the primary beam is not uniform, a single backscatter factor will not be able to describe the distribution of backscatter dose. For example, beam-shaping devices, such as region-of-interest (ROI) attenuators13–15 and compensation filters,16 modulate the photon fluence in the field of view (FOV) by attenuating the x-ray intensity outside the area-of-interest. ROI attenuators result in a dose saving to the patient’s skin, but with increased image noise in the periphery of the x-ray beam and improved image quality in the unattenuated ROI. Likewise, compensation filters in fluoroscopy are used for x-ray signal equalization by compensating for spatially varying patient transmission in the FOV.17 The compensation filter modifies the primary x-ray intensity and results in dose reduction under the filter. The presence of beam-shaping devices modulates the primary dose intensity and thereby the spatial backscatter distribution over the FOV. In this paper, we describe a method to determine the two-dimensional (2-D) skin dose distribution for fluoroscopically guided interventions (FGI) through convolution of a backscatter point spread function (PSFn) with the primary entrance x-ray field.18 Dose distribution determination through the 2-D and three-dimensional (3-D) convolution19,20 of invariant kernels has been reported in radiotherapy for treatment planning with irregular shaped photon beams.21 The current work expands the scope of skin dose determination for FGIs through the use of pencil beam backscatter kernel convolution for irregularly shaped beams resulting from beam-shaping devices.

Purpose: We have developed a Micro‐Solid State X‐ray Imaging Intensifier (μSSXII) detector with an electron multiplying (EM)CCD chip which allows for an ultra‐high pixel resolution of 8 × 8 microns. However, the EMCCD chip area of 1004 × 1002 pixels results in an 8.032 × 8.016 mm imaging area and restricts the ability of the user to orient and properly image larger pre‐clinical objects such as rat‐kidney‐vasculature casts. We propose a method to align and assemble images of such larger structures while preserving the ultra‐high spatial resolution of the μSSXII. Methods: An imaging platform of negligible attenuation was attached to a stepper motor giving the platform free movement in the plane normal to the fluoroscopic x‐ray beam. By alternating the detector'simage acquisition and the stepper movement, many adjacent overlapping views of the object are acquired. The user then identifies the two pixel coordinates in adjacent images which represent the same point in space. A custom Matlab computer code uses the pixel coordinates to aggregate and average the input images which results in a larger field of view consisting of many ultra‐high resolution images.Results: The images were able to be successfully combined to form a large image while preserving the ultra‐high resolution of the detector. In experimental tests, multiple portions of a mammography test object were imaged and virtually no spatial resolution degradation was found in the combined image. Additionally, when imaging a resin cast of rat‐kidney‐vasculature, vessels of less than 50 μm could be viewed in the combined image.Conclusions: The setup and method were found to preserve the ultra‐high resolution inherent to the μSSXII while allowing pre‐clinical imaging of objects larger than the detector's field of view. The large field of view was effective in orienting the user towards specific areas of interest in the objects imaged. Supported in part by: NIH Grants R01‐EB008425, R01‐EB002873, and an equipment grant from Toshiba Medical Systems Corp.

Purpose: A brain tissue surrogate material was needed to fill the anatomical cavity of a skull to create a phantom for use in simulated Neuro‐EIGI procedures. To enable diagnostic and interventional procedure simulation, the BIM must fit into and be congruous with the interior surface of the skull, be reusable, and allow the implantation of vascular phantoms. The material must reasonably reproduce the automatic technique parameter selections observed during Neuro‐EIGI procedures. Methods: We formulated a putty‐ like material to be used as the BIM. Its x‐ray attenuation properties were evaluated by comparison of the fluoroscopic and radiographic technique parameters automatically selected for a BIM‐filled skull on a Toshiba Infinix angiographic C‐arm unit to those of a solid anthropomorphic head phantom at various projection angles. The same comparison was made between the skull phantom without BIM in the cavity and the anthropomorphic head phantom. The BIM linear attenuation coefficient was calculated and compared to that of PMMA, a common tissue analog plastic. Results: The BIM keeps its shape, is moldable and reusable, and is congruent to the skull's interior surfaces. It allows for insertion and interchange of various custom vascular phantoms at proper anatomic locations. Addition of the BIM to the skull cavity improves the matching of the automatically selected parameters to those of the anthropomorphic phantom by an average of 96.3% for mAs and by 4.2% for kVp in fluoroscopy mode and by 88.6% and 9.0%, respectively, in DSA mode. The BIM's experimental and theoretical linear attenuation coefficient for the RQA5 spectrum differed from PMMA's by about 30%. Conclusions: Despite the difference in attenuation coefficients between the PMMA and BIM, the BIM is a good surrogate material for Neuro‐EIGI research as shown by its properties of congruity, reusability, and device implantation, along with the demonstrated improvement of automatically selected technique parameters. Supported in part by: NIH Grants R01‐EB008425, R01‐EB002873, and an equipment grant from Toshiba Medical Systems Corp.

This project assessed the effectiveness of using two different detectors to obtain dual-energy (DE) micro-CT data for the carrying out of material decomposition. A micro-CT coupled to either a complementary metal-oxide semiconductor (CMOS) or an electron multiplying CCD (EMCCD) detector was used to acquire image data of a 3D-printed phantom with channels filled with different materials. At any instance, materials such as iohexol contrast agent, water, and platinum were selected to make up the scanned object. DE micro-CT data was acquired, and slices of the scanned object were differentiated by material makeup. The success of the decomposition was assessed quantitatively through the computation of percentage normalized root-mean-square error (%NRMSE). Our results indicate a successful decomposition of iohexol for both detectors (%NRMSE values of 1.8 for EMCCD, 2.4 for CMOS), as well as platinum (%NRMSE value of 4.7). The CMOS detector performed material decomposition on air and water on average with 7 times more %NRMSE, possibly due to the decreased sensitivity of the CMOS system. Material decomposition showed the potential to differentiate between materials such as the iohexol and platinum, perhaps opening the door for its use in the neurovascular anatomical region. Work supported by Toshiba America Medical Systems, and partially supported by NIH grant 2R01EB002873.

In some medical-imaging procedures using CBCT and fluoroscopy, it may be needed to visualize only the center of the field-of-view with optimal quality. To reduce the dose to the patient as well as enable increased contrast in the region of interest (ROI) during CBCT and fluoroscopy procedures, a 0.7 mm thick Cu ROI attenuator with a circular aperture 12% of the FOV was used. The aim of this study was to quantify the dose-reduction benefit of ROI imaging during a typical CBCT and interventional fluoroscopy procedures in the head and torso. The Toshiba Infinix C-Arm System was modeled in BEAMnrc/EGSnrc with and without the ROI attenuator. Patient organ and effective doses were calculated in DOSXYZnrc/EGSnrc Monte-Carlo software for CBCT and interventional procedures. We first compared the entrance dose with and without the ROI attenuator on a 20 cm thick solid-water block. Then we simulated a CBCT scan and an interventional fluoroscopy procedure on the head and torso with and without an ROI attenuator. The results showed that the entrance-surface dose reduction in the solid water is about 85.7% outside the ROI opening and 10.5% in the ROI opening. The results showed a reduction in most organ doses of 45%-70% and in effective dose of 46%-66% compared to the dose in a CBCT scan and in an interventional procedure without the ROI attenuator. This work provides evidence of substantial reduction of organ and effective doses when using an ROI attenuator during CBCT and fluoroscopic procedures.

Modern 3D printing technology can fabricate vascular phantoms based on an actual human patient with a high degree of precision facilitating a realistic simulation environment for an intervention. We present two experimental setups using 3D printed patient-specific neurovasculature to simulate different disease anatomies. To simulate the human neurovasculature in the Circle of Willis, patient-based phantoms with aneurysms were 3D printed using a Objet Eden 260V printer. Anthropomorphic head phantoms and a human skull combined with acrylic plates simulated human head bone anatomy and x-ray attenuation. For dynamic studies the 3D printed phantom was connected to a pulsatile flow loop with the anthropomorphic phantom underneath. By combining different 3D printed phantoms and the anthropomorphic phantoms, different patient pathologies can be simulated. For static studies a 3D printed neurovascular phantom was embedded inside a human skull and used as a positional reference for treatment devices such as stents. To simulate tissue attenuation acrylic layers were added. Different combinations can simulate different patient treatment procedures. The Complementary-Metal-Oxide-Semiconductor (CMOS) based High Resolution Fluoroscope (HRF) with 75μm pixels offers an advantage over the state-of-the-art 200 μm pixel Flat Panel Detector (FPD) due to higher Nyquist frequency and better DQE performance. Whether this advantage is clinically useful during an actual clinical neurovascular intervention can be addressed by qualitatively evaluating images from a cohort of various cases performed using both detectors. The above-mentioned method can offer a realistic substitute for an actual clinical procedure. Also a large cohort of cases can be generated and used for a HRF clinical utility determination study.

Background Artificial stroke models can be used for testing various thrombectomy devices. Objective To determine the value of combined stentriever–aspiration thrombectomy compared with the stentriever-alone approach. Methods We designed an in vitro model of the intracranial circulation with a focus on the middle cerebral artery (MCA) that closely resembles the human intracranial circulation. After introducing fresh clot in the MCA, we used conventional biplane angiography and microangiographic fluoroscopy to compare recanalization rates and occurrence of emboli in new, unaffected territory for thrombectomy approaches in which a stentriever (Solitaire flow restoration stentriever, Covidien) was used alone or in combination with continuous manual aspiration through a Navien catheter (Covidien). Results In a total of 22 experiments (11 for each approach), successful clot delivery to the MCA was achieved in all cases. Successful angiographic recanalization (thrombolysis in cerebral infarction score of 2b–3) was achieved more frequently with the combined stentriever–aspiration approach than with the stentriever-alone approach (in 10 vs 4 experiments, p=0.023). Emboli in new territory occurred in three experiments with the stentriever-alone approach, and none were seen with the combined approach (p=0.21). Conclusions The combined stentriever–aspiration approach to thrombectomy leads to better angiographic recanalization rates than use of the stentriever alone. Further experiments are needed to test the value of balloon-guide catheters and aspiration performed using other types of catheters and modes of aspiration.

3D printing an anatomically accurate, functional flow loop phantom of a patient’s cardiac vasculature was used to assist in the surgical planning of one of the first native transcatheter mitral valve replacement (TMVR) procedures. CTA scans were acquired from a patient about to undergo the first minimally-invasive native TMVR procedure at the Gates Vascular Institute in Buffalo, NY. A python scripting library, the Vascular Modeling Toolkit (VMTK), was used to segment the 3D geometry of the patient’s cardiac chambers and mitral valve with severe stenosis, calcific in nature. A stereolithographic (STL) mesh was generated and AutoDesk Meshmixer was used to transform the vascular surface into a functioning closed flow loop. A Stratasys Objet 500 Connex3 multi-material printer was used to fabricate the phantom with distinguishable material features of the vasculature and calcified valve. The interventional team performed a mock procedure on the phantom, embedding valve cages in the model and imaging the phantom with a Toshiba Infinix INFX-8000V 5-axis Carm bi-Plane angiography system. Results: After performing the mock-procedure on the cardiac phantom, the cardiologists optimized their transapical surgical approach. The mitral valve stenosis and calcification were clearly visible. The phantom was used to inform the sizing of the valve to be implanted. Conclusion: With advances in image processing and 3D printing technology, it is possible to create realistic patientspecific phantoms which can act as a guide for the interventional team. Using 3D printed phantoms as a valve sizing method shows potential as a more informative technique than typical CTA reconstruction alone.

The skin-dose tracking system (DTS) provides a color-coded illustration of the cumulative skin-dose distribution on a closely-matching 3D graphic of the patient during fluoroscopic interventions in real-time for immediate feedback to the interventionist. The skin-dose tracking utility of DTS has been extended to include cone-beam computed tomography (CBCT) of neurointerventions. While the DTS was developed to track the entrance skin dose including backscatter, a significant part of the dose in CBCT is contributed by exit primary radiation and scatter due to the many overlapping projections during the rotational scan. The variation of backscatter inside and outside the collimated beam was measured with radiochromic film and a curve was fit to obtain a scatter spread function that could be applied in the DTS. Likewise, the exit dose distribution was measured with radiochromic film for a single projection and a correction factor was determined as a function of path length through the head. Both of these sources of skin dose are added for every projection in the CBCT scan to obtain a total dose mapping over the patient graphic. Results show the backscatter to follow a sigmoidal falloff near the edge of the beam, extending outside the beam as far as 8 cm. The exit dose measured for a cylindrical CTDI phantom was nearly 10 % of the entrance peak skin dose for the central ray. The dose mapping performed by the DTS for a CBCT scan was compared to that measured with radiochromic film and a CTDI-head phantom with good agreement.

An evaluation of the relation between parametric imaging results obtained from Digital Subtraction Angiography (DSA) images and blood-flow velocity measured using Doppler ultrasound in patient-specific neurovascular phantoms is provided. A silicone neurovascular phantom containing internal carotid artery, middle cerebral artery and anterior communicating artery was embedded in a tissue equivalent gel. The gel prevented movement of the vessels when blood mimicking fluid was pumped through it to obtain Colour Doppler images. The phantom was connected to a peristaltic pump, simulating physiological flow conditions. To obtain the parametric images, water was pumped through the phantom at various flow rates (100, 120 and 160 ml/min) and 10 ml contrast boluses were injected. DSA images were obtained at 10 frames/sec from the Toshiba C-arm and DSA image sequences were input into LabVIEW software to get parametric maps from time-density curves. The parametric maps were compared with velocities determined by Doppler ultrasound at the internal carotid artery. The velocities measured by the Doppler ultrasound were 38, 48 and 65 cm/s for flow rates of 100, 120 and 160 ml/min, respectively. For the 20% increase in flow rate, the percentage change of blood velocity measured by Doppler ultrasound was 26.3%. Correspondingly, there was a 20% decrease of Bolus Arrival Time (BAT) and 14.3% decrease of Mean Transit Time (MTT), showing strong inverse correlation with Doppler measured velocity. The parametric imaging parameters are quite sensitive to velocity changes and are well correlated to the velocities measured by Doppler ultrasound.

The skin dose tracking system (DTS) that we developed provides a color-coded mapping of the cumulative skin dose distribution on a 3D graphic of the patient during fluoroscopic procedures in real time. The DTS has now been modified to also calculate the kerma area product (KAP) and cumulative air kerma (CAK) for fluoroscopic interventions using data obtained in real-time from the digital bus on a Toshiba Infinix system. KAP is the integral of air kerma over the beam area and is typically measured with a large-area transmission ionization chamber incorporated into the collimator assembly. In this software, KAP is automatically determined for each x-ray pulse as the product of the air kerma/ mAs from a calibration file for the given kVp and beam filtration times the mAs per pulse times the length and width of the beam times a field nonuniformity correction factor. Field nonuniformity is primarily the result of the heel effect and the correction factor was determined from the beam profile measured using radio-chromic film. Dividing the KAP by the beam area at the interventional reference point provides the area-averaged CAK. The KAP and CAK per x-ray pulse are summed after each pulse to obtain the total procedure values in real-time. The calculated KAP and CAK were compared to the values displayed by the fluoroscopy machine with excellent agreement. The DTS now is able to automatically calculate both KAP and CAK without the need for measurement by an add-on transmission ionization chamber.

Because x-ray based image-guided vascular interventions are minimally invasive they are currently the most preferred method of treating disorders such as stroke, arterial stenosis, and aneurysms; however, the x-ray exposure to the patient during long image-guided interventional procedures could cause harmful effects such as cancer in the long run and even tissue damage in the short term. ROI fluoroscopy reduces patient dose by differentially attenuating the incident x-rays outside the region-of-interest. To reduce the noise in the dose-reduced regions previously recursive temporal filtering was successfully demonstrated for neurovascular interventions. However, in cardiac interventions, anatomical motion is significant and excessive recursive filtering could cause blur. In this work the effects of three noise-reduction schemes, including recursive temporal filtering, spatial mean filtering, and a combination of spatial and recursive temporal filtering, were investigated in a simulated ROI dose-reduced cardiac intervention. First a model to simulate the aortic arch and its movement was built. A coronary stent was used to simulate a bio-prosthetic valve used in TAVR procedures and was deployed under dose-reduced ROI fluoroscopy during the simulated heart motion. The images were then retrospectively processed for noise reduction in the periphery, using recursive temporal filtering, spatial filtering and a combination of both. Quantitative metrics for all three noise reduction schemes are calculated and are presented as results. From these it can be concluded that with significant anatomical motion, a combination of spatial and recursive temporal filtering scheme is best suited for reducing the excess quantum noise in the periphery. This new noise-reduction technique in combination with ROI fluoroscopy has the potential for substantial patient-dose savings in cardiac interventions.

X-ray and radionuclide imaging are widely popular medical imaging modalities. The requirements for each imaging modality are different, and different detectors are normally used. For this paper, we demonstrate our new detector capable of both fluoroscopy and angiography, to be used as an imager for both single-photon and integral-energy imaging and applied to the dual modalities of radionuclide imaging and X-ray imaging. This newly developed micro-angiographic fluoroscope (MAF) has 1024 × 1024 pixels of 35 μm effective size and is capable of real-time imaging at 30 fps. The large variable gain of its light image intensifier (LII) provides quantum-limited operation with essentially no additive instrumentation noise. We demonstrate that the MAF can be operated in single-photon counting (SPC) mode for X-ray imaging with substantially better resolution than in energy integration (EI) mode. We may use high LII gain with very low exposure (less than 1 X-ray photon/pixel) per frame for SPC mode (with X-ray and radionuclide) and higher exposure with lower gain for EI mode (transmission imaging with X-rays). For a demonstration of the operation in both EI and SPC mode, a heavily K-edge filtered X-ray beam (average energy of 31 keV) was used to provide a nearly monochromatic spectrum. The MTF measured using a standard slit method showed a dramatic improvement for the SPC mode over the EI mode at all spatial frequencies. Images of a line-pair phantom also showed improved spatial resolution in SPC mode compared to EI mode. In SPC mode, images of human distal and middle phalanges showed the trabecular structures of the bone with far better contrast and detail. We also show MAF operation in SPC mode for radionuclide imaging using a custom-built phantom filled with I-125. A 1-mm-diameter parallel hole, medium-energy gamma camera collimator was placed between the phantom, and the MAF and was moved multiple times at equal intervals in random directions to eliminate the pattern corresponding to the collimator septa. Data was acquired at 20 fps, and multiple signal-thresholded frames were summed in SPC mode to provide an integrated frame. The sharpness of the emission image is limited by the collimator resolution and could be improved by optimized collimator design. We demonstrate that the same MAF is capable of operating in both SPC and EI modes and can be used in both X-ray transmission imaging and radionuclide emission imaging.

The self-expanding V-POD is a second-generation flow-diverting device with a low-porosity PTFE patch on a self-expanding microstent. The authors evaluated this device for the treatment of elastase-induced aneurysms in rabbits.Three V-POD types (A, circumferential patch closed-cell stent [n = 9]; B, asymmetric patch closed-cell stent [n = 7]; and C, asymmetric patch open-cell stent [n = 4]) were evaluated by using angiography, conebeam micro-CT, histology, and SEM. Aneurysm flow modifications were expressed in terms of immediate poststent/prestent ratios of maximum CA volume entering the aneurysm dome tracked on procedural angiograms. Flow modifications were correlated with 4 weeks' follow-up angiographic, micro-CT, histologic, and SEM results.Mechanical stent-deployment difficulties in 4 aneurysms (1 type A; 3 type B) led to suboptimal results and exclusion from analysis. Of the remaining 16 aneurysms, 4-week post-treatment angiograms showed no aneurysm filling in 10 (63%), 3 (∼19%) had no filling with a small remnant neck, and 3 (∼19%) had0.25 filling. Successfully treated aneurysms (n = 16) demonstrated an immediate poststent/prestent CA maximum volume ratio of 0.13 ± 0.18% (0.0%-0.5%). Favorable contrast-flow modification on immediate angiography after deployment correlated significantly with aneurysm occlusion on follow-up angiography, micro-CT, and histology. The occlusion percentage derived from micro-CT was 96 ± 6.8%. Histology indicated advanced healing (grade ≥3) in the aneurysm dome in 13 of 16 cases. SEM revealed 15 of 16 stents in an advanced state of endothelialization.This study showed the feasibility and effectiveness of V-POD for aneurysm healing in a rabbit elastase model.

Partially polyurethane-covered stent (PPCS) is proposed for the treatment of cerebral aneurysms. The PPCSs were observed to substantially modify the flow entering the aneurysm in a patient-specific aneurysm phantom (PSAP). These stents can act as flow modulators and the polyurethane (PU) membrane can provide a smooth scaffold for restoring the structural integrity of the diseased vessel. Partial coating of the stent aids in sealing only the entrance to the aneurysm while keeping the perforators around the aneurysm open and patent. Biocompatibility of the PU membrane was monitored using contact angle measurements to show that critical surface tension (CST) values remained in the thromboresistant range of 20–30 mN/m. Stent flexibility, stiffness, and pressure–diameter relationship showed no significant change after asymmetric PU film application. No delamination of the PU membrane from the stent was observed within the working strains of the stent. The flow modulating capability of the PPCS was monitored by intentionally orienting the stent to cover either the proximal or the distal regions along the neck of the PSAP. Time density curves (TDCs) compared the relative metrics of input rate, washout rate, residence time, and influx in the aneurysm before and after the stent placement.

Minimally invasive interventions are rapidly replacing invasive surgical procedures for the most prevalent human disease conditions. X-ray image-guided interventions carried out using the insertion and navigation of catheters through the vasculature are increasing in number and sophistication. In this article, we offer our vision for the future of this dynamic field of endovascular image-guided interventions in the form of predictions about (1) improvements in high-resolution detectors for more accurate guidance, (2) the implementation of high-resolution region of interest computed tomography for evaluation and planning, (3) the implementation of dose tracking systems to control patient radiation risk, (4) the development of increasingly sophisticated interventional devices, (5) the use of quantitative treatment planning with patient-specific computer fluid dynamic simulations, and (6) the new expanding role of the medical physicist. We discuss how we envision our predictions will come to fruition and result in the universal goal of improved patient care.

Currently, a large number of endovascular interventions are performed for treatment of intracranial aneurysms. For these treatments, correct positioning of microcatheter tips, microguide wire tips, or coils is essential. Techniques to detect such devices may facilitate endovascular interventions. In this paper, we describe an algorithm for tracking of microcatheter tips during fluoroscopically guided neuroendovascular interventions. A sequence of fluoroscopic images (1,024 × 1,024 × 12 bits) was acquired using a C-arm angiography system as a microcatheter was passed through a carotid phantom which was on top of a head phantom. The carotid phantom was a silicone cylinder containing a simulated vessel with the shape and curvatures of the internal carotid artery. The head phantom consisted of a human skull and tissue-equivalent material. To detect the microcatheter in a given fluoroscopic frame, a background image consisting of an average of the four previous frames is subtracted from the current frame, the resulting image is filtered using a matched filter, and the position of maximum intensity in the filtered image is taken as the catheter tip position in the current frame. The distance between the tracked position and the correct position (error distance) was measured in each of the fluoroscopic images. The mean and standard deviation of the error distance values were 0.277 mm (1.59 pixels) and 0.26 mm (1.5 pixels), respectively. The error distance was less than 3 pixels in the 93.0% frames. Although the algorithm intermittently failed to correctly detect the catheter, the algorithm recovered the catheter in subsequent frames.

Minimally invasive endovascular image-guided interventions (EIGIs) are the preferred procedures for treatment of a wide range of vascular disorders. Despite benefits including reduced trauma and recovery time, EIGIs have their own challenges. Remote catheter actuation and challenging anatomical morphology may lead to erroneous endovascular device selections, delays or even complications such as vessel injury. EIGI planning using 3D phantoms would allow interventionists to become familiarized with the patient vessel anatomy by first performing the planned treatment on a phantom under standard operating protocols. In this study the optimal workflow to obtain such phantoms from 3D data for interventionist to practice on prior to an actual procedure was investigated. Patient-specific phantoms and phantoms presenting a wide range of challenging geometries were created. Computed Tomographic Angiography (CTA) data was uploaded into a Vitrea 3D station which allows segmentation and resulting stereo-lithographic files to be exported. The files were uploaded using processing software where preloaded vessel structures were included to create a closed-flow vasculature having structural support. The final file was printed, cleaned, connected to a flow loop and placed in an angiographic room for EIGI practice. Various Circle of Willis and cardiac arterial geometries were used. The phantoms were tested for ischemic stroke treatment, distal catheter navigation, aneurysm stenting and cardiac imaging under angiographic guidance. This method should allow for adjustments to treatment plans to be made before the patient is actually in the procedure room and enabling reduced risk of peri-operative complications or delays.

Cone beam computed tomography (CBCT) systems with rotational gantries that have standard flat panel detectors (FPD) are widely used for the 3D rendering of vascular structures using Feldkamp cone beam reconstruction algorithms. One of the inherent limitations of these systems is limited resolution (

Digital Subtraction Angiography (DSA) is the main diagnostic tool for intracranial aneurysms (IA) flow-diverter (FD) assisted treatment. Based on qualitative contrast flow evaluation, interventionists decide on subsequent steps. We developed a novel fully Retrievable Asymmetric Flow-Diverter (RAFD) which allows controlled deployment, repositioning and detachment achieve optimal flow diversion. The device has a small low porosity or solid region which is placed such that it would achieve maximum aneurysmal in-jet flow deflection with minimum impairment to adjacent vessels. We tested the new RAFD using a flow-loop with an idealized and a patient specific IA phantom in carotid-relevant physiological conditions. We positioned the deflection region at three locations: distally, center and proximally to the aneurysm orifice and analyzed aneurysm dome flow using DSA derived maps for mean transit time (MTT) and bolus arrival times (BAT). Comparison between treated and untreated (control) maps quantified the RAFD positioning effect. Average MTT, related to contrast presence in the aneurysm dome increased, indicating flow decoupling between the aneurysm and parent artery. Maximum effect was observed in the center and proximal position (~75%) of aneurysm models depending on their geometry. BAT maps, correlated well with inflow jet direction and magnitude. Reduction and jet dispersion as high as about 50% was observed for various treatments. We demonstrated the use of DSA data to guide the placement of the RAFD and showed that optimum flow diversion within the aneurysm dome is feasible. This could lead to more effective and a safer IA treatment using FDs.

Perfusion imaging is the most applied modality for the assessment of acute stroke. Parameters such as Cerebral Blood Flow (CBF), Cerebral Blood volume (CBV) and Mean Transit Time (MTT) are used to distinguish the tissue infarct core and ischemic penumbra. Due to lack of standardization these parameters vary significantly between vendors and software even when provided with the same data set. There is a critical need to standardize the systems and make them more reliable. We have designed a uniform phantom to test and verify the perfusion systems. We implemented a flow loop with different flow rates (250, 300, 350 ml/min) and injected the same amount of contrast. The images of the phantom were acquired using a Digital Angiographic system. Since this phantom is uniform, projection images obtained using DSA is sufficient for initial validation. To validate the phantom we measured the contrast concentration at three regions of interest (arterial input, venous output, perfused area) and derived time density curves (TDC). We then calculated the maximum slope, area under the TDCs and flow. The maximum slope calculations were linearly increasing with increase in flow rate, the area under the curve decreases with increase in flow rate. There was 25% error between the calculated flow and measured flow. The derived TDCs were clinically relevant and the calculated flow, maximum slope and areas under the curve were sensitive to the measured flow. We have created a systematic way to calibrate existing perfusion systems and assess their reliability.