Optimization of minimum volume constrained hyperspectral image unmixing on CPU–GPU heterogeneous platform

Hyperspectral unmixing is essential for efficient hyperspectral image processing. Nonnegative matrix factorization based on minimum volume constraint (MVC-NMF) is one of the most widely used methods for unsupervised unmixing for hyperspectral image without the pure-pixel assumption. But the model of MVC-NMF is unstable, and the traditional solution based on projected gradient algorithm (PG-MVC-NMF) converges slowly with low accuracy. In this paper, a novel parallel method is proposed for minimum volume constrained hyperspectral image unmixing on CPU–GPU Heterogeneous Platform. First, a optimized unmixing model of minimum logarithmic volume regularized NMF is introduced and solved based on the second-order approximation of function and alternating direction method of multipliers (SO-MVC-NMF). Then, the parallel algorithm for optimized MVC-NMF (PO-MVC-NMF) is proposed based on the CPU–GPU heterogeneous platform, taking advantage of the parallel processing capabilities of GPUs and logic control abilities of CPUs. Experimental results based on both simulated and real hyperspectral images indicate that the proposed algorithm is more accurate and robust than the traditional PG-MVC-NMF, and the total speedup of PO-MVC-NMF compared to PG-MVC-NMF is over 50 times.

[1]  Qian Du,et al.  High Performance Computing for Hyperspectral Remote Sensing , 2011, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[2]  José M. Bioucas-Dias,et al.  Does independent component analysis play a role in unmixing hyperspectral data? , 2005, IEEE Trans. Geosci. Remote. Sens..

[3]  Yin Zhang,et al.  An Alternating Direction Algorithm for Nonnegative Matrix Factorization , 2010 .

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

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

[6]  Yue Yu,et al.  Minimum distance constrained non-negative matrix factorization for the endmember extraction of hyperspectral images , 2007, International Symposium on Multispectral Image Processing and Pattern Recognition.

[7]  Wenfei Luo,et al.  Parallel implementation of N-FINDR algorithm for hyperspectral imagery on hybrid multiple-core CPU and GPU parallel platform , 2011, International Symposium on Multispectral Image Processing and Pattern Recognition.

[8]  Antonio J. Plaza,et al.  Parallel Hyperspectral Unmixing on GPUs , 2014, IEEE Geoscience and Remote Sensing Letters.

[9]  Chein-I Chang,et al.  A New Growing Method for Simplex-Based Endmember Extraction Algorithm , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[10]  José M. Bioucas-Dias,et al.  Hyperspectral unmixing algorithm via dependent component analysis , 2007, 2007 IEEE International Geoscience and Remote Sensing Symposium.

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

[12]  Wei Xia,et al.  An approach based on constrained nonnegative matrix factorization to unmix hyperspectral data , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[13]  Antonio J. Plaza,et al.  Recent Developments in High Performance Computing for Remote Sensing: A Review , 2011, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[14]  Chong-Yung Chi,et al.  A Convex Analysis-Based Minimum-Volume Enclosing Simplex Algorithm for Hyperspectral Unmixing , 2009, IEEE Transactions on Signal Processing.

[15]  Antonio J. Plaza,et al.  Real-time implementation of remotely sensed hyperspectral image unmixing on GPUs , 2012, Journal of Real-Time Image Processing.

[16]  Yuliya Tarabalka,et al.  Real-time anomaly detection in hyperspectral images using multivariate normal mixture models and GPU processing , 2009, Journal of Real-Time Image Processing.

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

[18]  Antonio J. Plaza,et al.  Real-Time Implementation of the Vertex Component Analysis Algorithm on GPUs , 2013, IEEE Geoscience and Remote Sensing Letters.

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

[20]  Jon Atli Benediktsson,et al.  Hyperspectral Unmixing on GPUs and Multi-Core Processors: A Comparison , 2013, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[21]  Antonio J. Plaza,et al.  Parallel implementation of endmember extraction algorithms using NVidia graphical processing units , 2009, 2009 IEEE International Geoscience and Remote Sensing Symposium.

[22]  Antonio J. Plaza,et al.  Special issue on architectures and techniques for real-time processing of remotely sensed images , 2009, Journal of Real-Time Image Processing.

[23]  Antonio J. Plaza,et al.  Real-Time Implementation of the Pixel Purity Index Algorithm for Endmember Identification on GPUs , 2014, IEEE Geoscience and Remote Sensing Letters.

[24]  Jun Zhou,et al.  Hyperspectral Unmixing via $L_{1/2}$ Sparsity-Constrained Nonnegative Matrix Factorization , 2011, IEEE Transactions on Geoscience and Remote Sensing.

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

[26]  Chih-Jen Lin,et al.  Projected Gradient Methods for Nonnegative Matrix Factorization , 2007, Neural Computation.

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

[28]  V. P. Pauca,et al.  Nonnegative matrix factorization for spectral data analysis , 2006 .

[29]  Roger N. Clark,et al.  The US Geological Survey, digital spectral reflectance library: version 1: 0.2 to 3.0 microns , 1993 .

[30]  Antonio J. Plaza,et al.  GPU Implementation of an Automatic Target Detection and Classification Algorithm for Hyperspectral Image Analysis , 2013, IEEE Geoscience and Remote Sensing Letters.

[31]  Julien Michel,et al.  Remote Sensing Processing: From Multicore to GPU , 2011, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

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

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