Parallel optimization of hyperspectral unmixing based on sparsity constrained nonnegative matrix factorization

Hyperspectral unmixing is a typical problem of blind source separation, which can be solved by nonnegative matrix factorization (NMF). Sparsity based NMF will increase the efficiency of unmixing, but its computational complexity limits the possibility of utilizing it in time-critical applications. In this paper, method of parallel hyperspectral unmixing based on sparsity constrained nonnegative matrix factorization on Graphics Processing Units (CSNMF-GPU) is investigated and compared in terms of both accuracy and speed. The realization of the proposed method using Compute Unified Device Architecture (CUDA) on GPU are described and evaluated. The experimental results comparing with the serial implementations based on both simulated and real hyperspectral data demonstrate the effectiveness of the proposed parallel optimization approach.

[1]  Antonio J. Plaza,et al.  Recent developments in sparse hyperspectral unmixing , 2010, 2010 IEEE International Geoscience and Remote Sensing Symposium.

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

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

[4]  H. Sebastian Seung,et al.  Learning the parts of objects by non-negative matrix factorization , 1999, Nature.

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

[6]  Liang Xiao,et al.  A novel sparsity constrained nonnegative matrix factorization for hyperspectral unmixing , 2012, 2012 IEEE International Geoscience and Remote Sensing Symposium.

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

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

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

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

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

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

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