Validation of statistical shape model based reconstruction of the proximal femur - A morphology study

Seventeen bones (sixteen cadaveric bones and one plastic bone) were used to validate a method for reconstructing a surface model of the proximal femur from 2D X-ray radiographs and a statistical shape model that was constructed from thirty training surface models. Unlike previously introduced validation studies, where surface-based distance errors were used to evaluate the reconstruction accuracy, here we propose to use errors measured based on clinically relevant morphometric parameters. For this purpose, a program was developed to robustly extract those morphometric parameters from the thirty training surface models (training population), from the seventeen surface models reconstructed from X-ray radiographs, and from the seventeen ground truth surface models obtained either by a CT-scan reconstruction method or by a laser-scan reconstruction method. A statistical analysis was then performed to classify the seventeen test bones into two categories: normal cases and outliers. This classification step depends on the measured parameters of the particular test bone. In case all parameters of a test bone were covered by the training population's parameter ranges, this bone is classified as normal bone, otherwise as outlier bone. Our experimental results showed that statistically there was no significant difference between the morphometric parameters extracted from the reconstructed surface models of the normal cases and those extracted from the reconstructed surface models of the outliers. Therefore, our statistical shape model based reconstruction technique can be used to reconstruct not only the surface model of a normal bone but also that of an outlier bone.

ANOVA kernels and RKHS of zero mean functions for model-based sensitivity analysis

Given a reproducing kernel Hilbert space (H,〈.,.〉)(H,〈.,.〉) of real-valued functions and a suitable measure μμ over the source space D⊂RD⊂R, we decompose HH as the sum of a subspace of centered functions for μμ and its orthogonal in HH. This decomposition leads to a special case of ANOVA kernels, for which the functional ANOVA representation of the best predictor can be elegantly derived, either in an interpolation or regularization framework. The proposed kernels appear to be particularly convenient for analyzing the effect of each (group of) variable(s) and computing sensitivity indices without recursivity.

Statistical model-based segmentation of the proximal femur in digital antero-posterior (AP) pelvic radiographs.

PURPOSE    Segmentation of the proximal femur in digital antero-posterior (AP) pelvic radiographs is required to create a three-dimensional model of the hip joint for use in planning and treatment. However, manually extracting the femoral contour is tedious and prone to subjective bias, while automatic segmentation must accommodate poor image quality, anatomical structure overlap, and femur deformity. A new method was developed for femur segmentation in AP pelvic radiographs. METHODS    Using manual annotations on 100 AP pelvic radiographs, a statistical shape model (SSM) and a statistical appearance model (SAM) of the femur contour were constructed. The SSM and SAM were used to segment new AP pelvic radiographs with a three-stage approach. At initialization, the mean SSM model is coarsely registered to the femur in the AP radiograph through a scaled rigid registration. Mahalanobis distance defined on the SAM is employed as the search criteria for each annotated suggested landmark location. Dynamic programming was used to eliminate ambiguities. After all landmarks are assigned, a regularized non-rigid registration method deforms the current mean shape of SSM to produce a new segmentation of proximal femur. The second and third stages are iteratively executed to convergence. RESULTS    A set of 100 clinical AP pelvic radiographs (not used for training) were evaluated. The mean segmentation error was [Formula: see text], requiring [Formula: see text] s per case when implemented with Matlab. The influence of the initialization on segmentation results was tested by six clinicians, demonstrating no significance difference. CONCLUSIONS    A fast, robust and accurate method for femur segmentation in digital AP pelvic radiographs was developed by combining SSM and SAM with dynamic programming. This method can be extended to segmentation of other bony structures such as the pelvis.

Early versions in medieval textual traditions. Wolfram’s Parzival as a test case

Wolfram von Eschenbach’s novel Parzival is a courtly romance composed in German language shortly after 1200. In a project, based at the University of Bern, a new critical edition of the poem is prepared in electronic and printed form. It visualizes parallel textual versions, which, depending on particular circumstances of oral performance, have developed in the early stage of the poem’s transmission. Philological research as well as phylogenetic techniques common in the natural sciences, e.g. in molecular biology, have been used to demonstrate the existence of these early textual versions. The article shows how both methods work and how they are applied to the ongoing edition. Exemplary passages to be presented include the text of some rare fragments written in the first decades of the 13th century, which might even go back to the author’s lifetime and which allow to date the existence of the versions they belong to.

Tambora 1815 as a test case for high impact volcanic eruptions: Earth system effects

The eruption of Tambora (Indonesia) in April 1815 had substantial effects on global climate and led to the ‘Year Without a Summer’ of 1816 in Europe and North America. Although a tragic event — tens of thousands of people lost their lives — the eruption also was an ‘experiment of nature’ from which science has learned until today. The aim of this study is to summarize our current understanding of the Tambora eruption and its effects on climate as expressed in early instrumental observations, climate proxies and geological evidence, climate reconstructions, and model simulations. Progress has been made with respect to our understanding of the eruption process and estimated amount of SO2 injected into the atmosphere, although large uncertainties still exist with respect to altitude and hemispheric distribution of Tambora aerosols. With respect to climate effects, the global and Northern Hemispheric cooling are well constrained by proxies whereas there is no strong signal in Southern Hemisphere proxies. Newly recovered early instrumental information for Western Europe and parts of North America, regions with particularly strong climate effects, allow Tambora’s effect on the weather systems to be addressed. Climate models respond to prescribed Tambora-like forcing with a strengthening of the wintertime stratospheric polar vortex, global cooling and a slowdown of the water cycle, weakening of the summer monsoon circulations, a strengthening of the Atlantic Meridional Overturning Circulation, and a decrease of atmospheric CO₂. Combining observations, climate proxies, and model simulations for the case of Tambora, a better understanding of climate processes has emerged.

Kernel density correlation based non-rigid point set matching for statistical model-based 2D/3D reconstruction

This paper presents a kernel density correlation based nonrigid point set matching method and shows its application in statistical model based 2D/3D reconstruction of a scaled, patient-specific model from an un-calibrated x-ray radiograph. In this method, both the reference point set and the floating point set are first represented using kernel density estimates. A correlation measure between these two kernel density estimates is then optimized to find a displacement field such that the floating point set is moved to the reference point set. Regularizations based on the overall deformation energy and the motion smoothness energy are used to constraint the displacement field for a robust point set matching. Incorporating this non-rigid point set matching method into a statistical model based 2D/3D reconstruction framework, we can reconstruct a scaled, patient-specific model from noisy edge points that are extracted directly from the x-ray radiograph by an edge detector. Our experiment conducted on datasets of two patients and six cadavers demonstrates a mean reconstruction error of 1.9 mm