Implementation of computer simulation to assess flow diversion treatment outcomes: systematic review and meta-analysis

Patient-specific cerebral 3D vessel model reconstruction using deep learning

Three-dimensional vessel model reconstruction from patient-specific magnetic resonance angiography (MRA) images often requires some manual maneuvers. This study aimed to establish the deep learning (DL)-based method for vessel model reconstruction. Time of flight MRA of 40 patients with internal carotid artery aneurysms was prepared, and three-dimensional vessel models were constructed using the threshold and region-growing method. Using those datasets, supervised deep learning using 2D U-net was performed to reconstruct 3D vessel models. The accuracy of the DL-based vessel segmentations was assessed using 20 MRA images outside the training dataset. The dice coefficient was used as the indicator of the model accuracy, and the blood flow simulation was performed using the DL-based vessel model. The created DL model could successfully reconstruct a three-dimensional model in all 60 cases. The dice coefficient in the test dataset was 0.859. Of note, the DL-generated model proved its efficacy even for large aneurysms (> 10 mm in their diameter). The reconstructed model was feasible in performing blood flow simulation to assist clinical decision-making. Our DL-based method could successfully reconstruct a three-dimensional vessel model with moderate accuracy. Future studies are warranted to exhibit that DL-based technology can promote medical image processing.

Quantitative analysis of hemodynamic changes induced by the discrepancy between the sizes of the flow diverter and parent artery

The efficacy of flow diverters is influenced by the strut configuration changes resulting from size discrepancies between the stent and the parent artery. This study aimed to quantitatively analyze the impact of size discrepancies between flow diverters and parent arteries on the flow diversion effects, using computational fluid dynamics. Four silicone models with varying parent artery sizes were developed. Real flow diverters were deployed in these models to assess stent configurations at the aneurysm neck. Virtual stents were generated based on these configurations for computational fluid dynamics analysis. The changes in the reduction rate of the hemodynamic parameters were quantified to evaluate the flow diversion effect. Implanting 4.0 mm flow diverters in aneurysm models with parent artery diameters of 3.0-4.5 mm, in 0.5 mm increments, revealed that a shift from oversized to undersized flow diverters led to an increase in the reduction rates of hemodynamic parameter, accompanied by enhanced metal coverage rate and pore density. However, the flow diversion effect observed transitioning from oversizing to matching was less pronounced when moving from matching to undersizing. This emphasizes the importance of proper sizing of flow diverters, considering the benefits of undersizing and not to exceed the threshold of advantages.

Experimental study using phantom models of cerebral aneurysms and 4D-DSA to measure blood flow on 3D-color-coded images

BACKGROUND

The current 3D-iFlow application can only measure the arrival time of contrast media through intensity values. If the flow rate could be estimated by 3D-iFlow, patient-specific hemodynamics could be determined within the scope of normal diagnostic management, eliminating the need for additional resources for blood flow rate estimation.

OBJECTIVE

The aim of this study is to develop and validate a method for measuring the flow rate by data obtained from 3D-iFlow images - a prototype application in Four-dimensional digital subtraction angiography (4D-DSA).

METHODS

Using phantom model and experimental circuit with circulating glycerin solution, an equation for the relationship between contrast media intensity and flow rate was developed. Applying the equation to the aneurysm phantom models, the derived flow rate was evaluated.

RESULTS

The average errors between the derived flow rate and setting flow rate became larger when the glycerin flow and the X-rays from the X-ray tube of the angiography system were parallel to each other or when the measurement point included overlaps with other contrast enhanced areas.

CONCLUSION

Although the error increases dependent on the imaging direction and overlap of contrast enhanced area, the developed equation can estimate the flow rate using the image intensity value measured on 3D-iFlow based on 4D-DSA.

Assessing a heterogeneous model for accounting for endovascular devices in hemodynamic simulations of cerebral aneurysms

The heterogeneous model developed by Berod et al [Int J Numer Method Biomed Eng 38, 2021] for representing the hemodynamic effects of endovascular prostheses is applied to a series of 10 patient specific cerebral aneurysms, 6 being treated by flow diverters, 4 being equipped with WEBs. Two markers correlated with the medical outcome of the treatment are used to assess the potential of the model, namely the saccular mean velocity and the inflow rate at the neck of the aneurysm. The comparison with the corresponding wire-resolved simulations is very favorable in both cases, and the model-based simulations also retrieve the jetting-type flows generated downstream of the struts. Noteworthy, the very same model was used for representing the flow diverters and the WEBs, showing the versatility and robustness of the heterogeneous modeling of the devices.

Effects of helical centerline stent vs. straight stent placement on blood flow velocity

As an approach to maintain patency in femoropopliteal stenting, a helical stent configuration was proposed, which showed improved patency in clinical trials. However, the effects of helical stent placement on the flow have not been quantitatively analyzed. The purpose of this study was to estimate flow velocities to quantify the influence of helical stent placement. Helical and straight stents were implanted in three healthy pigs, and the flow velocities were estimated using the time-intensity curve (TIC) in the angiography images. The angiographic images indicated thinning of the leading edge of the contrast medium through the helically deformed artery, which was not observed in the straight stent. The slower rise of the TIC peak in the helical stent indicated faster travel of this thinner edge. Arterial expansion due to stenting was observed in all cases, and the expansion rate varied according to location. All cases of helical stent implantation showed that velocity was maintained (55.0%-71.3% velocity retention), unlike for straight stent implantation (43.0%-68.0% velocity retention); however, no significant difference was observed.

Utility of virtual stenting in treatment of cerebral aneurysms by flow diverter devices

Successful endovascular treatment by stenting of intracranial aneurysms requires proper placement of the device and appropriate choice of its diameter and length. To date, several methods have been employed to achieve these goals, although each has inherent critical issues. Recently developed stent planning software applications can be used to assist interventional neuroradiologists. Based on a 3D-DSA image acquired before stenting, these applications simulate and visualize the final placement of the deployed stent. In this single-centre retrospective study, 27 patients undergoing an intravascular procedure for the treatment of intracranial aneurysms from June 2019 to July 2020 were evaluated according to strict inclusion criteria. Stent virtualization was performed with Syngo 3D Aneurysm Guidance Neuro software. We compared the software-generated stent measurement and measurements taken by the interventional radiologist. Statistical analysis was performed using the STAC web platform. Mean and standard deviations of absolute and relative discrepancies between predicted and implanted stents were recorded. Friedman's nonparametric test was used to refute the null hypotheses, i.e. (I) discrepancies between the size of virtual and implanted stents would occur, and (II) operator influence does not affect the outcome of the virtual stenting process. Based on these observations, it is believed that the virtual stenting process can validly assist interventional neuroradiologists in selecting the appropriate device and reducing peri- and post-procedural complications. The results of our study suggest that virtual reality simulation of devices used for endovascular treatment of intracranial aneurysms is a useful, rapid, and accurate tool for interventional procedure planning.

Evaluation of aneurysm rupture risk based upon flowrate-independent hemodynamic parameters: a multi-center pilot study

BACKGROUND

Specifying generic flow boundary conditions in aneurysm hemodynamic simulations yields a great degree of uncertainty for the evaluation of aneurysm rupture risk. Herein, we proposed the use of flowrate-independent parameters in discriminating unstable aneurysms and compared their prognostic performance against that of conventional absolute parameters.

METHODS

This retrospective study included 186 aneurysms collected from three international centers, with the stable aneurysms having a minimum follow-up period of 24 months. The flowrate-independent aneurysmal wall shear stress (WSS) and energy loss (EL) were defined as the coefficients of the second-order polynomials characterizing the relationships between the respective parameters and the parent-artery flows. Performance of the flowrate-independent parameters in discriminating unstable aneurysms with the logistic regression, Adaboost, and support-vector machine (SVM) methods was quantified and compared against that of the conventional parameters, in terms of sensitivity, specificity, and area under the curve (AUC).

RESULTS

In discriminating unstable aneurysms, the proposed flowrate-independent EL achieved the highest sensitivity (0.833, 95% CI 0.586 to 0.964) and specificity (0.833, 95% CI 0.672 to 0.936) on the SVM, with the AUC outperforming the conventional EL by 0.133 (95% CI 0.039 to 0.226, p=0.006). Likewise, the flowrate-independent WSS outperformed the conventional WSS in terms of the AUC (difference: 0.137, 95% CI 0.033 to 0.241, p=0.010).

CONCLUSION

The flowrate-independent hemodynamic parameters surpassed their conventional counterparts in predicting the stability of aneurysms, which may serve as a promising set of hemodynamic metrics to be used for the prediction of aneurysm rupture risk when physiologically real vascular boundary conditions are unavailable.

Critical Angiographic and Sonographic Analysis of Intra Aneurysmal and Downstream Hemodynamic Changes After Flow Diversion

Introduction

Successful treatment of intracranial aneurysms after flow diversion (FD) is dependent on the flow modulating effect of the device. We aimed to investigate the intra-aneurysmal and parent vessel hemodynamic changes, as well as the incidence of silent emboli following treatment with various FD devices.

Methods

We evaluated the appearance of the eclipse sign in nine distinct phases of cerebral angiography before and immediately after FD placement in correlation with aneurysm occlusion. Angiographic and clinical data of consecutive procedures were analyzed retrospectively. Patients who had successful FD procedure without adjunctive coiling, visible eclipse sign on post embolization angiography, and reliable follow-up angiographic data were included in the analysis. Detailed analysis of hemodynamic data from transcranial doppler after FD was performed in selected patients, such as monitoring for silent emboli.

Results

Among all patients (N = 65) who met inclusion criteria, complete aneurysm occlusion at 12 months was achieved in 89% (58/65). Eclipse sign prior to FD was observed in 42% (27/65) with unchanged appearance in 4.6% (3/65) of the treated patients. None of these three patients achieved complete aneurysm occlusion. Among all analyzed variables, such as aneurysm size, device type used, age, and appearance of the eclipse sign pre- and post-FD, the most reliable predictor of permanent aneurysm occlusion at 12 months was earlier, prolonged, and sustained eclipse sign visibility in more than three angiographic phases in comparison to the baseline (p < 0.001). Elevation in flow velocities within the ipsilateral vascular territory was noted in 70% (9/13), and bilaterally in 54% (7/13) of the treated patients. None of the patients had silent emboli.

Conclusions

Intra-aneurysmal and parent vessel hemodynamic changes after FD can be reliably assessed by the cerebral angiography and transcranial doppler with important implications for the prediction of successful treatment.

Incomplete stent expansion in flow-diversion treatment affects aneurysmal haemodynamics: a quantitative comparison of treatments affected by different severities of malapposition occurring in different segments of the parent artery

Incomplete stent expansion (IncSE) is occasionally seen in flow-diversion (FD) treatment of intracranial aneurysms; however, its haemodynamic consequences remain inconclusive. Through a parametric study, we quantify the aneurysmal haemodynamics subject to different severities of IncSE occurring in different portions of the stent. Two patient cases with IncSE confirmed in vivo were studied. To investigate a wider variety of IncSE scenarios, we modelled IncSE at two severity levels respectively located in the proximal, central, or distal segment of a stent, yielding a total of 14 treatment scenarios (including the ideal deployment). We examined stent wire configurations in 14 scenarios and resolved aneurysm haemodynamics through computational fluid dynamics (CFD). A considerable degradation of aneurysm flow-reduction performance was observed when central or distal IncSE occurred, with the maximal elevations of the inflow rate (IR) and energy loss (EL) being 10% and 15%. The underlying mechanism might be the increased resistance for flow to remain within the FD stent, which forces more blood to leak into the aneurysm sac. Counter-intuitively, a slight reduction of aneurysm inflow was associated with proximal IncSE, with the maximal further reduction of the IR and EL being 5% and 8%. This may be due to the disruption of the predominant parent-artery flow by the collapsed wires, which decreased the strength and altered the direction of aneurysmal inflow. The effects of IncSE vary greatly with the location of occurrence, revealing the importance of performing individualised, patient-specific risk assessment before treatment.

VICTORIA: VIrtual neck Curve and True Ostium Reconstruction of Intracranial Aneurysms

Purpose

For the status evaluation of intracranial aneurysms (IAs), morphological and hemodynamic parameters can provide valuable information. For their extraction, a separation of the aneurysm sac from its parent vessel is required that yields the neck curve and the ostium. However, manual and subjective neck curve and ostium definitions might lead to inaccurate IA assessments.

Methods

The research project VICTORIA was initiated, allowing users to interactively define the neck curve of five segmented IA models using a web application. The submitted results were qualitatively and quantitatively compared to identify the minimum, median and maximum aneurysm surface area. Finally, image-based blood flow simulations were carried out to assess the effect of variable neck curve definitions on relevant flow- and shear-related parameters.

Results

In total, 55 participants (20 physicians) from 18 countries participated in VICTORIA. For relatively simple aneurysms, a good agreement with respect to the neck curve definition was found. However, differences among the participants increased with increasing complexity of the aneurysm. Furthermore, it was observed that the majority of participants excluded any small arteries occurring in the vicinity of an aneurysm. This can lead to non-negligible deviations among the flow- and shear-related parameters, which need to be carefully evaluated, if quantitative analysis is desired. Finally, no differences between participants with medical and non-medical background could be observed.

Conclusions

VICTORIAs findings reveal the complexity of aneurysm neck curve definition, especially for bifurcation aneurysms. Standardization appears to be mandatory for future sac-vessel-separations. For hemodynamic simulations a careful neck curve definition is crucial to avoid inaccuracies during the quantitative flow analysis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13239-021-00535-w.

Association Between Aneurysmal Haemodynamics and Device Microstructural Characteristics After Flow-Diversion Treatments With Dual Stents of Different Sizes: A Numerical Study

Objectives

Treating intracranial aneurysms with flow-diverting stents sometimes requires deployment of a second device. Herein we quantify the sizing effects of devices in dual-stent treatments upon the final stent microstructure and the post-treatment aneurysmal haemodynamics.

Methods

Fifteen sidewall ICA aneurysm geometries were included. Using a virtual stenting technique, we implanted either one or two stents for each aneurysm treatment considered, with each stent specified as one of two different sizes, yielding a total of two single-stent and fouir dual-stent treatment scenarios for each aneurysm. Three stent microstructural parameters and nine aneurysmal haemodynamic parameters were quantified and systematically compared across the 90 treatment scenarios.

Results

Deployment of a second stent further reduced the aneurysmal inflow rate (IR) and energy loss (EL) by, respectively, 14 ± 11% (p = 0.001) and 9 ± 12% (p = 0.056), relative to the untreated condition. Sizing effects of the earlier-deployed stent led to largest differences of 6.9% for the final IR reduction and 11.1% for the EL, whereas sizing effects from the later-deployed stent were minor (≤2.1%). The change in stent pore size was the only microstructural parameter demonstrating a strong correlation with the reduction in the post-treatment aneurysmal haemodynamics, in terms of the IR (r = 0.50, p < 0.001) and pressure drop (r = 0.63, p < 0.001).

Conclusion

Size of the earlier-deployed stent has substantial effects on the final haemodynamic outcomes after dual-stent treatment. The average pore size of stent wires at the aneurysm orifice shows promise as a potential index for predicting the efficacy of flow-diversion treatments.