Advances in Mathematical Methods for Image and Signal Processing
View this Special IssueResearch Article  Open Access
A Multiresolution Image Completion Algorithm for Compressing Digital Color Images
Abstract
This paper introduces a new framework for image coding that uses image inpainting method. In the proposed algorithm, the input image is subjected to image analysis to remove some of the portions purposefully. At the same time, edges are extracted from the input image and they are passed to the decoder in the compressed manner. The edges which are transmitted to decoder act as assistant information and they help inpainting process fill the missing regions at the decoder. Textural synthesis and a new shearlet inpainting scheme based on the theory of pLaplacian operator are proposed for image restoration at the decoder. Shearlets have been mathematically proven to represent distributed discontinuities such as edges better than traditional wavelets and are a suitable tool for edge characterization. This novel shearlet pLaplacian inpainting model can effectively reduce the staircase effect in Total Variation (TV) inpainting model whereas it can still keep edges as well as TV model. In the proposed scheme, neural network is employed to enhance the value of compression ratio for image coding. Test results are compared with JPEG 2000 and H.264 Intracoding algorithms. The results show that the proposed algorithm works well.
1. Introduction
Image inpainting [1–5] is a method for recovering regions in images whose pixels are distorted or removed in some way. Inpainting methods are commonly based on Partial Differential Equations (PDEs) and Total Variation (TV) models. In PDEs technique [6], pixel values around the region to be inpainted are considered to be the boundary condition for a boundary value problem. Then, a proper equation for interpolating in that area will be solved. Image inpainting has a variety of applications such as text and object removal, denoising, superresolution, digital zooming, fillingin, and compression. Galić et al. [7] use an inpainting technique directly for compression whereas, in previous studies, image inpainting is treated only as a preprocessing step to increase other existing image compression standards. Bugeau et al. [8] offered a working algorithm for image inpainting trying to approximate the global minimum of an energy functional that combines the three fundamental concepts of selfsimilarity, coherence, and propagation. In this method, when the image does not have enough patches to copy from, either because the mask is too spread and the patch size is large or because the mask is placed on a singular location on the image, then the results are poor and although the presence of a geometry term seems to help, it is clearly not enough. The authors did not solve these problems and they did not optimize the search in the patch space.
Chan et al. [6] have proposed TV wavelet inpainting models. The main benefit of TV model is that it can keep the edges very well. But, the method has the drawback called staircase effect. To overcome this defect, we analyze the physical characteristics of TV model [9] and Laplacian operator [10] in local coordinates. At the same time, the traditional wavelets [11] are not very effective in dealing with multidimensional signals containing distributed discontinuities such as edges. To overcome this limitation, one has to use basis elements with much higher directional sensitivity and of various shapes to be able to capture the intrinsic geometrical features of multidimensional phenomena.
In this paper, a new discrete multiscale representation [12] called the Discrete Shearlet Transform (DST) is introduced to perform the inpainting based on Laplacian operator in wavelet domain. This approach, which is based on the shearlet transform, combines the power of multiscale methods with a unique ability to capture the geometry of multidimensional data and is optimally efficient in representing images containing edges. In the proposed algorithm, the correlated portions are identified and they are removed at encoder and filled by image inpainting at the decoder. A Laplacian based inpainting method which employs a DST is presented which can effectively reduce the staircase effect in TV model and can be used to achieve less computing time. Gradient descent back propagation [13] with adaptive learning rate is proposed for compression.
This paper is arranged as follows. First, in Section 2, the statement of the problem is described. In Section 3, the framework of proposed coding scheme is discussed. Specifically, Section 3 shows the ANNs based compression in shearlet domain and the shearlet image inpainting based on Laplacian operator. In Section 4, the experimental results are presented. Section 5 describes conclusion.
2. Statement of the Problem
A standard image model [6] is defined as where is original noise free model; is Gaussian white noise.
The standard wavelet transform of is given by where and is wavelet coefficients; is mother wavelet function.
Damages in the wavelet domain cause loss of wavelet coefficients of on the index region , with represent those wavelet components missing or damaged. The task of inpainting is to restore the missing coefficients in a proper manner, so that the image will have as much information being restored as possible.
For inpainting when we used the traditional wavelets [14–16], they do not deal with multidimensional discontinuities such as edges. Recently, a theory for multidimensional data called Multiscale Geometric Analysis (MGA) has been developed. Many new MGA tools have been proposed such as ridgelet, curvelet, bandlet, and contourlet, which provide higher directional sensitivity than wavelets. Shearlets [17–19], a new approach proposed in this paper, not only possess all the above properties but also are equipped with a rich mathematical structure similar to wavelets, which are associated with a multiresolution analysis. Shearlets form a tight frame at various scales and directions and are optimally sparse in representing image with edges.
The discrete shearlet transform of is given by and the corresponding image inpainting model in shearlet domain is where,
When this model is solved, the use of TV norm can retain sharp edges while reducing noise and other oscillations. But the corresponding EulerLagrange equation is not trivial to compute since it is highly nonlinear and illposed in strong sense. Furthermore, this model has the drawback that is called stair case effect. To overcome these deficiencies, a Laplacian operator [20] is introduced in this new shearlet inpainting model and is where can be adaptively selected based on the local gradient features of images. That is, away from edges, will be approached to 2 to overcome the staircase effect; on the contrary, will be approached to 1 to preserve edges. So this new model can effectively reduce the staircase effect in TV model whereas it can still retain the sharp edges as TV model.
The EulerLagrange equation of the above inpainting model is
The gradient descent flow of (7) is
The above equation is solved by the simple explicit finite difference algorithm. To simplify the formulation, we introduce the standard finite difference notations, such as
the forward differences:
and the backward differences:
We note that it is important to evaluate the nonlinear term, which we denote as
in (8). However, the Laplace operator is defined in the pixel domain. In this paper, we calculate it straightforwardly by transforming the shearlet domain to the pixel domain to compute the Laplace operator and then transform back to the shearlet domain. That is, we calculate the following:
For all, we compute the following: where is a small positive number which is used to prevent the numerical blowup when
Then we compute the curvature projection on the wavelet basis by
3. Framework of Proposed Scheme
The method proposed in this section is based on removing redundancy at the encoder and restoring the removed information using an inpainting method at the decoder. In this algorithm, redundancy removal is performed through detecting texture regions with similar statistical characteristics and dividing the image into homogeneous regions. The overall system with encoder and decoder diagrams is depicted in Figures 1 and 2. In the following subsections, encoder and decoder blocks are discussed separately.
3.1. Design of Encoder
3.1.1. Image Analysis
The input image is subjected to image analysis by extracting edges from the images in order to identify both the structural and textural regions. An input image is divided into blocks. Then they are identified as textural or structural blocks based on their distance from edges. The block is identified as a structural one when it contains more numbers of pixels having very small distance from edges. The remaining blocks are called textural blocks.
The important blocks from both textural regions and structural regions are selected by using different algorithms [21]. The remaining blocks are easily removed during encoding. In this way, the blocks are identified and removed. After the removal of specified blocks, in the original image, the regions are filled with the corresponding DC value or filled with pure red, or green, or blue component, in the case of color images.
3.1.2. ANN Based Image Compression
Gradient descent back propagation algorithm [13] is the widely used algorithm in Artificial Neural Networks (ANNs). The feedforward neural network architecture is capable of approximating most problems with high accuracy and generalization ability. This algorithm mainly focuses on the error correction learning rule. Error propagation consists of two passes through the different layer of the network: a forward pass and a backward pass. In the forward pass, the input vector is applied to the sensory nodes of the network and its effect propagates through the network layer by layer.
Finally, a set of outputs is produced as the actual response of the network. During the forward pass, the synaptic weights of the networks are all fixed. During the back pass, the synaptic weights are all adjusted in accordance with an error correction rule. The actual response of the network is subtracted from the desired response to produce an error signal. This error signal is then propagated through the network against the direction of synaptic conditions. The synaptic weights are adjusted to make the actual response of the network move closer to the desired response.
Procedure for Image Compression. For experiment, feedforward neural network with three layers is selected. Input layer, hidden layer with 16 neurons, and output layer with 64 neurons are introduced in that network. Back propagation algorithm is used for training process. For training the network, 256 × 256 Lena image is selected (see Algorithm 1).

3.2. Design of Decoder
After performing the ANN based compression at encoder, the decompression process is carried out at decoder. Final reconstructed image is obtained through Laplacian inpainting based on DST and texture synthesis. The theory of shearlets is discussed in addition to the two most important modules, namely, shearlet domain Laplacian inpainting and texture synthesis.
3.2.1. Continuous Shearlet Transform
The continuous shearlet transform [17, 18] is a nonisotropic version of the continuous wavelet transform with a superior directional sensitivity. For ,
Each analyzing element is termed shearlet. The frequency tiling of shearlets is shown in Figure 3.
3.2.2. Discrete Shearlet Transform
By sampling the continuous shearlet transform, we can get a discrete transform which is shown in Figure 4 (see Algorithm 2).

3.2.3. Shearlet Domain Laplacian Inpainting
From the pure inpainting perspective, the inpainting problem may be stated as follows. Let be the original image, which is composed by a source area, denoted by , whose pixel values are known, and a target area, denoted by , representing the damaged region to be repaired or a region to be filled in, and this means being inpainted. As shown in Figure 5, these are nonoverlapping areas, that is, and for the boundary between the source and target regions. The simplified architecture for image coding with inpainting is shown in Figure 5 and the pseudocode of proposed work is described as follows: the pseudocode of proposed framework with PDEs inpainting algorithm based onLaplacian operator (see Algorithm 3).

(a)
(b)
3.2.4. Texture Synthesis
The missing texture blocks are filled in with the texture from its surrounding [2]. Let the region to be filled be denoted by . The lost block will now be filled, pixel by pixel, in a raster fashion. Let be a representative template touching the left of a pixel . We proceed to find an estimate of from the available neighborhood, such that a given distance is minimized. As per [2], is a normalized Sum of Squared Differences (SSD) metric. Once such an estimate of is found, we choose the pixel to the immediate right of estimate of as our candidate for . For stochastic textures, the algorithm selects at random one of the pixels neighboring estimate of . The template can be a simple seed block of 3 × 3 pixels. Then, of all possible 3 × 3 blocks in the 8 neighbourhoods, the one with the minimum normalized SSD is found and a pixel to its right is copied into the current pixel in the lost block. This algorithm is considerably fast when using the improvements in [21–25].
4. Performance Evaluation
We illustrate the performance of the proposed algorithm for image compression with image inpainting in shearlet domain using ANNs and compare it with the image inpainting method proposed by Liu et al. [3], The codes are written in MATLAB 2008a.
4.1. Test Conditions and Parameters
The proposed algorithm is tested with color images from USCSIPI and Kodak image database. In all testes, we use shearlet base Laplacian inpainting and we set the parameter values , for noiseless images and for noisy images.
4.2. Overall Performance
Figure 6 shows test image Lena and corresponding results of proposed system. In this test, edges are extracted from the input image and shown in Figure 6(b) and the image with removed blocks (25% removal) is shown in Figure 6(c).
(a)
(b)
(c)
(d)
(e)
Based on the preserved blocks, the shearlet Laplacian inpainting gives results in Figure 6(d). When comparing the restored image in Figure 6(e) with JPEG 2000, proposed scheme saves 36.15% of bits with QP = 75. The comparison of standard images shows up to 50% bits saving which is achieved by the proposed scheme compared to JPEG 2000 and up to 31.61% bit saving compared to edge based image inpainting method proposed by Liu et al. [3]. The bitsaving results are shown in Table 1.

Figure 7 shows the reconstructed images by the proposed scheme with standard H.264/AVC Intracoding. The bitrate saving is also noticeable, shown in Table 2, but not as much as the comparison with JPEG 2000. The proposed scheme can acquire 48.07% bitrate saving compared to the stateoftheart H.264 Intracoding with QP = 24 for an image Kodim19. Bitrate saving of 33.01% is achieved for an image Kodim11 by the proposed scheme when compared with edge based image inpainting method.

(a)
(b)
(c)
(d)
Figure 8 shows the comparison results with JPEG 2000. The first row in Figure 8 shows results of the proposed method and the second row presents JPEG 2000 results. The bits saving of proposed system is indicated in Table 1. The proposed scheme averagely saves 43.54% bits with QP = 75 for the five images shown in Figure 8. Peak SignaltoNoise Ratio (PSNR) is also used to measure the quality of the restored images. It is defined as follows: where is the maximum possible pixel value of the image.
(a)
(b)
(c)
(d)
(e)
Figure 9 gives objective quality comparisons between the proposed scheme, JPEG 2000, and H.264/AVC Intracoding. It can be observed that the objective quality of the proposed scheme outperforms the two standard compression schemes (JPEG 2000 and H.264/AVC Intra) and edge based inpainting algorithm at both low and high bit rates. The proposed algorithm is adapted for subjective quality, and PSNR is not a good measure to evaluate the subjective quality, especially for the resulting images of Laplacian inpainting in shearlet domain.
(a)
(b)
(c)
(d)
(e)
(f)
The reconstructed images with same visual quality regardless of large PSNR difference are shown in Figure 10. Inpainting with edges produces the results which are difficult to differentiate compared with H.264/AVC Intra. But the inpainting method reduces the bit rate because much less number of bits is required for coding edges. The proposed algorithm averagely saves 21.07% bit rate for images having high percentage of texture regions (mandrill and kodim13) when QP = 24 and it is shown in Table 2.
(a)
(b)
The recent noreference quality assessment method is adopted in this paper that is proposed in [26] to evaluate the proposed scheme in comparison with JPEG 2000 and H.264/AVC Intra. Note that in [26], blocking artifacts are detected within a compressed image. From Figure 11, it can be observed that the reconstructed images by the proposed scheme contain less blocking than those by JPEG 2000 and H.264/AVC Intra at similar compression ratios. The result shows the better bitrate reduction when compared with previous methods as far as visual quality is concerned.
(a)
(b)
(c)
(d)
4.3. Computational Complexity
For the experiment, Intel 2 GHz CPU is used. The shearlet domain Laplacian inpainting with the help of assistant information module is realized in MATLAB 2008a. According to our empirical results, the shearlet domain Laplacian inpainting process at the decoder needs averagely 34 iterations to converge. We found that the proposed scheme needs 1–3 s to decode a 768 × 512 image at 40% removal. The time complexity of inpainting is in general proportional to the size of removed regions.
The proposed algorithm is simpler because of the presence of region removal and assistant information generation based on Discrete Shearlet Transform (DST).
5. Conclusion
In this paper, we develop an ANNs based image compression framework that adopts shearlet domain inpainting technique. In this proposed algorithm, the correlated regions are identified and removed automatically at the encoder. Then they are restored at the decoder by using inpainting scheme. The key techniques used for coding are gradient descent back propagation algorithm with adaptive learning rate and Laplacian image inpainting in shearlet domain. Experimental results show that the proposed scheme produces good results. The proposed scheme produces up to 50% and 48.07% bits saving when compared with JPEG 2000 and H.264/AVC Intra, respectively.
Edge extraction can be flexible and adaptable for compression. Finding the regions that can be eliminated is considered to be an open problem and it seems that solving this problem will lead to increase in compression ratios and output quality.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 V. Bastani, M. S. Helfroush, and K. Kasiri, “Image compression based on spatial redundancy removal and image inpainting,” Journal of Zhejiang University: Science C, vol. 11, no. 2, pp. 92–100, 2010. View at: Publisher Site  Google Scholar
 A. A. Efros and T. K. Leung, “Texture synthesis by nonparametric sampling,” in Proceedings of the 7th IEEE International Conference on Computer Vision (ICCV '99), pp. 1033–1038, September 1999. View at: Google Scholar
 D. Liu, X. Sun, F. Wu, S. Li, and Y.Q. Zhang, “Image compression with edgebased inpainting,” IEEE Transactions on Circuits and Systems for Video Technology, vol. 17, no. 10, pp. 1273–1286, 2007. View at: Publisher Site  Google Scholar
 S. D. Rane, G. Sapiro, and M. Bertalmio, “Structure and texture fillingin of missing image blocks in wireless transmission and compression applications,” IEEE Transactions on Image Processing, vol. 12, no. 3, pp. 296–303, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 M. Bertalmio, G. Sapiro, V. Caselles, and C. Ballester, “Image inpainting,” in Proceedings of the 27th ACM International Conference on Computer Graphics and Interactive Techniques (SIGGRAPH '00), pp. 417–424, July 2000. View at: Google Scholar
 T. F. Chan, J. Shen, and H.M. Zhou, “Total variation wavelet inpainting,” Journal of Mathematical Imaging and Vision, vol. 25, no. 1, pp. 107–125, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 I. Galić, J. Weickert, M. Welk, A. Bruhn, A. Belyaev, and H.P. Seidel, “Image compression with anisotropic diffusion,” Journal of Mathematical Imaging and Vision, vol. 31, no. 23, pp. 255–269, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 A. Bugeau, M. Bertalmío, V. Caselles, and G. Sapiro, “A comprehensive framework for image inpainting,” IEEE Transactions on Image Processing, vol. 19, no. 10, pp. 2634–2645, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 M. Lysaker and X.C. Tai, “Iterative image restoration combining total variation minimization and a secondorder functional,” International Journal of Computer Vision, vol. 66, no. 1, pp. 5–18, 2006. View at: Publisher Site  Google Scholar
 H.Y. Zhang, Q.C. Peng, and Y.D. Wu, “Wavelet inpainting based on pLaplace operator,” Acta Automatica Sinica, vol. 33, no. 5, pp. 546–549, 2007. View at: Publisher Site  Google Scholar
 U. A. Ignácio and C. R. Jung, “Blockbased image inpainting in the wavelet domain,” The Visual Computer, vol. 23, no. 9–11, pp. 733–741, 2007. View at: Publisher Site  Google Scholar
 G. R. Easley, D. Labate, and F. Colonna, “Shearletbased total variation diffusion for denoising,” IEEE Transactions on Image Processing, vol. 18, no. 2, pp. 260–268, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 A. Shaik, C. R. K. Reddy, and S. A. Ali, “Empirical analysis of image compression through wave transforms and NN,” International Journal of Computer Science and Information Technologies, vol. 2, pp. 924–931, 2011. View at: Google Scholar
 R. H. Chan, Y.W. Wen, and A. M. Yip, “A fast optimization transfer algorithm for image inpainting in wavelet domains,” IEEE Transactions on Image Processing, vol. 18, no. 7, pp. 1467–1476, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 M. S. Joshi, R. R. Manthalkar, and Y. V. Joshi, “Image compression using curvelet, ridgelet and wavelet transform: a comparative study,” ICGST International Journal on Graphics, Vision and Image Processing (GVIP), vol. 8, pp. 25–31, 2008. View at: Google Scholar
 L. Granai, F. Moschetti, and P. Vandergheynst, “Ridgelet transform applied to motion compensated images,” in Proceedings of the IEEE International Conference on Accoustics, Speech, and Signal Processing, vol. 3, pp. 381–384, April 2003. View at: Google Scholar
 G. Kutyniok and D. Labate, “Resolution of the wavefront set using continuous shearlets,” Transactions of the American Mathematical Society, vol. 361, no. 5, pp. 2719–2754, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. Guo, D. Labate, and W.Q. Lim, “Edge analysis and identification using the continuous shearlet transform,” Applied and Computational Harmonic Analysis, vol. 27, no. 1, pp. 24–46, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Easley, D. Labate, and W.Q. Lim, “Sparse directional image representations using the discrete shearlet transform,” Applied and Computational Harmonic Analysis, vol. 25, no. 1, pp. 25–46, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Li, “Existence of solutions to a superlinear pLaplacian equation,” Electronic Journal of Differential Equations, vol. 2001, no. 66, pp. 1–6, 2001. View at: Google Scholar  MathSciNet
 L. Luo, F. Wu, S. Li, Z. Xiong, and Z. Zhuang, “Advanced motion threading for 3D wavelet video coding,” Signal Processing: Image Communication, vol. 19, no. 7, pp. 601–616, 2004. View at: Publisher Site  Google Scholar
 P. J. Burt and E. H. Adelson, “The Laplacian pyramid as a compact image code,” IEEE Transactions on Communications, vol. 31, no. 4, pp. 532–540, 1983. View at: Google Scholar
 Ö. N. Gerek and A. E. Çetin, “A 2D orientationadaptive prediction filter in lifting structures for image coding,” IEEE Transactions on Image Processing, vol. 15, no. 1, pp. 106–111, 2006. View at: Publisher Site  Google Scholar
 C. Zhu, X. Sun, F. Wu, and H. Li, “Video coding with spatiotemporal texture synthesis and edgebased inpainting,” in Proceedings of the IEEE International Conference on Multimedia and Expo (ICME '08), pp. 813–816, June 2008. View at: Publisher Site  Google Scholar
 Y. Lu, W. Gao, and F. Wu, “Efficient background video coding with static sprite generation and arbitraryshape spatial prediction techniques,” IEEE Transactions on Circuits and Systems for Video Technology, vol. 13, no. 5, pp. 394–405, 2003. View at: Publisher Site  Google Scholar
 F. Pan, X. Lin, S. Rahardja, E. P. Ong, and W. S. Lin, “Measuring blocking artifacts using edge direction information,” in Proceedings of the IEEE International Conference on Multimedia and Expo (ICME '04), vol. 2, pp. 1491–1494, June 2004. View at: Google Scholar
Copyright
Copyright © 2014 R. Gomathi and A. Vincent Antony Kumar. 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.