Extraction of Endmembers From Hyperspectral Images Using A Weighted Fuzzy Purified-Means Clustering Model

Hyperspectral endmembers are the spectra of pure materials that are responsible for generating the mixed pixels in hyperspectral images (HSIs). Hyperspectral endmember extraction (HEE) is essentially an inverse problem, where the unknown endmembers are inferred from the spectral measurements. Efficient extraction of endmembers in HSI relies on a well-defined generative model that captures key factors in HSI generation process, such as the clustering effect in the spatial domain and the noise heterogeneity effect in the spectral domain. This paper presents a weighted fuzzy purified-means (WFP-means) clustering model for HEE, where the endmembers are modeled as mean vectors of individual classes, and the fractional contributions of individual endmembers, called abundances, are treated as soft class membership. Accordingly, an endmember is estimated as the weighted mean of purified pixels in HSI, while the abundances are estimated as the nonnegative regression coefficients. In contrast to a mixed pixel that consists of multiple endmembers, a “purified pixel” is due to a single endmember. The introduction of the concept of “purified pixels” into the fuzzy clustering model leads to an elegant optimization scheme. Moreover, the proposed model accounts for the noise variance heterogeneity issue, which is essential for achieving unbiased abundance estimation. The proposed method is tested on both simulated and real HSI, in comparison with several other HEE methods. The results demonstrate that the proposed method compares favorably with respect to the referenced methods in terms of both endmember and abundance estimation.

[1]  Derek M. Rogge,et al.  Iterative Spectral Unmixing for Optimizing Per-Pixel Endmember Sets , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[2]  Andreas T. Ernst,et al.  ICE: a statistical approach to identifying endmembers in hyperspectral images , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[3]  Antonio J. Plaza,et al.  Sparse Unmixing of Hyperspectral Data , 2011, IEEE Transactions on Geoscience and Remote Sensing.

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

[5]  Chein-I Chang,et al.  Constrained subpixel target detection for remotely sensed imagery , 2000, IEEE Trans. Geosci. Remote. Sens..

[6]  José M. Bioucas-Dias,et al.  Hyperspectral Unmixing Based on Mixtures of Dirichlet Components , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[7]  Vivek K. Goyal,et al.  Denoising Hyperspectral Imagery and Recovering Junk Bands using Wavelets and Sparse Approximation , 2006, 2006 IEEE International Symposium on Geoscience and Remote Sensing.

[8]  Qian Du,et al.  Foreword to the Special Issue on Spectral Unmixing of Remotely Sensed Data , 2011 .

[9]  Paul D. Gader,et al.  Multiple model endmember detection based on spectral and spatial information , 2010, WHISPERS.

[10]  Ali Mohammad-Djafari,et al.  Bayesian analysis of spectral mixture data using Markov Chain Monte Carlo Methods , 2006 .

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

[12]  J. Boardman Automating spectral unmixing of AVIRIS data using convex geometry concepts , 1993 .

[13]  Antonio J. Plaza,et al.  On the use of spectral libraries to perform sparse unmixing of hyperspectral data , 2010, 2010 2nd Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing.

[14]  J. Boardman,et al.  Mapping target signatures via partial unmixing of AVIRIS data: in Summaries , 1995 .

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

[16]  Paul D. Gader,et al.  PCE: Piecewise Convex Endmember Detection , 2010, IEEE Transactions on Geoscience and Remote Sensing.

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

[18]  R. Bro,et al.  A fast non‐negativity‐constrained least squares algorithm , 1997 .

[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]  Alexander Wong,et al.  K-P-Means: A Clustering Algorithm of K “Purified” Means for Hyperspectral Endmember Estimation , 2014, IEEE Geoscience and Remote Sensing Letters.

[21]  José M. Bioucas-Dias,et al.  Minimum Volume Simplex Analysis: A Fast Algorithm to Unmix Hyperspectral Data , 2008, IGARSS 2008 - 2008 IEEE International Geoscience and Remote Sensing Symposium.

[22]  G. Foody,et al.  Sub-pixel land cover composition estimation using a linear mixture model and fuzzy membership functions , 1994 .

[23]  Paul D. Gader,et al.  A comparison of deterministic and probabilistic approaches to endmember representation , 2010, 2010 2nd Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing.

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

[25]  Paul D. Gader,et al.  Piece-wise convex spatial-spectral unmixing of hyperspectral imagery using possibilistic and fuzzy clustering , 2011, 2011 IEEE International Conference on Fuzzy Systems (FUZZ-IEEE 2011).

[26]  Jean-Yves Tourneret,et al.  Bayesian separation of spectral sources under non-negativity and full additivity constraints , 2009, Signal Process..

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

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

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

[30]  Jan Larsen,et al.  Unmixing of Hyperspectral Images using Bayesian Non-negative Matrix Factorization with Volume Prior , 2011, J. Signal Process. Syst..

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

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

[33]  Donghui Chen,et al.  Nonnegativity constraints in numerical analysis , 2009, The Birth of Numerical Analysis.

[34]  Kanti V. Mardia,et al.  Spatial Classification Using Fuzzy Membership Models , 1988, IEEE Trans. Pattern Anal. Mach. Intell..