Hardware accelerated ray cast of volume data and volume gradient for an optimized splines-based multi-resolution 2D-3D registration

This paper describes a method for DRR generation as well as for volume gradients projection using hardware accelerated 2D texture mapping and accumulation buffering and demonstrates its application in 2D-3D registration of X-ray fluoroscopy to CT images. The robustness of the present registration scheme are guaranteed by taking advantage of a coarse-to-fine processing of the volume/image pyramids based on cubic B-splines. A human cadaveric spine specimen together with its ground truth was used to compare the present scheme with a purely software-based scheme in three aspects: accuracy, speed, and capture ranges. Our experiments revealed an equivalent accuracy and capture ranges but with much shorter registration time with the present scheme. More specifically, the results showed 0.8 mm average target registration error, 55 second average execution time per registration, and 10 mm and 10° capture ranges for the present scheme when tested on a 3.0 GHz Pentium 4 computer.

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and point distribution model

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and a point distribution model is discussed. We present a 2D/3D reconstruction scheme combining statistical extrapolation and regularized shape deformation with an iterative image-to-model correspondence establishing algorithm, and show its application to reconstruct the surface of proximal femur. The image-to-model correspondence is established using a non-rigid 2D point matching process, which iteratively uses a symmetric injective nearest-neighbor mapping operator and 2D thin-plate splines based deformation to find a fraction of best matched 2D point pairs between features detected from the fluoroscopic images and those extracted from the 3D model. The obtained 2D point pairs are then used to set up a set of 3D point pairs such that we turn a 2D/3D reconstruction problem to a 3D/3D one. We designed and conducted experiments on 11 cadaveric femurs to validate the present reconstruction scheme. An average mean reconstruction error of 1.2 mm was found when two fluoroscopic images were used for each bone. It decreased to 1.0 mm when three fluoroscopic images were used.

Impact of unlinked deaths and coding changes on mortality trends in the Swiss National Cohort

BACKGROUND Results of epidemiological studies linking census with mortality records may be affected by unlinked deaths and changes in cause of death classification. We examined these issues in the Swiss National Cohort (SNC). METHODS The SNC is a longitudinal study of the entire Swiss population, based on the 1990 (6.8 million persons) and 2000 (7.3 million persons) censuses. Among 1,053,393 deaths recorded 1991-2007 5.4% could not be linked using stringent probabilistic linkage. We included the unlinked deaths using pragmatic linkages and compared mortality rates for selected causes with official mortality rates. We also examined the impact of the 1995 change in cause of death coding from version 8 (with some additional rules) to version 10 of the International Classification of Diseases (ICD), using Poisson regression models with restricted cubic splines. Finally, we compared results from Cox models including and excluding unlinked deaths of the association of education, marital status, and nationality with selected causes of death. RESULTS SNC mortality rates underestimated all cause mortality by 9.6% (range 2.4%-17.9%) in the 85+ population. Underestimation was less pronounced in years nearer the censuses and in the 75-84 age group. After including 99.7% of unlinked deaths, annual all cause SNC mortality rates were reflecting official rates (relative difference between -1.4% and +1.8%). In the 85+ population the rates for prostate and breast cancer dropped, by 16% and 21% respectively, between 1994 and 1995 coincident with the change in cause of death coding policy. For suicide in males almost no change was observed. Hazard ratios were only negligibly affected by including the unlinked deaths. A sudden decrease in breast (21% less, 95% confidence interval: 12%-28%) and prostate (16% less, 95% confidence interval: 7%-23%) cancer mortality rates in the 85+ population coincided with the 1995 change in cause of death coding policy. CONCLUSIONS Unlinked deaths bias analyses of absolute mortality rates downwards but have little effect on relative mortality. To describe time trends of cause-specific mortality in the SNC, accounting for the unlinked deaths and for the possible effect of change in death certificate coding was necessary.

Blood lead level association with lower body weight in NHANES 1999-2006

BACKGROUND Lead exposure is associated with low birth-weight. The objective of this study is to determine whether lead exposure is associated with lower body weight in children, adolescents and adults. METHODS We analyzed data from NHANES 1999-2006 for participants aged ≥3 using multiple logistic and multivariate linear regression. Using age- and sex-standardized BMI Z-scores, overweight and obese children (ages 3-19) were classified by BMI ≥85 th and ≥95 th percentiles, respectively. The adult population (age ≥20) was classified as overweight and obese with BMI measures of 25-29.9 and ≥30, respectively. Blood lead level (BLL) was categorized by weighted quartiles. RESULTS Multivariate linear regressions revealed a lower BMI Z-score in children and adolescents when the highest lead quartile was compared to the lowest lead quartile (β (SE)=-0.33 (0.07), p<0.001), and a decreased BMI in adults (β (SE)=-2.58 (0.25), p<0.001). Multiple logistic analyses in children and adolescents found a negative association between BLL and the percentage of obese and overweight with BLL in the highest quartile compared to the lowest quartile (OR=0.42, 95% CI: 0.30-0.59; and OR=0.67, 95% CI: 0.52-0.88, respectively). Adults in the highest lead quartile were less likely to be obese (OR=0.42, 95% CI: 0.35-0.50) compared to those in the lowest lead quartile. Further analyses with blood lead as restricted cubic splines, confirmed the dose-relationship between blood lead and body weight outcomes. CONCLUSIONS BLLs are associated with lower body mass index and obesity in children, adolescents and adults.

Socioeconomic, Temporal and Regional Variation in Body Mass Index among 188,537 Swiss Male Conscripts Born between 1986 and 1992

BACKGROUND Rising levels of overweight and obesity are important public-health concerns worldwide. The purpose of this study is to elucidate their prevalence and trends in Switzerland by analyzing variations in Body Mass Index (BMI) of Swiss conscripts. METHODS The conscription records were provided by the Swiss Army. This study focussed on conscripts 18.5-20.5 years of age from the seven one-year birth cohorts spanning the period 1986-1992. BMI across professional status, area-based socioeconomic position (abSEP), urbanicity and regions was analyzed. Two piecewise quantile regression models with linear splines for three birth-cohort groups were used to examine the association of median BMI with explanatory variables and to determine the extent to which BMI has varied over time. RESULTS The study population consisted of 188,537 individuals. Median BMI was 22.51 kg/m2 (22.45-22.57 95% confidence interval (CI)). BMI was lower among conscripts of high professional status (-0.46 kg/m2; 95% CI: -0.50, -0.42, compared with low), living in areas of high abSEP (-0.11 kg/m2; 95% CI: -0.16, -0.07 compared to medium) and from urban communities (-0.07 kg/m2; 95% CI: -0.11, -0.03, compared with peri-urban). Comparing with Midland, median BMI was highest in the North-West (0.25 kg/m2; 95% CI: 0.19-0.30) and Central regions (0.11 kg/m2; 95% CI: 0.05-0.16) and lowest in the East (-0.19 kg/m2; 95% CI: -0.24, -0.14) and Lake Geneva regions (-0.15 kg/m2; 95% CI: -0.20, -0.09). Trajectories of regional BMI growth varied across birth cohorts, with median BMI remaining high in the Central and North-West regions, whereas stabilization and in some cases a decline were observed elsewhere. CONCLUSIONS BMI of Swiss conscripts is associated with individual and abSEP and urbanicity. Results show regional variation in the levels and temporal trajectories of BMI growth and signal their possible slowdown among recent birth cohorts.

Estimating HIV Incidence, Time to Diagnosis, and the Undiagnosed HIV Epidemic Using Routine Surveillance Data.

BACKGROUND Estimates of the size of the undiagnosed HIV-infected population are important to understand the HIV epidemic and to plan interventions, including "test-and-treat" strategies. METHODS We developed a multi-state back-calculation model to estimate HIV incidence, time between infection and diagnosis, and the undiagnosed population by CD4 count strata, using surveillance data on new HIV and AIDS diagnoses. The HIV incidence curve was modelled using cubic splines. The model was tested on simulated data and applied to surveillance data on men who have sex with men in The Netherlands. RESULTS The number of HIV infections could be estimated accurately using simulated data, with most values within the 95% confidence intervals of model predictions. When applying the model to Dutch surveillance data, 15,400 (95% confidence interval [CI] = 15,000, 16,000) men who have sex with men were estimated to have been infected between 1980 and 2011. HIV incidence showed a bimodal distribution, with peaks around 1985 and 2005 and a decline in recent years. Mean time to diagnosis was 6.1 (95% CI = 5.8, 6.4) years between 1984 and 1995 and decreased to 2.6 (2.3, 3.0) years in 2011. By the end of 2011, 11,500 (11,000, 12,000) men who have sex with men in The Netherlands were estimated to be living with HIV, of whom 1,750 (1,450, 2,200) were still undiagnosed. Of the undiagnosed men who have sex with men, 29% (22, 37) were infected for less than 1 year, and 16% (13, 20) for more than 5 years. CONCLUSIONS This multi-state back-calculation model will be useful to estimate HIV incidence, time to diagnosis, and the undiagnosed HIV epidemic based on routine surveillance data.

NLCHECK: Stata module to check linearity assumption after model estimation

nlcheck is a simple diagnostic tool that can be used after fitting a model to quickly check the linearity assumption for a given predictor. nlcheck categorizes the predictor into bins, refits the model including dummy variables for the bins, and then performs a joint Wald test for the added parameters. Alternative, nlcheck uses linear splines for the adaptive model. Support for discrete variables is also provided. Optionally, nlcheck also displays a graph of the adjusted linear predictions from the original model and the adaptive model

Lake-effect snow as the dominant control of mesic-forest distribution in Michigan, USA

1Recent studies demonstrated the sensitivity of northern forest ecosystems to changes in the amount and duration of snow cover at annual to decadal time scales. However, the consequences of snowfall variability remain uncertain for ecological variables operating at longer time scales, especially the distributions of forest communities. 2The Great Lakes region of North America offers a unique setting to examine the long-term effects of variable snowfall on forest communities. Lake-effect snow produces a three-fold gradient in annual snowfall over tens of kilometres, and dramatic edaphic variations occur among landform types resulting from Quaternary glaciations. We tested the hypothesis that these factors interact to control the distributions of mesic (dominated by Acer saccharum, Tsuga canadensis and Fagus grandifolia) and xeric forests (dominated by Pinus and Quercus spp.) in northern Lower Michigan. 3We compiled pre-European-settlement vegetation data and overlaid these data with records of climate, water balance and soil, onto Landtype Association polygons in a geographical information system. We then used multivariate adaptive regression splines to model the abundance of mesic vegetation in relation to environmental controls. 4Snowfall is the most predictive among five variables retained by our model, and it affects model performance 29% more than soil texture, the second most important variable. The abundance of mesic trees is high on fine-textured soils regardless of snowfall, but it increases with snowfall on coarse-textured substrates. Lake-effect snowfall also determines the species composition within mesic forests. The weighted importance of A. saccharum is significantly greater than of T. canadensis or F. grandifolia within the lake-effect snowbelt, whereas T. canadensis is more plentiful outside the snowbelt. These patterns are probably driven by the influence of snowfall on soil moisture, nutrient availability and fire return intervals. 5Our results imply that a key factor dictating the spatio-temporal patterns of forest communities in the vast region around the Great Lakes is how the lake-effect snowfall regime responds to global change. Snowfall reductions will probably cause a major decrease in the abundance of ecologically and economically important species, such as A. saccharum.