On derivatives and norms of generalized matrix functions and respective symmetric powers

In recent papers, the authors obtained formulas for directional derivatives of all orders, of the immanant and of the m-th xi-symmetric tensor power of an operator and a matrix, when xi is a character of the full symmetric group. The operator norm of these derivatives was also calculated. In this paper, similar results are established for generalized matrix functions and for every symmetric tensor power.

Sparse matrix multiplication on a reconfigurable many-core architecture

Sparse matrix-vector multiplication (SMVM) is a fundamental operation in many scientific and engineering applications. In many cases sparse matrices have thousands of rows and columns where most of the entries are zero, while non-zero data is spread over the matrix. This sparsity of data locality reduces the effectiveness of data cache in general-purpose processors quite reducing their performance efficiency when compared to what is achieved with dense matrix multiplication. In this paper, we propose a parallel processing solution for SMVM in a many-core architecture. The architecture is tested with known benchmarks using a ZYNQ-7020 FPGA. The architecture is scalable in the number of core elements and limited only by the available memory bandwidth. It achieves performance efficiencies up to almost 70% and better performances than previous FPGA designs.

Otimização da Amostragem e Estudo de Viabilidade da Monitorização de Sondas na ETAR de Rabada

Mestrado em Engenharia Química - Tecnologias de Protecção Ambiental

Unmixing hyperspectral data: independent and dependent component analysis

The development of high spatial resolution airborne and spaceborne sensors has improved the capability of ground-based data collection in the fields of agriculture, geography, geology, mineral identification, detection [2, 3], and classification [4–8]. The signal read by the sensor from a given spatial element of resolution and at a given spectral band is a mixing of components originated by the constituent substances, termed endmembers, located at that element of resolution. This chapter addresses hyperspectral unmixing, which is the decomposition of the pixel spectra into a collection of constituent spectra, or spectral signatures, and their corresponding fractional abundances indicating the proportion of each endmember present in the pixel [9, 10]. Depending on the mixing scales at each pixel, the observed mixture is either linear or nonlinear [11, 12]. The linear mixing model holds when the mixing scale is macroscopic [13]. The nonlinear model holds when the mixing scale is microscopic (i.e., intimate mixtures) [14, 15]. The linear model assumes negligible interaction among distinct endmembers [16, 17]. The nonlinear model assumes that incident solar radiation is scattered by the scene through multiple bounces involving several endmembers [18]. Under the linear mixing model and assuming that the number of endmembers and their spectral signatures are known, hyperspectral unmixing is a linear problem, which can be addressed, for example, under the maximum likelihood setup [19], the constrained least-squares approach [20], the spectral signature matching [21], the spectral angle mapper [22], and the subspace projection methods [20, 23, 24]. Orthogonal subspace projection [23] reduces the data dimensionality, suppresses undesired spectral signatures, and detects the presence of a spectral signature of interest. The basic concept is to project each pixel onto a subspace that is orthogonal to the undesired signatures. As shown in Settle [19], the orthogonal subspace projection technique is equivalent to the maximum likelihood estimator. This projection technique was extended by three unconstrained least-squares approaches [24] (signature space orthogonal projection, oblique subspace projection, target signature space orthogonal projection). Other works using maximum a posteriori probability (MAP) framework [25] and projection pursuit [26, 27] have also been applied to hyperspectral data. In most cases the number of endmembers and their signatures are not known. Independent component analysis (ICA) is an unsupervised source separation process that has been applied with success to blind source separation, to feature extraction, and to unsupervised recognition [28, 29]. ICA consists in finding a linear decomposition of observed data yielding statistically independent components. Given that hyperspectral data are, in given circumstances, linear mixtures, ICA comes to mind as a possible tool to unmix this class of data. In fact, the application of ICA to hyperspectral data has been proposed in reference 30, where endmember signatures are treated as sources and the mixing matrix is composed by the abundance fractions, and in references 9, 25, and 31–38, where sources are the abundance fractions of each endmember. In the first approach, we face two problems: (1) The number of samples are limited to the number of channels and (2) the process of pixel selection, playing the role of mixed sources, is not straightforward. In the second approach, ICA is based on the assumption of mutually independent sources, which is not the case of hyperspectral data, since the sum of the abundance fractions is constant, implying dependence among abundances. This dependence compromises ICA applicability to hyperspectral images. In addition, hyperspectral data are immersed in noise, which degrades the ICA performance. IFA [39] was introduced as a method for recovering independent hidden sources from their observed noisy mixtures. IFA implements two steps. First, source densities and noise covariance are estimated from the observed data by maximum likelihood. Second, sources are reconstructed by an optimal nonlinear estimator. Although IFA is a well-suited technique to unmix independent sources under noisy observations, the dependence among abundance fractions in hyperspectral imagery compromises, as in the ICA case, the IFA performance. Considering the linear mixing model, hyperspectral observations are in a simplex whose vertices correspond to the endmembers. Several approaches [40–43] have exploited this geometric feature of hyperspectral mixtures [42]. Minimum volume transform (MVT) algorithm [43] determines the simplex of minimum volume containing the data. The MVT-type approaches are complex from the computational point of view. Usually, these algorithms first find the convex hull defined by the observed data and then fit a minimum volume simplex to it. Aiming at a lower computational complexity, some algorithms such as the vertex component analysis (VCA) [44], the pixel purity index (PPI) [42], and the N-FINDR [45] still find the minimum volume simplex containing the data cloud, but they assume the presence in the data of at least one pure pixel of each endmember. This is a strong requisite that may not hold in some data sets. In any case, these algorithms find the set of most pure pixels in the data. Hyperspectral sensors collects spatial images over many narrow contiguous bands, yielding large amounts of data. For this reason, very often, the processing of hyperspectral data, included unmixing, is preceded by a dimensionality reduction step to reduce computational complexity and to improve the signal-to-noise ratio (SNR). Principal component analysis (PCA) [46], maximum noise fraction (MNF) [47], and singular value decomposition (SVD) [48] are three well-known projection techniques widely used in remote sensing in general and in unmixing in particular. The newly introduced method [49] exploits the structure of hyperspectral mixtures, namely the fact that spectral vectors are nonnegative. The computational complexity associated with these techniques is an obstacle to real-time implementations. To overcome this problem, band selection [50] and non-statistical [51] algorithms have been introduced. This chapter addresses hyperspectral data source dependence and its impact on ICA and IFA performances. The study consider simulated and real data and is based on mutual information minimization. Hyperspectral observations are described by a generative model. This model takes into account the degradation mechanisms normally found in hyperspectral applications—namely, signature variability [52–54], abundance constraints, topography modulation, and system noise. The computation of mutual information is based on fitting mixtures of Gaussians (MOG) to data. The MOG parameters (number of components, means, covariances, and weights) are inferred using the minimum description length (MDL) based algorithm [55]. We study the behavior of the mutual information as a function of the unmixing matrix. The conclusion is that the unmixing matrix minimizing the mutual information might be very far from the true one. Nevertheless, some abundance fractions might be well separated, mainly in the presence of strong signature variability, a large number of endmembers, and high SNR. We end this chapter by sketching a new methodology to blindly unmix hyperspectral data, where abundance fractions are modeled as a mixture of Dirichlet sources. This model enforces positivity and constant sum sources (full additivity) constraints. The mixing matrix is inferred by an expectation-maximization (EM)-type algorithm. This approach is in the vein of references 39 and 56, replacing independent sources represented by MOG with mixture of Dirichlet sources. Compared with the geometric-based approaches, the advantage of this model is that there is no need to have pure pixels in the observations. The chapter is organized as follows. Section 6.2 presents a spectral radiance model and formulates the spectral unmixing as a linear problem accounting for abundance constraints, signature variability, topography modulation, and system noise. Section 6.3 presents a brief resume of ICA and IFA algorithms. Section 6.4 illustrates the performance of IFA and of some well-known ICA algorithms with experimental data. Section 6.5 studies the ICA and IFA limitations in unmixing hyperspectral data. Section 6.6 presents results of ICA based on real data. Section 6.7 describes the new blind unmixing scheme and some illustrative examples. Section 6.8 concludes with some remarks.

Vertex component analysis: a geometric-based approach to unmix hyperspectral data

Hyperspectral remote sensing exploits the electromagnetic scattering patterns of the different materials at specific wavelengths [2, 3]. Hyperspectral sensors have been developed to sample the scattered portion of the electromagnetic spectrum extending from the visible region through the near-infrared and mid-infrared, in hundreds of narrow contiguous bands [4, 5]. The number and variety of potential civilian and military applications of hyperspectral remote sensing is enormous [6, 7]. Very often, the resolution cell corresponding to a single pixel in an image contains several substances (endmembers) [4]. In this situation, the scattered energy is a mixing of the endmember spectra. A challenging task underlying many hyperspectral imagery applications is then decomposing a mixed pixel into a collection of reflectance spectra, called endmember signatures, and the corresponding abundance fractions [8–10]. Depending on the mixing scales at each pixel, the observed mixture is either linear or nonlinear [11, 12]. Linear mixing model holds approximately when the mixing scale is macroscopic [13] and there is negligible interaction among distinct endmembers [3, 14]. If, however, the mixing scale is microscopic (or intimate mixtures) [15, 16] and the incident solar radiation is scattered by the scene through multiple bounces involving several endmembers [17], the linear model is no longer accurate. Linear spectral unmixing has been intensively researched in the last years [9, 10, 12, 18–21]. It considers that a mixed pixel is a linear combination of endmember signatures weighted by the correspondent abundance fractions. Under this model, and assuming that the number of substances and their reflectance spectra are known, hyperspectral unmixing is a linear problem for which many solutions have been proposed (e.g., maximum likelihood estimation [8], spectral signature matching [22], spectral angle mapper [23], subspace projection methods [24,25], and constrained least squares [26]). In most cases, the number of substances and their reflectances are not known and, then, hyperspectral unmixing falls into the class of blind source separation problems [27]. Independent component analysis (ICA) has recently been proposed as a tool to blindly unmix hyperspectral data [28–31]. ICA is based on the assumption of mutually independent sources (abundance fractions), which is not the case of hyperspectral data, since the sum of abundance fractions is constant, implying statistical dependence among them. This dependence compromises ICA applicability to hyperspectral images as shown in Refs. [21, 32]. In fact, ICA finds the endmember signatures by multiplying the spectral vectors with an unmixing matrix, which minimizes the mutual information among sources. If sources are independent, ICA provides the correct unmixing, since the minimum of the mutual information is obtained only when sources are independent. This is no longer true for dependent abundance fractions. Nevertheless, some endmembers may be approximately unmixed. These aspects are addressed in Ref. [33]. Under the linear mixing model, the observations from a scene are in a simplex whose vertices correspond to the endmembers. Several approaches [34–36] have exploited this geometric feature of hyperspectral mixtures [35]. Minimum volume transform (MVT) algorithm [36] determines the simplex of minimum volume containing the data. The method presented in Ref. [37] is also of MVT type but, by introducing the notion of bundles, it takes into account the endmember variability usually present in hyperspectral mixtures. The MVT type approaches are complex from the computational point of view. Usually, these algorithms find in the first place the convex hull defined by the observed data and then fit a minimum volume simplex to it. For example, the gift wrapping algorithm [38] computes the convex hull of n data points in a d-dimensional space with a computational complexity of O(nbd=2cþ1), where bxc is the highest integer lower or equal than x and n is the number of samples. The complexity of the method presented in Ref. [37] is even higher, since the temperature of the simulated annealing algorithm used shall follow a log( ) law [39] to assure convergence (in probability) to the desired solution. Aiming at a lower computational complexity, some algorithms such as the pixel purity index (PPI) [35] and the N-FINDR [40] still find the minimum volume simplex containing the data cloud, but they assume the presence of at least one pure pixel of each endmember in the data. This is a strong requisite that may not hold in some data sets. In any case, these algorithms find the set of most pure pixels in the data. PPI algorithm uses the minimum noise fraction (MNF) [41] as a preprocessing step to reduce dimensionality and to improve the signal-to-noise ratio (SNR). The algorithm then projects every spectral vector onto skewers (large number of random vectors) [35, 42,43]. The points corresponding to extremes, for each skewer direction, are stored. A cumulative account records the number of times each pixel (i.e., a given spectral vector) is found to be an extreme. The pixels with the highest scores are the purest ones. N-FINDR algorithm [40] is based on the fact that in p spectral dimensions, the p-volume defined by a simplex formed by the purest pixels is larger than any other volume defined by any other combination of pixels. This algorithm finds the set of pixels defining the largest volume by inflating a simplex inside the data. ORA SIS [44, 45] is a hyperspectral framework developed by the U.S. Naval Research Laboratory consisting of several algorithms organized in six modules: exemplar selector, adaptative learner, demixer, knowledge base or spectral library, and spatial postrocessor. The first step consists in flat-fielding the spectra. Next, the exemplar selection module is used to select spectral vectors that best represent the smaller convex cone containing the data. The other pixels are rejected when the spectral angle distance (SAD) is less than a given thresh old. The procedure finds the basis for a subspace of a lower dimension using a modified Gram–Schmidt orthogonalizati on. The selected vectors are then projected onto this subspace and a simplex is found by an MV T pro cess. ORA SIS is oriented to real-time target detection from uncrewed air vehicles using hyperspectral data [46]. In this chapter we develop a new algorithm to unmix linear mixtures of endmember spectra. First, the algorithm determines the number of endmembers and the signal subspace using a newly developed concept [47, 48]. Second, the algorithm extracts the most pure pixels present in the data. Unlike other methods, this algorithm is completely automatic and unsupervised. To estimate the number of endmembers and the signal subspace in hyperspectral linear mixtures, the proposed scheme begins by estimating sign al and noise correlation matrices. The latter is based on multiple regression theory. The signal subspace is then identified by selectin g the set of signal eigenvalue s that best represents the data, in the least-square sense [48,49 ], we note, however, that VCA works with projected and with unprojected data. The extraction of the end members exploits two facts: (1) the endmembers are the vertices of a simplex and (2) the affine transformation of a simplex is also a simplex. As PPI and N-FIND R algorithms, VCA also assumes the presence of pure pixels in the data. The algorithm iteratively projects data on to a direction orthogonal to the subspace spanned by the endmembers already determined. The new end member signature corresponds to the extreme of the projection. The algorithm iterates until all end members are exhausted. VCA performs much better than PPI and better than or comparable to N-FI NDR; yet it has a computational complexity between on e and two orders of magnitude lower than N-FINDR. The chapter is structure d as follows. Section 19.2 describes the fundamentals of the proposed method. Section 19.3 and Section 19.4 evaluate the proposed algorithm using simulated and real data, respectively. Section 19.5 presents some concluding remarks.

Comparison of Disposable Pipette Extraction and Dispersive Solid-Phase Extraction in the QuEChERS Method for Analysis of Pesticides in Strawberries

In this study, we sought to assess the applicability of GC–MS/MS for the identification and quantification of 36 pesticides in strawberry from integrated pest management (IPM) and organic farming (OF). Citrate versions of QuEChERS (quick, easy, cheap, effective, rugged and safe) using dispersive solid-phase extraction (d-SPE) and disposable pipette extraction (DPX) for cleanup were compared for pesticide extraction. For cleanup, a combination of MgSO4, primary secondary amine and C18 was used for both the versions. Significant differences were observed in recovery results between the two sample preparation versions (DPX and d-SPE). Overall, 86% of the pesticides achieved recoveries (three spiking levels 10, 50 and 200 µg/kg) in the range of 70–120%, with <13% RSD. The matrix effects were also evaluated in both the versions and in strawberries from different crop types. Although not evidencing significant differences between the two methodologies were observed, however, the DPX cleanup proved to be a faster technique and easy to execute. The results indicate that QuEChERS with d-SPE and DPX and GC–MS/MS analysis achieved reliable quantification and identification of 36 pesticide residues in strawberries from OF and IPM.

Paracoccidioides brasiliensis: a mycologic and immunochemical study of a sample isolated from an armadillo (Dasipus novencinctus)

A sample of P. brasiliensis isolated from the spleen and the liver of an armadillo (Dasipus novencinctus) has been analysed under a mycological and immunochemical viewpoint. The armadillo was captured in an area of Tucuruí (State of Pará, Brazil), the animal being already established as an enzootic reservoir of P. brasiliensis at that region of the country. This sample maintained in the fungal collection of the Tropical Medicine Institute of São Paulo (Brazil) numbered 135, has got all the characteristics of P. brasiliensis, with a strong antigenic power and low virulence for guinea-pigs and Wistar rats. The specific exoantigen of P. brasiliensis - the glycoprotein with a molecular weight of 43 kDa - was easily demonstrated with double immunodiffusion, immunoelectrophoresis, SDS-PAGE and immunobloting techniques.

New analytical methodologies for doping control – detection of anabolic androgenic steroids in human urine

Dissertation submitted to Faculdade de Ciências e Tecnologia - Universidade Nova de Lisboa in fulfilment of the requirements for the degree of Doctor of Philosophy (Biochemistry - Biotechnology)

DETECTION OF HTLV-I ANTIBODIES AND DNA IN BLOOD SAMPLE OF A PATIENT WITH MYELOPATHY IN NIGERIA

We describe a case of human T-lymphotropic virus type I associated myelopathy in a 50-year old woman in Nigeria. The patient presented with progressive loss of tone to the two lower limbs and later inability to walk. The HTLV-I antibody presence in the plasma collected from the patient was repeatedly detected by enzyme immunoassays (Abbott HTLV-I EIA and Coulter SELECT-HTLV I/II) and confirmed by Western blot technique. In addition, HTLV-I DNA was amplified from the genomic DNA isolated from the peripheral blood mononuclear cells of the patient by the polymerase chain reaction technique. This finding is significant being the first report of association of HTLV-I with myelopathy in Nigeria.

PREVALENCE OF ANTI- Toxocara ANTIBODIES IN A RANDOM SAMPLE OF INPATIENTS AT A CHILDREN'S HOSPITAL IN VITÓRIA, ESPÍRITO SANTO, BRAZIL

In the streets of Vitória, in the State of Espírito Santo, Brazil, are large number of stray dogs, many of which are infected with Toxocara canis, suggesting a high risk for human infection. In order to investigate the prevalence of Toxocara infection in children in Espírito Santo we studied the prevalence of anti-Toxocara antibodies in 100 random inpatients over one year of age, at the Children's Hospital N.S. da Glória, the reference children's hospital for the State.All the sera were collected during the period between October 1996 and January 1997. The mean age was 6.6±4.1 yrs. (1 to 14 yrs., median 6yrs.) and there were patients from all of the different wards of the hospital. Sixty-eigth patients came from the metropolitan area of Vitória and the other 32 from 17 other municipalities. The anti-Toxocara antibodies were investigated by ELISA-IgG using a secretory-excretory antigen obtained from second stage larvae. All sera were adsorbed with Ascaris suum antigen before the test. Thirty-nine sera (39%) were positive, predominantly from boys, but the gender difference was not statistically significant (boys:25/56 or 44.6%; girls:14/44 or 31.8%; p=0.311). The prevalence of positive sera was higher, but not statistically significant, in children from the urban periphery of metropolitan Vitória (formed by the cities of Vitória, Cariacica, Vila Velha, Serra and Viana) than in children from 17 other municipalities (44.1% and 28.1% respectively, p=0.190). Although the samples studied do not represent all children living in the State of Espírito Santo, since the Children's Hospital N.S. da Glória admits only patients from the state health system, it is probable that these results indicate a high frequency of Toxocara infection in children living in Espírito Santo. Further studies of population samples are necessary to ascertain the prevalence of Toxocara infection in our country.

The effect of method, standard and sample components on the total antioxidant capacity of commercial waters assessed by optical conventional assays

The total antioxidant capacity (TAC) of 28 flavoured water samples was assessed by ferric reducing antioxidant potential (FRAP), oxygen radical absorbance capacity (ORAC), trolox equivalent antioxidant capacity (TEAC) and total reactive antioxidant potential (TRAP) methods. It was observed that flavoured waters had higher antioxidant activity than the corresponding natural ones. The observed differences were attributed to flavours, juice and vitamins. Generally, higher TAC contents were obtained on lemon waters and lower values on guava and raspberry flavoured waters. Lower and higher TACs were obtained by TRAP and ORAC method, respectively. Statistical analysis suggested that vitamins and flavours increased the antioxidant content of the commercial waters.

A novel antibody-like material for breast cancer antigen CA15-3, used to track breast cancer by potentiometric transduction

This work presents the development of a low cost sensor device for the diagnosis of breast cancer in point-of-care, made with new synthetic biomimetic materials inside plasticized poly(vinyl chloride), PVC, membranes, for subsequent potentiometric detection. This concept was applied to target a conventional biomarker in breast cancer: Breast Cancer Antigen (CA15-3). The new biomimetic material was obtained by molecularly-imprinted technology. In this, a plastic antibody was obtained by polymerizing around the biomarker that acted as an obstacle to the growth of the polymeric matrix. The imprinted polymer was specifically synthetized by electropolymerization on an FTO conductive glass, by using cyclic voltammetry, including 40 cycles within -0.2 and 1.0 V. The reaction used for the polymerization included monomer (pyrrol, 5.0×10-3 mol/L) and protein (CA15-3, 100U/mL), all prepared in phosphate buffer saline (PBS), with a pH of 7.2 and 1% of ethylene glycol. The biomarker was removed from the imprinted sites by proteolytic action of proteinase K. The biomimetic material was employed in the construction of potentiometric sensors and tested with regard to its affinity and selectivity for binding CA15-3, by checking the analytical performance of the obtained electrodes. For this purpose, the biomimetic material was dispersed in plasticized PVC membranes, including or not a lipophilic ionic additive, and applied on a solid conductive support of graphite. The analytical behaviour was evaluated in buffer and in synthetic serum, with regard to linear range, limit of detection, repeatability, and reproducibility. This antibody-like material was tested in synthetic serum, and good results were obtained. The best devices were able to detect 5 times less CA15-3 than that required in clinical use. Selectivity assays were also performed, showing that the various serum components did not interfere with this biomarker. Overall, the potentiometric-based methods showed several advantages compared to other methods reported in the literature. The analytical process was simple, providing fast responses for a reduced amount of analyte, with low cost and feasible miniaturization. It also allowed the detection of a wide range of concentrations, diminishing the required efforts in previous sample pre-treating stages.

The effect of water inorganic matrix in ibuprofen adsorption onto activated carbon for water and wastewater treatment

Dissertação para obtenção do Grau de Mestre em Engenharia Química e Bioquímica