CN102279969A - 基于contourlet和商空间的抗打印扫描数字水印方法 - Google Patents
基于contourlet和商空间的抗打印扫描数字水印方法 Download PDFInfo
- Publication number
- CN102279969A CN102279969A CN2011102354331A CN201110235433A CN102279969A CN 102279969 A CN102279969 A CN 102279969A CN 2011102354331 A CN2011102354331 A CN 2011102354331A CN 201110235433 A CN201110235433 A CN 201110235433A CN 102279969 A CN102279969 A CN 102279969A
- Authority
- CN
- China
- Prior art keywords
- watermark
- image
- coefficient
- contourlet
- conversion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Editing Of Facsimile Originals (AREA)
- Image Processing (AREA)
Abstract
一种基于contourlet和商空间的抗打印扫描数字水印方法,属于数字水印技术领域,其水印嵌入步骤如下:一、读入水印图像,采用Arnold变换对二值水印图像进行置乱;二、将置乱后的水印通过运算将其转换为二值序列;三、对原始载体图形进行三层Contourlet分解,得到等子带;四、找出三层低频区域系数,确定水印的嵌入位置、实现水印信息的嵌入;五、对嵌入水印后的图像进行逆Contourlet变换,重构出含有水印信息的图像。本方法能有效地平衡水印的嵌入容量和不可见性,根据PSNR(峰值信噪比)和SSIM(结构相似度索引)的值来判断水印在打印扫描过程中失真不大,算法的鲁棒性得到了一定的保证。具有兼顾数字水印的鲁棒性和不可见性两方面的优点。该方法可以有效抵抗打印扫描过程中的各种攻击。
Description
技术领域
本发明涉及一种基于contourlet和商空间的抗打印扫描数字水印方法,属于数字水印技术领域。
背景技术
由于图像在打印和扫描处理的过程所受的攻击太强,针对该现象,目前使用的抗打印扫描数字水印的方法主要集中在基于DCT变换和DFT变换。如有的做DCT变换后,再选取一些符合特定条件的变换系数,对这些系数进行线性运算变换,从而实现了水印信息的嵌入。水印是通过多次反复操作实现嵌入的,嵌入后的水印还能用来确定图像所经受的扫描次数,这种水印算法提高了水印的鲁棒性,但稳健性不是太好而且水印的容量小;还有的根据傅立叶—梅林(FM)变换,并利用了DFT能抵抗各种几何变形,对平移操作具有不变性等特性,提出了一种健壮性不太理想的水印算法,因为图像经过此种DFT后系数会发生较大的变化,使得提取方法无通用性。以上方法都存在一定的简单化问题。
发明内容
为了解决上述存在的问题,本发明提供一种能够利用contourlet变换能有效有效地捕捉图像边缘轮廓等性,有效地提高了水印算法的鲁棒性的抗打印扫描数字水印方法。
本发明的目的是通过下述技术方案实现的:一种基于contourlet和商空间的抗打印扫描数字水印方法,包括水印嵌入和水印提取过程,其中,水印嵌入步骤如下:
Step1:读入水印图像,若水印图像是灰度图像则将其转化为二值图像,然后采用Arnold变换对二值水印图像进行置乱;
Step2:将置乱后的水印通过运算将其转换为二值序列:
Step3:对图像I进行三层Contourlet分解,得到等子带;
Step4:找出由第三层低频区域系数所构成的商空间中的粗粒度X1;根据X1中系数的位置分别选取第二级子带和第一级子带中对应的细粒度,找出水印的嵌入位置;记录这三层商空间中符合条件的重要系数的位置并统计已确定的重要系数的个数k,并且记住这些系数的位置;
Step5:根据所要嵌入的水印信息的值来修改Step4中涉及系数的值,具体修改方法如下:
其中a为一个参数,它可以控制水印信息的嵌入强度,其值是通过具体实验结果进行调整;
Step6:按照同样方式来处理其他两层中嵌入位置的低频系数值,从而实现水印信息的嵌入;
Step7:对嵌入水印后的图像进行逆Contourlet变换,重构出含有水印信息的图像。
基于contourlet和商空间的抗打印扫描数字水印方法,其水印提取过程如下:
Step1:对打印扫描后的图像做三级Contourlet变换,得到三层Contourlet分解系数;
Step2:解开密钥,获取水印的嵌入位置和嵌入个数,从而提取Contourlet变换后的图像的三层商空间;
Step4:按照同样的方法分别处理其他两层的嵌入位置的Contourlet变换系数,由于商空间为三层结构,因此我们就能够提取出三组水印信息 S1、S2、S3,对M进行二值化处理就获得了经过Arnold变换后的水印信息S’;
Step5:对水印信息进行逆Arnold变换就获得了提取出的水印图像S。
本发明的有益效果:本发明采用上述方法与现有技术相比,能有效地平衡水印的嵌入容量和不可见性,根据PSNR(峰值信噪比)和SSIM(结构相似度索引)的值来判断水印在打印扫描过程中失真不大,算法的鲁棒性得到了一定的保证。具有兼顾数字水印的鲁棒性和不可见性两方面的优点。该方法可以有效抵抗打印扫描过程中的各种攻击:半色调、D/A转换、A/D转换、缩放、旋转、剪切、扭曲、各种随机噪声等。
附图说明
图1是contourlet变换的示意图。
图2是Lena图像分解的示例。
图3是商空间理论中的关系图。
图4粒度层次结构关系图。
图5水印嵌入过程示意图。
图6 水印提取过程示意图。
图7粒度层次分解矩阵示意图。
图8实验鲁棒性效果图。
具体实施方式
本发明涉及的基于contourlet和商空间的抗打印扫描数字水印方法,其中Contourlet变换是第二代小波,它不仅继承了小波变换的多分辨率时频分析特征,而且拥有良好的各向异性特征,它能用比小波变换更少的系数来表达光滑的曲线,并且克服了小波变换中方向性的不足(只有水平、竖直和对角线方向)。这种变换可以看成是由拉普拉斯金字塔(1aplacian pyramid,LP)分解和方向滤波器组(pyranmidal directional filter bank,PDFB)两部分构成的双重方向滤波器组(参见文献M.N. Do and M. Vetterli. Framing Pyramids. IEEE Transaction. Signal Processing,no.9,pp.2329-2342,Sep.2003.)。如图1和图2分别是contourlet变换的示意图和Lena图像分解的示例。
它提供了丰富的方向和形状集合,故其在捕获光滑轮廓和几何结构上更为有效。利用这种变换,根据Contourlet系数之间的父子关系,构造出Contourlet系数重要子树关系。
本发明方法中还引进了商空间理论和粒度计算的概念,在这个理论中,将Contourlet系数划分为不同的粒度。最后根据粒度计算的保真特性和保假特性,将水印信息快速地嵌入到这些粒度中。
我们用一个三元组(X,F,T)来表示商空间,商空间的定义为(参见刘仁金,黄贤武.图像分割的商空间粒度原理[J].计算机学报,2005,Vol.28 No.10):
商空间理论与模糊集和粗糙集的关系如图3所示,它们同属于近几年新提出来的粒度计算范畴。引入商空间理论可以解决如何尽可能找到所感兴趣的区域而又不会因传统计算技术的范围过于精确而有所遗漏的问题。
因而,引入了商空间理论和粒度计算概念,可以快速地找到嵌入位置,有效地降低了算法的时间复杂度。
其中粒度计算的保假原理:假如存在一个问题在粗粒度空间中找不到解,那么这个问题在比它细的粒度空间中肯定也找不到解。
粒度计算的保真原理:假如存在一个问题在两个粗粒度空间中能找到问题的解,那么这个问题在这两个粗粒度空间的合成空间上肯定也能找到解,并且合成空间上的求解复杂度是这两个粗粒度空间上的求解复杂度之和。
在商空间理论中,应该使用分层方法对各层变换系数进行粗分割,再考虑向更高层次系数继续划分,例如在三级Contourlet分解后的低频系数中,在第三层系数中如果一个系数可以划归到粗粒度里面,根据粒度计算的保真原理,那么和这个系数对应的第二层子系数就可以划归到第二层粒度中,同理可以迅速找到第一层细粒度中的元素。这种分层思路有助于实现多层嵌入水印信息,可以全面分析问题,减少运算步骤,降低算法的计算复杂度。
Step1:读入水印图像,(若水印图像是灰度图像则要将其转化为二值图像)然后采用Arnold变换对二值水印图像进行置乱,这样可以提高水印的保密性,并且能打乱水印信号的自身相关性;具体到N的数字图像的Arnold变换为:
(1)
Step2:将置乱后的水印通过相应的运算将其转换为二值序列:
Step4:找出由第三层低频区域系数所构成的商空间中的粗粒度X1,根据X1中系数的位置分别选取第二级子带和第一级子带中对应的细粒度,找出水印的嵌入位置。记录这三层商空间中符合条件的重要系数的位置并统计已确定的重要系数的个数k,并且记住这些系数的位置,在水印的提取过程中还要用到;
Step5:根据所要嵌入的水印信息的值来修改上面选取的系数的值,具体修改方法如下:
其中a为一个参数,它可以控制水印信息的嵌入强度,其值是通过具体实验结果进行调整;
Step6:按照类似方式来处理其他两层中嵌入位置的低频系数值,从而实现水印信息的嵌入;
Step7:对嵌入水印后的图像进行逆Contourlet变换,重构出含有水印信息的图像。
即先对载体图像Lena图像进行三级Contourlet分解(如图2所示),接着要实现在三层分解系数中嵌入水印信息,对于这个复杂的嵌入工作一般先将这项工作放在粗粒度空间中(第三层粗粒度空间系数,如图7所示)图7是粒度层次分解矩阵图,其中的“第三层粒度(粗粒度)”按下述的公式(2)和(3)来选取确定,由于进行的是三层Contourlet分解,直接在对应Contourlet分解的下一层(即第二层粒度(细粒度))的4个像素即为对应的第二层细粒度空间系数,同理可得最后一层的细粒度空间系数。
进行初步宏观上的分析,找出水印在粗粒度空间中的嵌入位置,对于不能作为嵌入位置的系数再结合上述的“保假原理”推出水印在细粒度空间中的对应位置也不能嵌入水印。这样在计算代价很小的情况下,我们依然能够实现水印信息的多层嵌入,大大提高了水印的嵌入容量,从而提高了算法的性能。二者的结合实现了在contourlet变换各层系数中嵌入水印,不仅增加了水印信息的嵌入量,提高了水印的鲁棒性,还降低了算法的时间复杂度,可用于印刷作品的版权保护。
上述Step3对图像进行三级Contourlet分解步骤:
a由于在Contourlet变换后的低频区域嵌入水印信息不易被人发现,水印的安全性比较好,因此在这里考察由它的第一级子带的低频区域(LL1)所构成的商空间中的系数,第一级中所有满足的 就构成一个集合—粗粒度X1;
(2)
这里的是一个阈值,其定义为:
其中,代表该子图像的能量,表示子带图像,和分别表示子带图像的宽度和高度。我们知道人眼对纹理细节变化不太敏感,能量的大小代表图像在该方向上的纹理细节的丰富程度,因为能量大的图像在该方向上的纹理细节就越丰富,所以在能量大的系数上嵌入水印比较安全。因此可以用平均能量作为阈值,来找出该层上的重要系数。
c:根据细粒度X2中的元素的位置迅速找到其所对应的由第三级子带中低频系数所构成的商空间中的系数,使用第一步的方法找到满足条件的系数形成细粒度X3。经过上述三步的划分后,就可以得到粒度的层次结构图,如图4所示。
水印的提取是水印嵌入的一个逆过程,在水印的检测中无需原始图像,只需记录重要系数位置的密钥,因此此算法是盲水印算法。
具体过程如下:
Step1:对打印扫描后的图像做三级Contourlet变换,得到三层Contourlet分解系数;
Step2:解开密钥,获取水印的嵌入位置和嵌入个数,从而提取Contourlet变换后的图像的三层商空间;
Step4:按照同样的方法分别处理其他两层的嵌入位置的Contourlet变换系数,由于商空间为三层结构,因此我们就能够提取出三组水印信息 S1、S2、S3,对M进行二值化处理就获得了经过Arnold变换后的水印信息S’;
Step5:对水印信息进行逆Arnold变换就获得了提取出的水印图像S。,其提取过程示意图如图6所示。
实验结果的鲁棒性效果,如图8所示。将不同的载体图像嵌入不同的水印后,分别经过HP4VC激光打印机打印,扫描仪为紫光B6210(光学分辨率为400dpi)的扫描作为攻击,所提取的水印表明该方法的鲁棒性较好。其主观效果如图8所示,而其客观效果可以由PSNR和SSIM作为参考。峰值信噪比(PNSR,peak signal to noise ratio),单位是dB,它的值越大,那么图像的重建质量就越好;另一个是SSIM(结构相似度索引),在文献WANG Z,BOVIK A C,SHEIKH H R,SIMONCELLI E P.Image quality assessment:from error visibility to structural similarity[J].IEEE Transactions on Image Processing,2004,13(4):600-612. 中详细介绍了这个指标的相关内容。它考虑了人类视觉模型(HVS,human visual system),从而能够更准确地衡量原始水印和提取水印两幅图像的相似程度,其值越大,两幅图像越相似,目前大多算法常采用这两个指标。
Claims (2)
1.一种基于contourlet和商空间的抗打印扫描数字水印方法,包括水印嵌入和水印提取过程,其中,水印嵌入步骤如下:
Step1:读入水印图像,若水印图像是灰度图像则将其转化为二值图像,然后采用Arnold变换对二值水印图像进行置乱;
Step2:将置乱后的水印通过运算将其转换为二值序列:
Step4:找出由第三层低频区域系数所构成的商空间中的粗粒度X1;根据X1中系数的位置分别选取第二级子带和第一级子带中对应的细粒度,找出水印的嵌入位置;记录这三层商空间中符合条件的重要系数的位置并统计已确定的重要系数的个数k,并且记住这些系数的位置;
Step5:根据所要嵌入的水印信息的值来修改Step4中涉及系数的值,具体修改方法如下:
其中a为一个参数,它可以控制水印信息的嵌入强度,其值是通过具体实验结果进行调整;
Step6:按照同样方式来处理其他两层中嵌入位置的低频系数值,从而实现水印信息的嵌入;
Step7:对嵌入水印后的图像进行逆Contourlet变换,重构出含有水印信息的图像。
2.根据权利要求1所述的基于contourlet和商空间的抗打印扫描数字水印方法,水印提取过程如下:
Step1:对打印扫描后的图像做三级Contourlet变换,得到三层Contourlet分解系数;
Step2:解开密钥,获取水印的嵌入位置和嵌入个数,从而提取Contourlet变换后的图像的三层商空间;
Step4:按照同样的方法分别处理其他两层的嵌入位置的Contourlet变换系数,由于商空间为三层结构,因此我们就能够提取出三组水印信息 S1、S2、S3,对M进行二值化处理就获得了经过Arnold变换后的水印信息S’;
Step5:对水印信息进行逆Arnold变换就获得了提取出的水印图像S。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110235433 CN102279969B (zh) | 2011-08-17 | 2011-08-17 | 基于contourlet和商空间的抗打印扫描数字水印方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110235433 CN102279969B (zh) | 2011-08-17 | 2011-08-17 | 基于contourlet和商空间的抗打印扫描数字水印方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102279969A true CN102279969A (zh) | 2011-12-14 |
CN102279969B CN102279969B (zh) | 2013-04-10 |
Family
ID=45105398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110235433 Expired - Fee Related CN102279969B (zh) | 2011-08-17 | 2011-08-17 | 基于contourlet和商空间的抗打印扫描数字水印方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102279969B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102843493A (zh) * | 2012-07-17 | 2012-12-26 | 辽宁大学 | 一种基于局部平衡三叉树理论的快速图像置乱方法 |
CN102938133A (zh) * | 2012-11-20 | 2013-02-20 | 海南大学 | 一种基于Arnold置乱变换和DWT-DFT的医学图像鲁棒水印方法 |
CN103426142A (zh) * | 2013-08-14 | 2013-12-04 | 中山大学深圳研究院 | 一种基于Arnold变换的Patchwork数字水印编解码方法 |
CN108492238A (zh) * | 2018-03-07 | 2018-09-04 | 小明太极(湖北)国漫文化有限公司 | 盲水印嵌入方法和提取方法 |
CN109102449A (zh) * | 2018-07-19 | 2018-12-28 | 郭伟 | 一种基于第二代Bandelet变换的数字水印的生成及提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010020270A1 (en) * | 1998-06-29 | 2001-09-06 | Minerva M. Yeung | Fragile watermarking for objects |
CN101308566A (zh) * | 2008-06-02 | 2008-11-19 | 西安电子科技大学 | 基于contourlet变换抗几何攻击数字图像水印方法 |
-
2011
- 2011-08-17 CN CN 201110235433 patent/CN102279969B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010020270A1 (en) * | 1998-06-29 | 2001-09-06 | Minerva M. Yeung | Fragile watermarking for objects |
CN101308566A (zh) * | 2008-06-02 | 2008-11-19 | 西安电子科技大学 | 基于contourlet变换抗几何攻击数字图像水印方法 |
Non-Patent Citations (3)
Title |
---|
张楠,范铁生: "基于小波变换的数字水印算法研究", 《科技创新导报》, no. 20085, 29 February 2008 (2008-02-29), pages 2 - 3 * |
王丹华,范铁生,魏晓蕊: "基于不定算法的频域音频数字水印的改进方案", 《自动化技术与应用》, vol. 26, no. 1, 31 January 2007 (2007-01-31), pages 99 - 102 * |
韦升文: "基于Contourlet和SVD的抗打印扫描零水印算法", 《福建电脑》, no. 20116, 30 June 2011 (2011-06-30), pages 7 - 10 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102843493A (zh) * | 2012-07-17 | 2012-12-26 | 辽宁大学 | 一种基于局部平衡三叉树理论的快速图像置乱方法 |
CN102843493B (zh) * | 2012-07-17 | 2015-02-25 | 辽宁大学 | 一种基于局部平衡三叉树理论的快速图像置乱方法 |
CN102938133A (zh) * | 2012-11-20 | 2013-02-20 | 海南大学 | 一种基于Arnold置乱变换和DWT-DFT的医学图像鲁棒水印方法 |
CN103426142A (zh) * | 2013-08-14 | 2013-12-04 | 中山大学深圳研究院 | 一种基于Arnold变换的Patchwork数字水印编解码方法 |
CN108492238A (zh) * | 2018-03-07 | 2018-09-04 | 小明太极(湖北)国漫文化有限公司 | 盲水印嵌入方法和提取方法 |
CN109102449A (zh) * | 2018-07-19 | 2018-12-28 | 郭伟 | 一种基于第二代Bandelet变换的数字水印的生成及提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102279969B (zh) | 2013-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Su et al. | A novel blind color image watermarking based on Contourlet transform and Hessenberg decomposition | |
Kumar et al. | DWT based color image watermarking using maximum entropy | |
Khan et al. | Reduced time complexity for detection of copy-move forgery using discrete wavelet transform | |
Keshavarzian et al. | ROI based robust and secure image watermarking using DWT and Arnold map | |
Tian et al. | Robust and blind watermarking algorithm based on DCT and SVD in the contourlet domain | |
Roy et al. | A hybrid domain color image watermarking based on DWT–SVD | |
CN102279969B (zh) | 基于contourlet和商空间的抗打印扫描数字水印方法 | |
Liu et al. | Robust blind image watermarking based on chaotic mixtures | |
Gul et al. | A novel triple recovery information embedding approach for self-embedded digital image watermarking | |
CN105512999A (zh) | 一种双变换的彩色图像全息水印方法 | |
Subhedar et al. | Image steganography using contourlet transform and matrix decomposition techniques | |
Niu et al. | A novel robust image watermarking in quaternion wavelet domain based on superpixel segmentation | |
Lu et al. | Wavelet-based CNN for robust and high-capacity image watermarking | |
Jamali et al. | Adaptive blind image watermarking using fuzzy inference system based on human visual perception | |
Luo et al. | Stereo image watermarking scheme for authentication with self-recovery capability using inter-view reference sharing | |
Saneie et al. | Introducing a new method of robust digital image watermarking against cropping and salt & pepper noise using Sudoku | |
Motwani et al. | Fuzzy Perceptual Watermarking For Ownership Verification. | |
Aberna et al. | Optimal Semi-Fragile Watermarking based on Maximum Entropy Random Walk and Swin Transformer for Tamper Localization | |
Rachmawanto et al. | A non-blind robust and impercept watermarking using discrete cosine transform and discrete wavelet transform | |
Jamali et al. | Adaptive image watermarking using human perception based fuzzy inference system | |
Hosam et al. | A hybrid ROI-embedding based watermarking technique using DWT and DCT transforms | |
Agarwal et al. | Enhanced data hiding method using DWT based on Saliency model | |
Xiaolei et al. | Robust Digital Watermarking Algorithm Based on DCT-SVD and QR Code | |
Umaamaheshvari et al. | A robust digital watermarking technique based on feature and transform method | |
Allwadhi et al. | A comprehensive survey of state-of-art techniques in digital watermarking |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130410 Termination date: 20150817 |
|
EXPY | Termination of patent right or utility model |