Fully Constrained Least Squares Spectral Unmixing by Simplex Projection

We present a new algorithm for linear spectral mixture analysis, which is capable of supervised unmixing of hyperspectral data while respecting the constraints on the abundance coefficients. This simplex-projection unmixing algorithm is based upon the equivalence of the fully constrained least squares problem and the problem of projecting a point onto a simplex. We introduce several geometrical properties of high-dimensional simplices and combine them to yield a recursive algorithm for solving the simplex-projection problem. A concrete implementation of the algorithm for large data sets is provided, and the algorithm is benchmarked against well-known fully constrained least squares unmixing (FCLSU) techniques, on both artificial data sets and real hyperspectral data collected over the Cuprite mining region. Unlike previous algorithms for FCLSU, the presented algorithm possesses no optimization steps and is completely analytical, severely reducing the required processing power.

[1]  Chein-I Chang,et al.  Constrained band selection for hyperspectral imagery , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[2]  Yixiong Lin,et al.  Geometric linear algebra , 2005 .

[3]  Paul D. Gader,et al.  Sparsity Promoting Iterated Constrained Endmember Detection in Hyperspectral Imagery , 2007, IEEE Geoscience and Remote Sensing Letters.

[4]  C. Michelot.  A finite algorithm for finding the projection of a point onto the canonical simplex of ∝n , 1986 .

[5]  Chein-I Chang,et al.  Estimation of number of spectrally distinct signal sources in hyperspectral imagery , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[6]  Antonio J. Plaza,et al.  Survey of geometric and statistical unmixing algorithms for hyperspectral images , 2010, 2010 2nd Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing.

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

[8]  John F. Mustard,et al.  Spectral unmixing , 2002, IEEE Signal Process. Mag..

[9]  Sen Jia,et al.  Constrained Nonnegative Matrix Factorization for Hyperspectral Unmixing , 2009, IEEE Transactions on Geoscience and Remote Sensing.

[10]  Yoram Singer,et al.  Efficient projections onto the l1-ball for learning in high dimensions , 2008, ICML '08.

[11]  Stefan A. Robila,et al.  Considerations on Parallelizing Nonnegative Matrix Factorization for Hyperspectral Data Unmixing , 2009, IEEE Geoscience and Remote Sensing Letters.

[12]  Yosio Edemir Shimabukuro,et al.  The least-squares mixing models to generate fraction images derived from remote sensing multispectral data , 1991, IEEE Trans. Geosci. Remote. Sens..

[13]  Manuel Graña,et al.  Endmember Extraction Methods: A Short Review , 2008, KES.

[14]  Fred A. Kruse,et al.  Comparison of airborne hyperspectral data and EO-1 Hyperion for mineral mapping , 2003, IEEE Trans. Geosci. Remote. Sens..

[15]  Chein-I Chang,et al.  Estimation of subpixel target size for remotely sensed imagery , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[16]  David Gillis,et al.  Using endmembers as a coordinate system in hyperspectral imagery , 2002, SPIE Optics + Photonics.

[17]  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..

[18]  Jin Chen,et al.  Generalization of Subpixel Analysis for Hyperspectral Data With Flexibility in Spectral Similarity Measures , 2009, IEEE Transactions on Geoscience and Remote Sensing.

[19]  Antonio J. Plaza,et al.  A quantitative and comparative analysis of endmember extraction algorithms from hyperspectral data , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[20]  Chein-I Chang,et al.  Hyperspectral image classification and dimensionality reduction: an orthogonal subspace projection approach , 1994, IEEE Trans. Geosci. Remote. Sens..

[21]  D. Brie,et al.  Separation of Non-Negative Mixture of Non-Negative Sources Using a Bayesian Approach and MCMC Sampling , 2006, IEEE Transactions on Signal Processing.

[22]  Chein-I Chang,et al.  A New Growing Method for Simplex-Based Endmember Extraction Algorithm , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[23]  I-Hsiung Lin.  Geometric Linear Algebra volume 2 , 2008 .

[24]  J. Boardman,et al.  Geometric mixture analysis of imaging spectrometry data , 1994, Proceedings of IGARSS '94 - 1994 IEEE International Geoscience and Remote Sensing Symposium.

[25]  Jon Atli Benediktsson,et al.  On the decomposition of Mars hyperspectral data by ICA and Bayesian positive source separation , 2008, Neurocomputing.

[26]  C. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[27]  Alfred O. Hero,et al.  Joint Bayesian Endmember Extraction and Linear Unmixing for Hyperspectral Imagery , 2009, IEEE Transactions on Signal Processing.

[28]  Peter Bajorski,et al.  Simplex projection methods for selection of endmembers in hyperspectral imagery , 2004, IGARSS 2004. 2004 IEEE International Geoscience and Remote Sensing Symposium.

[29]  Saïd Moussaoui,et al.  Implementation Strategies for Hyperspectral Unmixing Using Bayesian Source Separation , 2010, IEEE Transactions on Geoscience and Remote Sensing.

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

[31]  Yoram Singer,et al.  Efficient Learning of Label Ranking by Soft Projections onto Polyhedra , 2006, J. Mach. Learn. Res..

[32]  Miguel Velez-Reyes,et al.  Solving adundance estimation in hyperspectral unmixing as a least distance problem , 2004, IGARSS 2004. 2004 IEEE International Geoscience and Remote Sensing Symposium.

[33]  Le Wang,et al.  Fully Constrained Linear Spectral Unmixing: Analytic Solution Using Fuzzy Sets , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[34]  Chein-I Chang,et al.  Multispectral and hyperspectral image analysis with convex cones , 1999, IEEE Trans. Geosci. Remote. Sens..

[35]  D. Stein,et al.  Application of the normal compositional model to the analysis of hyperspectral imagery , 2003, IEEE Workshop on Advances in Techniques for Analysis of Remotely Sensed Data, 2003.

[36]  Mark L. G. Althouse,et al.  Least squares subspace projection approach to mixed pixel classification for hyperspectral images , 1998, IEEE Trans. Geosci. Remote. Sens..

[37]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[38]  J. Settle,et al.  Linear mixing and the estimation of ground cover proportions , 1993 .

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