Densities with Spherical Level Sets in the Gauss-exponential Domain

2000 Mathematics Subject Classification: 60F05, 60B10.

Modeling and Model Predictive Control of an Industrial Hydraulic System

This thesis aims to illustrate the construction of a mathematical model of a hydraulic system, oriented to the design of a model predictive control (MPC) algorithm. The modeling procedure starts with the basic formulation of a piston-servovalve system. The latter is a complex non linear system with some unknown and not measurable effects that constitute a challenging problem for the modeling procedure. The first level of approximation for system parameters is obtained basing on datasheet informations, provided workbench tests and other data from the company. Then, to validate and refine the model, open-loop simulations have been made for data matching with the characteristics obtained from real acquisitions. The final developed set of ODEs captures all the main peculiarities of the system despite some characteristics due to highly varying and unknown hydraulic effects, like the unmodeled resistive elements of the pipes. After an accurate analysis, since the model presents many internal complexities, a simplified version is presented. The latter is used to linearize and discretize correctly the non linear model. Basing on that, a MPC algorithm for reference tracking with linear constraints is implemented. The results obtained show the potential of MPC in this kind of industrial applications, thus a high quality tracking performances while satisfying state and input constraints. The increased robustness and flexibility are evident with respect to the standard control techniques, such as PID controllers, adopted for these systems. The simulations for model validation and the controlled system have been carried out in a Python code environment.

Mercúrio total em cabelo de crianças de uma população costeira, Cananéia, São Paulo, Brasil

Os níveis de mercúrio (Hg) total em cabelos estão diretamente relacionados à alimentação, particularmente ao consumo de peixes por populações costeiras com grande representação caiçara. No presente estudo foram avaliados os níveis de mercúrio total em cabelos de crianças com idade entre 4 e 12 anos, pertencente a três escolas públicas da cidade de Cananéia, São Paulo, Brasil. Os resultados obtidos (mediana e intervalo) para mercúrio total foram de: 0,04mg.kg-1 (0,01-0,77mg.kg-1), 0,39mg.kg-1 (< 0,01-3,33mg.kg-1) e 0,39mg.kg-1 (< 0,01-2,81mg.kg-1) considerando as escolas ES1, ES2 e ES3, respectivamente. Em geral, os valores encontrados estiveram bem abaixo do valor preconizado pela Organização Mundial da Saúde para uma população adulta não exposta ao mercúrio (2,0mg.kg-1). Os baixos valores observados e a inexistência de valores de referência para mercúrio total em cabelos de crianças brasileiras possibilitam a consideração desses valores como possível referência nacional em cabelos de populações costeiras, uma vez que foram obtidos em região de baixo impacto ambiental.

Qualidade nutricional das refeições servidas em uma unidade de alimentação e nutrição de uma indústria da região metropolitana de São Paulo

OBJETIVO: Avaliar a qualidade nutricional das refeições servidas em uma Unidade de Alimentação e Nutrição de uma fábrica da região metropolitana da cidade de São Paulo. MÉTODOS: Dentre os cardápios praticados no período de um ano (242 dias) na unidade mencionada, foram selecionados 30% por sorteio sistemático, e avaliados utilizando-se o Índice de Qualidade da Refeição, com base nas recomendações da Organização Mundial da Saúde e do Ministério da Saúde brasileiro. Esse índice compõe-se de cinco itens que variam entre zero e 20 pontos cada um: adequação na oferta de hortaliças e frutas; oferta de carboidratos; oferta de gordura total; oferta de gordura saturada e variabilidade do cardápio. No período analisado, foram servidas 367 preparações, agrupadas em 30 categorias, segundo composição e forma de preparo. A correlação de Spearman foi utilizada para investigar a correlação do índice com os nutrientes da refeição. As análises foram realizadas no pacote estatístico STATA, considerando-se o nível de significância de 5%. RESULTADOS: O valor médio do Índice de Qualidade da Refeição foi de 64,60 (DP=21,18) pontos, sendo 44% das refeições classificadas como "refeição que necessita de melhora" e apenas 25% como "adequadas". Além do arroz e do feijão, servidos diariamente, as preparações mais frequentes foram: legumes e frutas (30%), massas e cremes (12%), frituras (9%) e sobremesas com creme (8%). Encontrou-se correlação positiva entre o Índice de Qualidade da Refeição e a vitamina C (r=0,32). CONCLUSÃO: Apesar da presença constante de frutas, legumes e verduras, há a necessidade de adequar a oferta das preparações às recomendações para uma alimentação saudável, que efetivamente colaborem na promoção da saúde.

Effects of phonophoresis with Arnica montana onto acute inflammatory process in rat skeletal muscles: An experimental study

This study aimed at verifying the effects of phonophoresis associated with Arnica montana on the acute phase of an inflammatory muscle lesion. Forty Wistar male rats (300 +/- 50 g), of which the Tibialis Anterior muscle was surgically lesioned, were divided into four groups (n = 10 each): control group received no treatment; the ultrasound group (US) was treated in pulsed mode with 1-MHz frequency, 0.5 W/cm(2) intensity (spatial and temporal average - SATA), duty cycle of 1: 2 (2 ms on, 4 ms off, 50%), time of application 3 min per session, one session per day, for 3 days; the phonophoresis or ultrasound plus arnica (US+A) group was treated with arnica with the same US parameters plus arnica gel; and the arnica group (A) was submitted to massage with arnica gel, also for 3 min, once a day, for 3 days. Treatment started 24 h after the surgical lesion. On the 4th day after lesion creation, animals were sacrificed and sections of the lesioned, inflamed muscle were removed for quantitative (mononuclear and polymorphonuclear cell count) and qualitative histological analysis. Collected data from the 4 groups were statistically analyzed and the significance level set at p < 0.05. Results show higher mononuclear cell density in all three treated groups with no significant difference between them, but values were significantly different (p < 0.0001) when compared to control group`s. As to polymorphonuclear cell density, significant differences were found between control group (p = 0.0134) and US, US+A and A groups; the arnica group presented lesser density of polymorphonuclear cells when compared (p = 0.0134) to the other groups. No significant difference was found between US and US+A groups. While the massage with arnica gel proved to be an effective anti-inflammatory on acute muscle lesion in topic use, these results point to ineffectiveness of Arnica montana phonophoresis, US having seemingly checked or minimized its anti-inflammatory effect. (C) 2008 Elsevier B. V. All rights reserved.

Exploring the knowledge contained in neuroimages: Statistical discriminant analysis and automatic segmentation of the most significant changes

Objective: The aim of this article is to propose an integrated framework for extracting and describing patterns of disorders from medical images using a combination of linear discriminant analysis and active contour models. Methods: A multivariate statistical methodology was first used to identify the most discriminating hyperplane separating two groups of images (from healthy controls and patients with schizophrenia) contained in the input data. After this, the present work makes explicit the differences found by the multivariate statistical method by subtracting the discriminant models of controls and patients, weighted by the pooled variance between the two groups. A variational level-set technique was used to segment clusters of these differences. We obtain a label of each anatomical change using the Talairach atlas. Results: In this work all the data was analysed simultaneously rather than assuming a priori regions of interest. As a consequence of this, by using active contour models, we were able to obtain regions of interest that were emergent from the data. The results were evaluated using, as gold standard, well-known facts about the neuroanatomical changes related to schizophrenia. Most of the items in the gold standard was covered in our result set. Conclusions: We argue that such investigation provides a suitable framework for characterising the high complexity of magnetic resonance images in schizophrenia as the results obtained indicate a high sensitivity rate with respect to the gold standard. (C) 2010 Elsevier B.V. All rights reserved.

Pulmonary function testing in quadriplegic subjects

Study design: Cross-sectional study. Objective: Pulmonary functional capacity in 23 Brazilian quadriplegic subjects (ASIA A), aged 30 (9.5) years, weight 66 (10.75) kg, height 176 (7) cm, was investigated at 42 ( 64) months postinjury. Setting: University Hospital-UNICAMP, Campinas, Brazil. Method: Subjects performed forced vital capacity ( FVC) and maximal voluntary ventilation (MVV) tests while seated in their standard wheelchairs. Forced Expired Volume after 1 s (FEV1) and FVC/FEV1 ratio were calculated from these tests. Values obtained were compared to three prediction equations from the literature that are used specifically for spinal cord subjects and include different variables in their formulae, such as age, gender, height, postinjury time and injury level. Data are expressed as median (interquartile interval). Differences between values were demonstrated by median confidence interval with significance level set at a 0.05. Results: Obtained data were statistically different from prediction equation results, with FVC 3.11 ( 0.81), 4.46 (0.28), 4.16 (0.33), 4.26 (0.42); FEV1 2.77 (1.03), 3.67 (0.21), 3.66 (0.30), 3.45 (0.39) and MVV 92 (27), 154.2 (11.9), 156.6 (14),157.3 (16.8), where the first value is obtained experimentally and the second, third and fourth values correspond to predicted values. The results obtained from spirometry test in this study differed significantly from the results obtained when prediction equations were used. Conclusion: The use of prediction equations developed to estimate pulmonary function in wheelchair users significantly overestimates pulmonary function of quadriplegic individuals with complete lesions (ASIA group A), in comparison to measured values.

Human Fatalities from Cyanobacteria: Chemical and Biological Evidence for Cyanotoxins

An outbreak of acute liver failure occurred at a dialysis center in Caruaru, Brazil (8 degrees 17 'S, 35 degrees 58 'W), 134 km from Recife, the state capital of Pernambuco. At the clinic, 116 (89%) of 131 patients experienced visual disturbances, nausea, and vomiting after routine hemodialysis treatment on 13-20 February 1996. Subsequently, 100 patients developed acute liver failure, and of these 76 died. As of December 1996, 52 of the deaths could be attributed to a common syndrome now called Caruaru syndrome. Examination of phytoplankton from the dialysis clinic's water source, analyses of the clinic's water treatment system, plus serum and liver tissue of clinic patients led to the identification of two groups of cyanobacterial toxins, the hepatotoxic cyclic peptide microcystins and the hepatotoxic alkaloid cylindrospermopsin. Comparison of victims' symptoms and pathology using animal studies of these two cyanotoxins leads us to conclude that the major contributing factor to death of the dialyses patients was intravenous exposure to microcystins, specifically microcystin-YR, -LR, and -AR. From liver concentrations and exposure volumes, it was estimated that 19.5 mug/L microcystin was in the water used for dialysis treatments. This is 19.5 times the level set as a guideline for safe drinking water supplies by the World. Health Organization.

Brief note on the variation of constants formula for fuzzy differential equations

This note gives a theory of state transition matrices for linear systems of fuzzy differential equations. This is used to give a fuzzy version of the classical variation of constants formula. A simple example of a time-independent control system is used to illustrate the methods. While similar problems to the crisp case arise for time-dependent systems, in time-independent cases the calculations are elementary solutions of eigenvalue-eigenvector problems. In particular, for nonnegative or nonpositive matrices, the problems at each level set, can easily be solved in MATLAB to give the level sets of the fuzzy solution. (C) 2002 Elsevier Science B.V. All rights reserved.

PAH air pollution at a Portuguese urban area: carcinogenic risks and sources identification

This study aimed to characterize air pollution and the associated carcinogenic risks of polycyclic aromatic hydrocarbon (PAHs) at an urban site, to identify possible emission sources of PAHs using several statistical methodologies, and to analyze the influence of other air pollutants and meteorological variables on PAH concentrations.The air quality and meteorological data were collected in Oporto, the second largest city of Portugal. Eighteen PAHs (the 16 PAHs considered by United States Environment Protection Agency (USEPA) as priority pollutants, dibenzo[a,l]pyrene, and benzo[j]fluoranthene) were collected daily for 24 h in air (gas phase and in particles) during 40 consecutive days in November and December 2008 by constant low-flow samplers and using polytetrafluoroethylene (PTFE) membrane filters for particulate (PM10 and PM2.5 bound) PAHs and pre-cleaned polyurethane foam plugs for gaseous compounds. The other monitored air pollutants were SO2, PM10, NO2, CO, and O3; the meteorological variables were temperature, relative humidity, wind speed, total precipitation, and solar radiation. Benzo[a]pyrene reached a mean concentration of 2.02 ngm−3, surpassing the EU annual limit value. The target carcinogenic risks were equal than the health-based guideline level set by USEPA (10−6) at the studied site, with the cancer risks of eight PAHs reaching senior levels of 9.98×10−7 in PM10 and 1.06×10−6 in air. The applied statistical methods, correlation matrix, cluster analysis, and principal component analysis, were in agreement in the grouping of the PAHs. The groups were formed according to their chemical structure (number of rings), phase distribution, and emission sources. PAH diagnostic ratios were also calculated to evaluate the main emission sources. Diesel vehicular emissions were the major source of PAHs at the studied site. Besides that source, emissions from residential heating and oil refinery were identified to contribute to PAH levels at the respective area. Additionally, principal component regression indicated that SO2, NO2, PM10, CO, and solar radiation had positive correlation with PAHs concentrations, while O3, temperature, relative humidity, and wind speed were negatively correlated.

Numerische Simulation dendritischen Kristallwachstums

Dendritic Growth, Stefan-Problem, Finite-Element-Method, Level-Set-Method

Blood Pressure and Hemodynamic Adaptations after a Training Program in Young Individuals with Down Syndrome

Background:Cardiovascular diseases affect people worldwide. Individuals with Down Syndrome (DS) have an up to sixteen-time greater risk of mortality from cardiovascular diseases.Objective:To evaluate the effects of aerobic and resistance exercises on blood pressure and hemodynamic variables of young individuals with DS.Methods:A total of 29 young individuals with DS participated in the study. They were divided into two groups: aerobic training (AT) (n = 14), and resistance training (TR) (n = 15). Their mean age was 15.7 ± 2.82 years. The training program lasted 12 weeks, and had a frequency of three times a week for AT and twice a week for RT. AT was performed in treadmill/ bicycle ergometer, at an intensity between 50%-70% of the HR reserve. RT comprised nine exercises with three sets of 12 repetition-maximum. Systolic blood pressure (SBP), diastolic blood pressure (DBP), mean blood pressure (MBP) and hemodynamic variables were assessed beat-to-beat using the Finometer device before/after the training program. Descriptive analysis, the Shapiro-Wilk test to check the normality of data, and the two-way ANOVA for repeated measures were used to compare pre- and post-training variables. The Pearson’s correlation coefficient was calculated to correlate hemodynamic variables. The SPSS version 18.0 was used with the significance level set at p < 0.05.Results:After twelve weeks of aerobic and/or resistance training, significant reductions in variables SBP, DBP and MBP were observed.Conclusion:This study suggests a chronic hypotensive effect of moderate aerobic and resistance exercises on young individuals with DS.

Central and Cortical Gray Mater Segmentation of Magnetic Resonance Images of the Fetal Brain

Motivation. The study of human brain development in itsearly stage is today possible thanks to in vivo fetalmagnetic resonance imaging (MRI) techniques. Aquantitative analysis of fetal cortical surfacerepresents a new approach which can be used as a markerof the cerebral maturation (as gyration) and also forstudying central nervous system pathologies [1]. However,this quantitative approach is a major challenge forseveral reasons. First, movement of the fetus inside theamniotic cavity requires very fast MRI sequences tominimize motion artifacts, resulting in a poor spatialresolution and/or lower SNR. Second, due to the ongoingmyelination and cortical maturation, the appearance ofthe developing brain differs very much from thehomogenous tissue types found in adults. Third, due tolow resolution, fetal MR images considerably suffer ofpartial volume (PV) effect, sometimes in large areas.Today extensive efforts are made to deal with thereconstruction of high resolution 3D fetal volumes[2,3,4] to cope with intra-volume motion and low SNR.However, few studies exist related to the automatedsegmentation of MR fetal imaging. [5] and [6] work on thesegmentation of specific areas of the fetal brain such asposterior fossa, brainstem or germinal matrix. Firstattempt for automated brain tissue segmentation has beenpresented in [7] and in our previous work [8]. Bothmethods apply the Expectation-Maximization Markov RandomField (EM-MRF) framework but contrary to [7] we do notneed from any anatomical atlas prior. Data set &Methods. Prenatal MR imaging was performed with a 1-Tsystem (GE Medical Systems, Milwaukee) using single shotfast spin echo (ssFSE) sequences (TR 7000 ms, TE 180 ms,FOV 40 x 40 cm, slice thickness 5.4mm, in plane spatialresolution 1.09mm). Each fetus has 6 axial volumes(around 15 slices per volume), each of them acquired inabout 1 min. Each volume is shifted by 1 mm with respectto the previous one. Gestational age (GA) ranges from 29to 32 weeks. Mother is under sedation. Each volume ismanually segmented to extract fetal brain fromsurrounding maternal tissues. Then, in-homogeneityintensity correction is performed using [9] and linearintensity normalization is performed to have intensityvalues that range from 0 to 255. Note that due tointra-tissue variability of developing brain someintensity variability still remains. For each fetus, ahigh spatial resolution image of isotropic voxel size of1.09 mm is created applying [2] and using B-splines forthe scattered data interpolation [10] (see Fig. 1). Then,basal ganglia (BS) segmentation is performed on thissuper reconstructed volume. Active contour framework witha Level Set (LS) implementation is used. Our LS follows aslightly different formulation from well-known Chan-Vese[11] formulation. In our case, the LS evolves forcing themean of the inside of the curve to be the mean intensityof basal ganglia. Moreover, we add local spatial priorthrough a probabilistic map created by fitting anellipsoid onto the basal ganglia region. Some userinteraction is needed to set the mean intensity of BG(green dots in Fig. 2) and the initial fitting points forthe probabilistic prior map (blue points in Fig. 2). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed as described in [8]. Results.The case study presented here has 29 weeks of GA. Thehigh resolution reconstructed volume is presented in Fig.1. The steps of BG segmentation are shown in Fig. 2.Overlap in comparison with manual segmentation isquantified by the Dice similarity index (DSI) equal to0.829 (values above 0.7 are considered a very goodagreement). Such BG segmentation has been applied on 3other subjects ranging for 29 to 32 GA and the DSI hasbeen of 0.856, 0.794 and 0.785. Our segmentation of theinner (red and blue contours) and outer cortical surface(green contour) is presented in Fig. 3. Finally, torefine the results we include our WM segmentation in theFreesurfer software [12] and some manual corrections toobtain Fig.4. Discussion. Precise cortical surfaceextraction of fetal brain is needed for quantitativestudies of early human brain development. Our workcombines the well known statistical classificationframework with the active contour segmentation forcentral gray mater extraction. A main advantage of thepresented procedure for fetal brain surface extraction isthat we do not include any spatial prior coming fromanatomical atlases. The results presented here arepreliminary but promising. Our efforts are now in testingsuch approach on a wider range of gestational ages thatwe will include in the final version of this work andstudying as well its generalization to different scannersand different type of MRI sequences. References. [1]Guibaud, Prenatal Diagnosis 29(4) (2009). [2] Rousseau,Acad. Rad. 13(9), 2006, [3] Jiang, IEEE TMI 2007. [4]Warfield IADB, MICCAI 2009. [5] Claude, IEEE Trans. Bio.Eng. 51(4) (2004). [6] Habas, MICCAI (Pt. 1) 2008. [7]Bertelsen, ISMRM 2009 [8] Bach Cuadra, IADB, MICCAI 2009.[9] Styner, IEEE TMI 19(39 (2000). [10] Lee, IEEE Trans.Visual. And Comp. Graph. 3(3), 1997, [11] Chan, IEEETrans. Img. Proc, 10(2), 2001 [12] Freesurfer,http://surfer.nmr.mgh.harvard.edu.

Neurally adjusted ventilatory assist (NAVA) improves patient-ventilator interaction during non-invasive ventilation delivered by face mask.

PURPOSE: To determine if, compared to pressure support (PS), neurally adjusted ventilatory assist (NAVA) reduces patient-ventilator asynchrony in intensive care patients undergoing noninvasive ventilation with an oronasal face mask. METHODS: In this prospective interventional study we compared patient-ventilator synchrony between PS (with ventilator settings determined by the clinician) and NAVA (with the level set so as to obtain the same maximal airway pressure as in PS). Two 20-min recordings of airway pressure, flow and electrical activity of the diaphragm during PS and NAVA were acquired in a randomized order. Trigger delay (T(d)), the patient's neural inspiratory time (T(in)), ventilator pressurization duration (T(iv)), inspiratory time in excess (T(iex)), number of asynchrony events per minute and asynchrony index (AI) were determined. RESULTS: The study included 13 patients, six with COPD, and two with mixed pulmonary disease. T(d) was reduced with NAVA: median 35 ms (IQR 31-53 ms) versus 181 ms (122-208 ms); p = 0.0002. NAVA reduced both premature and delayed cyclings in the majority of patients, but not the median T(iex) value. The total number of asynchrony events tended to be reduced with NAVA: 1.0 events/min (0.5-3.1 events/min) versus 4.4 events/min (0.9-12.1 events/min); p = 0.08. AI was lower with NAVA: 4.9 % (2.5-10.5 %) versus 15.8 % (5.5-49.6 %); p = 0.03. During NAVA, there were no ineffective efforts, or late or premature cyclings. PaO(2) and PaCO(2) were not different between ventilatory modes. CONCLUSION: Compared to PS, NAVA improved patient ventilator synchrony during noninvasive ventilation by reducing T(d) and AI. Moreover, with NAVA, ineffective efforts, and late and premature cyclings were absent.

Harmonic active contours.

We propose a segmentation method based on the geometric representation of images as 2-D manifolds embedded in a higher dimensional space. The segmentation is formulated as a minimization problem, where the contours are described by a level set function and the objective functional corresponds to the surface of the image manifold. In this geometric framework, both data-fidelity and regularity terms of the segmentation are represented by a single functional that intrinsically aligns the gradients of the level set function with the gradients of the image and results in a segmentation criterion that exploits the directional information of image gradients to overcome image inhomogeneities and fragmented contours. The proposed formulation combines this robust alignment of gradients with attractive properties of previous methods developed in the same geometric framework: 1) the natural coupling of image channels proposed for anisotropic diffusion and 2) the ability of subjective surfaces to detect weak edges and close fragmented boundaries. The potential of such a geometric approach lies in the general definition of Riemannian manifolds, which naturally generalizes existing segmentation methods (the geodesic active contours, the active contours without edges, and the robust edge integrator) to higher dimensional spaces, non-flat images, and feature spaces. Our experiments show that the proposed technique improves the segmentation of multi-channel images, images subject to inhomogeneities, and images characterized by geometric structures like ridges or valleys.