An integrated approach for reconstructing surface models of the proximal femur from sparse input data for surgical navigation

A patient-specific surface model of the proximal femur plays an important role in planning and supporting various computer-assisted surgical procedures including total hip replacement, hip resurfacing, and osteotomy of the proximal femur. The common approach to derive 3D models of the proximal femur is to use imaging techniques such as computed tomography (CT) or magnetic resonance imaging (MRI). However, the high logistic effort, the extra radiation (CT-imaging), and the large quantity of data to be acquired and processed make them less functional. In this paper, we present an integrated approach using a multi-level point distribution model (ML-PDM) to reconstruct a patient-specific model of the proximal femur from intra-operatively available sparse data. Results of experiments performed on dry cadaveric bones using dozens of 3D points are presented, as well as experiments using a limited number of 2D X-ray images, which demonstrate promising accuracy of the present approach.

Evaluating the Quality and Accuracy of TanDEM-X Digital Elevation Models at Archaeological Sites in the Cilician Plain, Turkey

Satellite remote sensing provides a powerful instrument for mapping and monitoring traces of historical settlements and infrastructure, not only in distant areas and crisis regions. It helps archaeologists to embed their findings from field surveys into the broader context of the landscape. With the start of the TanDEM-X mission, spatially explicit 3D-information is available to researchers at an unprecedented resolution worldwide. We examined different experimental TanDEM-X digital elevation models (DEM) that were processed from two different imaging modes (Stripmap/High Resolution Spotlight) using the operational alternating bistatic acquisition mode. The quality and accuracy of the experimental DEM products was compared to other available DEM products and a high precision archaeological field survey. The results indicate the potential of TanDEM-X Stripmap (SM) data for mapping surface elements at regional scale. For the alluvial plain of Cilicia, a suspected palaeochannel could be reconstructed. At the local scale, DEM products from TanDEM-X High Resolution Spotlight (HS) mode were processed at 2 m spatial resolution using a merge of two monostatic/bistatic interferograms. The absolute and relative vertical accuracy of the outcome meet the specification of high resolution elevation data (HRE) standards from the National System for Geospatial Intelligence (NSG) at the HRE20 level.

Sevoflurane remifentanil interaction: comparison of different response surface models

Various pharmacodynamic response surface models have been developed to quantitatively describe the relationship between two or more drug concentrations with their combined clinical effect. We examined the interaction of remifentanil and sevoflurane on the probability of tolerance to shake and shout, tetanic stimulation, laryngeal mask airway insertion, and laryngoscopy in patients to compare the performance of five different response surface models.

Evaluation of 3-dimensional superimposition techniques on various skeletal structures of the head using surface models.

OBJECTIVES To test the applicability, accuracy, precision, and reproducibility of various 3D superimposition techniques for radiographic data, transformed to triangulated surface data. METHODS Five superimposition techniques (3P: three-point registration; AC: anterior cranial base; AC + F: anterior cranial base + foramen magnum; BZ: both zygomatic arches; 1Z: one zygomatic arch) were tested using eight pairs of pre-existing CT data (pre- and post-treatment). These were obtained from non-growing orthodontic patients treated with rapid maxillary expansion. All datasets were superimposed by three operators independently, who repeated the whole procedure one month later. Accuracy was assessed by the distance (D) between superimposed datasets on three form-stable anatomical areas, located on the anterior cranial base and the foramen magnum. Precision and reproducibility were assessed using the distances between models at four specific landmarks. Non parametric multivariate models and Bland-Altman difference plots were used for analyses. RESULTS There was no difference among operators or between time points on the accuracy of each superimposition technique (p>0.05). The AC + F technique was the most accurate (D<0.17 mm), as expected, followed by AC and BZ superimpositions that presented similar level of accuracy (D<0.5 mm). 3P and 1Z were the least accurate superimpositions (0.790.05), the detected structural changes differed significantly between different techniques (p<0.05). Bland-Altman difference plots showed that BZ superimposition was comparable to AC, though it presented slightly higher random error. CONCLUSIONS Superimposition of 3D datasets using surface models created from voxel data can provide accurate, precise, and reproducible results, offering also high efficiency and increased post-processing capabilities. In the present study population, the BZ superimposition was comparable to AC, with the added advantage of being applicable to scans with a smaller field of view.

Achieving a more realistic assessment of rockfall hazards by coupling three-dimensional process models and field-based tree-ring data

Sound knowledge of the spatial and temporal patterns of rockfalls is fundamental for the management of this very common hazard in mountain environments. Process-based, three-dimensional simulation models are nowadays capable of reproducing the spatial distribution of rockfall occurrences with reasonable accuracy through the simulation of numerous individual trajectories on highly-resolved digital terrain models. At the same time, however, simulation models typically fail to quantify the ‘real’ frequency of rockfalls (in terms of return intervals). The analysis of impact scars on trees, in contrast, yields real rockfall frequencies, but trees may not be present at the location of interest and rare trajectories may not necessarily be captured due to the limited age of forest stands. In this article, we demonstrate that the coupling of modeling with tree-ring techniques may overcome the limitations inherent to both approaches. Based on the analysis of 64 cells (40 m × 40 m) of a rockfall slope located above a 1631-m long road section in the Swiss Alps, we illustrate results from 488 rockfalls detected in 1260 trees. We illustrate that tree impact data cannot only be used (i) to reconstruct the real frequency of rockfalls for individual cells, but that they also serve (ii) the calibration of the rockfall model Rockyfor3D, as well as (iii) the transformation of simulated trajectories into real frequencies. Calibrated simulation results are in good agreement with real rockfall frequencies and exhibit significant differences in rockfall activity between the cells (zones) along the road section. Real frequencies, expressed as rock passages per meter road section, also enable quantification and direct comparison of the hazard potential between the zones. The contribution provides an approach for hazard zoning procedures that complements traditional methods with a quantification of rockfall frequencies in terms of return intervals through a systematic inclusion of impact records in trees.

Reconstruction of 3D Vertebral Models from a Single 2D Lateral Fluoroscopic Image

Accurate three-dimensional (3D) models of lumbar vertebrae are required for image-based 3D kinematics analysis. MRI or CT datasets are frequently used to derive 3D models but have the disadvantages that they are expensive, time-consuming or involving ionizing radiation (e.g., CT acquisition). In this chapter, we present an alternative technique that can reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image and a statistical shape model. Cadaveric studies are conducted to verify the reconstruction accuracy by comparing the surface models reconstructed from a single lateral fluoroscopic image to the ground truth data from 3D CT segmentation. A mean reconstruction error between 0.7 and 1.4 mm was found.

An ice-rich flow origin for the banded terrain in the Hellas basin, Mars

The interior of Hellas Basin displays a complex landscape and a variety of geomorphological domains. One of these domains, the enigmatic banded terrain covers much of the northwestern part of the basin. We use high-resolution (CTX and HiRISE) Digital Terrain Models to show that most of the complex viscous flowing behavior exhibited by the banded terrain is controlled by topography and flow-like interactions between neighboring banded terrain. Furthermore, the interior of the basin hosts several landforms suggestive of the presence of near-surface ice, which include polygonal patterns with elongated pits, scalloped depressions, isolated mounds and collapse structures. We suggest that thermal contraction cracking and sublimation of near-surface ice are responsible for the formation and the development of most of the ice-related landforms documented in Hellas. The relatively pristine form, lack of superposed craters, and strong association with the banded terrain, suggest an Amazonian (<3 Ga) age of formation for these landforms. Finally, relatively high surface pressures (above the triple point of water) expected in Hellas and summer-time temperatures often exceeding the melting point of water ice suggest that the basin may have recorded relatively “temperate” climatic conditions compared to other places on Mars. Therefore, the potentially ice-rich banded terrain may have deformed with lower viscosity and stresses compared to other locations on Mars, which may account for its unique morphology.

Accuracy of linear measurements using three imaging modalities: two lateral cephalograms and one 3D model from CBCT data

BACKGROUND The aim of this study was to evaluate the accuracy of linear measurements on three imaging modalities: lateral cephalograms from a cephalometric machine with a 3 m source-to-mid-sagittal-plane distance (SMD), from a machine with 1.5 m SMD and 3D models from cone-beam computed tomography (CBCT) data. METHODS Twenty-one dry human skulls were used. Lateral cephalograms were taken, using two cephalometric devices: one with a 3 m SMD and one with a 1.5 m SMD. CBCT scans were taken by 3D Accuitomo® 170, and 3D surface models were created in Maxilim® software. Thirteen linear measurements were completed twice by two observers with a 4 week interval. Direct physical measurements by a digital calliper were defined as the gold standard. Statistical analysis was performed. RESULTS Nasion-Point A was significantly different from the gold standard in all methods. More statistically significant differences were found on the measurements of the 3 m SMD cephalograms in comparison to the other methods. Intra- and inter-observer agreement based on 3D measurements was slightly better than others. LIMITATIONS Dry human skulls without soft tissues were used. Therefore, the results have to be interpreted with caution, as they do not fully represent clinical conditions. CONCLUSIONS 3D measurements resulted in a better observer agreement. The accuracy of the measurements based on CBCT and 1.5 m SMD cephalogram was better than a 3 m SMD cephalogram. These findings demonstrated the linear measurements accuracy and reliability of 3D measurements based on CBCT data when compared to 2D techniques. Future studies should focus on the implementation of 3D cephalometry in clinical practice.

Patterns of landscape form in the upper Rhône basin, Central Swiss Alps, predominantly show lithologic controls despite multiple glaciations and variations in rock uplift rates

The development of topography depends mainly on the interplay between uplift and erosion. These processes are controlled by various factors including climate, glaciers, lithology, seismic activity and short-term variables, such as anthropogenic impact. Many studies in orogens all over the world have shown how these controlling variables may affect the landscape's topography. In particular, it has been hypothesized that lithology exerts a dominant control on erosion rates and landscape morphology. However, clear demonstrations of this influence are rare and difficult to disentangle from the overprint of other signals such as climate or tectonics. In this study we focus on the upper Rhône Basin situated in the Central Swiss Alps in order to explore the relation between topography, possible controlling variables and lithology in particular. The Rhône Basin has been affected by spatially variable uplift, high orographically driven rainfalls and multiple glaciations. Furthermore, lithology and erodibility vary substantially within the basin. Thanks to high-resolution geological, climatic and topographic data, the Rhône Basin is a suitable laboratory to explore these complexities. Elevation, relief, slope and hypsometric data as well as river profile information from digital elevation models are used to characterize the landscape's topography of around 50 tributary basins. Additionally, uplift over different timescales, glacial inheritance, precipitation patterns and erodibility of the underlying bedrock are quantified for each basin. Results show that the chosen topographic and controlling variables vary remarkably between different tributary basins. We investigate the link between observed topographic differences and the possible controlling variables through statistical analyses. Variations of elevation, slope and relief seem to be linked to differences in long-term uplift rate, whereas elevation distributions (hypsometry) and river profile shapes may be related to glacial imprint. This confirms that the landscape of the Rhône Basin has been highly preconditioned by (past) uplift and glaciation. Linear discriminant analyses (LDAs), however, suggest a stronger link between observed topographic variations and differences in erodibility. We therefore conclude that despite evident glacial and tectonic conditioning, a lithologic control is still preserved and measurable in the landscape of the Rhône tributary basins.

Geomorphology of the Imhotep region on comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Context. Since August 2014, the OSIRIS Narrow Angle Camera (NAC) onboard the Rosetta spacecraft has acquired high spatial resolution images of the nucleus of comet 67P/Churyumov-Gerasimenko, down to the decimeter scale. This paper focuses on the Imhotep region, located on the largest lobe of the nucleus, near the equator. Aims. We map, inventory, and describe the geomorphology of the Imhotep region. We propose and discuss some processes to explain the formation and ongoing evolution of this region. Methods. We used OSIRIS NAC images, gravitational heights and slopes, and digital terrain models to map and measure the morphologies of Imhotep. Results. The Imhotep region presents a wide variety of terrains and morphologies: smooth and rocky terrains, bright areas, linear features, roundish features, and boulders. Gravity processes such as mass wasting and collapse play a significant role in the geomorphological evolution of this region. Cometary processes initiate erosion and are responsible for the formation of degassing conduits that are revealed by elevated roundish features on the surface. We also propose a scenario for the formation and evolution of the Imhotep region; this implies the presence of large primordial voids inside the nucleus, resulting from its formation process.

Single fluoroscopic image based 2D/3D reconstruction of a scaled, patient-specific vertebral model

Purpose Accurate three-dimensional (3D) models of lumbar vertebrae can enable image-based 3D kinematic analysis. The common approach to derive 3D models is by direct segmentation of CT or MRI datasets. However, these have the disadvantages that they are expensive, timeconsuming and/or induce high-radiation doses to the patient. In this study, we present a technique to automatically reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image. Methods Our technique is based on a hybrid 2D/3D deformable registration strategy combining a landmark-to-ray registration with a statistical shape model-based 2D/3D reconstruction scheme. Fig. 1 shows different stages of the reconstruction process. Four cadaveric lumbar spine segments (total twelve lumbar vertebrae) were used to validate the technique. To evaluate the reconstruction accuracy, the surface models reconstructed from the lateral fluoroscopic images were compared to the associated ground truth data derived from a 3D CT-scan reconstruction technique. For each case, a surface-based matching was first used to recover the scale and the rigid transformation between the reconstructed surface model Results Our technique could successfully reconstruct 3D surface models of all twelve vertebrae. After recovering the scale and the rigid transformation between the reconstructed surface models and the ground truth models, the average error of the 2D/3D surface model reconstruction over the twelve lumbar vertebrae was found to be 1.0 mm. The errors of reconstructing surface models of all twelve vertebrae are shown in Fig. 2. It was found that the mean errors of the reconstructed surface models in comparison to their associated ground truths after iterative scaled rigid registrations ranged from 0.7 mm to 1.3 mm and the rootmean squared (RMS) errors ranged from 1.0 mm to 1.7 mm. The average mean reconstruction error was found to be 1.0 mm. Conclusion An accurate, scaled 3D reconstruction of the lumbar vertebra can be obtained from a single lateral fluoroscopic image using a statistical shape model based 2D/3D reconstruction technique. Future work will focus on applying the reconstructed model for 3D kinematic analysis of lumbar vertebrae, an extension of our previously-reported imagebased kinematic analysis. The developed method also has potential applications in surgical planning and navigation.

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.

Accurate and robust reconstruction of proximal femur from sparse intraoperative data and dense point distribution model for surgical navigation

Constructing a 3D surface model from sparse-point data is a nontrivial task. Here, we report an accurate and robust approach for reconstructing a surface model of the proximal femur from sparse-point data and a dense-point distribution model (DPDM). The problem is formulated as a three-stage optimal estimation process. The first stage, affine registration, is to iteratively estimate a scale and a rigid transformation between the mean surface model of the DPDM and the sparse input points. The estimation results of the first stage are used to establish point correspondences for the second stage, statistical instantiation, which stably instantiates a surface model from the DPDM using a statistical approach. This surface model is then fed to the third stage, kernel-based deformation, which further refines the surface model. Handling outliers is achieved by consistently employing the least trimmed squares (LTS) approach with a roughly estimated outlier rate in all three stages. If an optimal value of the outlier rate is preferred, we propose a hypothesis testing procedure to automatically estimate it. We present here our validations using four experiments, which include 1 leave-one-out experiment, 2 experiment on evaluating the present approach for handling pathology, 3 experiment on evaluating the present approach for handling outliers, and 4 experiment on reconstructing surface models of seven dry cadaver femurs using clinically relevant data without noise and with noise added. Our validation results demonstrate the robust performance of the present approach in handling outliers, pathology, and noise. An average 95-percentile error of 1.7-2.3 mm was found when the present approach was used to reconstruct surface models of the cadaver femurs from sparse-point data with noise added.

Valgus Hip With High Antetorsion Causes Pain Through Posterior Extraarticular FAI

BACKGROUND: Valgus hips with increased antetorsion present with lack of external rotation and posterior hip pain that is aggravated with hip extension and external rotation. This may be the result of posterior femoroacetabular impingement (FAI). QUESTIONS/PURPOSES: We asked whether (1) the range of motion (ROM); (2) the location of anterior and posterior bony collision zones; and (3) the prevalence of extraarticular impingement differ between valgus hips with increased antetorsion compared with normal hips and hips with idiopathic FAI. METHODS: Surface models based on CT scan reconstructions of 13 valgus hips with increased antetorsion, 22 hips with FAI, and 27 normal hips were included. Validated three-dimensional collision detection software was used to quantify the simulated hip ROM and the location of impingement on the acetabular and the femoral sides. RESULTS: Hips with coxa valga and antetorsion showed decreased extension, external rotation, and adduction, whereas internal rotation in 90° of flexion was increased. Impingement zones were more anteroinferior on the femur and posteroinferior on the acetabular (pelvic) side; and the zones were more frequently extraarticular, posterior, or to a lesser degree anterior against the inferior iliac spine. We found a higher prevalence of extraarticular impingement for valgus hips with increased antetorsion. CONCLUSIONS: Valgus hips with increased antetorsion predispose to posterior extraarticular FAI and to a lesser degree anteroinferior spine impingement. LEVEL OF EVIDENCE: Level II, prognostic study. See Guidelines for Authors for a complete description of levels of evidence.