Border figure detection using a phase oscillator network with dynamical coupling

Oscillator networks have been developed in order to perform specific tasks related to image processing. Here we analytically investigate the existence of synchronism in a pair of phase oscillators that are short-range dynamically coupled. Then, we use these analytical results to design a network able of detecting border of black-and-white figures. Each unit composing this network is a pair of such phase oscillators and is assigned to a pixel in the image. The couplings among the units forming the network are also dynamical. Border detection emerges from the network activity.

First stars XII. Abundances in extremely metal-poor turnoff stars, and comparison with the giants

Context. The detailed chemical abundances of extremely metal-poor (EMP) stars are key guides to understanding the early chemical evolution of the Galaxy. Most existing data, however, treat giant stars that may have experienced internal mixing later. Aims. We aim to compare the results for giants with new, accurate abundances for all observable elements in 18 EMP turno. stars. Methods. VLT/UVES spectra at R similar to 45 000 and S/N similar to 130 per pixel (lambda lambda 330-1000 nm) are analysed with OSMARCS model atmospheres and the TURBOSPECTRUM code to derive abundances for C, Mg, Si, Ca, Sc, Ti, Cr, Mn, Co, Ni, Zn, Sr, and Ba. Results. For Ca, Ni, Sr, and Ba, we find excellent consistency with our earlier sample of EMP giants, at all metallicities. However, our abundances of C, Sc, Ti, Cr, Mn and Co are similar to 0.2 dex larger than in giants of similar metallicity. Mg and Si abundances are similar to 0.2 dex lower (the giant [Mg/Fe] values are slightly revised), while Zn is again similar to 0.4 dex higher than in giants of similar [Fe/H] (6 stars only). Conclusions. For C, the dwarf/giant discrepancy could possibly have an astrophysical cause, but for the other elements it must arise from shortcomings in the analysis. Approximate computations of granulation (3D) effects yield smaller corrections for giants than for dwarfs, but suggest that this is an unlikely explanation, except perhaps for C, Cr, and Mn. NLTE computations for Na and Al provide consistent abundances between dwarfs and giants, unlike the LTE results, and would be highly desirable for the other discrepant elements as well. Meanwhile, we recommend using the giant abundances as reference data for Galactic chemical evolution models.

The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 2: Model sensitivity to the biomass burning

We describe an estimation technique for biomass burning emissions in South America based on a combination of remote-sensing fire products and field observations, the Brazilian Biomass Burning Emission Model (3BEM). For each fire pixel detected by remote sensing, the mass of the emitted tracer is calculated based on field observations of fire properties related to the type of vegetation burning. The burnt area is estimated from the instantaneous fire size retrieved by remote sensing, when available, or from statistical properties of the burn scars. The sources are then spatially and temporally distributed and assimilated daily by the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) in order to perform the prognosis of related tracer concentrations. Three other biomass burning inventories, including GFEDv2 and EDGAR, are simultaneously used to compare the emission strength in terms of the resultant tracer distribution. We also assess the effect of using the daily time resolution of fire emissions by including runs with monthly-averaged emissions. We evaluate the performance of the model using the different emission estimation techniques by comparing the model results with direct measurements of carbon monoxide both near-surface and airborne, as well as remote sensing derived products. The model results obtained using the 3BEM methodology of estimation introduced in this paper show relatively good agreement with the direct measurements and MOPITT data product, suggesting the reliability of the model at local to regional scales.

Estimation of the thickness and the optical parameters of several stacked thin films using optimization

The reverse engineering problem addressed in the present research consists of estimating the thicknesses and the optical constants of two thin films deposited on a transparent substrate using only transmittance data through the whole stack. No functional dispersion relation assumptions are made on the complex refractive index. Instead, minimal physical constraints are employed, as in previous works of some of the authors where only one film was considered in the retrieval algorithm. To our knowledge this is the first report on the retrieval of the optical constants and the thickness of multiple film structures using only transmittance data that does not make use of dispersion relations. The same methodology may be used if the available data correspond to normal reflectance. The software used in this work is freely available through the PUMA Project web page (http://www.ime.usp.br/similar to egbirgin/puma/). (C) 2008 Optical Society of America

Improvements on ICA mixture models for image pre-processing and segmentation

Today several different unsupervised classification algorithms are commonly used to cluster similar patterns in a data set based only on its statistical properties. Specially in image data applications, self-organizing methods for unsupervised classification have been successfully applied for clustering pixels or group of pixels in order to perform segmentation tasks. The first important contribution of this paper refers to the development of a self-organizing method for data classification, named Enhanced Independent Component Analysis Mixture Model (EICAMM), which was built by proposing some modifications in the Independent Component Analysis Mixture Model (ICAMM). Such improvements were proposed by considering some of the model limitations as well as by analyzing how it should be improved in order to become more efficient. Moreover, a pre-processing methodology was also proposed, which is based on combining the Sparse Code Shrinkage (SCS) for image denoising and the Sobel edge detector. In the experiments of this work, the EICAMM and other self-organizing models were applied for segmenting images in their original and pre-processed versions. A comparative analysis showed satisfactory and competitive image segmentation results obtained by the proposals presented herein. (C) 2008 Published by Elsevier B.V.

Registration of temporal sequences of coronal and sagittal MR images through respiratory patterns

This work discusses the determination of the breathing patterns in time sequence of images obtained from magnetic resonance (MR) and their use in the temporal registration of coronal and sagittal images. The registration is made without the use of any triggering information and any special gas to enhance the contrast. The temporal sequences of images are acquired in free breathing. The real movement of the lung has never been seen directly, as it is totally dependent on its surrounding muscles and collapses without them. The visualization of the lung in motion is an actual topic of research in medicine. The lung movement is not periodic and it is susceptible to variations in the degree of respiration. Compared to computerized tomography (CT), MR imaging involves longer acquisition times and it is preferable because it does not involve radiation. As coronal and sagittal sequences of images are orthogonal to each other, their intersection corresponds to a segment in the three-dimensional space. The registration is based on the analysis of this intersection segment. A time sequence of this intersection segment can be stacked, defining a two-dimension spatio-temporal (2DST) image. The algorithm proposed in this work can detect asynchronous movements of the internal lung structures and lung surrounding organs. It is assumed that the diaphragmatic movement is the principal movement and all the lung structures move almost synchronously. The synchronization is performed through a pattern named respiratory function. This pattern is obtained by processing a 2DST image. An interval Hough transform algorithm searches for synchronized movements with the respiratory function. A greedy active contour algorithm adjusts small discrepancies originated by asynchronous movements in the respiratory patterns. The output is a set of respiratory patterns. Finally, the composition of coronal and sagittal image pairs that are in the same breathing phase is realized by comparing of respiratory patterns originated from diaphragmatic and upper boundary surfaces. When available, the respiratory patterns associated to lung internal structures are also used. The results of the proposed method are compared with the pixel-by-pixel comparison method. The proposed method increases the number of registered pairs representing composed images and allows an easy check of the breathing phase. (C) 2010 Elsevier Ltd. All rights reserved.

Compact tiled display with uniform illumination

The demand for more pixels is beginning to be met as manufacturers increase the native resolution of projector chips. Tiling several projectors still offers a solution to augment the pixel capacity of a display. However, problems of color and illumination uniformity across projectors need to be addressed as well as the computer software required to drive such devices. We present the results obtained on a desktop-size tiled projector array of three D-ILA projectors sharing a common illumination source. A short throw lens (0.8:1) on each projector yields a 21-in. diagonal for each image tile; the composite image on a 3×1 array is 3840×1024 pixels with a resolution of about 80 dpi. The system preserves desktop resolution, is compact, and can fit in a normal room or laboratory. The projectors are mounted on precision six-axis positioners, which allow pixel level alignment. A fiber optic beamsplitting system and a single set of red, green, and blue dichroic filters are the key to color and illumination uniformity. The D-ILA chips inside each projector can be adjusted separately to set or change characteristics such as contrast, brightness, or gamma curves. The projectors were then matched carefully: photometric variations were corrected, leading to a seamless image. Photometric measurements were performed to characterize the display and are reported here. This system is driven by a small PC cluster fitted with graphics cards and running Linux. It can be scaled to accommodate an array of 2×3 or 3×3 projectors, thus increasing the number of pixels of the final image. Finally, we present current uses of the display in fields such as astrophysics and archaeology (remote sensing).

The Reflectance Signature of Canopy Components: Implications for the Interpretation of Remotely Sensed Images

The reflectance signatures of plantation pine canopy and understorey components were measured using a spectro-radiometer. The aim was to establish whether differences observed in the reflectance signature of stressed and unstressed pine needles were consistent with observed differences in the reflectance of multispectral Landsat Thematic Mapper (TM) images of healthy and stressed forest. Because overall scene reflectance includes the contribution of each scene component, needle reflectance may not be representative of canopy reflectance. In this investigation, a limited dataset of reflectance signatures from stressed and unstressed needles confirmed the negative relationship between pine needle health and reflectance which was observed in visible red wavelengths. However, the reflectance contribution from bushes, pine needle litter and bare soil tended to reinforce this relationship suggesting that in this instance, overall scene reflectance is comprised of the proportional reflectance of each scene component. In near infrared wavelengths, differences between healthy and stressed needle reflectance suggested a strong positive relationship between reflectance and tree health. For Landsat TM images, previous research had only observed a weak positive relationship between stand health and near infrared reflectance in these pine canopies. This suggests that for multispectral Landsat TM images, reflectance of near infrared light from pine canopies may be affected by other factors which may include the scattering of light within canopies. These results are seen as promising for the use of hyperspectral images to detect stand health, provided that pixel reflectance is not influenced by other scene components.

An improved seeded region growing algorithm

Recently Adams and Bischof (1994) proposed a novel region growing algorithm for segmenting intensity images. The inputs to the algorithm are the intensity image and a set of seeds - individual points or connected components - that identify the individual regions to be segmented. The algorithm grows these seed regions until all of the image pixels have been assimilated. Unfortunately the algorithm is inherently dependent on the order of pixel processing. This means, for example, that raster order processing and anti-raster order processing do not, in general, lead to the same tessellation. In this paper we propose an improved seeded region growing algorithm that retains the advantages of the Adams and Bischof algorithm fast execution, robust segmentation, and no tuning parameters - but is pixel order independent. (C) 1997 Elsevier Science B.V.

Centrifugal processing of cell debris and inclusion bodies from recombinant Escherichia coli

The settling characteristics of cell debris and inclusion bodies prior to, and following, fractionation in a disc-stack centrifuge were measured using Cumulative Sedimentation Analysis (CSA) and Centrifugal Disc photosedimentation (CDS). The impact of centrifuge feedrate and repeated homogenisation on both cell debris and inclusion body collection efficiency was investigated. Increasing the normalised centrifuge feedrate (Q/Sigma) from 1.32 x 10(-9) m s(-1) to 3.97 x 10(-9) m s(-1) leads to a 36% increase in inclusion body paste purity. Purity may also be improved by repeated homogenisation. Increasing the number of homogeniser passes results in smaller cell debris size whilst leaves inclusion body size unaltered. At a normalised centrifuge feedrate of 2.65 x 10(-9) m s(-1), increasing the number of homogeniser passes from two (2) to ten (10) improved overall inclusion body paste purity by 58%. Grade-efficiency curves for both the cell debris and inclusion bodies have also been generated in this study. The data are described using an equation developed by Mannweiler (1989) with parameters of k = 0.15-0.26 and n = 2.5-2.6 for inclusion bodies, and k = 0.12-0.14 and n = 2.0-2.2 for cell debris. This is the first accurate experimentally-determined grade efficiency curve for cell debris. Previous studies have simply estimated debris grade efficiency curves using an approximate debris size distribution and grade efficiency curves determined with 'ideal particles' (e.g. spherical PVA particles). The findings of this study may be used to simulate and optimise the centrifugal fractionation of inclusion bodies from cell debris.

Quantifying subclinical central nervous lesions in primary antiphospholipid syndrome: The role of magnetization transfer imaging

Purpose: To define the role of magnetization transfer imaging (MTI) in detecting subclinical central nervous system (CNS) lesions in primary antiphospholipid syndrome (PAPS). Materials and Methods: Ten non-CNS PAPS patients were compared to 10 CNS PAPS patients and 10 age- and sex-matched controls. All PAPS patients met Sapporo criteria. All Subjects underwent conventional MRI and complementary MTI analysis to compose histograms. CNS viability was determined according to the magnetization transfer ratio (MTR) by mean pixel intensity (MPI) and the mean peak height (MPH). Volumetric cerebral measurements were assessed by brain parenchyma factor (BPF) and total/cerebral volume. Results: MTR histograms analysis revealed that MPI was significantly different among groups (P < 0.0001). Non-CNS PAPS had a higher MPI than CNS PAPS, (30.5 +/- 1.01 vs. 25.1 +/- 3.17 percent unit (pu); P < 0.05) although lower than controls (30.5 +/- 1.01 vs. 31.20 < 0.50 pu; P < 0.05). MPH in non-CNS PAPS (5.57 +/- 0.20% (1/pu)} was similar to controls (5.63 +/- 0.20% (1/pu), P > 0.05) and higher than CNS PAPS (4.71 +/- 0.30% (1/pu), P < 0.05). A higher peak location (PL) was also observed in the CNS PAPS group in comparison with the other groups (P < 0.0001). In addition, a lower BPF was found in non-CNS PAPS compared to controls (0.80 +/- 0.03 vs. 0.84 +/- 0.02 units; P < 0.05) but similar to CNS PAPS (0.80 +/- 0.03 vs. 0.79 +/- 0.05 units; P > 0.05). Conclusion: Our findings suggest that non-CNS PAPS patients have subclinical cerebral damage. The long-term-clinical relevance of MTI analysis in these patients needs to be defined by prospective studies.

Extraction of electrical characteristics from pixels of multifrequency EIT images

Computer modelling has shown that electrical characteristics of individual pixels may be extracted from within multiple-frequency electrical impedance tomography (MFEIT) images formed using a reference data set obtained from a purely resistive, homogeneous medium. In some applications it is desirable to extract the electrical characteristics of individual pixels from images where a purely resistive, homogeneous reference data set is not available. One such application of the technique of MFEIT is to allow the acquisition of in vivo images using reference data sets obtained from a non-homogeneous medium with a reactive component. However, the reactive component of the reference data set introduces difficulties with the extraction of the true electrical characteristics from the image pixels. This study was a preliminary investigation of a technique to extract electrical parameters from multifrequency images when the reference data set has a reactive component. Unlike the situation in which a homogenous, resistive data set is available, it is not possible to obtain the impedance and phase information directly from the image pixel values of the MFEIT images data set, as the phase of the reactive reference is not known. The method reported here to extract the electrical characteristics (the Cole-Cole plot) initially assumes that this phase angle is zero. With this assumption, an impedance spectrum can be directly extracted from the image set. To obtain the true Cole-Cole plot a correction must be applied to account for the inherent rotation of the extracted impedance spectrum about the origin, which is a result of the assumption. This work shows that the angle of rotation associated with the reactive component of the reference data set may be determined using a priori knowledge of the distribution of frequencies of the Cole-Cole plot. Using this angle of rotation, the true Cole-Cole plot can be obtained from the impedance spectrum extracted from the MFEIT image data set. The method was investigated using simulated data, both with and without noise, and also for image data obtained in vitro. The in vitro studies involved 32 logarithmically spaced frequencies from 4 kHz up to 1 MHz and demonstrated that differences between the true characteristics and those of the impedance spectrum were reduced significantly after application of the correction technique. The differences between the extracted parameters and the true values prior to correction were in the range from 16% to 70%. Following application of the correction technique the differences were reduced to less than 5%. The parameters obtained from the Cole-Cole plot may be useful as a characterization of the nature and health of the imaged tissues.

A mathematical model for Escherichia coli debris size reduction during high pressure homogenisation based on grinding theory

Experimental data for E. coli debris size reduction during high-pressure homogenisation at 55 MPa are presented. A mathematical model based on grinding theory is developed to describe the data. The model is based on first-order breakage and compensation conditions. It does not require any assumption of a specified distribution for debris size and can be used given information on the initial size distribution of whole cells and the disruption efficiency during homogenisation. The number of homogeniser passes is incorporated into the model and used to describe the size reduction of non-induced stationary and induced E. coil cells during homogenisation. Regressing the results to the model equations gave an excellent fit to experimental data ( > 98.7% of variance explained for both fermentations), confirming the model's potential for predicting size reduction during high-pressure homogenisation. This study provides a means to optimise both homogenisation and disc-stack centrifugation conditions for recombinant product recovery. (C) 1997 Elsevier Science Ltd.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

Kallmann syndrome and mirror movements: White matter quantitative evaluation with magnetic resonance imaging

Kallmann syndrome (KS), characterized by the association of hypogonadotropic hypogonadism and anosmia, may present many other phenotypic abnormalities, including neurologic features as involuntary movements, called mirror movements (MM). MM etiology probably involves a complex mechanism comprising corticospinal tract abnormal development associated with deficient contralateral motor cortex inhibitory system. In this study, in order to address previous hypotheses concerning MM etiology, we identified and quantified white matter (WM) alterations in 21 KS patients, comparing subjects with and without MM and 16 control subjects, using magnetization transfer ratio (MTR) and T2 relaxometry (R2). Magnetization transfer and 12 double-echo images were acquired in a 1.5 T system. MTR and R2 were calculated pixel by pixel to initially create individual maps, and then, group average maps, co-registered with MNI305 stereotaxic coordinate system. After analysis of selected regions of interest, we demonstrated areas with higher 12 relaxation time and lower MTR values in KS patients, with and without MM, differently involving corticospinal tract projection, frontal lobes and corpus callosum. Higher MTR was observed only in pyramidal decussation when compared in both groups of patients with controls. In conclusion, we demonstrated that patients with KS have altered WM areas, presenting in a different manner in patients with and without MM. These data suggest axonal loss or disorganization involving abnormal pyramidal tracts and other associative/connective areas, relating to the presence or absence of MM. We also found a different pattern of alteration in pyramidal decussation, which can represent the primary area of neuronal disarrangement. (C) 2010 Elsevier B.V. All rights reserved.