HYCA: A New Technique for Hyperspectral Compressive Sensing

Hyperspectral imaging relies on sophisticated acquisition and data processing systems able to acquire, process, store, and transmit hundreds or thousands of image bands from a given area of interest. In this paper, we exploit the high correlation existing among the components of the hyperspectral data sets to introduce a new compressive sensing methodology, termed hyperspectral coded aperture (HYCA), which largely reduces the number of measurements necessary to correctly reconstruct the original data. HYCA relies on two central properties of most hyperspectral images, usually termed data cubes: 1) the spectral vectors live on a low-dimensional subspace; and 2) the spectral bands present high correlation in both the spatial and the spectral domain. The former property allows to represent the data vectors using a small number of coordinates. In this paper, we particularly exploit the high spatial correlation mentioned in the latter property, which implies that each coordinate is piecewise smooth and thus compressible using local differences. The measurement matrix computes a small number of random projections for every spectral vector, which is connected with coded aperture schemes. The reconstruction of the data cube is obtained by solving a convex optimization problem containing a data term linked to the measurement matrix and a total variation regularizer. The solution of this optimization problem is obtained by an instance of the alternating direction method of multipliers that decomposes very hard problems into a cyclic sequence of simpler problems. In order to address the need to set up the parameters involved in the HYCA algorithm, we also develop a constrained version of HYCA (C-HYCA), in which all the parameters can be automatically estimated, which is an important aspect for practical application of the algorithm. A series of experiments with simulated and real data shows the effectiveness of HYCA and C-HYCA, indicating their potential in real-world applications.

[1]  J. Tropp,et al.  CoSaMP: Iterative signal recovery from incomplete and inaccurate samples , 2008, Commun. ACM.

[2]  Liangpei Zhang,et al.  Hyperspectral Image Denoising Employing a Spectral–Spatial Adaptive Total Variation Model , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[3]  Stephen J. Wright,et al.  Sparse reconstruction by separable approximation , 2009, IEEE Trans. Signal Process..

[4]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[5]  Michael Elad,et al.  The Cosparse Analysis Model and Algorithms , 2011, ArXiv.

[6]  Javier Portilla,et al.  Image restoration through l0 analysis-based sparse optimization in tight frames , 2009, 2009 16th IEEE International Conference on Image Processing (ICIP).

[7]  Qian Du,et al.  Low-Complexity Principal Component Analysis for Hyperspectral Image Compression , 2008, Int. J. High Perform. Comput. Appl..

[8]  Bormin Huang,et al.  GPU Acceleration of Predictive Partitioned Vector Quantization for Ultraspectral Sounder Data Compression , 2011, IEEE J. Sel. Top. Appl. Earth Obs. Remote. Sens..

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

[10]  Simon Foucart,et al.  Hard Thresholding Pursuit: An Algorithm for Compressive Sensing , 2011, SIAM J. Numer. Anal..

[11]  Pierre Vandergheynst,et al.  Joint trace/TV norm minimization: A new efficient approach for spectral compressive imaging , 2012, 2012 19th IEEE International Conference on Image Processing.

[12]  Balas K. Natarajan,et al.  Sparse Approximate Solutions to Linear Systems , 1995, SIAM J. Comput..

[13]  M. Nikolova An Algorithm for Total Variation Minimization and Applications , 2004 .

[14]  D. L. Donoho,et al.  Compressed sensing , 2006, IEEE Trans. Inf. Theory.

[15]  M. Lustig,et al.  Compressed Sensing MRI , 2008, IEEE Signal Processing Magazine.

[16]  Pierre Vandergheynst,et al.  Compressive Source Separation: Theory and Methods for Hyperspectral Imaging , 2012, IEEE Transactions on Image Processing.

[17]  Rahul Garg,et al.  Gradient descent with sparsification: an iterative algorithm for sparse recovery with restricted isometry property , 2009, ICML '09.

[18]  Antonio J. Plaza,et al.  GPU implementation of a constrained hyperspectral coded aperture algorithm for compressive sensing , 2015, 2015 7th Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing (WHISPERS).

[19]  Mike E. Davies,et al.  Iterative Hard Thresholding for Compressed Sensing , 2008, ArXiv.

[20]  Michael Elad,et al.  Stable recovery of sparse overcomplete representations in the presence of noise , 2006, IEEE Transactions on Information Theory.

[21]  D. Donoho,et al.  Atomic Decomposition by Basis Pursuit , 2001 .

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

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

[24]  J. Tropp JUST RELAX: CONVEX PROGRAMMING METHODS FOR SUBSET SELECTION AND SPARSE APPROXIMATION , 2004 .

[25]  Jessica A. Faust,et al.  Imaging Spectroscopy and the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) , 1998 .

[26]  S. Foucart,et al.  Sparsest solutions of underdetermined linear systems via ℓq-minimization for 0 , 2009 .

[27]  T. Chan,et al.  Fast dual minimization of the vectorial total variation norm and applications to color image processing , 2008 .

[28]  Giovanni Motta Hyperspectral Data Compression , 2006 .

[29]  Emmanuel J. Candès,et al.  Decoding by linear programming , 2005, IEEE Transactions on Information Theory.

[30]  R. DeVore,et al.  A Simple Proof of the Restricted Isometry Property for Random Matrices , 2008 .

[31]  Mário A. T. Figueiredo,et al.  Signal restoration with overcomplete wavelet transforms: comparison of analysis and synthesis priors , 2009, Optical Engineering + Applications.

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

[33]  Yunsong Li,et al.  A GPU-Accelerated Wavelet Decompression System With SPIHT and Reed-Solomon Decoding for Satellite Images , 2011, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[34]  Deanna Needell,et al.  Stable Image Reconstruction Using Total Variation Minimization , 2012, SIAM J. Imaging Sci..

[35]  Antonio J. Plaza,et al.  A new technique for hyperspectral compressive sensing using spectral unmixing , 2012, Optics & Photonics - Optical Engineering + Applications.

[36]  Richard G. Baraniuk,et al.  A new compressive imaging camera architecture using optical-domain compression , 2006, Electronic Imaging.

[37]  Bormin Huang Satellite Data Compression , 2011 .

[38]  Bhaskar D. Rao,et al.  Sparse signal reconstruction from limited data using FOCUSS: a re-weighted minimum norm algorithm , 1997, IEEE Trans. Signal Process..

[39]  A F Goetz,et al.  Imaging Spectrometry for Earth Remote Sensing , 1985, Science.

[40]  Lawrence Carin,et al.  Bayesian Compressive Sensing , 2008, IEEE Transactions on Signal Processing.

[41]  Yin Zhang,et al.  A Compressive Sensing and Unmixing Scheme for Hyperspectral Data Processing , 2012, IEEE Transactions on Image Processing.

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

[43]  R. Tibshirani Regression Shrinkage and Selection via the Lasso , 1996 .

[44]  Michael Elad,et al.  Optimally sparse representation in general (nonorthogonal) dictionaries via ℓ1 minimization , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Ashwin A. Wagadarikar,et al.  Single disperser design for coded aperture snapshot spectral imaging. , 2008, Applied optics.

[46]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[47]  M E Gehm,et al.  Single-shot compressive spectral imaging with a dual-disperser architecture. , 2007, Optics express.

[48]  Qiang Zhang,et al.  Joint segmentation and reconstruction of hyperspectral data with compressed measurements. , 2011, Applied optics.