Factors Associated With Decreased Accuracy of Modified Thrombolysis in Cerebral Infarct Scoring Among Neurointerventionalists During Thrombectomy.

Background and Purpose: The modified Thrombolysis in Cerebral Infarct (mTICI) score is used to grade angiographic outcome after endovascular thrombectomy. We sought to identify factors that decrease the accuracy of intraprocedural mTICI. Methods: We performed a 2-center retrospective cohort study comparing operator (n=6) mTICI scores to consensus scores from blinded adjudicators. Groups were also assessed by dichotomizing mTICI scores to 0­2a versus 2b­3. Results: One hundred thirty endovascular thrombectomy procedures were included. Operators and adjudicators had a pairwise agreement in 96 cases (73.8%). Krippendorff α was 0.712. Multivariate analysis showed endovascular thrombectomy overnight (odds ratio [OR]=3.84 [95% CI, 1.22­12.1]), lacking frontal (OR, 5.66 [95 CI, 1.36­23.6]), or occipital (OR, 7.18 [95 CI, 2.12­24.3]) region reperfusion, and higher operator mTICI scores (OR, 2.16 [95 CI, 1.16­4.01]) were predictive of incorrectly scoring mTICI intraprocedurally. With dichotomized mTICI scores, increasing number of passes was associated with increased risk of operator error (OR, 1.93 [95 CI, 1.22­3.05]). Conclusions: In our study, mTICI disagreement between operator and adjudicators was observed in 26.2% of cases. Interventions that took place between 22:30 and 4:00, featured frontal or occipital region nonperfusion, higher operator mTICI scores, and increased number of passes had higher odds of intraprocedural mTICI inaccuracy.

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.

Semi-automated measurement of vascular tortuosity and its implications for mechanical thrombectomy performance.

PURPOSE: Few studies have examined the geometry of endovascular mechanical thrombectomy pathways. Here we examine the tortuosity and angulations of catheter pathways from the aortic arch to the termination of the internal carotid artery (ICA) and its association with thrombectomy performance. METHODS: We included 100 consecutive anterior circulation large vessel occlusion thrombectomy patients over 12 months. Computed tomography angiograms (CTA) were used for 3D segmentation of catheter pathway from the aortic arch to ICA termination. Tortuosity index (TI) and angulations of the catheter pathway were measured in a semi-automated fashion. TI and angulation degree were compared between sides and correlated with age and procedural measures. RESULTS: We analyzed 188 catheter pathways in 100 patients. Severe angulation (≤ 30°) was present in 5.8% and 39.4% of common carotid artery (CCA) and extracranial ICA segments, respectively. Five pathways (2.6%) had 360° loop. CCA and extracranial ICA tortuosity had a weak but significant correlation with age (r = 0.17, 0.21, p value = 0.05, 0.02 respectively), time from groin puncture to the site of occlusion (r = 0.29, 0.25, p values = 0.008, 0.026 respectively), and fluoroscopy time (r = 0.022, 0.31, p values = 0.016, 0.001 respectively). There was a significant difference in the pattern of angulation (p value = 0.04) and tortuosity between right and left side in CCA segment (TI = 0.20 ± 0.086 vs. 0.15 ± 0.82, p value < 0.001). CONCLUSIONS: There was a significant difference in CCA angulation between right and left sides. TI of extracranial CCA and ICA correlated with age and influenced time from groin puncture to the occlusion site and total fluoroscopy time.

Antimicrobial Peptide Combined with BMP2-Modified Mesenchymal Stem Cells Promotes Calvarial Repair in an Osteolytic Model.

Repair and regeneration of inflammation-induced bone loss remains a clinical challenge. LL37, an antimicrobial peptide, plays critical roles in cell migration, cytokine production, apoptosis, and angiogenesis. Migration of stem cells to the affected site and promotion of vascularization are essential for tissue engineering therapy, including bone regeneration. However, it is largely unknown whether LL37 affects mesenchymal stem cell (MSC) behavior and bone morphogenetic protein 2 (BMP2)-mediated bone repair during the bone pathologic remodeling process. By performing in vitro and in vivo studies with MSCs and a lipopolysaccharide (LPS)-induced mouse calvarial osteolytic bone defect model, we found that LL37 significantly promotes cell differentiation, migration, and proliferation in both unmodified MSCs and BMP2 gene-modified MSCs. Additionally, LL37 inhibited LPS-induced osteoclast formation and bacterial activity in vitro. Furthermore, the combination of LL37 and BMP2 markedly promoted MSC-mediated angiogenesis and bone repair and regeneration in LPS-induced osteolytic defects in mouse calvaria. These findings demonstrate for the first time that LL37 can be a potential candidate drug for promoting osteogenesis and for inhibiting bacterial growth and osteoclastogenesis, and that the combination of BMP2 and LL37 is ideal for MSC-mediated bone regeneration, especially for inflammation-induced bone loss.

Three-dimensional printing of MRI-visible phantoms and MR image-guided therapy simulation.

PURPOSE: To demonstrate the use of anatomic MRI-visible three-dimensional (3D)-printed phantoms and to assess process accuracy and material MR signal properties. METHODS: A cervical spine model was generated from computed tomography (CT) data and 3D-printed using an MR signal-generating material. Printed phantom accuracy and signal characteristics were assessed using 120 kVp CT and 3 Tesla (T) MR imaging. The MR relaxation rates and diffusion coefficient of the fabricated phantom were measured and 1 H spectra were acquired to provide insight into the nature of the proton signal. Finally, T2 -weighted imaging was performed during cryoablation of the model. RESULTS: The printed model produced a CT signal of 102 ± 8 Hounsfield unit, and an MR signal roughly 1/3rd that of saline in short echo time/short repetition time GRE MRI (456 ± 36 versus 1526 ± 121 arbitrary signal units). Compared with the model designed from the in vivo CT scan, the printed model differed by 0.13 ± 0.11 mm in CT, and 0.62 ± 0.28 mm in MR. The printed material had T2 ∼32 ms, T2*∼7 ms, T1 ∼193 ms, and a very small diffusion coefficient less than olive oil. MRI monitoring of the cryoablation demonstrated iceball formation similar to an in vivo procedure. CONCLUSION: Current 3D printing technology can be used to print anatomically accurate phantoms that can be imaged by both CT and MRI. Such models can be used to simulate MRI-guided interventions such as cryosurgeries. Future development of the proposed technique can potentially lead to printed models that depict different tissues and anatomical structures with different MR signal characteristics. Magn Reson Med 77:613-622, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Lightweight preprocessing and template matching facilitate streamlined ischemic myocardial scar classification.

Purpose: Ischemic myocardial scarring (IMS) is a common outcome of coronary artery disease that potentially leads to lethal arrythmias and heart failure. Late-gadolinium-enhanced cardiac magnetic resonance (CMR) imaging scans have served as the diagnostic bedrock for IMS, with recent advancements in machine learning enabling enhanced scar classification. However, the trade-off for these improvements is intensive computational and time demands. As a solution, we propose a combination of lightweight preprocessing (LWP) and template matching (TM) to streamline IMS classification. Approach: CMR images from 279 patients (151 IMS, 128 control) were classified for IMS presence using two convolutional neural networks (CNNs) and TM, both with and without LWP. Evaluation metrics included accuracy, sensitivity, specificity, F1-score, area under the receiver operating characteristic curve (AUROC), and processing time. External testing dataset analysis encompassed patient-level classifications (PLCs) and a CNN versus TM classification comparison (CVTCC). Results: LWP enhanced the speed of both CNNs (4.9x) and TM (21.9x). Furthermore, in the absence of LWP, TM outpaced CNNs by over 10x, while with LWP, TM was more than 100x faster. Additionally, TM performed similarly to the CNNs in accuracy, sensitivity, specificity, F1-score, and AUROC, with PLCs demonstrating improvements across all five metrics. Moreover, the CVTCC revealed a substantial 90.9% agreement. Conclusions: Our results highlight the effectiveness of LWP and TM in streamlining IMS classification. Anticipated enhancements to LWP's region of interest (ROI) isolation and TM's ROI targeting are expected to boost accuracy, positioning them as a potential alternative to CNNs for IMS classification, supporting the need for further research.

Enhancing cerebral vasculature analysis with pathlength-corrected 2D angiographic parametric imaging: A feasibility study.

BACKGROUND: 2D angiographic parametric imaging (API) quantitatively extracts imaging biomarkers related to contrast flow and is conventionally applied to 2D digitally subtracted angiograms (DSA's). In the interventional suite, API is typically performed using 1-2 projection views and is limited by vessel overlap, foreshortening, and depth-integration of contrast motion. PURPOSE: This work explores the use of a pathlength-correction metric to overcome the limitations of 2D-API: the primary objective was to study the effect of converting 3D contrast flow to projected contrast flow using a simulated angiographic framework created with computational fluid dynamics (CFD) simulations, thereby removing acquisition variability. METHODS: The pathlength-correction framework was applied to in-silico angiograms, generating a reference (i.e., ground-truth) volumetric contrast distribution in four patient-specific intracranial aneurysm geometries. Biplane projections of contrast flow were created from the reference volumetric contrast distributions, assuming a cone-beam geometry. A Parker-weighted reconstruction was performed to obtain a binary representation of the vessel structure in 3D. Standard ray tracing techniques were then used to track the intersection of a ray from the focal spot with each voxel of the reconstructed vessel wall to a pixel in the detector plane. The lengths of each ray through the 3D vessel lumen were then projected along each ray-path to create a pathlength-correction map, where the pixel intensity in the detector plane corresponds to the vessel width along each source-detector ray. By dividing the projection sequences with this correction map, 2D pathlength-corrected in-silico angiograms were obtained. We then performed voxel-wise (3D) API on the ground-truth contrast distribution and compared it to pixel-wise (2D) API, both with and without pathlength correction for each biplane view. The percentage difference (PD) between the resultant API biomarkers in each dataset were calculated within the aneurysm region of interest (ROI). RESULTS: Intensity-based API parameters, such as the area under the curve (AUC) and peak height (PH), exhibited notable changes in magnitude and spatial distribution following pathlength correction: these now accurately represent conservation of mass of injected contrast media within each arterial geometry and accurately reflect regions of stagnation and recirculation in each aneurysm ROI. Improved agreement was observed between these biomarkers in the pathlength-corrected biplane maps: the maximum PD within the aneurysm ROI is 3.3% with pathlength correction and 47.7% without pathlength correction. As expected, improved agreement with ROI-averaged ground-truth 3D counterparts was observed for all aneurysm geometries, particularly large aneurysms: the maximum PD for both AUC and PH was 5.8%. Temporal parameters (mean transit time, MTT, time-to-peak, TTP, time-to-arrival, TTA) remained unaffected after pathlength correction. CONCLUSIONS: This study indicates that the values of intensity-based API parameters obtained with conventional 2D-API, without pathlength correction, are highly dependent on the projection orientation, and uncorrected API should be avoided for hemodynamic analysis. The proposed metric can standardize 2D API-derived biomarkers independent of projection orientation, potentially improving the diagnostic value of all acquired 2D-DSA's. Integration of a pathlength correction map into the imaging process can allow for improved interpretation of biomarkers in 2D space, which may lead to improved diagnostic accuracy during procedures involving the cerebral vasculature.

Effect of singular value decomposition on removing injection variability in 2D quantitative angiography: An in silico and in vitro phantoms study.

BACKGROUND: Intraoperative 2D quantitative angiography (QA) for intracranial aneurysms (IAs) has accuracy challenges due to the variability of hand injections. Despite the success of singular value decomposition (SVD) algorithms in reducing biases in computed tomography perfusion (CTP), their application in 2D QA has not been extensively explored. This study seeks to bridge this gap by investigating the potential of SVD-based deconvolution methods in 2D QA, particularly in addressing the variability of injection durations. PURPOSE: Building on the identified limitations in QA, the study aims to adapt SVD-based deconvolution techniques from CTP to QA for IAs. This adaptation seeks to capitalize on the high temporal resolution of QA, despite its two-dimensional nature, to enhance the consistency and accuracy of hemodynamic parameter assessment. The goal is to develop a method that can reliably assess hemodynamic conditions in IAs, independent of injection variables, for improved neurovascular diagnostics. MATERIALS AND METHODS: The study included three internal carotid aneurysm (ICA) cases. Virtual angiograms were generated using computational fluid dynamics (CFD) for three physiologically relevant inlet velocities to simulate contrast media injection durations. Time-density curves (TDCs) were produced for both the inlet and aneurysm dome. Various SVD variants, including standard SVD (sSVD) with and without classical Tikhonov regularization, block-circulant SVD (bSVD), and oscillation index SVD (oSVD), were applied to virtual angiograms. The method was applied on virtual angiograms to recover the aneurysmal dome impulse response function (IRF) and extract flow related parameters such as Peak Height PHIRF, Area Under the Curve AUCIRF, and Mean transit time MTT. Next, correlations between QA parameters, injection duration, and inlet velocity were assessed for unconvolved and deconvolved data for all SVD methods. Additionally, we performed an in vitro study, to complement our in silico investigation. We generated a 2D DSA using a flow circuit design for a patient-specific internal carotid artery phantom. The DSA showcases factors like x-ray artifacts, noise, and patient motion. We evaluated QA parameters for the in vitro phantoms using different SVD variants and established correlations between QA parameters, injection duration, and velocity for unconvolved and deconvolved data. RESULTS: The different SVD algorithm variants showed strong correlations between flow and deconvolution-adjusted QA parameters. Furthermore, we found that SVD can effectively reduce QA parameter variability across various injection durations, enhancing the potential of QA analysis parameters in neurovascular disease diagnosis and treatment. CONCLUSION: Implementing SVD-based deconvolution techniques in QA analysis can enhance the precision and reliability of neurovascular diagnostics by effectively reducing the impact of injection duration on hemodynamic parameters.

Clinical situations for which 3D printing is considered an appropriate representation or extension of data contained in a medical imaging examination: pediatric congenital heart disease conditions.

BACKGROUND: The use of medical 3D printing (focusing on anatomical modeling) has continued to grow since the Radiological Society of North America's (RSNA) 3D Printing Special Interest Group (3DPSIG) released its initial guideline and appropriateness rating document in 2018. The 3DPSIG formed a focused writing group to provide updated appropriateness ratings for 3D printing anatomical models across a variety of congenital heart disease. Evidence-based- (where available) and expert-consensus-driven appropriateness ratings are provided for twenty-eight congenital heart lesion categories. METHODS: A structured literature search was conducted to identify all relevant articles using 3D printing technology associated with pediatric congenital heart disease indications. Each study was vetted by the authors and strength of evidence was assessed according to published appropriateness ratings. RESULTS: Evidence-based recommendations for when 3D printing is appropriate are provided for pediatric congenital heart lesions. Recommendations are provided in accordance with strength of evidence of publications corresponding to each cardiac clinical scenario combined with expert opinion from members of the 3DPSIG. CONCLUSIONS: This consensus appropriateness ratings document, created by the members of the RSNA 3DPSIG, provides a reference for clinical standards of 3D printing for pediatric congenital heart disease clinical scenarios.

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

Purpose: Physics-informed neural networks (PINNs) and computational fluid dynamics (CFD) have both demonstrated an ability to derive accurate hemodynamics if boundary conditions (BCs) are known. Unfortunately, patient-specific BCs are often unknown, and assumptions based upon previous investigations are used instead. High speed angiography (HSA) may allow extraction of these BCs due to the high temporal fidelity of the modality. We propose to investigate whether PINNs using convection and Navier-Stokes equations with BCs derived from HSA data may allow for extraction of accurate hemodynamics in the vasculature. Materials and Methods: Imaging data generated from in vitro 1000 fps HSA, as well as simulated 1000 fps angiograms generated using CFD were utilized for this study. Calculations were performed on a 3D lattice comprised of 2D projections temporally stacked over the angiographic sequence. A PINN based on an objective function comprised of the Navier-Stokes equation, the convection equation, and angiography-based BCs was used for estimation of velocity, pressure and contrast flow at every point in the lattice. Results: Imaging-based PINNs show an ability to capture such hemodynamic phenomena as vortices in aneurysms and regions of rapid transience, such as outlet vessel blood flow within a carotid artery bifurcation phantom. These networks work best with small solution spaces and high temporal resolution of the input angiographic data, meaning HSA image sequences represent an ideal medium for such solution spaces. Conclusions: The study shows the feasibility of obtaining patient-specific velocity and pressure fields using an assumption-free data driven approach based purely on governing physical equations and imaging data.

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

Purpose: Contrast dilution gradient (CDG) analysis is a quantitative method allowing blood velocity estimation using angiographic acquisitions. Currently, CDG is restricted to peripheral vasculature due to the suboptimal temporal resolution of current imaging systems. We investigate extension of CDG methods to the flow conditions of proximal vasculature using 1000 frames per second (fps) high-speed angiographic (HSA) imaging. Approach: We performed in-vitro HSA acquisitions using the XC-Actaeon detector and 3D-printed patient-specific phantoms. The CDG approach was used for blood velocity estimation expressed as the ratio of temporal and spatial contrast gradients. The gradients were extracted from 2D contrast intensity maps synthesized by plotting intensity profiles along the arterial centerline at each frame. In-vitro results obtained at various frame rates via temporal binning of 1000 fps data were retrospectively compared to computational fluid dynamics (CFD) velocimetry. Full-vessel velocity distributions were estimated at 1000 fps via parallel line expansion of the arterial centerline analysis. Results: Using HSA, the CDG method displayed agreement with CFD at or above 250 fps [mean-absolute error (MAE): 2.6±6.3 cm/s, p=0.05]. Relative velocity distributions correlated well with CFD at 1000 fps with universal underapproximation due to effects of pulsatile contrast injection (MAE: 4.3 cm/s). Conclusions: Using 1000 fps HSA, CDG-based extraction of velocities across large arteries is possible. The method is sensitive to noise; however, image processing techniques and a contrast injection, which adequately fills the vessel assist algorithm accuracy. The CDG method provides high resolution quantitative information for rapidly transient flow patterns observed in arterial circulation.

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.

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

Purpose: Previous studies have demonstrated the efficacy of contrast dilution gradient (CDG) analysis in determining large vessel velocity distributions from 1000 fps high-speed angiography (HSA). However, the method required vessel centerline extraction, which made it applicable only to non-tortuous geometries using a highly specific contrast injection technique. This study seeks to remove the need for a priori knowledge regarding the direction of flow and modify the vessel sampling method to make the algorithm more robust to non-linear geometries. Materials and Methods: 1000 fps HSA acquisitions were obtained in vitro with a benchtop flow loop using the XC-Actaeon (Varex Inc.) photon-counting detector, and in silico using a passive-scalar transport model within a computational fluid dynamics (CFD) simulation. CDG analyses were obtained using gridline sampling across the vessel, and subsequent 1D velocity measurement in both the x- and y-directions. The velocity magnitudes derived from the component CDG velocity vectors were aligned with CFD results via co-registration of the resulting velocity maps and compared using mean absolute percent error (MAPE) between pixels values in each method after temporal averaging of the 1-ms velocity distributions. Results: Regions well-saturated with contrast throughout the acquisition showed agreement when compared to CFD (MAPE of 18% for the carotid bifurcation inlet and MAPE of 27% for the internal carotid aneurysm), with respective completion times of 137 seconds and 5.8 seconds. Conclusions: CDG may be used to obtain velocity distributions in and surrounding vascular pathologies provided the contrast injection is sufficient to provide a gradient, and diffusion of contrast through the system is negligible.

Clinical situations for which 3D Printing is considered an appropriate representation or extension of data contained in a medical imaging examination: vascular conditions.

BACKGROUND: Medical three-dimensional (3D) printing has demonstrated utility and value in anatomic models for vascular conditions. A writing group composed of the Radiological Society of North America (RSNA) Special Interest Group on 3D Printing (3DPSIG) provides appropriateness recommendations for vascular 3D printing indications. METHODS: A structured literature search was conducted to identify all relevant articles using 3D printing technology associated with vascular indications. Each study was vetted by the authors and strength of evidence was assessed according to published appropriateness ratings. RESULTS: Evidence-based recommendations for when 3D printing is appropriate are provided for the following areas: aneurysm, dissection, extremity vascular disease, other arterial diseases, acute venous thromboembolic disease, venous disorders, lymphedema, congenital vascular malformations, vascular trauma, vascular tumors, visceral vasculature for surgical planning, dialysis access, vascular research/development and modeling, and other vasculopathy. Recommendations are provided in accordance with strength of evidence of publications corresponding to each vascular condition combined with expert opinion from members of the 3DPSIG. CONCLUSION: This consensus appropriateness ratings document, created by the members of the 3DPSIG, provides an updated reference for clinical standards of 3D printing for the care of patients with vascular conditions.

Improving the accuracy of computational fluid dynamics simulations of coiled cerebral aneurysms using finite element modeling.

Cerebral aneurysms are a serious clinical challenge, with ∼half resulting in death or disability. Treatment via endovascular coiling significantly reduces the chances of rupture, but the techniquehas failure rates of ∼20 %. This presents a pressing need to develop a method fordetermining optimal coildeploymentstrategies. Quantification of the hemodynamics of coiled aneurysms using computational fluid dynamics (CFD) has the potential to predict post-treatment outcomes, but representing the coil mass in CFD simulations remains a challenge. We use the Finite Element Method (FEM) for simulating patient-specific coil deployment for n = 4 ICA aneurysms for which 3D printed in vitro models were also generated, coiled, and scanned using ultra-high resolution synchrotron micro-CT. The physical and virtual coil geometries were voxelized onto a binary structured grid and porosity maps were generated for geometric comparison. The average binary accuracy score is 0.8623 and the average error in porosity map is 4.94 %. We then conduct patient-specific CFD simulations of the aneurysm hemodynamics using virtual coils geometries, micro-CT generated oil geometries, and using the porous medium method to represent the coil mass. Hemodynamic parameters including Neck Inflow Rate (Qneck) and Wall Shear Stress (WSS) were calculated for each of the CFD simulations. The average relative error in Qneck and WSS from CFD using FEM geometry were 6.6 % and 21.8 % respectively, while the error from CFD using a porous media approximation resulted in errors of 55.1 % and 36.3 % respectively; demonstrating a marked improvement in the accuracy of CFD simulations using FEM generated coil geometries.

Identification of infarct core and ischemic penumbra using computed tomography perfusion and deep learning.

Purpose: The size and location of infarct and penumbra are key to decision-making for acute ischemic stroke (AIS) management. CT perfusion (CTP) software estimate infarct and penumbra volume using contralateral hemisphere relative thresholding. This approach is not robust and widely contested by the scientific community. In this study, we investigate the use of deep learning-based algorithms to efficiently locate infarct and penumbra tissue on CTP hemodynamic maps. Approach: CTP scans were retrospectively collected for 60 and 59 patients in the infarct only and infarct + penumbra substudies respectively. Commercial CTP software was used to generate cerebral blood flow, cerebral blood volume, mean transit time, time to peak, and delay time maps. U-Net-shaped architectures were trained to segment infarct or infarct + penumbra. Test-time-augmentation, ensembling, and watershed segmentation were used as postprocessing techniques. Segmentation performance was evaluated using Dice coefficients (DC) and mean absolute volume errors (MAVE). Results: The algorithm segmented infarct tissue resulted in DC of 0.64 ± 0.03 (0.63, 0.65), and MAVE of 4.91 ± 0.94 (4.5, 5.32) mL. In comparison, the commercial software predicted infarct with a DC of 0.31 ± 0.17 (0.26, 0.36) and MAVE of 9.77 ± 8.35 (7.12, 12.42) mL. The algorithm was able to segment infarct + penumbra with a DC of 0.61 ± 0.04 (0.6, 0.63), and MAVE of 6.51 ± 1.37 (5.91, 7.11) mL. In comparison, the commercial software predicted infarct + penumbra with a DC of 0.3 ± 0.19 (0.25, 0.35) and MAVE of 9.18 ± 7.55 (7.25, 11.11) mL. Conclusions: Use of deep learning algorithms to assess severity of AIS in terms of infarct and penumbra volume is precise and outperforms current relative thresholding methods. Such an algorithm would enhance the reliability of CTP in guiding treatment decisions.

Radiological Society of North America (RSNA) 3D Printing Special Interest Group (SIG) clinical situations for which 3D printing is considered an appropriate representation or extension of data contained in a medical imaging examination: breast conditions.

The use of medical 3D printing has expanded dramatically for breast diseases. A writing group composed of the Radiological Society of North America (RSNA) Special Interest Group on 3D Printing (SIG) provides updated appropriateness criteria for breast 3D printing in various clinical scenarios. Evidence-based appropriateness criteria are provided for the following clinical scenarios: benign breast lesions and high-risk breast lesions, breast cancer, breast reconstruction, and breast radiation (treatment planning and radiation delivery).

Recovery of complete time density curves from incomplete angiographic data using recurrent neural networks.

Quantitative angiography is a 2D/3D x-ray imaging modality that summarizes hemodynamic information using time density curve (TDC) based parameters. Estimation of the TDC parameters are susceptible to errors due to various factors including, patient motion, incomplete temporal data, imaging trigger errors etc. In this study, we tested the feasibility of using recurrent neural networks (RNN) to recover complete TDC temporal information from incomplete sequences and evaluate quantitative parameters generated from the corrected TDCs. Digital subtraction angiograms (DSAs) were collected from patients undergoing endovascular treatments and angiographic parametric imaging (API) parameters were calculated from each DSA. Each set of API parameters was used to simulate a TDC resulting in a dataset of 760 TDCs. One-third of each TDC was continuously masked from pseudo-random points past the peak height (PH) point to simulate missing/artifact information. An RNN was developed, trained and tested to generate completed/corrected TDCs. The RNN recovered complete TDC temporal information with an average mean squared error of 0.0086±0.002. Average mean absolute errors were calculated between each API parameter generated from the ground truth TDCs and RNN corrected TDCs, these were 11.02%±0.91 for time to peak, 10.97%±0.69 for mean transit time, 5.65%±0.76 for PH, and 15.08%±0.98 for area under the TDC. The change in API parameters was not clinically significant and the predictive power of the API parameters was retained. This study proved the feasibility of using RNNs to mitigate motion artifacts and incomplete angiographic acquisitions to extract accurate quantitative parameters.

Automation of ischemic myocardial scar detection in cardiac magnetic resonance imaging of the left ventricle using machine learning.

Purpose: Machine learning techniques can be applied to cardiac magnetic resonance imaging (CMR) scans in order to differentiate patients with and without ischemic myocardial scarring (IMS). However, processing the image data in the CMR scans requires manual work that takes a significant amount of time and expertise. We propose to develop and test an AI method to automatically identify IMS in CMR scans to streamline processing and reduce time costs. Materials and Methods: CMR scans from 170 patients (138 IMS & 32 without IMS as identified by a clinical expert) were processed using a multistep automatic image data selection algorithm. This algorithm consisted of cropping, circle detection, and supervised machine learning to isolate focused left ventricle image data. We used a ResNet-50 convolutional neural network to evaluate manual vs. automatic selection of left ventricle image data through calculating accuracy, sensitivity, specificity, F1 score, and area under the receiver operating characteristic curve (AUROC). Results: The algorithm accuracy, sensitivity, specificity, F1 score, and AUROC were 80.6%, 85.6%, 73.7%, 83.0%, and 0.837, respectively, when identifying IMS using manually selected left ventricle image data. With automatic selection of left ventricle image data, the same parameters were 78.5%, 86.0%, 70.7%, 79.7%, and 0.848, respectively. Conclusion: Our proposed automatic image data selection algorithm provides a promising alternative to manual selection when there are time and expertise limitations. Automatic image data selection may also prove to be an important and necessary step toward integration of machine learning diagnosis and prognosis in clinical workflows.

Comparison of fluid dynamics changes due to physical activity in 3D printed patient specific coronary phantoms with the Windkessel equivalent model of coronary flow.

BACKGROUND: 3D printing (3DP) used to replicate the geometry of normal and abnormal vascular pathologies has been demonstrated in many publications; however, reproduction of hemodynamic changes due to physical activities, such as rest versus moderate exercise, need to be investigated. We developed a new design for patient specific coronary phantoms, which allow adjustable physiological variables such as coronary distal resistance and coronary compliance in patients with coronary artery disease. The new design was tested in precise benchtop experiments and compared with a theoretical Windkessel electrical circuit equivalent, that models coronary flow and pressure using arterial resistance and compliance. METHODS: Five phantoms from patients who underwent clinically indicated elective invasive coronary angiography were built from CCTA scans using multi-material 3D printing. Each phantom was used in a controlled flow system where patient specific flow conditions were simulated by a programmable cardiac pump. To simulate the arteriole and capillary beds flow resistance and the compliance for various physical activities, we designed a three-chamber outlet system which controls the outflow dynamics of each coronary tree. Benchtop pressure measurements were recorded using sensors embedded in each of the main coronary arteries. Using the Windkessel model, patient specific flow equivalent electrical circuit models were designed for each coronary tree branch, and flow in each artery was determined for known inflow conditions. Local flow resistances were calculated through Poiseuille's Law derived from the radii and lengths of the coronary arteries using CT angiography based multi-planar reconstructions. The coronary stenosis flow rates from the benchtop and the electrical models were compared to the localized flow rates calculated from invasive pressure measurements recorded in the angio-suites. RESULTS: The average Pearson correlations of the localized flow rates at the location of the stenosis between each of the models (Benchtop/Electrical, Benchtop/Angio, Electrical/Angio) are 0.970, 0.981, and 0.958 respectively. CONCLUSIONS: 3D printed coronary phantoms can be used to replicate the human arterial anatomy as well as blood flow conditions. It displays high levels of correlation when compared to hemodynamics calculated in electrically-equivalent coronary Windkessel models as well as invasive angio-suite pressure measurements.