Research Article  Open Access
Ran Li, Hongbing Liu, Yu Zeng, Yanling Li, "Block Compressed Sensing of Images Using Adaptive Granular Reconstruction", Advances in Multimedia, vol. 2016, Article ID 1280690, 9 pages, 2016. https://doi.org/10.1155/2016/1280690
Block Compressed Sensing of Images Using Adaptive Granular Reconstruction
Abstract
In the framework of block Compressed Sensing (CS), the reconstruction algorithm based on the Smoothed Projected Landweber (SPL) iteration can achieve the better ratedistortion performance with a low computational complexity, especially for using the Principle Components Analysis (PCA) to perform the adaptive hardthresholding shrinkage. However, during learning the PCA matrix, it affects the reconstruction performance of Landweber iteration to neglect the stationary local structural characteristic of image. To solve the above problem, this paper firstly uses the Granular Computing (GrC) to decompose an image into several granules depending on the structural features of patches. Then, we perform the PCA to learn the sparse representation basis corresponding to each granule. Finally, the hardthresholding shrinkage is employed to remove the noises in patches. The patches in granule have the stationary local structural characteristic, so that our method can effectively improve the performance of hardthresholding shrinkage. Experimental results indicate that the reconstructed image by the proposed algorithm has better objective quality when compared with several traditional ones. The edge and texture details in the reconstructed image are better preserved, which guarantees the better visual quality. Besides, our method has still a low computational complexity of reconstruction.
1. Introduction
Nyquist frequency sampling theorem is the theoretical basis of traditional image coding (such as JPEG and JPEG 2000), which requests that the number of image transformations is the total number of pixels at least. The full transformation of image results in a large amount of calculations, but most of transformation coefficients are discarded during encoding, which causes the waste of energy consumption. Due to a high computational complexity and a low utilization rate of energy consumption, the traditional image coding is not suitable for the application requiring a light load, for example, wireless sensor network node with the limited energy consumption [1] and remote sensing imaging [2]. In addition, since a few transformation coefficients contain most of useful information, the quality of image reconstruction can be degraded severely once these important coefficients lose; thus it is also challenging for the faulttolerant mechanism in the wireless communication [3]. Compressed Sensing (CS) [4, 5] transforms the signal with the subNyquist rate, and it can still accurately recover the signal, which motivates the image CS as a new image coding scheme [6]. The CS random sampling can be regarded as the partial image transformation, and it compresses image by dimensionality reduction while transforming image. Due to the advantage of saving mostly the costs of image coding, the image CS attracts lots of academic interests.
Many researchers are dedicated to improving the ratedistortion (RD) performance of image CS; the popular method explores the adaptive sparse representation model to improve the convergence performance of minimum norm reconstruction [7, 8]; for example, Chen et al. [9] use the spatial correlation to establish MultipleHypotheses (MH) prediction model and recover the more sparse residual to improve the reconstruction quality; Becker et al. [10] proposed the NESTA algorithm based on the firstorder analysis to solve the minimum Total Variation (TV) model, which ensures the speediness and robustness of the sparse decomposition; Zhang et al. [11] exploit the nonlocal selfsimilarity existing in the group of blocks to express the concise sparse representation model; Wu et al. [12] introduce the local autoregressive (AR) model to trace the nonstationary statistical property of image. For abovementioned methods, the improvement of RD performance is at the expense of the high computational complexity, which results in the fact that the reconstruction time rises rapidly as the dimensionality of image increases; for example, the AR algorithm proposed by Wu et al. [12] requires 1 hour to reconstruct an image with size of 512 × 512 pixels, which makes this method lose the value of application. Compared with the above methods which ignore the computational complexity but pursue the high reconstruction quality, Gan [13] and Mun and Fowler [14] proposed the Smoothed Projected Landweber (SPL) algorithm with a low computational complexity to ensure the better RD performance. However, the SPL algorithm adopts a fixed sparse representation basis (e.g., DCT and wavelet basis) which cannot change adaptively according to the image content in the process of iteration; therefore the potential to improve the reconstruction quality is not developed completely. The Principal Component Analysis (PCA) [15] is the optimal orthogonal transformation matrix to remove spatial redundancy of image; our previous work [16] uses PCA to update continually the sparse representation basis at each SPL iteration. Adapted by the image statistical property, the PCAbased SPL algorithm guarantees the performance improvement. The defect of our previous work is that we ignore the stationary local structural characteristic of image when learning PCA matrix. This defect suppresses the performance improvement of PCA decorrelation, which results in the degradation of reconstructed image by the SPL algorithm.
Aimed at the problem that PCA cannot exploit the stationary local structural characteristic of image to update the sparse representation basis at each SPL iteration, this paper proposes the adaptive reconstruction algorithm based on Granular Computing (GrC) theory [17]. The proposed method divides the image into several granules, in which any granule is a set of image patches with the similar structural characteristic. We use PCA to learn the corresponding optimal representation basis of each granule so as to obtain the efficient hardthresholding shrinkage. Experimental results show that the proposed algorithm can improve the RD performance of image CS and achieve the better subjective visual quality.
This paper will explore the adaptive granular reconstruction in block CS of images. First of all, we briefly introduce the background of block CS and GrCbased clustering, and then we describe the adaptive reconstruction algorithm in detail. Afterward, the experimental results are shown and discussed, and finally we make a conclusion on the paper.
2. Background
2.1. Image Block CS
When applied in image acquisition, CS has high time and space complexity, which may be computationally expensive for some realworld systems. Hence, most of existing image CS approaches [9, 13, 14, 16] split the images into nonoverlapping blocks, and, at the decoder, each block is recovered independently in the recovery algorithm. To the best of our knowledge, the pioneering work of block CS was proposed by Gan [13]. This work was further extended by Gan [13] and Mun and Fowler [14] to improve the performance of block CS. Block CS framework is described as follows. First of all, an image with pixels is divided into small blocks with size of pixels each. Let represent the vectorized signal of the th block (; ) through raster scanning. Construct random Gaussian matrix, and get orthogonal matrix by orthonormalizing the Gaussian matrix. Then, we randomly pick rows from to produce the block measurement matrix . Finally, the corresponding observation vector of each block is obtained by measuring with as follows:There are CS observed values for the original image (), and the total measurement matrix is a block diagonal matrix, in which each diagonal element is ; that is,The block CS can compute an initial image to accelerate the reconstruction process by adopting the linear Minimum Mean Square Error (MMSE) criterion. The MMSE solution of image block can be computed as follows:in which ranges from 0.9 to 1. We set and to be 0.95 and 32, respectively, by experience.
When reconstructing image, we can construct the relationship model between the block observation vector and a whole image : in which is the elementary matrix to rearrange the column vectors block by block to a raster scanning column vector of image, and . The minimum norm reconstruction model can be constructed by using (4) as follows: in which and denote and norm, respectively, is the sparse representation basis, and is the fixed regularization factor. Equation (5) can be solved by using the SPL algorithm with a low computational complexity; in particular can be updated flexibly at each SPL iteration. Therefore, the PCA can be used to learn the adaptive sparse representation basis.
2.2. GrCBased Clustering
GrC is proposed by Zadeh in 1979, and he said that the information granule exists widely in our daily life and it can be viewed as an abstraction of actual object. Granulating, organization, and causality are three important concepts of GrC [18]; that is, the granulating divides a whole of the object into several parts, the organization combines several parts into a whole by some specific operators, and the causality is used to model the relationship between the cause and its result. In the above, there is a natural relation between GrC and clustering. The training set can be regarded as the object, and each sample in training set is defined as a single granule. From the above views, we can see that the granulating and organization correspond to the clustering and its inverse, respectively, and the causality describes the internal relation between samples. Different from the traditional clustering method (e.g., means [19]), the GrCbased clustering chooses a variety of shapes to flexibly represent a granule, and the relationship between granules can be modeled by using the more mature algebraic system; therefore the GrCbased clustering has better generalization ability.
Although the group of points in training set has irregular shape, their borders can be better distinguished by GrCbased clustering, since GrC can express the granule as various shapes, for example, hyperdiamond, hypersphere, and hyperbox. The vector form of granule is denoted as , in which is the center of granule and is the granularity. Denote as a point contained in the granule , and the granularity of is computed byin which is norm to measure the distance. The different distance measure represents the different shape of granule; for example, expresses the corresponding granule as hyperdiamond granule, expresses the corresponding granule as hypersphere, and expresses the corresponding granule as hyperbox granule. The operators between granules include the merge operator and the decomposing operator . The merge operator is to merge two smaller granules into a big one, and the decomposition operator is to decompose a big granule into two smaller ones, in which the merge operator is used to cluster the samples in the training set. The merge granules of two granules and are merged as a new granule:in whichFigure 1 illustrates the merging of two spherical granules 1.0, 2.0, 2.0) and 2.5, 2.5, 2.5) in the twodimensional space, and their merging result is 1.75, 2.25, 1.75).
According to the vectorized representation of granule, the relation between granules inevitably involves the relation between vectors, but the partial order relation between vectors is not consistent with inclusion relation between granules. Therefore, we regard that the fuzzy inclusion relation exists between granules. When merging granules, we measure the fuzzy inclusion relation as follows:in which is a mapping from the space of granule to the space of real numbers, and this mapping uses the following formula:
Aimed to the training set in the clustering problem, we construct an algebra system , which is made up of the granule set , the fuzzy inclusion relation , and merge operator . By the fuzzy inclusion relation between granules, we control the granule merging and divide gradually into a number of granules in order for each granule to be corresponding to a category. Figure 2 shows the GrCbased clustering results of spherical and square granules in the twodimensional space; we can see that the GrC performs the granule merging by using the fuzzy inclusion relation and the granularity, which guarantees the better clustering performance; in particular the control of inclusion relation is very beneficial to distinguish the ownership of the training samples in the border. In a natural image signal, some relations exist between image blocks in the spatial domain. Depending on this experience, the GrCbased clustering divides automatically a whole image into several regions including blocks according to the image content. In each cluster, the structure feature of any block is very similar to others; thus it is more conducive to find the sparse representation basis adapted to the fixed data pattern.
(a) Original training set
(b) Clustering results
3. Adaptive Granular Reconstruction
3.1. SPL Iteration
In the framework of block CS, the convex optimization algorithm (e.g., NESTA [10]) can be only used to independently reconstruct each block. However, the varying block sparsity results in the uneven qualities of reconstructed blocks due to the different local structure characteristics of image. The uneven reconstructed qualities bring about lots of blocking artifacts. In addition, the small block size restricts the performance of block reconstruction, which makes the reconstructed image contain a lot of noises. Aimed to the above problem, Gan [13] proposed the image reconstruction algorithm based on SPL iteration by combining Projection Onto Convex Set (POCS), hardthresholding shrinkage, and wiener filtering. The SPL algorithm guarantees the better reconstruction performance with a low computational complexity. The flow of SPL algorithm is presented as shown in the following
The Flow of SPL Iteration
Task. Find the optimal solution of model (5). Input: block CS observation vector , ; block measurement matrix ; sparse representation basis ; regularization factor for hard thresholding. Initialization: set the initial MMSE linear estimator , . Main iteration: increment by 1, and apply these steps:(i)Weiner filtering: the Winner filter with 3 × 3 sliding window is performed to process as follows:(ii)POCS: the filtered image is projected into the convex set ; that is,(iii)Hardthresholding shrinkage: transform by using the sparse representation basis , and perform the hardthresholding shrinkage as follows:(iv)POCS: perform the POCS operator again:(v)Stopping rule: stop when , where ; is the stopping tolerance. Output: the reconstructed image .
Due to the presence of blocking artifacts and noise in the MMSE linear estimator, we firstly use the 3 × 3 Wiener filter to reduce the blocking artifacts and smoothen the image. Afterward, the filtered image is projected to the convex set again, and the transformation coefficients on the sparse representation basis are shrunk by the hard thresholding to eliminate further noises. Finally, the POCS is performed again to enforce the image into the convex set . Repeat the above flow until meeting the stopping criterion. In the process of hardthresholding shrinkage, the solution is not sparse due to the existence of noises. We can get the more sparse transformation coefficients by enforcing the noise components to 0.
3.2. GrCBased PCA HardThresholding Shrinkage
Whether the hardthresholding shrinkage can eliminate the noises effectively is important for the reconstruction performance of SPL iteration. Gan [13] performs the hardthresholding operator in the overlapping block DCT and the undecimated wavelet domain, but they cannot better capture the directional characteristics (e.g., point, line, and edge) so as to degrade the reconstructed quality of image. To overcome the above defects, Mun and Fowler [14] proposed the bivariate shrinkage based on directional transform (e.g., contour and doubletree wavelet), which achieves some performance improvements. The abovementioned methods focus the useful information as much as possible on a few transformation coefficients by constructing some specific wavelets and remove lots of transform coefficients with little information to suppress noises. According to [20], when performing the hardthresholding shrinkage, it is more effective to capture the important information by learning the dictionary adapted by image content than by constructing the specific wavelet. However, the dictionary learning algorithms (e.g., KSVD algorithm [21]) have a high computational complexity; therefore the dictionary learning algorithm is not suited to the SPL iteration requiring a low computational complexity. Our previous work [16] proposed the PCAbased hardthresholding shrinkage to suppress noises. These methods extract the image patches from the noisy image as samples and perform the PCA to train the orthogonal transformation matrix which removes the spatial correlation. Since the PCA matrix can better separate useful components and useless noises existing in image, the hardthresholding shrinkage gets some performance improvements. In addition, the PCAbased hardthresholding shrinkage has a low computational complexity; thus it is so suitable to SPL iteration. Due to the fact that the image has only locally stationary statistics and varying structural characteristics globally, the training set composed by global patches will generate a single PCA matrix with the poor ability to extract the sparse coefficients. Depending on the local stationary statistics of image, it is more practical to divide the global patches into several subsets which contain the patches with similar structural characteristics, and PCA is used to generate the sparse representation basis of each subset. These PCA matrices are specialized to some particular structure characteristics; thus they guarantee that all samples in subset have the more sparse representations.
According to the abovementioned viewpoints, we firstly cluster image patches before using PCA to learn the sparse representation basis. means [19] is a popular clustering algorithm, but it requires setting the number of clusters in advance, so that the training set cannot be classified adaptively in terms of the data pattern. Due to the nonstationary statistical characteristics of image, a large difference exists between the structure features of patches; thus means is not appropriate for the classification of image patches. The GrC does not need to set a fixed number of clusters, and it can achieve the suitable number of clusters by controlling the threshold of granularity to merge the granules depending on the fuzzy inclusion relation between granules. Therefore, we use the GrC to cluster image patches according to the data pattern of training set. The steps of GrCbased PCA hardthresholding shrinkage are described as follows.
Step 1. We extract all patches with the size of from the noisy image and denote these patches as . The upperleft pixel of patch is corresponding to the th pixel; thus each pixel is corresponding to a patch. There are overlapping pixels between patches, and those pixels beyond the border of image cannot be extracted; that is, .
Step 2. We use as the finest atom to constitute the training set and set the control parameter of granularity. The merge operator is used to merge the granule conditionally until all training samples are contained in some granule, and finally each granule corresponds to a category. The flow of clustering is shown in the following.
The Flow of Dividing the Training Set Using GrC Clustering
Task. Merge the samples in into several granules according to the specified conditions. Input: the training set and the threshold of granularity. Initialization: use the samples in to induce granule atoms; that is, the sample corresponds to the hyperbox granule atom , in which the first entry is the starting point, the second entry is the terminal point, and the third entry is the granularity. In the above, we obtain the initial hyperbox granule set . Main iteration: increment by 1 to , and apply these steps:(i)Compute the fuzzy inclusions between and other granules of ; that is, .(ii)Select the index of granule with the maximum fuzzy inclusion; that is, .(iii)Compute the granularity , and judge whether to merge granules and : if , and are merged into a new granule ; otherwise is still a member in .(iv)Update the training set , that is, to remove from .(v)Stop the above iteration once is empty. Output: hyperbox granule set .
By this flow, we can get the granule set composed of granules.
Step 3. We compute the covariance matrix of the th granule as follows:in which is the number of samples in .
Step 4. We compute the eigenvalues of and their corresponding normalized eigenvectors .
Step 5. The eigenvectors are used to construct the orthogonal transformation matrix , and we transform each patch into domain as follows:
Step 6. The hardthresholding shrinkage is performed as follows:
Step 7. The thresholding coefficients are transformed inversely into the pixel domain; that is,
Step 8. We combine all patches into a whole image after shrinking all patches in . Because of the existence of overlapped area between patches, these pixels in overlapped area are computed by averaging all pixel values at the same position.
By the above flow, each granule corresponds to a PCA matrix, and the additional computational complexity is from the GrC clustering. Therefore, the GrCbased PCA learning algorithm maintains a low computational complexity when compared with the computations of learning global PCA matrix.
4. Experimental Results and Analysis
The performances of proposed algorithm are evaluated on the five 512 × 512 test images Lenna, Barbara, Parrot, Leaves, and House including the varying smooth, edge, and texture components. Its parameter settings are described as follows: the regularization factor is 2.5, the stopping tolerance of SPL iteration is 0.01, the patch size is 7, and the threshold of granularity is 0.03. On one hand, we evaluate the performance of hardthresholding shrinkage based on the different sparse representation bases. The proposed algorithm uses the GrCbased PCA basis, which is abbreviated to GrCPCA, and we select the comparing bases including DCT basis, DoubleTree Discrete Wavelet Transform (DDWT) basis, and Global PCA basis (GPCA). On the other hand, we compare the reconstruction performance of proposed algorithm with the MH based SPL (MH_SPL) algorithm [9], the NESTA algorithm with the minimum TV (NESTA_TV) [10], and the Groupbased Sparse Recovery (GSR) algorithm [11]. In all experiments, the block size is 32; the preset total measurement rate is set from 0.2 to 0.6. The Peak SignaltoNoise Ratio (PSNR) is used to evaluate the objective performance, and all PSNR values are averaged over 5 trials, since the quality of reconstruction can vary due to the randomness of random transformation matrix . All experiments run under the following computer configuration: Intel(R) Core(TM) i7 @ 3.40 GHz CPU, 8 GB RAM, Microsoft Windows 7 32 bits, and MATLAB Version 7.6.0.324 (R2008a).
4.1. Performance Comparison of HardThresholding Shrinkage
Figure 3 presents the average RD performance on the five test images of the SPL algorithm with the hardthresholding shrinkages based on DCT, DDWT, GPCA, and GrCPCA. We can see that our GrCPCA basis obtains the highest PSNR value among all bases at any measurement rate, and particularly it has about 1 dB gains compared to the DCT basis. Our GrCPCA basis exploits the local stationary structural characteristics of image; thus the PSNR of GrCPCA basis is 0.2 dB larger than GPCA basis. Figure 4 shows the performance of hardthreshold shrinkage using different clustering algorithms. We can see that the RD performance degrades once using the traditional means algorithm. However, the PSNR value obtains the obvious improvement at any measurement rate when using the GrCbased clustering method; particularly when the measurement rate is 0.2, the PSNR gain of GrCPCA is around 0.6 dB.
4.2. Comparison of Reconstruction Performance
To evaluate the RD performance of our GrCPCAbased SPL algorithm, we select the popular reconstruction algorithms as benchmarks, and they include MH_SPL [9], NESTA_TV [10], and GSR [11]. Table 1 lists the average PSNR value on the five test images when the measurement rate ranges from 0.2 to 0.4. We can see that our method obviously outperforms MH_SPL and the NESTA_TV, and it obtains the PSNR gain of 0.48 dB and 3.27 dB, respectively. However, when compared with the GSR algorithm, the proposed algorithm in this paper deteriorates 0.64 dB on average, which results from the fact that the GSR uses both the local stationary statistics of image and nonlocal similarities of patches. Therefore, the GSR obtains more sparse representation compared to GrCPCA. However, it can be observed from Table 2 that the GSR algorithm has a high computational complexity, and it requires 2542.65 s to reconstruct Lenna on average, but our algorithm requires only 30.20 s. In the above, when considering both the reconstruction quality and computational complexity, the proposed algorithm is more practical than GSR. Compared with the DDWTbased MH_SPL algorithm, our algorithm increases 18.54 s on average, which results from the fact that PCA needs to relearn the sparse representation basis at each iteration, so that the proposed algorithm obtains the PSNR improvement at the cost of computational complexity. Figure 5 shows the visual qualities of reconstructed Lenna using the different algorithms when the measurement rate is 0.3; we can see that the reconstructed image by NESTA_TV contains lots of blocking artifacts, and MH_SPL algorithm eliminates the blocking artifacts but introduces the noises. However, our algorithm gets the subjective visual quality similar to the GSR algorithm. Figure 6 shows the visual qualities of reconstructed remote sensing image Kirtland using the different algorithms when the measurement rate is 0.3. It can be observed that our method provides better visual quality for the remote sensing application when compared the other reconstruction algorithms.


(a) NESTA_TV
(b) MH_SPL
(c) GSR
(d) Proposed
(a) NESTA_TV
(b) MH_SPL
(c) GSR
(d) Proposed
5. Conclusion
In this paper, we use the GrCbased clustering to improve the performance of PCAbased hardthresholding shrinkage in the SPL iteration. Considering the local stationary statistical property, the image is firstly decomposed into some granules. Since GrC can divide adaptively the training set according to the data pattern, the samples in granule have the similar structure. Then, the PCA is performed to learn the sparse representation basis dedicated to the granule, and the patches in granule are shrunk by the hard thresholding using PCA matrix. Finally, all patches merge into a whole image. The experimental results show that the proposed algorithm has better RD performance and guarantees the better subjective visual quality of reconstructed image with a low computational complexity.
Competing Interests
The authors declare that there are no competing interests regarding the publication of this paper.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China under Grants nos. 61501393 and 61572417, in part by the Key Scientific Research Project of Colleges and Universities in Henan Province of China under Grant no. 16A520069, and in part by Youth Sustentation Fund of Xinyang Normal University under Grant no. 2015QN043.
References
 P. Guo, J. Wang, B. Li, and S. Lee, “A variable thresholdvalue authentication architecture for wireless mesh networks,” Journal of Internet Technology, vol. 15, no. 6, pp. 929–936, 2014. View at: Google Scholar
 J. Wu, F. Liu, L. C. Jiao, and X. Wang, “Compressive sensing SAR image reconstruction based on Bayesian framework and evolutionary computation,” IEEE Transactions on Image Processing, vol. 20, no. 7, pp. 1904–1911, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 S. Xie and Y. Wang, “Construction of tree network with limited delivery latency in homogeneous wireless sensor networks,” Wireless Personal Communications, vol. 78, no. 1, pp. 231–246, 2014. View at: Publisher Site  Google Scholar
 E. J. Candes, 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. View at: Publisher Site  Google Scholar  MathSciNet
 D. L. Donoho, “Compressed sensing,” IEEE Transactions on Information Theory, vol. 52, no. 4, pp. 1289–1306, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 C. Deng, W. Lin, B.S. Lee, and C. T. Lau, “Robust image coding based upon compressive sensing,” IEEE Transactions on Multimedia, vol. 14, no. 2, pp. 278–290, 2012. View at: Publisher Site  Google Scholar
 X. Gao, M. Yu, J. Wang, and J. Wei, “l_{0} sparsity for image denoising with local and global priors,” Advances in Multimedia, vol. 2015, Article ID 386134, 9 pages, 2015. View at: Publisher Site  Google Scholar
 S. Xiao, “Deblurring by solving a TV^{p}Regularized optimization problem using split bregman method,” Advances in Multimedia, vol. 2014, Article ID 906464, 11 pages, 2014. View at: Publisher Site  Google Scholar
 C. Chen, E. W. Tramel, and J. E. Fowler, “Compressedsensing recovery of images and video using multihypothesis predictions,” in Proceedings of the 45th Asilomar Conference on Signals, Systems and Computers (ASILOMAR '11), pp. 1193–1198, Pacific Grove, Calif, USA, November 2011. View at: Publisher Site  Google Scholar
 S. Becker, J. Bobin, and E. J. Candes, “NESTA: a fast and accurate firstorder method for sparse recovery,” SIAM Journal on Imaging Sciences, vol. 4, no. 1, pp. 1–39, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 J. Zhang, D. Zhao, and W. Gao, “Groupbased sparse representation for image restoration,” IEEE Transactions on Image Processing, vol. 23, no. 8, pp. 3336–3351, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 X. Wu, W. Dong, X. Zhang, and G. Shi, “Modelassisted adaptive recovery of compressed sensing with imaging applications,” IEEE Transactions on Image Processing, vol. 21, no. 2, pp. 451–458, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 L. Gan, “Block compressed sensing of natural images,” in Proceedings of the International Conference on Digital Signal Processing, pp. 403–406, 2007. View at: Google Scholar
 S. Mun and J. E. Fowler, “Block compressed sensing of images using directional transforms,” in Proceedings of the 2009 16th IEEE International Conference on Image Processing (ICIP '09), pp. 3021–3024, Cairo, Egypt, November 2009. View at: Publisher Site  Google Scholar
 C. A. Deledalle, J. Salmon, and A. Dalalyan, “Image denoising with patch based PCA: local versus global,” in Proceedings of the 22nd British Machine Vision Conference, pp. 1–10, 2011. View at: Google Scholar
 R. Li and X. Zhu, “A PCAbased smoothed projected Landweber algorithm for block compressed sensing image reconstruction,” in Proceedings of the 4th International Conference on Image Analysis and Signal Processing (IASP '12), pp. 1–6, November 2012. View at: Publisher Site  Google Scholar
 L. A. Zadeh, Fuzzy Sets and Information Granulation, North Holland, Amsterdam, The Netherlands, 1979.
 L. A. Zadeh, “Toward a theory of fuzzy information granulation and its centrality in human reasoning and fuzzy logic,” Fuzzy Sets and Systems, vol. 90, no. 2, pp. 111–127, 1997. View at: Publisher Site  Google Scholar  MathSciNet
 C. M. Bishop, Pattern Recognition and Machine Learning, Springer, New York, NY, USA, 2006. View at: Publisher Site  MathSciNet
 M. Elad and M. Aharon, “Image denoising via sparse and redundant representations over learned dictionaries,” IEEE Transactions on Image Processing, vol. 15, no. 12, pp. 3736–3745, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 M. Aharon, M. Elad, and A. Bruckstein, “KSVD: an algorithm for designing overcomplete dictionaries for sparse representation,” IEEE Transactions on Signal Processing, vol. 54, no. 11, pp. 4311–4322, 2006. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Ran Li 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.