Skip to main content

Reproducibility of image-based computational models of intracranial aneurysm: a comparison between 3D rotational angiography, CT angiography and MR angiography

Abstract

Background

Reconstruction of patient-specific biomechanical model of intracranial aneurysm has been based on different imaging modalities. However, different imaging techniques may influence the model geometry and the computational fluid dynamics (CFD) simulation. The aim of this study is to evaluate the differences of the morphological and hemodynamic parameters in the computational models reconstructed from computed tomography angiography (CTA), magnetic resonance angiography (MRA) and 3D rotational angiography (3DRA).

Methods

Ten patients with cerebral aneurysms were enrolled in the study. MRA, CTA and 3DRA were performed on all patients. For each patient, three patient-specific models were reconstructed respectively based on the three sets of imaging data of the patient. CFD simulations were performed on each model. Model geometry and hemodynamic parameters were compared between the three models.

Results

In terms of morphological parameters, by comparing CTA based models (CM) and 3DRA based models (DM) which were treated as the “standard models”, the aspect ratio had the minimum difference (Δ = 8.3 ± 1.72 %, P = 0.953) and the surface distance was 0.25 ± 0.07 mm. Meanwhile, by comparing MRA based models (MM) and DM, the size had the minimum difference (Δ = 6.6 ± 1.85 %, P = 0.683) and the surface distance was 0.36 ± 0.1 mm. In respect of hemodynamic parameters, all three models showed a similar distribution: low average WSS at the sack, high OSI at the body and high average WSSG at the neck. However, there was a large variation in the average WSS (Δ = 34 ± 5.13 % for CM, Δ = 40.6 ± 9.21 % for MM).

Conclusion

CTA and MRA have no significant differences in reproducing intracranial aneurysm geometry. The CFD results suggests there might be some significant differences in hemodynamic parameters between the three imaging-based models and this needs to be considered when interpreting the CFD results of different imaging-based models. If we only need to study the main flow patterns, three types of image-based model might be all suitable for patient-specific computational modeling studies.

Background

Intracranial aneurysms are common neurovascular disease and the most serious consequences are rupture leading to subarachnoid hemorrhage [1]. Hemodynamics may play an important role in the process of aneurysm formation, progression and rupture [25]. It has been shown that the size of intracranial aneurysms and other morphological parameters are related to their risk of rupture [6]. Moreover, the change of the hemodynamic parameters, such as wall shear stress (WSS) or wall shear stress gradient (WSSG),may influence the behavior of endothelial cells and smooth muscle cells [7], resulting in flow-mediated vasodilatation and vascular remodeling.

With the development of CT angiography (CTA), magnetic resonance angiography (MRA) and 3D rotational angiography (3DRA) techniques, the image-based patient-specific analysis of the hemodynamics of intracranial aneurysm is becoming important for diagnosis and treatment plan. Image-based analysis of intracranial aneurysms can help with the assessment of the potential risk of aneurysm rupture [8, 9] and provide guidance for treatment choice [10, 11]. Compared with other image technologies, 3DRA can depict considerably more small additional angiographic aneurysms [12, 13]. In addition, the reconstructed 3DRA images can show only the enhanced vascular lumina which allow observing any desired region without hindering over projecting bony structures [14]. Therefore, the diagnosis and measurement of aneurysms can be performed more accurately by using 3DRA technology. However, because 3DRA is invasive and expensive, most clinical imaging for cerebral aneurysm prefer to CTA and MRA rather than 3DRA.

Recently, computational fluid dynamics (CFD) plays an important role in researching cerebral aneurysm. Different imaging characteristics of CTA, MRA and 3DRA may have an impact on the reconstruction of the aneurysm models in the CFD simulation. Comparison studies of the model generation from different images have been presented in previous work [15, 16]. A previous study pointed out that even if the reconstructed aneurysm model has minor changes, there might still be a large difference in the hemodynamic results [17]. The purpose of this study is to evaluate the difference in both model geometry and fluid mechanics parameters between the reconstructed models from CTA, MRA and 3DRA.

Methods

Patients and imaging

Clinical study was performed at the Department of Radiology in Nanjing General Hospital. Ten patients (three males and seven females; mean age 56.4 years) with saccular cerebral aneurysms were enrolled in the study.

3DSA image data were obtained from Axiom Aritistd TA digital subtraction angiography (Siemens Medical Systems, Germany). The reconstruction of artery was performed in the 3D-image post processing workstation (Siemens). CTA image data were obtained from Somatom Definition AS 128 spiral CT scanner (Siemens Medical Systems, Germany). MRA image data were obtained from Achieva_1.5T Nova-Dual magnetic resonance (Philips Healthcare, The Netherlands), using maximum intensity projection (MIP) and multi-planner reformation (MPR) method to process vascular image reconstruction.

All patients with incidental findings of cerebral aneurysms at MRA and CTA underwent further evaluation of 3DRA within a 5-day interval. In this way we obtained the image data of 3DRA, CTA and MRA for each patient. Three of the ten aneurysms were found to be ruptured. The patient characteristics and aneurysm locations are shown in Table 1. This study was approved by the internal review board (Research Committee for Clinical Research of Zhongda Hospital, Affiliated to Southeast University), and informed written consent was obtained.

Table 1 Patient characteristics and aneurysm location

Model reconstruction

Three different patient-specific models were reconstructed for each patient based on the three different imaging modalities. In brief, 2D cross-sectional images of the cerebral aneurysms and its parent arteries were imported into ScanIP, Version 6.0 (Simpleware Ltd) for image segmentation and reconstruction. The segmentation procedure mainly detected the optimal boundary of the artery surface with different segmentation method. Firstly, a rough model was obtained by using the threshold segmentation, and then the algorithm of region-growing was used to remove the unconnected regions, conserving the aneurysm and its parent vessel. The small branches of the parent artery were removed artificially, based on the published studies [18]. Surface smoothing was performed with recursive Gaussian filter to reduce the sharp corners. Figure 1 shows the procedure of the patient-specific model development.

Fig. 1
figure 1

Model reconstruction of patient-specific model. a A original image data, b a rough model by using the threshold segmentation and artificially adjusting, c the final model after surface smoothing, d showing different areas in the model. Arrows represent the aneurysm

The vascular models were converted to a stereolithography (STL) format and exported to Workbench, Version 15.0 (ANSYS Inc.). In order to visualize and analyze the region of interests (ROI), for example, the aneurysm or the parent vessel, we divided the model into different regions. The ROI extent of the parent vessel applied to all models was 1.5 times of the diameter of the proximal vessel neck. In addition, an extended vessel was added to the inlet to allow the fully developed flow formation before entering the vascular model (Fig. 1d).

CFD simulation

Unstructured meshes were created with ICEM CFD, Version 15.0 (ANSYS Inc.). The maximum element size was 0.25 mm with a minimum size of 0.1 mm for the high curvature regions and the surface of the aneurysms. On average, meshes of each model consisted of 0.15 million nodes and 0.85 million elements.

Blood was treated as incompressible, viscous Newtonian fluid (ρ = 1056 kg/m3, μ = 0.0035 Pa·s), and the vessel wall was assumed to be rigid with no-slip boundary condition. The range of Reynolds numbers at the inlets was 500–750. Transient CFD simulation was performed with Fluent, Version 15.0 (ANSYS Inc.), which uses a finite volume approach to solve the Navier-Stokes equations. Based on phase contrast (PC) MRA [19], a pulsatile volume flow waveform of a healthy volunteer in the ICA was imposed at the inlet. According to the inlet volume flow rate, a simple relation assuming constant radius was employed to obtain the mean blood velocity at the inlet section as:

$$v_{mean} \left( t \right) = \frac{q(t)}{{\pi (d/2)^{2} }}$$
(1)

where d is the inlet diameter, t is time and v mean is the mean velocity. If the inlet was not located on the ICA, the waveform was scaled to simulate the physiologically realistic flow rate. The same inlet flow rate of the three models for each patient was used in the simulation. An outflow boundary with weight of 1 was used when the model had only one outlet. If there were two outlets, an appropriate flow rate relation was set based on the normal physiological blood flow distribution rate of different arteries. In this study, two outlets were the middle cerebral artery and anterior cerebral artery. The flow distribution ratio between MCA and ACA was 0.65:0.35 [20].

In the simulation, the convergence criterion was satisfied when the residual of continuity was less than 10−4 and the residual of velocity component was less than 5.0 × 10−5.

Parameter comparison and analysis

Seven morphological parameters were investigated, namely size [21], aspect ratio (AR) [22, 23], neck area (NA), parent vessel diameter (PVD) [24], size ratio (SR) [24], aneurysm angle (AA) [24] and surface distance (Table 2).

Table 2 Definition of the morphological parameters

When defining AR, the average neck diameter was calculated by the maximum diameter and the minimum diameter. Besides, the PVD was obtained by measuring two representative vessel cross sections of the parent vessel (D1 at the proximal neck and D2 at 1.5 × D1 upstream). The diameter of the parent vessel was calculated in the same way as the neck diameter. It must be emphasized that, when defining SR, the maximum aneurysm height was the distance from the centroid of the aneurysm neck to any point of the dome, rather than the maximum perpendicular height. The illustration of the morphological parameters and formulas are presented in Fig. 2. The surface distance, which was different compared with other six parameters, was calculated with a tool available in Amira, Version 5.4 (Visage Imaging Inc.). When calculating the surface distance, the two models must be approximately overlapped by using automatically aligned function. In this paper, 3DRA based models (DM) were treated as the “standard models”, CTA based models (CM) and MRA based models (MM) were compared with the “standard model” and the surface distances were calculated, respectively.

Fig. 2
figure 2

The illustration of the morphological parameters. a Definition of size and AR, b definition of AA and SR

Nine hemodynamic parameters were studied, namely average wall shear stress (AWSSA, AWSSP), low wall shear stress (LWSSA), maximum wall shear stress (MWSSA), 90th percentile value of wall shear stress (90WSSA), oscillatory shear index (OSIA), maximum oscillatory shear index (MOSIA), high oscillatory shear index (HOSIA) and average wall shear stress gradient (AWSSGA). The subscript “A” and “P” represents aneurysm and parent vessel, respectively. These nine hemodynamic parameters were derived from three main parameters: AWSS, OSI and AWSSG. They were defined as:

$${\text{AWSS}} = \left| {\frac{1}{T}\int\nolimits_{0}^{T} {\overrightarrow {WSS} dt} } \right|$$
(2)
$${\text{OSI}} = \frac{1}{2}\left( {1 - \frac{{\left| {\frac{1}{T}\int_{0}^{T} {\overrightarrow {WSS} } } \right|dt}}{{\frac{1}{T}\int_{0}^{T} {\left| {\overrightarrow {WSS} } \right|dt} }}} \right)$$
(3)
$${\text{AWSSG}} = \frac{1}{T}\int\nolimits_{0}^{T} {\sqrt {\left| {grad(WSS_{x} )} \right|^{2} \,+\, \left| {grad(WSS_{y} )} \right|^{2} + \left| {grad(WSS_{z} )} \right|^{2} d}}$$
(4)

The LWSSA was defined as the portion of aneurysm wall under low WSS (<0.4 Pa) at the end diastole [25] and the MWSSA represented the maximum WSS on the aneurysm wall at the peak systole [26]. These two parameters were believed to be important in the rupture of the aneurysm. The 90WSSA calculated the 90th percentile value of the WSS on the aneurysm at peak systole [27]. OSI is a non-dimensional parameter with a value between 0–0.5, and it is the measurement of the directional change of WSS during the cardiac cycle [28]. MOSIA and HOSIA refer to the maximum OSI and the portion of aneurysm wall above the high OSI (>0.2) at the peak systole, respectively.

According to the classification proposed by Cebral et al, three flow characteristics were compared in this paper, naming flow pattern (FP), size of the impingement region (IS) and size of the inflow jet (JS), respectively. Basedon the streamlines of velocity and the distribution of WSS and OSI, two independent observers assessed the three flow characteristics and a third observer re-examined when the first two had different results.

Two observers independently performed the measurement of the morphological parameters of different patient-specific models and the mean values were used for analysis. The hemodynamic parameters were obtained by Fluent and Tecplot 360, Version 2014 (Tecplot Inc.). CFD-Post, Version 15.0 (ANSYS Inc.) was used for visualization and the assessment of flow characteristics.

3DRA was considered as the gold standard in detecting the intracranial aneurysms. Patients were imaged more than once with each imaging modality and that the models based on 3DRA showed the least variance. Therefore, in this study, we treated the DM as the “standard model” and other two models (CM and MM) were compared with DM, respectively. The difference of each parameter was calculated as: (|CM (or MM) − DM|)/DM × 100 %. The means and standard errors of all morphological and hemodynamic parameters were calculated. The differences between DM and CM or DM and MM were analyzed by a paired nonparametric Wilcoxon test. Differences were considered statistically signification for P ≤ 0.05. As for the assessment of flow characteristic, agreement between DM and CM or DM and MM were tested with Kappa test. Agreement was categorized as poor (k < 0.4), moderate (0.4 < k < 0.75), or good (k > 0.75). Statistical analyses were performed using SPSS 19.0.

Results

Comparison of morphological parameters

All values for means and SEs for each morphological parameter are displayed in Table 3. Figure 3a presents the boxplot about difference of the morphology parameters between CM and DM. Wherein, AR had the minimum difference (Δ = 8.3 ± 1.39 %, P = 0.953), followed by size (Δ = 8.7 ± 1.39 %, P = 0.333). However, NA had the maximum difference (Δ = 21.5 ± 5.5 %, P = 0.799) with the biggest difference over 40 %. Generally, CM was similar with DM in the above morphological parameters with the maximum difference less than 25 %, except for NA. It may be the partial volume effect in CT imaging which had the biggest influence of measuring the area. As a result, the NA may be overestimated in modeling.

Table 3 Differences between models of morphologic parameters (mean ± SE)
Fig. 3
figure 3

Differences distribution for morphologic parameters. a Difference between CM and MM, b difference between MM and DM

Figure 3b shows the boxplot about difference between MM and DM. In these parameters, size had the minimum difference (Δ = 6.6 ± 1.85 %, P = 0.683). On the whole, the distribution range of differences was greater, suggesting that MRA had more disturbances in imaging the intracranial aneurysm. Nevertheless, the maximum difference was less than 30 % and the measuring results could also reflect the basic geometry of the aneurysm.

Surface differences were presented in Fig. 4, which shows the distribution of the surface distance of a cerebral aneurysm located in ICA between CM (left), MM (right) and DM, respectively (patient No. 1). Figure 5 shows the surface distance (mean and standard deviation) of the ten patients. As for the surface distance between DM and CM (grey column), the maximum mean distance occurred in patient 1 (0.37 ± 0.17 mm) and the minimum mean distance occurred in patient 3 (0.15 ± 0.15 mm). As for the surface distance between DM and MM (red column), the maximum and the minimum mean distance occurred in patient 3 (0.49 ± 0.29 mm) and patient 9 (0.19 ± 0.14 mm). DM was more closely approximated by CM than by MM. The mean surface distance of the ten patients was (0.25 ± 0.07 mm) for CM and (0.36 ± 0.1 mm) for MM.

Fig. 4
figure 4

Surface distance distribution of a cerebral aneurysm located in ICA

Fig. 5
figure 5

Surface distance data (mean and standard deviation) of the patients

Comparison of hemodynamic parameters

Figure 6 shows the AWSS, OSI and AWSSG maps of a cerebral aneurysm located in ACA-A1 and its parent artery (patient No. 5). We can see all of the three models (DM, CM and MM) shared a similar distribution of the three hemodynamic parameters: low AWSS at the sack, high OSI at the body and high AWSSG at the neck. Similar phenomena can be found in the other nine cases.

Fig. 6
figure 6

The distribution of AWSS, OSI and AWSSG in the three image based models (DM, CM and MM)

The differences information boxplot of nine hemodynamic parameters is presented in Fig. 7. All quantitative values for means and SEs for each hemodynamic parameter were displayed in Table 4. Overall, the differences of hemodynamic parameters appeared to be relatively significant, especially for the parameter of HOSI. AWSSA was usually used when studying the hemodynamics of aneurysm. Here, there was a large variation in the average WSS, either the difference between CM and DM (Δ = 34 ± 5.13 %, P = 0.646) or between MM and DM (Δ = 40.6 ± 9.21 %, P = 0.575). Compared with other parameters, the difference of MOSIA was smaller relatively (Δ = 7.3 ± 1.79 % for CM, Δ = 13.8 ± 3.13 % for MM). Meanwhile, in some cases, the difference about HOSI has the biggest value which is larger than 100 % in both two charts.

Fig. 7
figure 7

Differences distribution for hemodynamic parameters. a Difference between CM and DM, b difference between MM and DM

Table 4 Differences between models of hemodynamic parameters (mean ± SE)

Comparison of flow characteristics

Figure 8 shows the streamline colored with velocity magnitude at the peak systole of two aneurysms (patient No. 4, No. 5). As for the case located on ACA-A1 (top row in Fig. 8), they had the same FP (unchanging direction of inflow jet with single stable vortex), same IS (large) and same JS (large). As for the case located on ACoA (bottom row in Fig. 8), all the three models shared the same FP (unchanging direction of inflow jet with multiple stable vortices), but the IS and JS were small for CM while they were lager for DM and MM. Other aneurysms were also assessed in this way and the results are shown in Table 5. Both CM and MM had a good agreement in assessing the flow characteristics when comparing with DM.

Fig. 8
figure 8

Velocity-colored streamline at the peak systole of two aneurysms located on ACA-A1 (top row) and ACoA (bottom row). Models based on different imaging techniques were present from left to right (DM, CM and MM).The black arrow shows the flow direction

Table 5 Agreement judgment of FP, IS and JS between different image based models

Discussion and conclusion

As for intracranial aneurysm, both morphological and hemodynamic parameters are closely associated with the process of aneurysm formation, progression and rupture. Each imaging modality has its own unique features which may affect the measurement of morphological parameter and consequently hemodynamic parameters [29, 30]. 3DRA has greater patient risk and discomfort because it is more invasive than the others. Moreover, 3DRA is highly sensitive to the inhomogeneity caused by the contrast agent filling process during acquisition [31]. In CTA images, it may be hard to segment near-bone vasculature because the contrast agent and the intensity values overlap which consequently affects the accuracy of model reconstruction from CTA images. In MRA images, there is some signal loss intra large aneurysms with disturbed or slow flow and therefore cannot display the entire structure of the aneurysm [32]. Although 3DRA is the gold standard, CTA and MRA are still widely used in the diagnosis and treatment of aneurysms in clinical. In general, CTA is usually used in diagnosis and follow-up studies, and 3DRA is usually used before and during treatment [33]. MRA is also used because of its advantages such as cheap, non-invasive and non- radioactive. In our study, 3DRA was used as the reference standard because it has the highest spatial resolution and lowest visibility [14, 34, 35] of bone which result in better anatomic accuracy among the three modalities [12, 36].

Comparison studies of the model generation from different images have been presented in previous work. Goubergrits et al. have compared the difference of 3DRA, CTA and MRA in cerebral aneurysm geometry and hemodynamics based on an in-vitro study. They found that all three imaging modalities adequately reproduce aneurysm geometry and allow meaningful CFD analyses [12]. Piotin et al. [26] have assessed the impact of CTA, MRA and 3DRA for measuring the volume of an in-vitro aneurysm model and found that 3DRA is more accurate than CTA, while CTA is more accurate than MRA. However, most of the studies above were based on in-vitro experiment and had a phantom model to be compared. Based on ten aneurysms imaged with 3DRA and CTA, Geers et al. [20] have evaluated the impact of imaging technique on model reconstruction and CFD simulation. Geers’ results show that, although differences of quantitative hemodynamic variables were relatively large, the main flow characteristics were in excellent agreement. This is consistent with our findings. Meanwhile, in the recent research about hemodynamics impact related to differences in geometry, Schneider et al. [15] evaluated the neck size overestimation’ effect on hemodynamic features. They found that the lager neck size can result in a different impingement region and the area of low WSS. The study of Hoiet al. [16] also showed that differences in neck size can result in different flow patterns. In this paper, we find patient No. 3 has the largest surface distance difference, meanwhile its hemodynamic parameters (like LWSSA, OSIA, HOSI) have similar largest differences. This is suggest that differences in geometry are related to hemodynamics results.

The purpose of this study is to evaluate the difference in model geometry and fluid mechanics parameters of reconstructed models from 3DRA, CTA, and MRA. To our knowledge, this is the first in-vivo study to assess the impact of the three imaging modalities on model reconstruction and CFD simulation. Ten aneurysms, imaged with 3DRA, CTA and MRA, were reconstructed to computational model DM, CM and MM, respectively. Seven morphologic parameters and nine hemodynamic parameters were chosen to compare and analyze the effects of different imaging method on the modeling and CFD simulation.

In the comparison of morphology, CM has a better reproducibility than MM for most parameters. Difference between CM and DM mainly lies within 15 % with biggest difference of nearly 25 %, except for NA. Although MM has a bigger difference range, it has also an acceptable reproducibility with biggest difference less than 35 %. These differences may be resulted from different imaging technique or artificial measurement errors. However, all of the three modalities has no significant differences and allow meaningful measurement of morphologic parameters. In the comparison of hemodynamics, DM, CM, MM share a similar distribution of the main hemodynamic parameters, low AWSS at the sack, high OSI at the body and high AWSSG at the neck. However, quantitative comparison reveals that both CM and MM have significant difference from DM. Especially, the difference about HOSI has the biggest value which is larger than 100 % in some cases. This is because even though there are only small differences in morphology of some individuals, those small differences can result in obvious differences in hemodynamic parameters.

In conclusion, 3DRA, CTA and MRA have no significant differences in reproducing intracranial aneurysm geometry. The CFD results suggests there might be some significant differences in hemodynamic parameters between the three imaging-based models and this needs to be considered when interpreting the CFD results of different imaging-based models. If we only need to study the main intra-aneurysmal flow patterns, three types of image-based model might be all suitable for patient-specific computational modeling studies.

References

  1. Schievink WI. Intracranial aneurysms. N Engl J Med. 1997;336(1):28–40.

    Article  Google Scholar 

  2. Hashimoto T, Meng H, Young WL. Intracranial aneurysms: links among inflammation, hemodynamics and vascular remodeling. Neurol Res. 2006;28(4):372–80.

    Article  Google Scholar 

  3. Ho A, Lin N, Charoenvimolphan N, Stanley M, Frerichs KU, Day AL, Du R. Morphological parameters associated with ruptured posterior communicating aneurysms. PLoS ONE. 2014;9(4):e94837.

    Article  Google Scholar 

  4. Kaspera W, Ladzinski P, Larysz P, Hebda A, Ptaszkiewicz K, Kopera M, Larysz D. Morphological, hemodynamic, and clinical independent risk factors for anterior communicating artery aneurysms. Stroke. 2014;45(10):2906–11.

    Article  Google Scholar 

  5. Sforza DM, Putman CM, Cebral JR. Hemodynamics of cerebral aneurysms. Annu Rev Fluid Mech. 2009;41:91–107.

    Article  MATH  Google Scholar 

  6. Weir B, Disney L, Karrison T. Sizes of ruptured and unruptured aneurysms in relation to their sites and the ages of patients. J Neurosurg. 2002;96(1):64–70.

    Article  Google Scholar 

  7. Luscher TF, Tanner FC. Endothelial regulation of vascular tone and growth. Am J Hypertens. 1993;6(7 Pt 2):283S–93S.

    Google Scholar 

  8. Cebral JR, Castro MA, Burgess JE, Pergolizzi RS, Sheridan MJ, Putman CM. Characterization of cerebral aneurysms for assessing risk of rupture by using patient-specific computational hemodynamics models. AJNR Am J Neuroradiol. 2005;26(10):2550–9.

    Google Scholar 

  9. Cebral JR, Hendrickson S, Putman CM. Hemodynamics in a lethal basilar artery aneurysm just before its rupture. AJNR Am J Neuroradiol. 2009;30(1):95–8.

    Article  Google Scholar 

  10. Cha KS, Balaras E, Lieber BB, Sadasivan C, Wakhloo AK. Modeling the interaction of coils with the local blood flow after coil embolization of intracranial aneurysms. J Biomech Eng. 2007;129(6):873–9.

    Article  Google Scholar 

  11. Kakalis NM, Mitsos AP, Byrne JV, Ventikos Y. The haemodynamics of endovascular aneurysm treatment: a computational modelling approach for estimating the influence of multiple coil deployment. IEEE Trans Med Imaging. 2008;27(6):814–24.

    Article  Google Scholar 

  12. van Rooij WJ, Sprengers ME, de Gast AN, Peluso JP, Sluzewski M. 3D rotational angiography: the new gold standard in the detection of additional intracranial aneurysms. AJNR Am J Neuroradiol. 2008;29(5):976–9.

    Article  Google Scholar 

  13. Hernandez M, Frangi AF. Non-parametric geodesic active regions: method and evaluation for cerebral aneurysms segmentation in 3DRA and CTA. Med Image Anal. 2007;11(3):224–41.

    Article  Google Scholar 

  14. Sugahara T, Korogi Y, Nakashima K, Hamatake S, Honda S, Takahashi M. Comparison of 2D and 3D digital subtraction angiography in evaluation of intracranial aneurysms. AJNR Am J Neuroradiol. 2002;23(9):1545–52.

    Google Scholar 

  15. Schneiders JJ, Marquering HA, Antiga L, van den Berg R, VanBavel E, Majoie CB. Intracranial aneurysm neck size overestimation with 3D rotational angiography: the impact on intra-aneurysmal hemodynamics simulated with computational fluid dynamics. AJNR Am J Neuroradiol. 2013;34(1):121–8.

    Article  Google Scholar 

  16. Hoi Y, Meng H, Woodward SH, Bendok BR, Hanel RA, Guterman LR, Hopkins LN. Effects of arterial geometry on aneurysm growth: three-dimensional computational fluid dynamics study. J Neurosurg. 2004;101(4):676–81.

    Article  Google Scholar 

  17. Goubergrits L, Schaller J, Kertzscher U, Petz C, Hege HC, 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(5):294–302.

    Article  Google Scholar 

  18. Cebral JR, Castro MA, Appanaboyina S, Putman CM, Millan D, Frangi AF. Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: technique and sensitivity. IEEE Trans Med Imaging. 2005;24(4):457–67.

    Article  Google Scholar 

  19. Zhao M, Amin-Hanjani S, Ruland S, Curcio AP, Ostergren L, Charbel FT. Regional cerebral blood flow using quantitative MR angiography. AJNR Am J Neuroradiol. 2007;28(8):1470–3.

    Article  Google Scholar 

  20. Jing F, Per O. Intracerebral regional distribution of blood flow in response to uterine contractions in growth-restricted human fetuses. Early Hum Dev. 2007;83(9):607–12.

    Article  Google Scholar 

  21. Raghavan ML, Ma B, Harbaugh RE. Quantified aneurysm shape and rupture risk. J Neurosurg. 2005;102(2):355–62.

    Article  Google Scholar 

  22. Ujiie H, Tachibana H, Hiramatsu O, Hazel AL, Matsumoto T, Ogasawara Y, Nakajima H, Hori T, Takakura K, Kajiya F. Effects of size and shape (aspect ratio) on the hemodynamics of saccular aneurysms: a possible index for surgical treatment of intracranial aneurysms. Neurosurgery. 1999;45(1):119–29 (discussion 129–30).

    Article  Google Scholar 

  23. 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(3):495–502 (discussion 502–3).

    Article  Google Scholar 

  24. Dhar S, Tremmel M, Mocco J, Kim M, Yamamoto J, Siddiqui AH, Hopkins LN, Meng H. Morphology parameters for intracranial aneurysm rupture risk assessment. Neurosurgery. 2008;63(2):185–96 (discussion 196–7).

    Article  Google Scholar 

  25. Jou LD, Lee DH, Morsi H, Mawad ME. Wall shear stress on ruptured and unruptured intracranial aneurysms at the internal carotid artery. AJNR Am J Neuroradiol. 2008;29(9):1761–7.

    Article  Google Scholar 

  26. Castro M, Putman C, Radaelli A, Frangi A, Cebral J. Hemodynamics and rupture of terminal cerebral aneurysms. Academ Radiol. 2009;16(10):1201–7.

    Article  Google Scholar 

  27. Geers AJ, Larrabide I, Radaelli AG, Bogunovic H, Kim M, Gratama van Andel HA, Majoie CB, VanBavel E, Frangi AF. Patient-specific computational hemodynamics of intracranial aneurysms from 3D rotational angiography and CT angiography: an in vivo reproducibility study. AJNR Am J Neuroradiol. 2011;32(3):581–6.

    Article  Google Scholar 

  28. He X, Ku DN. Pulsatile flow in the human left coronary artery bifurcation: average conditions. J Biomech Eng. 1996;118(1):74–82.

    Article  Google Scholar 

  29. Majoie CB, Sprengers ME, van Rooij WJ, Lavini C, Sluzewski M, van Rijn JC, den Heeten GJ. MR angiography at 3T versus digital subtraction angiography in the follow-up of intracranial aneurysms treated with detachable coils. AJNR Am J Neuroradiol. 2005;26(6):1349–56.

    Google Scholar 

  30. Nijenhuis RJ, Jacobs MJ, Jaspers K, Reinders M, van Engelshoven JM, Leiner T, Backes WH. Comparison of magnetic resonance with computed tomography angiography for preoperative localization of the Adamkiewicz artery in thoracoabdominal aortic aneurysm patients. J Vasc Surg. 2007;45(4):677–85.

    Article  Google Scholar 

  31. Ernemann UU, Gronewaller E, Duffner FB, Guervit O, Claassen J, Skalej MD. Influence of geometric and hemodynamic parameters on aneurysm visualization during three-dimensional rotational angiography: an in vitro study. AJNR Am J Neuroradiol. 2003;24(4):597–603.

    Google Scholar 

  32. Piotin M, Gailloud P, Bidaut L, Mandai S, Muster M, Moret J, Rufenacht DA. CT angiography, MR angiography and rotational digital subtraction angiography for volumetric assessment of intracranial aneurysms. An experimental study. Neuroradiology. 2003;45(6):404–9.

    Article  Google Scholar 

  33. Brisman JL, Song JK, Newell DW. Cerebral aneurysms. N Engl J Med. 2006;355(9):928–39.

    Article  Google Scholar 

  34. Anxionnat R, Bracard S, Ducrocq X, Trousset Y, Launay L, Kerrien E, Braun M, Vaillant R, Scomazzoni F, Lebedinsky A, Picard L. Intracranial aneurysms: clinical value of 3D digital subtraction angiography in the therapeutic decision and endovascular treatment. Radiology. 2001;218(3):799–808.

    Article  Google Scholar 

  35. White PM, Wardlaw JM, Easton V. Can noninvasive imaging accurately depict intracranial aneurysms? A systematic review. Radiology. 2000;217(2):361–70.

    Article  Google Scholar 

  36. Tanoue S, Kiyosue H, Kenai H, Nakamura T, Yamashita M, Mori H. Three-dimensional reconstructed images after rotational angiography in the evaluation of intracranial aneurysms: surgical correlation. Neurosurgery. 2000;47(4):866–71.

    Article  Google Scholar 

Download references

Authors’ contributions

YR carried out the imaging analysis, computational analysis, participated in data analysis and drafted the manuscript. GC carried outpatient recruitment and imaging analysis. ZL participated in the computational simulation. YC participated in the design of the study and helped to draft the manuscript. GL participated in its design and coordination. ZL participated in the design of the study, data analysis and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This study was partially supported by the National 973 Basic Research Program of China [No. 2013CB733800], the National Natural Science Foundation of China (NSFC) (Nos. 11272091, 11422222, 31470043) and ARC (FT140101152).

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Yan Cai, Guang-Ming Lu or Zhi-Yong Li.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ren, Y., Chen, GZ., Liu, Z. et al. Reproducibility of image-based computational models of intracranial aneurysm: a comparison between 3D rotational angiography, CT angiography and MR angiography. BioMed Eng OnLine 15, 50 (2016). https://doi.org/10.1186/s12938-016-0163-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12938-016-0163-4

Keywords