Research Article  Open Access
Zhao Jian, Sun Meiling, Jia Jian, Huang Luxi, Han Fan, Liu Shan, "Image Watermark Based on Extended Shearlet and Insertion Using the Largest Information Entropy on Horizontal Cone", Mathematical Problems in Engineering, vol. 2015, Article ID 450819, 10 pages, 2015. https://doi.org/10.1155/2015/450819
Image Watermark Based on Extended Shearlet and Insertion Using the Largest Information Entropy on Horizontal Cone
Abstract
Extended discrete shearlet provides a directional multiresolution decomposition. It has been mathematically shown that extended discrete shearlet is a more efficient representation for the signals containing distributed discontinuities such as edges, compared to discrete wavelet. Multiresolution analyses such as curvelet and ridgelet share similar properties, yet their directional representations are significantly different from that of extended discrete shearlet. Taking advantage of the unique properties of directional representation of extended discrete shearlet, we develop an image watermark algorithm based on the largest information entropy. In proposed algorithm, firstly, 1level extended discrete shearlet transform decomposes the test image into directional components on horizontal cone; each directional component reflects directional features and textured features differently. Next, the directional component whose information entropy is the highest is selected to carry watermark. Compared with related algorithms based on DWT and DCT, the proposed algorithm tends to obtain preferable invisibility when it is robust against common attacks.
1. Introduction
Image watermark is one of the most active and challenging subjects in the information hiding research because it is an efficient solution to protect the copyright of the digital media. Wavelets have been widely applied in image watermark owning to their perfect performance for piecewise smooth signals in one dimension [1, 2].
However, wavelets are not actually optimal representation for signals in two dimensions or higher dimensions such as images [3–7]. In fact, wavelets not only fail to preserve the detailed features but they also fail to capture the directional features [3, 7]. With the development of multiresolution analysis, various multiresolution methods have been proposed to overcome this limitation of wavelets, and it has been experimentally proved that these multiresolution methods perform well in image watermark [4, 5, 8]. Among various multiresolution methods, discrete shearlet [6, 9] is a relatively new multiresolution method, which has drawn considerable attention and has been applied in some fields of image processing such as image denoising, image approximations, and image inpainting for its superior performance in dealing with signals in high dimensions [6, 9–11]. Extended discrete shearlet [9, 11] is derived from discrete shearlet transform. It is based on a simple but rigorous mathematical framework which is composed of directional orthonormal bases; each orthonormal base exhibits high directional sensitivity and is highly anisotropic. The extended discrete shearlet transform attempts to obtain a sparse image representation which contains more directional information and preserves more detailed information than discrete wavelet.
Taking advantage of these unique properties of extended discrete shearlet, we develop an image watermark algorithm based on extended discrete shearlet transform and the largest information entropy [12, 13]. In proposed algorithm, by using extended discrete shearlet transform, an image is composed into directional components on the horizontal cone; each directional component is transformed into a set of shearlet coefficients. The watermark is embedded into the largest coefficients of the directional component whose information entropy is the highest. The proposed algorithm is compared with DCTbased algorithm and DWTbased algorithm.
2. Correlation with Basic Theory
2.1. Correlation with Shearlet Transform
A 2D affine system with composite dilations [6, 9, 14] is described bywhere the parabolic scaling matrices are associated with scale transform and the shear matrices are related to areapreserving geometrical transform. The generating function is chosen like a tensor product bywhere is continuous wavelet which is supported with and is a bump function with certain weak additional properties and is supported with .
The associated continuous shearlet transform [6, 9–11, 14] is defined byThe shearlet transform is a function associated with , , and , whereas wavelet transform is a function which is merely associated with and . The main idea of the shear parameter is to capture directional features, such as orientations of curves [6, 9, 11]. As a consequence, the shearlet is able to capture directional information. Figure 1 illustrates the ability of capturing directions of the shearlet transform.
(a)
(b)
(c)
(d)
2.2. Correlation with Discrete Shearlet Transform
By sampling on an appropriate discrete set of , , and , a discrete shearlet transform is obtained [6, 9, 10, 14]. The discrete shearlet transform is described bywhere , , , and . If satisfy the propertythe discrete shearlet transform forms a tight frame for .
For any , , supposeSpecifically, define where is a discrete wavelet function with and is a bump function with . From the supporting condition of and , in frequency domain is supported with
Equation (8) implies that each element is supported on a pair of trapezoids, each of which is contained in a box of approximate size and oriented along line of slope [6, 9, 10, 14].
2.3. Correlation with Extended Discrete Shearlet
Extended discrete shearlet [6, 11] is obtained by directly extending discrete shearlet. Yet its framework is slightly different from that of the discrete shearlet. Framework of the extended discrete shearlet is composed of directional orthonormal bases; each orthonormal basis provides anisotropic window functions elongated along the direction determined by shear matrix. In addition, the extended discrete shearlet framework uses a MRA approach to compute shearlet coefficients and uses an iterative method to approximate the inverse. Based on this construction, a given image is explicitly decomposed into variously directional components for different cones when extended discrete shearlet is associated with different shear matrices and direction numbers; each directional component reflects directional features and textured features differently (see Figure 2); besides, a stable reconstruction of a given image is obtained.
(a)
(b)
(c)
Actually, the extended discrete shearlet framework exhibits highly directional sensitivity; and, furthermore, the extended discrete shearlet framework is optimally sparse representation and is well localized in both the spatial and the frequency domains. Based on these properties, the extended discrete shearlet can be considered as an appropriate transform for image watermark.
3. A Novel Watermark Algorithm Based on Extended Discrete Shearlet
The most important issue in extended discrete shearletbased watermark algorithm is how to choose the watermark embedding location or the effective coefficients to be embedded, since the choice of watermark embedding location has an effect on both invisibility and robustness of watermark algorithm to some extent [1, 2, 4, 5, 8]. In this section, a novel algorithm for image watermark is proposed, which is based on combining the analysis of extended discrete shearlet coefficients [6, 11] and the basic theory of the largest information entropy to search the watermark embedding location for different images.
3.1. Analysis of Extended Discrete Shearlet Coefficients
After 1level extended discrete shearlet decomposition associated with shear matrix and direction numbers , the test image of size is decomposed into directional components on horizontal cone (see Figure 2(b)); each of directional components is transformed into a set of shearlet coefficients.
Each directional component reflects directional features and preserved textured features differently; therefore, the occurring probabilities of different pixel values in each directional component are different from those of all others and the probability distribution of each directional component’s information is different from that of all others in confusion degree. The higher information entropy [12, 13] of the directional component is, the more uncertain occurring probabilities of its different pixel values are, and probability distribution of its information is more disordered and more chaotic; that means it can hide more watermark bits. Based on this analysis, the directional component whose information entropy is the largest is selected to carry watermark. The information entropy of the th directional component on horizontal cone is computed as follows:(1)Mark pixel values of the th directional component by .(2)Count occurring probability of different pixel values of the th directional component; here is the total number of different pixel values.(3)Calculate the information entropy of the th directional component by where .
Generally, watermark insertion is realized by modifying some special values of pixels or transforming domain coefficients [1, 2, 4, 5]. Here watermark insertion is realized by modifying the largest coefficients of selected directional component in extended discrete shearlet domain.
3.2. Watermark Insertion Using the Largest Information Entropy
For an image used as watermark, in order to uniformly and randomly spread bits over the host image, it needs to be transformed into a binary image by Here threshold is selected by OSTA method. After processing the watermark image , a binary image is obtained and then reshaped to binary sequence with length .
According to the above analysis in Section 3.1, the directional component, whose information entropy is the highest, can hide more watermark bits and is selected to carry watermark. To obtain directional components on horizontal cone, a llevel extended discrete shearlet decomposition is applied to test image with size , which is associated with shear matrix and direction numbers . Then, to find out the directional component whose information entropy is the highest, information entropy of each directional component on horizontal cone is computed by (9).
Here we suppose the th directional component’s information entropy is the highest. Then the binary sequence is embedded into the largest coefficients of the th directional component bywhere , is used to control embedding strength, and is the threshold.
Finally, the watermarked image is obtained after 1level extended inverse discrete shearlet transform by using these modified shearlet coefficients.
3.3. Watermark Extraction
The watermark can be extracted correctly from the watermarked image without original test image. The extracting process is as follows, which is mirroring embedding process.
To extract embedded watermark bits from the watermarked image, 1level extended discrete shearlet decomposition is firstly applied to the watermarked image, where the shear matrix is selected as and the number of direction is ; then shearlet coefficients of every directional component on horizontal cone are obtained.
Then, the watermarked coefficients are found out according to watermark size and the location key of embedding watermark. The extracted watermark bits are recovered into a sequence bywhere , , , and is the threshold.
Finally, the extracted watermark is recovered by reshaping .
4. Experimental Results
4.1. Invisibility of Watermark
Invisibility is an evaluative measure of perceptual quality of the watermarked image. In a satisfactory image watermark algorithm, watermark should not cause much degradation of perceptual quality of the watermarked image. In the proposed algorithm, a watermark image is embedded into different test images (“Lena,” “Boat,” and “Barbara” with ) to test invisibility, where the number of directions is and the embedding strength is 0.1.
As shown in Table 1, there are not much visual differences between original test images and their corresponding watermarked images. In addition, the PSNR value of each watermarked image is above 53 dB, and NC value of corresponding extracted watermark image is above 0.98. The extracted watermarks are all easily distinguishable. Whether from a subjective view or from the objective view, the proposed algorithm has preferable invisibility.

Furthermore, by analyzing the absolute difference between the test image and the watermarked image, it is not difficult to find that the watermark tends to spread on these regions with strong edges and extreme brightness. In particular, it is noted that the watermarked image is magnified a hundred times.
4.2. Robustness of the Proposed Algorithm
Robustness is a default measure which is used to evaluate the watermarked image’s resistance against attacks. In a satisfactory algorithm for image watermark, the watermark would not be easily removed from the watermarked image after common and deliberate attacks. In this paper, we simulate the real situation of disguised attacks and check out robustness for the proposed algorithm, where the number of orientations is 9 and the embedding strength is 0.1.
4.2.1. Robustness to Scaling Attack
For the “Lena” image and “Barbara” image, when scaling parameters are 1/2, 3/4, 3/2, 2, and 4, the extracted watermarks are still distinguishable and NC values are all above 0.90. For the “Boat” image, the extracted watermarks are still distinguishable and NC values are all above 0.90 when scaling parameters are 3/4, 3/2, 2, and 4. With the increase of scaling attack parameter, PSNR values are increasing; MSE values are decreasing. Therefore, the proposed algorithm can resist scaling attack.
4.2.2. Robustness to JPEG Compression Attack
JPEG lossy compression is the most common attack during the transformation of the watermarked image. In order to evaluate robustness of proposed algorithm, every watermarked image is compressed using JPEG with different quality factors. From extracted watermarks in Table 3, it is not hard for anyone to observe that extracted watermark is still distinguishable until the quality factor is set to 30. In addition, for all test images, PSNR values are all above 30 dB; NC values are all above 0.8. Therefore, the proposed algorithm is very robust against JPEG compression.
4.2.3. Robustness to Filtering Attack
As shown in Table 4, whether Gaussian lowpass filtering, average filtering, or median filtering, the extracted watermarks are still distinguishable for the “Lena” image and the “Barbara” image. In addition, NC values are all above 0.855 for the “Lena” image, 0.77 for the “Boat” image, and 0.86 for the “Barbara” image. However, the extracted watermarks are imperfect for the “Boat” image.
4.2.4. Robustness under Various Noises Attack
The watermarked image is easily disturbed by various noises. In the test, we add “Salt and Pepper” (Table 5) noise and “Gaussian” white noise (Table 6) to the watermarked images.
For the watermarked image with “Salt and Pepper” noise, with the increase of noise attack strength, more and more noise points appear in the extracted watermark, but the extracted watermark is still distinguishable and complete until the attack strength is set to 0.05. When the attack strength is set to 0.05, the extracted watermark is not perfect but distinguishable; the NC value is 0.8262 for the Lena image, 0.7866 for the Boat image, and 0.7861 for the Barbara image, respectively.
Similarly, for the watermarked image under the attack of “Gaussian” white noise, when the noise mean is set to zero, noise points increasingly appear in extracted watermark with the increase of noise variance. The extracted watermarks are still distinguishable and perfect for all test images until the noise variance is set to 0.0005. When the noise variance is set to 0.0005 and 0.001, the extracted watermark is incomplete but distinguishable for all test images; the NC values are all above 0.8. Therefore, the proposed algorithm can resist noise attack to some extent.
4.2.5. The Comparison Experiments of Related Algorithm
To verify the effectiveness of the proposed algorithm, the proposed algorithm is compared with related algorithms based on DWT and DCT. DCTbased watermark algorithm is performed by embedding watermark in the largest coefficients of the firstlevel DCT subband. Watermarking based on DWT is performed by embedding watermark in the largest coefficients of the firstlevel subband, where a onelevel DWT decomposition is performed using “Haar” wavelet.
To be fair, both the watermark image and the test images are the same to all algorithms. Table 7 shows the results of comparison on “Lena” image. As shown in Table 7, the proposed algorithm outperforms other algorithms in invisibility. Meanwhile, in robustness against these attacks, such as scaling attack and cropping attack, the proposed algorithm gets advantage over other algorithms.
In Table 2, although the NC values of proposed algorithm are not higher than other algorithms except no attack and scaling attack, it is noted that PSNR values of proposed algorithm are higher than other algorithms except cropping middle (50%). In other words, watermark which is embedded in host image cause much degradation of perceptual quality of the watermarked image, yet the quality of the watermarked image for proposed algorithm is better.






5. Conclusion
This paper tries to implement an image watermark by taking advantage of the unique advantages of extended discrete shearlet and selecting watermark embedding location based on the highest information entropy. The unique advantages of extended discrete shearlet made the watermark system more invisible. Experimental results demonstrate that visual performance of the watermarked images is preferable to that of related algorithms based on DWT and DCT; meanwhile, the proposed watermark algorithm can resist common attack. Therefore, the extended discrete shearlet can be considered as an appropriate transform for digital watermark.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by National Natural Science Foundation of China (no. 61379010) and Natural Science Basic Research Plan in Shaanxi Province of China (nos. 2012JQ1012 and 2015JM6293).
References
 G. Y. Liu, H. J. Liu, and A. Kadir, “Waveletbased color pathological image watermark through dynamically adjusting the embedding intensity,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 406349, 10 pages, 2012. View at: Publisher Site  Google Scholar
 G. Bhatnagar, Q. M. J. Wu, and P. K. Atrey, “Secure randomized image watermarking based on singular value decomposition,” ACM Transactions on Multimedia Computing, Communications and Applications, vol. 10, no. 1, article 4, 2013. View at: Publisher Site  Google Scholar
 Y. X. Ma, J. C. Xie, and J. G. Luo, “Image enhancement based on nonsubsampled contourlet transform,” in Proceedings of the 5th International Conference on Information Assurance and Security (IAS '09), pp. 31–34, September 2009. View at: Publisher Site  Google Scholar
 H. Sadreazami and M. Amini, “A robust spread spectrum based image watermarking in ridgelet domain,” AEU—International Journal of Electronics and Communications, vol. 66, no. 5, pp. 364–371, 2012. View at: Publisher Site  Google Scholar
 P. Mangaiyarkarasi and S. Arulselvi, “A new digital image watermarking based on finite ridgelet transform and extraction using lCA,” in Proceedings of the International Conference on Emerging Trends in Electrical and Computer Technology, pp. 837–841, March 2011. View at: Publisher Site  Google Scholar
 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  MathSciNet
 V. M. Patel, G. R. Easley, and J. Healy, “Shearletbased deconvolution,” IEEE Transactions on Image Processing, vol. 18, no. 12, pp. 2673–2685, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. Ranjbar, F. Zargari, and M. Ghanbari, “A highly robust twostage Contourletbased digital image watermarking method,” Signal Processing: Image Communication, vol. 28, no. 10, pp. 1526–1536, 2013. View at: Publisher Site  Google Scholar
 W.Q. Lim, “The discrete shearlet transform: a new directional transform and compactly supported shearlet frames,” IEEE Transactions on Image Processing, vol. 19, no. 5, pp. 1166–1180, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 R. Gomathi and A. V. A. Kumar, “An efficient GEM model for image inpainting using a new directional sparse representation: Discrete Shearlet Transform,” in Proceedings of the IEEE International Conference on Computational Intelligence and Computing Research (ICCIC '10), pp. 1–4, Coimbatore, India, December 2010. View at: Publisher Site  Google Scholar
 X.Y. Wang, Y.C. Liu, and H.Y. Yang, “An efficient remote sensing image denoising method in extended discrete shearlet domain,” Journal of Mathematical Imaging and Vision, vol. 49, no. 2, pp. 434–453, 2014. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. E. Shanonn and W. Weaver, The Mathematical Theory of Communication, University of Illinois Press, 1946.
 C. E. Shannon and E. Claude, “A mathematical theory of communication,” Bell System Technical Journal, vol. 27, no. 3, pp. 379–423, 623–656, 1948. View at: Google Scholar
 K. Guo and D. Labate, “Optimally sparse multidimensional representation using shearlets,” SIAM Journal on Mathematical Analysis, vol. 39, no. 1, pp. 298–318, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2015 Zhao Jian 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.