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.

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.

The effect of pulsed jet lavage in vertebroplasty on injection forces of polymethylmethacrylate bone cement, material distribution, and potential fat embolism: a cadaver study

STUDY DESIGN: In vitro testing of vertebroplasty techniques including pulsed jet-lavage for fat and marrow removal in human cadaveric lumbar and thoracic vertebrae. OBJECTIVE: To develop jet-lavage techniques for vertebroplasty and investigate their effect on cement distribution, injection forces, and fat embolism. SUMMARY OF BACKGROUND DATA: The main complications of cement vertebroplasty are cement leakage and pulmonary fat embolism, which can have fatal consequences and are difficult to prevent reliably by current vertebroplasty techniques. METHODS: Twenty-four vertebrae (Th8-L04) from 5 osteoporotic cadaver spines were grouped in triplets depending on bone mineral density (BMD). Before polymethylmethacrylate (PMMA) vertebroplasty, a pulsatile jet-lavage for removal of intertrabecular fat and bone marrow was performed in 2 groups with 8 specimens each, performing radial and axial irrigation from the biopsy needles. One hundred mL of Ringer solution were injected through 1 pedicle and regained by low vacuum via the contralateral pedicle. Eight control vertebrae were not irrigated. All specimens underwent standardized PMMA cement augmentation injecting 20% of the vertebral volume. Injection forces, cement distribution, and extravasations were quantified. RESULTS: All irrigation solution could be retrieved with the vacuum applied. A Kruskal-Wallis test revealed significantly higher injection forces of the control group as compared with the irrigated groups (P = 0.021). Dilatation of the syringe at forces above 300 N occurred in 75% of the untreated compared with 12.5% of the lavaged specimens. CT distribution analysis showed more homogenous cement distribution of the cement and significantly less extravasation in the irrigated specimens. CONCLUSION: The developed lavage technique for vertebroplasty showed to be feasible and reproducible. The reduction of injection forces would allow the use of more viscous PMMA cement lowering the risk for cement embolization and results in a safer procedure. The wash-out of bone marrow and the possible reduction of pulmonary fat embolism have to be verified with in vivo models.

Cement distribution, volume, and compliance in vertebroplasty: some answers from an anatomy-based nonlinear finite element study

STUDY DESIGN: The biomechanics of vertebral bodies augmented with real distributions of cement were investigated using nonlinear finite element (FE) analysis. OBJECTIVES: To compare stiffness, strength, and stress transfer of augmented versus nonaugmented osteoporotic vertebral bodies under compressive loading. Specifically, to examine how cement distribution, volume, and compliance affect these biomechanical variables. SUMMARY OF BACKGROUND DATA: Previous FE studies suggested that vertebroplasty might alter vertebral stress transfer, leading to adjacent vertebral failure. However, no FE study so far accounted for real cement distributions and bone damage accumulation. METHODS: Twelve vertebral bodies scanned with high-resolution pQCT and tested in compression were augmented with various volumes of cements and scanned again. Nonaugmented and augmented pQCT datasets were converted to FE models, with bone properties modeled with an elastic, plastic and damage constitutive law that was previously calibrated for the nonaugmented models. The cement-bone composite was modeled with a rule of mixture. The nonaugmented and augmented FE models were subjected to compression and their stiffness, strength, and stress map calculated for different cement compliances. RESULTS: Cement distribution dominated the stiffening and strengthening effects of augmentation. Models with cement connecting either the superior or inferior endplate (S/I fillings) were only up to 2 times stiffer than the nonaugmented models with minimal strengthening, whereas those with cement connecting both endplates (S + I fillings) were 1 to 8 times stiffer and 1 to 12 times stronger. Stress increases above and below the cement, which was higher for the S + I cases and was significantly reduced by increasing cement compliance. CONCLUSION: The developed FE approach, which accounts for real cement distributions and bone damage accumulation, provides a refined insight into the mechanics of augmented vertebral bodies. In particular, augmentation with compliant cement bridging both endplates would reduce stress transfer while providing sufficient strengthening.

Theoretical models of planetary system formation: mass vs. semi-major axis

Context. Planet formation models have been developed during the past years to try to reproduce what has been observed of both the solar system and the extrasolar planets. Some of these models have partially succeeded, but they focus on massive planets and, for the sake of simplicity, exclude planets belonging to planetary systems. However, more and more planets are now found in planetary systems. This tendency, which is a result of radial velocity, transit, and direct imaging surveys, seems to be even more pronounced for low-mass planets. These new observations require improving planet formation models, including new physics, and considering the formation of systems. Aims: In a recent series of papers, we have presented some improvements in the physics of our models, focussing in particular on the internal structure of forming planets, and on the computation of the excitation state of planetesimals and their resulting accretion rate. In this paper, we focus on the concurrent effect of the formation of more than one planet in the same protoplanetary disc and show the effect, in terms of architecture and composition of this multiplicity. Methods: We used an N-body calculation including collision detection to compute the orbital evolution of a planetary system. Moreover, we describe the effect of competition for accretion of gas and solids, as well as the effect of gravitational interactions between planets. Results: We show that the masses and semi-major axes of planets are modified by both the effect of competition and gravitational interactions. We also present the effect of the assumed number of forming planets in the same system (a free parameter of the model), as well as the effect of the inclination and eccentricity damping. We find that the fraction of ejected planets increases from nearly 0 to 8% as we change the number of embryos we seed the system with from 2 to 20 planetary embryos. Moreover, our calculations show that, when considering planets more massive than ~5 M⊕, simulations with 10 or 20 planetary embryos statistically give the same results in terms of mass function and period distribution.

Clinical versus pre-clinical FE models for vertebral body strength predictions

The finite element analysis is an accepted method to predict vertebral body compressive strength. This study compares measurements obtained from in vitro tests with the ones from two different simulation models: clinical quantitative computer tomography (QCT) based homogenized finite element (hFE) models and pre-clinical high-resolution peripheral QCT-based (HR-pQCT) hFE models. About 37 vertebral body sections were prepared by removing end-plates and posterior elements, scanned with QCT (390/450μm voxel size) as well as HR-pQCT (82μm voxel size), and tested in compression up to failure. Non-linear viscous damage hFE models were created from QCT/HT-pQCT images and compared to experimental results based on stiffness and ultimate load. As expected, the predictability of QCT/HR-pQCT-based hFE models for both apparent stiffness (r2=0.685/0.801r2=0.685/0.801) and strength (r2=0.774/0.924r2=0.774/0.924) increased if a better image resolution was used. An analysis of the damage distribution showed similar damage locations for all cases. In conclusion, HR-pQCT-based hFE models increased the predictability considerably and do not need any tuning of input parameters. In contrast, QCT-based hFE models usually need some tuning but are clinically the only possible choice at the moment.

A critical evaluation of secondary cancer risk models applied to Monte Carlo dose distributions of 2-dimensional, 3-dimensional conformal and hybrid intensity-modulated radiation therapy for breast cancer

The comparison of radiotherapy techniques regarding secondary cancer risk has yielded contradictory results possibly stemming from the many different approaches used to estimate risk. The purpose of this study was to make a comprehensive evaluation of different available risk models applied to detailed whole-body dose distributions computed by Monte Carlo for various breast radiotherapy techniques including conventional open tangents, 3D conformal wedged tangents and hybrid intensity modulated radiation therapy (IMRT). First, organ-specific linear risk models developed by the International Commission on Radiological Protection (ICRP) and the Biological Effects of Ionizing Radiation (BEIR) VII committee were applied to mean doses for remote organs only and all solid organs. Then, different general non-linear risk models were applied to the whole body dose distribution. Finally, organ-specific non-linear risk models for the lung and breast were used to assess the secondary cancer risk for these two specific organs. A total of 32 different calculated absolute risks resulted in a broad range of values (between 0.1% and 48.5%) underlying the large uncertainties in absolute risk calculation. The ratio of risk between two techniques has often been proposed as a more robust assessment of risk than the absolute risk. We found that the ratio of risk between two techniques could also vary substantially considering the different approaches to risk estimation. Sometimes the ratio of risk between two techniques would range between values smaller and larger than one, which then translates into inconsistent results on the potential higher risk of one technique compared to another. We found however that the hybrid IMRT technique resulted in a systematic reduction of risk compared to the other techniques investigated even though the magnitude of this reduction varied substantially with the different approaches investigated. Based on the epidemiological data available, a reasonable approach to risk estimation would be to use organ-specific non-linear risk models applied to the dose distributions of organs within or near the treatment fields (lungs and contralateral breast in the case of breast radiotherapy) as the majority of radiation-induced secondary cancers are found in the beam-bordering regions.

Impact of synchronous metastasis distribution on cancer specific survival in renal cell carcinoma after radical nephrectomy with tumor thrombectomy

PURPOSE Metastatic renal cell carcinoma can be clinically diverse in terms of the pattern of metastatic disease and response to treatment. We studied the impact of metastasis and location on cancer specific survival. MATERIALS AND METHODS The records of 2,017 patients with renal cell cancer and tumor thrombus who underwent radical nephrectomy and tumor thrombectomy from 1971 to 2012 at 22 centers in the United States and Europe were analyzed. Number and location of synchronous metastases were compared with respect to patient cancer specific survival. Multivariable Cox regression models were used to quantify the impact of covariates. RESULTS Lymph node metastasis (155) or distant metastasis (725) was present in 880 (44%) patients. Of the patients with distant disease 385 (53%) had an isolated metastasis. The 5-year cancer specific survival was 51.3% (95% CI 48.6-53.9) for the entire group. On univariable analysis patients with isolated lymph node metastasis had a significantly worse cancer specific survival than those with a solitary distant metastasis. The location of distant metastasis did not have any significant effect on cancer specific survival. On multivariable analysis the presence of lymph node metastasis, isolated distant metastasis and multiple distant metastases were independently associated with cancer specific survival. Moreover higher tumor thrombus level, papillary histology and the use of postoperative systemic therapy were independently associated with worse cancer specific survival. CONCLUSIONS In our multi-institutional series of patients with renal cell cancer who underwent radical nephrectomy and tumor thrombectomy, almost half of the patients had synchronous lymph node or distant organ metastasis. Survival was superior in patients with solitary distant metastasis compared to isolated lymph node disease.

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.

Associations of obesity and body fat distribution with incident atrial fibrillation in the biracial health aging and body composition cohort of older adults

UNLABELLED Obesity is a well-recognized risk factor for atrial fibrillation (AF), yet adiposity measures other than body mass index (BMI) have had limited assessment in relation to AF risk. We examined the associations of adiposity measures with AF in a biracial cohort of older adults. Given established racial differences in obesity and AF, we assessed for differences by black and white race in relating adiposity and AF. METHODS We analyzed data from 2,717 participants of the Health, Aging, and Body Composition Study. Adiposity measures were BMI, abdominal circumference, subcutaneous and visceral fat area, and total and percent fat mass. We determined the associations between the adiposity measures and 10-year incidence of AF using Cox proportional hazards models and assessed for their racial differences in these estimates. RESULTS In multivariable-adjusted models, 1-SD increases in BMI, abdominal circumference, and total fat mass were associated with a 13% to 16% increased AF risk (hazard ratio [HR] 1.14, 95% CI 1.02-1.28; HR 1.16, 95% CI 1.04-1.28; and HR 1.13, 95% CI 1.002-1.27). Subcutaneous and visceral fat areas were not significantly associated with incident AF. We did not identify racial differences in the associations between the adiposity measures and AF. CONCLUSION Body mass index, abdominal circumference, and total fat mass are associated with risk of AF for 10years among white and black older adults. Obesity is one of a limited number of modifiable risk factors for AF; future studies are essential to evaluate how obesity reduction can modify the incidence of AF.

Variability of daily winter wind speed distribution over Northern Europe during the past millennium in regional and global climate simulations

We analyse the variability of the probability distribution of daily wind speed in wintertime over Northern and Central Europe in a series of global and regional climate simulations covering the last centuries, and in reanalysis products covering approximately the last 60 years. The focus of the study lies on identifying the link of the variations in the wind speed distribution to the regional near-surface temperature, to the meridional temperature gradient and to the North Atlantic Oscillation. Our main result is that the link between the daily wind distribution and the regional climate drivers is strongly model dependent. The global models tend to behave similarly, although they show some discrepancies. The two regional models also tend to behave similarly to each other, but surprisingly the results derived from each regional model strongly deviates from the results derived from its driving global model. In addition, considering multi-centennial timescales, we find in two global simulations a long-term tendency for the probability distribution of daily wind speed to widen through the last centuries. The cause for this widening is likely the effect of the deforestation prescribed in these simulations. We conclude that no clear systematic relationship between the mean temperature, the temperature gradient and/or the North Atlantic Oscillation, with the daily wind speed statistics can be inferred from these simulations. The understand- ing of past and future changes in the distribution of wind speeds, and thus of wind speed extremes, will require a detailed analysis of the representation of the interaction between large-scale and small-scale dynamics.

Theoretical models of planetary system formation II. Post-formation evolution

Aims. We extend the results of planetary formation synthesis by computing the long-term evolution of synthetic systems from the clearing of the gas disk into the dynamical evolution phase. Methods. We use the symplectic integrator SyMBA to numerically integrate the orbits of planets for 100 Myr, using populations from previous studies as initial conditions. Results. We show that within the populations studied, mass and semimajor axis distributions experience only minor changes from post-formation evolution. We also show that, depending upon their initial distribution, planetary eccentricities can statistically increase or decrease as a result of gravitational interactions. We find that planetary masses and orbital spacings provided by planet formation models do not result in eccentricity distributions comparable to observed exoplanet eccentricities, requiring other phenomena, such as stellar fly-bys, to account for observed eccentricities.