A Fast Hyperplane-Based Minimum-Volume Enclosing Simplex Algorithm for Blind Hyperspectral Unmixing

• Published in: IEEE transactions on signal processing Vol. 64; no. 8; pp. 1946 - 1961
• Main Authors: Lin, Chia-Hsiang, Chi, Chong-Yung, Wang, Yu-Hsiang, Chan, Tsung-Han
• Format: Journal Article
• Language: English
• Published: New York IEEE 15.04.2016; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Algorithm design and analysis; Algorithms; Blinds; Computation; Computational efficiency; Computer simulation; convex geometry; Craig's criterion; Estimation; Geometry; hyperplane; Hyperplanes; Hyperspectral imaging; Hyperspectral unmixing; Linear matrix inequalities; minimum-volume enclosing simplex; Monte Carlo simulation; Pixels; Signal processing; Signal processing algorithms
• ISSN: 1053-587X; 1941-0476
• DOI: 10.1109/TSP.2015.2508778

Hyperspectral unmixing (HU) is a crucial signal processing procedure to identify the underlying materials (or endmembers) and their corresponding proportions (or abundances) from an observed hyperspectral scene. A well-known blind HU criterion, advocated by Craig during the early 1990s, considers the vertices of the minimum-volume enclosing simplex of the data cloud as good endmember estimates, and it has been empirically and theoretically found effective even in the scenario of no pure pixels. However, such kinds of algorithms may suffer from heavy simplex volume computations in numerical optimization, etc. In this paper, without involving any simplex volume computations, by exploiting a convex geometry fact that a simplest simplex of N vertices can be defined by N associated hyperplanes, we propose a fast blind HU algorithm, for which each of the N hyperplanes associated with the Craig's simplex of N vertices is constructed from N-1 affinely independent data pixels, together with an endmember identifiability analysis for its performance support. Without resorting to numerical optimization, the devised algorithm searches for the N(N-1) active data pixels via simple linear algebraic computations, accounting for its computational efficiency. Monte Carlo simulations and real data experiments are provided to demonstrate its superior efficacy over some benchmark Craig-criterion-based algorithms in both computational efficiency and estimation accuracy.