Maria-Cruz Villa-Uriol

Fast virtual deployment of self-expandable stents: Method and in vitro evaluation for intracranial aneurysmal stenting

INTRODUCTION Minimally invasive treatment approaches, like the implantation of percutaneous stents, are becoming more popular every day for the treatment of intracranial aneurysms. The outcome of such treatments is related to factors like vessel and aneurysm geometry, hemodynamic conditions and device design. For this reason, having a tool for assessing stenting alternatives beforehand is crucial. METHODOLOGY The Fast Virtual Stenting (FVS) method, which provides an estimation of the configuration of intracranial stents when released in realistic geometries, is proposed in this paper. This method is based on constrained simplex deformable models. The constraints are used to account for the stent design. An algorithm for its computational implementation is also proposed. The performance of the proposed methodology was contrasted with real stents released in a silicone phantom. RESULTS In vitro experiments were performed on the phantom where a contrast injection was performed. Subsequently, corresponding Computational Fluid Dynamics (CFD) analyzes were carried out on a digital replica of the phantom with the virtually released stent. Virtual angiographies are used to compare in vitro experiments and CFD analysis. Contrast time-density curves for in vitro and CFD data were generated and used to compare them. CONCLUSIONS Results of both experiments resemble very well, especially when comparing the contrast density curves. The use of FVS methodology in the clinical environment could provide additional information to clinicians before the treatment to choose the therapy that best fits the patient.

How do coil configuration and packing density influence intra-aneurysmal hemodynamics?

BACKGROUND AND PURPOSE: Endovascular coiling is a well-established therapy for treating intracranial aneurysms. Nonetheless, postoperative hemodynamic changes induced by this therapy remain not fully understood. The purpose of this work is to assess the influence of coil configuration and packing density on intra-aneurysmal hemodynamics. MATERIALS AND METHODS: Three 3D rotational angiography images of 3 intracranial aneurysms before and after endovascular coiling were used. For each aneurysm, a 3D representation of the vasculature was obtained after the segmentation of the images. Afterward, a virtual coiling technique was used to treat the aneurysm geometries with coil models. The aneurysms were coiled with 5 packing densities, and each was generated by using 3 coil configurations. Computational fluid dynamics analyses were carried out in both untreated and treated aneurysm geometries. Statistical tests were performed to evaluate the relative effect of coil configuration on local hemodynamics. RESULTS: The intra-aneurysmal blood flow velocity and wall shear stress were diminished as packing density increased. Aneurysmal flow velocity was reduced >50% due to the first inserted coils (packing density <12%) but with a high dependency on coil configuration. Nonsignificant differences (P > .01) were found in the hemodynamics due to coil configuration for high packing densities (near 30%). A damping effect was observed on the intra-aneurysmal blood flow waveform after coiling. CONCLUSIONS: Intra-aneurysmal hemodynamics are altered by coils. Coil configuration might reduce its influence on intra-aneurysmal hemodynamics as the packing density increases until an insignificant influence could be achieved for high packing densities.

Automated segmentation of cerebral vasculature with aneurysms in 3DRA and TOF-MRA using geodesic active regions: an evaluation study.

PURPOSE To evaluate the suitability of an improved version of an automatic segmentation method based on geodesic active regions (GAR) for segmenting cerebral vasculature with aneurysms from 3D x-ray reconstruction angiography (3DRA) and time of flight magnetic resonance angiography (TOF-MRA) images available in the clinical routine. METHODS Three aspects of the GAR method have been improved: execution time, robustness to variability in imaging protocols, and robustness to variability in image spatial resolutions. The improved GAR was retrospectively evaluated on images from patients containing intracranial aneurysms in the area of the Circle of Willis and imaged with two modalities: 3DRA and TOF-MRA. Images were obtained from two clinical centers, each using different imaging equipment. Evaluation included qualitative and quantitative analyses of the segmentation results on 20 images from 10 patients. The gold standard was built from 660 cross-sections (33 per image) of vessels and aneurysms, manually measured by interventional neuroradiologists. GAR has also been compared to an interactive segmentation method: isointensity surface extraction (ISE). In addition, since patients had been imaged with the two modalities, we performed an intermodality agreement analysis with respect to both the manual measurements and each of the two segmentation methods. RESULTS Both GAR and ISE differed from the gold standard within acceptable limits compared to the imaging resolution. GAR (ISE) had an average accuracy of 0.20 (0.24) mm for 3DRA and 0.27 (0.30) mm for TOF-MRA, and had a repeatability of 0.05 (0.20) mm. Compared to ISE, GAR had a lower qualitative error in the vessel region and a lower quantitative error in the aneurysm region. The repeatability of GAR was superior to manual measurements and ISE. The intermodality agreement was similar between GAR and the manual measurements. CONCLUSIONS The improved GAR method outperformed ISE qualitatively as well as quantitatively and is suitable for segmenting 3DRA and TOF-MRA images from clinical routine.

Biomechanical wall properties of human intracranial aneurysms resected following surgical clipping (IRRAs Project)

BACKGROUND AND PURPOSE Individual rupture risk assessment of intracranial aneurysms is a major issue in the clinical management of asymptomatic aneurysms. Aneurysm rupture occurs when wall tension exceeds the strength limit of the wall tissue. At present, aneurysmal wall mechanics are poorly understood and thus, risk assessment involving mechanical properties is inexistent. Aneurysm computational hemodynamics studies make the assumption of rigid walls, an arguable simplification. We therefore aim to assess mechanical properties of ruptured and unruptured intracranial aneurysms in order to provide the foundation for future patient-specific aneurysmal risk assessment. This work also challenges some of the currently held hypotheses in computational flow hemodynamics research. METHODS A specific conservation protocol was applied to aneurysmal tissues following clipping and resection in order to preserve their mechanical properties. Sixteen intracranial aneurysms (11 female, 5 male) underwent mechanical uniaxial stress tests under physiological conditions, temperature, and saline isotonic solution. These represented 11 unruptured and 5 ruptured aneurysms. Stress/strain curves were then obtained for each sample, and a fitting algorithm was applied following a 3-parameter (C(10), C(01), C(11)) Mooney-Rivlin hyperelastic model. Each aneurysm was classified according to its biomechanical properties and (un)rupture status. RESULTS Tissue testing demonstrated three main tissue classes: Soft, Rigid, and Intermediate. All unruptured aneurysms presented a more Rigid tissue than ruptured or pre-ruptured aneurysms within each gender subgroup. Wall thickness was not correlated to aneurysmal status (ruptured/unruptured). An Intermediate subgroup of unruptured aneurysms with softer tissue characteristic was identified and correlated with multiple documented risk factors of rupture. CONCLUSION There is a significant modification in biomechanical properties between ruptured aneurysm, presenting a soft tissue and unruptured aneurysms, presenting a rigid material. This finding strongly supports the idea that a biomechanical risk factor based assessment should be utilized in the to improve the therapeutic decision making.

neurIST complex information processing toolchain for the integrated management of cerebral aneurysms

Cerebral aneurysms are a multi-factorial disease with severe consequences. A core part of the European project @neurIST was the physical characterization of aneurysms to find candidate risk factors associated with aneurysm rupture. The project investigated measures based on morphological, haemodynamic and aneurysm wall structure analyses for more than 300 cases of ruptured and unruptured aneurysms, extracting descriptors suitable for statistical studies. This paper deals with the unique challenges associated with this task, and the implemented solutions. The consistency of results required by the subsequent statistical analyses, given the heterogeneous image data sources and multiple human operators, was met by a highly automated toolchain combined with training. A testimonial of the successful automation is the positive evaluation of the toolchain by over 260 clinicians during various hands-on workshops. The specification of the analyses required thorough investigations of modelling and processing choices, discussed in a detailed analysis protocol. Finally, an abstract data model governing the management of the simulation-related data provides a framework for data provenance and supports future use of data and toolchain. This is achieved by enabling the easy modification of the modelling approaches and solution details through abstract problem descriptions, removing the need of repetition of manual processing work.

Morphodynamic Analysis of Cerebral Aneurysm Pulsation From Time-Resolved Rotational Angiography

This paper presents a technique to estimate and model patient-specific pulsatility of cerebral aneurysms over one cardiac cycle, using 3D rotational X-ray angiography (3DRA) acquisitions. Aneurysm pulsation is modeled as a time varying B-spline tensor field representing the deformation applied to a reference volume image, thus producing the instantaneous morphology at each time point in the cardiac cycle. The estimated deformation is obtained by matching multiple simulated projections of the deforming volume to their corresponding original projections. A weighting scheme is introduced to account for the relevance of each original projection for the selected time point. The wide coverage of the projections, together with the weighting scheme, ensures motion consistency in all directions. The technique has been tested on digital and physical phantoms that are realistic and clinically relevant in terms of geometry, pulsation and imaging conditions. Results from digital phantom experiments demonstrate that the proposed technique is able to recover subvoxel pulsation with an error lower than 10% of the maximum pulsation in most cases. The experiments with the physical phantom allowed demonstrating the feasibility of pulsation estimation as well as identifying different pulsation regions under clinical conditions.

Automated landmarking and geometric characterization of the carotid siphon

The geometry of the carotid siphon has a large variability between subjects, which has prompted its study as a potential geometric risk factor for the onset of vascular pathologies on and off the internal carotid artery (ICA). In this work, we present a methodology for an objective and extensive geometric characterization of carotid siphon parameterized by a set of anatomical landmarks. We introduce a complete and automated characterization pipeline. Starting from the segmentation of vasculature from angiographic image and its centerline extraction, we first identify ICA by characterizing vessel tree bifurcations and training a support vector machine classifier to detect ICA terminal bifurcation. On ICA centerline curve, we detect anatomical landmarks of carotid siphon by modeling it as a sequence of four bends and selecting their centers and interfaces between them. Bends are detected from the trajectory of the curvature vector expressed in the parallel transport frame of the curve. Finally, using the detected landmarks, we characterize the geometry in two complementary ways. First, with a set of local and global geometric features, known to affect hemodynamics. Second, using large deformation diffeomorphic metric curve mapping (LDDMCM) to quantify pairwise shape similarity. We processed 96 images acquired with 3D rotational angiography. ICA identification had a cross-validation success rate of 99%. Automated landmarking was validated by computing limits of agreement with the reference taken to be the locations of the manually placed landmarks averaged across multiple observers. For all but one landmark, either the bias was not statistically significant or the variability was within 50% of the inter-observer one. The subsequently computed values of geometric features and LDDMCM were commensurate to the ones obtained with manual landmarking. The characterization based on pair-wise LDDMCM proved better in classifying the carotid siphon shape classes than the one based on geometric features. The proposed characterization provides a rich description of geometry and is ready to be applied in the search for geometric risk factors of the carotid siphon.

Efficient 3D Geometric and Zernike Moments Computation from Unstructured Surface Meshes

This paper introduces and evaluates a fast exact algorithm and a series of faster approximate algorithms for the computation of 3D geometric moments from an unstructured surface mesh of triangles. Being based on the object surface reduces the computational complexity of these algorithms with respect to volumetric grid-based algorithms. In contrast, it can only be applied for the computation of geometric moments of homogeneous objects. This advantage and restriction is shared with other proposed algorithms based on the object boundary. The proposed exact algorithm reduces the computational complexity for computing geometric moments up to order N with respect to previously proposed exact algorithms, from N9 to N6. The approximate series algorithm appears as a power series on the rate between triangle size and object size, which can be truncated at any desired degree. The higher the number and quality of the triangles, the better the approximation. This approximate algorithm reduces the computational complexity to N3. In addition, the paper introduces a fast algorithm for the computation of 3D Zernike moments from the computed geometric moments, with a computational complexity N4, while the previously proposed algorithm is of order N6. The error introduced by the proposed approximate algorithms is evaluated in different shapes and the cost-benefit ratio in terms of error, and computational time is analyzed for different moment orders.

Toward integrated management of cerebral aneurysms

In the last few years, some of the visionary concepts behind the virtual physiological human began to be demonstrated on various clinical domains, showing great promise for improving healthcare management. In the current work, we provide an overview of image- and biomechanics-based techniques that, when put together, provide a patient-specific pipeline for the management of intracranial aneurysms. The derivation and subsequent integration of morphological, morphodynamic, haemodynamic and structural analyses allow us to extract patient-specific models and information from which diagnostic and prognostic descriptors can be obtained. Linking such new indices with relevant clinical events should bring new insights into the processes behind aneurysm genesis, growth and rupture. The development of techniques for modelling endovascular devices such as stents and coils allows the evaluation of alternative treatment scenarios before the intervention takes place and could also contribute to the understanding and improved design of more effective devices. A key element to facilitate the clinical take-up of all these developments is their comprehensive validation. Although a number of previously published results have shown the accuracy and robustness of individual components, further efforts should be directed to demonstrate the diagnostic and prognostic efficacy of these advanced tools through large-scale clinical trials.

AngioLab-A software tool for morphological analysis and endovascular treatment planning of intracranial aneurysms

Determining whether and how an intracranial aneurysm should be treated is a tough decision that clinicians face everyday. Emerging computational tools could help clinicians analyze clinical data and make these decisions. AngioLab is a single graphical user interface, developed on top of the open source framework GIMIAS, that integrates some of the latest image analysis and computational modeling tools for intracranial aneurysms. Two workflows are available: Advanced Morphological Analysis (AMA) and Endovascular Treatment Planning (ETP). AngioLab has been evaluated by a total of 62 clinicians, who considered the information provided by AngioLab relevant and meaningful. They acknowledged the emerging need of these type of tools and the potential impact they might have on the clinical decision-making process.