CN112488904A - 基于矩阵奇异值分解的四元数彩色数字图像盲水印方法 - Google Patents
基于矩阵奇异值分解的四元数彩色数字图像盲水印方法 Download PDFInfo
- Publication number
- CN112488904A CN112488904A CN202110059404.8A CN202110059404A CN112488904A CN 112488904 A CN112488904 A CN 112488904A CN 202110059404 A CN202110059404 A CN 202110059404A CN 112488904 A CN112488904 A CN 112488904A
- Authority
- CN
- China
- Prior art keywords
- watermark
- image
- quaternion
- blue
- block
- 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
- 239000011159 matrix material Substances 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 20
- 238000000605 extraction Methods 0.000 claims abstract description 5
- 239000003086 colorant Substances 0.000 claims description 9
- 238000013139 quantization Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 239000006002 Pepper Substances 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000006399 behavior Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/0021—Image watermarking
- G06T1/005—Robust watermarking, e.g. average attack or collusion attack resistant
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2201/00—General purpose image data processing
- G06T2201/005—Image watermarking
- G06T2201/0083—Image watermarking whereby only watermarked image required at decoder, e.g. source-based, blind, oblivious
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
Abstract
本发明利用频域数字水印算法鲁棒性强的优点,公开了一种基于矩阵奇异值分解(Singular Value Decomposition,SVD)的四元数彩色数字图像盲水印方法。本发明根据四元数表示的图像在奇异值分解后所得对角矩阵的特点,通过量化对角矩阵的最大特征值来完成彩色数字水印的嵌入与盲提取。本发明将彩色数字水印嵌入到彩色宿主图像中,不但具有较强的鲁棒性,而且具有较好的水印隐蔽性和较大的嵌入容量,解决了大容量彩色图像数字水印算法鲁棒性差的难题,适用于强鲁棒性的大容量数字媒体版权保护的场合。
Description
技术领域
本发明属于信息安全技术领域,涉及强鲁棒和大容量彩色数字图像的版权保护。
背景技术
随着计算机网络和通信技术的飞速发展,越来越多的彩色数字图像在网络上涌现、传播,随之而来的盗版、侵权等非法行为屡见不鲜,其版权保护问题变得日益突出,急需解决。为此,一方面要求版权保护的标识趋向于美观实用、大容量的彩色图像数字水印,而不再局限于容量较小的伪随机序列、二值图像或灰度图像;另一方面,在网络传输过程中,彩色数字图像会受到各种各样未知的客观因素影响以及人为因素的改变,需要进一步提高其算法的鲁棒性。
因此,如何充分利用频域数字水印算法鲁棒性强的优点设计一种水印容量大、隐蔽性好和鲁棒性强的彩色图像数字水印算法成为目前亟待解决的问题之一。
发明内容
本发明的目的是提供一种基于矩阵奇异值分解的四元数彩色数字图像盲水印方法,其特征在于通过具体的水印嵌入过程和提取过程来实现的,其水印嵌入过程描述如下:
第一步:彩色图像数字水印的预处理:首先,将一幅大小为N×N的彩色图像数字水印W依照红、绿、蓝三基色的顺序分成3个分层水印图像W i ;然后,将每个分层水印图像进行基于密钥Ka i 的Arnold置乱变换以提高水印的安全性;最后,将置乱后的分层水印图像中的每个十进制数表示的像素用8位二进制数表示,并依次连接形成长度为8N 2 的分层水印位序列SW i ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第二步:获取宿主图像的嵌入块:将一幅大小为M×M的原始彩色宿主图像H分成大小为m×m的图像块,每一个图像块纵向包含红、绿、蓝三层;根据分层水印信息长度8N 2 ,利用Matlab系统内置函数randperm生成的伪随机序列在宿主图像H中选择合适位置的图像块,以提高水印抗剪切攻击的鲁棒性,其中3×8N 2<=(M×M)/(m×m);
第三步:选取一个图像块A,进行是否调整极端像素值的判断,如果图像块A蓝色通道的(1, 1)位置上的像素值为0,那么将其像素值改为1;然后,将彩色宿主图像H依照红、绿、蓝三基色的顺序分成3个分层宿主图像H i ,图像块A的每个像素值都通过H i 用四元数表示出来,得到四元数像素块qA,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第四步:根据公式(1),对四元数像素块qA进行奇异值分解(Singular ValueDecomposition, SVD),得到左奇异矩阵U、右奇异矩阵V和对角矩阵S,其中S仅在主对角线上有奇异值;
[U,S,V]=svd(qA) (1)
第五步:从分层水印序列SW i 中按先后顺序取出一个嵌入水印信息w,依据该嵌入水印信息及公式(2)求得量化值sq;
其中,mval=mod(sm,T),sm是对角矩阵S中的最大值,T为量化步长,mod(.)为求余函数,T1=0.25×T,T2=0.75×T,judge1(mval)是当mval<=T/2+T/8或mval>=T-T/8时返回true,否则返回false;judge0(mval)是当mval<=T/8或mval>=T/2-T/8时返回true,否则返回false;
第六步:用量化值sq替换sm,得到S * ;再利用公式(3)进行逆奇异值分解,得到含水印的四元数像素块qA * ;
qA * =U×S * ×V T (3)
第七步:将含水印的四元数像素块qA * 转换为含水印的像素块A * ,然后进行溢出判断:如果sq≠sm,就按照公式(4)进行溢出判断得到sq * ,并用其替换sm * ;之后再利用公式(3)进行逆奇异值分解,得到最终的含水印的四元数像素块qA * ,再进一步转换为最终的含水印的像素块A * ;
第八步:重复执行上述第三步到第七步,直到所有的水印信息都被嵌入完成为止;最后,获得含水印图像H * ;
其水印提取过程描述如下:
第一步:将含水印图像H * 分成m×m的非重叠图像块,每一个图像块纵向包含红、绿、蓝三层;
第二步:利用上述水印嵌入过程中所提到的Matlab系统内置函数randperm生成的伪随机序列选择图像块;
第三步:选取一个图像块A * ,将含水印图像H * 依照红、绿、蓝三基色的顺序分成3个分层含水印图像H i * ,图像块A * 的每个像素值都通过H i * 用四元数表示出来,得到含水印的四元数像素块为qA * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第四步:根据公式(5),对含水印的四元数像素块qA * 进行奇异值分解,得到对角矩阵S * ;
S * =svd(qA * ) (5)
第五步:利用公式(6),提取含水印的四元数像素块qA * 中所含有的水印w * ;
其中,sm是对角矩阵S * 中的最大值,mod(.)是求余函数,T为量化步长;
第六步:重复执行第三步到第五步,提取每层二进制水印序列SW i * ,然后每8位二进制信息为一组转换成十进制的像素值,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第七步:对转化后的每层十进制像素执行基于密钥Ka i 的逆Arnold变换并获得分层水印W i * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第八步:组合获得的分层水印W i * 形成最终的提取水印W * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层。
该方法利用矩阵奇异值分解后对角矩阵S中的系数分布规律,通过对其量化完成了数字水印嵌入与盲提取;该方法既具有较强的鲁棒性,又具有较好的水印不可见性和较大的嵌入容量。
附图说明
图1(a)、图1(b)是两幅原始彩色宿主图像。
图2(a)、图2(b)是两幅原始彩色水印图像。
图3(a)、图3(b)是将图2(a)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像,其结构相似度SSIM值依次是0.9488、0.9365,其峰值信噪比PSNR值依次是39.2818dB、38.8812dB。
图4(a)、图4(b)是依次从图3(a)、图3(b)中提取的水印,其归一化互相关系数NC值分别是1.0000、1.0000。
图5(a)、图5(b)、图5(c)、图5(d)、图5(e)、图5(f)、图5(g)是将图3(a)所示的含水印图像依次进行JPEG70、JPEG2000压缩(5:1)、高斯噪声(0,0.001)、椒盐噪声(2%)、旋转(45°)、缩放(400%)、高斯低通滤波(3×3)等攻击后所提取的水印,其归一化互相关系数NC值分别是0.9849、0.9992、0.9040、0.9899、0.9541、0.9994、0.9666。
图6(a)、图6(b)是将图2(b)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像,其结构相似度SSIM值依次是0.9447、0.9309,其峰值信噪比PSNR值依次是38.8930dB、38.5167dB。
图7(a)、图7(b)是依次从图6(a)、图6(b)中提取的水印,其归一化互相关系数NC值分别是1.0000、1.0000。
图8(a)、图8(b)、图8(c)、图8(d)、图8(e)、图8(f)、图8(g)是将图6(a)所示的含水印图像依次进行JPEG70、JPEG2000压缩(5:1)、高斯噪声(0,0.001)、椒盐噪声(2%)、旋转(45°)、缩放(400%)、高斯低通滤波(3×3)等攻击后所提取的水印,其归一化互相关系数NC值分别是0.9896、0.9997、0.9364、0.9936、0.9674、0.9993、0.9746。
具体实施方式
本发明的目的是提供一种基于矩阵奇异值分解的四元数彩色数字图像盲水印方法,其特征在于通过具体的水印嵌入过程和提取过程来实现的,其水印嵌入过程描述如下:
第一步:彩色图像数字水印的预处理:首先,将一幅大小为32×32的彩色图像数字水印W依照红、绿、蓝三基色的顺序分成3个分层水印图像W i ;然后,将每个分层水印图像进行基于密钥Ka i 的Arnold置乱变换以提高水印的安全性;最后,将置乱后的分层水印图像中的每个十进制数表示的像素用8位二进制数表示(例如:可将205转换成二进制数11001101),并依次连接形成长度为8×322 =8192的分层水印位序列SW i ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第二步:获取宿主图像的嵌入块:将一幅大小为512×512的原始彩色宿主图像H分成大小为2×2的图像块,每一个图像块纵向包含红、绿、蓝三层;根据分层水印信息长度,利用基于Matlab系统内置函数randperm生成的伪随机序列在宿主图像H中选择合适位置的图像块,以提高水印抗剪切攻击的鲁棒性,其中3×8192<=(512×512)/(2×2);
第三步:选取一个图像块A,进行是否调整极端像素值的判断,如果图像块A蓝色通道的(1,1)位置上的像素值为0,那么将其像素值改为1;然后,将彩色宿主图像H依照红、绿、蓝三基色的顺序分成3个分层宿主图像H i ,图像块A的每个像素值都通过H i 用四元数表示出来,得到四元数像素块qA,其中i=1, 2, 3, 分别表示红、绿、蓝三层;此处,设选取的图像块A的红色通道的像素值为,绿色通道的像素值为,蓝色通道的像素值为,可知A的蓝色通道的(1,1)位置上的像素值不为0,所以无需调整其像素值,然后得到的四元数像素块qA(1,1)为w=0, x=194, y=81, z=87;qA(1,2)为w=0, x=203,y=88, z=91;qA(2,1)为w=0, x=195, y=82, z=88;qA(2,2)为w=0, x=198, y=83, z=86;
第四步:根据公式(1),对四元数像素块qA进行奇异值分解(Singular ValueDecomposition, SVD),得到左奇异矩阵U、右奇异矩阵V和对角矩阵S,其中S仅在主对角线上有奇异值;
[U,S,V]=svd(qA) (1)
此处,根据公式(1),得到的左奇异矩阵U(1,1)为w=-2.5290e-05, x=0.6055, y=0.2555, z=0.2741;U(1,2)为w=-0.0022, x=0.4034, y=0.4028, z=0.4099;U(2,1)为w=2.5290e-05, x=0.5992, y=0.2493, z=0.2678;U(2,2)为w=0.0022, x=-0.4126, y=-0.4067, z=-0.4141;右奇异矩阵V(1,1)为w=0.6964, x=0, y=0, z=0;V(1,2)为w=-0.7177, x=0, y=0, z=0;V(2,1)为w=0.7176, x=-0.0036, y=0.0044, z=0.0039;V(2,2)为w=0.6964, x=-0.0035, y=0.0043, z=0.0037;对角矩阵S= ;
第五步:从分层水印序列SW i 中按先后顺序取出一个嵌入水印信息w,依据该嵌入水印信息及公式(2)求得量化值sq;
其中,mval=mod(sm,T),sm是对角矩阵S中的最大值,T为量化步长,mod(.)为求余函数,T1=0.25×T,T2=0.75×T,judge1(mval)是当mval<=T/2+T/8或mval>=T-T/8时返回true,否则返回false;judge0(mval)是当mval<=T/8或mval>=T/2-T/8时返回true,否则返回false;此时,sm=463.6330,设w=‘1’,T=54,则根据公式(2),得sq=472.5000;
第六步:用量化值sq替换sm,得到S * ;再利用公式(3)进行逆奇异值分解,得到含水印的四元数像素块qA * ;
qA * =U×S * ×V T (3)
此处得到S * = ,则根据公式(3),得qA * (1,1)为w=-1.5617e-04, x=197.7387, y=82.5776, z=88.5776;qA * (1,2)为w=-3.1233e-05, x=206.8547, y=89.6551, z=92.7126;qA * (2,1)为w=1.5617e-04, x=198.7003, y=83.5395, z=89.6538;qA * (2,2)为w=3.1233e-05, x=201.8151, y=84.6155, z=87.6728;
第七步:将含水印的四元数像素块qA * 转换为含水印的像素块A * ,然后进行溢出判断:如果sq≠sm,就按照公式(4)进行溢出判断得到sq * ,并用其替换sm * ;之后再利用公式(3)进行逆奇异值分解,得到最终的含水印的四元数像素块qA * ,再进一步转换为最终的含水印的像素块A * ;
其中,sm * 是对角矩阵S * 中的最大值,T为量化步长,是含水印的像素块A * 在第p行第q列的像素值,1≤p,q≤m,m是图像块的尺寸大小;此处,A * 的红色通道的像素值为,绿色通道的像素值为,蓝色通道的像素值为,又由第五步知sq=472.5000,sm=463.6330,得sq≠sm,则根据公式(4),得sq值保持不变,所以A * 就是真正含水印的图像块。
第八步:重复执行上述第三步到第七步,直到所有的水印信息都被嵌入完成为止;最后,获得含水印图像H * ;
其水印提取过程描述如下:
第一步:将含水印图像H * 分成2×2的非重叠图像块,每一个图像块纵向包含红、绿、蓝三层;
第二步:利用上述水印嵌入过程中所提到的Matlab系统内置函数randperm生成的伪随机序列选择图像块;
第三步:选取一个图像块A * ,将含水印图像H * 依照红、绿、蓝三基色的顺序分成3个分层含水印图像H i * ,图像块A * 的每个像素值都通过H i * 用四元数表示出来,得到含水印的四元数像素块为qA * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;此处,设选取的图像块A * 的红色通道的像素值为,绿色通道的像素值为,蓝色通道的像素值为,然后得到的四元数像素块qA * (1,1)为w=0, x=198, y=83, z=89;qA * (1,2)为w=0, x=207, y=90, z=93;qA * (2,1)为w=0, x=199, y=84, z=90;qA * (2,2)为w=0, x=202, y=85, z=88;
第四步:根据公式(5),对含水印的四元数像素块qA * 进行奇异值分解,得到对角矩阵S * ;
S * =svd(qA * ) (5)
第五步:利用公式(6),提取含水印的四元数像素块qA * 中所含有的水印w * ;
其中,sm是对角矩阵S * 中的最大值,mod(.)是求余函数,T为量化步长;此时,sm=473.4084,T=54,利用公式(6)提取所含有的水印w * =‘1’;
第六步:重复执行第三步到第五步,提取每层二进制水印序列SW i * ,然后每8位二进制信息为一组转换成十进制的像素值,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第七步:对转化后的每层十进制像素执行基于密钥Ka i 的逆Arnold变换并获得分层水印W i * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第八步:组合获得的分层水印W i * 形成最终的提取水印W * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层。
该方法具有频域数字水印算法的优点,鲁棒性强,水印不可见性好;算法简单,适用于彩色图像作为数字水印的高效版权保护。
本发明有效性验证
为了证明本发明的有效性,选择如图1(a)、图1(b)所示的两幅大小为512×512的24位标准图像作为宿主图像,并分别用如图2(a)、图2(b)所示的两幅大小为32×32的24位彩色图像作为数字水印进行验证。
图3(a)、图3(b)是将图2(a)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像,其结构相似度SSIM值依次是0.9488、0.9365
,其峰值信噪比PSNR值依次是39.2818dB、38.8812dB;图4(a)、图4(b)是依次从图3(a)、图3(b)中提取的水印,其归一化互相关系数NC值分别是1.0000、1.0000;图5(a)、图5(b)、图5(c)、图5(d)、图5(e)、图5(f)、图5(g)是将图3(a)所示的含水印图像依次进行JPEG70、JPEG2000压缩(5:1)、高斯噪声(0,0.001)、椒盐噪声(2%)、旋转(45°)、缩放(400%)、高斯低通滤波(3×3)等攻击后所提取的水印,其归一化互相关系数NC值分别是0.9849、0.9992、0.9040、0.9899、0.9541、0.9994、0.9666。
图6(a)、图6(b)是将图2(b)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像,其结构相似度SSIM值依次是0.9447、0.9309,其峰值信噪比PSNR值依次是38.8930dB、38.5167dB;图7(a)、图7(b)是依次从图6(a)、图6(b)中提取的水印,其归一化互相关系数NC值分别是1.0000、1.0000;图8(a)、图8(b)、图8(c)、图8(d)、图8(e)、图8(f)、图8(g)是将图6(a)所示的含水印图像依次进行JPEG70、JPEG2000压缩(5:1)、高斯噪声(0,0.001)、椒盐噪声(2%)、旋转(45°)、缩放(400%)、高斯低通滤波(3×3)等攻击后所提取的水印,其归一化互相关系数NC值分别是0.9896、0.9997、0.9364、0.9936、0.9674、0.9993、0.9746。
综上所述,所嵌入的彩色图像数字水印具有较好的不可见性,满足了水印算法的不可见性要求;同时,从各种受攻击图像中所提取的彩色图像数字水印具有较好的可鉴别性和较高的NC值,说明该方法具有较强的鲁棒性;另外,实际嵌入彩色数字图像中的二进制信息为3×213位,具有较大的水印容量,满足了大容量的彩色图像数字水印版权保护的需要。
Claims (1)
1.一种基于矩阵奇异值分解的四元数彩色数字图像盲水印方法,其特征在于通过具体的水印嵌入过程和提取过程来实现的,其水印嵌入过程描述如下:
第一步:彩色图像数字水印的预处理:首先,将一幅大小为N×N的彩色图像数字水印W依照红、绿、蓝三基色的顺序分成3个分层水印图像W i ;然后,将每个分层水印图像进行基于密钥Ka i 的Arnold置乱变换以提高水印的安全性;最后,将置乱后的分层水印图像中的每个十进制数表示的像素用8位二进制数表示,并依次连接形成长度为8N 2 的分层水印位序列SW i ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第二步:获取宿主图像的嵌入块:将一幅大小为M×M的原始彩色宿主图像H分成大小为m×m的图像块,每一个图像块纵向包含红、绿、蓝三层;根据分层水印信息长度8N 2 ,利用Matlab系统内置函数randperm生成的伪随机序列在宿主图像H中选择合适位置的图像块,以提高水印抗剪切攻击的鲁棒性,其中3×8N 2<=(M×M)/(m×m);
第三步:选取一个图像块A,进行是否调整极端像素值的判断,如果图像块A蓝色通道的(1, 1)位置上的像素值为0,那么将其像素值改为1;然后,将彩色宿主图像H依照红、绿、蓝三基色的顺序分成3个分层宿主图像H i ,图像块A的每个像素值都通过H i 用四元数表示出来,得到四元数像素块qA,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第四步:根据公式(1),对四元数像素块qA进行奇异值分解(Singular ValueDecomposition, SVD),得到左奇异矩阵U、右奇异矩阵V和对角矩阵S,其中S仅在主对角线上有奇异值;
[U,S,V]=svd(qA) (1)
第五步:从分层水印序列SW i 中按先后顺序取出一个嵌入水印信息w,依据该嵌入水印信息及公式(2)求得量化值sq;
其中,mval=mod(sm,T),sm是对角矩阵S中的最大值,T为量化步长,mod(.)为求余函数,T1=0.25×T,T2=0.75×T,judge1(mval)是当mval<=T/2+T/8或mval>=T-T/8时返回true,否则返回false;judge0(mval)是当mval<=T/8或mval>=T/2-T/8时返回true,否则返回false;
第六步:用量化值sq替换sm,得到S * ;再利用公式(3)进行逆奇异值分解,得到含水印的四元数像素块qA * ;
qA * =U×S * ×V T (3)
第七步:将含水印的四元数像素块qA * 转换为含水印的像素块A * ,然后进行溢出判断:如果sq≠sm,就按照公式(4)进行溢出判断得到sq * ,并用其替换sm * ;之后再利用公式(3)进行逆奇异值分解,得到最终的含水印的四元数像素块qA * ,再进一步转换为最终的含水印的像素块A * ;
第八步:重复执行上述第三步到第七步,直到所有的水印信息都被嵌入完成为止;最后,获得含水印图像H * ;
其水印提取过程描述如下:
第一步:将含水印图像H * 分成m×m的非重叠图像块,每一个图像块纵向包含红、绿、蓝三层;
第二步:利用上述水印嵌入过程中所提到的Matlab系统内置函数randperm生成的伪随机序列选择图像块;
第三步:选取一个图像块A * ,将含水印图像H * 依照红、绿、蓝三基色的顺序分成3个分层含水印图像H i * ,图像块A * 的每个像素值都通过H i * 用四元数表示出来,得到含水印的四元数像素块为qA * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第四步:根据公式(5),对含水印的四元数像素块qA * 进行奇异值分解,得到对角矩阵S * ;
S * =svd(qA * ) (5)
第五步:利用公式(6),提取含水印的四元数像素块qA * 中所含有的水印w * ;
其中,sm是对角矩阵S * 中的最大值,mod(.)是求余函数,T为量化步长;
第六步:重复执行第三步到第五步,提取每层二进制水印序列SW i * ,然后每8位二进制信息为一组转换成十进制的像素值,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第七步:对转化后的每层十进制像素执行基于密钥Ka i 的逆Arnold变换并获得分层水印W i * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层;
第八步:组合获得的分层水印W i * 形成最终的提取水印W * ,其中i=1, 2, 3, 分别表示红、绿、蓝三层。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110059404.8A CN112488904B (zh) | 2021-01-18 | 2021-01-18 | 基于矩阵奇异值分解的四元数彩色数字图像盲水印方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110059404.8A CN112488904B (zh) | 2021-01-18 | 2021-01-18 | 基于矩阵奇异值分解的四元数彩色数字图像盲水印方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112488904A true CN112488904A (zh) | 2021-03-12 |
CN112488904B CN112488904B (zh) | 2023-09-29 |
Family
ID=74912308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110059404.8A Active CN112488904B (zh) | 2021-01-18 | 2021-01-18 | 基于矩阵奇异值分解的四元数彩色数字图像盲水印方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112488904B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113139890A (zh) * | 2021-05-18 | 2021-07-20 | 辽宁工程技术大学 | 一种基于子块均值的空域盲水印算法 |
CN113222803A (zh) * | 2021-05-31 | 2021-08-06 | 鲁东大学 | 基于纯四元数Schur分解的彩色数字图像盲水印方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999065241A1 (en) * | 1998-06-08 | 1999-12-16 | Digital Video Express, L.P. | System and methodology for tracing to a source of unauthorized copying of prerecorded proprietary material, such as movies |
CN102880999A (zh) * | 2012-09-21 | 2013-01-16 | 吉林大学 | 抗打印、印刷和扫描过程的彩色图像数字水印方法 |
CN107274333A (zh) * | 2017-06-21 | 2017-10-20 | 鲁东大学 | 一种融合奇异值分解的空域彩色图像盲水印方法 |
CN107895340A (zh) * | 2017-11-29 | 2018-04-10 | 鲁东大学 | 一种融合qr分解的空频域彩色数字图像盲水印方法 |
CN109035129A (zh) * | 2018-08-15 | 2018-12-18 | 鲁东大学 | 一种基于二维离散正弦变换的彩色数字图像盲水印方法 |
CN109829845A (zh) * | 2019-01-29 | 2019-05-31 | 鲁东大学 | 基于矩阵Schur分解的可变步长彩色图像盲水印方法 |
CN110349073A (zh) * | 2019-07-24 | 2019-10-18 | 鲁东大学 | 基于Schur分解的四进制彩色数字图像盲水印方法 |
-
2021
- 2021-01-18 CN CN202110059404.8A patent/CN112488904B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999065241A1 (en) * | 1998-06-08 | 1999-12-16 | Digital Video Express, L.P. | System and methodology for tracing to a source of unauthorized copying of prerecorded proprietary material, such as movies |
CN102880999A (zh) * | 2012-09-21 | 2013-01-16 | 吉林大学 | 抗打印、印刷和扫描过程的彩色图像数字水印方法 |
CN107274333A (zh) * | 2017-06-21 | 2017-10-20 | 鲁东大学 | 一种融合奇异值分解的空域彩色图像盲水印方法 |
CN107895340A (zh) * | 2017-11-29 | 2018-04-10 | 鲁东大学 | 一种融合qr分解的空频域彩色数字图像盲水印方法 |
CN109035129A (zh) * | 2018-08-15 | 2018-12-18 | 鲁东大学 | 一种基于二维离散正弦变换的彩色数字图像盲水印方法 |
CN109829845A (zh) * | 2019-01-29 | 2019-05-31 | 鲁东大学 | 基于矩阵Schur分解的可变步长彩色图像盲水印方法 |
CN110349073A (zh) * | 2019-07-24 | 2019-10-18 | 鲁东大学 | 基于Schur分解的四进制彩色数字图像盲水印方法 |
Non-Patent Citations (3)
Title |
---|
褚静;徐安成;张美凤;: "DWT和SVD相融合的彩色图像水印算法", 电视技术, no. 17 * |
郭巧丹;吴锡生;: "基于SVD的彩色图像盲水印算法", 计算机工程与科学, no. 11 * |
陈北京;田翠翠;戴慧;王定成;舒华忠;: "四元数调制重叠变换及其在双彩色图像数字水印中的应用", 东南大学学报(自然科学版), no. 06 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113139890A (zh) * | 2021-05-18 | 2021-07-20 | 辽宁工程技术大学 | 一种基于子块均值的空域盲水印算法 |
CN113139890B (zh) * | 2021-05-18 | 2024-04-19 | 辽宁工程技术大学 | 一种基于子块均值的空域盲水印算法 |
CN113222803A (zh) * | 2021-05-31 | 2021-08-06 | 鲁东大学 | 基于纯四元数Schur分解的彩色数字图像盲水印方法 |
CN113222803B (zh) * | 2021-05-31 | 2022-03-29 | 鲁东大学 | 基于纯四元数Schur分解的彩色数字图像盲水印方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112488904B (zh) | 2023-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108648134B (zh) | 一种融合离散小波变换的空域彩色数字图像盲水印方法 | |
CN109035129B (zh) | 一种基于二维离散正弦变换的彩色数字图像盲水印方法 | |
CN107895340B (zh) | 一种融合qr分解的空频域彩色数字图像盲水印方法 | |
CN111199508B (zh) | 一种融合dct和dht的空域彩色数字图像盲水印方法 | |
CN109829845B (zh) | 基于矩阵Schur分解的可变步长彩色图像盲水印方法 | |
CN110349073B (zh) | 基于Schur分解的四进制彩色数字图像盲水印方法 | |
CN110390621A (zh) | 一种基于可变步长的dct域彩色数字图像盲水印方法 | |
CN112488904B (zh) | 基于矩阵奇异值分解的四元数彩色数字图像盲水印方法 | |
CN109102454A (zh) | 一种融合快速傅里叶变换的彩色qr码数字盲水印方法 | |
CN104050622B (zh) | 基于三进制编码的彩色图像盲水印方法 | |
CN106204410B (zh) | 一种新型的基于矩阵Schur分解的数字水印方法 | |
Muñoz-Ramirez et al. | A robust watermarking scheme to JPEG compression for embedding a color watermark into digital images | |
CN112508765B (zh) | 一种基于沃尔什-哈达玛变换的频域彩色数字图像盲水印方法 | |
CN106991636B (zh) | 一种融合近似Schur分解的空域彩色图像盲水印方法 | |
CN110570345B (zh) | 一种融合离散余弦变换的空域彩色数字图像盲水印方法 | |
CN105023236A (zh) | 一种基于矩阵lu分解的双彩色图像盲水印方法 | |
CN109829846B (zh) | 一种基于二维离散余弦变换的数字图像盲水印方法 | |
CN116993567A (zh) | 一种基于哈达玛变换和教学优化算法的频域盲水印方法 | |
CN106169171B (zh) | 基于Hessenberg分解的高性能数字水印方法 | |
CN111242828B (zh) | 一种融合离散傅里叶变换的空域彩色数字图像盲水印方法 | |
CN110415154B (zh) | 一种基于哈尔变换的四进制彩色数字图像盲水印方法 | |
CN113191932B (zh) | 一种融合离散切比雪夫变换的空域彩色数字图像盲水印方法 | |
CN113222803B (zh) | 基于纯四元数Schur分解的彩色数字图像盲水印方法 | |
CN112017098B (zh) | 一种自适应彩色数字图像盲水印方法 | |
CN112488903B (zh) | 一种融合多级离散傅里叶变换的空域彩色数字图像盲水印方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |