Blind Spectral Unmixing Based on Sparse Nonnegative Matrix Factorization

Nonnegative matrix factorization (NMF) is a widely used method for blind spectral unmixing (SU), which aims at obtaining the endmembers and corresponding fractional abundances, knowing only the collected mixing spectral data. It is noted that the abundance may be sparse (i.e., the endmembers may be with sparse distributions) and sparse NMF tends to lead to a unique result, so it is intuitive and meaningful to constrain NMF with sparseness for solving SU. However, due to the abundance sum-to-one constraint in SU, the traditional sparseness measured by L0/L1-norm is not an effective constraint any more. A novel measure (termed as S-measure) of sparseness using higher order norms of the signal vector is proposed in this paper. It features the physical significance. By using the S-measure constraint (SMC), a gradient-based sparse NMF algorithm (termed as NMF-SMC) is proposed for solving the SU problem, where the learning rate is adaptively selected, and the endmembers and abundances are simultaneously estimated. In the proposed NMF-SMC, there is no pure index assumption and no need to know the exact sparseness degree of the abundance in prior. Yet, it does not require the preprocessing of dimension reduction in which some useful information may be lost. Experiments based on synthetic mixtures and real-world images collected by AVIRIS and HYDICE sensors are performed to evaluate the validity of the proposed method.

Spectral unmixing based on nonnegative matrix factorization with local smoothness constraint

Spectral unmixing (SU) is an emerging problem in the remote sensing image processing. Since both the endmember signatures and their abundances have nonnegative values, it is a natural choice to employ the attractive nonnegative matrix factorization (NMF) methods to solve this problem. Motivated by that the abundances are sparse, the NMF with local smoothness constraint (NMF-LSC) is proposed in this paper. In the proposed method, the smoothness constraint is utilized to impose the sparseness, instead of the traditional L1-norm which is restricted by the underlying column-sum-to-one requirement of the to the abundance matrix. Simulations show the advantages of our algorithm over the compared methods.

Fluctuating dimension in a discrete model for quantum gravity based on the spectral principle

The spectral principle of Connes and Chamseddine is used as a starting point to define a discrete model for Euclidean quantum gravity. Instead of summing over ordinary geometries, we consider the sum over generalized geometries where topology, metric, and dimension can fluctuate. The model describes the geometry of spaces with a countable number n of points, and is related to the Gaussian unitary ensemble of Hermitian matrices. We show that this simple model has two phases. The expectation value , the average number of points in the Universe, is finite in one phase and diverges in the other. We compute the critical point as well as the critical exponent of . Moreover, the space-time dimension delta is a dynamical observable in our model, and plays the role of an order parameter. The computation of is discussed and an upper bound is found, < 2.

SZEGO and PARA-ORTHOGONAL POLYNOMIALS on THE REAL LINE: ZEROS and CANONICAL SPECTRAL TRANSFORMATIONS

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Szego{double acute} and para-orthogonal polynomials on the real line: Zeros and canonical spectral transformations

We study polynomials which satisfy the same recurrence relation as the Szego{double acute} polynomials, however, with the restriction that the (reflection) coefficients in the recurrence are larger than one in modulus. Para-orthogonal polynomials that follow from these Szego{double acute} polynomials are also considered. With positive values for the reflection coefficients, zeros of the Szego{double acute} polynomials, para-orthogonal polynomials and associated quadrature rules are also studied. Finally, again with positive values for the reflection coefficients, interlacing properties of the Szego{double acute} polynomials and polynomials arising from canonical spectral transformations are obtained. © 2012 American Mathematical Society.

Thermal 2-loop master spectral function at finite momentum

When considering NLO corrections to thermal particle production in the “relativistic” regime, in which the invariant mass squared of the produced particle is K2 ~ (πT)2, then the production rate can be expressed as a sum of a few universal “master” spectral functions. Taking the most complicated 2-loop master as an example, a general strategy for obtaining a convergent 2-dimensional integral representation is suggested. The analysis applies both to bosonic and fermionic statistics, and shows that for this master the non-relativistic approximation is only accurate for K2 ~(8πT)2, whereas the zero-momentum approximation works surprisingly well. Once the simpler masters have been similarly resolved, NLO results for quantities such as the right-handed neutrino production rate from a Standard Model plasma or the dilepton production rate from a QCD plasma can be assembled for K2 ~ (πT)2.

Classification of Satellite Images by means of Fuzzy Rules generated by a Genetic Algorithm

The data acquired by Remote Sensing systems allow obtaining thematic maps of the earth's surface, by means of the registered image classification. This implies the identification and categorization of all pixels into land cover classes. Traditionally, methods based on statistical parameters have been widely used, although they show some disadvantages. Nevertheless, some authors indicate that those methods based on artificial intelligence, may be a good alternative. Thus, fuzzy classifiers, which are based on Fuzzy Logic, include additional information in the classification process through based-rule systems. In this work, we propose the use of a genetic algorithm (GA) to select the optimal and minimum set of fuzzy rules to classify remotely sensed images. Input information of GA has been obtained through the training space determined by two uncorrelated spectral bands (2D scatter diagrams), which has been irregularly divided by five linguistic terms defined in each band. The proposed methodology has been applied to Landsat-TM images and it has showed that this set of rules provides a higher accuracy level in the classification process

Short-term spectral analysis of ODP Hole 114-704B

Short-term spectral analysis was carried out on geochemical logging data from ODP Site 704. The FFT was used to compute the amplitude spectra of short-term overlapping segments to produce depth-period-amplitude spectrograms of the logging data. The spectrograms provided a means of evaluating the significance of the observed periodic components. The periodic components that were consistently present and prominent across a given record interval were considered to be significant. Changes in the spectrogram characteristics seem to reflect changes in either lithology, sedimentation rates, or hiatuses and may therefore provide useful information to aid in stratigraphic and paleoenvironmental studies. The dominant periodicity during the late Pleistocene and Brunhes Chron (0.97 to 0.47 Ma) was determined to be > 100,000 yr whereas the upper Matuyama Chron was dominated by the 41,000-yr periodicity. These periodicities suggest that the sedimentation patterns within the upper Matuyama Chron (0.98-1.78 Ma) were influenced by the Milankovitch obliquity cycle and those within the latest Matuyama-Brunhes Chron (<0.98 Ma) by the eccentricity cycle. The Brunhes/Matuyama boundary therefore represents a major discontinuity. Periodicities observed within the lower Matuyama and the upper Gauss Chron did not correlate with any of the periodicities within the Milankovitch frequency bands.

Downwelling spectral irradiance as measured in different depths, remote sensing reflectance and phytoplankton pigment concentration during POLARSTERN cruises ANT-XXIV/1, ANT-XXIV/4, ANT-XXVI/4 and Maria S. Merian cruise MSM18/3

The composition and abundance of algal pigments provide information on phytoplankton community characteristics such as photoacclimation, overall biomass and taxonomic composition. In particular, pigments play a major role in photoprotection and in the light-driven part of photosynthesis. Most phytoplankton pigments can be measured by high-performance liquid chromatography (HPLC) techniques applied to filtered water samples. This method, as well as other laboratory analyses, is time consuming and therefore limits the number of samples that can be processed in a given time. In order to receive information on phytoplankton pigment composition with a higher temporal and spatial resolution, we have developed a method to assess pigment concentrations from continuous optical measurements. The method applies an empirical orthogonal function (EOF) analysis to remote-sensing reflectance data derived from ship-based hyperspectral underwater radiometry and from multispectral satellite data (using the Medium Resolution Imaging Spectrometer - MERIS - Polymer product developed by Steinmetz et al., 2011, doi:10.1364/OE.19.009783) measured in the Atlantic Ocean. Subsequently we developed multiple linear regression models with measured (collocated) pigment concentrations as the response variable and EOF loadings as predictor variables. The model results show that surface concentrations of a suite of pigments and pigment groups can be well predicted from the ship-based reflectance measurements, even when only a multispectral resolution is chosen (i.e., eight bands, similar to those used by MERIS). Based on the MERIS reflectance data, concentrations of total and monovinyl chlorophyll a and the groups of photoprotective and photosynthetic carotenoids can be predicted with high quality. As a demonstration of the utility of the approach, the fitted model based on satellite reflectance data as input was applied to 1 month of MERIS Polymer data to predict the concentration of those pigment groups for the whole eastern tropical Atlantic area. Bootstrapping explorations of cross-validation error indicate that the method can produce reliable predictions with relatively small data sets (e.g., < 50 collocated values of reflectance and pigment concentration). The method allows for the derivation of time series from continuous reflectance data of various pigment groups at various regions, which can be used to study variability and change of phytoplankton composition and photophysiology.