Blind and Fully Constrained Unmixing of Hyperspectral Images

This paper addresses the problem of blind and fully constrained unmixing of hyperspectral images. Unmixing is performed without the use of any dictionary, and assumes that the number of constituent materials in the scene and their spectral signatures are unknown. The estimated abundances satisfy the desired sum-to-one and nonnegativity constraints. Two models with increasing complexity are developed to achieve this challenging task, depending on how noise interacts with hyperspectral data. The first one leads to a convex optimization problem and is solved with the alternating direction method of multipliers. The second one accounts for signal-dependent noise and is addressed with a reweighted least squares algorithm. Experiments on synthetic and real data demonstrate the effectiveness of our approach.

[1]  Richard H. Bartels,et al.  Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4] , 1972, Commun. ACM.

[2]  W. Fuller,et al.  Estimation for a Linear Regression Model with Unknown Diagonal Covariance Matrix , 1978 .

[3]  Solution of the matrix equation AX−XB=C , 1986 .

[4]  R. Carroll,et al.  An asymptotic theory for weighted least-squares with weights estimated by replication , 1988 .

[5]  Dimitri P. Bertsekas,et al.  On the Douglas—Rachford splitting method and the proximal point algorithm for maximal monotone operators , 1992, Math. Program..

[6]  Peter M. Hooper.  Iterative weighted least squares estimation in heteroscedastic linear models , 1993 .

[7]  Stéphane Mallat,et al.  Matching pursuits with time-frequency dictionaries , 1993, IEEE Trans. Signal Process..

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

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

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

[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]  Tongxing Lu,et al.  Solution of the matrix equation AX−XB=C , 2005, Computing.

[13]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[14]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[15]  M. Yuan,et al.  Model selection and estimation in regression with grouped variables , 2006 .

[16]  Hairong Qi,et al.  A Maximum Entropy Approach to Unsupervised Mixed-Pixel Decomposition , 2007, IEEE Transactions on Image Processing.

[17]  I. Daubechies,et al.  Iteratively reweighted least squares minimization for sparse recovery , 2008, 0807.0575.

[18]  Michael Elad,et al.  From Sparse Solutions of Systems of Equations to Sparse Modeling of Signals and Images , 2009, SIAM Rev..

[19]  Jean-Yves Tourneret,et al.  Bayesian Estimation of Linear Mixtures Using the Normal Compositional Model. Application to Hyperspectral Imagery , 2010, IEEE Transactions on Image Processing.

[20]  Alfred O. Hero,et al.  Multidimensional Shrinkage-Thresholding Operator and Group LASSO Penalties , 2011, IEEE Signal Processing Letters.

[21]  Patrick L. Combettes,et al.  Proximal Splitting Methods in Signal Processing , 2009, Fixed-Point Algorithms for Inverse Problems in Science and Engineering.

[22]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[23]  Shengli Xie,et al.  Blind Spectral Unmixing Based on Sparse Nonnegative Matrix Factorization , 2011, IEEE Transactions on Image Processing.

[24]  Michael Möller,et al.  A Convex Model for Nonnegative Matrix Factorization and Dimensionality Reduction on Physical Space , 2011, IEEE Transactions on Image Processing.

[25]  Alfonso Fernández-Manso,et al.  Spectral unmixing , 2012 .

[26]  Paul Honeine,et al.  Geometric Unmixing of Large Hyperspectral Images: A Barycentric Coordinate Approach , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[27]  Marian-Daniel Iordache,et al.  Greedy algorithms for pure pixels identification in hyperspectral unmixing: A multiple-measurement vector viewpoint , 2013, 21st European Signal Processing Conference (EUSIPCO 2013).

[28]  Jean-Yves Tourneret,et al.  Nonlinear Spectral Unmixing of Hyperspectral Images Using Gaussian Processes , 2012, IEEE Transactions on Signal Processing.

[29]  Ferréol Soulez,et al.  Exploiting spatial sparsity for multi-wavelength imaging in optical interferometry , 2012, Journal of the Optical Society of America. A, Optics, image science, and vision.

[30]  Antonio J. Plaza,et al.  MUSIC-CSR: Hyperspectral Unmixing via Multiple Signal Classification and Collaborative Sparse Regression , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[31]  Antonio J. Plaza,et al.  Collaborative Sparse Regression for Hyperspectral Unmixing , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[32]  Ben Somers,et al.  A multi-measurement vector approach for endmember extraction in urban environments , 2014 .