Steady-state equilibrium phase inversion recovery ON-resonant water suppression (IRON) MR angiography in conjunction with superparamagnetic nanoparticles. A robust technique for imaging within a wide range of contrast agent dosages.

PURPOSE: To investigate the ability of inversion recovery ON-resonant water suppression (IRON) in conjunction with P904 (superparamagnetic nanoparticles which consisting of a maghemite core coated with a low-molecular-weight amino-alcohol derivative of glucose) to perform steady-state equilibrium phase MR angiography (MRA) over a wide dose range. MATERIALS AND METHODS: Experiments were approved by the institutional animal care committee. Rabbits (n = 12) were imaged at baseline and serially after the administration of 10 incremental dosages of 0.57-5.7 mgFe/Kg P904. Conventional T1-weighted and IRON MRA were obtained on a clinical 1.5 Tesla (T) scanner to image the thoracic and abdominal aorta, and peripheral vessels. Contrast-to-noise ratios (CNR) and vessel sharpness were quantified. RESULTS: Using IRON MRA, CNR and vessel sharpness progressively increased with incremental dosages of the contrast agent P904, exhibiting constantly higher contrast values than T1 -weighted MRA over a very wide range of contrast agent doses (CNR of 18.8 ± 5.6 for IRON versus 11.1 ± 2.8 for T1 -weighted MRA at 1.71 mgFe/kg, P = 0.02 and 19.8 ± 5.9 for IRON versus -0.8 ± 1.4 for T1-weighted MRA at 3.99 mgFe/kg, P = 0.0002). Similar results were obtained for vessel sharpness in peripheral vessels, (Vessel sharpness of 46.76 ± 6.48% for IRON versus 33.20 ± 3.53% for T1-weighted MRA at 1.71 mgFe/Kg, P = 0.002, and of 48.66 ± 5.50% for IRON versus 19.00 ± 7.41% for T1-weighted MRA at 3.99 mgFe/Kg, P = 0.003). CONCLUSION: Our study suggests that quantitative CNR and vessel sharpness after the injection of P904 are consistently higher for IRON MRA when compared with conventional T1-weighted MRA. These findings apply for a wide range of contrast agent dosages.

A noninvasive method for measuring the velocity of diffuse hydrothermal flow by tracking moving refractive index anomalies

Diffuse flow velocimetry (DFV) is introduced as a new, noninvasive, optical technique for measuring the velocity of diffuse hydrothermal flow. The technique uses images of a motionless, random medium (e.g.,rocks) obtained through the lens of a moving refraction index anomaly (e.g., a hot upwelling). The method works in two stages. First, the changes in apparent background deformation are calculated using particle image velocimetry (PIV). The deformation vectors are determined by a cross correlation of pixel intensities across consecutive images. Second, the 2-D velocity field is calculated by cross correlating the deformation vectors between consecutive PIV calculations. The accuracy of the method is tested with laboratory and numerical experiments of a laminar, axisymmetric plume in fluids with both constant and temperaturedependent viscosity. Results show that average RMS errors are ∼5%–7% and are most accurate in regions of pervasive apparent background deformation which is commonly encountered in regions of diffuse hydrothermal flow. The method is applied to a 25 s video sequence of diffuse flow from a small fracture captured during the Bathyluck’09 cruise to the Lucky Strike hydrothermal field (September 2009). The velocities of the ∼10°C–15°C effluent reach ∼5.5 cm/s, in strong agreement with previous measurements of diffuse flow. DFV is found to be most accurate for approximately 2‐D flows where background objects have a small spatial scale, such as sand or gravel

A fast Bayesian reconstruction algorithm for emission tomography with entropy prior converging to feasible images

The development and tests of an iterative reconstruction algorithm for emission tomography based on Bayesian statistical concepts are described. The algorithm uses the entropy of the generated image as a prior distribution, can be accelerated by the choice of an exponent, and converges uniformly to feasible images by the choice of one adjustable parameter. A feasible image has been defined as one that is consistent with the initial data (i.e. it is an image that, if truly a source of radiation in a patient, could have generated the initial data by the Poisson process that governs radioactive disintegration). The fundamental ideas of Bayesian reconstruction are discussed, along with the use of an entropy prior with an adjustable contrast parameter, the use of likelihood with data increment parameters as conditional probability, and the development of the new fast maximum a posteriori with entropy (FMAPE) Algorithm by the successive substitution method. It is shown that in the maximum likelihood estimator (MLE) and FMAPE algorithms, the only correct choice of initial image for the iterative procedure in the absence of a priori knowledge about the image configuration is a uniform field.

A new maximum likelihood method for luminosity calibrations

A new statistical parallax method using the Maximum Likelihood principle is presented, allowing the simultaneous determination of a luminosity calibration, kinematic characteristics and spatial distribution of a given sample. This method has been developed for the exploitation of the Hipparcos data and presents several improvements with respect to the previous ones: the effects of the selection of the sample, the observational errors, the galactic rotation and the interstellar absorption are taken into account as an intrinsic part of the formulation (as opposed to external corrections). Furthermore, the method is able to identify and characterize physically distinct groups in inhomogeneous samples, thus avoiding biases due to unidentified components. Moreover, the implementation used by the authors is based on the extensive use of numerical methods, so avoiding the need for simplification of the equations and thus the bias they could introduce. Several examples of application using simulated samples are presented, to be followed by applications to real samples in forthcoming articles.

An Efficient Segmentation Method for Ultrasound Images based on a Semi-supervised Approach and Patch-based Features

Segmenting ultrasound images is a challenging problemwhere standard unsupervised segmentation methods such asthe well-known Chan-Vese method fail. We propose in thispaper an efficient segmentation method for this class ofimages. Our proposed algorithm is based on asemi-supervised approach (user labels) and the use ofimage patches as data features. We also consider thePearson distance between patches, which has been shown tobe robust w.r.t speckle noise present in ultrasoundimages. Our results on phantom and clinical data show avery high similarity agreement with the ground truthprovided by a medical expert.

Modelling of weather radar echoes from anomalous propagation using a hybrid parabolic equation method and NWP model data

Contamination of weather radar echoes by anomalous propagation (anaprop) mechanisms remains a serious issue in quality control of radar precipitation estimates. Although significant progress has been made identifying clutter due to anaprop there is no unique method that solves the question of data reliability without removing genuine data. The work described here relates to the development of a software application that uses a numerical weather prediction (NWP) model to obtain the temperature, humidity and pressure fields to calculate the three dimensional structure of the atmospheric refractive index structure, from which a physically based prediction of the incidence of clutter can be made. This technique can be used in conjunction with existing methods for clutter removal by modifying parameters of detectors or filters according to the physical evidence for anomalous propagation conditions. The parabolic equation method (PEM) is a well established technique for solving the equations for beam propagation in a non-uniformly stratified atmosphere, but although intrinsically very efficient, is not sufficiently fast to be practicable for near real-time modelling of clutter over the entire area observed by a typical weather radar. We demonstrate a fast hybrid PEM technique that is capable of providing acceptable results in conjunction with a high-resolution terrain elevation model, using a standard desktop personal computer. We discuss the performance of the method and approaches for the improvement of the model profiles in the lowest levels of the troposphere.

Analysis of warm season thunderstorms using an object-oriented tracking method based on radar and total lightning data

Monitoring thunderstorms activity is an essential part of operational weather surveillance given their potential hazards, including lightning, hail, heavy rainfall, strong winds or even tornadoes. This study has two main objectives: firstly, the description of a methodology, based on radar and total lightning data to characterise thunderstorms in real-time; secondly, the application of this methodology to 66 thunderstorms that affected Catalonia (NE Spain) in the summer of 2006. An object-oriented tracking procedure is employed, where different observation data types generate four different types of objects (radar 1-km CAPPI reflectivity composites, radar reflectivity volumetric data, cloud-to-ground lightning data and intra-cloud lightning data). In the framework proposed, these objects are the building blocks of a higher level object, the thunderstorm. The methodology is demonstrated with a dataset of thunderstorms whose main characteristics, along the complete life cycle of the convective structures (development, maturity and dissipation), are described statistically. The development and dissipation stages present similar durations in most cases examined. On the contrary, the duration of the maturity phase is much more variable and related to the thunderstorm intensity, defined here in terms of lightning flash rate. Most of the activity of IC and CG flashes is registered in the maturity stage. In the development stage little CG flashes are observed (2% to 5%), while for the dissipation phase is possible to observe a few more CG flashes (10% to 15%). Additionally, a selection of thunderstorms is used to examine general life cycle patterns, obtained from the analysis of normalized (with respect to thunderstorm total duration and maximum value of variables considered) thunderstorm parameters. Among other findings, the study indicates that the normalized duration of the three stages of thunderstorm life cycle is similar in most thunderstorms, with the longest duration corresponding to the maturity stage (approximately 80% of the total time).

A wavelet-based method for the determination of the relative resolution between remotely sensed images

Spatial resolution is a key parameter of all remote sensing satellites and platforms. The nominal spatial resolution of satellites is a well-known characteristic because it is directly related to the area in ground that represents a pixel in the detector. Nevertheless, in practice, the actual resolution of a specific image obtained from a satellite is difficult to know precisely because it depends on many other factors such as atmospheric conditions. However, if one has two or more images of the same region, it is possible to compare their relative resolutions. In this paper, a wavelet-decomposition-based method for the determination of the relative resolution between two remotely sensed images of the same area is proposed. The method can be applied to panchromatic, multispectral, and mixed (one panchromatic and one multispectral) images. As an example, the method was applied to compute the relative resolution between SPOT-3, Landsat-5, and Landsat-7 panchromatic and multispectral images taken under similar as well as under very different conditions. On the other hand, if the true absolute resolution of one of the images of the pair is known, the resolution of the other can be computed. Thus, in the last part of this paper, a spatial calibrator that is designed and constructed to help compute the absolute resolution of a single remotely sensed image is described, and an example of its use is presented.

Electroencephalogram paroxysmal theta characterizes cataplexy in mice and children.

Astute control of brain activity states is critical for adaptive behaviours and survival. In mammals and birds, electroencephalographic recordings reveal alternating states of wakefulness, slow wave sleep and paradoxical sleep (or rapid eye movement sleep). This control is profoundly impaired in narcolepsy with cataplexy, a disease resulting from the loss of orexin/hypocretin neurotransmitter signalling in the brain. Narcolepsy with cataplexy is characterized by irresistible bouts of sleep during the day, sleep fragmentation during the night and episodes of cataplexy, a sudden loss of muscle tone while awake and experiencing emotions. The neural mechanisms underlying cataplexy are unknown, but commonly thought to involve those of rapid eye movement-sleep atonia, and cataplexy typically is considered as a rapid eye movement sleep disorder. Here we reassess cataplexy in hypocretin (Hcrt, also known as orexin) gene knockout mice. Using a novel video/electroencephalogram double-blind scoring method, we show that cataplexy is not a state per se, as believed previously, but a dynamic, multi-phased process involving a reproducible progression of states. A knockout-specific state and a stereotypical paroxysmal event were introduced to account for signals and electroencephalogram spectral characteristics not seen in wild-type littermates. Cataplexy almost invariably started with a brief phase of wake-like electroencephalogram, followed by a phase featuring high-amplitude irregular theta oscillations, defining an activity profile distinct from paradoxical sleep, referred to as cataplexy-associated state and in the course of which 1.5-2 s high-amplitude, highly regular, hypersynchronous paroxysmal theta bursts (∼7 Hz) occurred. In contrast to cataplexy onset, exit from cataplexy did not show a predictable sequence of activities. Altogether, these data contradict the hypothesis that cataplexy is a state similar to paradoxical sleep, even if long cataplexies may evolve into paradoxical sleep. Although not exclusive to overt cataplexy, cataplexy-associated state and hypersynchronous paroxysmal theta activities are highly enriched during cataplexy in hypocretin/orexin knockout mice. Their occurrence in an independent narcolepsy mouse model, the orexin/ataxin 3 transgenic mouse, undergoing loss of orexin neurons, was confirmed. Importantly, we document for the first time similar paroxysmal theta hypersynchronies (∼4 Hz) during cataplexy in narcoleptic children. Lastly, we show by deep recordings in mice that the cataplexy-associated state and hypersynchronous paroxysmal theta activities are independent of hippocampal theta and involve the frontal cortex. Cataplexy hypersynchronous paroxysmal theta bursts may represent medial prefrontal activity, associated in humans and rodents with reward-driven motor impulse, planning and conflict monitoring.

New trapping method to survey for presence of the Etruscan shrew Suncus etruscus, the smallest mammal

The presence of the Etruscan shrew Suncus etruscus is hard to prove where its predator, the barn owl Tyto alba, is absent, because most live traps are not triggered by it. I therefore developed a new trapping method involving a feeding period of 1 week followed by one night of trapping using modified Trip Trap traps. I show here in detail how I caught four Etruscan shrews in 2010 with 24 traps in the Valley of Dora Baltea (Piemonte, Italy). In 2011, another 11 Etruscan shrews were caught in Piemonte and Lombardia, Italy, and Ticino, Switzerland. The proposed new method is useful for establishing the presence of the species.

A selective nociceptive block is not sufficient to prevent central changes after peripheral nerve injury

Background: Providing analgesia without suppressing motor or sensory function is a challenge for regional anesthesia and postoperative pain management. Resiniferatoxin (RTX), an ultrapotent agonist for transient receptor potential subtype-1 (TRPV1) can produce this selective blockade, as TRPV1 is selectively expressed on nociceptors. Futhermore, after peripheral nerve injury, spontaneous ectopic activity arises from all types of nerve fibers that can affect spinal neurons and glial cells. The goal of the present experiment is to determine whether spontaneous activity generated in C-fibers or in both A&C-fibers is required for microglia activation. Method: We applied RTX (0.01%) or bupivacaine microspheres to the sciatic nerve of rats to block the conduction of C-fibers or A&C-fibers, respectively, before spared nerve injury (SNI). Behavior was tested and all the rats were sacrificed 2 days later; immunohistochemistry was performed on their spinal cord for mitogen-activated protein kinase (MAPK) p38, bromodeoxyuridine (BrdU, marker of proliferation) and Iba1 (microglial marker). Result: At day 2 after SNI robust mechanical allodynia and p38 activation in spinal microglia were documented. There was also a substantial cell proliferation in the spinal cord, all proliferating cells (BrdU+) being microglia (Iba1+). RTX blocked heat sensitivity and produced heat hypoalgesia without affecting mechanical allodynia and motor function. Microglial proliferation and p38 activation in the spinal cord were not affected by RTX (p >0.05). In contrast, a complete sensory and motor blockade was seen with bupivacaine which also significantly inhibited p38 activation and microglial proliferation in the spinal cord (p <0.05). Conclusion: We conclude that (1) RTX can provide a selective nociceptive blockade but that (2) blocking only nociceptive fibers does not impair the development of mechanical allodynia and microglia activation. Therefore (3) if microglia activation is important for chronic pain development then specific nociceptive blockade won't be sufficient to prevent it.

A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

Low income, community poverty and risk of end stage renal disease

BACKGROUND: The risk of end stage renal disease (ESRD) is increased among individuals with low income and in low income communities. However, few studies have examined the relation of both individual and community socioeconomic status (SES) with incident ESRD. METHODS: Among 23,314 U.S. adults in the population-based Reasons for Geographic and Racial Differences in Stroke study, we assessed participant differences across geospatially-linked categories of county poverty [outlier poverty, extremely high poverty, very high poverty, high poverty, neither (reference), high affluence and outlier affluence]. Multivariable Cox proportional hazards models were used to examine associations of annual household income and geospatially-linked county poverty measures with incident ESRD, while accounting for death as a competing event using the Fine and Gray method. RESULTS: There were 158 ESRD cases during follow-up. Incident ESRD rates were 178.8 per 100,000 person-years (105 py) in high poverty outlier counties and were 76.3 /105 py in affluent outlier counties, p trend = 0.06. In unadjusted competing risk models, persons residing in high poverty outlier counties had higher incidence of ESRD (which was not statistically significant) when compared to those persons residing in counties with neither high poverty nor affluence [hazard ratio (HR) 1.54, 95% Confidence Interval (CI) 0.75-3.20]. This association was markedly attenuated following adjustment for socio-demographic factors (age, sex, race, education, and income); HR 0.96, 95% CI 0.46-2.00. However, in the same adjusted model, income was independently associated with risk of ESRD [HR 3.75, 95% CI 1.62-8.64, comparing the < $20,000 income group to the > $75,000 group]. There were no statistically significant associations of county measures of poverty with incident ESRD, and no evidence of effect modification. CONCLUSIONS: In contrast to annual family income, geospatially-linked measures of county poverty have little relation with risk of ESRD. Efforts to mitigate socioeconomic disparities in kidney disease may be best appropriated at the individual level.

A new method to measure rolling resistance in treadmill cycling.

The purpose of this study was to evaluate a new method of measuring rolling resistance in treadmill cycling and to establish its sensitivity and reproducibility. One participant was asked to keep a bicycle in equilibrium on a treadmill without pedalling at a constant speed of 5.56 m x s(-1), which was held in place in the front by a dynamometer. For each condition, the method consisted of 11 measurements of the force required to hold the cycle at different treadmill slopes (0-10%, increment 1%). The coefficient of rolling resistance was calculated based on the forces applied to the bicycle in equilibrium. To test the sensitivity of the method, the bicycle was successively equipped with three tyre types (700 x 28, 700 x 23, 700 x 22) and inflation pressure was set at 150, 300, 600, 900, and 1100 kPa. To test the reproducibility of the method, a second experimenter repeated all measurements done with the 700 x 23 tyres. The method was sensitive enough to detect an effect of both tyre type and inflation pressure (P < 0.001: two-way ANOVA). The measurement of the coefficient of rolling resistance by two separate experimenters resulted in a small bias of 0.00029 (95% CI, -0.00011 to 0.00068). In conclusion, the new method is sensitive and reliable, as well as being simple and affordable.