CN108564519A - 一种基于多变换域的彩色图像数字水印方法 - Google Patents
一种基于多变换域的彩色图像数字水印方法 Download PDFInfo
- Publication number
- CN108564519A CN108564519A CN201711401354.7A CN201711401354A CN108564519A CN 108564519 A CN108564519 A CN 108564519A CN 201711401354 A CN201711401354 A CN 201711401354A CN 108564519 A CN108564519 A CN 108564519A
- Authority
- CN
- China
- Prior art keywords
- matrix
- watermark
- component
- key
- sub
- 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
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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/0065—Extraction of an embedded watermark; Reliable detection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Editing Of Facsimile Originals (AREA)
- Image Processing (AREA)
Abstract
本发明提出一种基于多变换域的彩色图像数字水印方法,结合离散小波变换(DWT)、四元数离散余弦变换(QDCT)和奇异值分解(SVD)。首先将彩色图像分通道DWT,再将三个通道同类频率子带用四元数描述,然后根据水印尺寸大小进行分块,并将每个子块进行DQCT,并组成一个系数矩阵,最后,在保证视觉可感知质量下,将水印嵌入到系数矩阵中。本发明提供一种新的彩色图像水印方法,针对现有彩色图像数字水印算法在平衡鲁棒性与透明性关系上还具有提升空间,充分利用各个变换方法的优势,在保证透明性前提下,提高水印算法的鲁棒性和嵌入容量。
Description
技术领域
本发明属于图像信号处理领域,具体涉及一种基于多变换域的彩色图像数字水印方法。
背景技术
随着互联网技术的迅猛发展,数字多媒体数据传输越发迅速、便捷,引发越益严重的个人版权侵害问题。因此,加强对数字作品的个人版权保护问题的需要,其中,数字水印技术作为一种有效的版权保护方式也备受重视。数字水印技术是将版权信息直接嵌入到载图中,需要验证版权时,直接通过与嵌入算法相应的提取算法提取出版权信息,以达到对作品版权保护的目的。
灰度图像数字水印算法主要分为空域算法和变换域算法两类。其中,空域算法是最早的数字水印算法,主要有LSB、Patchwork以及纹理映射等方法。此类算法嵌入容量虽大,但通过改变图像的灰度值嵌入水印,导致算法的鲁棒性较差。因此,为提高算法的鲁棒性,研究人员提出变换算法,例如离散小波变换(Discrete Wavelet Transform,DWT)、四元数离散余弦变换(Discrete Cosine Transform,DCT)和奇异值分解(Singular ValueDecomposition,SVD)等。此类算法是将水印嵌入到某些频率值,一般具有较好的鲁棒性和透明性。但是不同变换域具有各自的优势和缺点,例如DWT水印算法具有良好抗JPEG压缩和抗剪切能力,透明性良好,但是不能抵抗几何攻击能;DCT水印算法则对低通滤波和有损压缩能力,但透明性较差;SVD水印算法因其具有良好的几何不变性,也备受关注。为获得更好的鲁棒性和透明性,研究人员提出了双变换域水印算法的思想,即将各种变换域方法的优势结合,缺陷弥补,例如DWT_DCT、DCT_SVD、DWT_DFT等。
目前,针对图像数字水印算法多是针对于灰度图像,而彩色图像的应用日益广泛,为切合实际需要,彩色图像数字水印算法的研究意义加大。现有彩色图像水印算法,主要分为三大类:1)基于颜色空间的单通道处理方式,将水印嵌入到某一个色彩通道,实质还是灰度图像处理方式,没有最大化利用彩色图像嵌入空间;2)基于色彩空间的多通道处理方式,将水印分散到多个通道进行嵌入,这种方式将各个通道独立处理,嵌入空间较大,但没有考虑各个通道间的联系,且计算复杂度较高;3)基于四元数的整体处理方式,利用一个四元数描述多个色彩通道信息,水印通过调节四元数嵌入到载体当中,这种整体处理方式,保留各个通道间的联系的同时降低计算复杂度,并且具有良好的透明性,但在嵌入容量方面还有待提升。综上所述,现有针对彩色图像水印算法存在一些问题:2)处理方式比较单一,没有将分通道处理和整体处理方式结合,不能充分利用彩色图像巨大的嵌入空间问题。2)在平衡鲁棒性、透明性以及嵌入容量方面还有提升空间。
发明内容
本发明针对现有技术中存在的上述问题,结合DWT、四元数离散余弦变换(Quaternion Discrete Cosine Transform,QDCT)和奇异值变换提出一种混合变换域数字水印嵌入方法。通过分通道小波变换和四元数离散余弦变换相结合的方法,利用小波变换的多尺度频率分解特性和四元数余弦变换整体式计算的特点,将数字水印嵌入其中,该方法有较好的鲁棒性和透明性,有效提高嵌入容量。
本发明解决上述技术问题的技术方案是提出一种多变换域的彩色图像数字水印处理方法,包括水印嵌入和水印提取两个阶段,水印嵌入阶段包括:将彩色图像分通道处理后对各个色彩分量进行DWT,提取各色彩分量上同类频率子带,将其描述为一个四元数矩阵,根据水印尺寸对四元数矩阵进行分块,再将每个分块QDCT,提取QDCT的实数矩阵的第一个分量,组成一个与水印尺寸大小一致的系数矩阵,将水印嵌入到系数矩阵中。水印提取阶段包括:将彩色图像进行分通道的DWT,将嵌入阶段选定的各色彩分量的频率子带描述为一个四元数矩阵,将四元数矩阵进行分块并进行QDCT得到子块,将各子块的实数部分的第一个分量提取出组成一个系数矩阵,根据嵌入密钥提取数字水印。
提取各色彩分量上同类频率子带进一步包括:读取RGB彩色图像F,分通道处理,得红色分量(R)、绿色分量(G)、蓝色分量(B),将各个通道进行一级DWT获得四个频率子带,分别为低频(LLq)、水平(LHq)、垂直(HLq)、对角线(HHq)子带。其中,q代表所属色彩分量,即q∈{R,G,B}。例如,R分量对应的低频为LLR,G分量对应的低频为LLG,B分量对应的低频为LLB。
将水印嵌入到系数矩阵中进一步包括:先将水印图像W用猫映射(Arnold)置乱得加密水印W1,其置乱次数为密钥key1。再根据公式E=USVT将系数矩阵进行奇异值分解,提取对角矩阵S和两个酉矩阵U与V,将水印根据公式E′=S+aW1嵌入对角矩阵S中得过渡矩阵E′,根据公式E′=U2S2V2 T将过渡水印E′分解提取对角矩阵S2,最后根据公式E″=US2VT得到重构的系数矩阵E″。其中,S、U2、V2依次保存作嵌入密钥key2、key3、key4,a为控制嵌入强度参数,T表示转置。
水印的提取中将各子块的实数部分的第一个分量提取出组成一个系数矩阵具体包括,获取含水印的彩色图像F*,分通道处理得R、G、B分量。将各个通道进行一级DWT获得四个频率子带,分别为低频(LLq *)、水平(LHq *)、垂直(HLq *)、对角线(HHq *)子带,将水印嵌入区域的三个色彩分量的频率子带用四元数描述为一个四元数矩阵再将四元数矩阵进行大小为n×n的分块,共分为l×l块子块。对每个子块进行QDCT,提取坐标为(o,p)的子块的实部矩阵Drel* O,P的第一个分量Drel* O,P(1,1),用其构成系数矩阵E*,
即E*(o,p)=Drel* O,P(1,1),
一个四元数可表示为为Q=Q1+Q2i+Q3j+Q4k,其中,Q1、Q2、Q3和Q4为任意实数,i,j,k为四元数的虚数单位。本发明中将任意彩色图像的各个色彩分量的同类频率子带用四元数描述为一个纯四元数,即实数为0,虚数部分用于表征各色分量的频率子带。
低频子带四元数矩阵为LL=LLRi+LLGj+LLBk,水平子带四元数矩阵LH=LHRi+LHGj+LHBk,垂直子带四元数矩阵为HL=HLRi+HLGj+HLBk,高频子带四元数矩阵为HH=HHRi+HHGj+HHBk,选择LL、LH、HL、HH四个矩阵中的一个矩阵对作为水印待嵌入区域矩阵Qemb。
根据公式确定矩阵Qemb分块尺寸为n×n,其中,彩色图像F大小M×N,u=min(M,N),ceil表示正无穷方向取最小整数,按照分块位置(o,p)标记子块为Qemb o,p,对子块Qemb o,p进行QDCT得四元数矩阵DO,P,DO,P=Drel O,P+Dimg1 O,Pi+Dimg2 O,Pj+Dimg3 O,Pk,其中,Dimg1 O,P、Dimg2 O,P、Dimg3 O,P为虚部矩阵,Drel O,P为实部矩阵,且提取所有子块实部矩阵的第一个分量Drel O,P(1,1)构造一个系数矩阵E,即该系数矩阵为
根据嵌入密钥提取数字水印进一步包括,按照水印嵌入阶段的方法提取的系数矩阵E*,并将E*进行奇异值分解,提取对角阵S*,计算过渡矩阵E**=U2S*V2 T,根据公式W*=(E**-S)/a提取加密水印W*,其中,酉矩阵U2、V2分别为嵌入过程中的密钥key3、key4,S为嵌入过程中的密钥key2,a嵌入过程中的控制嵌入强度参数。将W*进行Arnold映射,其置乱次数为密钥key1,获得提取的水印W**。
本发明还提出一种多变换域的彩色图像数字水印水印嵌入方法,将彩色图像,分通道处理后对各个色彩分量进行DWT,提取各色彩分量上同类频率子带,将其描述为一个四元数矩阵,根据水印尺寸对四元数矩阵进行分块,再将每个分块QDCT,提取QDCT的实数矩阵的第一个分量,组成一个与水印尺寸大小一致的系数矩阵,将水印嵌入到系数矩阵中。
本发明还提出一种多变换域的彩色图像数字水印提取方法,将彩色图像进行分通道的DWT,将嵌入阶段选定的各色彩分量的频率子带描述为一个四元数矩阵,将四元数矩阵进行分块并进行离散余弦变换QDCT得到子块,将各子块的实数部分的第一个分量提取出组成一个系数矩阵,根据嵌入密钥提取数字水印。
本发明将RGB彩色图像先分通道DWT处理后采用QDCT的整体处理方式,前者提高水印系统鲁棒性,后者保留彩色图像各个色彩分量间的关联性,提高透明性的同时简化了多通道处理存在的计算复杂度。
附图说明
图1为本发明的水印嵌入流程图;
图2为本发明的水印提取流程图。
具体实施方式
为了使本发明的技术方案更加清楚明白,将结合附图和实施例进行进一步阐明本发明的具体内容,但本发明的内容不仅仅局限于下面的实施例。本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样在本申请所列权利要求书限定范围之内。
本发明提出一种基于多变换域的彩色图像数字水印处理方法。主要包括水印嵌入和水印提取两个部分。
水印嵌入部分包括:对彩色图像分通道置乱后作DWT,将同类频率子带用四元数描述;根据水印信息大小,确定所需嵌入的频率子带,将该频率子带对应的四元数频率矩阵进行分块及DQCT得到子块;提取所有子块四元数余弦系数的实数部分的第一个分量,构造一个系数矩阵;将水印奇异值分解(SVD)的特征值嵌入到系数矩阵中,最后进行逆变换得到含水印的彩色图像。提取部分包括:将含水印彩色图像进行分通道的DWT,将嵌入阶段选定的各色彩分量的频率子带描述为一个四元数矩阵;然后,将四元数矩阵进行分块DQCT,再将各子块的实数部分的第一个分量提取出组成一个系数矩阵;最后根据嵌入阶段的嵌入方法对应的提取方法,提取数字水印。
水印嵌入部分具体包括:读取大小为l×l的灰度水印图像W,进行Arnold映射,得置乱水印W1,记置乱次数AW为密钥key1。读取大小为M×N的RGB彩色图像F,先将R、G、B三个色彩分量分进行一次小波分解,分别获取各分量的低频(LLq)、水平(LHq)、垂直(HLq)、对角线(HHq)子带。四元数描述公式为Q=Q1+Q2i+Q3j+Q4k,其中,Q1、Q2、Q3和Q4为任意实数,i,j,k为四元数的虚数单位。本发明中将任意彩色图像的各个色彩分量的同类频率子带用四元数描述为一个纯四元数,即实数为0,虚数部分用于表征各色分量的频率子带。则R、G、B三个色彩分量的低频(LLR,LLG,LLB)、水平(LHR,LHG,LHB)、垂直(HLR,HLG,HLB)、对角线(HHR,HHG,HHB)子带,获得低频子带四元数矩阵为LL=LLRi+LLGj+LLBk,水平子带四元数矩阵LH=LHRi+LHGj+LHBk,垂直子带四元数矩阵为HL=HLRi+HLGj+HLBk,高频子带四元数矩阵为HH=HHRi+HHGj+HHBk。
选择LL、LH、HL、HH四个矩阵中的一个矩阵对作为水印待嵌入区域矩阵,记该待嵌入区域矩阵为Qemb,Qemb∈{LL,LH,HL,HH}。设水印尺寸l×l,彩色图像F大小M×N,根据公式确定矩阵Qemb分块尺寸为n×n,共计l×l块,其中,u=min(M,N),ceil表示正无穷方向取最小整数。
u=min(M,N)。按照分块位置标记子块为Qemb o,p,(o,p)为子块分块位置,1≤o,p≤l。然后,对子块Qemb o,p进行DQCT得四元数DO,P,DO,P=Drel O,P+Dimg1 O,Pi+Dimg2 O,Pj+Dimg3 O,Pk,其中,Dimg1 O,P、Dimg2 O,P、Dimg3 O,P为虚部矩阵,Drel O,P为实部矩阵,且提取所有子块实部矩阵的第一个分量Drel O,P(1,1)构造一个新的系数矩阵E,其中E(o,p)=Drel O,P(1,1),即该系数矩阵为
将水印W嵌入系数矩阵E中。更进一步地,为有更好的安全性和透明性在具体实施时可将水印W嵌入到矩阵E奇异值分解的对角矩阵中。具体包括:首先将水印图像W用猫映射(Arnold)置乱得加密水印W1,其置乱次数为密钥key1将。然后根据奇异值分解公式E=USVT将系数矩阵进行奇异值分解,提取对角矩阵S和两个酉矩阵U与V,其中T表示转置,再将水印嵌入到对角矩阵S中得过渡矩阵E′,即E′=S+aW1,此时置乱水印W1通过控制嵌入强度参数a嵌入到系数矩阵的对角矩阵中。之后根据公式E′=U2S2V2 T将过渡矩阵E′分解提取对角矩阵S2。最后根据公式E″=US2VT得到重构系数矩阵E″。其中,S、U2、V2依次保存作嵌入密钥key2、key3、key4。
令Drel O,P(1,1)=E″(O,P),将每个子块对应的实数矩阵Drel O,P的第一个分量用重构系数矩阵E″(O,P)中对应的系数值替代,得到新的实数矩阵Drel O,P′。将新的实数矩阵与未经更改的虚部系数Dimg1 O,P、Dimg2 O,P和Dimg3 O,P重新用四元数表示成一个四元数矩阵DO,P′=Drel O,P′+Dimg1 O,Pi+Dimg2 O,Pj+Dimg2 O,Pk,再将DO,P′进行四元数离散余弦的逆变换,得各子块嵌入水印后的频率子带四元数矩阵Qemb o,p′;将所有子块得到的频率子带四元数矩阵按原有位置组合,即可得嵌入水印后的频率子带四元数矩阵Qemb′,该四元数矩阵的三个虚部矩阵分别为红色分量、绿色分量和蓝色分量含水印信息的频率子带;最后将同一色彩分量进行小波重构,得到各自含水印信息的色彩分量,最后进行通道合并,即可得含水印的彩色图像F*。
以下以一具体实例说明本发明水印方法的嵌入步骤。如图1所示为本发明的水印嵌入流程图。先读取大小为64×64的灰度水印图像W,进行猫映射(Arnold)置乱水印为W1,记置乱次数32为密钥key1;然后读取大小为512×512的RGB彩色图像F,分离R、G、B三个色彩分量,将每个色彩分量进行一次小波分解,分别获取各分量的四个频率子带:低频子带(LLq)、水平子带(LHq)、垂直子带(HLq)、对角线子带(HHq),q代表频率子带所属色彩分量,即q∈{R,G,B};其次,根据四元数描述公式:Q=Q1+Q2i+Q3j+Q4k将在不同色彩分量上的同类频率子带用四元数描述为一个纯四元数,获得低频子带四元数矩阵为LL,水平子带四元数矩阵LH,垂直子带四元数矩阵为HL,高频子带四元数矩阵为HH;考虑子带兼具透明性和鲁棒性特性,以下选择垂直子带四元数矩阵HL作为待嵌入水印区域为例进一步具体说明水印嵌入过程。
将HL进行分块,由公式确定每子块分块尺寸为n×n,其中ceil为正无穷方向取最小整数函数,u为彩色图像F的最小边长,l为水印的边长,在此实例中u=512,l=64,因此将HL以4×4的分块大小分为64×64块。同时,标记每个字块的坐标为(o,p),1≤o,p≤64。对每个小块进行DQCT,得一个实部矩阵和三个虚部矩阵,三个虚部矩阵分别为Dimg1 O,P、Dimg2 O,P、Dimg3 O,P,记坐标为(o,p)子块对应的实部矩阵为提取每个子块实数部分的第一个分量DO,P(1,1),根据其子块所在位置,组成一个大小为64×64的系数矩阵E,其中E(o,p)=Drel O,P(1,1),即该系数为将加密水印W1嵌入系数矩阵E中,但为有更好的安全性和透明性在具体实施时可将W1嵌入到E奇异值分解的对角矩阵中。具体步骤包括:首先将系数矩阵进行奇异值分解E=USVT,提取对角矩阵S和两个酉矩阵U与V,再用加法原则将水印嵌入到对角矩阵S得过渡矩阵E′,即E′=S+aW1,此时水印W1通过控制嵌入强度参数a嵌入到系数矩阵的对角矩阵中。然后,根据公式E′=U2S2V2 T将过渡矩阵E′分解,提取对角矩阵S2;最后根据公式E″=US2VT得到重构系数矩阵E″。其中,将S、U2、V2依次保存作嵌入密钥key2、key3、key4。
令Drel O,P(1,1)=E″(O,P),将坐标为(O,P)的子块的实数矩阵Drel O,P的第一个分量Drel O,P(1,1)用重构系数矩阵E″中(O,P)处的系数值替代,得到新的实数矩阵Drel O,P′。将新的实部矩阵与未经更改的虚部系数Dimg1 O,P、Dimg2 O,P和Dimg3 O,P重新用四元数表示成一个四元数矩阵DO,P=Drel O,P+Dimg1 O,Pi+Dimg2 O,Pj+Dimg2 O,Pk,再将DO,P进行四元数离散余弦的逆变换,得各子块嵌入水印后的频率子带四元数矩阵。将所有子块得到的频率子带四元数矩阵按原有位置组合,即可得嵌入水印后的频率子带四元数矩阵HL′,HL′的三个虚部矩阵分别为红色分量、绿色分量和蓝色分量含水印信息的垂直子带。最后将同一色彩分量进行小波重构,得到各自含水印信息的色彩分量,进行通道合并,即可得含水印的彩色图像F*。
以下以具体的实例说明水印的提取步骤。图2为本发明的水印提取流程图。
获取含水印的彩色图像F*,分通道处理,得红色分量(R)、绿色分量(G)、蓝色分量(B),再将各个通道进行一级DWT获得四个频率子带,分别为低频(LLq *)、水平(LHq *)、垂直(HLq *)、对角线(HHq *)子带,q代表频率子带所属色彩分量,记q∈{R,G,B}。如将垂直子带作为水印待嵌入区域为例进行说明,三个色彩分量的垂直子带用四元数描述为一个四元数矩阵,则垂直子带四元数矩阵为HL=HLRi+HLGj+HLBk。再将四元数矩阵HL进行4×4大小分块,共计64×64块,将每个子块坐标记为(o,p),1≤o,p≤l。对每个小块进行DQCT,提取坐标为(o,p)的子块的实部矩阵Drel* O,P的第一个分量Drel* O,P(1,1),用其构成系数矩阵即E*(o,p)=Drel* O,P(1,1)。
对E*进行奇异值分解(SVD)分解,提取对角阵S*。计算过渡矩阵E**=U2S*V2 T,其中酉矩阵U2、V2分别为嵌入过程中的密钥key3、key4。根据公式W*=(E**-S)/a计算提取加密水印W*,S为嵌入过程中的密钥key2,a为嵌入过程中的控制嵌入强度参数,将W*进行Arnold映射,其置乱次数为32,获得提取的水印W**。
目前,针对数字图像水印方法性能优劣主要判别两个指标:透明性和鲁棒性。透明性度量常用峰值信噪比(PSNR)方法计算,主要用于判断含水印载图与原始图像的失真情况,PSNR值越大,则说明该水印算法的透明性越高。尺寸为M×N的RGB彩色图像的PSNR定义如下:
式中:F(i,j)为RGB彩色图像的分量像素值;F′(i,j)为嵌入水印后RGB彩色图像的分量像素值。
除了对透明性进行度量,还需要对评价算法的鲁棒性,一般采用计算检测水印和初始水印的归一化相关系数NC来进行评价,取值0到1,NC值越高,说明算法的鲁棒性越强。尺寸为l×l水印的NC值定义如下:
式中,W(i,j)为原始水印,W′(i,j)为提取水印。
本方案使用的实例选用尺寸大小为512×512标准彩色BMP真彩图像做载图,包括经典实验图像Airplane、Lena、Baboon和Pepper,选用64×64大小的灰度图像Lena作水印图像。均在MATLAB2014a软件平台下进行仿真,这些源数据被绝大多数研究使用,实验结果具有代表性。
参数a是强度控制参数,其取值将直接影响水印算法的透明性和鲁棒性。表1为不同参数a下,从未受攻击的水印载图中提取的水印与原始水印的归一化相关系数(NC)。如表所示随着a的取值从0.1到1的增加,NC值也逐渐增加,虽然a≤0.3时,NC值不为1,但NC值在大于0.75的情况下,有意义水印的提取可视为有效,依旧可以实现版权认证,因此本发明算法是有效算法。
表1未受攻击下提取水印与原始水印的NC值
表2为不同参数a下,不同载图对应的峰值信噪比(PSNR),随着a的取值从0.1到1的增加,各载图对应的峰值信噪比(PSNR)也逐渐降低,但PSNR取值均在35以上,说明本发明具有良好的透明性。
表2含水印载图与载图的PSNR值
为验证本发明的抗攻击能力,将对含水印载图进行多种攻击测试,包括椒盐噪声、高斯噪声、中值滤波、对比度攻击、JPEG压缩、剪切和旋转。通过计算受攻击含水印载图中提取的水印和原始水印的NC值判断算法的鲁棒性。如表3所示,在参数a=0.5时,不同载图应对不同攻击下,提取的水印与原始水印的NC值多数都接近1。本发明算法在抗击噪声、中值滤波、JPEG压缩以及对比度攻击的普通信号处理攻击具有良好的鲁棒性,在抗击剪切和旋转这种去同步攻击,NC值接近0.9,仍然具有一定的抗攻击性。
表3攻击测试下提取水印与原水印的NC值
综上所述,本发明采用双变换域的思想,结合DWT、DQCT和奇异值分解(SVD)变换方法,即将分通道处理和整体处理相结合,简化分通道的双变换域算法的计算复杂度,同时因四元数的引入,保持各个通道间的关联,提升水印算法的透明性。根据实验表明本发明具有良好抗一般信号处理攻击能力,同时,具有一定抗几何攻击能力。目前针对彩色图像的水印算法较少,本发明方法可为彩色图像水印算法的研究提供一种新思路。
Claims (15)
1.一种多变换域的彩色图像数字水印处理方法,包括水印嵌入和水印提取两个阶段,其特征在于,将彩色图像分通道处理后,对各个色彩分量进行DWT,提取各色彩分量上同类频率子带,将其描述为一个四元数矩阵,根据水印尺寸对四元数矩阵进行分块,再将每个分块QDCT并提取实数矩阵的第一个分量,组成一个与水印尺寸大小一致的系数矩阵,将水印置乱后嵌入到系数矩阵中;水印提取阶段包括:将彩色图像进行分通道的离散小波变换,将嵌入阶段选定的各色彩分量的频率子带描述为一个四元数矩阵,将四元数矩阵进行分块并进行离散余弦变换得到子块,将各子块的实数部分的第一个分量提取出组成一个系数矩阵,根据嵌入密钥提取数字水印。
2.根据权利要求1所述的方法,其特征在于,提取各色彩分量上同类频率子带进一步包括:读取RGB彩色图像F,分通道处理,得红色分量(R)、绿色分量(G)、蓝色分量(B),将各个通道进行一次DWT获得四个频率子带,分别为低频(LLq)、水平(LHq)、垂直(HLq)、对角线(HHq)子带,其中,q代表频率子带所属色彩分量,即q∈{R,G,B}。
3.根据权利要求1所述的方法,其特征在于,将水印置乱后嵌入到系数矩阵中进一步包括:先将水印W用猫映射(Arnold)置乱得置乱水印W1,记置乱次数为密钥key1。根据公式E=USVT将系数矩阵进行奇异值分解,提取对角矩阵S和两个酉矩阵U与V,将水印W1根据公式E′=S+aW1嵌入对角矩阵S中得过渡矩阵E′,根据公式E′=U2S2V2 T将过渡水印E′分解提取对角矩阵S2,最后根据公式E″=US2VT得到重构系数矩阵E″,其中,S、U2、V2依次保存作嵌入密钥key1、key2、key3,a为控制嵌入强度参数,T表示转置。
4.根据权利要求1所述的方法,其特征在于,水印的提取中将各子块的实数部分的第一个分量提取出组成一个系数矩阵具体包括,获取含水印的彩色图像F*,分通道处理得RGB分量,将各个通道进行一次离散小波变换获得四个频率子带,分别为低频(LLq *)、水平(LHq *)、垂直(HLq *)、对角线(HHq *)子带,q代表频率子带所属色彩分量,记q∈{R,G,B},将水印嵌入区域的三个色彩分量的频率子带用四元数描述为一个四元数矩阵再将四元数矩阵进行n×n大小分块,对每个子块进行DQCT,提取坐标为(o,p)的子块的实部矩阵的第一个分量用其构成系数矩阵E*,
即
5.根据权利要求2所述的方法,其特征在于,根据四元数描述公式Q=Q1+Q2i+Q3j+Q4k将在不同色彩分量上的同类频率子带描述为一个纯四元数,Q1为四元数的实数分量,Q2、Q3和Q4为四元数的虚数分量,其中,虚数分量用于表征各色分量的频率子带,i,j,k为四元数的虚数单位。
6.根据权利要求2所述的方法,其特征在于,低频子带四元数矩阵为LL=LLRi+LLGj+LLBk,水平子带四元数矩阵LH=LHRi+LHGj+LHBk,垂直子带四元数矩阵为HL=HLRi+HLGj+HLBk,高频子带四元数矩阵为HH=HHRi+HHGj+HHBk,选择LL、LH、HL、HH四个矩阵中的一个矩阵对作为水印待嵌入区域矩阵Qemb。
7.根据权利要求5所述的方法,其特征在于,根据公式确定矩阵Qemb分块尺寸为n×n,其中,彩色图像F大小M×N,u=min(M,N),ceil为正无穷方向取最小整数函数,按照分块位置(o,p)标记子块为Qemb o,p,对子块Qemb o,p进行QDCT得DO,P,DO,P=Drel O,P+Dimg1 O,Pi+Dimg2 O,Pj+Dimg3 O,Pk,其中,Dimg1 O,P、Dimg2 O,P、Dimg3 O,P为虚部矩阵,Drel O,P为实部矩阵,且提取所有子块实部矩阵的第一个分量Drel O,P(1,1)构造一个新的系数矩阵E,即该系数矩阵为
8.根据权利要求4所述的方法,其特征在于,根据嵌入密钥提取数字水印进一步包括,对E*进行奇异值分解并提取对角阵S*,计算过渡矩阵E**=U2S*V2 T,根据公式W*=(E**-S)/a计算加密水印W*,其中,酉矩阵U2、V2分别为嵌入过程中的密钥key3、key4,S为嵌入过程中的密钥key2,a嵌入过程中的控制嵌入强度参数,利用密钥key1将W*进行Arnold置乱,获得提取的水印W**。
9.一种多变换域的彩色图像数字水印水印嵌入方法,其特征在于,将彩色图像,分通道处理后对各个色彩分量进行DWT,提取各色彩分量上同类频率子带,将其描述为一个四元数矩阵,根据水印尺寸对四元数矩阵进行分块,再将每个分块QDCT,提取QDCT的实数矩阵的第一个分量,组成一个与水印尺寸大小一致的系数矩阵,将水印嵌入到系数矩阵中。
10.一种多变换域的彩色图像数字水印提取方法,其特征在于,将彩色图像进行分通道的DWT,将嵌入阶段选定的各色彩分量的频率子带描述为一个四元数矩阵,将四元数矩阵进行分块并进行离散余弦变换QDCT得到子块,将各子块的实数部分的第一个分量提取出组成一个系数矩阵,根据嵌入密钥提取数字水印。
11.根据权利要求9所述的方法,其特征在于,将水印置乱后嵌入到系数矩阵中进一步包括:先将水印W用猫映射(Arnold)置乱得置乱水印W1,记置乱次数为密钥key1。根据公式E=USVT将系数矩阵进行奇异值分解,提取对角矩阵S和两个酉矩阵U与V,将水印W1根据公式E′=S+aW1嵌入对角矩阵S中得过渡矩阵E′,根据公式E′=U2S2V2 T将过渡水印E′分解提取对角矩阵S2,最后根据公式E″=US2VT得到重构系数矩阵E″,其中,S、U2、V2依次保存作嵌入密钥key1、key2、key3,a为控制嵌入强度参数,T表示转置。
12.根据权利要求9所述的方法,其特征在于,根据四元数描述公式Q=Q1+Q2i+Q3j+Q4k将在不同色彩分量上的同类频率子带描述为一个纯四元数,Q1为四元数的实数分量,Q2、Q3和Q4为四元数的虚数分量,其中,虚数分量用于表征各色分量的频率子带,i,j,k为四元数的虚数单位,低频子带四元数矩阵为LL=LLRi+LLGj+LLBk,水平子带四元数矩阵LH=LHRi+LHGj+LHBk,垂直子带四元数矩阵为HL=HLRi+HLGj+HLBk,高频子带四元数矩阵为HH=HHRi+HHGj+HHBk,选择LL、LH、HL、HH四个矩阵中的一个矩阵对作为水印待嵌入区域矩阵Qemb。
13.根据权利要求12所述的方法,其特征在于,根据公式确定矩阵Qemb分块尺寸为n×n,其中,彩色图像F大小M×N,u=min(M,N),ceil为正无穷方向取最小整数函数,按照分块位置(o,p)标记子块为Qemb o,p,对子块Qemb o,p进行QDCT得DO,P,DO,P=Drel O,P+Dimg1 O,Pi+Dimg2 O,Pj+Dimg3 O,Pk,其中,Dimg1 O,P、Dimg2 O,P、Dimg3 O,P为虚部矩阵,Drel O,P为实部矩阵,且提取所有子块实部矩阵的第一个分量Drel O,P(1,1)构造一个新的系数矩阵E,即该系数矩阵为
14.根据权利要求10所述的方法,其特征在于,水印的提取中将各子块的实数部分的第一个分量提取出组成一个系数矩阵具体包括,获取含水印的彩色图像F*,分通道处理得RGB分量,将各个通道进行一次离散小波变换获得四个频率子带,分别为低频(LLq *)、水平(LHq *)、垂直(HLq *)、对角线(HHq *)子带,q代表频率子带所属色彩分量,记q∈{R,G,B},将水印嵌入区域的三个色彩分量的频率子带用四元数描述为一个四元数矩阵再将四元数矩阵进行n×n大小分块,对每个子块进行DQCT,提取坐标为(o,p)的子块的实部矩阵的第一个分量用其构成系数矩阵E*,即
15.根据权利要求14所述的方法,其特征在于,根据嵌入密钥提取数字水印进一步包括,对E*进行奇异值分解并提取对角阵S*,计算过渡矩阵E**=U2S*V2 T,根据公式W*=(E**-S)/a计算加密水印W*,其中,酉矩阵U2、V2分别为嵌入过程中的密钥key3、key4,S为嵌入过程中的密钥key2,a嵌入过程中的控制嵌入强度参数,利用密钥key1将W*进行Arnold置乱,最终获得提取的水印W**。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711401354.7A CN108564519B (zh) | 2017-12-22 | 2017-12-22 | 一种基于多变换域的彩色图像数字水印方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711401354.7A CN108564519B (zh) | 2017-12-22 | 2017-12-22 | 一种基于多变换域的彩色图像数字水印方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108564519A true CN108564519A (zh) | 2018-09-21 |
CN108564519B CN108564519B (zh) | 2021-11-16 |
Family
ID=63529437
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711401354.7A Active CN108564519B (zh) | 2017-12-22 | 2017-12-22 | 一种基于多变换域的彩色图像数字水印方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108564519B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109544434A (zh) * | 2018-10-16 | 2019-03-29 | 江苏信实云安全技术有限公司 | 用于raw格式图像的cwt-qim数字盲水印算法 |
CN110097490A (zh) * | 2019-04-26 | 2019-08-06 | 上海理工大学 | 基于IWT-Schur的感兴趣区域可逆水印嵌入及提取方法 |
CN110163266A (zh) * | 2019-05-08 | 2019-08-23 | 南京邮电大学 | 基于四元数矩阵奇异值分解的图像哈希算法及其应用 |
CN110189243A (zh) * | 2019-05-13 | 2019-08-30 | 杭州电子科技大学上虞科学与工程研究院有限公司 | 一种基于张量奇异值分解的彩色图像鲁棒水印方法 |
CN110992236A (zh) * | 2019-11-28 | 2020-04-10 | 陕西科技大学 | 数字水印嵌入环境确定方法、装置、设备及可读存储介质 |
CN111242827A (zh) * | 2020-01-15 | 2020-06-05 | 复旦大学 | 基于dt-cwt和svd变换的鲁棒彩色数字图像水印方法 |
CN111583085A (zh) * | 2020-04-16 | 2020-08-25 | 天津大学 | 一种基于变换域的彩色图像数字水印方法 |
CN111598197A (zh) * | 2020-05-18 | 2020-08-28 | 燕山大学 | 一种基于小波与幻方变换的qr码加密、解密方法及系统 |
CN111654592A (zh) * | 2020-06-05 | 2020-09-11 | 中国矿业大学 | 基于秘密分享的多图像隐藏方法 |
CN111681153A (zh) * | 2020-06-08 | 2020-09-18 | 重庆邮电大学 | 基于非下采样轮廓波变换的cs-svd的鲁棒数字水印算法 |
CN112149062A (zh) * | 2019-06-26 | 2020-12-29 | 中电万维信息技术有限责任公司 | 基于水印的电子印章生成方法、装置、系统及存储介质 |
CN112529756A (zh) * | 2020-12-01 | 2021-03-19 | 湖南长城信息金融设备有限责任公司 | 一种数字图像水印嵌入方法、提取方法、装置、存储介质 |
CN113822786A (zh) * | 2021-11-25 | 2021-12-21 | 成都盛思睿信息技术有限公司 | 一种强鲁棒性的隐形数字图像水印处理方法及装置 |
CN112529756B (zh) * | 2020-12-01 | 2024-05-31 | 长城信息股份有限公司 | 一种数字图像水印嵌入方法、提取方法、装置、存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030095682A1 (en) * | 2001-11-20 | 2003-05-22 | Sanghyun Joo | Apparatus and method for embedding and extracting digital watermarks based on wavelets |
US20030099373A1 (en) * | 2001-11-27 | 2003-05-29 | Sanghyun Joo | Robust blind watermarking method in wavelet DC components |
EP1316917A2 (en) * | 2001-12-03 | 2003-06-04 | Samsung Electronics Co., Ltd. | Apparatus and method for embedding watermark into original information, transmitting watermarked information, and reconstructing the watermark |
CN101042769A (zh) * | 2007-01-12 | 2007-09-26 | 中国人民解放军国防科学技术大学 | 一种基于小波与dct双重域的主动式数字图像内容鉴别方法 |
CN102236882A (zh) * | 2011-07-21 | 2011-11-09 | 上海理工大学 | 基于版权保护的变换域加密全息水印算法 |
CN102496135A (zh) * | 2011-12-06 | 2012-06-13 | 银江股份有限公司 | 一种基于dwt域的数字水印方法及系统 |
CN106447590A (zh) * | 2016-10-11 | 2017-02-22 | 华南理工大学 | 数字图像中水印的加入及提取方法 |
-
2017
- 2017-12-22 CN CN201711401354.7A patent/CN108564519B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030095682A1 (en) * | 2001-11-20 | 2003-05-22 | Sanghyun Joo | Apparatus and method for embedding and extracting digital watermarks based on wavelets |
US20030099373A1 (en) * | 2001-11-27 | 2003-05-29 | Sanghyun Joo | Robust blind watermarking method in wavelet DC components |
EP1316917A2 (en) * | 2001-12-03 | 2003-06-04 | Samsung Electronics Co., Ltd. | Apparatus and method for embedding watermark into original information, transmitting watermarked information, and reconstructing the watermark |
CN101042769A (zh) * | 2007-01-12 | 2007-09-26 | 中国人民解放军国防科学技术大学 | 一种基于小波与dct双重域的主动式数字图像内容鉴别方法 |
CN102236882A (zh) * | 2011-07-21 | 2011-11-09 | 上海理工大学 | 基于版权保护的变换域加密全息水印算法 |
CN102496135A (zh) * | 2011-12-06 | 2012-06-13 | 银江股份有限公司 | 一种基于dwt域的数字水印方法及系统 |
CN106447590A (zh) * | 2016-10-11 | 2017-02-22 | 华南理工大学 | 数字图像中水印的加入及提取方法 |
Non-Patent Citations (2)
Title |
---|
KAUSHIK DEB ET AL.: "Combined DWT-DCT Based Digital Image Watermarking Technique for Copyright Protection", 《2012 7TH INTERNATIONAL CONFERENCE ON ELECTRICAL AND COMPUTER ENGINEERING》 * |
孔玲君 等: "一种基于DWT-DCT变换域的全息水印技术", 《光电子·激光》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109544434A (zh) * | 2018-10-16 | 2019-03-29 | 江苏信实云安全技术有限公司 | 用于raw格式图像的cwt-qim数字盲水印算法 |
CN109544434B (zh) * | 2018-10-16 | 2023-06-09 | 江苏水印科技有限公司 | 用于raw格式图像的cwt-qim数字盲水印算法 |
CN110097490A (zh) * | 2019-04-26 | 2019-08-06 | 上海理工大学 | 基于IWT-Schur的感兴趣区域可逆水印嵌入及提取方法 |
CN110163266A (zh) * | 2019-05-08 | 2019-08-23 | 南京邮电大学 | 基于四元数矩阵奇异值分解的图像哈希算法及其应用 |
CN110163266B (zh) * | 2019-05-08 | 2022-09-06 | 南京邮电大学 | 基于四元数矩阵奇异值分解的图像哈希算法及其应用 |
CN110189243A (zh) * | 2019-05-13 | 2019-08-30 | 杭州电子科技大学上虞科学与工程研究院有限公司 | 一种基于张量奇异值分解的彩色图像鲁棒水印方法 |
CN110189243B (zh) * | 2019-05-13 | 2023-03-24 | 杭州电子科技大学上虞科学与工程研究院有限公司 | 一种基于张量奇异值分解的彩色图像鲁棒水印方法 |
CN112149062A (zh) * | 2019-06-26 | 2020-12-29 | 中电万维信息技术有限责任公司 | 基于水印的电子印章生成方法、装置、系统及存储介质 |
CN110992236A (zh) * | 2019-11-28 | 2020-04-10 | 陕西科技大学 | 数字水印嵌入环境确定方法、装置、设备及可读存储介质 |
CN110992236B (zh) * | 2019-11-28 | 2023-05-26 | 北京奇良海德印刷股份有限公司 | 数字水印嵌入环境确定方法、装置、设备及可读存储介质 |
CN111242827A (zh) * | 2020-01-15 | 2020-06-05 | 复旦大学 | 基于dt-cwt和svd变换的鲁棒彩色数字图像水印方法 |
CN111242827B (zh) * | 2020-01-15 | 2022-12-20 | 复旦大学 | 基于dt-cwt和svd变换的鲁棒彩色数字图像水印方法 |
CN111583085A (zh) * | 2020-04-16 | 2020-08-25 | 天津大学 | 一种基于变换域的彩色图像数字水印方法 |
CN111583085B (zh) * | 2020-04-16 | 2023-07-04 | 天津大学 | 一种基于变换域的彩色图像数字水印方法 |
CN111598197A (zh) * | 2020-05-18 | 2020-08-28 | 燕山大学 | 一种基于小波与幻方变换的qr码加密、解密方法及系统 |
CN111654592B (zh) * | 2020-06-05 | 2022-02-25 | 中国矿业大学 | 基于秘密分享的多图像隐藏方法 |
CN111654592A (zh) * | 2020-06-05 | 2020-09-11 | 中国矿业大学 | 基于秘密分享的多图像隐藏方法 |
CN111681153A (zh) * | 2020-06-08 | 2020-09-18 | 重庆邮电大学 | 基于非下采样轮廓波变换的cs-svd的鲁棒数字水印算法 |
CN112529756A (zh) * | 2020-12-01 | 2021-03-19 | 湖南长城信息金融设备有限责任公司 | 一种数字图像水印嵌入方法、提取方法、装置、存储介质 |
CN112529756B (zh) * | 2020-12-01 | 2024-05-31 | 长城信息股份有限公司 | 一种数字图像水印嵌入方法、提取方法、装置、存储介质 |
CN113822786B (zh) * | 2021-11-25 | 2022-02-22 | 成都盛思睿信息技术有限公司 | 一种强鲁棒性的隐形数字图像水印处理方法及装置 |
CN113822786A (zh) * | 2021-11-25 | 2021-12-21 | 成都盛思睿信息技术有限公司 | 一种强鲁棒性的隐形数字图像水印处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108564519B (zh) | 2021-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108564519A (zh) | 一种基于多变换域的彩色图像数字水印方法 | |
CN105512999B (zh) | 一种双变换的彩色图像全息水印方法 | |
Kumar et al. | DWT based color image watermarking using maximum entropy | |
Divecha et al. | Implementation and performance analysis of DCT-DWT-SVD based watermarking algorithms for color images | |
Patil et al. | DWT based invisible watermarking technique for digital images | |
CN106803229A (zh) | 基于相位奇异值校正的图像水印方法 | |
Ernawan et al. | Image watermarking based on integer wavelet transform-singular value decomposition with variance pixels | |
Muñoz-Ramirez et al. | A robust watermarking scheme to JPEG compression for embedding a color watermark into digital images | |
CN109754358B (zh) | 一种基于显著性检测和轮廓波变换的图像水印方法与系统 | |
Loukhaoukha et al. | Hybrid watermarking algorithm based on SVD and lifting wavelet transform for ownership verification | |
CN104766263A (zh) | 基于四元数Legendre矩校正的彩色图像水印嵌入及检测方法 | |
KR100337954B1 (ko) | 텍스트 워터마크 삽입 방법 및 장치 | |
Rahim et al. | Exploiting de-noising convolutional neural networks DnCNNs for an efficient watermarking scheme: a case for information retrieval | |
Zhiwei et al. | Steganography based on wavelet transform and modulus function | |
CN108765249B (zh) | 基于方向各向异性和奇异值分解相结合的彩色水印方法 | |
Kekre et al. | Robust Watermaking Technique Using Hybrid Wavelet Transform Generated From Kekre Transforms and Discrete Cosine Transform | |
CN110428355A (zh) | 一种基于压缩感知与nsct的彩色图像数字水印方法 | |
CN114037593B (zh) | 基于逆向直方图平移的可逆图像水印算法 | |
Chaitanya et al. | Digital color image watermarking using DWTDCT coefficients in RGB planes | |
Kaur et al. | Comparative analyses of YCbCr color space and CIELab color space based on DWT and SVD | |
Chang et al. | A multipurpose wavelet-based image watermarking | |
CN113034333A (zh) | Zernike矩和DWT-NSCT的数字水印算法 | |
Hu et al. | A blind watermarking algorithm for color image based on wavelet transform and Fourier transform | |
Tomar et al. | A statistical comparison of digital image watermarking techniques | |
CN113139890B (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 |