CN107274333B - 一种融合奇异值分解的空域彩色图像盲水印方法 - Google Patents

一种融合奇异值分解的空域彩色图像盲水印方法 Download PDF

Info

Publication number
CN107274333B
CN107274333B CN201710472958.4A CN201710472958A CN107274333B CN 107274333 B CN107274333 B CN 107274333B CN 201710472958 A CN201710472958 A CN 201710472958A CN 107274333 B CN107274333 B CN 107274333B
Authority
CN
China
Prior art keywords
watermark
image
singular value
image block
layered
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.)
Active
Application number
CN201710472958.4A
Other languages
English (en)
Other versions
CN107274333A (zh
Inventor
苏庆堂
苏琳
王环英
马艳萍
薛梅
巩艳华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SHANDONG DAOWAN ELECTRIC Co.,Ltd.
Original Assignee
Ludong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ludong University filed Critical Ludong University
Priority to CN201710472958.4A priority Critical patent/CN107274333B/zh
Publication of CN107274333A publication Critical patent/CN107274333A/zh
Application granted granted Critical
Publication of CN107274333B publication Critical patent/CN107274333B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • G06T1/0021Image watermarking
    • G06T1/005Robust watermarking, e.g. average attack or collusion attack resistant
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2201/00General purpose image data processing
    • G06T2201/005Image watermarking
    • G06T2201/0083Image 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个二维分层水印图像Wi,并将分层水印图像Wi进行基于Kai的Arnold置乱变换以提高其安全性;然后,将每个十进制表示的像素用8位二进制数表示,并依次组合形成分层水印序列Wi,其中i∈{1,2,3}分别表示红、绿、蓝三层;
第二步:将大小为M×M原始宿主图像H分成3个分层宿主图像Hi,其中i∈{1,2,3}分别表示红、绿、蓝三层;同时,把每一个分层宿主图像Hi分成大小为n×n的图像块;利用基于密钥K的图像块间隔选择法在分层宿主图像Hi中选择嵌入水印位置的图像块:即在奇数行依序选取第qK-1个图像块,在偶数行依序选取第qK个图像块,其中q∈[1,floor(M/(n*K))],K是块间隔值并作为密钥,M,n分别是宿主图像和图像块的尺寸,floor(.)是向下取整函数;
第三步:选取一个图像块A,并按照公式(1)在空域中直接求出图像块A经过奇异值分解后对应的最大奇异值λ1
Figure GDA0002386666650000021
其中,||·||2表示矩阵的二范数,ap,q∈A,1≤p,q≤n,p,q分别表示像素在n×n的图像块A中所处的行号和列号,图像块A经过奇异值分解后的最大奇异值λ1在空域中利用公式(1)直接求得而不需通过真正的奇异值分解后获得;
第四步:在空域中,依据水印W及公式(2)或(3),计算当将水印W量化嵌入到奇异值分解的最大奇异值λ1上时,将得到两个极限值
Figure GDA0002386666650000022
Figure GDA0002386666650000023
如果w=0,则
Figure GDA0002386666650000024
如果w=1,则
Figure GDA0002386666650000025
其中T是量化步长;
第五步:依据公式(4),理论上获得最大奇异值λ1在嵌入水印W后的最终结果
Figure GDA0002386666650000026
Figure GDA0002386666650000027
其中abs(.)是绝对值函数;
第六步:依据公式(5),计算最大奇异值λ1的变化量为Δλ1
Figure GDA0002386666650000028
第七步:依据公式(6),计算图像块在嵌入水印后每个像素的平均修改量ΔEaverage
ΔEaverage=Δλ1/n (6)
该公式推导如下:
对于n×n的图像块A,其奇异值分解可表示为:
Figure GDA0002386666650000031
其中U,V是正交矩阵,D是对角矩阵,[U1,U2,…,Un],[V1,V2,…,Vn]分别是形成矩阵U,V的特征向量,且对角矩阵中的奇异值满足λ1≥λ2≥…≥λr>λr+1=…=λn=0;当最大奇异值λ1的变化量是Δλ1时,则整个图像块的变化量ΔA=Δλ1U1V1 T
当图像块中像素值变化时,则U1V1 T=n,故ΔA=nΔλ1
因为图像块中共有n×n个像素,所以每个像素的平均变化量
ΔEaverage=ΔA/(n×n)=nΔλ1/(n×n)=Δλ1/n;
第八步:在空域中将每个像素的平均修改量ΔEaverage依次加到图像块A的所有像素上,这样在空域中完成对奇异值分解最大奇异值的修改并嵌入了水印W,此时在空域中得到真正含水印的图像块A*
第九步:重复执行步骤第三步到第八步,直到所有的水印信息都被嵌入完成为止;最后,将含水印的红、绿、蓝分层图像重新组合并获得含水印图像H*
其水印提取过程描述如下:
第一步:将含水印图像H*分成3个分层含水印图像
Figure GDA0002386666650000032
其中i∈{1,2,3}分别表示红、绿、蓝三层,并将每一分层含水印图像
Figure GDA0002386666650000033
进一步分成n×n的非重叠图像块;
第二步:在分层图像
Figure GDA0002386666650000034
中,利用上述水印嵌入过程中所提到的基于密钥K的图像块间隔选择法选择图像块;
第三步:选取一个图像块A*,在空域中利用公式(7)计算其奇异值分解中的最大奇异值
Figure GDA0002386666650000035
Figure GDA0002386666650000036
其中,||·||2表示矩阵的二范数,
Figure GDA0002386666650000041
1≤p,q≤n,p,q分别表示像素在n×n的图像块A*中所处的行号和列号;
第四步:利用公式(8),提取图像块A*中所含有的水印W*
Figure GDA0002386666650000042
其中,mod(.)是求余函数,round(.)是四舍五入函数,T为量化步长;
第五步:重复执行第三步、第四步,直到提取所有的水印信息,并将所提取的二进制水印信息每8位一组转换成十进制的像素值;
第六步:执行基于密钥Kai的逆Arnold变换并获得分层水印,其中i∈{1,2,3};
第七步:组合分层水印形成最终的提取水印W*
该方法利用在空域求得奇异值分解中最大奇异值的原理及最大奇异值的变化量在空域像素上的分布规律,在空域中直接完成了需要进行奇异值分解进行的数字水印嵌入与盲提取,简单快捷;该方法既具有较好的水印鲁棒性,又具有较高的算法实时性。
附图说明
图1(a)、图1(b)是两幅原始彩色宿主图像。
图2(a)、图2(b)是两幅彩色水印图像。
图3(a)、图3(b)是将图2(a)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像,其结构相似度SSIM值依次是0.9701、0.9657,其峰值信噪比PSNR值依次是40.6976dB、40.5906dB。
图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)是将图3(a)所示的含水印图像依次进行JPEG2000压缩(5:1)、椒盐噪声(2%)、中值滤波(3×3)、锐化(1.0)、缩放(400%)、剪切(25%)等攻击后所提取的水印,其归一化互相关系数NC值分别是1.0000、0.9843、0.9661、0.9998、0.9997、0.7523。
图6(a)、图6(b)是将图2(b)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像,其结构相似度SSIM值依次是0.9677、0.9616,其峰值信噪比PSNR值依次是40.7046dB、40.7368dB。
图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)是将图6(a)所示的含水印图像依次进行JPEG2000压缩(5:1)、椒盐噪声(2%)、中值滤波(3×3)、锐化(1.0)、缩放(400%)、剪切(25%)等攻击后所提取的水印,其归一化互相关系数NC值分别是1.0000、0.9819、0.9636、0.9999、0.9997、0.7550。
具体实施方式
本发明的目的是提供一种融合奇异值分解的空域彩色图像盲水印方法,其特征在于通过具体的水印嵌入过程和提取过程来实现的,其水印嵌入过程描述如下:
第一步:数字水印预处理:将大小为32×32三维彩色图像数字水印W按红、绿、蓝的三基色顺序分成3个二维分层水印图像Wi,并将分层水印图像进行基于Kai的Arnold置乱变换以提高其安全性;然后,将每个十进制表示的像素用8位二进制数表示,并依次组合形成分层水印序列Wi,其中i∈{1,2,3}分别表示红、绿、蓝三层;例如可将225,224分别转换成二进制数11100001,11100000,并连接为水印序列1110000111100000;
第二步:将大小为512×512原始宿主图像H分成3个分层宿主图像Hi,其中i∈{1,2,3}分别表示红、绿、蓝三层;同时,把每一个分层宿主图像Hi分成大小为4×4的图像块;利用基于密钥K的图像块间隔选择法在分层宿主图像Hi中选择嵌入水印位置的图像块:即在奇数行依序选取第qK-1个图像块,在偶数行依序选取第qK个图像块,其中q∈[1,floorM/(n+K))],K是块间隔值并作为密钥,M,n分别是宿主图像和图像块的尺寸,floor(.)是向下取整函数;例如:当K=2,M=512,n=4时,则q∈[1,64],在奇数行选取第1,3,…,127个图像块,在偶数行选取第2,4,…,128个图像块;
第三步:选取一个图像块A,并按照公式(1)在空域中直接求出图像块A经过奇异值分解后对应的最大奇异值λ1
Figure GDA0002386666650000051
其中,||·||2表示矩阵的二范数,ap,q∈A,1≤p,q≤n,p,q分别表示像素在n×n的图像块A中所处的行号和列号,图像块A经过奇异值分解后的最大奇异值λ1在空域中利用公式(1)直接求得而不需通过真正的奇异值分解后获得;
此处,设选取的图像块A为
Figure GDA0002386666650000061
利用公式(1)可在空域中得该图像块的最大奇异值为λ1=896.2539;
第四步:在空域中,依据水印W及公式(2)或(3),计算当将水印W量化嵌入到奇异值分解的最大奇异值上λ1时,将得到的两个极限值
Figure GDA0002386666650000062
Figure GDA0002386666650000063
如果w=0,则
Figure GDA0002386666650000064
如果w=1,则
Figure GDA0002386666650000065
其中T是量化步长;此时,设w=0,T=46,则根据公式(2),得
Figure GDA0002386666650000066
Figure GDA0002386666650000067
第五步:依据公式(4),理论上获得最大奇异值λ1在嵌入水印W后的最终结果
Figure GDA0002386666650000068
Figure GDA0002386666650000069
其中abs(.)是绝对值函数;此时,
Figure GDA00023866666500000610
所以根据公式(4),得
Figure GDA00023866666500000611
第六步:依据公式(5),计算最大奇异值λ1的变化量为Δλ1
Figure GDA00023866666500000612
此时,
Figure GDA00023866666500000613
λ1=896.2539,所以根据公式(5),得Δλ1=-10.7539;
第七步:依据公式(6),计算图像块在嵌入水印后每个像素的平均修改量ΔEaverage
ΔEaverage=Δλ1/n (6)
此时,Δλ1=-10.7539,n=4,所以根据公式(6),得ΔEaverage=-2.6885;
第八步:在空域中将每个像素的平均修改量ΔEaverage=-2.6885依次加到图像块A的所有像素上,这样在空域中完成对奇异值分解最大奇异值的修改并嵌入了水印0,此时在空域中得到真正含水印的图像块A*
Figure GDA0002386666650000071
第九步:重复执行步骤第三步到第八步,直到所有的水印信息都被嵌入完成为止;最后,将含水印的红、绿、蓝分层图像重新组合并获得含水印图像H*
其水印提取过程描述如下:
第一步:将含水印图像H*分成3个分层含水印图像
Figure GDA0002386666650000072
其中i∈{1,2,3}分别表示红、绿、蓝三层,并将每一分层含水印图像
Figure GDA0002386666650000073
进一步分成4×4的非重叠图像块;
第二步:在分层图像
Figure GDA0002386666650000074
中,利用上述水印嵌入过程中所提到的基于密钥K的图像块间隔选择法选择图像块;
第三步:选取一个图像块A*,在空域中利用公式(7)计算其奇异值分解中的最大奇异值
Figure GDA0002386666650000075
Figure GDA0002386666650000076
其中,||·||2表示矩阵的二范数,
Figure GDA0002386666650000077
1≤p,q≤n,p,q分别表示像素在n×n的图像块A*中所处的行号和列号;
此时,设选取的图像块A*
Figure GDA0002386666650000078
则利用公式(7)得其最大奇异值
Figure GDA0002386666650000079
第四步:利用公式(8),提取图像块A*中所含有的水印w*
Figure GDA00023866666500000710
其中,mod(.)是求余函数,round(.)是四舍五入函数,T为量化步长;此时,
Figure GDA00023866666500000711
T=46,利用公式(8)提取所含有的水印w*=0;
第五步:重复执行第三步、第四步,直到提取所有的水印信息,并将所提取的二进制水印信息每8位一组转换成十进制的像素值;
第六步:执行基于密钥Kai的逆Arnold变换并获得分层水印,其中i∈{1,2,3};
第七步:组合分层水印形成最终的提取水印W*
该方法不但具有较强的水印鲁棒性,而且具有较高的算法实时性,简单快捷,适用于彩色图像作为数字水印的高效版权保护。
本发明有效性验证
为了证明本发明的有效性,选择如图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.9701、0.9657,其峰值信噪比PSNR值依次是40.6976dB、40.5906dB;图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)是将图3(a)所示的含水印图像依次进行JPEG2000压缩(5:1)、椒盐噪声(2%)、中值滤波(3×3)、锐化(1.0)、缩放(400%)、剪切(25%)等攻击后所提取的水印,其归一化互相关系数NC值分别是1.0000、0.9843、0.9661、0.9998、0.9997、0.7523。
图6(a)、图6(b)是将图2(b)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像,其结构相似度SSIM值依次是0.9677、0.9616,其峰值信噪比PSNR值依次是40.7046dB、40.7368dB;图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)是将图6(a)所示的含水印图像依次进行JPEG2000压缩(5:1)、椒盐噪声(2%)、中值滤波(3×3)、锐化(1.0)、缩放(400%)、剪切(25%)等攻击后所提取的水印,其归一化互相关系数NC值分别是1.0000、0.9819、0.9636、0.9999、0.9997、0.7550。
该算法在平台2.27GHZ CPU,2.00GB RAM,Win 7,MATLAB 7.10.0(R2010a)上进行过近万次运行,其数字水印的平均嵌入时间是0.406780秒,平均提取时间是0.252644秒,总计时间为0.659424秒。
由此可见,所嵌入的彩色图像数字水印具有较好的不可见性;同时,从各种受攻击图像中所提取的数字水印图像具有较好的可鉴别性,说明该方法具有较强的鲁棒性;同时,该算法的平均运行总时间少于1秒,满足了彩色数字图像快速版权保护的需要。

Claims (1)

1.一种融合奇异值分解的空域彩色图像盲水印方法,其特征在于通过具体的水印嵌入过程和提取过程来实现的,其水印嵌入过程描述如下:
第一步:数字水印预处理:将大小为N×N三维彩色图像数字水印W按红、绿、蓝的三基色顺序分成3个二维分层水印图像Wi,并将分层水印图像Wi进行基于Kai的Arnold置乱变换以提高其安全性;然后,将每个十进制表示的像素用8位二进制数表示,并依次组合形成分层水印序列Wi,其中i∈{1,2,3}分别表示红、绿、蓝三层;
第二步:将大小为M×M原始宿主图像H分成3个分层宿主图像Hi,其中i∈{1,2,3}分别表示红、绿、蓝三层;同时,把每一个分层宿主图像Hi分成大小为n×n的图像块;利用基于密钥K的图像块间隔选择法在分层宿主图像Hi中选择嵌入水印位置的图像块:即在奇数行依序选取第qK-1个图像块,在偶数行依序选取第qK个图像块,其中q∈[1,floor(M/(n*K))],K是块间隔值并作为密钥,M,n分别是宿主图像和图像块的尺寸,floor(.)是向下取整函数;
第三步:选取一个图像块A,并按照公式(1)在空域中直接求出图像块A经过奇异值分解后对应的最大奇异值λ1
Figure FDA0002386666640000011
其中,||·||2表示矩阵的二范数,ap,q∈A,1≤p,q≤n,p,q分别表示像素在n×n的图像块A中所处的行号和列号,图像块A经过奇异值分解后的最大奇异值λ1在空域中利用公式(1)直接求得而不需通过真正的奇异值分解后获得;
第四步:在空域中,依据水印W及公式(2)或(3),计算当将水印W量化嵌入到奇异值分解的最大奇异值λ1上时,将得到两个极限值
Figure FDA0002386666640000012
Figure FDA0002386666640000013
如果w=0,则
Figure FDA0002386666640000014
如果w=1,则
Figure FDA0002386666640000015
其中T是量化步长;
第五步:依据公式(4),理论上获得最大奇异值λ1在嵌入水印W后的最终结果
Figure FDA0002386666640000016
Figure FDA0002386666640000021
其中abs(.)是绝对值函数;
第六步:依据公式(5),计算最大奇异值λ1的变化量为Δλ1
Figure FDA0002386666640000022
第七步:依据公式(6),计算图像块在嵌入水印后每个像素的平均修改量ΔEaverage
ΔEaverage=Δλ1/n (6)
该公式推导如下:
对于n×n的图像块A,其奇异值分解可表示为:
Figure FDA0002386666640000023
其中U,V是正交矩阵,D是对角矩阵,[U1,U2,…,Un],[V1,V2,…,Vn]分别是形成矩阵U,V的特征向量,且对角矩阵中的奇异值满足λ1≥λ2≥…≥λr>λr+1=…=λn=0;
当最大奇异值λ1的变化量是Δλ1时,则整个图像块的变化量ΔA=Δλ1U1V1 T
当图像块中像素值变化时,则U1V1 T=n,故ΔA=nΔλ1
因为图像块中共有n×n个像素,所以每个像素的平均变化量
ΔEaverage=ΔA/(n×n)=nΔλ1/(n×n)=Δλ1/n;
第八步:在空域中将每个像素的平均修改量ΔEaverage依次加到图像块A的所有像素上,这样在空域中完成对奇异值分解最大奇异值的修改并嵌入了水印W,此时在空域中得到真正含水印的图像块A*
第九步:重复执行步骤第三步到第八步,直到所有的水印信息都被嵌入完成为止;最后,将含水印的红、绿、蓝分层图像重新组合并获得含水印图像H*
其水印提取过程描述如下:
第一步:将含水印图像H*分成3个分层含水印图像
Figure FDA0002386666640000031
其中i∈{1,2,3}分别表示红、绿、蓝三层,并将每一分层含水印图像
Figure FDA0002386666640000032
进一步分成n×n的非重叠图像块;
第二步:在分层图像
Figure FDA0002386666640000033
中,利用上述水印嵌入过程中所提到的基于密钥K的图像块间隔选择法选择图像块;
第三步:选取一个图像块A*,在空域中利用公式(7)计算其奇异值分解中的最大奇异值
Figure FDA0002386666640000034
Figure FDA0002386666640000035
其中,||·||2表示矩阵的二范数,
Figure FDA0002386666640000036
1≤p,q≤n,p,q分别表示像素在n×n的图像块A*中所处的行号和列号;
第四步:利用公式(8),提取图像块A*中所含有的水印w*
Figure FDA0002386666640000037
其中,mod(.)是求余函数,round(.)是四舍五入函数,T为量化步长;
第五步:重复执行第三步、第四步,直到提取所有的水印信息,并将所提取的二进制水印信息每8位一组转换成十进制的像素值;
第六步:执行基于密钥Kai的逆Arnold变换并获得分层水印,其中i∈{1,2,3};
第七步:组合分层水印形成最终的提取水印W*
CN201710472958.4A 2017-06-21 2017-06-21 一种融合奇异值分解的空域彩色图像盲水印方法 Active CN107274333B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710472958.4A CN107274333B (zh) 2017-06-21 2017-06-21 一种融合奇异值分解的空域彩色图像盲水印方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710472958.4A CN107274333B (zh) 2017-06-21 2017-06-21 一种融合奇异值分解的空域彩色图像盲水印方法

Publications (2)

Publication Number Publication Date
CN107274333A CN107274333A (zh) 2017-10-20
CN107274333B true CN107274333B (zh) 2020-05-05

Family

ID=60068565

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710472958.4A Active CN107274333B (zh) 2017-06-21 2017-06-21 一种融合奇异值分解的空域彩色图像盲水印方法

Country Status (1)

Country Link
CN (1) CN107274333B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110211016B (zh) * 2018-02-28 2022-11-01 佛山科学技术学院 一种基于卷积特征的水印嵌入方法
CN109543425B (zh) * 2018-10-26 2023-03-10 浙江师范大学 一种基于张量分解的图像数据隐藏方法
CN112488904B (zh) * 2021-01-18 2023-09-29 鲁东大学 基于矩阵奇异值分解的四元数彩色数字图像盲水印方法
CN116228507B (zh) * 2023-03-14 2023-10-27 上海阅文信息技术有限公司 一种数字文本水印溯源方法及系统
CN116563081B (zh) * 2023-07-06 2023-09-26 云南日报报业集团 一种空间域盲图像水印方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103810668A (zh) * 2014-02-19 2014-05-21 鲁东大学 一种基于奇异值分解的双彩色图像盲水印方法
CN105761199A (zh) * 2016-03-25 2016-07-13 鲁东大学 一种基于组合域的鲁棒彩色图像水印方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103810668A (zh) * 2014-02-19 2014-05-21 鲁东大学 一种基于奇异值分解的双彩色图像盲水印方法
CN105761199A (zh) * 2016-03-25 2016-07-13 鲁东大学 一种基于组合域的鲁棒彩色图像水印方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
" A dual color images watermarking scheme based on the optimized compensation of singular value decomposition";Qingtang Su.et al;《AEU - International Journal of Electronics and Communications 》;20130831;第652-664页 *
"A blind color image watermarking based on DC component in the spatial domain";Qingtang Su .et al;《Optik》;20131231;第6255-6260页 *
"A blind dual color images watermarking based on singular value decomposition";Qingtang Su.et al;《Applied Mathematics and Computation》;20130415;第219卷;第8455-8466页 *

Also Published As

Publication number Publication date
CN107274333A (zh) 2017-10-20

Similar Documents

Publication Publication Date Title
CN107274333B (zh) 一种融合奇异值分解的空域彩色图像盲水印方法
CN107895340B (zh) 一种融合qr分解的空频域彩色数字图像盲水印方法
CN108648134B (zh) 一种融合离散小波变换的空域彩色数字图像盲水印方法
Narayana et al. Visual cryptography for gray scale images with enhanced security mechanisms
CN110390621B (zh) 一种基于可变步长的dct域彩色数字图像盲水印方法
CN112381701B (zh) 一种基于FFST和Hessenberg分解的抗RST攻击立体图像零水印方法
CN109102454B (zh) 一种融合快速傅里叶变换的彩色qr码数字盲水印方法
CN104053074B (zh) 基于深度图像和Otsu分割的视频水印方法
Patel et al. Secure and authentic DCT image steganography through DWT-SVD based digital watermarking with RSA encryption
CN111199508B (zh) 一种融合dct和dht的空域彩色数字图像盲水印方法
CN106204410B (zh) 一种新型的基于矩阵Schur分解的数字水印方法
CN106991636B (zh) 一种融合近似Schur分解的空域彩色图像盲水印方法
Muñoz-Ramirez et al. A robust watermarking scheme to JPEG compression for embedding a color watermark into digital images
Singh High PSNR based image steganography
CN110570345B (zh) 一种融合离散余弦变换的空域彩色数字图像盲水印方法
CN112508765A (zh) 一种基于沃尔什-哈达玛变换的频域彩色数字图像盲水印方法
CN112488904A (zh) 基于矩阵奇异值分解的四元数彩色数字图像盲水印方法
CN111242828B (zh) 一种融合离散傅里叶变换的空域彩色数字图像盲水印方法
Kaur et al. Comparative analyses of YCbCr color space and CIELab color space based on DWT and SVD
Veličković et al. Digital video protection in the DWT-SVD domain using scrambled watermark by GMSAT algorithm
CN113222803B (zh) 基于纯四元数Schur分解的彩色数字图像盲水印方法
Maheshwari et al. Robust multi-modal watermarking using visually encrypted watermark
CN113191933B (zh) 一种基于近似哈达玛变换的彩色数字图像盲水印方法
Venugopala et al. Evaluation of video watermarking algorithms on mobile device
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
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A spatial color image blind watermarking method based on singular value decomposition

Effective date of registration: 20211216

Granted publication date: 20200505

Pledgee: Yantai financing guarantee Group Co.,Ltd.

Pledgor: LUDONG University

Registration number: Y2021980015152

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220130

Address after: 264000 Yantai Beihang science and Technology Park, No. 77, Hangtian Road, high tech Zone, Yantai City, Shandong Province

Patentee after: SHANDONG DAOWAN ELECTRIC Co.,Ltd.

Address before: 264025 No. 186 Hongqi Middle Road, Zhifu District, Shandong, Yantai

Patentee before: LUDONG University

PC01 Cancellation of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20220311

Granted publication date: 20200505

Pledgee: Yantai financing guarantee Group Co.,Ltd.

Pledgor: LUDONG University

Registration number: Y2021980015152