Philippe Büchler

The virtual skeleton database: an open access repository for biomedical research and collaboration.

BACKGROUND Statistical shape models are widely used in biomedical research. They are routinely implemented for automatic image segmentation or object identification in medical images. In these fields, however, the acquisition of the large training datasets, required to develop these models, is usually a time-consuming process. Even after this effort, the collections of datasets are often lost or mishandled resulting in replication of work. OBJECTIVE To solve these problems, the Virtual Skeleton Database (VSD) is proposed as a centralized storage system where the data necessary to build statistical shape models can be stored and shared. METHODS The VSD provides an online repository system tailored to the needs of the medical research community. The processing of the most common image file types, a statistical shape model framework, and an ontology-based search provide the generic tools to store, exchange, and retrieve digital medical datasets. The hosted data are accessible to the community, and collaborative research catalyzes their productivity. RESULTS To illustrate the need for an online repository for medical research, three exemplary projects of the VSD are presented: (1) an international collaboration to achieve improvement in cochlear surgery and implant optimization, (2) a population-based analysis of femoral fracture risk between genders, and (3) an online application developed for the evaluation and comparison of the segmentation of brain tumors. CONCLUSIONS The VSD is a novel system for scientific collaboration for the medical image community with a data-centric concept and semantically driven search option for anatomical structures. The repository has been proven to be a useful tool for collaborative model building, as a resource for biomechanical population studies, or to enhance segmentation algorithms.

Biomechanical model of human cornea based on stromal microstructure.

The optical characteristics of the human cornea depends on the mechanical balance between the intra-ocular pressure and intrinsic tissue stiffness. A wide range of ophthalmic surgical procedures alter corneal biomechanics to induce local or global curvature changes for the correction of visual acuity. Due to the large number of surgical interventions performed every day, a deeper understanding of corneal biomechanics is needed to improve the safety of these procedures and medical devices. The aim of this study is to propose a biomechanical model of the human cornea, based on stromal microstructure. The constitutive mechanical law includes collagen fiber distribution based on X-ray scattering analysis, collagen cross-linking, and fiber uncrimping. Our results showed that the proposed model reproduced inflation and extensiometry experimental data [Elsheikh et al., Curr. Eye Res., 2007; Elsheikh et al., Exp. Eye Res., 2008] successfully. The mechanical properties obtained for different age groups demonstrated an increase in collagen cross-linking for older specimens. In future work such a model could be used to simulate non-symmetric interventions, and provide better surgical planning.

Optimisation of orthopaedic implant design using statistical shape space analysis based on level sets.

Statistical shape analysis techniques have shown to be efficient tools to build population specific models of anatomical variability. Their use is commonplace as prior models for segmentation, in which case the instance from the shape model that best fits the image data is sought. In certain cases, however, it is not just the most likely instance that must be searched, but rather the whole set of shape instances that meet certain criterion. In this paper we develop a method for the assessment of specific anatomical/morphological criteria across the shape variability found in a population. The method is based on a level set segmentation approach, and used on the parametric space of the statistical shape model of the target population, solved via a multi-level narrow-band approach for computational efficiency. Based on this technique, we develop a framework for evidence-based orthopaedic implant design. To date, implants are commonly designed and validated by evaluating implant bone fitting on a limited set of cadaver bones, which not necessarily span the whole variability in the population. Based on our framework, we can virtually fit a proposed implant design to samples drawn from the statistical model, and assess which range of the population is suitable for the implant. The method highlights which patterns of bone variability are more important for implant fitting, allowing and easing implant design improvements, as to fit a maximum of the target population. Results are presented for the optimisation of implant design of proximal human tibia, used for internal fracture fixation.

Statistical finite element model for bone shape and biomechanical properties

We present a framework for statistical finite element analysis combining shape and material properties, and allowing performing statistical statements of biomechanical performance across a given population. In this paper, we focus on the design of orthopaedic implants that fit a maximum percentage of the target population, both in terms of geometry and biomechanical stability. CT scans of the bone under consideration are registered non-rigidly to obtain correspondences in position and intensity between them. A statistical model of shape and intensity (bone density) is computed by means of principal component analysis. Afterwards, finite element analysis (FEA) is performed to analyse the biomechanical performance of the bones. Realistic forces are applied on the bones and the resulting displacement and bone stress distribution are calculated. The mechanical behaviour of different PCA bone instances is compared.

Multiscale Modeling for Image Analysis of Brain Tumor Studies

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multiscale, multiphysics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlas-based segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

Oxygen and glucose distribution after intracorneal lens implantation.

Purpose. Insertion of an implant in the cornea to achieve corneal multifocality has been suggested as a solution for presbyopia. However, unresolved issues related to nutrient transport need to be resolved. Our aim was to find the best lens position and influence lens transport properties in order to optimize nutrient supply to corneal cells. Method. An axisymmetric corneal model was built to simulate the nutrient transport in the cornea. Oxygen and glucose concentrations were calculated for normal cornea and intracorneal lens wearing conditions. The simulation considers the different tissue layers (epithelium, stroma, and endothelium) as well as layer and solute concentration dependent consumption. Results. The minimum oxygen tension in the cornea was found to be higher when the lens was placed at 3/4 of the corneal thickness. Moreover, in this position, the influence of the inlay diffusivity was smaller than at more anterior or posterior placements. The diffusivity of the inlay affects the way nutrients will be transported through the cornea. The threshold where glucose may diffuse through or around the implant was found to be 1/100th of the stromal diffusivity. Conclusions. Computational methods are especially attractive to study nutrient transport in the cornea due to the difficulties associated with in vivo or in vitro measurements. The exact parameters that dictate the corneal metabolism are not known. However, the combined analysis of oxygen and glucose distribution is valuable in order to predict the complex physiological changes that arise under intracorneal lens implantation.

Stress free configuration of the human eye

Numerical simulations of eye globes often rely on topographies that have been measured in vivo using devices such as the Pentacam or OCT. The topographies, which represent the form of the already stressed eye under the existing intraocular pressure, introduce approximations in the analysis. The accuracy of the simulations could be improved if either the stress state of the eye under the effect of intraocular pressure is determined, or the stress-free form of the eye estimated prior to conducting the analysis. This study reviews earlier attempts to address this problem and assesses the performance of an iterative technique proposed by Pandolfi and Holzapfel [1], which is both simple to implement and promises high accuracy in estimating the eyes stress-free form. A parametric study has been conducted and demonstrated reliance of the error level on the level of flexibility of the eye model, especially in the cornea region. However, in all cases considered 3-4 analysis iterations were sufficient to produce a stress-free form with average errors in node location <10(-6)mm and a maximal error <10(-4)mm. This error level, which is similar to what has been achieved with other methods and orders of magnitude lower than the accuracy of current clinical topography systems, justifies the use of the technique as a pre-processing step in ocular numerical simulations.

Specimen specific parameter identification of ovine lumbar intervertebral discs: On the influence of fibre–matrix and fibre–fibre shear interactions

Numerical models of the intervertebral disc, which address mechanical questions commonly make use of the difference in water content between annulus and nucleus, and thus fluid and solid parts are separated. Despite this simplification, models remain complex due to the anisotropy and nonlinearity of the annulus and regional variations of the collagen fibre density. Additionally, it has been shown that cross-links make a large contribution to the stiffness of the annulus. Because of this complex composite structure, it is difficult to reproduce several sets of experimental data with one single set of material parameters. This study addresses the question to which extent the ultrastructure of the intervertebral disc should be modelled so that its moment-angle behaviour can be adequately described. Therefore, a hyperelastic constitutive law, based on continuum mechanical principles was derived, which does not only consider the anisotropy from the collagen fibres, but also interactions among the fibres and between the fibres and the ground substance. Eight ovine lumbar intervertebral discs were tested on a custom made spinal loading simulator in flexion/extension, lateral bending and axial rotation. Specimen-specific geometrical models were generated using CT images and T2 maps to distinguish between annulus fibrosus and nucleus pulposus. For the identification of the material parameters the annulus fibrosus was described with two scenarios: with and without fibre-matrix and fibre-fibre interactions. Both scenarios showed a similar behaviour on a load displacement level. Comparing model predictions to the experimental data, the mean RMS of all specimens and all load cases was 0.54±0.15° without the interaction and 0.54±0.19° when the fibre-matrix and fibre-fibre interactions were included. However, due to the increased stiffness when cross-links effects were included, this scenario showed more physiological stress-strain relations in uniaxial and biaxial stress states. Thus, the present study suggests that fibre-matrix and fibre-fibre interactions should be considered in the constitutive law when the model addresses questions concerning the stress field of the annulus fibrosus.

Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

Artificial Muscle Devices: Innovations and Prospects for Fecal Incontinence Treatment

Fecal incontinence describes the involuntary loss of bowel content, which is responsible for stigmatization and social exclusion. It affects about 45% of retirement home residents and overall more than 12% of the adult population. Severe fecal incontinence can be treated by the implantation of an artificial sphincter. Currently available implants, however, are not part of everyday surgery due to long-term re-operation rates of 95% and definitive explantation rates of 40%. Such figures suggest that the implants fail to reproduce the capabilities of the natural sphincter. This article reviews the artificial sphincters on the market and under development, presents their physical principles of operation and critically analyzes their performance. We highlight the geometrical and mechanical parameters crucial for the design of an artificial fecal sphincter and propose more advanced mechanisms of action for a biomimetic device with sensory feedback. Dielectric electro-active polymer actuators are especially attractive because of their versatility, response time, reaction forces, and energy consumption. The availability of such technology will enable fast pressure adaption comparable to the natural feedback mechanism, so that tissue atrophy and erosion can be avoided while maintaining continence during daily activities.