Skip to main content

Advertisement

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH): uncertainty quantification of geometric rupture risk parameters

Abstract

Background

Geometric parameters have been proposed for prediction of cerebral aneurysm rupture risk. Predicting the rupture risk for incidentally detected unruptured aneurysms could help clinicians in their treatment decision. However, assessment of geometric parameters depends on several factors, including the spatial resolution of the imaging modality used and the chosen reconstruction procedure. The aim of this study was to investigate the uncertainty of a variety of previously proposed geometric parameters for rupture risk assessment, caused by variability of reconstruction procedures.

Materials

26 research groups provided segmentations and surface reconstructions of five cerebral aneurysms as part of the Multiple Aneurysms AnaTomy CHallenge (MATCH) 2018. 40 dimensional and non-dimensional geometric parameters, describing aneurysm size, neck size, and irregularity of aneurysm shape, were computed. The medians as well as the absolute and relative uncertainties of the parameters were calculated. Additionally, linear regression analysis was performed on the absolute uncertainties and the median parameter values.

Results

A large variability of relative uncertainties in the range between 3.9 and 179.8% was found. Linear regression analysis indicates that some parameters capture similar geometric aspects. The lowest uncertainties < 6% were found for the non-dimensional parameters isoperimetric ratio, convexity ratio, and ellipticity index. Uncertainty of 2D and 3D size parameters was significantly higher than uncertainty of 1D parameters. The most extreme uncertainties > 80% were found for some curvature parameters.

Conclusions

Uncertainty analysis is essential on the road to clinical translation and use of rupture risk prediction models. Uncertainty quantification of geometric rupture risk parameters provided by this study may help support development of future rupture risk prediction models.

Background

Intracranial aneurysms (IAs) are complexly shaped local dilatations of the cerebral vasculature. They occur at different locations in the anterior and posterior cerebral circulation and are unique in size and phenotype [1, 2]. Unruptured IAs are found in approximately 3.2% of the adult population worldwide [3] and are being discovered with an increasing frequency due to the widespread use of high-resolution magnetic resonance imaging or computed tomography [3].

However, according to published data of two large studies (the International Study of Unruptured Intracranial Aneurysms (ISUIA) [4] and a Japanese study of unruptured aneurysms [5]), the large majority of IAs never rupture. As aneurysm treatment is associated with significant risks, physicians have to weigh the risk of treatment against the rupture risk of yet unruptured aneurysms.

Rupture risk of unruptured IAs increases with increasing aneurysmal diameter [4] and current clinical guidelines therefore recommend to treat aneurysms larger than 7 mm, it was shown that smaller aneurysms can also lead to subarachnoid hemorrhage [6, 7]. To integrate other clinical factors for the assessment of IA rupture risk, (e.g., specific population, arterial hypertension, age, sex, previous rupture of an independent intracranial aneurysm, site), the PHASES score was developed [8, 9]. However, it has been demonstrated that in practice PHASES underestimates individual rupture risk and may lead to an inappropriate treatment [10, 11].

Therefore, advanced techniques were developed to improve assessment of aneurysm morphology. Raghavan et al. [12] defined several metrics (e.g., nonsphericity index, undulation index, ellipticity index) to enable accurate description of complex aneurysm shapes. In the last decade, multiple geometric parameters were added [13,14,15,16]. Nevertheless, none of these complex quantifications are introduced into clinics yet. Physicians still prefer assessment of rupture risk using one-dimensional or normalized features such as aneurysm diameter or aspect ratio. One reason for this is, that these parameters can be measured directly from imaging data, while complex geometric risk parameters require segmentation and reconstruction of the aneurysm’s 3D geometry.

To provide a clinically applicable workflow, Saalfeld et al. [17] analysed 100 IAs by means of a semi-automatic neck curve reconstruction algorithm. The authors identified the characteristic dome point angle as a potential morphologic candidate to assess rupture risk. A very recent literature review by Liang et al. [18], included 46 studies containing 2791 aneurysms and found the simple metrics aspect ratio and aneurysm size as the most relevant parameters for rupture prediction. As these studies depend strongly on specific processing techniques, inter-study comparability might not be given concerning the lack of a gold standard for segmentation. Multiple algorithms are available for general vessel segmentation, but are not specifically tailored to the reconstruction of IAs [19]. Furthermore, custom-made smoothing algorithms, the subjective impression of the corresponding investigator, and different levels of experience might result in further inaccuracies.

Models for discriminating between ruptured and unruptured aneurysms proposed for a cohort available to one research group often fail when applied to other patient cohorts. While the segmentation procedure might be standardized fairly well in one group, other research groups might use different methods as no consensus on image segmentation of cerebral aneurysms exists. Uncertainties caused by different image segmentation methods will result in uncertainties in geometric and hemodynamic parameters used for rupture risk prediction. These uncertainties are major factors preventing successful introduction of image-based methods in clinical practice [20, 21].

Accuracy of segmentation and the impact of its uncertainties were investigated by several studies. Berg et al. investigated the impact of voxel size and reconstruction kernel on morphologic and hemodynamic parameters. They found a higher susceptibility of parameters associated with neck representation and wall shear stress [22, 23]. Klepaczko et al. found a variation of vessel size between 7 and 79% depending on the vessel size, image acquisition and segmentation method [24]. However, most studies focus on the impact of segmentation on hemodynamic parameters [25,26,27,28]. As of yet, no information is available about uncertainties of morphometric rupture risk parameters caused by inter-group variability in segmentation.

Therefore, the international Multiple Aneurysms AnaTomy CHallenge (MATCH) was announced to directly compare segmentation capabilities based on identical aneurysm datasets. Within MATCH, 26 groups from 13 countries contributed the segmentation results of five IAs with different location, size, and shape. Initial results revealed that remarkable differences with respect to vessel and aneurysm representation occurred [29].

To further address this situation, the current study is a follow-up analysis of reconstructed IAs contributed within MATCH. The aim is to quantify the segmentation-related uncertainty of multiple geometric parameters previously proposed for rupture risk prediction. This analysis of variability is intended to enable physicians and biomedical engineers to critically assess those quantities. Finally, the knowledge of segmentation uncertainty could affect the performance of rupture risk prediction models both directly and indirectly. This is because geometric parameters and hemodynamic parameters (based on computational fluid dynamics simulations), both of which are affected by geometry reconstruction, have been proposed as rupture risk predictors [3, 12, 30, 31]. Furthermore, parameters with low uncertainty can be identified for preferential use in future studies.

Methods

In order to quantify the uncertainties of a set of geometric parameters proposed in the literature as possible rupture risk predictors, reconstructions of five IAs provided in the MATCH challenge were analyzed.

Multiple Aneurysms AnaTomy CHallenge (MATCH)

The international challenge was announced on November 3rd, 2017. Interested groups were provided with anonymized clinical DICOM datasets, acquired with an interventional angio suite (Artis Q, Siemens Healthineers, Forchheim, Germany). Rotational angiographies were extracted within a 5 s run comprising 133 slices with a voxel resolution of 1240 × 960. From the rotational angiographies, a 3D digital subtraction angiography CT scan (512 × 512 × 512) was reconstructed exhibiting an isotropic voxel size of 0.28 mm. The datasets comprised three 3D digital subtraction angiography volumes containing the left anterior, right anterior, and posterior cerebral circulation. In total, the corresponding patient harbored five IAs: three in the left and right M1-segment of the middle cerebral artery (MCA), one on the left MCA bifurcation, and one on the left posterior inferior cerebellar artery.

To assess the real-world variability of existing segmentation techniques, each team was requested to apply their own algorithms and strategies. Research groups were allowed to submit their contributions until the end of January 2018 and in total 26 research teams from 13 countries participated. However, group 5 had to be excluded completely.

Overall, 18 different software tools were used for segmentation. Furthermore, four different segmentation algorithms were applied and large variation with respect to the required processing time was reported. For further details regarding the design and initial outcome of MATCH, please refer to Berg et al. [29].

Geometric rupture risk parameters

Altogether, 40 geometric parameters were calculated. These comprise twelve aneurysm size (#1–12), six neck size (#13–18), twelve non-dimensional (#19–30), as well as ten curvature parameters (#31–40), respectively. Table 1 summarizes the corresponding geometric rupture risk parameters, including a short description and references to the definitions used.

Table 1 Complete list of all geometric parameters

Note that this set of parameters does not represent all geometric parameters previously proposed in literature. For example, we excluded parameters whose definition requires manual, subjective processing steps (e.g. [32]), as well as parameters whose description in the literature did not allow us to accurately reproduce parameter calculation (e.g. [33]).

Calculation of parameters

Prior to parameter calculation, all five aneurysms were manually extracted from the segmentations provided by the 26 research groups participating in the challenge. Extraction was done using ZIBAmira (v. 2015.28, Zuse Institute Berlin, Germany). Since not every research group segmented every aneurysm, this resulted in a total of 121 aneurysm geometries. The aneurysm geometries were then checked for manifoldness, self-intersections, consistent triangle orientation, and lack of holes in the aneurysm surface and any such geometric errors were manually corrected.

Subsequently, the geometric parameters were automatically calculated for all aneurysm geometries using ZIBAmira, MATLAB (R2017b, MathWorks, Natick, USA), and Python (v. 3.7.1, Python Software Foundation, Delaware, USA) scripts.

All final aneurysm geometries as well as the table containing all parameter values for all aneurysms are provided as supplemental material, see “Availability of data and materials” section as well as Additional file 1.

Uncertainty quantification

First, median and the middle 68.3% range (i.e. the difference between the 84.13th and the 15.87th percentile) were calculated for each geometric parameter and aneurysm using SPSS (v. 23, IBM, Armonk, USA). This absolute uncertainty range corresponds to the range of mean ± one standard deviation for normally distributed data. We chose this range to allow comparison of reported uncertainty ranges against standard deviations, which are commonly reported for normally distributed data. The relative uncertainty of a parameter was then defined as the uncertainty range divided by the median value, in percent. Thus, for each parameter five relative uncertainties were calculated, one for each aneurysm.

Finally, these relative uncertainties were averaged across all five aneurysms. The average and variation of the relative uncertainty for all aneurysms were then specified using the median and the interquartile range (IQR), respectively. We decided to specify the IQR rather than the range, as the range is more sensitive to outliers. Additionally, linear regression between the absolute uncertainty range (as defined above) and the median calculated for each aneurysm was investigated.

Furthermore, the absolute relative deviations were calculated for each parameter. Here, the absolute deviation between each parameter value and the median of the parameter for the respective aneurysm was calculated and subsequently normalized using the same median value. These absolute relative deviations (121 values) were then visualized using a boxplot for easier comparison of relative parameter uncertainties (see Fig. 1).

Fig. 1
figure1

Boxplots of absolute relative deviations between measured values and the median of the respective aneurysm for all investigated geometric parameters. The plot was cropped at 100%. The number of outliers omitted due to cropping is specified on the right side of the boxplot (red numbers). For two parameters (MSD and GSD) the 75% percentile whisker is omitted as well and these values are specified in parentheses next to the plot

Statistical analysis

Statistical analysis was performed using SPSS. All data was tested for normality using the Shapiro–Wilk test. Since a large subset of parameters was not normally distributed, we used median and IQR to describe central tendency and variability, respectively. Correlation between parameters was quantified using Kendall’s tau τ. To determine whether a parameter’s uncertainty is affected by the parameter’s average value, correlations were calculated between median values and relative uncertainties calculated for each parameter and aneurysm. Furthermore, correlations were calculated between all 120 parameter values to identify similar parameters. Here, correlations of parameters describing similar geometric aspects are reported for τ > 0.7. Additionally, Wilcoxon signed-rank test was performed on the absolute relative deviations to assess differences in relative uncertainties between parameters. The significance level of all tests was set at 0.05.

Results

Uncertainty of size parameters

Table 2 summarizes averaged values and relative uncertainties of 12 aneurysm size parameters calculated for the five investigated IAs.

Table 2 Evaluation of size parameters

Among 1D size parameters, the relative uncertainty of Hb (#5) was considerably higher than the relative uncertainties of the other parameters. Ignoring Hb, the relative uncertainties of 3D size parameters were higher than the relative uncertainties of 2D size parameters, whose relative uncertainties were in turn higher than the 1D size parameters’ relative uncertainties.

Uncertainty of neck size parameters

Table 3 summarizes averaged values and relative uncertainties of six parameters (four 1D and two 2D) describing neck size. Relative uncertainties of all 1D neck size parameters were approximately twice as high as the relative uncertainties of 1D size parameters (#1–4), excluding Hb (#5). The relative uncertainties of both 2D neck size parameters were considerably higher than the relative uncertainties of the 1D neck size parameters.

Table 3 Evaluation of neck parameters

Figure 2 shows scatter plots between the four 1D size (excluding Hb) and the four 1D neck size parameters (#13–16). No correlation with τ > 0.7 was found between 1D size and neck size parameters.

Fig. 2
figure2

Scatter plot matrix of the four 1D size and the four 1D neck size parameters, which showed low uncertainties. Histograms of the respective parameters are shown along the diagonal entries

Uncertainty of non-dimensional aneurysm shape parameters

Table 4 summarizes averaged values and relative uncertainties for twelve non-dimensional parameters aiming to describe aneurysm shape, deviation from a (hemi)spherical shape or aneurysm size relative to parent vessel size (#19).

Table 4 Evaluation of non-dimensional parameters

Figure 3 shows scatter plots for these twelve parameters. Four pairs of parameters seem to capture similar geometric aspects: AR and BF (#20, #22) correlate with τ = 0.80, NSI and IPR (#24, #25) correlate with τ = 1, AVSV and AASA (#26, #27) correlate with τ = 0.81, and UI and CR (#28, #29) correlate with τ = − 1 (since UI = 1-CR by definition). At the same time, relative uncertainty of BF is significantly lower than AR, relative uncertainty of IPR is significantly lower than NSI, relative uncertainty of AASA is significantly lower than AVSV, and relative uncertainty of CR is significantly lower than UI (all p < 0.001).

Fig. 3
figure3

Scatter plot matrix of twelve non-dimensional shape parameters. Histograms of the respective parameters are shown along the diagonal entries

Uncertainty of aneurysm surface curvature parameters

Table 5 summarizes averaged values and relative uncertainties for eight parameters characterizing aneurysm surface curvature.

Table 5 Evaluation of curvature parameters

Two mean curvature parameters (#31 and #32) show a relatively low and similar relative uncertainty of around 15%. Both parameters are significantly correlated at τ = 0.8. They are also correlated with two 1D size parameters (#2 to #4) at τ > 0.7, as well as all 2D and 3D size parameters (#6 to #12) at τ > 0.8.

Among all other curvature parameters, only MLN (#35) has low relative uncertainty. However, only a weak correlation between MLN and the two mean curvature parameters mentioned above (#31 and #32) was observed (τ < 0.5).

Discussion

Uncertainty associated with inter-operator or inter-group variability during segmentation was quantified for a large set of geometric parameters. A large variability in uncertainty for different parameters was found, with uncertainties ranging from 3.9 to 179.8%. This information is valuable since advanced image-based methodologies are increasingly applied to assess the rupture risk of individual IAs, aiming to support clinicians in making treatment decisions for individual patients. In these research efforts, several groups have attempted to develop a model predicting rupture risk of cerebral aneurysms based on geometric [12,13,14,15,16, 34] or hemodynamic parameters [30, 31], or a combination of both parameter types, as recently proposed by Detmer et al. [35]. However, the acceptance among physicians remains limited. One reason for this limited acceptance might be the difficulty in addressing and controlling all possible sources of error as well as often controversial findings of different research groups and differences in significance between univariate and multivariate analysis [16, 36]. Therefore, this study focuses on the quantification of uncertainty introduced during the crucial step of vessel wall reconstruction (segmentation and surface post processing, e.g. smoothing).

The uncertainties and their variability found in our study were caused by different segmentation and surface reconstruction approaches from medical image data. Considering the median values of the relative uncertainty ranges of all investigated parameters, the following key findings can be summarized:

  • Non-dimensional shape parameters constructed from surface areas and/or volumes tended to exhibit the lowest uncertainties: IPR (#25), CR (#29), EI (#30). This is not surprising, since the parameters used for the calculation of these non-dimensional parameters are either overestimated or underestimated concordantly, thus reducing difference between non-dimensional values.

  • Four out of five 1D size parameters also showed comparatively low relative uncertainty values: H, Lmax, Hmax, Dmax (#1–4).

  • Relative uncertainties of size parameters grow with dimension. Relative uncertainties for 2D size parameters (i.e. areas) were larger than for 1D size parameters and relative uncertainties for 3D size parameters (i.e. volumes) were larger than for 2D size parameters. Relative uncertainties for 2D and 3D size parameters were at least twice as high as for 1D size parameters. This increase in relative uncertainty with dimension seems reasonable, given that the 2D and 3D size parameters can be expected to correlate with the square and cube of the 1D size parameters and hence the same can be expected from the relative uncertainties. This result was anticipated, since a small earlier study of the reproducibility of aneurysm segmentation found a segmentation accuracy of around one voxel [25]. Thus, the relative uncertainty of a 1D size parameter can be approximated as the voxel size divided by the 1D parameter (e.g. diameter). Here, the voxel resolution of the image data was 0.28 mm and the average aneurysm diameter was 4 mm. Thus, an estimate for the expected relative uncertainty of 1D size parameters would be 7%. Since 2D and 3D size parameters correlated with the square and cube of 1D size parameters, respectively, the relative uncertainties can also be expected to be approximately the square/cube of the 1D size parameter uncertainties.

  • Calculated uncertainties for 1D neck size parameters were approximately twice as high as for 1D size parameters. This can be explained by the aneurysm neck being particularly sensitive to segmentation uncertainty, as was shown by Berg et al. [22]. For example, different segmentation approaches can lead to the in-/exclusion of small branching vessels, which can considerably affect the neck shape and size. This may result in bifurcation aneurysms being classified as sidewall aneurysms, an example is shown in Fig. 4.

    Fig. 4
    figure4

    Example of differences in the aneurysm neck definition, depending on whether branching vessels were included in the image segmentation or not

  • Curvature parameters generally exhibit the highest relative uncertainties, with the exception of three parameters, two of which (MAA and absMAA) are measures of aneurysm shape and size, while the third parameter (MLN) is a measure of aneurysm shape alone.

Uncertainty analysis based recommendations

Based on the described findings, several recommendations can be formulated, which could be taken into account in related future studies:

  • Non-dimensional parameters seem to be preferable, since over- or undersizing during segmentation is partially compensated by the normalization, which results in lower relative uncertainty.

  • In a recent study by Liang et al. [18], which reviewed 46 morphologic and hemodynamic studies, aspect ratio was identified as the most relevant parameter for rupture risk assessment among geometric rupture risk parameters. However, our result suggest that the relative uncertainty of the aspect ratio is approximately 30%. Aneurysm height was another frequently used parameter, which exhibits a relative uncertainty of 14% in our study. Thus, processing of image data should be done particularly cautiously, if these parameters are to be evaluated.

  • Most curvature parameters exhibit high relative uncertainties in our study. This likely results from the curvature being particularly affected by both major steps of the reconstruction procedure, segmentation and smoothing. From the viewpoint of uncertainty, the inclusion of curvature parameters, with the exception of MAA, absMAA, and MLN, for rupture risk estimation warrants caution. If curvature parameters are to be used, careful standardization of image reconstruction should be aimed for.

  • In addition to relative uncertainty, the differences in parameter value between ruptured and unruptured aneurysms should also be considered, when choosing a rupture risk parameter. For example, Weir et al. found an average AR of 1.8 for unruptured and 3.4 for ruptured aneurysms [37], while we found a relative uncertainty of 29.3% for AR. This suggests that AR is able to discriminate between the two groups fairly well, even though it has only average relative uncertainty. On the other hand, Weir et al. also found an average maximum dimension of 7 mm for unruptured and 8 mm for ruptured aneurysms, while we found a relative uncertainty of 15.7% for the maximum dimension. Thus, the resulting absolute uncertainty can be expected to be similar to the difference between both groups, making the parameter less suited for discriminating between the two groups. The example shows that relative uncertainty alone is insufficient to assess whether a parameter is suited for rupture risk estimation.

Limitations

The calculation of geometric parameters requires the separation of the aneurysm head from the parent vessel, which was done manually by one operator for all data sets. Although a small inter-operator dependency analysis found minor impact on the calculation of geometric parameters [31], this procedure could be replaced by an automated approach in the future [17].

Another aspect possibly affecting the results of the uncertainty analysis is the accuracy of the geometric parameters’ measurements. Beside the voxel resolution of the acquisition discussed above, the resolution of the surface meshes provided by MATCH participants was important. The average and standard deviation of the mesh resolution was 0.199 ± 0.096 mm with a range between 0.010 mm and 0.453 mm. Note, that the accuracy of the used measurement tools is much higher than both the mesh and voxel resolution, with at least 32 bit accuracy for mesh node coordinates. Since accumulation of numerical errors might reduce this accuracy, we performed a numerical experiment on a synthetic aneurysm, which had the shape of a spherical cap (diameter 5 mm, height 4 mm, triangular mesh resolution 0.2 mm). Diameter, surface and volume of the meshed synthetic aneurysm were calculated as usual and compared to the analytical solution, the differences were 0.016%, 0.17% and 0.31%, respectively. Hence, the uncertainty introduced by the measurement tools used in our study can be neglected.

Further limitations of this study are linked to the design of MATCH: The major limitation of this study is the small number of segmented IAs. This impedes comparison of inter-aneurysm and inter-operator variability, because the former can only be assessed on the basis of five aneurysms.

Additionally, MATCH only evaluated limited information on the reconstruction procedure used by each group. Examples of missing information are the segmentation thresholds as well as the surface smoothing algorithms used.

The present study investigated the uncertainty of individual geometric rupture risk parameters only. Uncertainty of complex rupture risk models including interactions between two or more geometric risk parameters, such as proposed by Prestigiacomo et al. [38] or Detmer et al. [39], for example, was not investigated. However, all dimensionless parameters investigated here are combinations of two single rupture risk parameters and all data necessary to assess the uncertainty of other complex geometric rupture risk models is provided. Finally, uncertainty of the reconstruction procedure also affects calculation of hemodynamic parameters proposed as rupture risk parameters, which are not part of this study. Note, that the majority of complex rupture risk models include both geometric and hemodynamic risk parameters. Within the MATCH challenge, 17 groups tried to predict ruptured aneurysm. Among these groups, 12 proposed models combining geometric and hemodynamic parameters, whereas 5 used only hemodynamic rupture risk predictors.

Conclusions

Uncertainty analysis is mandatory on the road to clinical translation and use of rupture risk prediction models. The presented uncertainty analysis shows, that accurate rupture risk estimation requires reliable and reproducible assessment of geometric parameters. Thus, developing standards for reconstruction of aneurysm image data seems warranted. The data provided by this study might be able to help support development of future rupture risk prediction models by providing estimates of the uncertainties of individual parameters. In the future, an uncertainty analysis based on a larger aneurysm cohort, representing a more representative range of geometric rupture risk parameter values, is desirable.

References

  1. 1.

    Keedy A. An overview of intracranial aneurysms. Mcgill J Med. 2006;9:141–6.

  2. 2.

    Jeong Y-G, Jung Y-T, Kim M-S, Eun C-K, Jang S-H. Size and location of ruptured intracranial aneurysms. J Korean Neurosurg Soc. 2009;45:11–5. https://doi.org/10.3340/jkns.2009.45.1.11.

  3. 3.

    Thompson BG, Brown RD, Amin-Hanjani S, Broderick JP, Cockroft KM, Connolly ES, et al. Guidelines for the management of patients with unruptured intracranial aneurysms: a guideline for healthcare professionals from the american heart Association/American stroke association. Stroke. 2015;46:2368–400. https://doi.org/10.1161/STR.0000000000000070.

  4. 4.

    Wiebers DO, Whisnant JP, Huston J, Meissner I, Brown RD, Piepgras DG, et al. Unruptured intracranial aneurysms: natural history, clinical outcome, and risks of surgical and endovascular treatment. Lancet. 2003;362:103–10.

  5. 5.

    Morita A, Kirino T, Hashi K, Aoki N, Fukuhara S, Hashimoto N, et al. The natural course of unruptured cerebral aneurysms in a Japanese cohort. N Engl J Med. 2012;366:2474–82. https://doi.org/10.1056/NEJMoa1113260.

  6. 6.

    Lee G-J, Eom K-S, Lee C, Kim D-W, Kang S-D. Rupture of very small intracranial aneurysms: incidence and clinical characteristics. J Cerebrovasc Endovasc Neurosurg. 2015;17:217–22. https://doi.org/10.7461/jcen.2015.17.3.217.

  7. 7.

    Ishibashi T, Murayama Y, Urashima M, Saguchi T, Ebara M, Arakawa H, et al. Unruptured intracranial aneurysms: incidence of rupture and risk factors. Stroke. 2009;40:313–6. https://doi.org/10.1161/STROKEAHA.108.521674.

  8. 8.

    Greving JP, Wermer MJH, Brown RD, Morita A, Juvela S, Yonekura M, et al. Development of the PHASES score for prediction of risk of rupture of intracranial aneurysms: a pooled analysis of six prospective cohort studies. Lancet Neurol. 2014;13:59–66. https://doi.org/10.1016/S1474-4422(13)70263-1.

  9. 9.

    Backes D, Vergouwen MDI, Tiel Groenestege AT, Bor ASE, Velthuis BK, Greving JP, et al. PHASES score for prediction of intracranial aneurysm growth. Stroke. 2015;46:1221–6. https://doi.org/10.1161/STROKEAHA.114.008198.

  10. 10.

    Bijlenga P, Gondar R, Schilling S, Morel S, Hirsch S, Cuony J, et al. PHASES score for the management of intracranial aneurysm: a cross-sectional population-based retrospective study. Stroke. 2017;48:2105–12. https://doi.org/10.1161/STROKEAHA.117.017391.

  11. 11.

    Darsaut T, Fahed R, Raymond J. PHASES and the natural history of unruptured aneurysms: science or pseudoscience? J Neurointerv Surg. 2017;9:527–8. https://doi.org/10.1136/neurintsurg-2016-012295.

  12. 12.

    Raghavan ML, Ma B, Harbaugh RE. Quantified aneurysm shape and rupture risk. J Neurosurg. 2005;102:355–62. https://doi.org/10.3171/jns.2005.102.2.0355.

  13. 13.

    Hoh BL, Sistrom CL, Firment CS, Fautheree GL, Velat GJ, Whiting JH, et al. Bottleneck factor and height-width ratio: association with ruptured aneurysms in patients with multiple cerebral aneurysms. Neurosurgery. 2007;61:716–22. https://doi.org/10.1227/01.neu.0000298899.77097.bf.

  14. 14.

    Dhar S, Tremmel M, Mocco J, Kim M, Yamamoto J, Siddiqui AH, et al. Morphology parameters for intracranial aneurysm rupture risk assessment. Neurosurgery. 2008;63:185. https://doi.org/10.1227/01.NEU.0000316847.64140.81.

  15. 15.

    Xiang J, Natarajan SK, Tremmel M, Ma D, Mocco J, Hopkins LN, et al. Hemodynamic-morphologic discriminants for intracranial aneurysm rupture. Stroke. 2011;42:144–52. https://doi.org/10.1161/STROKEAHA.110.592923.

  16. 16.

    Lall RR, Eddleman CS, Bendok BR, Batjer HH. Unruptured intracranial aneurysms and the assessment of rupture risk based on anatomical and morphological factors: sifting through the sands of data. Neurosurg Focus. 2009;26:E2. https://doi.org/10.3171/2009.2.FOCUS0921.

  17. 17.

    Saalfeld S, Berg P, Niemann A, Luz M, Preim B, Beuing O. Semiautomatic neck curve reconstruction for intracranial aneurysm rupture risk assessment based on morphological parameters. Int J Comput Assist Radiol Surg. 2018;13:1781–93. https://doi.org/10.1007/s11548-018-1848-x.

  18. 18.

    Liang L, Steinman DA, Brina O, Chnafa C, Cancelliere NM, Pereira VM. Towards the Clinical utility of CFD for assessment of intracranial aneurysm rupture—a systematic review and novel parameter-ranking tool. J Neurointerv Surg. 2018. https://doi.org/10.1136/neurintsurg-2018-014246.

  19. 19.

    Young PG, Beresford-West TBH, Coward SRL, Notarberardino B, Walker B, Abdul-Aziz A. An efficient approach to converting three-dimensional image data into highly accurate computational models. Philos Trans A Math Phys Eng Sci. 2008;366:3155–73. https://doi.org/10.1098/rsta.2008.0090.

  20. 20.

    Sarrami-Foroushani A, Lassila T, Frangi AF. Virtual endovascular treatment of intracranial aneurysms: models and uncertainty. Wiley Interdiscip Rev Syst Biol Med. 2017. https://doi.org/10.1002/wsbm.1385.

  21. 21.

    Fossan FE, Sturdy J, Müller LO, Strand A, Bråten AT, Jørgensen A, et al. Uncertainty quantification and sensitivity analysis for computational FFR estimation in stable coronary artery disease. Cardiovasc Eng Technol. 2018;9:597–622. https://doi.org/10.1007/s13239-018-00388-w.

  22. 22.

    Berg P, Radtke L, Vos S, Serowy S, Janiga G, Preim B, Beuing O, Saalfeld S. 3DRA reconstruction of intracranial aneurysms—how does voxel size influences morphologic and hemodynamic parameters. Conf Proc IEEE Eng Med Biol Soc. 2018. https://doi.org/10.1109/embc.2018.8512524.

  23. 23.

    Berg P, Saalfeld S, Voß S, Redel T, Preim B, Janiga G, Beuing O. Does the DSA reconstruction kernel affect hemodynamic predictions in intracranial aneurysms? An analysis of geometry and blood flow variations. J Neurointerv Surg. 2018;10(3):290–6. https://doi.org/10.1136/neurintsurg-2017-012996.

  24. 24.

    Klepaczko A, Szczypiński P, Deistung A, Reichenbach JR, Materka A. Simulation of MR angiography imaging for validation of cerebral arteries segmentation algorithms. Comput Methods Programs Biomed. 2016;137:293–309. https://doi.org/10.1016/j.cmpb.2016.09.020.

  25. 25.

    Goubergrits L, Schaller J, Kertzscher U, Petz C, Hege H-C, Spuler A. Reproducibility of image-based analysis of cerebral aneurysm geometry and hemodynamics: an in vitro study of magnetic resonance imaging, computed tomography, and three-dimensional rotational angiography. J Neurol Surg A Cent Eur Neurosurg. 2013;74:294–302. https://doi.org/10.1055/s-0033-1342937.

  26. 26.

    Hua Y, Oh JH, Kim YB. Influence of parent artery segmentation and boundary conditions on hemodynamic characteristics of intracranial aneurysms. Yonsei Med J. 2015;56(5):1328–37. https://doi.org/10.3349/ymj.2015.56.5.1328.

  27. 27.

    Sen Y, Qian Y, Avolio A, Morgan M. Image segmentation methods for intracranial aneurysm haemodynamic research. J Biomech. 2014;47(5):1014–9. https://doi.org/10.1016/j.jbiomech.2013.12.035.

  28. 28.

    Castro MA, Putman CM, Cebral JR. Computational fluid dynamics modeling of intracranial aneurysms: effects of parent artery segmentation on intra-aneurysmal hemodynamics. AJNR Am J Neuroradiol. 2006;27(8):1703–9.

  29. 29.

    Berg P, Voß S, Saalfeld S, Janiga G, Bergersen AW, Valen-Sendstad K, et al. Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH): phase I: segmentation. Cardiovasc Eng Technol. 2018. https://doi.org/10.1007/s13239-018-00376-0.

  30. 30.

    Berg P, Beuing O. Multiple intracranial aneurysms: a direct hemodynamic comparison between ruptured and unruptured vessel malformations. Int J Comput Assist Radiol Surg. 2018;13:83–93. https://doi.org/10.1007/s11548-017-1643-0.

  31. 31.

    Goubergrits L, Schaller J, Kertzscher U, van den Bruck N, Poethkow K, Petz C, et al. Statistical wall shear stress maps of ruptured and unruptured middle cerebral artery aneurysms. J R Soc Interface. 2012;9:677–88. https://doi.org/10.1098/rsif.2011.0490.

  32. 32.

    Rohde S, Lahmann K, Beck J, Nafe R, Yan B, Raabe A, Berkefeld J. Fourier analysis of intracranial aneurysms: towards an objective and quantitative evaluation of the shape of aneurysms. Neuroradiology. 2005;47:121–6. https://doi.org/10.1007/s00234-004-1324-x.

  33. 33.

    Lauric A, Miller EL, Baharoglu MI, Malek AM. Rupture status discrimination in intracranial aneurysms using the centroid-radii model. IEEE Trans Biomed Eng. 2011;58:2895–903. https://doi.org/10.1109/TBME.2011.2162410.

  34. 34.

    Ma B, Harbaugh RE, Raghavan ML. Three-dimensional geometrical characterization of cerebral aneurysms. Ann Biomed Eng. 2004;32:264–73. https://doi.org/10.1023/B:ABME.0000012746.31343.92.

  35. 35.

    Detmer FJ, Chung BJ, Mut F, Slawki M, Hamzei-Sichani F, Putman C, Jiménez C, Cebral JR. Development and internal validation of an aneurysm rupture probability model based on patient characteristics and aneurysm location, morphology, and hemodynamics. Int J Comput Assist Radiol Surg. 2018;13:1767–79. https://doi.org/10.1007/s11548-018-1837-0.

  36. 36.

    Miura Y, Ishida F, Umeda Y, Tanemura H, Suzuki H, Matsushima S, Shimosaka S, Taki W. Low wall shear stress is independently associated with the rupture status of middle cerebral artery aneurysms. Stroke. 2013;44:519–21. https://doi.org/10.1161/STROKEAHA.112.675306.

  37. 37.

    Weir B, Amidei C, Kongable G, Findlay JM, Kassell NF, Kelly J, et al. The aspect ratio (dome/neck) of ruptured and unruptured aneurysms. J Neurosurg. 2003;99:447–51. https://doi.org/10.3171/jns.2003.99.3.0447.

  38. 38.

    Prestigiacomo CJ, He W, Catrambone J, Chung S, Kasper L, Pasupuleti L, Mittal N. Predicting aneurysm rupture probabilities through the application of a computed tomography angiography-derived binary logistic regression model. J Neurosurg. 2009;110(1):1–6. https://doi.org/10.3171/2008.5.17558.

  39. 39.

    Detmer FJ, Chung BJ, Mut F, Pritz M, Slawski M, Hamzei-Sichani F, Kallmes D, Putman C, Jimenez C, Cebral JR. Development of a statistical model for discrimination of rupture status in posterior communicating artery aneurysms. Acta Neurochir. 2018;160(8):1643–52. https://doi.org/10.1007/s00701-018-3595-8.

  40. 40.

    Lauric A, Miller EL, Baharoglu MI, Malek AM. 3D shape analysis of intracranial aneurysms using the writhe number as a discriminant for rupture. Ann Biomed Eng. 2011;39:1457–69. https://doi.org/10.1007/s10439-010-0241-x.

  41. 41.

    Chien A, Sayre J, Viñuela F. Comparative morphological analysis of the geometry of ruptured and unruptured aneurysms. Neurosurgery. 2011;69:349–56. https://doi.org/10.1227/NEU.0b013e31821661c3.

  42. 42.

    Ujiie H, Tamano Y, Sasaki K, Hori T. Is the aspect ratio a reliable index for predicting the rupture of a saccular aneurysm? Neurosurgery. 2001;48:495–502.

Download references

Authors’ contributions

Development of the analysis strategy and the studies concept was done by LG and PB. LG separated all aneurysm dome geometries from the parent vessel. FH calculated all geometric parameters analyzed in this study. Statistical Analysis was performed by LG, PB, FH and JB. Those four authors also wrote major parts of the manuscript. SV, GJ, SS, OB were crucially involved in MATCH and processed required data. All authors critically reviewed the manuscript, discussed findings and their relevance. Expertise of HCH And AS was especially necessary for mathematical and clinical interpretation of this study’s findings respectively. All authors read and approved the final manuscript.

Acknowledgements

We would like to thank all participants of MATCH. Without their participation this analysis would not have been possible.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The dataset supporting the conclusions of this article (i.e. all aneurysm dome geometries) is available in the figshare repository, https://doi.org/10.6084/m9.figshare.7485035 (https://figshare.com/s/a98c587dc35b79c27384).

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

This study was funded by the German Federal Ministry of Education and Research within the research campus STIMULATE (Grant Number 13GW0095A) and the German Research Foundation (Grant Number 399581926).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Correspondence to Leonid Goubergrits.

Additional file

12938_2019_657_MOESM1_ESM.xlsx

Additional file 1. All geometric parameters calculated for the aneurysm geometries provided in frames of MATCH.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Uncertainty
  • Cerebral aneurysms
  • Segmentation
  • Rupture risk