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

A model for radiofrequency electromagnetic field predictions at outdoor and indoor locations in the context of epidemiological research

We present a geospatial model to predict the radiofrequency electromagnetic field from fixed site transmitters for use in epidemiological exposure assessment. The proposed model extends an existing model toward the prediction of indoor exposure, that is, at the homes of potential study participants. The model is based on accurate operation parameters of all stationary transmitters of mobile communication base stations, and radio broadcast and television transmitters for an extended urban and suburban region in the Basel area (Switzerland). The model was evaluated by calculating Spearman rank correlations and weighted Cohen's kappa (kappa) statistics between the model predictions and measurements obtained at street level, in the homes of volunteers, and in front of the windows of these homes. The correlation coefficients of the numerical predictions with street level measurements were 0.64, with indoor measurements 0.66, and with window measurements 0.67. The kappa coefficients were 0.48 (95%-confidence interval: 0.35-0.61) for street level measurements, 0.44 (95%-CI: 0.32-0.57) for indoor measurements, and 0.53 (95%-CI: 0.42-0.65) for window measurements. Although the modeling of shielding effects by walls and roofs requires considerable simplifications of a complex environment, we found a comparable accuracy of the model for indoor and outdoor points.

Histological evidence of delayed ischemic brain tissue damage in the rat double-hemorrhage model

The rat double-SAH model is one of the standard models to simulate delayed cerebral vasospasm (CVS) in humans. However, the proof of delayed ischemic brain damage is missing so far. Our objective was, therefore, to determine histological changes in correlation with the development of symptomatic and perfusion weighted imaging (PWI) proven CVS in this animal model. CVS was induced by injection of autologous blood in the cisterna magna of 22 Sprague-Dawley rats. Histological changes were analyzed on day 3 and day 5. Cerebral blood flow (CBF) was assessed by PWI at 3 tesla magnetic resonance (MR) tomography. Neuronal cell count did not differ between sham operated and SAH rats in the hippocampus and the cerebral cortex on day 3. In contrast, on day 5 after SAH the neuronal cell count was significantly reduced in the hippocampus (p<0.001) and the inner cortical layer (p=0.03). The present investigation provides quantitative data on brain tissue damage in association with delayed CVS for the first time in a rat SAH model. Accordingly, our data suggest that the rat double-SAH model may be suitable to mimic delayed ischemic brain damage due to CVS and to investigate the neuroprotective effects of drugs.

Genome-wide investigation reveals high evolutionary rates in annual model plants

Background: ;Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast;genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.;Results: ;According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall;evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.;Conclusions: ;The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the;mechanism determining the molecular evolutionary rate at the genomic level.

A prediction model for assessing residential radon concentration in Switzerland

Indoor radon is regularly measured in Switzerland. However, a nationwide model to predict residential radon levels has not been developed. The aim of this study was to develop a prediction model to assess indoor radon concentrations in Switzerland. The model was based on 44,631 measurements from the nationwide Swiss radon database collected between 1994 and 2004. Of these, 80% randomly selected measurements were used for model development and the remaining 20% for an independent model validation. A multivariable log-linear regression model was fitted and relevant predictors selected according to evidence from the literature, the adjusted R², the Akaike's information criterion (AIC), and the Bayesian information criterion (BIC). The prediction model was evaluated by calculating Spearman rank correlation between measured and predicted values. Additionally, the predicted values were categorised into three categories (50th, 50th-90th and 90th percentile) and compared with measured categories using a weighted Kappa statistic. The most relevant predictors for indoor radon levels were tectonic units and year of construction of the building, followed by soil texture, degree of urbanisation, floor of the building where the measurement was taken and housing type (P-values <0.001 for all). Mean predicted radon values (geometric mean) were 66 Bq/m³ (interquartile range 40-111 Bq/m³) in the lowest exposure category, 126 Bq/m³ (69-215 Bq/m³) in the medium category, and 219 Bq/m³ (108-427 Bq/m³) in the highest category. Spearman correlation between predictions and measurements was 0.45 (95%-CI: 0.44; 0.46) for the development dataset and 0.44 (95%-CI: 0.42; 0.46) for the validation dataset. Kappa coefficients were 0.31 for the development and 0.30 for the validation dataset, respectively. The model explained 20% overall variability (adjusted R²). In conclusion, this residential radon prediction model, based on a large number of measurements, was demonstrated to be robust through validation with an independent dataset. The model is appropriate for predicting radon level exposure of the Swiss population in epidemiological research. Nevertheless, some exposure misclassification and regression to the mean is unavoidable and should be taken into account in future applications of the model.

Cardiac transplantation with hearts from donors after circulatory declaration of death: haemodynamic and biochemical parameters at procurement predict recovery following cardioplegic storage in a rat model

OBJECTIVES: Donation after circulatory declaration of death (DCDD) could significantly improve the number of cardiac grafts for transplantation. Graft evaluation is particularly important in the setting of DCDD given that conditions of cardio-circulatory arrest and warm ischaemia differ, leading to variable tissue injury. The aim of this study was to identify, at the time of heart procurement, means to predict contractile recovery following cardioplegic storage and reperfusion using an isolated rat heart model. Identification of reliable approaches to evaluate cardiac grafts is key in the development of protocols for heart transplantation with DCDD. METHODS: Hearts isolated from anaesthetized male Wistar rats (n = 34) were exposed to various perfusion protocols. To simulate DCDD conditions, rats were exsanguinated and maintained at 37°C for 15-25 min (warm ischaemia). Isolated hearts were perfused with modified Krebs-Henseleit buffer for 10 min (unloaded), arrested with cardioplegia, stored for 3 h at 4°C and then reperfused for 120 min (unloaded for 60 min, then loaded for 60 min). Left ventricular (LV) function was assessed using an intraventricular micro-tip pressure catheter. Statistical significance was determined using the non-parametric Spearman rho correlation analysis. RESULTS: After 120 min of reperfusion, recovery of LV work measured as developed pressure (DP)-heart rate (HR) product ranged from 0 to 15 ± 6.1 mmHg beats min(-1) 10(-3) following warm ischaemia of 15-25 min. Several haemodynamic parameters measured during early, unloaded perfusion at the time of heart procurement, including HR and the peak systolic pressure-HR product, correlated significantly with contractile recovery after cardioplegic storage and 120 min of reperfusion (P < 0.001). Coronary flow, oxygen consumption and lactate dehydrogenase release also correlated significantly with contractile recovery following cardioplegic storage and 120 min of reperfusion (P < 0.05). CONCLUSIONS: Haemodynamic and biochemical parameters measured at the time of organ procurement could serve as predictive indicators of contractile recovery. We believe that evaluation of graft suitability is feasible prior to transplantation with DCDD, and may, consequently, increase donor heart availability.

Noninvasive Doppler-derived myocardial performance index in rats with myocardial infarction: validation and correlation by conductance catheter

The rodent model of myocardial infarction (MI) is extensively used in heart failure studies. However, long-term follow-up of echocardiographic left ventricular (LV) function parameters such as the myocardial performance index (MPI) and its ratio with the fractional shortening (LVFS/MPI) has not been validated in conjunction with invasive indexes, such as those derived from the conductance catheter (CC). Sprague-Dawley rats with left anterior descending coronary artery ligation (MI group, n = 9) were compared with a sham-operated control group (n = 10) without MI. Transthoracic echocardiography (TTE) was performed every 2 wk over an 8-wk period, after which classic TTE parameters, especially MPI and LVFS/MPI, were compared with invasive indexes obtained by using a CC. Serial TTE data showed significant alterations in the majority of the noninvasive functional and structural parameters (classic and novel) studied in the presence of MI. Both MPI and LVFS/MPI significantly (P < 0.05 for all reported values) correlated with body weight (r = -0.58 and 0.76 for MPI and LVFS/MPI, respectively), preload recruitable stroke work (r = -0.61 and 0.63), LV end-diastolic pressure (LVEDP) (r = 0.82 and -0.80), end-diastolic volume (r = 0.61 and -0.58), and end-systolic volume (r = 0.46 and -0.48). Forward stepwise linear regression analysis revealed that, of all variables tested, LVEDP was the only independent determinant of MPI (r = 0.84) and LVFS/MPI (r = 0.83). We conclude that MPI and LVFS/MPI correlate strongly and better than the classic noninvasive TTE parameters with established, invasively assessed indexes of contractility, preload, and volumetry. These findings support the use of these two new noninvasive indexes for long-term analysis of the post-MI LV remodeling.

Experimental model to evaluate in vivo and in vitro cartilage MR imaging by means of histological analyses

OBJECTIVES: Implementation of an experimental model to compare cartilage MR imaging by means of histological analyses. MATERIAL AND METHODS: MRI was obtained from 4 patients expecting total knee replacement at 1.5 and/or 3T prior surgery. The timeframe between pre-op MRI and knee replacement was within two days. Resected cartilage-bone samples were tagged with Ethi((R))-pins to reproduce the histological cutting course. Pre-operative scanning at 1.5T included following parameters for fast low angle shot (FLASH: TR/TE/FA=33ms/6ms/30 degrees , BW=110kHz, 120mmx120mm FOV, 256x256 matrix, 0.65mm slice-thickness) and double echo steady state (DESS: TR/TE/FA=23.7ms/6.9ms/40 degrees , BW=130kHz, 120x120mm FOV, 256x256 matrix, 0.65mm slice-thickness). At 3T, scan parameters were: FLASH (TR/TE/FA=12.2ms/5.1ms/10 degrees , BW=130kHz, 170x170mm FOV, 320x320, 0.5mm slice-thickness) and DESS (TR/TE/FA=15.6ms/4.5ms/25 degrees , BW=200kHz, 135mmx150mm FOV, 288x320matrix, 0.5mm slice-thickness). Imaging of the specimens was done the same day at 1.5T. MRI (Noyes) and histological (Mankin) score scales were correlated using the paired t-test. Sensitivity and specificity for the detection of different grades of cartilage degeneration were assessed. Inter-reader and intra-reader reliability was determined using Kappa analysis. RESULTS: Low correlation (sensitivity, specificity) was found for both sequences in normal to mild Mankin grades. Only moderate to severe changes were diagnosed with higher significance and specificity. The use of higher field-strengths was advantageous for both protocols with sensitivity values ranging from 13.6% to 93.3% (FLASH) and 20.5% to 96.2% (DESS). Kappa values ranged from 0.488 to 0.944. CONCLUSIONS: Correlating MR images with continuous histological slices was feasible by using three-dimensional imaging, multi-planar-reformat and marker pins. The capability of diagnosing early cartilage changes with high accuracy could not be proven for both FLASH and DESS.

Correlation between skeletal trauma and energy in falls from great height detected by post-mortem multislice computed tomography (MSCT)

Fatal falls from great height are a frequently encountered setting in forensic pathology. They present--by virtue of a calculable energy transmission to the body--an ideal model for the assessment of the effects of blunt trauma to a human body. As multislice computed tomography (MSCT) has proven not only to be invaluable in clinical examinations, but also to be a viable tool in post-mortem imaging, especially in the field of osseous injuries, we performed a MSCT scan on 20 victims of falls from great height. We hereby detected fractures and their distributions were compared with the impact energy. Our study suggests a marked increase of extensive damage to different body regions at about 20 kJ and more. The thorax was most often affected, regardless of the amount of impacting energy and the primary impact site. Cranial fracture frequency displayed a biphasic distribution with regard to the impacting energy; they were more frequent in energies of less than 10, and more than 20 kJ, but rarer in the intermediate energy group, namely that of 10-20 kJ.

Near-infrared digital photography to estimate snow correlation length for microwave emission modeling

The study is based on experimental work conducted in alpine snow. We made microwave radiometric and near-infrared reflectance measurements of snow slabs under different experimental conditions. We used an empirical relation to link near-infrared reflectance of snow to the specific surface area (SSA), and converted the SSA into the correlation length. From the measurements of snow radiances at 21 and 35 GHz , we derived the microwave scattering coefficient by inverting two coupled radiative transfer models (the sandwich and six-flux model). The correlation lengths found are in the same range as those determined in the literature using cold laboratory work. The technique shows great potential in the determination of the snow correlation length under field conditions.

Bone mineral density (BMD) and computer tomographic measurements of the equine proximal phalanx in correlation with breaking strength

Despite the fact that bone mineral density (BMD) is an important fracture risk predictor in human medicine, studies in equine orthopedic research are still lacking. We hypothesized that BMD correlates with bone failure and fatigue fractures of this bone. Thus, the objectives of this study were to measure the structural and mechanical properties of the proximal phalanx with dual energy X-ray absorptiometry (DXA), to correlate the data obtained from DXA and computer tomography (CT) measurements to those obtained by loading pressure examination and to establish representative region of interest (ROI) for in vitro BMD measurements of the equine proximal phalanx for predicting bone failure force. DXA was used to measure the whole bone BMD and additional three ROI sites in 14 equine proximal phalanges. Following evaluation of the bone density, whole bone, cortical width and area in the mid-diaphyseal plane were measured on CT images. Bones were broken using a manually controlled universal bone crusher to measure bone failure force and reevaluated for the site of fractures on follow-up CT images. Compressive load was applied at a constant displacement rate of 2 mm/min until failure, defined as the first clear drop in the load measurement. The lowest BMD was measured at the trabecular region (mean +/- SD: 1.52 +/- 0.12 g/cm2; median: 1.48 g/cm2; range: 1.38-1.83 g/cm2). There was a significant positive linear correlation between trabelcular BMD and the breaking strength (P = 0.023, r = 0.62). The trabecular region of the proximal phalanx appears to be the only significant indicator of failure of strength in vitro. This finding should be reassessed to further reveal the prognostic value of trabecular BMD in an in vivo fracture risk model.

Daily ozone cycle in the stratosphere: global, regional and seasonal behaviour modelled with the Whole Atmosphere Community Climate Model

The Whole Atmosphere Community Climate Model (WACCM) is utilised to study the daily ozone cycle and underlying photochemical and dynamical processes. The analysis is focused on the daily ozone cycle in the middle stratosphere at 5 hPa where satellite-based trend estimates of stratospheric ozone are most biased by diurnal sampling effects and drifting satellite orbits. The simulated ozone cycle shows a minimum after sunrise and a maximum in the late afternoon. Further, a seasonal variation of the daily ozone cycle in the stratosphere was found. Depending on season and latitude, the peak-to-valley difference of the daily ozone cycle varies mostly between 3 and 5% (0.4 ppmv) with respect to the midnight ozone volume mixing ratio. The maximal variation of 15% (0.8 ppmv) is found at the polar circle in summer. The global pattern of the strength of the daily ozone cycle is mainly governed by the solar zenith angle and the sunshine duration. In addition, we find synoptic-scale variations in the strength of the daily ozone cycle. These variations are often anti-correlated to regional temperature anomalies and are due to the temperature dependence of the rate coefficients k2 and k3 of the Chapman cycle reactions. Further, the NOx catalytic cycle counteracts the accumulation of ozone during daytime and leads to an anti-correlation between anomalies in NOx and the strength of the daily ozone cycle. Similarly, ozone recombines with atomic oxygen which leads to an anti-correlation between anomalies in ozone abundance and the strength of the daily ozone cycle. At higher latitudes, an increase of the westerly (easterly) wind cause a decrease (increase) in the sunshine duration of an air parcel leading to a weaker (stronger) daily ozone cycle.

Reconstructing the subsurface ocean decadal variability using surface nudging in a perfect model framework

Initialising the ocean internal variability for decadal predictability studies is a new area of research and a variety of ad hoc methods are currently proposed. In this study, we explore how nudging with sea surface temperature (SST) and salinity (SSS) can reconstruct the threedimensional variability of the ocean in a perfect model framework. This approach builds on the hypothesis that oceanic processes themselves will transport the surface information into the ocean interior as seen in ocean-only simulations. Five nudged simulations are designed to reconstruct a 150 years ‘‘target’’ simulation, defined as a portion of a long control simulation. The nudged simulations differ by the variables restored to, SST or SST + SSS, and by the area where the nudging is applied. The strength of the heat flux feedback is diagnosed from observations and the restoring coefficients for SSS use the same time-scale. We observed that this choice prevents spurious convection at high latitudes and near sea-ice border when nudging both SST and SSS. In the tropics, nudging the SST is enough to reconstruct the tropical atmosphere circulation and the associated dynamical and thermodynamical impacts on the underlying ocean. In the tropical Pacific Ocean, the profiles for temperature show a significant correlation from the surface down to 2,000 m, due to dynamical adjustment of the isopycnals. At mid-tohigh latitudes, SSS nudging is required to reconstruct both the temperature and the salinity below the seasonal thermocline. This is particularly true in the North Atlantic where adding SSS nudging enables to reconstruct the deep convection regions of the target. By initiating a previously documented 20-year cycle of the model, the SST + SSS nudging is also able to reproduce most of the AMOC variations, a key source of decadal predictability. Reconstruction at depth does not significantly improve with amount of time spent nudging and the efficiency of the surface nudging rather depends on the period/events considered. The joint SST + SSS nudging applied verywhere is the most efficient approach. It ensures that the right water masses are formed at the right surface density, the subsequent circulation, subduction and deep convection further transporting them at depth. The results of this study underline the potential key role of SSS for decadal predictability and further make the case for sustained largescale observations of this field.

Use of activity monitors for assessment of pruritus in an acute model of canine atopic dermatitis

BACKGROUND: We developed a canine model of acute atopic dermatitis to evaluate the potential of compounds to treat pruritus and skin lesions induced in Dermatophagoides farinae (Df)-sensitized dogs. HYPOTHESIS/OBJECTIVES: The aim was to investigate the effectiveness of long-term recording activity monitors to assess pruritus induced by allergen challenges. ANIMALS: Thirty-two Df-sensitized laboratory dogs. METHODS: In two blinded crossover studies, 28 Df-sensitized dogs were challenged on 3 days with a Df slurry applied to clipped abdominal skin. Dogs were treated with a positive control (prednisolone 1 mg/kg once daily for 5 days, starting 1 day before challenge) or left untreated; all were fitted with activity monitors. To confirm pruritus, a parallel study with four dogs was conducted, filming the dogs before and during challenge and assessing the film for pruritic behaviour. RESULTS: The activity of dogs treated with prednisolone was significantly lower between 00.00 and 03.00 h and between 03.00 and 06.00 h compared with untreated dogs (repeated-measures ANCOVA; P < 0.0001). To determine whether the recorded night-time activity corresponded to pruritic manifestations, we compared activity monitor and video recordings of four dogs for two periods (16.30-20.30 and 24.00-03.00 h) before and during a Df challenge. The correlation between night-time activity monitor activity and observed pruritic behaviour was highly significant (test of correlation coefficient versus zero: r = 0.57, P < 0.0001). CONCLUSIONS AND CLINICAL IMPORTANCE: Determination of night-time activity with activity monitors after allergen challenge appears to be an objective and practical way to assess pruritus in this experimental model of canine atopic dermatitis.