Mesh-based vs. Image-based Statistical Appearance Model of the Human Femur: a Preliminary Comparison Study for the Creation of Finite Element Meshes

Statistical models have been recently introduced in computational orthopaedics to investigate the bone mechanical properties across several populations. A fundamental aspect for the construction of statistical models concerns the establishment of accurate anatomical correspondences among the objects of the training dataset. Various methods have been proposed to solve this problem such as mesh morphing or image registration algorithms. The objective of this study is to compare a mesh-based and an image-based statistical appearance model approaches for the creation of nite element(FE) meshes. A computer tomography (CT) dataset of 157 human left femurs was used for the comparison. For each approach, 30 finite element meshes were generated with the models. The quality of the obtained FE meshes was evaluated in terms of volume, size and shape of the elements. Results showed that the quality of the meshes obtained with the image-based approach was higher than the quality of the mesh-based approach. Future studies are required to evaluate the impact of this finding on the final mechanical simulations.

hp-dGFEM for Second-Order Elliptic Problems in Plyhedra I: Stability and Quasioptimality on Geometric Meshes

We introduce and analyze hp-version discontinuous Galerkin (dG) finite element methods for the numerical approximation of linear second-order elliptic boundary-value problems in three-dimensional polyhedral domains. To resolve possible corner-, edge- and corner-edge singularities, we consider hexahedral meshes that are geometrically and anisotropically refined toward the corresponding neighborhoods. Similarly, the local polynomial degrees are increased linearly and possibly anisotropically away from singularities. We design interior penalty hp-dG methods and prove that they are well-defined for problems with singular solutions and stable under the proposed hp-refinements. We establish (abstract) error bounds that will allow us to prove exponential rates of convergence in the second part of this work.

Getting in touch-3D printing in Forensic Imaging

With the increasing use of medical imaging in forensics, as well as the technological advances in rapid prototyping, we suggest combining these techniques to generate displays of forensic findings. We used computed tomography (CT), CT angiography, magnetic resonance imaging (MRI) and surface scanning with photogrammetry in conjunction with segmentation techniques to generate 3D polygon meshes. Based on these data sets, a 3D printer created colored models of the anatomical structures. Using this technique, we could create models of bone fractures, vessels, cardiac infarctions, ruptured organs as well as bitemark wounds. The final models are anatomically accurate, fully colored representations of bones, vessels and soft tissue, and they demonstrate radiologically visible pathologies. The models are more easily understood by laypersons than volume rendering or 2D reconstructions. Therefore, they are suitable for presentations in courtrooms and for educational purposes.

Influence of Smoothing on Voxel-Based Mesh Accuracy in Micro-Finite Element

The interest in automatic volume meshing for finite element analysis (FEA) has grown more since the appearance of microfocus CT (μCT), due to its high resolution, which allows for the assessment of mechanical behaviour at a high precision. Nevertheless, the basic meshing approach of generating one hexahedron per voxel produces jagged edges. To prevent this effect, smoothing algorithms have been introduced to enhance the topology of the mesh. However, whether smoothing also improves the accuracy of voxel-based meshes in clinical applications is still under question. There is a trade-off between smoothing and quality of elements in the mesh. Distorted elements may be produced by excessive smoothing and reduce accuracy of the mesh. In the present work, influence of smoothing on the accuracy of voxel-based meshes in micro-FE was assessed. An accurate 3D model of a trabecular structure with known apparent mechanical properties was used as a reference model. Virtual CT scans of this reference model (with resolutions of 16, 32 and 64 μm) were then created and used to build voxel-based meshes of the microarchitecture. Effects of smoothing on the apparent mechanical properties of the voxel-based meshes as compared to the reference model were evaluated. Apparent Young’s moduli of the smooth voxel-based mesh were significantly closer to those of the reference model for the 16 and 32 μm resolutions. Improvements were not significant for the 64 μm, due to loss of trabecular connectivity in the model. This study shows that smoothing offers a real benefit to voxel-based meshes used in micro-FE. It might also broaden voxel-based meshing to other biomechanical domains where it was not used previously due to lack of accuracy. As an example, this work will be used in the framework of the European project ContraCancrum, which aims at providing a patient-specific simulation of tumour development in brain and lungs for oncologists. For this type of clinical application, such a fast, automatic, and accurate generation of the mesh is of great benefit.

Feature-invariant image registration method for quantification of surgical outcomes in patients with craniosynostosis: a preliminary study

Craniosynostosis consists of a premature fusion of the sutures in an infant skull that restricts skull and brain growth. During the last decades, there has been a rapid increase of fundamentally diverse surgical treatment methods. At present, the surgical outcome has been assessed using global variables such as cephalic index, head circumference, and intracranial volume. However, these variables have failed in describing the local deformations and morphological changes that may have a role in the neurologic disorders observed in the patients. This report describes a rigid image registration-based method to evaluate outcomes of craniosynostosis surgical treatments, local quantification of head growth, and indirect intracranial volume change measurements. The developed semiautomatic analysis method was applied to computed tomography data sets of a 5-month-old boy with sagittal craniosynostosis who underwent expansion of the posterior skull with cranioplasty. Quantification of the local changes between pre- and postoperative images was quantified by mapping the minimum distance of individual points from the preoperative to the postoperative surface meshes, and indirect intracranial volume changes were estimated. The proposed methodology can provide the surgeon a tool for the quantitative evaluation of surgical procedures and detection of abnormalities of the infant skull and its development.

The examination and identification of bite marks in foods using 3D scanning and 3D comparison methods

Bite mark analysis offers the opportunity to identify the biter based on the individual characteristics of the dentitions. Normally, the main focus is on analysing bite mark injuries on human bodies, but also, bite marks in food may play an important role in the forensic investigation of a crime. This study presents a comparison of simulated bite marks in different kinds of food with the dentitions of the presumed biter. Bite marks were produced by six adults in slices of buttered bread, apples, different kinds of Swiss chocolate and Swiss cheese. The time-lapse influence of the bite mark in food, under room temperature conditions, was also examined. For the documentation of the bite marks and the dentitions of the biters, 3D optical surface scanning technology was used. The comparison was performed using two different software packages: the ATOS modelling and analysing software and the 3D studio max animation software. The ATOS software enables an automatic computation of the deviation between the two meshes. In the present study, the bite marks and the dentitions were compared, as well as the meshes of each bite mark which were recorded in the different stages of time lapse. In the 3D studio max software, the act of biting was animated to compare the dentitions with the bite mark. The examined food recorded the individual characteristics of the dentitions very well. In all cases, the biter could be identified, and the dentitions of the other presumed biters could be excluded. The influence of the time lapse on the food depends on the kind of food and is shown on the diagrams. However, the identification of the biter could still be performed after a period of time, based on the recorded individual characteristics of the dentitions.

Cartilage tissue engineering by expanded goat articular chondrocytes

In this study we investigated whether expanded goat chondrocytes have the capacity to generate cartilaginous tissues with biochemical and biomechanical properties improving with time in culture. Goat chondrocytes were expanded in monolayer with or without combinations of FGF-2, TGF-beta1, and PDGFbb, and the postexpansion chondrogenic capacity assessed in pellet cultures. Expanded chondrocytes were also cultured for up to 6 weeks in HYAFF-M nonwoven meshes or Polyactive foams, and the resulting cartilaginous tissues were assessed histologically, biochemically, and biomechanically. Supplementation of the expansion medium with FGF-2 increased the proliferation rate of goat chondrocytes and enhanced their postexpansion chondrogenic capacity. FGF-2-expanded chondrocytes seeded in HYAFF-M or Polyactive scaffolds formed cartilaginous tissues with wet weight, glycosaminoglycan, and collagen content, increasing from 2 days to 6 weeks culture (up to respectively 2-, 8-, and 41-fold). Equilibrium and dynamic stiffness measured in HYAFF M-based constructs also increased with time, up to, respectively, 1.3- and 16-fold. This study demonstrates the feasibility to engineer goat cartilaginous tissues at different stages of development by varying culture time, and thus opens the possibility to test the effect of maturation stage of engineered cartilage on the outcome of cartilage repair in orthotopic goat models.

Laparoscopic ventral hernia repair is safe and cost effective

BACKGROUND: Ventral hernia repair is increasingly performed by laparoscopic means since the introduction of dual-layer meshes. This study aimed to compare the early complications and cost effectiveness of open hernia repair with those associated with laparoscopic repair. METHODS: Open ventral hernia repair was performed for 92 consecutive patients using a Vypro mesh, followed by laparoscopic repair for 49 consecutive patients using a Parietene composite mesh. RESULTS: The rate of surgical-site infections was significantly higher with open ventral hernia repair (13 vs 1; p = 0.03). The median length of hospital stay was significantly shorter with laparoscopic surgery (7 vs 6 days; p = 0.02). For laparoscopic repair, the direct operative costs were higher (2,314 vs 2,853 euros; p = 0.03), and the overall hospital costs were lower (9,787 vs 7,654 euros; p = 0.02). CONCLUSIONS: Laparoscopic ventral hernia repair leads to fewer surgical-site infections and a shorter hospital stay than open repair. Despite increased operative costs, overall hospital costs are lowered by laparoscopic ventral hernia repair.

hp-dGFEM for Second-Order Elliptic Problems in Plyhedra II: Exponential Convergence

The goal of this paper is to establish exponential convergence of $hp$-version interior penalty (IP) discontinuous Galerkin (dG) finite element methods for the numerical approximation of linear second-order elliptic boundary-value problems with homogeneous Dirichlet boundary conditions and piecewise analytic data in three-dimensional polyhedral domains. More precisely, we shall analyze the convergence of the $hp$-IP dG methods considered in [D. Schötzau, C. Schwab, T. P. Wihler, SIAM J. Numer. Anal., 51 (2013), pp. 1610--1633] based on axiparallel $\sigma$-geometric anisotropic meshes and $\bm{s}$-linear anisotropic polynomial degree distributions.

Intussusceptive microvascular growth: a common alternative to capillary sprouting.

Intussusceptive capillary growth represents a new principle for microvascular growth as described in the lungs of growing rats. According to this concept, the capillary network expands by the formation of slender transcapillary tissue pillars, which give rise to new vascular meshes. The process was first observed in Mercox casts of the lung microvasculature, which revealed the existence of multiple tiny holes with diameters around 1.5 microns. Consecutive transmission electron microscopic investigation of serial sections demonstrated that the holes corresponded to slender tissue pillars (Burri and Tarek, 1990). The corrosion cast technique thus appears to be an adequate screening method for intussusceptive growth. In the present investigation, Mercox casts of various vascular systems, namely, those of the eye, submandibular gland, heart, liver, stomach, small and large intestine, trachea, kidney, uterus and ovary were prepared from rats aged between 4 and 9 weeks in order to screen them for the existence of the typical tiny holes representing tissue pillars. In all organs investigated, these structures were observed in various locations to a variable degree. They were mainly encountered within dilated vascular segments or at triple or quadruple branching points of the circulation. Even in capillary networks with a three-dimensional arrangement could these pillars be detected. Intussusception thus appears to be a principle of growth appertaining to many vascular systems.

Scanning electron microscope study of the developing microvasculature in the postnatal rat lung.

During postnatal growth the parenchymal septa of rat lung undergo an impressive restructuring. While immature septa are thick and contain two capillary layers, mature septa are slender and contain a single microvascular network. Using the Mercox casting technique and scanning electron microscopy, we investigated the mode and the timing of the transformation of the pulmonary capillary bed. During the third postnatal week the parenchymal septa rapidly mature to match adult morphology. Even in adult lungs, however, remnants of the immature status are present: A capillary bilayer is regularly found at the base and the tip of the septa. Our observations support the concept that reduction of intervening tissue, partial fusion of the two capillary networks, and preferential growth lead to the mature vascular arrangement. The fact that true mature interalveolar septa show a denser capillary network than alveolar walls abutting onto pleura, bronchi, or larger vessels is consonant with the fusion theory. Towards the nonparenchyma, the capillary network surrounding every airspace had no counterpart to fuse with. From quantitative data it can be calculated that owing to lung growth, mesh size should increase more than four times between birth and adult age. The adult lung network, however, is denser than the one in young animals. This means that new meshes must be added during growth. We propose that small holes observed in sheet-like regions of the microvasculature enlarge to form new capillary meshes. With this mechanism of in-itself or intussusceptional growth, sprouting of individual capillary segments to increase network size is no longer needed.

An image-based method to automatically propagate bony landmarks: application to computational spine biomechanics.

In attempts to elucidate the underlying mechanisms of spinal injuries and spinal deformities, several experimental and numerical studies have been conducted to understand the biomechanical behavior of the spine. However, numerical biomechanical studies suffer from uncertainties associated with hard- and soft-tissue anatomies. Currently, these parameters are identified manually on each mesh model prior to simulations. The determination of soft connective tissues on finite element meshes can be a tedious procedure, which limits the number of models used in the numerical studies to a few instances. In order to address these limitations, an image-based method for automatic morphing of soft connective tissues has been proposed. Results showed that the proposed method is capable to accurately determine the spatial locations of predetermined bony landmarks. The present method can be used to automatically generate patient-specific models, which may be helpful in designing studies involving a large number of instances and to understand the mechanical behavior of biomechanical structures across a given population.

Image-based vs. mesh-based statistical appearance models of the human femur: Implications for finite element simulations.

Statistical appearance models have recently been introduced in bone mechanics to investigate bone geometry and mechanical properties in population studies. The establishment of accurate anatomical correspondences is a critical aspect for the construction of reliable models. Depending on the representation of a bone as an image or a mesh, correspondences are detected using image registration or mesh morphing. The objective of this study was to compare image-based and mesh-based statistical appearance models of the femur for finite element (FE) simulations. To this aim, (i) we compared correspondence detection methods on bone surface and in bone volume; (ii) we created an image-based and a mesh-based statistical appearance models from 130 images, which we validated using compactness, representation and generalization, and we analyzed the FE results on 50 recreated bones vs. original bones; (iii) we created 1000 new instances, and we compared the quality of the FE meshes. Results showed that the image-based approach was more accurate in volume correspondence detection and quality of FE meshes, whereas the mesh-based approach was more accurate for surface correspondence detection and model compactness. Based on our results, we recommend the use of image-based statistical appearance models for FE simulations of the femur.

Fluid-structure interaction of an elastic pipe immersed inside a coaxial rigid pipe: a model for the Tullio phenomenon

An axisymmetric, elastic pipe is filled with an incompressible fluid and is immersed in a second, coaxial rigid pipe which contains the same fluid. A pressure pulse in the outer fluid annulus deforms the elastic pipe which invokes a fluid motion in the fluid core. It is the aim of this study to investigate streaming phenomena in the core which may originate from such a fluid-structure interaction. This work presents a numerical solver for such a configuration. It was developed in the OpenFOAM software environment and is based on the Arbitrary Lagrangian Eulerian (ALE) approach for moving meshes. The solver features a monolithic integration of the one-dimensional, coupled system between the elastic structure and the outer fluid annulus into a dynamic boundary condition for the moving surface of the fluid core. Results indicate that our configuration may serve as a mechanical model of the Tullio Phenomenon (sound-induced vertigo).