Linear unmixing of hyperspectral images using a scaled gradient method

This paper addresses the problem of linear unmixing for hyperspectral imagery. This problem can be formulated as a linear regression problem whose regression coefficients (abundances) satisfy sum-toone and positivity constraints. Two scaled gradient iterative methods are proposed for estimating the abundances of the linear mixing model. The first method is obtained by including a normalization step in the scaled gradient method. The second method inspired by the fully constrained least squares algorithm includes the sum-to-one constraint in the observation model with an appropriate weighting parameter. Simulations on synthetic data illustrate the performance of these algorithms.

[1]  Chein-I Chang,et al.  Fully constrained least squares linear spectral mixture analysis method for material quantification in hyperspectral imagery , 2001, IEEE Trans. Geosci. Remote. Sens..

[2]  David A. Landgrebe,et al.  Signal Theory Methods in Multispectral Remote Sensing , 2003 .

[3]  H. Sebastian Seung,et al.  Learning the parts of objects by non-negative matrix factorization , 1999, Nature.

[4]  Chih-Jen Lin,et al.  Projected Gradient Methods for Nonnegative Matrix Factorization , 2007, Neural Computation.

[5]  Mario Winter,et al.  N-FINDR: an algorithm for fast autonomous spectral end-member determination in hyperspectral data , 1999, Optics & Photonics.

[6]  Chein-I Chang,et al.  Semi-Supervised Linear Spectral Unmixing Using a Hierarchical Bayesian Model for Hyperspectral Imagery , 2008, IEEE Transactions on Signal Processing.

[7]  Jing Wang,et al.  Applications of Independent Component Analysis in Endmember Extraction and Abundance Quantification for Hyperspectral Imagery , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[8]  M. Daube-Witherspoon,et al.  An Iterative Image Space Reconstruction Algorthm Suitable for Volume ECT , 1986, IEEE Transactions on Medical Imaging.

[9]  H. Lantéri,et al.  Penalized maximum likelihood image restoration with positivity constraints:multiplicative algorithms , 2002 .

[10]  Jorge Igual,et al.  Nonnegative Matrix Factorization of Laboratory Astrophysical Ice Mixtures , 2008, IEEE Journal of Selected Topics in Signal Processing.

[11]  José M. Bioucas-Dias,et al.  Vertex component analysis: a fast algorithm to unmix hyperspectral data , 2005, IEEE Transactions on Geoscience and Remote Sensing.

[12]  Christian P. Robert,et al.  Monte Carlo Statistical Methods (Springer Texts in Statistics) , 2005 .