Borja Rodríguez-Vila

A method for incorporating three-dimensional residual stretches/stresses into patient-specific finite element simulations of arteries.

The existence of residual stresses in human arteries has long been shown experimentally. Researchers have also demonstrated that residual stresses have a significant effect on the distribution of physiological stresses within arterial tissues, and hence on their development, e.g., stress-modulated remodeling. Through progress in medical imaging, image analysis and finite element (FE) meshing tools it is now possible to construct in vivo patient-specific geometries and thus to study specific, clinically relevant problems in arterial mechanics via FE simulations. Classical continuum mechanics and FE methods assume that constitutive models and the corresponding simulations start from unloaded, stress-free reference configurations while the boundary-value problem of interest represents a loaded geometry and includes residual stresses. We present a pragmatic methodology to simultaneously account for both (i) the three-dimensional (3-D) residual stress distributions in the arterial tissue layers, and (ii) the equilibrium of the in vivo patient-specific geometry with the known boundary conditions. We base our methodology on analytically determined residual stress distributions (Holzapfel and Ogden, 2010, J. R. Soc. Interface 7, 787-799) and calibrate it using data on residual deformations (Holzapfel et al., 2007, Ann. Biomed. Eng. 35, 530-545). We demonstrate our methodology on three patient-specific FE simulations calibrated using experimental data. All data employed here are generated from human tissues - both the aorta and thrombus, and their respective layers - including the geometries determined from magnetic resonance images, and material properties and 3-D residual stretches determined from mechanical experiments. We study the effect of 3-D residual stresses on the distribution of physiological stresses in the aortic layers (intima, media, adventitia) and the layers of the intraluminal thrombus (luminal, medial, abluminal) by comparing three types of FE simulations: (i) conventional calculations; (ii) calculations accounting only for prestresses; (iii) calculations including both 3-D residual stresses and prestresses. Our results show that including residual stresses in patient-specific simulations of arterial tissues significantly impacts both the global (organ-level) deformations and the stress distributions within the arterial tissue (and its layers). Our method produces circumferential Cauchy stress distributions that are more uniform through the tissue thickness (i.e., smaller stress gradients in the local radial directions) compared to both the conventional and prestressing calculations. Such methods, combined with appropriate experimental data, aim at increasing the accuracy of classical FE analyses for patient-specific studies in computational biomechanics and may lead to increased clinical application of simulation tools.

Methodology for registration of distended rectums in pelvic CT studies

PURPOSE Accurate delineation of the rectum is of high importance in off-line adaptive radiation therapy since it is a major dose-limiting organ in prostate cancer radiotherapy. The intensity-based deformable image registration (DIR) methods cannot create a correct spatial transformation if there is no correspondence between the template and the target images. The variation of rectal filling, gas, or feces, creates a non correspondence in image intensities that becomes a great obstacle for intensity-based DIR. METHODS In this study the authors have designed and implemented a semiautomatic method to create a rectum mask in pelvic computed tomography (CT) images. The method, that includes a DIR based on the demons algorithm, has been tested in 13 prostate cancer cases, each comprising of two CT scans, for a total of 26 CT scans. RESULTS The use of the manual segmentation in the planning image and the proposed rectum mask method (RMM) method in the daily image leads to an improvement in the DIR performance in pelvic CT images, obtaining a mean value of overlap volume index = 0.89, close to the values obtained using the manual segmentations in both images. CONCLUSIONS The application of the RMM method in the daily image and the manual segmentations in the planning image during prostate cancer treatments increases the performance of the registration in presence of rectal fillings, obtaining very good agreement with a physicians manual contours.

3D deformable registration for monitoring radiotherapy treatment in prostate cancer

Two deformable registration methods, the Demons and the Morphon algorithms, have been used for registration of CT datasets to evaluate their usability in radiotherapy planning for prostate cancer. These methods were chosen because they can perform deformable registration in a fully automated way. The experiments show that for intrapatient registration both of the methods give useful results, although some differences exist in the way they deform the template. The Morphon method has, however, some advantageous compared to the Demons method. It is invariant to the image intensity and it does not distort the deformed data. The conclusion is therefore to recommend the Morphon method as a registration tool for this application. A more flexible regularization model is needed, though, in order to be able to catch the full range of deformations required to match the datasets.

Rupture risk in abdominal aortic aneurysms: A realistic assessment of the explicit GPU approach

Accurate estimation of peak wall stress (PWS) is the crux of biomechanically motivated rupture risk assessment for abdominal aortic aneurysms aimed to improve clinical outcomes. Such assessments often use the finite element (FE) method to obtain PWS, albeit at a high computational cost, motivating simplifications in material or element formulations. These simplifications, while useful, come at a cost of reliability and accuracy. We achieve research-standard accuracy and maintain clinically applicable speeds by using novel computational technologies. We present a solution using our custom finite element code based on graphics processing unit (GPU) technology that is able to account for added complexities involved with more physiologically relevant solutions, e.g. strong anisotropy and heterogeneity. We present solutions up to 17× faster relative to an established finite element code using state-of-the-art nonlinear, anisotropic and nearly-incompressible material descriptions. We show a realistic assessment of the explicit GPU FE approach by using complex problem geometry, biofidelic material law, double-precision floating point computation and full element integration. Due to the increased solution speed without loss of accuracy, shown on five clinical cases of abdominal aortic aneurysms, the method shows promise for clinical use in determining rupture risk of abdominal aortic aneurysms.

Automated Delineation of Vessel Wall and Thrombus Boundaries of Abdominal Aortic Aneurysms Using Multispectral MR Images.

A correct patient-specific identification of the abdominal aortic aneurysm is useful for both diagnosis and treatment stages, as it locates the disease and represents its geometry. The actual thickness and shape of the arterial wall and the intraluminal thrombus are of great importance when predicting the rupture of the abdominal aortic aneurysms. The authors describe a novel method for delineating both the internal and external contours of the aortic wall, which allows distinguishing between vessel wall and intraluminal thrombus. The method is based on active shape model and texture statistical information. The method was validated with eight MR patient studies. There was high correspondence between automatic and manual measurements for the vessel wall area. Resulting segmented images presented a mean Dice coefficient with respect to manual segmentations of 0.88 and a mean modified Hausdorff distance of 1.14 mm for the internal face and 0.86 and 1.33 mm for the external face of the arterial wall. Preliminary results of the segmentation show high correspondence between automatic and manual measurements for the vessel wall and thrombus areas. However, since the dataset is small the conclusions cannot be generalized.

Interpretation of motion analysis of laparoscopic instruments based on principal component analysis in box trainer settings

BackgroundMotion analysis parameters (MAPs) have been extensively validated for assessment of minimally invasive surgical skills. However, there are discrepancies on how specific MAPs, tasks, and skills match with each other, reflecting that motion analysis cannot be generalized independently of the learning outcomes of a task. Additionally, there is a lack of knowledge on the meaning of motion analysis in terms of surgical skills, making difficult the provision of meaningful, didactic feedback. In this study, new higher significance MAPs (HSMAPs) are proposed, validated, and discussed for the assessment of technical skills in box trainers, based on principal component analysis (PCA).MethodsMotion analysis data were collected from 25 volunteers performing three box trainer tasks (peg grasping/PG, pattern cutting/PC, knot suturing/KS) using the EVA tracking system. PCA was applied on 10 MAPs for each task and hand. Principal components were trimmed to those accounting for an explained variance > 80% to define the HSMAPs. Individual contributions of MAPs to HSMAPs were obtained by loading analysis and varimax rotation. Construct validity of the new HSMAPs was carried out at two levels of experience based on number of surgeries.ResultsThree new HSMAPs per hand were defined for PG and PC tasks, and two per hand for KS task. PG presented validity for HSMAPs related to insecurity and economy of space. PC showed validity for HSMAPs related to cutting efficacy, peripheral unawareness, and confidence. Finally, KS presented validity for HSMAPs related with economy of space and knotting security.ConclusionsPCA-defined HSMAPs can be used for technical skills’ assessment. Construct validation and expert knowledge can be combined to infer how competences are acquired in box trainer tasks. These findings can be exploited to provide residents with meaningful feedback on performance. Future works will compare the new HSMAPs with valid scoring systems such as GOALS.

Laparoscopic Video and Ultrasounds Image Processing in Minimally Invasive Pancreatic Surgeries

Due to limitations in conventional medical imaging and the restrictions imposed by both the anatomy and the surgical approach in pancreatic cancer, there is a need for methods to support intraoperative imaging in order to improve their accurate anatomical localization and the characterization of their nature. Laparoscopic ultrasounds (LUS) images and endoscopic videos can be used to extract useful information during the surgical procedures. A fast approach for acquiring an estimation of the tumor positioning and size through laparoscopic ultrasounds images has been developed. Based on the surgical video, endoscope 3D tracking is achieved by means of a Shape-from-Motion technique. Intraoperative imaging algorithms’ validation has been carried out in an ex vivo porcine model and results have shown the viability of exploiting them for structures characterization and their 3D reconstruction.

Quantitative evaluation of patient-specific conforming hexahedral meshes of abdominal aortic aneurysms and intraluminal thrombus generated from MRI

A novel method for generating patient-specific high quality conforming hexahedral meshes is presented. The meshes are directly obtained from the segmentation of patient magnetic resonance (MR) images of abdominal aortic aneurysms (AAA). The MRI permits distinguishing between structures of interest in soft tissue. Being so, the contours of the lumen, the aortic wall and the intraluminal thrombus (ILT) are available and thus the meshes represent the actual anatomy of the patient’s aneurysm, including the layered morphologies of these structures. Most AAAs are located in the lower part of the aorta and the upper section of the iliac arteries, where the inherent tortuosity of the anatomy and the presence of the ILT makes the generation of high-quality elements at the bifurcation is a challenging task. In this work we propose a novel approach for building quadrilateral meshes for each surface of the sectioned geometry, and generating conforming hexahedral meshes by combining the quadrilateral meshes. Conforming hexahedral meshes are created for the wall and the ILT. The resulting elements are evaluated on four patients’ datasets using the Scaled Jacobian metric. Hexahedral meshes of 25,000 elements with 94.8% of elements well-suited for FE analysis are generated.

Method for Incorporating Three-Dimensional Residual Stresses Into Patient-Specific Simulations of Arteries

METHOD FOR INCORPORATING THREE-DIMENSIONAL RESIDUAL STRESSES INTO PATIENT-SPECIFIC SIMULATIONS OF ARTERIES