A nonparametric ensemble transform method for Bayesian inference

Many applications, such as intermittent data assimilation, lead to a recursive application of Bayesian inference within a Monte Carlo context. Popular data assimilation algorithms include sequential Monte Carlo methods and ensemble Kalman filters (EnKFs). These methods differ in the way Bayesian inference is implemented. Sequential Monte Carlo methods rely on importance sampling combined with a resampling step, while EnKFs utilize a linear transformation of Monte Carlo samples based on the classic Kalman filter. While EnKFs have proven to be quite robust even for small ensemble sizes, they are not consistent since their derivation relies on a linear regression ansatz. In this paper, we propose another transform method, which does not rely on any a priori assumptions on the underlying prior and posterior distributions. The new method is based on solving an optimal transportation problem for discrete random variables. © 2013, Society for Industrial and Applied Mathematics

Aerosol inorganic composition at a tropical site: Discrepancies between filter-based sampling and a semi-continuous method

The concentrations of the water-soluble inorganic aerosol species, ammonium (NH4+), nitrate (NO3-), chloride (Cl-), and sulfate (SO42-), were measured from September to November 2002 at a pasture site in the Amazon Basin (Rondnia, Brazil) (LBA-SMOCC). Measurements were conducted using a semi-continuous technique (Wet-annular denuder/Steam-Jet Aerosol Collector: WAD/SJAC) and three integrating filter-based methods, namely (1) a denuder-filter pack (DFP: Teflon and impregnated Whatman filters), (2) a stacked-filter unit (SFU: polycarbonate filters), and (3) a High Volume dichotomous sampler (HiVol: quartz fiber filters). Measurements covered the late dry season (biomass burning), a transition period, and the onset of the wet season (clean conditions). Analyses of the particles collected on filters were performed using ion chromatography (IC) and Particle-Induced X-ray Emission spectrometry (PIXE). Season-dependent discrepancies were observed between the WAD/SJAC system and the filter-based samplers. During the dry season, when PM2.5 (D-p <= 2.5 mu m) concentrations were similar to 100 mu g m(-3), aerosol NH4+ and SO42- measured by the filter-based samplers were on average two times higher than those determined by the WAD/SJAC. Concentrations of aerosol NO3- and Cl- measured with the HiVol during daytime, and with the DFP during day- and nighttime also exceeded those of the WAD/SJAC by a factor of two. In contrast, aerosol NO3- and Cl- measured with the SFU during the dry season were nearly two times lower than those measured by the WAD/SJAC. These differences declined markedly during the transition period and towards the cleaner conditions during the onset of the wet season (PM2.5 similar to 5 mu g m(-3)); when filter-based samplers measured on average 40-90% less than the WAD/SJAC. The differences were not due to consistent systematic biases of the analytical techniques, but were apparently a result of prevailing environmental conditions and different sampling procedures. For the transition period and wet season, the significance of our results is reduced by a low number of data points. We argue that the observed differences are mainly attributable to (a) positive and negative filter sampling artifacts, (b) presence of organic compounds and organosulfates on filter substrates, and (c) a SJAC sampling efficiency of less than 100%.

A non-parametric vessel detection method for complex vascular structures

Modern medical imaging techniques enable the acquisition of in vivo high resolution images of the vascular system. Most common methods for the detection of vessels in these images, such as multiscale Hessian-based operators and matched filters, rely on the assumption that at each voxel there is a single cylinder. Such an assumption is clearly violated at the multitude of branching points that are easily observed in all, but the Most focused vascular image studies. In this paper, we propose a novel method for detecting vessels in medical images that relaxes this single cylinder assumption. We directly exploit local neighborhood intensities and extract characteristics of the local intensity profile (in a spherical polar coordinate system) which we term as the polar neighborhood intensity profile. We present a new method to capture the common properties shared by polar neighborhood intensity profiles for all the types of vascular points belonging to the vascular system. The new method enables us to detect vessels even near complex extreme points, including branching points. Our method demonstrates improved performance over standard methods on both 2D synthetic images and 3D animal and clinical vascular images, particularly close to vessel branching regions. (C) 2008 Elsevier B.V. All rights reserved.

Design and optimization of LNA topologies with image rejection filters

An important problem in designing RFIC in CMOS technology is the parasitic elements of passive and active devices that complicate design calculations. This article presents three LNA topologies including cascode, folded cascade, and differential cascode and then introduces image rejection filters for low-side and high-side injection. Then, a new method for design and optimization of the circuits based on a Pareto-based multiobjective genetic algorithm is proposed. A set of optimum device values and dimensions that best match design specifications are obtained. The optimization method is layout aware, parasitic aware, and simulation based. Circuit simulations are carried out based on TSMC 0.18 um CMOS technology by using Hspice.

Achieving illumination invariance using image filters

In this chapter we described a novel framework for automatic face recognition in the presence of varying illumination, primarily applicable to matching face sets or sequences. The framework is based on simple image processing filters that compete with unprocessed greyscale input to yield a single matching score between individuals. By performing all numerically consuming computation offline, our method both (i) retains the matching efficiency of simple image filters, but (ii) with a greatly increased robustness, as all online processing is performed in closed-form. Evaluated on a large, real-world data corpus, the proposed framework was shown to be successful in video-based recognition across a wide range of illumination, pose and face motion pattern changes

Hydrograph separation and development of empirical relationships using single-parameter digital filters

A reliable hydrograph separation method is necessary for surface runoff modeling and hydrological studies. This paper investigates and compares the separation characteristics of two single-parameter digital filters, which are herein referred to as the one-parameter algorithm and the conceptual method. The application of the one-parameter algorithm was found to be restricted to low and medium baseflow separations, with a maximum separation limit of 50% of the total runoff hydrograph. The one-parameter algorithm was also observed to produce unrealistic sharp peaks under the peaks of the measured hydrograph when recession constant is smaller than 0.96. On the other hand, the conceptual method is applicable even for catchments fed largely by groundwater discharge. However, a reliable estimation of recession constant is a prerequisite for applying the conceptual method for large baseflow separations. Based on the hydrograph separation results, useful empirical relationships were developed for a partially urbanized watershed to estimate total runoff and direct runoff from the measured rainfall depth. The relationships between rainfall depth and total runoff depth and rainfall depth and direct runoff depth were found to be well represented by linear equations. The empirical relationships were then applied to estimate the long-term contribution of baseflow and surface runoff to total runoff at the study site. Baseflow was found to contribute about 58–61% of the annual total runoff.

Generation of total runoff hydrographs using a method derived from a digital filter algorithm

A reliable prediction of the total runoff hydrograph is necessary for water resources management. This study investigates two approaches to generate total runoff hydrograph by adding baseflow to direct runoff hydrographs. The first approach uses a method, derived from a digital filter algorithm for hydrograph separation, to generate baseflow hydrographs from direct runoff hydrographs. The method appears to perform well in producing the overall shape of the total runoff hydrographs and the acceptable mass balance errors for a year of water cycle. For application, the recession baseflow constant needs to be estimated reliably and the initial baseflow could be approximated to the long-term mean dry weather flow. The second approach assumes a constant baseflow rate. Although this approach is still capable of producing the overall hydrograph shape, it yields high mass balance errors in the total runoff hydrographs for both monthly and long-term periods. Further analysis shows that two-third of the mass balance errors are contributed from periods with direct runoff, implying that the constant baseflow assumption could introduce significant errors into the computations of total runoff hydrograph

A new fluorescence method for determination of ozone in ambient air

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Measurements of ambient ozone using indigo blue-coated filters

Ozone monitoring techniques utilize expensive instruments that are often large and heavy. These instruments are not easy to handle in the field, and their size also limits some sampling schemes, principally for indoor ozone determination. We have developed a lightweight, inexpensive, and sensitive method that offers flexibility to undertake measurements of ambient ozone in many environments, both indoor and outdoor. The method is based on the reaction of ozone with indigo blue dye. The indigo molecule contains 1 carbon double bond (C = C) that reacts with ozone and results in nearly colorless reaction products. During sample collection, 2 cellulose filters coated with 40 mu L of 1.0 x 10(-3) M indigo blue were used. The determinations were done spectrophotometrically at 250 and 600 nm. The analytical parameters studied were sampling time and flow rate. Analytical curves were constructed with concentrations ranging from 37 to 123 parts per billion by volume (ppbv) of standard ozone, at 0.4 L/min and 15 min sampling time. The detection limits achieved were 6 and 9 ppbv, respectively, at 250 and 600 nm. Considering interferences, measurements made at 250 nm gave more reliable and specific values for ozone.

A method for the digital image analysis of ceramic grains based on shape factor segmentation

A very simple and robust method for ceramics grains quantitative image analysis is presented. Based on the use of optimal imaging conditions for reflective light microscopy of bulk samples, a digital image processing routine was developed for shading correction, noise suppressing and contours enhancement. Image analysis was done for grains selected according to their concavities, evaluated by perimeter ratio shape factor, to avoid consider the effects of breakouts and ghost boundaries due to ceramographic preparation limitations. As an example, the method was applied for two ceramics, to compare grain size and morphology distributions. In this case, most of artefacts introduced by ceramographic preparation could be discarded due to the use of perimeter ratio exclusion range.

A neural networks-based method for single-phase harmonic content identification

A neural method is presented in this paper to identify the harmonic components of an ac controller. The components are identified by analyzing the single-phase current waveform. The method effectiveness is verified by applying it to an active power filter (APF) model dedicated to the selective harmonic compensation. Simulation results using theoretical and experimental data are presented to validate the proposed approach. © 2008 IEEE.

UV filters in water samples: Experimental design on the SPE optimization followed by GC-MS/MS analysis

Esse estudo descreve o desenvolvimento e otimização de um método de extração em fase solida (SPE) para análise dos filtros ultravioletas (UV): benzofenona-3 (BP-3), etilhexil salicilato (ES), etilhexil metoxinamato (EHMC) e octocrileno (OC) em matrizes ambientais. Um planejamento fatorial fracionário (PFF) 25-1 foi empregado na avaliação das variáveis significativas do método de extração. As condições experimentais otimizadas da avaliação estatística foram: capacidade do cartucho de 500 mL, eluente acetato de etila, metanol como solvente de lavagem (10% em água, v/v) and volume do eluente de 3 × 2 mL e pH 3. Os parâmetros analíticos avaliados foram satisfatõrios, apresentando linearidade de 100 a 4000 ng L -1, recuperaç ões para os quatro níveis de fortificação (Limite de Quantificação do Método, 200, 1000 e 2000 ng L-1) entre 62 e 107% com desvio padrão relativo menor que 14%. Os limites de quantificação foram encontrados na faixa de ng L-1, variando entre 10 e 100 ng L-1. O método proposto foi aplicado para a determinação dos quatro filtros UV em amostras de águas naturais. This study describes the development and optimization of a solid-phase extraction (SPE) method for analysis of ultraviolet (UV) filters, benzophenone-3 (BP-3), ethylhexyl methoxycinnamate (EHMC), ethylhexyl salicylate (ES) and octocrylene (OC), in environmental matrices. A 25-1 fractional factorial design (FFD) was used to evaluate the significant variables for the extraction method. The optimized experimental conditions determined from the statistical evaluation were: breakthrough volume of 500 mL, eluent of ethyl acetate, wash solvent of methanol (10% in water, v/v), eluent volume of 3 × 2 mL and pH 3. The evaluated analytical parameters were satisfactory for the analytes and showed linearity between 100 and 4000 ng L-1, recoveries for four fortification levels (Method Quantification Limit, 200, 1000 and 2000 ng L-1) were between 62 and 107% with relative standard deviations less than 14%. Limits of quantification were in the ng L-1 range and were between 10 and 100 ng L-1. The proposed method was used to analyze four UV filters in natural water samples. ©2013 Sociedade Brasileira de Química.

Processing of high resolution magic angle spinning spectra of breast cancer cells by the filter diagonalization method

Proton nuclear magnetic resonance (H-1 NMR) spectroscopy for detection of biochemical changes in biological samples is a successful technique. However, the achieved NMR resolution is not sufficiently high when the analysis is performed with intact cells. To improve spectral resolution, high resolution magic angle spinning (HR-MAS) is used and the broad signals are separated by a T-2 filter based on the CPMG pulse sequence. Additionally, HR-MAS experiments with a T-2 filter are preceded by a water suppression procedure. The goal of this work is to demonstrate that the experimental procedures of water suppression and T-2 or diffusing filters are unnecessary steps when the filter diagonalization method (FDM) is used to process the time domain HR-MAS signals. Manipulation of the FDM results, represented as a tabular list of peak positions, widths, amplitudes and phases, allows the removal of water signals without the disturbing overlapping or nearby signals. Additionally, the FDM can also be used for phase correction and noise suppression, and to discriminate between sharp and broad lines. Results demonstrate the applicability of the FDM post-acquisition processing to obtain high quality HR-MAS spectra of heterogeneous biological materials.

Vibration as a method for muscular stimulation

Human reactions to vibration have been extensively investigated in the past. Vibration, as well as whole-body vibration (WBV), has been commonly considered as an occupational hazard for its detrimental effects on human condition and comfort. Although long term exposure to vibrations may produce undesirable side-effects, a great part of the literature is dedicated to the positive effects of WBV when used as method for muscular stimulation and as an exercise intervention. Whole body vibration training (WBVT) aims to mechanically activate muscles by eliciting neuromuscular activity (muscle reflexes) via the use of vibrations delivered to the whole body. The most mentioned mechanism to explain the neuromuscular outcomes of vibration is the elicited neuromuscular activation. Local tendon vibrations induce activity of the muscle spindle Ia fibers, mediated by monosynaptic and polysynaptic pathways: a reflex muscle contraction known as the Tonic Vibration Reflex (TVR) arises in response to such vibratory stimulus. In WBVT mechanical vibrations, in a range from 10 to 80 Hz and peak to peak displacements from 1 to 10 mm, are usually transmitted to the patient body by the use of oscillating platforms. Vibrations are then transferred from the platform to a specific muscle group through the subject body. To customize WBV treatments, surface electromyography (SEMG) signals are often used to reveal the best stimulation frequency for each subject. Use of SEMG concise parameters, such as root mean square values of the recordings, is also a common practice; frequently a preliminary session can take place in order to discover the more appropriate stimulation frequency. Soft tissues act as wobbling masses vibrating in a damped manner in response to mechanical excitation; Muscle Tuning hypothesis suggest that neuromuscular system works to damp the soft tissue oscillation that occurs in response to vibrations; muscles alters their activity to dampen the vibrations, preventing any resonance phenomenon. Muscle response to vibration is however a complex phenomenon as it depends on different parameters, like muscle-tension, muscle or segment-stiffness, amplitude and frequency of the mechanical vibration. Additionally, while in the TVR study the applied vibratory stimulus and the muscle conditions are completely characterised (a known vibration source is applied directly to a stretched/shortened muscle or tendon), in WBV study only the stimulus applied to a distal part of the body is known. Moreover, mechanical response changes in relation to the posture. The transmissibility of vibratory stimulus along the body segment strongly depends on the position held by the subject. The aim of this work was the investigation on the effects that the use of vibrations, in particular the effects of whole body vibrations, may have on muscular activity. A new approach to discover the more appropriate stimulus frequency, by the use of accelerometers, was also explored. Different subjects, not affected by any known neurological or musculoskeletal disorders, were voluntarily involved in the study and gave their informed, written consent to participate. The device used to deliver vibration to the subjects was a vibrating platform. Vibrations impressed by the platform were exclusively vertical; platform displacement was sinusoidal with an intensity (peak-to-peak displacement) set to 1.2 mm and with a frequency ranging from 10 to 80 Hz. All the subjects familiarized with the device and the proper positioning. Two different posture were explored in this study: position 1 - hack squat; position 2 - subject standing on toes with heels raised. SEMG signals from the Rectus Femoris (RF), Vastus Lateralis (VL) and Vastus medialis (VM) were recorded. SEMG signals were amplified using a multi-channel, isolated biomedical signal amplifier The gain was set to 1000 V/V and a band pass filter (-3dB frequency 10 - 500 Hz) was applied; no notch filters were used to suppress line interference. Tiny and lightweight (less than 10 g) three-axial MEMS accelerometers (Freescale semiconductors) were used to measure accelerations of onto patient’s skin, at EMG electrodes level. Accelerations signals provided information related to individuals’ RF, Biceps Femoris (BF) and Gastrocnemius Lateralis (GL) muscle belly oscillation; they were pre-processed in order to exclude influence of gravity. As demonstrated by our results, vibrations generate peculiar, not negligible motion artifact on skin electrodes. Artifact amplitude is generally unpredictable; it appeared in all the quadriceps muscles analysed, but in different amounts. Artifact harmonics extend throughout the EMG spectrum, making classic high-pass filters ineffective; however, their contribution was easy to filter out from the raw EMG signal with a series of sharp notch filters centred at the vibration frequency and its superior harmonics (1.5 Hz wide). However, use of these simple filters prevents the revelation of EMG power potential variation in the mentioned filtered bands. Moreover our experience suggests that the possibility of reducing motion artefact, by using particular electrodes and by accurately preparing the subject’s skin, is not easily viable; even though some small improvements were obtained, it was not possible to substantially decrease the artifact. Anyway, getting rid of those artifacts lead to some true EMG signal loss. Nevertheless, our preliminary results suggest that the use of notch filters at vibration frequency and its harmonics is suitable for motion artifacts filtering. In RF SEMG recordings during vibratory stimulation only a little EMG power increment should be contained in the mentioned filtered bands due to synchronous electromyographic activity of the muscle. Moreover, it is better to remove the artifact that, in our experience, was found to be more than 40% of the total signal power. In summary, many variables have to be taken into account: in addition to amplitude, frequency and duration of vibration treatment, other fundamental variables were found to be subject anatomy, individual physiological condition and subject’s positioning on the platform. Studies on WBV treatments that include surface EMG analysis to asses muscular activity during vibratory stimulation should take into account the presence of motion artifacts. Appropriate filtering of artifacts, to reveal the actual effect on muscle contraction elicited by vibration stimulus, is mandatory. However as a result of our preliminary study, a simple multi-band notch filtering may help to reduce randomness of the results. Muscle tuning hypothesis seemed to be confirmed. Our results suggested that the effects of WBV are linked to the actual muscle motion (displacement). The greater was the muscle belly displacement the higher was found the muscle activity. The maximum muscle activity has been found in correspondence with the local mechanical resonance, suggesting a more effective stimulation at the specific system resonance frequency. Holding the hypothesis that muscle activation is proportional to muscle displacement, treatment optimization could be obtained by simply monitoring local acceleration (resonance). However, our study revealed some short term effects of vibratory stimulus; prolonged studies should be assembled in order to consider the long term effectiveness of these results. Since local stimulus depends on the kinematic chain involved, WBV muscle stimulation has to take into account the transmissibility of the stimulus along the body segment in order to ensure that vibratory stimulation effectively reaches the target muscle. Combination of local resonance and muscle response should also be further investigated to prevent hazards to individuals undergoing WBV treatments.

Pascal Filters

We introduce a new type of filter approximation method and call it the Pascal filter, which we construct from the Pascal polynomials. The roll-off characteristics of the Pascal, Butterworth, and the Chebyshev filters are compared.