New Theory for Unmixing ILL-Conditioned Hyperspectral Mixtures

Hyperspectral unmixing (HU), a blind source separation problem, aims at unambiguiously identifying the spectral signatures of the materials, as well as their abundances, from the measured hyperspectral mixtures. In real hyperspectral scenes, high correlation between the spectral signatures is commonly observed, making HU quite challenging. Although such ill-conditioning is critical for effective HU, it is often ignored in existing HU literature. To the best of our knowledge, existing preconditioning techniques, for reducing the condition number of the signature matrix, were developed based on the pure-pixel assumption, which can, however, be seriously violated in remote sensing. Under a relaxed purity assumption, with respect to the pure-pixel one, this paper proposes novel theory for unmixing ill-conditioned hyperspectral mixtures. Specifically, we exactly identify the John's ellipsoid (i.e., the maximum ellipsoid inscribed in the convex hull of the hyperspectral data vectors) via split augmented Lagrangian shrinkage algorithm (SALSA), and transform this ellipsoid into an Euclidean ball. This transformation brings the data vectors into a new space wherein the corresponding material signature vectors form a regular simplex, which is a very strong prior information. Based on this prior, we design an HU criterion, and prove its perfect identifiability under a very mild sufficient condition. Then, we demonstrate the feasibility of realizing our criterion via non-convex optimization and guarantee a stationary point solution.

[1]  Franz E Schuster,et al.  An arithmetic proof of John’s ellipsoid theorem , 2005, 1207.7246.

[2]  José M. Bioucas-Dias,et al.  An Augmented Lagrangian Approach to the Constrained Optimization Formulation of Imaging Inverse Problems , 2009, IEEE Transactions on Image Processing.

[3]  Antonio J. Plaza,et al.  A Signal Processing Perspective on Hyperspectral Unmixing: Insights from Remote Sensing , 2014, IEEE Signal Processing Magazine.

[4]  Wei-Chiang Li,et al.  Convex Optimization for Signal Processing and Communications: From Fundamentals to Applications , 2017 .

[5]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[6]  José M. Bioucas-Dias,et al.  A variable splitting augmented Lagrangian approach to linear spectral unmixing , 2009, 2009 First Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing.

[7]  Shuzhong Zhang,et al.  Maximum Block Improvement and Polynomial Optimization , 2012, SIAM J. Optim..

[8]  José M. Bioucas-Dias,et al.  Hyperspectral Subspace Identification , 2008, IEEE Transactions on Geoscience and Remote Sensing.

[9]  Wei-Chiang Li,et al.  Identifiability of the Simplex Volume Minimization Criterion for Blind Hyperspectral Unmixing: The No-Pure-Pixel Case , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[10]  Mário A. T. Figueiredo,et al.  Near-infrared hyperspectral unmixing based on a minimum volume criterion for fast and accurate chemometric characterization of counterfeit tablets. , 2010, Analytical chemistry.

[11]  Stephen P. Boyd,et al.  Disciplined Convex Programming , 2006 .

[12]  Sanjeev Arora,et al.  A Practical Algorithm for Topic Modeling with Provable Guarantees , 2012, ICML.

[13]  Nicolas Gillis,et al.  Enhancing pure-pixel identification performance via preconditioning , 2014, 2014 6th Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing (WHISPERS).

[14]  José M. Bioucas-Dias,et al.  Alternating direction algorithms for constrained sparse regression: Application to hyperspectral unmixing , 2010, 2010 2nd Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing.

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

[16]  Yue Wang,et al.  Convex-Optimization-Based Compartmental Pharmacokinetic Analysis for Prostate Tumor Characterization Using DCE-MRI , 2016, IEEE Transactions on Biomedical Engineering.

[17]  Yue Wang,et al.  Detection of Sources in Non-Negative Blind Source Separation by Minimum Description Length Criterion , 2018, IEEE Transactions on Neural Networks and Learning Systems.

[18]  Yue Wang,et al.  Maximum Volume Inscribed Ellipsoid: A New Simplex-Structured Matrix Factorization Framework via Facet Enumeration and Convex Optimization , 2017, SIAM J. Imaging Sci..

[19]  Antonio J. Plaza,et al.  Hyperspectral Unmixing Overview: Geometrical, Statistical, and Sparse Regression-Based Approaches , 2012, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[20]  Chong-Yung Chi,et al.  A Fast Hyperplane-Based Minimum-Volume Enclosing Simplex Algorithm for Blind Hyperspectral Unmixing , 2015, IEEE Transactions on Signal Processing.

[21]  Chia-Hsiang Lin,et al.  A Convex Optimization-Based Coupled Nonnegative Matrix Factorization Algorithm for Hyperspectral and Multispectral Data Fusion , 2018, IEEE Transactions on Geoscience and Remote Sensing.

[22]  David P. Dobkin,et al.  The quickhull algorithm for convex hulls , 1996, TOMS.