Barium sulfate and pigment admixture for photoacoustic and x-ray contrast imaging of the gut.

Significance: X-ray imaging is frequently used for gastrointestinal imaging. Photoacoustic imaging (PAI) of the gastrointestinal tract is an emerging approach that has been demonstrated for preclinical imaging of small animals. A contrast agent active in both modalities could be useful for imaging applications. Aim: We aimed to develop a dual-modality contrast agent comprising an admixture of barium sulfate with pigments that absorb light in the second near-infrared region (NIR-II), for preclinical imaging with both x-ray and PAI modalities. Approach: Eleven different NIR-II dyes were evaluated after admixture with a 40% w/v barium sulfate mixture. The resulting NIR-II absorption in the soluble fraction and in the total mixture was characterized. Proof-of-principle imaging studies in mice were carried out. Results: Pigments that produced more uniform suspensions were assessed further for photoacoustic contrast signal at a wavelength of 1064 nm that corresponds to the output of the Nd:YAG laser used. Phantom imaging studies demonstrated that the pigment-barium sulfate mixture generated imaging contrast in both x-ray and PAI modalities. The optimal pigment selected for further study was a cyanine tetrafluoroborate salt. Ex-vivo and whole-body mouse imaging demonstrated that photoacoustic and x-ray contrast signals co-localized in the intestines for both imaging modalities. Conclusion: These data demonstrate that commercially-available NIR-II pigments can simply be admixed with barium sulfate to generate a dual-modality contrast agent appropriate for small animal gastrointestinal imaging.

Initial investigation of predicting hematoma expansion for intracerebral hemorrhage using imaging biomarkers and machine learning.

Purpose: Intracerebral Hemorrhage (ICH) is one of the most devastating types of strokes with mortality and morbidity rates ranging from about 51%-65% one year after diagnosis. Early hematoma expansion (HE) is a known cause of worsening neurological status of ICH patients. The goal of this study was to investigate whether non-contrast computed tomography imaging biomarkers (NCCT-IB) acquired at initial presentation can predict ICH growth in the acute stage. Materials and Methods: We retrospectively collected NCCT data from 200 patients with acute (<6 hours) ICH. Four NCCT-IBs (blending region, dark hole, island, and edema) were identified for each hematoma, respectively. HE status was recorded based on the clinical observation reported in the patient chart. Supervised machine learning models were developed, trained, and tested for 15 different input combinations of the NCCT-IBs to predict HE. Model performance was assessed using area under the receiver operating characteristic curve and probability for accurate diagnosis (PAD) was calculated. A 20-fold Monte-Carlo cross validation was implemented to ensure model reliability on a limited sample size of data, by running a myriad of random training/testing splits. Results: The developed algorithm was able to predict expansion utilizing all four inputs with an accuracy of 70.17%. Further testing of all biomarker combinations yielded P AD ranging from 0.57, to 0.70. Conclusion: Specific attributes of ICHs may influence the likelihood of HE and can be evaluated via a machine learning algorithm. However, certain parameters may differ in importance to reach accurate conclusions about potential expansion.

Quantitative angiography prognosis of intracranial aneurysm treatment failure using parametric imaging and distal vessel analysis.

Purpose: Data-driven methods based on x-ray angiographic parametric imaging (API) have been successfully used to provide prognosis for intracranial aneurysm (IA) treatment outcome. Previous studies have mainly focused on embolization devices where the flow pattern visualization is in the aneurysm dome; however, this is not possible in IAs treated with endovascular coils due to high x-ray attenuation of the devices. To circumvent this challenge, we propose to investigate whether flow changes in the parent artery distal to the coil-embolized IAs could be used to achieve the same accuracy of surgical outcome prognosis. Methods: Eighty digital subtraction angiography sequences were acquired from patients with IA embolized with coils. Five API parameters were recorded from a region of interest (ROI) placed distal to the IA neck in the main artery. Average API values were recorded and pre-treatment values. A supervised machine learning algorithm was trained to provide a six-month post procedure binary outcome (occluded/not occluded). Receiver operating characteristic (ROC) analysis was used to assess the accuracy of the method. Results: Use of API parameters with data driven methods yielded an area under the ROC curve of 0.77 ±0.11 and accuracy of 78.6%. Single parameter-based analysis yielded accuracies which were suboptimal for clinical acceptance. Conclusions: We determined that data-driven method based on API analysis of flow in the parent artery of IA treated with coils provide clinically acceptable accuracy for the prognosis of six months occlusion outcome.

Prognosis of ischemia recurrence in patients with moderate intracranial atherosclerotic disease using quantitative MRA measurements.

Purpose: To investigate the relation between delayed ischemic stroke and the intracranial atherosclerotic disease (ICAD) hemodynamics as determined by Non-invasive Optimal Vessel Analysis (NOVA) MRI measurements. Materials and Methods: Thirty-three patients with ICAD were enrolled in this study. All patients underwent clinically indicated angioplasty followed by 2-dimensional phase contrast MR (2D PCMR) performed on a 3.0 Tesla MRI scanner using either a 16-channel neurovascular coil or 32-channel head coil. The volumetric flow rate measurements were calculated from 2D PCMR with Non-invasive Optimal Vessel Analysis (NOVA) software (VasSol, Chicago, IL, USA). Flow rate measurements were obtained in 20 major arteries distal, proximal and within the Circle of Willis. Patients were followed up for six month, and ischemia reoccurrence and location were recorded. Receiver operating characteristic (ROC) analysis was performed using flow rates measurements in the ipsilateral side of the ischemic event occurrence. Results: Complete set of measurements was achieved in n=34. Left and right hemisphere ischemia recurrence was observed in seven and three cases respectively. Best predictor of ischemic event reoccurrence was flow rate in the middle cerebral artery with area under the ROC of 0.821±0.109. Conclusions: This is an effectiveness study to determine whether blood flow measurements in the intracranial vasculature may be predictive of future ischemic events. Our results demonstrated significant correlation between the blood flow measurements using 2D PCMR processed with the NOVA software and the reoccurrence of ischemia. These results support further investigation for using this method for risk stratification of ICAD patients.

Investigation of the efficacy of a data-driven CT artifact correction scheme for sparse and truncated projection data for intracranial hemorrhage diagnosis.

Data-driven CT-image reconstruction techniques for truncated or sparsely acquired projection data to reduce radiation dose, iodine volume, and patient motion artifacts have been proposed. These approaches have shown good performance and preservation of image quality metrics. To continue these efforts, we investigated whether these techniques affect the performance of a machine-learning algorithm to identify the presence of intracranial hemorrhage (ICH). Ten-thousand head CT scans were collected from the 2019 RSNA Intracranial Hemorrhage Detection and Classification Challenge dataset. Sinograms were simulated and then resampled in both a one-third truncated and one-third sparse manner. GANs were tasked with correcting the incomplete projection data in two ways. Firstly, in the sinogram domain, where the incomplete sinogram was filled by the GAN and then reconstructed. Secondly, in the reconstruction domain, where the incomplete data were first reconstructed and the sparse or truncation artifacts were corrected by the GAN. Eighty-five hundred images were used for artifact correction network training, and 1500 were withheld for network assessment via an already trained machine-learning algorithm tasked with diagnosis of ICH presence. Fully-sampled reconstructions were compared with the sparse and truncated reconstructions for classification accuracy. Difference in classification accuracy between the fully sampled (83.4%), sparse (82.0%), and truncated (82.3%) reconstructions was minimal, demonstrating that the network diagnosis performance is unaffected by 2/3 reduction of projection data. This work indicates that data-driven reconstructions for a sparse or truncated projection dataset can provide high diagnostic performance for ICH detection at a fraction of the typical radiation dose.

Use of quantitative angiographic methods with a data-driven model to evaluate reperfusion status (mTICI) during thrombectomy.

PURPOSE: Intra-procedural assessment of reperfusion during mechanical thrombectomy (MT) for emergent large vessel occlusion (LVO) stroke is traditionally based on subjective evaluation of digital subtraction angiography (DSA). However, semi-quantitative diagnostic tools which encode hemodynamic properties in DSAs, such as angiographic parametric imaging (API), exist and may be used for evaluation of reperfusion during MT. The objective of this study was to use data-driven approaches, such as convolutional neural networks (CNNs) with API maps, to automatically assess reperfusion in the neuro-vasculature during MT procedures based on the modified thrombolysis in cerebral infarction (mTICI) scale. METHODS: DSAs from patients undergoing MTs of anterior circulation LVOs were collected, temporally cropped to isolate late arterial and capillary phases, and quantified using API peak height (PH) maps. PH maps were normalized to reduce injection variability. A CNN was developed, trained, and tested to classify PH maps into 2 outcomes (mTICI 0,1,2a/mTICI 2b,2c,3) or 3 outcomes (mTICI 0,1,2a/mTICI 2b/mTICI 2c,3), respectively. Ensembled networks were used to combine information from multiple views (anteroposterior and lateral). RESULTS: The study included 383 DSAs. For the 2-outcome classification, average accuracy was 81.0% (95% CI, 79.0-82.9%), and the area under the receiver operating characteristic curve (AUROC) was 0.86 (0.84-0.88). For the 3-outcome classification, average accuracy was 64.0% (62.0-66.0), and AUROC values were 0.85 (0.83-0.87), 0.74 (0.71-0.77), and 0.78 (0.76-0.81) for the mTICI 0,1,2a, mTICI 2b, and mTICI 2c,3 classes, respectively. CONCLUSION: This study demonstrated the feasibility of using hemodynamic information in API maps with data-driven models to autonomously assess intra-procedural reperfusion during MT.

Assessment of computed tomography perfusion software in predicting spatial location and volume of infarct in acute ischemic stroke patients: a comparison of Sphere, Vitrea, and RAPID.

BACKGROUND: CT perfusion (CTP) infarct and penumbra estimations determine the eligibility of patients with acute ischemic stroke (AIS) for endovascular intervention. This study aimed to determine volumetric and spatial agreement of predicted RAPID, Vitrea, and Sphere CTP infarct with follow-up fluid attenuation inversion recovery (FLAIR) MRI infarct. METHODS: 108 consecutive patients with AIS and large vessel occlusion were included in the study between April 2019 and January 2020 . Patients were divided into two groups: endovascular intervention (n=58) and conservative treatment (n=50). Intervention patients were treated with mechanical thrombectomy and achieved successful reperfusion (Thrombolysis in Cerebral Infarction 2b/2 c/3) while patients in the conservative treatment group did not receive mechanical thrombectomy or intravenous thrombolysis. Intervention and conservative treatment patients were included to assess infarct and penumbra estimations, respectively. It was assumed that in all patients treated conservatively, penumbra converted to infarct. CTP infarct and penumbra volumes were segmented from RAPID, Vitrea, and Sphere to assess volumetric and spatial agreement with follow-up FLAIR MRI. RESULTS: Mean infarct differences (95% CIs) between each CTP software and FLAIR MRI for each cohort were: intervention cohort: RAPID=9.0±7.7 mL, Sphere=-0.2±8.7 mL, Vitrea=-7.9±8.9 mL; conservative treatment cohort: RAPID=-31.9±21.6 mL, Sphere=-26.8±17.4 mL, Vitrea=-15.3±13.7 mL. Overlap and Dice coefficients for predicted infarct were (overlap, Dice): intervention cohort: RAPID=(0.57, 0.44), Sphere=(0.68, 0.60), Vitrea=(0.70, 0.60); conservative treatment cohort: RAPID=(0.71, 0.56), Sphere=(0.73, 0.60), Vitrea=(0.72, 0.64). CONCLUSIONS: Sphere proved the most accurate in patients who had intervention infarct assessment as Vitrea and RAPID overestimated and underestimated infarct, respectively. Vitrea proved the most accurate in penumbra assessment for patients treated conservatively although all software overestimated penumbra.

CT artifact correction for sparse and truncated projection data using generative adversarial networks.

PURPOSE: Computed tomography image reconstruction using truncated or sparsely acquired projection data to reduce radiation dose, iodine volume, and patient motion artifacts has been widely investigated. To continue these efforts, we investigated the use of machine learning-based reconstruction techniques using deep convolutional generative adversarial networks (DCGANs) and evaluated its effect using standard imaging metrics. METHODS: Ten thousand head computed tomography (CT) scans were collected from the 2019 RSNA Intracranial Hemorrhage Detection and Classification Challenge dataset. Sinograms were simulated and then resampled in both a one-third truncated and one-third sparse manner. DCGANs were tasked with correcting the incomplete projection data, either in the sinogram domain where the full sinogram was recovered by the DCGAN and then reconstructed, or the reconstruction domain where the incomplete data were first reconstructed and the sparse or truncation artifacts were corrected by the DCGAN. Seventy-five hundred images were used for network training and 2500 were withheld for network assessment using mean absolute error (MAE), structural similarity index measure (SSIM), and peak signal-to-noise ratio (PSNR) between results of different correction techniques. Image data from a quality-assurance phantom were also resampled in the two manners and corrected and reconstructed for network performance assessment using line profiles across high-contrast features, the modulation transfer function (MTF), noise power spectrum (NPS), and Hounsfield Unit (HU) linearity analysis. RESULTS: Better agreement with the fully sampled reconstructions were achieved from sparse acquisition corrected in the sinogram domain and the truncated acquisition corrected in the reconstruction domain. MAE, SSIM, and PSNR showed quantitative improvement from the DCGAN correction techniques. HU linearity of the reconstructions was maintained by the correction techniques for the sparse and truncated acquisitions. MTF curves reached the 10% modulation cutoff frequency at 5.86 lp/cm for the truncated corrected reconstruction compared with 2.98 lp/cm for the truncated uncorrected reconstruction, and 5.36 lp/cm for the sparse corrected reconstruction compared with around 2.91 lp/cm for the sparse uncorrected reconstruction. NPS analyses yielded better agreement across a range of frequencies between the resampled corrected phantom and truth reconstructions. CONCLUSIONS: We demonstrated the use of DCGANs for CT-image correction from sparse and truncated simulated projection data, while preserving imaging quality of the fully sampled projection data.

Automatic radiomic feature extraction using deep learning for angiographic parametric imaging of intracranial aneurysms.

BACKGROUND: Angiographic parametric imaging (API) is an imaging method that uses digital subtraction angiography (DSA) to characterize contrast media dynamics throughout the vasculature. This requires manual placement of a region of interest over a lesion (eg, an aneurysm sac) by an operator. OBJECTIVE: The purpose of our work was to determine if a convolutional neural network (CNN) was able to identify and segment the intracranial aneurysm (IA) sac in a DSA and extract API radiomic features with minimal errors compared with human user results. METHODS: Three hundred and fifty angiographic images of IAs were retrospectively collected. The IAs and surrounding vasculature were manually contoured and the masks put to a CNN tasked with semantic segmentation. The CNN segmentations were assessed for accuracy using the Dice similarity coefficient (DSC) and Jaccard index (JI). Area under the receiver operating characteristic curve (AUROC) was computed. API features based on the CNN segmentation were compared with the human user results. RESULTS: The mean JI was 0.823 (95% CI 0.783 to 0.863) for the IA and 0.737 (95% CI 0.682 to 0.792) for the vasculature. The mean DSC was 0.903 (95% CI 0.867 to 0.937) for the IA and 0.849 (95% CI 0.811 to 0.887) for the vasculature. The mean AUROC was 0.791 (95% CI 0.740 to 0.817) for the IA and 0.715 (95% CI 0.678 to 0.733) for the vasculature. All five API features measured inside the predicted masks were within 18% of those measured inside manually contoured masks. CONCLUSIONS: CNN segmentation of IAs and surrounding vasculature from DSA images is non-inferior to manual contours of aneurysms and can be used in parametric imaging procedures.

Feasibility study for use of angiographic parametric imaging and deep neural networks for intracranial aneurysm occlusion prediction.

BACKGROUND: Angiographic parametric imaging (API), based on digital subtraction angiography (DSA), is a quantitative imaging tool that may be used to extract contrast flow parameters related to hemodynamic conditions in abnormal pathologies such as intracranial aneurysms (IAs). OBJECTIVE: To investigate the feasibility of using deep neural networks (DNNs) and API to predict IA occlusion using pre- and post-intervention DSAs. METHODS: We analyzed DSA images of IAs pre- and post-treatment to extract API parameters in the IA dome and the corresponding main artery (un-normalized data). We implemented a two-step correction to account for injection variability (normalized data) and projection foreshortening (relative data). A DNN was trained to predict a binary IA occlusion outcome: occluded/unoccluded. Network performance was assessed with area under the receiver operating characteristic curve (AUROC) and classification accuracy. To evaluate the effect of the proposed corrections, prediction accuracy analysis was performed after each normalization step. RESULTS: The study included 190 IAs. The mean and median duration between treatment and follow-up was 9.8 and 8.0 months, respectively. For the un-normalized, normalized, and relative subgroups, the DNN average prediction accuracies for IA occlusion were 62.5% (95% CI 60.5% to 64.4%), 70.8% (95% CI 68.2% to 73.4%), and 77.9% (95% CI 76.2% to 79.6%). The average AUROCs for the same subgroups were 0.48 (0.44-0.52), 0.67 (0.61-0.73), and 0.77 (0.74-0.80). CONCLUSIONS: The study demonstrated the feasibility of using API and DNNs to predict IA occlusion using only pre- and post-intervention angiographic information.