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

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.

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.

Long-term effects of watershed management on surface runoff and sediment loss in the Ethiopian Highlands

Data on rainfall, runoff and sediment loss from different land use types have been collected by the Soil Conservation Research Programme in seven small catchments (73-673 hectares) throughout the Ethiopian Highlands since the early 1980s. Monitoring was carried out on a storm-to-storm basis for extended periods of 10-20 years, and the data are analysed here to assess long-term effects of changes. Soil and water conservation technologies were introduced in the early years in the catchments in view of their capacity to reduce runoff and sediment yield. Results indicate that rainfall did not substantially change over the observation periods. Land use changes and land degradation, however, altered runoff, as shown by the data from small test plots (30 m2), which were not altered by conservation measures during the monitoring periods. Sediment delivery from the catchments may have decreased due to soil and water conservation, while runoff rates did not change significantly. Extrapolation of the results in the highlands, however, showed that expansion of cultivated and grazing land induced by population growth may have increased the overall surface runoff. Watershed management in the catchments, finally, had beneficial effects on ecosystem services by reducing soil erosion, restoring soil fertility, enhancing agricultural production, and maintaining overall runoff to the benefit of lowland areas and neighbouring countries.

A system for 3-D reconstruction of a patient-specific surface model from calibrated X-ray images

This paper presents a system for 3-D reconstruction of a patient-specific surface model from calibrated X-ray images. Our system requires two X-ray images of a patient with one acquired from the anterior-posterior direction and the other from the axial direction. A custom-designed cage is utilized in our system to calibrate both images. Starting from bone contours that are interactively identified from the X-ray images, our system constructs a patient-specific surface model of the proximal femur based on a statistical model based 2D/3D reconstruction algorithm. In this paper, we present the design and validation of the system with 25 bones. An average reconstruction error of 0.95 mm was observed.

Ecosystem functions are resistant to extreme changes to rainfall regimes in a mesotrophic grassland

Major changes to rainfall regimes are predicted for the future but the effect of such changes on terrestrial ecosystem function is largely unknown. We created a rainfall manipulation experiment to investigate the effects of extreme changes in rainfall regimes on ecosystem functioning in a grassland system. We applied two rainfall regimes; a prolonged drought treatment (30 % reduction over spring and summer) and drought/downpour treatment (long periods of no rainfall interspersed with downpours), with an ambient control. Both rainfall manipulations included increased winter rainfall. We measured plant community composition, CO2 fluxes and soil nutrient availability. Plant species richness and cover were lower in the drought/downpour treatment, and showed little recovery after the treatment ceased. Ecosystem processes were less affected, possibly due to winter rainfall additions buffering reduced summer rainfall, which saw relatively small soil moisture changes. However, soil extractable P and ecosystem respiration were significantly higher in rainfall change treatments than in the control. This grassland appears fairly resistant, in the short term, to even the more extreme rainfall changes that are predicted for the region, although prolonged study is needed to measure longer-term impacts. Differences in ecosystem responses between the two treatments emphasise the variety of ecosystem responses to changes in both the size and frequency of rainfall events. Given that model predictions are inconsistent there is therefore a need to assess ecosystem function under a range of potential climate change scenarios.