CN103778592A - 一种基于进制余码的双彩色图像盲水印方法 - Google Patents
一种基于进制余码的双彩色图像盲水印方法 Download PDFInfo
- Publication number
- CN103778592A CN103778592A CN201410014632.3A CN201410014632A CN103778592A CN 103778592 A CN103778592 A CN 103778592A CN 201410014632 A CN201410014632 A CN 201410014632A CN 103778592 A CN103778592 A CN 103778592A
- Authority
- CN
- China
- Prior art keywords
- watermark
- coefficient
- image
- color
- ary
- 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
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
Abstract
本发明公开了一种基于进制余码的双彩色图像盲水印方法,将大容量彩色图像作为数字水印嵌入到彩色宿主图像,满足了基于彩色图像标识的版权保护需要。本发明提出了进制余码的定义,通过调整宿主图像像素数据的进制余码来嵌入非二进制水印信息;提取水印时,利用含水印图像的进制余码即可提取水印,不要原始水印或原始宿主图像的帮助,达到了盲检测的目的。该方法简单快捷,具有较高的水印嵌入量和水印不可见性,本发明适用于彩色数字图像作为数字水印的版权保护。
Description
技术领域
本发明属于信息安全技术领域,涉及大容量彩色数字图像作为数字水印的版权保护。
背景技术
彩色数字图像作为网络信息的主要载体,很容易发生伪造、滥用、侵权等不良行为,面对多媒体版权保护的现实状况,需要加强对嵌入水印是高信息量彩色图像的数字水印技术研究。如今的图像数字水印方法多数是将二值或灰度图像作为数字水印,而将彩色数字图像作为数字水印的比较少,其主要原因是将彩色数字图像作为水印时,其含有的信息量是相同尺寸灰度图像的3倍,是二值图像的24倍,故增加了水印嵌入的难度,降低了水印不可见性,因此,如何将彩色数字图像作为数字水印成为亟待解决的问题之一。
在近几年的数字水印技术研究中,虽然盲检测数字水印方法逐渐成为数字水印技术发展的主流,但是由于彩色数字图像所包含的版权保护信息非常大,用非盲水印方法可以比较方便地嵌入或提取水印,现有的很多彩色图像数字水印方法研究的是非盲水印,因此,如何改变现状实现彩色水印图像的盲提取是当前图像数字水印技术的难点之一。
发明内容
本发明的目的是提供一种基于进制余码的双彩色图像盲水印方法,包含具体的水印嵌入过程和提取过程,其水印嵌入的具体过程描述如下:
第一步:彩色水印图像的预处理:将大小为N×N的24位原始彩色水印图像W的R、G、B三个水印分量分别降维处理得到一维数据,同时将其每一数值型数据换为长度为3的等长字符型数据,并进行字符连接得到字符型水印信息Wr、Wg、Wb;
第二步:将彩色宿主图像I的R、G、B三个分量依次进行一级整数小波变换,获得各分量的小波低频系数Hr、Hg、Hb;
第四步:利用公式(1),计算该系数单元的初始进制余码s0;
其中,mod(.)表示求余函数,r表示数值a i 所具有的进制数;
第五步:从字符型水印信息中截取当前水印信息w,并利用公式(2)计算w与初始进制余码s0的差值余数e;
第六步:按照下面的规则来嵌入水印w,并按公式(3)计算新的进制余码s1,直到新的进制余码s1与嵌入水印w一致,完成水印w嵌入;
规则一:如果e=0,则该系数单元的值不作任何修改;
规则二:如果1≤e≤5,则该系数单元中a e 值加1,当a e >255时,则令a e =255;
规则三:如果e>5,则该系数单元中a 10-e 值减1,当a 10-e <0时,则令a 10-e =1;
第七步:重复执行第四步到第六步,直到所有水印信息被嵌入到宿主图像;
第八步:对修改后的整型小波系数进行逆小波变换得到含水印的R、G、B分量,然后合并各分量图像得到含水印的图像。
本发明所述水印提取的具体过程如下:
第一步:将含水印的彩色图像按三原色分成R、G、B三个分量,并分别执行一级整数小波变换,获得低频系数Tr、Tg、Tb;
第二步:将每一低频系数中每5个系数值组成一个系数单元a i * ;
第三步:依据公式(4)从系数单元a i * 中提取水印信息w*;
第四步:重复执行第二步、第三步,直到提取所有的水印信息;
第五步:按先后顺序将提取的字符型水印信息w*进行3位组合,并进行数据类型转换,最终提取R、G、B分量中含有的水印;
第六步:组合分量水印形成最终水印W*。
该方法简单快捷,嵌入的高信息量水印具有较好的水印不可见性,本发明达到如下效果:将彩色数字图像作为数字水印并“毫无觉察”地嵌入到宿主图像,提取水印时不需要原始宿主图像或原始水印图像的帮助,能从各种受攻击图像中快速提取所嵌入的水印,该发明适用于彩色数字图像作为数字水印的版权保护。
附图说明
图1(a)、图1(b)、图1(c)、图1(d)是四幅原始彩色宿主图像。
图2是一幅彩色水印图像。
图3(a)、图3(b)、图3(c)、图3(d)是将图2所示的水印依次嵌入到宿主图像图1(a)、图1(b)、图1(c)、图1(d)后所得到的含水印图像,其结构相似度SSIM值依次是0.99891、0.99962、0.99880、0.99874,其峰值信噪比PSNR值依次是55.8543dB、55.8482dB、55.8415dB、55.8774dB。
图4(a)、图4(b)、图4(c)、图4(d)是依次从图3(a)、图3(b)、图3(c)、图3(d)中提取的水印,其归一化互相关系数NC值分别是1.00000、1.00000、0.99998、1.00000。
图5(a)、图5(b)、图5(c)、图5(d)、图5(e)、图5(f)是将图3(a)所示的含水印图像依次进行JPEG压缩、加噪、中值滤波、巴特沃斯低通滤波、Mosaic 拼贴、旋转等攻击后所提取的水印,其归一化互相关系数NC值分别是0.89854、0.98491、0.90453、0.89365、0.89473、0.92016。
具体实施方式
本发明的目的是提供一种基于进制余码的双彩色图像盲水印方法,包含具体的水印嵌入过程和提取过程,通过实例描述其水印嵌入过程如下:
第一步:彩色水印图像的预处理:将大小为64×64的24位原始彩色水印图像W的R、G、B三个水印分量分别降维处理得到一维数据,同时将其每一数值型数据转换为长度为3的等长字符型数据,并进行字符连接得到字符型水印信息Wr、Wg、Wb;例如,将三个像素值206,66,5分别转换为“206”,“066”,“005”,依次连接后得到字符型水印信息“206066005”。
第二步:将512×512的彩色宿主图像I的R、G、B三个分量依次进行一级整数小波变换,获得各分量的小波低频系数Hr、Hg、Hb;
第四步:设r=10,利用公式(1),计算该系数单元的初始进制余码s0;
第五步:从字符型水印信息中截取当前水印信息w=“2”,并利用公式(2)计算w与初始进制余码s0的差值余数e;
第六步:因此时1≤e≤5,则按照规则二来嵌入水印w,将系数单元中的第4个元素修改为155,并利用公式(3)计算新的进制余码s1,此时新的进制余码s1与嵌入水印w一致,完成了水印嵌入;
第七步:重复执行第四步到第六步,直到所有水印信息被嵌入到宿主图像;
第八步:对修改后的整型小波系数进行逆小波变换得到含水印的R、G、B分量,然后合并各分量图像得到含水印的图像。
通过实例描述其水印提取的具体过程如下:
第一步:将含水印的彩色图像按三原色分成R、G、B三个分量,并分别执行一级整数小波变换,获得低频系数Tr、Tg、Tb;
第二步:每一低频系数中每5个系数值组成一个系数单元a i * ,设当前系数单元是{156,155,157,155,149};
第三步:设r=10,利用公式(4),从系数单元a i * 中提取水印信息w*;
第四步:重复执行第二步、第三步,直到提取所有的水印信息;
第五步:按先后顺序将提取的字符型水印信息w*进行3位组合,并进行数据类型转换,最终提取R、G、B分量中含有的水印;
第六步:组合分量水印形成最终水印W*。
该方法简单快捷,具有较高的水印嵌入量和水印不可见性,适用于彩色图像作为数字水印的版权保护。
本发明有效性验证
为了证明本发明的有效性,选择如图1(a)、图1(b)、图1(c)、图1(d)所示的四幅大小为512×512的24位标准图像作为宿主图像,并将如图2所示的一幅大小为64×64的24位彩色图像作为数字水印进行验证。
图3(a)、图3(b)、图3(c)、图3(d)是将图2所示的水印依次嵌入到宿主图像图1(a)、图1(b)、图1(c)、图1(d)后所得到的含水印图像,其结构相似度SSIM值依次是0.99891、0.99962、0.99880、0.99874,其峰值信噪比PSNR值依次是55.8543dB、55.8482dB、55.8415dB、55.8774dB,可以看出所嵌入的水印具有良好的不可见性;
图4(a)、图4(b)、图4(c)、图4(d)是依次从图3(a)、图3(b)、图3(c)、图3(d)中提取的水印,其归一化互相关系数NC值分别是1.00000、1.00000、0.99998、1.00000。图5(a)、图5(b)、图5(c)、图5(d)、图5(e)、图5(f)是将图3(a)所示的含水印图像依次进行JPEG压缩、加噪、中值滤波、巴特沃斯低通滤波、Mosaic 拼贴、旋转等攻击后所提取的水印,其归一化互相关系数NC值分别是0.89854、0.98491、0.90453、0.89365、0.89473、0.92016,由此可见,所提取的数字水印图像具有良好的可鉴别性,说明该方法具有较强的鲁棒性,能够很好地提取所嵌入的彩色水印。
Claims (1)
1.一种基于进制余码的双彩色图像盲水印方法,其特征在于:是通过具体的水印嵌入过程和提取过程的,其水印嵌入的具体过程描述如下:
第一步:彩色水印图像的预处理:将大小为N×N的24位原始彩色水印图像W的R、G、B三个水印分量分别降维处理得到一维数据,同时将其每一数值型数据换为长度为3的等长字符型数据,并进行字符连接得到字符型水印信息Wr、Wg、Wb;
第二步:将彩色宿主图像I的R、G、B三个分量依次进行一级整数小波变换,获得各分量的小波低频系数Hr、Hg、Hb;
第四步:利用公式(1),计算该系数单元的初始进制余码s0;
其中,mod(.)表示求余函数,r表示数值a i 所具有的进制数;
第五步:从字符型水印信息中截取当前水印信息w,并利用公式(2)计算w与初始进制余码s0的差值余数e;
第六步:按照下面的规则来嵌入水印w,并按公式(3)计算新的进制余码s1,直到新的进制余码s1与嵌入水印w一致,完成水印w嵌入;
规则一:如果e=0,则该系数单元的值不作任何修改;
规则二:如果1≤e≤5,则该系数单元中a e 值加1,当a e >255时,则令a e =255;
规则三:如果e>5,则该系数单元中a 10-e 值减1,当a 10-e <0时,则令a 10-e =1;
第七步:重复执行第四步到第六步,直到所有水印信息被嵌入到宿主图像;
第八步:对修改后的整型小波系数进行逆小波变换得到含水印的R、G、B分量,然后合并各分量图像得到含水印的图像;
本发明所述水印提取的具体过程如下:
第一步:将含水印的彩色图像按三原色分成R、G、B三个分量,并分别执行一级整数小波变换,获得低频系数Tr、Tg、Tb;
第二步:将每一低频系数中每5个系数值组成一个系数单元a i * ;
第三步:依据公式(4)从系数单元a i * 中提取水印信息w*;
第四步:重复执行第二步、第三步,直到提取所有的水印信息;
第五步:按先后顺序将提取的字符型水印信息w*进行3位组合,并进行数据类型转换,最终提取R、G、B分量中含有的水印;
第六步:组合分量水印形成最终水印W*。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410014632.3A CN103778592B (zh) | 2014-01-14 | 2014-01-14 | 一种基于进制余码的双彩色图像盲水印方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410014632.3A CN103778592B (zh) | 2014-01-14 | 2014-01-14 | 一种基于进制余码的双彩色图像盲水印方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103778592A true CN103778592A (zh) | 2014-05-07 |
CN103778592B CN103778592B (zh) | 2016-07-06 |
Family
ID=50570791
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410014632.3A Expired - Fee Related CN103778592B (zh) | 2014-01-14 | 2014-01-14 | 一种基于进制余码的双彩色图像盲水印方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103778592B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106157233A (zh) * | 2016-07-11 | 2016-11-23 | 鲁东大学 | 一种基于矩阵Schur分解的高性能数字水印方法 |
CN106169171A (zh) * | 2016-07-15 | 2016-11-30 | 鲁东大学 | 基于Hessenberg分解的高性能数字水印方法 |
CN107135395A (zh) * | 2017-05-11 | 2017-09-05 | 宁波大学 | 一种在数字彩色图像中实现秘密信息密写及提取的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101835049A (zh) * | 2010-05-24 | 2010-09-15 | 上海师范大学 | 一种jpeg图像自嵌入数字水印的生成及认证方法 |
US20110228972A1 (en) * | 2005-08-04 | 2011-09-22 | Nippon Telegraph And Telephone Corporation | Digital watermark padding method, digital watermark padding device, digital watermark detecting method, digital watermark detecting device, and program |
-
2014
- 2014-01-14 CN CN201410014632.3A patent/CN103778592B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110228972A1 (en) * | 2005-08-04 | 2011-09-22 | Nippon Telegraph And Telephone Corporation | Digital watermark padding method, digital watermark padding device, digital watermark detecting method, digital watermark detecting device, and program |
CN101835049A (zh) * | 2010-05-24 | 2010-09-15 | 上海师范大学 | 一种jpeg图像自嵌入数字水印的生成及认证方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106157233A (zh) * | 2016-07-11 | 2016-11-23 | 鲁东大学 | 一种基于矩阵Schur分解的高性能数字水印方法 |
CN106157233B (zh) * | 2016-07-11 | 2019-11-05 | 鲁东大学 | 一种基于矩阵Schur分解的高性能数字水印方法 |
CN106169171A (zh) * | 2016-07-15 | 2016-11-30 | 鲁东大学 | 基于Hessenberg分解的高性能数字水印方法 |
CN106169171B (zh) * | 2016-07-15 | 2019-08-13 | 鲁东大学 | 基于Hessenberg分解的高性能数字水印方法 |
CN107135395A (zh) * | 2017-05-11 | 2017-09-05 | 宁波大学 | 一种在数字彩色图像中实现秘密信息密写及提取的方法 |
CN107135395B (zh) * | 2017-05-11 | 2019-08-20 | 宁波大学 | 一种在数字彩色图像中实现秘密信息密写及提取的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103778592B (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105761199B (zh) | 一种基于组合域的鲁棒彩色图像水印方法 | |
CN108648134B (zh) | 一种融合离散小波变换的空域彩色数字图像盲水印方法 | |
CN106817586B (zh) | 一种h.264/avc视频可逆数据隐藏方法 | |
CN104168484A (zh) | 一种基于视觉注意机制的视频水印方法 | |
CN109102454B (zh) | 一种融合快速傅里叶变换的彩色qr码数字盲水印方法 | |
CN109829845A (zh) | 基于矩阵Schur分解的可变步长彩色图像盲水印方法 | |
CN104050622B (zh) | 基于三进制编码的彩色图像盲水印方法 | |
CN102663669B (zh) | 一种基于Schur分解的双彩色数字图像盲水印方法 | |
CN106204410B (zh) | 一种新型的基于矩阵Schur分解的数字水印方法 | |
CN104899823B (zh) | 基于Hessenberg分解的双彩色图像盲水印方法 | |
CN103810668A (zh) | 一种基于奇异值分解的双彩色图像盲水印方法 | |
CN103617591A (zh) | 一种基于qr分解与补偿的双彩色图像盲水印方法 | |
CN103778592B (zh) | 一种基于进制余码的双彩色图像盲水印方法 | |
CN105023236B (zh) | 一种基于矩阵lu分解的双彩色图像盲水印方法 | |
CN104123693A (zh) | 一种多功能立体图像数字水印方法 | |
CN106991636B (zh) | 一种融合近似Schur分解的空域彩色图像盲水印方法 | |
CN116579908B (zh) | 一种在图像中植入加密隐藏信息的方法及装置 | |
CN106157233B (zh) | 一种基于矩阵Schur分解的高性能数字水印方法 | |
CN105141962B (zh) | 一种面向h.264/avc视频码流的大容量数据隐藏方法 | |
CN106169171B (zh) | 基于Hessenberg分解的高性能数字水印方法 | |
CN110415155A (zh) | 一种融合哈尔变换的空域彩色图像盲水印方法 | |
CN102314668B (zh) | 一种用于提高嵌入水印后图像质量的差值扩展嵌入数字水印改进方法 | |
CN109829846A (zh) | 一种基于二维离散余弦变换的数字图像盲水印方法 | |
CN105208462B (zh) | 视频数字水印处理方法及执行该方法的视频数字水印系统 | |
CN110415154B (zh) | 一种基于哈尔变换的四进制彩色数字图像盲水印方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160706 Termination date: 20170114 |