Journal of Sensors
Volume 2015 (2015), Article ID 256391, 11 pages
http://dx.doi.org/10.1155/2015/256391

Research Article

Undersampled Hyperspectral Image Reconstruction Based on Surfacelet Transform

Lei Liu,¹ Jingwen Yan,² Di Guo,³ Yunsong Liu,⁴ and Xiaobo Qu⁴

¹Department of Mathematics, Shantou University, Shantou 515063, China
²Guangdong Provincial Key Laboratory of Digital Signal and Image Processing, Shantou University, Shantou 515063, China
³School of Computer and Information Engineering, Xiamen University of Technology, Xiamen 361005, China
⁴Department of Electronic Science, Xiamen University, Xiamen 361005, China

Received 7 July 2014; Revised 23 September 2014; Accepted 24 September 2014

Academic Editor: Yongqiang Zhao

Copyright © 2015 Lei Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Hyperspectral imaging is a crucial technique for military and environmental monitoring. However, limited equipment hardware resources severely affect the transmission and storage of a huge amount of data for hyperspectral images. This limitation has the potentials to be solved by compressive sensing (CS), which allows reconstructing images from undersampled measurements with low error. Sparsity and incoherence are two essential requirements for CS. In this paper, we introduce surfacelet, a directional multiresolution transform for 3D data, to sparsify the hyperspectral images. Besides, a Gram-Schmidt orthogonalization is used in CS random encoding matrix, two-dimensional and three-dimensional orthogonal CS random encoding matrixes and a patch-based CS encoding scheme are designed. The proposed surfacelet-based hyperspectral images reconstruction problem is solved by a fast iterative shrinkage-thresholding algorithm. Experiments demonstrate that reconstruction of spectral lines and spatial images is significantly improved using the proposed method than using conventional three-dimensional wavelets, and growing randomness of encoding matrix can further improve the quality of hyperspectral data. Patch-based CS encoding strategy can be used to deal with large data because data in different patches can be independently sampled.

1. Introduction

Typical hyperspectral imaging (HSI) is acquired on satellites and aerospace probes and then transmitted to grounds. The imaging spectrometer can provide tens to hundreds of narrow-band spectral information for each spatial location. The huge amounts of data but scarce equipment hardware resources on satellites and aerospace severely limit the transmission and storage of hyperspectral images [1, 2]. However, the ground receiving side holds very strong processing capability. If we transfer the system complexity from the satellites and aerospace probes to the ground, it is expected to potentially solve the limitation in traditional high spectral sampling.

Traditionally, super-resolution reconstruction is used to improve the spatial resolution of HSI [3, 4] while compression can improve the transmission efficiency [5–8]. We acquire “all” data and then “throw away” most of it in transmission process. Can we just directly measure the part that we need to save the cost of storage? Compressed sensing (CS) has been proposed to solve these contradictions in remote sensing. The theory of CS shows that a sparse signal can be recovered from a relatively small number of linear measurements [9, 10]. The difference between the conventional method and CS based transmission method is shown in Figure 1.

(a)

(b)

Figure 1: Conventional HSI transmission and CS HSI transmission. (a) Conventional sampling, compression, and decompression. (b) Compressed sensing and sparse reconstruction.

CS has been widely used in medical imaging [11, 12] and wireless sensor network [13, 14] to reduce the sampled data in recent years which has been applied in remote sensing. Duarte et al. [15] proposed a single-pixel imaging method. Ma [2] applied CS in single frame imaging. Aravind et al. [16] used ten spectral bands to compare orthogonal matching pursuit with simultaneous orthogonal matching pursuit in reconstruction. Fowler [17] presented compressive-projection PCA to effectively shift the computational burden of PCA from the encoder to the decoder.

According to the CS theory, the reconstruction error is bounded by the sparse approximation error [18]. For a signal , let be the solution to obey CS conditions, a forward transform, and the sparse approximation of in transform domain by retaining only the -largest entries of ; then, the reconstruction error is

where , are small positive constants and is noise power [19]. Equation () implies that a sparser representation will reduce the reconstruction error. An optimal sparsifying transform is always important for sparse image reconstruction to reduce the reconstruction error.

In this study, we introduce the surfacelet transform (ST) to sparsely represent hyperspectral images by making use of the spatial and spectral information. The advantage of ST over wavelet on sparse HSI data reconstruction is demonstrated, where the data are reconstructed by using a fast iterative shrinkage-thresholding algorithm [20]. To further improve the reconstruction performance, a 3-dimensional (3D) random encoding is designed and the Gram-Schmidt orthogonalization is adopted. Finally, a patch-based CS encoding scheme is designed to deal with large size data.

The remainder of this paper is organized as follows. First, ST-based compressive sensing HSI reconstruction method is introduced in Section 2. The experimental results are presented in Section 3. Finally, conclusions are given in Section 4.

2. ST-Based Compressive Sensing HSI Reconstruction

The CS theory comes up with two principles: sparsity, which asserts that the reconstructed signal is sparse with a transform , and incoherence, which requires that the encoding matrix is incoherent with [10].

2.1. Sparse Representation of HSI Using Surfacelet Transform

Wavelet is commonly used as a typical sparse transform for HSI [21–24]. However, wavelet sometimes fails in sparsely representing HSI because it has only 7 directions [25–27]. Surfacelet transform (ST), proposed by Lu and Do [25], can efficiently capture the surface intrinsic geometrical structure within -dimensional signals [25]. It offers directional subbands with decomposition level by combining the multiscale pyramid with the 3-dimensional directional filter banks (3D-DFB). Thus, surfacelet with more directions may help reducing the blocky artifacts caused by orthogonal wavelets [24, 25].

ST is a multiscale version of the 3D-DFB. The input signal first goes through the 3D hourglass filter , which is a three-channel undecimated filter bank. One branch of the three-channel structure of 3D-DFB is given by Figure 2. The output is then fed into a 2D filter bank, denoted by , which operates on the (, ) planes. The tree-structured filter bank produces output subbands, denoted by for . Each output is then fed into another 2D filter bank operating on the (, ) planes. In the end, we get outputs, represented by for and [25]. spatial domain basis images of ST with 2 levels of decomposition are shown in Figure 3.

Figure 2: One branch of the three-channel structure of 3D-DFB.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 3: basis images of ST in 3D spatial domains with 2-level decomposition.

ST is optimal for the singularities, thus providing sparser representation of smooth curves and surface singularities than wavelet. Theoretically, for the singularities of an image , the best -term approximation error using ST has error decay rate of , while this error decay rate is for typical wavelets [28, 29]. Figure 4 shows that a sparser representation is achieved using ST. According to the CS theory [19], the reconstruction error bound is proportional to , where is the forward transform and means the number of preserved coefficients in transform domain. Therefore, ST is expected to reduce the image reconstruction error of HSI.

(a)

(b)

(c)

(d)

(e)

Figure 4: The -term approximation errors and reconstructed images when , (a) Lunar Lake HSI data, (b) the -term approximation error, (c) the original 5th band, (d) sparsely represented image using wavelet ( dB), and (e) sparsely represented image using surfacelet ( dB).

2.2. Incoherence and Gaussian Random Encoding Matrix

Incoherence means a column of the CS encoding matrix must be trying the proliferation in the corresponding sparsity basis . Since this paper focuses on investigating the spatial-spectral sparsity of HSI, a Gaussian random matrix is chosen as because it is incoherent with the entire existing basis [9, 10]. If each spatial image is encoded separately, can be represented by block matrix:

where is an encoding matrix for the band of HSI. When , that is, each 2D image is encoded with the same encoding matrix, then is called 2D random encoding matrix; when , that is, each 2D image is encoded differently, then is called 3D random encoding matrix. Physically, 3D random encoding means that digital micromirror device (DMD) arrays are different for different bands. The performance of both encoding schemes will be discussed in Section 3.

In order to improve the performance of the recovered signal, a Gram-Schmidt orthogonalization (GSO) is used in CS encoding matrix. Given one band random encoding matrix and the th column denote as , the GSO is expressed by

The new matrix is adopted as the CS encoding matrix of one spectral band. The randomness of the matrix is not affected by GSO. In the following, columns of each are processed with GSO.

2.3. Reconstruction Algorithm

For HSI , where represents the th spectral band image, represent the spatial dimensions, and represents the spectral depth of HSI. Let , where is a column vector of the th band HSI with size . The data acquisition model for CS is given by

where is a random encoding matrix for the th spectral band image and is the acquired undersampled data. The sampling ratio is defined as

where , which means the HSI are undersampled.

In this paper, ST is adopted to provide sparser representation of HSI data and is expected to reduce the reconstruction error. Let represent inverse ST and let denote forward ST; CS recovers by solving [18, 19]

where , stands for -norm, and is the regularization parameter which decides the tradeoff between the sparsity and the data fidelity. Many researchers seek for a simple and fast algorithm to solve (), such as the conjugate gradient [30], Bregman iteration [31, 32], low rank reconstruction [33], and other methods. In this paper, we choose FISTA to solve () because of its simplicity and fast convergence, whose convergence rate is , where is an iteration counter [20].

To solve (), the smallest Lipschitz constant [20] of the gradient is . The Lipschitz constant ensures the convergence of algorithm [34]. Taking and as the initial values, the threshold is . For , solutions are found by iterating from () as follows:

The maximum number of iterations of FISTA is set as 200 to achieve stable solutions. The final output is .

2.4. Patch-Based Compressed Sensing with Surfacelet Reconstruction (PCSST)

For a larger HSI data , bigger encoding matrix may exceed the memory of computer or leads long computation time. In this case, a patch-based sampling operation can be used. Let be a patch operator and let divide the th band of into patches, and then the data acquisition model for CS is given by

where denotes measurements of patches in th band of and satisfies , where means random encoding on each block. In our scheme, the encoding on different patches of each band is different and these patches are nonoverlapped.

PCSST allows implementing the proposed method on a larger dataset. It has the potential to allow for parallel computing as shown in Figure 5 and also provides the possibility to reduce the complexity of sensor arrays since the sensing detectors independently sample the data in different patches.

Figure 5: Parallel processing of compressive sampling on HIS.

3. Experimental Results

Experiments are conducted in three aspects. First, better reconstructed signal using ST than using wavelet is demonstrated. Second, improving the randomness of encoding matrix for each spectral band is shown to improve the recovery quality and the reconstruction performance of two encoding schemes is compared. Third, a patch-based CS encoding scheme is designed to deal with large data.

The HSI data is obtained from U.S. AVRIS website, including Moffet Field and Lunar Lake [35]. These data contain 224 spectral bands with spatial size , and every pixel is encoded with 16 bits. Linear interpolation is employed to fix the junk bands [36] and maintain the consistency of spectral lines, and all data are normalized.

To evaluate the performance, the mean-square-error (MSE), peak signal-to-noise ratio (PSNR), and structural similarity (SSIM) [37] are adopted as criteria to measure the reconstruction error. Their definitions are

where denotes the original image, and are mean and standard deviation of , stands for the recovered image, and are mean and standard deviation of ; , are small constant, and is adopted in our work, which is used to avoid instability when either or is very close to zero [37]. Mean SSIM index is adopted in this paper. Simulations run on a dual core 2.5 GHz CPU laptop with 4 GB RAM. A fast ST implementation in C language is adopted.

3.1. Surfacelet and Wavelet with 2D Encoding

To simulate the CS data acquisition, the HSI are undersampled by random encoding matrix . 50% sampled data means for . In order to improve the quality of the recovered signal, Gram-Schmidt orthogonalization is used to obtain the CS encoding matrix. With the 2D encoding matrix (i.e., each band has the same encoding matrix), the ST-based reconstruction is compared with the 3D wavelets-based reconstruction. Daubechies filters “db4” with 2 levels of decomposition is used, and there are 7 directional subbands for 3D wavelet. The multiscale pyramid with 2 decomposition levels is chosen, and there are 12 directional subbands for ST. The complexity of ST is and wavelet is .

In order to well discuss the relationship between and PSNR, we compared the PSNR performance from to as shown in Figure 6. is empirically chosen to give optimal PSNR in our work. As shown in Figure 7, ST significantly improves the PSNR of each spectral band more than 3D wavelet. Edges and curves are better reconstructed using ST than 3D wavelets as shown in Figure 8. At a given spatial location, the spectral line reconstructed using ST is much more consistent with the ground truth than wavelets as shown in Figure 9. Under different sampling ratios, ST achieves much better PSNRs than wavelet as shown in Figure 10, and more advantage of ST is seen at low ratios. Besides, ST shows better performance for the Lunar Lake dataset which has richer geometric textures than Moffet Field dataset. Running time (unit: second) based on ST and WT with different sampled rate is shown in Table 1.

Table 1: Running time (unit: second) based on ST and WT with different sampled rate.

Sampling ratio
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8

ST
458.9
629.7
833.8
1002.9
1224.2
1381.8
1544.9
1758.7

WT
989.8
1115.5
1359.7
1639.9
1871.6
2089.6
2324.1
2452.9

Note: the ST is implemented in C and WT is implemented in MATLAB.

Figure 6: PSNR performance versus different . Note: test is performed on Lunar Lake dataset at sampling ratio .

(a)

(b)

Figure 7: PSNR performance using wavelet and surfacelet when sampling ratio is 0.40. (a) Moffet Field dataset and (b) Lunar Lake dataset.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 8: The reconstructed 199th band of Moffet Field and 5th band of Lunar Lake using wavelets or surfacelet when sampling ratio . (a) and (d) are original Moffet Field and Lunar Lake bands, (b)  dB, and (e)  dB, are recovered bands using wavelet, and (c)  dB, and (f)  dB, are reconstructed bands using surfacelet.

(a)

(b)

Figure 9: Reconstructed spectral line when sampling ratio . (a) Moffet Field dataset and (b) Lunar Lake dataset.

(a)

(b)

(c)

(d)

Figure 10: PSNR and SSIM at different sampling ratios. (a), (c) PSNR and SSIM for Moffet Field dataset and (b), (d) PSNR and SSIM for Lunar Lake dataset.

3.2. 2D and 3D Encoding Matrix

The performance of 2D encoding and 3D encoding schemes is compared in this section. The recovered spatial images with two encoding schemes are shown in Figure 11. Edges and curves are better reconstructed using 3D encoding than 2D encoding. The recovered spectral lines using 3D encoding are more consistent with ground truth than using 2D encoding as shown in Figure 12. The improvement is more obvious for Moffet Field dataset with textures in low intensities as shown in Figures 11(a) and 12(a). These results imply that increasing the encoding randomness among different bands will achieve better reconstruction.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 11: The reconstructed 199th band of Moffet Field and 5th band of Lunar Lake using 2D encoding and 3D encoding when sampling ratio . (a) and (d) are original Moffet Field and Lunar Lake bands, (b)  dB, and (e)  dB, are recovered bands using 2D encoding matrix for each band and (c)  dB, and (f)  dB, are reconstructed bands using different encoding matrix for each band.

(a)

(b)

Figure 12: Reconstructed spectral line when sampling ratio . (a) Moffet Field dataset and (b) Lunar Lake dataset.

3.3. Reconstruction with PCSST

We tested the methods on Lunar Lake data with size . A reconstructed spectral band is shown in Figure 13. The image structures recovered by surfacelet are much sharper than wavelet.

(a)

(b)

(c)

Figure 13: Patch-based compressed sensing ST reconstruction. (a) The original 5th band of Lunar Lake dataset, (b) the patch-based reconstructed image using wavelet when undersampled rate and the  dB, , and (c) the patch-based reconstructed image using ST when undersampled rate and the  dB, .

4. Conclusions

Compressive sensing is a new sampling theorem. In this paper, the surfacelet transform is introduced into hyperspectral image reconstruction from compressive sampled data. The surfacelet transform is a directional multiresolution transform for 3D data, which is applied to sparsify the hyperspectral images for the first time. Simulations are conducted in three aspects. First, better reconstructed signal using surfacelet than using wavelet is demonstrated. Second, improving the randomness of encoding matrix for each spectral band is shown to improve the recovery quality. Third, a patch-based CS encoding scheme is designed to deal with large data. It provides the possibility to reduce the complexity of sensor arrays because sensing detectors in different patches independently sample data. Experiments demonstrate that reconstruction of spectral lines and spatial images is significantly improved using the proposed method than using conventional three-dimensional wavelets.

In the future, our work includes the following two aspects.

The first aspect is optimizing sparse representation of hyperspectral imaging, for example, combining adaptive sparse representation [38–40] and surfacelet transform. An adaptive dictionary may provide a sparser representation leading to a lower reconstruction error. Therefore, it is meaningful to try a low-complexity training method in the future.

The second aspect is overlapping patches compressed sensing reconstruction. Overlapping patches can reduce the “block artifacts.” But overlapping also introduces more encoded data and higher computations on image reconstructions. Trading the compressive sampling rate with image qualities using overlapping and speeding up the reconstruction algorithm [41] will be discussed in the future.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors wish to express their thanks to all reviewers for their comments. This work was supported in part by the National Natural Science Foundation of China (no. 61201045, no. 40971206, no. 61302174), Fundamental Research Funds for the Central Universities (no. 2013SH002), Open Fund from Key Lab of Digital Signal and Image Processing of Guangdong Province (no. 54600321, no. 2013GDDSIPL-07), Academic Innovation Team Building Project of Shantou University (no. ITC12002), and Scientific Research Foundation for the Introduction of Talent at Xiamen University of Technology (no. 90030606).

References

A. F. H. Goetz, “Three decades of hyperspectral remote sensing of the Earth: a personal view,” Remote Sensing of Environment, vol. 113, supplement 1, pp. S5–S16, 2009.

J. Ma, “Single-Pixel remote sensing,” IEEE Geoscience and Remote Sensing Letters, vol. 6, no. 2, pp. 199–203, 2009.

Y. Zhao, J. Yang, and J. C. W. Chan, “Hyperspectral imagery super-resolution by spatial-spectral joint nonlocal similarity,” IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, vol. 99, pp. 2671–2679, 2013.

Y. Zhao, J. Yang, Q. Zhang, L. Song, Y. Cheng, and Q. Pan, “Hyperspectral imagery super-resolution by sparse representation and spectral regularization,” EURASIP Journal on Advances in Signal Processing, vol. 2011, article 87, 2011.

B. Aiazzi, L. Alparone, S. Baronti, C. Lastri, and M. Selva, “Spectral distortion in lossy compression of hyperspectral data,” Journal of Electrical and Computer Engineering, vol. 2012, Article ID 850637, 8 pages, 2012.

J. G. Kolo, S. A. Shanmugam, D. W. G. Lim, L.-M. Ang, and K. P. Seng, “An adaptive lossless data compression scheme for wireless sensor networks,” Journal of Sensors, vol. 2012, Article ID 539638, 20 pages, 2012.

Q. Du and J. E. Fowler, “Hyperspectral image compression using JPEG2000 and principal component analysis,” IEEE Geoscience and Remote Sensing Letters, vol. 4, no. 2, pp. 201–205, 2007.

C.-I. Chang, B. Ramakishna, J. Wang, and A. Plaza, “Exploitation-based hyperspectral image compression,” Journal of Applied Remote Sensing, vol. 4, Article ID 041760, 2010.

E. J. Candès, J. Romberg, and T. Tao, “Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information,” IEEE Transactions on Information Theory, vol. 52, no. 2, pp. 489–509, 2006.

D. L. Donoho, “Compressed sensing,” IEEE Transactions on Information Theory, vol. 52, no. 4, pp. 1289–1306, 2006.

X. Qu, Y. Hou, F. Lam, D. Guo, J. Zhong, and Z. Chen, “Magnetic resonance image reconstruction from undersampled measurements using a patch-based nonlocal operator,” Medical Image Analysis, vol. 18, no. 6, pp. 843–856, 2014.

W. Hao, J. Li, X. Qu, and Z. Dong, “Fast iterative contourlet thresholding for compressed sensing MRI,” Electronics Letters, vol. 49, no. 19, pp. 1206–1208, 2013.

J. Xiong and Q. Tang, “1-bit compressive data gathering for wireless sensor networks,” Journal of Sensors, vol. 2014, Article ID 805423, 8 pages, 2014.

C. Caione, D. Brunelli, and L. Benini, “Compressive sensing optimization for signal ensembles in WSNs,” IEEE Transactions on Industrial Informatics, vol. 10, no. 1, pp. 382–392, 2014.

M. F. Duarte, M. A. Davenport, D. Takbar et al., “Single-pixel imaging via compressive sampling: building simpler, smaller, and less-expensive digital cameras,” IEEE Signal Processing Magazine, vol. 25, no. 2, pp. 83–91, 2008.

N. V. Aravind, K. Abhinandan, V. V. Acharya, and D. S. Sumam, “Comparison of OMP and SOMP in the reconstruction of compressively sensed hyperspectral images,” in Proceedings of the International Conference on Communications and Signal Processing (ICCSP '11), pp. 188–192, Hamirpur, India, February 2011.

J. E. Fowler, “Compressive-projection principal component analysis,” IEEE Transactions on Image Processing, vol. 18, no. 10, pp. 2230–2242, 2009.

E. J. Candès, “The restricted isometry property and its implications for compressed sensing,” Comptes Rendus Mathematique, vol. 346, no. 9-10, pp. 589–592, 2008.

E. J. Candès, Y. C. Eldar, D. Needell, and P. Randall, “Compressed sensing with coherent and redundant dictionaries,” Applied and Computational Harmonic Analysis, vol. 31, no. 1, pp. 59–73, 2011.

A. Beck and M. Teboulle, “A fast iterative shrinkage-thresholding algorithm for linear inverse problems,” SIAM Journal on Imaging Sciences, vol. 2, no. 1, pp. 183–202, 2009.

B. Rasti, J. R. Sveinsson, M. O. Ulfarsson, and J. A. Benediktsson, “Hyperspectral image denoising using 3D wavelets,” in Proceedings of the 32nd IEEE International Geoscience and Remote Sensing Symposium (IGARSS '12), pp. 1349–1352, July 2012.

X. Tang and W. A. Pearlman, “Three-dimensional wavelet-based compression of hyperspectral images,” Hyperspectral Data Compression, pp. 273–308, 2006.

S. Prasad, W. Li, J. E. Fowler, and L. M. Bruce, “Information fusion in the redundant-wavelet-transform domain for noise-robust hyperspectral classification,” IEEE Transactions on Geoscience and Remote Sensing, vol. 50, no. 9, pp. 3474–3486, 2012.

A. Karami, M. Yazdi, and G. Mercier, “Compression of hyperspectral images using discerete wavelet transform and tucker decomposition,” IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, vol. 5, no. 2, pp. 444–450, 2012.

Y. M. Lu and M. N. Do, “Multidimensional directional filter banks and surfacelets,” IEEE Transactions on Image Processing, vol. 16, no. 4, pp. 918–931, 2007.

X. Qu, W. Zhang, D. Guo, C. Cai, S. Cai, and Z. Chen, “Iterative thresholding compressed sensing MRI based on contourlet transform,” Inverse Problems in Science and Engineering, vol. 18, no. 6, pp. 737–758, 2010.

P. Feng, Y. Pan, B. Wei, W. Jin, and D. Mi, “Enhancing retinal image by the contourlet transform,” Pattern Recognition Letters, vol. 28, no. 4, pp. 516–522, 2007.

M. N. Do and M. Vetterli, “The contourlet transform: An efficient directional multiresolution image representation,” IEEE Transactions on Image Processing, vol. 14, no. 12, pp. 2091–2106, 2005.

E. J. Candès and D. L. Donoho, “New tight frames of curvelets and optimal representations of objects with piecewise C ² singularities,” Communications on Pure and Applied Mathematics, vol. 57, no. 2, pp. 219–266, 2004.

M. Lustig, D. Donoho, and J. M. Pauly, “Sparse MRI: the application of compressed sensing for rapid MR imaging,” Magnetic Resonance in Medicine, vol. 58, no. 6, pp. 1182–1195, 2007.

W. Yin, S. Osher, D. Goldfarb, and J. Darbon, “Bregman iterative algorithms for -minimization with applications to compressed sensing,” SIAM Journal on Imaging Sciences, vol. 1, no. 1, pp. 143–168, 2008.

J.-F. Cai, S. Osher, and Z. Shen, “Linearized Bregman iterations for compressed sensing,” Mathematics of Computation, vol. 78, no. 267, pp. 1515–1536, 2009.

Y. Q. Zhao and J. Yang, “Hyperspectral image denoising via sparse representation and low-rank constraint,” IEEE Transactions on Geoscience and Remote Sensing, vol. 53, pp. 296–308, 2015.

K. Goebel and W. A. Kirk, “A fixed point theorem for transformations whose iterates have uniform Lipschitz constant,” Studia Mathematica, vol. 47, pp. 135–140, 1973.

AVIRIS Homepage, JPL, NASA, December 2010, http://www.nasa.gov/

P. Zhong and R. Wang, “Multiple-spectral-band CRFs for denoising junk bands of hyperspectral imagery,” IEEE Transactions on Geoscience and Remote Sensing, vol. 51, no. 4, pp. 2260–2275, 2013.

Z. Wang, A. C. Bovik, H. R. Sheikh, and E. P. Simoncelli, “Image quality assessment: From error visibility to structural similarity,” IEEE Transactions on Image Processing, vol. 13, no. 4, pp. 600–612, 2004.

E. le Pennec and S. Mallat, “Sparse geometric image representations with bandelets,” IEEE Transactions on Image Processing, vol. 14, no. 4, pp. 423–438, 2005.

X. Qu, D. Guo, B. Ning et al., “Undersampled MRI reconstruction with patch-based directional wavelets,” Magnetic Resonance Imaging, vol. 30, no. 7, pp. 964–977, 2012.

M. Aharon, M. Elad, and A. Bruckstein, “K-SVD: an algorithm for designing overcomplete dictionaries for sparse representation,” IEEE Transactions on Signal Processing, vol. 54, no. 11, pp. 4311–4322, 2006.

J. Yang, Y. Yamaguchi, and W.-M. Boerner, “Numerical methods for solving the optimal problem of contrast enhancement,” IEEE Transactions on Geoscience and Remote Sensing, vol. 38, no. 2, pp. 965–971, 2000.