Abstract:An increasing number of astronomical instruments (on Earth and space-based) provide hyperspectral images, that is three-dimensional data cubes with two spatial dimensions and one spectral dimension. The intrinsic limitation in spatial resolution of these instruments implies that the spectra associated with pixels of such images are most often mixtures of the spectra of the “pure” components that exist in the considered region. In order to estimate the spectra and spatial abundances of these pure components, we here propose an original blind signal separation (BSS), that is to say an unsupervised unmixing method. Our approach is based on extensions and combinations of linear BSS methods that belong to two major classes of methods, namely nonnegative matrix factorization (NMF) and sparse component analysis (SCA). The former performs the decomposition of hyperspectral images, as a set of pure spectra and abundance maps, by using nonnegativity constraints, but the estimated solution is not unique: It highly depends on the initialization of the algorithm. The considered SCA methods are based on the assumption of the existence of points or tiny spatial zones where only one source is active (i.e., one pure component is present). These points or zones are then used to estimate the mixture and perform the decomposition. In real conditions, the assumption of perfect single-source points or zones is not always realistic. In such conditions, SCA yields approximate versions of the unknown sources and mixing coefficients. We propose to use part of these preliminary estimates from the SCA to initialize several runs of the NMF in order to refine these estimates and further constrain the convergence of the NMF algorithm. The proposed methods also estimate the number of pure components involved in the data and they provide error bars associated with the obtained solution. Detailed tests with synthetic data show that the decomposition achieved with such hybrid methods is nearly unique and provides good performance, illustrating the potential of applications to real data.

Decomposition Of Mixed Pixels Using Bayesian Self-Organizing Map (Bsom) Neural Networks

Decomposition of Mixed Pixels Based on Bayesian Self-Organizing Map and Gaussian Mixture Model.

Decomposition of Mixed Pixels in MODIS Data Using Bernstein Basis Functions

An approach based on self-organizing map and fuzzy membership for decomposition of mixed pixels in hyperspectral imagery

Blind Spectral Unmixing Based on Sparse Component Analysis for Hyperspectral Remote Sensing Imagery

A New Scheme for Decomposition of Mixed Pixels Based on Modified Nonnegative Matrix Factorization and Genetic-Algorithm

Subpixel Mapping Method of Hyperspectral Images Based on Modified Binary Quantum Particle Swarm Optimization

SUBPIXEL MAPPING OF HYPERSPECTRAL IMAGE BASED ON LINEAR SUBPIXEL FEATURE DETECTION AND OBJECT OPTIMIZATION

Unsupervised Hyperspectral Unmixing Based on Constrained Nonnegative Matrix Factorization and Particle Swarm Optimization

Hopfield Neural Network Approach for Supervised Nonlinear Spectral Unmixing

Unsupervised Nonlinear Hyperspectral Unmixing Based on the Generalized Bilinear Model

A Fast Fully Constrained Geometric Unmixing of Hyperspectral Images

New Scheme for Decomposition of Mixed Pixels of Remote Sensing Images Based on Nonnegative Matrix Factorization

Mixed Pixel Decomposition Based on Extended Fuzzy Clustering for Single Spectral Value Remote Sensing Images

A new scheme for blind decomposition of mixed pixels based on non-negative matrix factorization

Unmixing methods based on nonnegativity and weakly mixed pixels for astronomical hyperspectral datasets

A Kernel Spatial Complexity-Based Nonlinear Unmixing Method of Hyperspectral Imagery.

Nonlinear Spectral Mixture Analysis by Determining Per-Pixel Endmember Sets

Hyperspectral Unmixing Network Accounting for Spectral Variability Based on a Modified Scaled and a Perturbed Linear Mixing Model

Spatial-Spectral Nonlinear Hyperspectral Unmixing under Complex Noise

A maximum entropy approach to unsupervised mixed-pixel decomposition