专利名称：一种数字水印的嵌入方法
技术领域：
本发明涉及一种基于小波变换的数字水印技术，是一种多媒体数据保护方法，属于多媒体信号处理领域。
数字媒体有一些模拟媒体不可比拟的优点，如数字信号高质量、编辑加工容易、拷贝不失真、易于通过电子的(网络)或物理的(CD-ROM)系统低价高效地迅速传输和分配等。这些优点使得最近几年数字媒体(数字音频、数字图象、数字视频)的技术开发和应用迅猛发展。然而，也正是由于这些优点，给媒体原始拥有者的权益造成了潜在的威胁其艰苦劳动的成果有可能在一夜之间被无偿地大批完美复制并传遍世界的每个角落。因而数字媒体的版权保护成为一个迫切需要解决的问题。数字水印则为解决数字媒体版权保护问题的一个有效办法。近几年来，数字水印技术在许多应用领域体现了它的重要性，并得到了广泛的重视。
目前大多数数据隐藏的研究和文献集中于图象水印。根据水印嵌入的方式而言，目前提出的数字水印技术主要可分为二类空域和变换域技术。前者通过直接改变主图象某些象素值来嵌入水印。后者是在变换域中实现，先将图象做某种变换，然后把水印嵌入到图象的变换域中。由于DWT(Discrete Wavelet Transform)良好的空间一频率局部特性和与人眼视觉特性相符的变换机制，在新一代静止图象压缩标准JPEG 2000占据了核心位置，并且逐渐代替DCT成为变换域数字水印算法的主要工具。
著名并被人们广泛接受的是Cox等人的观点，DCT域水印应嵌入幅值最大前k个AC系数中。黄等进一步推广了这一思想，提出了利用DC分量来嵌入水印。但在DWT域，这方面的研究尚未见报道，并且嵌入公式和小波分解的级数的研究也尚属空白。
本发明的目的是提出一种基于小波变换的数字水印嵌入方法，用于提高水印的稳健性，保护数字媒体版权。
为了实现上述目的，本发明方法采用如下三个步骤1)根据所要嵌入水印数据量的多少尽量提高小波分解的级数，确定分解级数后，对原始数据进行小波分解；2)嵌入公式是vi′＝vi+αxi，其中，vi是小波系数，α是拉伸因子，xi是水印分量；3)按嵌入公式水印首先嵌入小波系数的低频带，若有剩余，再按小波系数频带重要性的排序顺序嵌入高频带。
步骤1)根据水印数据量的多少确定小波分解的级数，小波分解的级数要尽量高，然后进行小波分解。嵌入公式vi′＝vi+αxi中的拉伸因子α在保证水印不可见的前提下，尽可能大。
小波系数频带重要性的排序顺序选择可分离的滤波器组，对输入图象进行小波分解，产生LHj、HLj、HHj三个高频带系列，一个LL3低频带(三级分解时)(如
图1所示)。其中低频带表示由小波变换分解级数决定的最大尺度、最小分辨率下对原始图象的最佳逼近。它的统计特征和原图象相似，图象的大部分能量集中在此。高频带系列则分别是图象在不同尺度、不同分辨率下的细节信息。分辨率越低，其中有用信息的比例越高。也就是说，经过小波分解把一个图象分成了若干级。对于同一级图象，低频子图象LLj最重要，其次是HLj与LHj，高频子图象HHj相对最不重要。对于不同级来说，级高者重要，级低者不重要。所以，小波图象子频带按其重要性总体的排序为LLk、HLk、LHk、HHk、HLk-1、LHk-1、HHk-1、…、HL1、LH1、HH1(如图1所示)。
本发明与现有水印技术比较有以下优点通过研究DWT域的嵌入方法(包括嵌入公式、嵌入位置和小波分解的级数)与稳健性的关系，找到了比较好的嵌入方法。与现有水印技术相比，该技术显著提高了水印的稳健性。
以下我们从理论和实验数据说明本发明带来的积极效果。
1)小波分解的级数Mallat在1988年提出多分辨率分析(MRA)，它是信号分析和处理的有效工具。在MRA的基础上，S.Mallat提出了用子带结构实现DWT的算法，统一了子带滤波器与小波变换的计算。
给出平方可和的二维信号{xm，n}m，n∈Z，令c0，m，n＝xm，n，m，n∈Z则二维信号的小波分解公式为cj,m,n=2Σk,lhl-2mhl-2ncj-1,k,l,dj,m,n1=2Σk,lhk-2mgl-2ncj-1,k,l,]]>dj,m,n2=2Σk,lgk-2mhl-2ncj-1,k,l,dj,m,n3=2Σk,lgk-2mgl-2ncj-1,k,l.]]>其中 gk＝(-1)nh1-n， 从上面公式可知，小波分解是每级以2加权，并注意到低通滤波器系数和为1，从而随着小波分解级数的增加，低频系数的幅值以近似2的倍数增长。而水印编码可看为在强背景(原始图象)下迭加一个弱信号(水印)，只要迭加的信号低于对比度门限，视觉系统就无法感觉到信号的存在。而根据Weber定律，对比度门限和背景信号的幅值成比例。这就意味着随着小波分解级数的增加，嵌入水印的强度大幅度增加，从而水印的稳健性增强。同时，小波分解的级数越多，水印分量可以更好地扩散。因此在水印算法中，应根据水印数据量的多少，尽可能提高小波分解的级数。
2)嵌入位置从多分辨分析的角度考虑小波图象的各个频带时，这些频带之间并不是纯粹无关的。对于各个高频带，由于它们是图象同一个边缘、轮廓和纹理信息在不同方向、不同尺度和不同分辨率下由粗到细的描述，它们之间存在着一定的关系小波系数的幅值随着尺度的减小而衰减，即使图象函数有奇异点，只要在奇异点的邻域里振荡不激烈。若图象函数f(x，y)在(x0，y0)邻域内具有Lipschitz α连续，则尺度为2j、位于(p，q)的小波系数的幅值有|αj(p，q)|≤A2j(α+1)其中A为一正常数。因此小波系数的幅值随着尺度的减小而衰减。振荡模型不具有这个性质，幸而振荡模型在实际图象中不常遇到。基于零树的编码正是利用了小波系数的幅值随着尺度的减小而衰减这个性质而获得了极大成功。
随机序列、有意义的文本、图象等都可以作为水印。不同类型的水印包含的数据量不同。不同数据量的水印应有不同的嵌入对策。
水印应放在哪里才足够稳健？在DWT域，不同的小波系数作为水印载体，水印有不同的稳健性。
水印长度较短当水印长度不大于小波图象低频系数个数时，水印应嵌入哪里？首先，Cox等提出水印应放在HVS感觉上最重要的分量上。这个观点现在已被人们广泛接受。感觉上重要的分量是图象信号的主要成分，携带较多的信号能量，在图象有一定失真的情况下，仍能保留主要成分。但Cox等把DC系数排除在外，原因在于避免加水印的图象出现方块效应。由于小波变换是全局变换，把水印加在低频系数的水印图象不会出现方块效应。同时由前面可知低频带是对原始图象的低通逼近，图象的大部分能量集中在此，而高频带系列只是图象的细节信息。因此水印应当首先嵌入小波图象低频系数。
其次，同第一步相同的道理，由于低频系数的幅值一般远大于高频系数，从而具有较大的感觉容量，嵌入一定强度的水印后不会引起原始图象视觉质量的明显改变。
再次，根据信号处理理论，嵌入水印的图象最有可能遭遇到的信号处理过程，如数据压缩、低通滤波、次抽样、插值、D/A和A/D转换等，对低频系数的保护比高频系数好。即这些小波系数在经过常见的信号处理和噪声干扰后仍能很好地保留，不过多地为信号处理和噪声干扰所改变。
因此，水印分量嵌入到DWT域低频系数具有足够的稳健性。
水印长度较长当水印长度大于小波图象低频系数个数时，除了把水印嵌入低频带系数外，首先利用HL3(三级分解)带的系数来嵌入水印。这是因为1)HL3带在高频带系列中最重要。2)HL3带的系数通常比其它尺度的系数大(见本部分开头)。同理，接下来应利用LH3带的系数。以此类推，我们得到如下结论当水印长度较长时，按小波图象频带重要性的排序顺序嵌入水印稳健性最好(图1)。
综上所述，我们通过对小波图象系数的分布特点和振幅进行定性、定量分析，得到了一个新的嵌入对策水印应当首先嵌入小波图象低频系数，若有剩余，再按小波图象频带重要性的排序顺序嵌入高频带。
3)嵌入公式嵌入公式也影响水印的稳健性。目前常用的嵌入公式有两个(1)vi′＝vi+axi；(2)vi′＝vi(1+αxi)，其中α为拉伸因子。
由于嵌入公式(1)对每个小波系数叠加相同的强度，水印图象受攻击后，水印分量可能都以一定的强度存在。但对嵌入公式(2)来说，由于大的系数叠加大的强度，小的系数叠加小的强度，水印图象受攻击后，小系数的水印分量可能不存在。只有在攻击强度很强的情况下，基于嵌入公式(1)的水印分量不存在了，基于嵌入公式(2)的大系数的水印分量仍存在。但这种情况很少存在。
因此，嵌入公式(1)是比较好的选择。
我们用纹理比较简单和比较复杂的图象“Lena”(256×256×8bits)和“Baboon”(256×256×8bits)为测试图象。实验数据表明新方法是非常稳健的。
图1是图象小波分解图。
图2是水印的不可性比较图。
图3为条件完全相同的二种情况水印图象在JPEG压缩和Gaussian噪声干扰下的稳健性能比较图。
图1中，对输入图象进行小波分解，产生LHj、HLj、HHj三个高频带系列，一个LL3低频带(三级分解时)。
图2中，(a)为利用LL3嵌入水印的水印图象；(b)为利用HL3嵌入水印的水印图象。在完全相同的条件下，当PSNR都为44.4dB时，“Lena”第一种情况水印不可见，而第二种情况水印明显可见，说明低频系数具有较大的感觉容量。
图3中，(a)为抗JPEG压缩性能比较；(b)为抗Gaussian噪声性能比较。纵轴表示从失真的水印图象中抽取的水印W*与原始水印W的相似度。图3表明在低频带系数嵌入水印比在HL3带系数嵌入水印稳健好。把水印嵌入其它高频子带上的实验也可得到类似结果。
表1为基于嵌入公式(1)的稳健性。表中，(a)Lena水印图象抗JPEG(一级分解)(b)Baboon水印图象抗noise(一级分解)(c))Baboon水印图象抗JPEG(三级分解)(d)Lena水印图象抗noise(三级分解)。
表2为基于嵌入公式(2)的稳健性。表2中，(a)Lena水印图象抗JPEG(一级分解)(b)Baboon水印图象抗noise(一级分解)(c))Baboon水印图象抗JPEG(三级分解)(d)Lena水印图象抗noise(三级分解)。
表1-2中小波基选Daubechies小波DbN(1≤N≤10)，所有Baboon抗JPEG都表示Baboon水印图象受JPEG压缩攻击后，PSNR分别为20db、20.9db、21 8db、22.4db、22.9db时，抽取的水印W*与原始水印W的相似度；表中所有Lena抗JPEG都表示Lena水印图象受JPEG压缩攻击后，PSNR分别为23.1db、25db、26.9db、28.1db、29db时，抽取的水印W*与原始水印W的相似度。表中所有Baboon抗noise都表示Baboon水印图象受加性Gaussian噪声攻击后，PSNR分别为18.7db、16.3db、14.6db、13.2db、12.2db时，抽取的水印W*与原始水印W的相似度；表中所有Lena抗noise都表示Lena水印图象受加性Gaussian噪声攻击后，PSNR分别为18.7db、16.3db、14.6db、13.1db、12.1db时，抽取的水印W*与原始水印W的相似度。
表1基于嵌入公式(1)的稳健性(a)(b)
(c) (d)
表2基于嵌入公式(2)的稳健性(a) (b)
(c) (d)
权利要求
1.一种数字水印的嵌入方法，其特征在于该方法的步骤为1)根据所要嵌入水印数据量的多少确定小波分解的级数；2)嵌入公式是vi′＝vi+αxi，其中，vi是小波系数，α是拉伸因子，xi是水印分量；3)水印按嵌入公式首先嵌入小波系数的低频带，若有剩余，再按小波系数频带重要性的排序顺序嵌入高频带。
2.根据权利要求1所述的一种数字水印的嵌入方法，其特征是步骤1)根据水印数据量的多少确定小波分解的级数，小波分解的级数要尽量高，然后进行小波分解。
3.根据权利要求1所述的一种数字水印的嵌入方法，其特征是嵌入公式Vi′＝vi+αxi中的拉伸因子α在保证水印不可见的前提下，尽可能大。
全文摘要
本发明涉及一种基于小波变换的数字水印技术,是一种多媒体视频数据保护方法,属于多媒体信号处理领域。本发明首先选定小波分解的级数,利用小波变换将多媒体数据进行分解,然后选择一些小波系数按一定的规则嵌入水印。最后通过小波逆变换得到嵌入水印的多媒体数据。与原始媒体相比,人眼或耳无法区别两者的差别。而数字水印则为这个问题提供了一种区别的手段。水印可以是任何有意义的数字文档或无意义的随机序列。本发明可使通过网络传播的多媒体数据或文件获得保护。
文档编号G06T1/00GK1377184SQ0211517
公开日2002年10月30日 申请日期2002年4月30日 优先权日2002年4月30日
发明者刘九芬, 黄达人, 黄继武 申请人:中山大学