Sparse nonnegative matrix underapproximation and its application to hyperspectral image analysis

Dimensionality reduction techniques such as principal component analysis (PCA) are powerful tools for the analysis of high-dimensional data. In hyperspectral image analysis, nonnegativity of the data can be taken into account, leading to an additive linear model called nonnegative matrix factorization (NMF), which improves interpretability of the decomposition. Recently, another technique based on under-approximations (NMU) has been introduced, which allows the extraction of features in a recursive way, such as PCA, but preserving nonnegativity, such as NMF. However, for difficult hyperspectral datasets, even NMU can mix some materials together, and is therefore not able to separate of all materials properly for accurate target identification. In this paper we introduce sparse NMU by adding a sparsity constraint on the abundance matrix and use it to extract materials individually in a more efficient way than NMU. This is experimentally demonstrated on a HYDICE image of the San Diego airport.

[1]  Naum Zuselevich Shor,et al.  Minimization Methods for Non-Differentiable Functions , 1985, Springer Series in Computational Mathematics.

[2]  Luigi Grippo,et al.  On the convergence of the block nonlinear Gauss-Seidel method under convex constraints , 2000, Oper. Res. Lett..

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

[4]  Stanley Osher,et al.  L1 unmixing and its application to hyperspectral image enhancement , 2009, Defense + Commercial Sensing.

[5]  Nicolas Gillis Nonnegative matrix factorization : complexity, algorithms and applications , 2011 .

[6]  Pablo A. Parrilo,et al.  Guaranteed Minimum-Rank Solutions of Linear Matrix Equations via Nuclear Norm Minimization , 2007, SIAM Rev..

[7]  Hualiang Li,et al.  Non-negative Matrix Factorization with Orthogonality Constraints and its Application to Raman Spectroscopy , 2007, J. VLSI Signal Process..

[8]  Hairong Qi,et al.  Endmember Extraction From Highly Mixed Data Using Minimum Volume Constrained Nonnegative Matrix Factorization , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[9]  M. J. D. Powell,et al.  On search directions for minimization algorithms , 1973, Math. Program..

[10]  Laurence A. Wolsey,et al.  Two “well-known” properties of subgradient optimization , 2009, Math. Program..

[11]  Guillermo Sapiro,et al.  Learning discriminative sparse models for source separation and mapping of hyperspectral imagery , 2010 .

[12]  Hyunsoo Kim,et al.  Sparse Non-negative Matrix Factorizations via Alternating Non-negativity-constrained Least Squares , 2006 .

[13]  Xinjian Chen,et al.  Nonlinear band expansion and nonnegative matrix underapproximation for unsupervised segmentation of a liver from a multi-phase CT image , 2011, Medical Imaging.

[14]  Nicolas Gillis,et al.  Using underapproximations for sparse nonnegative matrix factorization , 2009, Pattern Recognit..

[15]  Nicolas Gillis,et al.  Dimensionality reduction, classification, and spectral mixture analysis using non-negative underapproximation , 2011 .

[16]  Qiang Zhang,et al.  Joint segmentation and reconstruction of hyperspectral data with compressed measurements. , 2011, Applied optics.

[17]  Yahya M. Masalmah,et al.  A full algorithm to compute the constrained positive matrix factorization and its application in unsupervised unmixing of hyperspectral imagery , 2008, SPIE Defense + Commercial Sensing.

[18]  Mark D. Plumbley,et al.  Theorems on Positive Data: On the Uniqueness of NMF , 2008, Comput. Intell. Neurosci..

[19]  Nicolas Gillis,et al.  Dimensionality reduction, classification, and spectral mixture analysis using nonnegative underapproximation , 2010, Defense + Commercial Sensing.

[20]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

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