Rubén Cárdenes

Patient-specific simulations of stenting procedures in coronary bifurcations: Two clinical cases

Computational simulations of stenting procedures in idealized geometries can only provide general guidelines and their use in the patient-specific planning of percutaneous treatments is inadequate. Conversely, image-based patient-specific tools that are able to realistically simulate different interventional options might facilitate clinical decision-making and provide useful insights on the treatment for each individual patient. The aim of this work is the implementation of a patient-specific model that uses image-based reconstructions of coronary bifurcations and is able to replicate real stenting procedures following clinical indications. Two clinical cases are investigated focusing the attention on the open problems of coronary bifurcations and their main treatment, the provisional side branch approach. Image-based reconstructions are created combining the information from conventional coronary angiography and computed tomography angiography while structural finite element models are implemented to replicate the real procedure performed in the patients. First, numerical results show the biomechanical influence of stents deployment in the coronary bifurcations during and after the procedures. In particular, the straightening of the arterial wall and the influence of two overlapping stents on stress fields are investigated here. Results show that a sensible decrease of the vessel tortuosity occurs after stent implantation and that overlapping devices result in an increased stress state of both the artery and the stents. Lastly, the comparison between numerical and image-based post-stenting configurations proved the reliability of such models while replicating stent deployment in coronary arteries.

Computational fluid dynamic simulations of image-based stented coronary bifurcation models.

One of the relevant phenomenon associated with in-stent restenosis in coronary arteries is an altered haemodynamics in the stented region. Computational fluid dynamics (CFD) offers the possibility to investigate the haemodynamics at a level of detail not always accessible within experimental techniques. CFD can quantify and correlate the local haemodynamics structures which might lead to in-stent restenosis. The aim of this work is to study the fluid dynamics of realistic stented coronary artery models which replicate the complete clinical procedure of stent implantation. Two cases of pathologic left anterior descending coronary arteries with their bifurcations are reconstructed from computed tomography angiography and conventional coronary angiography images. Results of wall shear stress and relative residence time show that the wall regions more prone to the risk of restenosis are located next to stent struts, to the bifurcations and to the stent overlapping zone for both investigated cases. Considering a bulk flow analysis, helical flow structures are generated by the curvature of the zone upstream from the stent and by the bifurcation regions. Helical recirculating microstructures are also visible downstream from the stent struts. This study demonstrates the feasibility to virtually investigate the haemodynamics of patient-specific coronary bifurcation geometries.

A multidimensional segmentation evaluation for medical image data

Evaluation of segmentation methods is a crucial aspect in image processing, especially in the medical imaging field, where small differences between segmented regions in the anatomy can be of paramount importance. Usually, segmentation evaluation is based on a measure that depends on the number of segmented voxels inside and outside of some reference regions that are called gold standards. Although some other measures have been also used, in this work we propose a set of new similarity measures, based on different features, such as the location and intensity values of the misclassified voxels, and the connectivity and the boundaries of the segmented data. Using the multidimensional information provided by these measures, we propose a new evaluation method whose results are visualized applying a Principal Component Analysis of the data, obtaining a simplified graphical method to compare different segmentation results. We have carried out an intensive study using several classic segmentation methods applied to a set of MRI simulated data of the brain with several noise and RF inhomogeneity levels, and also to real data, showing that the new measures proposed here and the results that we have obtained from the multidimensional evaluation, improve the robustness of the evaluation and provides better understanding about the difference between segmentation methods.

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.

Anatomical Labeling of the Circle of Willis Using Maximum A Posteriori Probability Estimation

Anatomical labeling of the cerebral arteries forming the Circle of Willis (CoW) enables inter-subject comparison, which is required for geometric characterization and discovering risk factors associated with cerebrovascular pathologies. We present a method for automated anatomical labeling of the CoW by detecting its main bifurcations. The CoW is modeled as rooted attributed relational graph, with bifurcations as its vertices, whose attributes are characterized as points on a Riemannian manifold. The method is first trained on a set of pre-labeled examples, where it learns the variability of local bifurcation features as well as the variability in the topology. Then, the labeling of the target vasculature is obtained as maximum a posteriori probability (MAP) estimate where the likelihood of labeling individual bifurcations is regularized by the prior structural knowledge of the graph they span. The method was evaluated by cross-validation on 50 subjects, imaged with magnetic resonance angiography, and showed a mean detection accuracy of 95%. In addition, besides providing the MAP, the method can rank the labelings. The proposed method naturally handles anatomical structural variability and is demonstrated to be suitable for labeling arterial segments of the CoW.

Automatic Aneurysm Neck Detection Using Surface Voronoi Diagrams

A new automatic approach for saccular intracranial aneurysm isolation is proposed in this work. Due to the inter- and intra-observer variability in manual delineation of the aneurysm neck, a definition based on a minimum cost path around the aneurysm sac is proposed that copes with this variability and is able to make consistent measurements along different data sets, as well as to automate and speedup the analysis of cerebral aneurysms. The method is based on the computation of a minimal path along a scalar field obtained on the vessel surface, to find the aneurysm neck in a robust and fast manner. The computation of the scalar field on the surface is obtained using a fast marching approach with a speed function based on the exponential of the distance from the centerline bifurcation between the aneurysm dome and the parent vessels. In order to assure a correct topology of the aneurysm sac, the neck computation is constrained to a region defined by a surface Voronoi diagram obtained from the branches of the vessel centerline. We validate this method comparing our results in 26 real cases with manual aneurysm isolation obtained using a cut-plane, and also with results obtained using manual delineations from three different observers by comparing typical morphological measures.

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.

Fast and accurate geodesic distance transform by ordered propagation

In this paper, we present a new geodesic distance transform that uses a non-Euclidean metric suitable for non-convex discrete 2D domains. The geodesic metric used is defined as the shortest path length through a set of pixels called Locally Nearest Hidden Pixels, and manages visibility zones using bounding angles. The algorithm is designed using ordered propagation, which makes it extremely efficient and linear in the number of pixels in the domain. We have compared our algorithm with the four most similar geodesic distance transform techniques, and we show that our approach has higher accuracy and lower computational complexity.

An Efficient Algorithm for Multiple Sclerosis Lesion Segmentation from Brain MRI

We propose a novel method for the segmentation of Multiple Sclerosis (MS) lesions in MRI. The method is based on a three-step approach: first a conventional k-NN classifier is applied to pre-classify gray matter (GM), white matter (WM), cerebro-spinal fluid (CSF) and MS lesions from a set of prototypes selected by an expert. Second, the classification of problematic patterns is resolved computing a fast distance transformation (DT) algorithm from the set of prototypes in the Euclidean space defined by the MRI dataset. Finally, a connected component filtering algorithm is used to remove lesion voxels not connected to the real lesions. This method uses distance information together with intensity information to improve the accuracy of lesion segmentation and, thus, it is specially useful when MS lesions have similar intensity values than other tissues. It is also well suited for interactive segmentations due to its efficiency. Results are shown on real MRI data as wall as on a standard database of synthetic images.

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.