专利名称：结合多尺度特征与中国余数定理的大容量数字水印方法
技术领域：
本发明涉及一种用于隐秘通信和数字图像的大容量数字水印方法，特别是一种结 合多尺度特征与中国余数定理的大容量数字水印方法。
背景技术：
互联网的开放性给人们的生活带来了很大的便利。借助于网络，人们可以很容易 地实现数据共享。但互联网的开放性也带来一个严峻的问题，即未经授权的多媒体产品被 任意地上传到网络并被无限制的下载，严重侵犯了作者的合法权益。如何对数字多媒体产 品进行有效的保护，已经成为学术界和产业界密切关注的问题。 数字图像水印技术通过向图像数据中加入具有标识意义的信息，在发生版权纠 纷的时候将其检测出来作为证据。在实际使用时，数字图像经常会受到无意或有意的变 换。一方面，为了达到理想的视觉效果或适应特殊应用的需要，图像经常会受到压缩、滤 波等信号处理操作；另一方面，恶意破坏者为了去除图像中的水印信息可能对图像进行有 意的变换。这两种操作都会对水印的正确检测造成破坏。因此，水印对各种常见变换的 鲁棒性是其实际应用的前提要求，尤其是对旋转、缩放、平移等几何变换的鲁棒性。由于 图像在受到几何变形时，图像像素重新分布，水印检测器如果没有针对性的设计，很难自 动地定位含水印区域，从而导致检测失败。目前能够抵抗几何变换的数字水印方法主要 分为以下几种(1)基于几何校正的方法。主要思想是对受到几何变换的图像在水印检 测之前进行校正，然后在校正后的图像中检测水印。几何校正的过程可以借助于图像配 准技术求出几何变换的参数然后进行逆变换；也可以在嵌入水印信号的同时嵌入模板，并 在检测时根据模板确定几何变换参数进行相应的逆变换。基于几何校正的方法在水印检 测时经常需要原始图像的参与，因此很难实现盲检测(水印检测时不需要原始图像)，目 前对其研究较少。(2)基于几何不变域的水印嵌入。主要思想是借助于对几何操作具有 不变性的图像变换将水印信号嵌入其中。常用的几何不变域有Fourier-MelIin变换域、 图像归一化域等。基于 Fourier-Mellin 变换的算法C. Lin et al. "Rotation, Scale, and Translation Resilient Watermarking of Images，，，IEEE Transactions on Image Processing, vol. 10，no. 5，pp. 767-782，2001.基于图像归一化的算法:M. Alghoniemy et al. "Geometric Invariance in Image Watermarking", IEEE Transactions on Image Processing, vol. 13，no. 2，pp. 145-153，2004.均是此类方法中的经典方法。此类算法的 缺点在于水印嵌入过程中往往伴随有插值运算，造成加水印图像的质量下降；此外，插值 运算还会对水印检测造成不利影响。(3)基于图像特征的水印嵌入。从图像中提取的局 部特征点一般具有几何不变性，以这些特征点为参考确定一定的区域用于嵌入水印，则水 印信号具有对几何变换的不变性，这正是此类算法的基本出发点。基于图像特征的水印嵌 入方法能同时抵抗常规的图像处理操作和几何变换，且加水印图像质量好，因此是近年来 该令页域石if究的热点。文献 H. Y. Lee et al. “Robust image watermarking using local invariant features, ” Optical Engineering, vol.45, no.3, pp.037002(1-11),2006.禾口文献 Χ· Y. Wang et al. "A new digital image watermarking algorithm resilient to desynchronization attacks,，，IEEE Transactions on Information Forensics and Security, vol. 2, no. 4, pp. 655-663，2007.是近年来此类算法中较为经典的算法。Lee等 提出的算法首先利用尺度不变特征变换技术检测特征点，然后依据特征点在图像中确定局 部区域，最后将水印信号在空间域叠加到局部图像块上。检测时，通过Wiener滤波的方法 提取水印信号。该算法嵌入的水印信号为具有正态分布的随机二值图像，且水印嵌入在空 间域进行。由于水印嵌入过程在空域进行，类似于叠加随机噪声，水印鲁棒性不理想。Wang 等提出的方法利用Harris-Laplace特征提取算法检测特征点，确定局部区域，然后将二值 水印序列在离散傅里叶变换域(DFT)嵌入。水印提取时，通过比较局部区域中对应位置的 傅里叶系数进行。该算法由于在变换域嵌入水印，水印鲁棒性较好，但水印的容量有限，仅 有32比特，而且当图像受到旋转的时候需要预先进行校正。近年来出现的各种算法基本思 想和以上两种方法类似。目前，基于图像特征的水印算法中大多都将水印在空间域嵌入，导 致水印鲁棒性不理想；在变换域嵌入的水印算法容量均为几十比特，水印容量低。如何将大 容量的水印信号(二值图像)在变换域嵌入到由特征点确定的局部区域中，是目前研究的 难点。
发明内容
本发明的目的是要提供一种在嵌入大容量水印信号的同时，使水印能够同时抵 抗常规的信号处理操作和几何变换，对数字图像进行有效保护的结合多尺度特征与中国余 数定理的大容量数字水印方法。本发明的目的是这样实现的该大容量数字水印方法包括水印嵌入和水印提取二 个过程；一、水印嵌入过程包括以下步骤(1)选取圆形的二值黑白图像作为水印图像；当需要嵌入的水印图像为其他形状 时，选其外接圆图像作为水印图像，记水印信号为W。(2)对原始图像，利用尺度不变特征变换算法提取多尺度特征点，并以特征点为中 心，特征尺度的整数倍作为半径确定圆形区域。(3)以尺度不变特征变换特征点检测过程中对应的高斯差函数值作为特征点稳定 性的定量度量，结合步骤(2)中确定的区域半径，从图像中筛选出局部最稳定的特征点，并 确定出其所对应的圆形区域。(4)以高斯差函数值作为特征点稳定性强弱的定量度量，对步骤(3)中相互之间 存在重叠的区域进行优化，确定互不重叠的圆形区域。(5)对提取的每一个区域，进行一层非下采样轮廓波变换，并取低频子带内接圆内 的系数作为水印嵌入所要修改的系数。(6)将水印图像W进行缩放，缩放后的二值水印图像与非下采样轮廓波变换的低 频子带大小相同。(7)水印嵌入时利用中国余数定理，根据水印比特对与其具有相同坐标处的低频 子带系数进行修改，对修改后的系数进行非下采样轮廓波重构，得到含水印的图像块。 (8)用含水印的图像块替换原始图像中的对应区域，得到整幅含水印的图像。
二、所述的水印提取过程包括如下步骤(1)对待检测图像，利用尺度不变特征变换算法提取多尺度特征点，并以特征点为中心，特征尺度的整数倍作为半径确定圆形区域，特征尺度的放大倍数与水印嵌入时相同。(2)以尺度不变特征变换特征点检测过程中对应的高斯差函数值作为特征点稳定 性的定量度量，结合步骤(1)中确定的区域半径，从图像中筛选出局部最稳定的特征点，并 确定出其所对应的圆形区域。(3)以步骤⑵中特征点的稳定性程度作为筛选依据，对步骤(2)中相互间存在重 叠的区域进行优化，得到互不重叠的圆形区域。(4)对每一个局部区域，进行一层非下采样轮廓波变换，并取低频子带内接圆内的 系数作为水印提取所采用的系数。(5)利用中国余数定理进行水印提取，并将提取的水印图像变换到与原始水印图 像相同的大小。(6)计算提取的水印图像与原始水印图像的篡改评价函数，作为水印的相似度度 量。也可以根据提取的水印图像对水印存在与否做出判断。所述的水印嵌入过程的步骤(3)中所述筛选出局部最稳定的特征点为初步确定 水印嵌入的圆形区域，按如下步骤进行设从原始图像中提取的特征点中特征尺度处于 [4,8]之间的特征点集合为Q1 = {ρι，ρ2，ρ3，…，pm}，m为特征点的个数，对Q1中的特征点 P1，以其坐标作为圆心，特征尺度的整数倍为半径生成一个圆形区域R1，如果R1中包含有其 它特征点，则将Pji应的高斯差函数值与R1中其它特征点对应的高斯差函数值进行比较， 若为局部极大值，则保留P1，并从Q1中去掉存在于R1中的其它特征点，对Ωι中的所有特 征点进行上述操作，最后便得到局部最稳定的特征点，记为Ω2 = {qi, q2，q3，…，qn}，η为 局部最稳定特征点个数。所述的水印嵌入过程的步骤(4)中所述的以高斯差函数值作为特征点稳定性强 弱的定量度量，为对相互之间存在重叠的区域进行优化，筛选出互不重叠的区域，按如下步 骤进行(1)对局部最稳定特征点集合Ω2中的特征点qi，q2，q3，…，qn，首先筛选出稳定性 最好的点，设为由Qi与其特征尺度共同决定的圆形区域记为Ri,并将Qi从Ω2中去除， 记入集合Ω 3 ；(2)判断Ω 2中剩余特征点对应的圆形区域与Ri是否重叠，若重叠则从Ω2中去掉 相应的特征点，否则保留；(3)对集合Ω2中的特征点，重复步骤⑴和步骤(2)，直到处理完所有的特征点， 得到新的特征点集合ω3;(4)对Ω3中的每一个特征点，结合其特征尺度，最终确定互不重叠的圆形区域，用 于水印嵌入。所述的水印嵌入过程的步骤(5)中所述的对每一个区域进行一层非下采样轮廓 波变换，并选取低频子带内接圆内的系数是用于水印嵌入，水印嵌入前，首先将水印图像变 换到与低频子带相同的大小。所述的水印嵌入过程的步骤(7)中所述的水印嵌入时利用中国余数定理，且在水 印嵌入时，选取互质的两个整数分别为M1 = 13，M2 = 20，具体嵌入时根据所要嵌入的水印比特对非下采样轮廓波变换低频子带内接圆内的系数进行修改，强制使得系数对礼、M2的 余数Rei、Re2满足特定的大小关系，即待嵌入水印位为1时，Re1 ^ Re2 ；待嵌入水印位为O 时，Re1 < Re20有益效果该方法在自然图像的局部 区域中嵌入大容量的二值水印图像，提取的 水印信号既能抵抗常规的压缩、滤波等图像处理操作，又能抵抗旋转、缩放、平移等几何变 换，且提取的水印图像与含水印图像所受的几何变形相一致。能用于隐秘通信和数字图像 的版权保护，具体优点如下(1)水印方法具有对图像旋转、缩放、平移等几何变换的不变性。换句话说，当含 水印图像受到上述几何变换后，本发明中的方法仍能够有效地进行水印提取。这主要是因 为本发明中水印信号是嵌入在由多尺度特征点所确定的局部区域中。与传统的特征点不同 之处在于多尺度特征点是从图像的尺度空间中检测出的，每一个特征点除了具有位置信息 夕卜，还对应一个特征尺度。该特征尺度跟随图像的尺度变换而变化，因此本发明中以特征点 位置作为中心，特征尺度的整数倍作为半径生成圆形区域作为水印嵌入和提取的区域，这 些区域本身具有对图像几何变换的不变性。因此，将水印信号嵌入到这些区域中，就使得水 印信号具有几何不变性。(2)加水印后的图像质量好。相比于常规的水印算法，本发明将水印嵌入到由多尺 度特征点确定的局部区域中，而非整幅图像中。水印嵌入过程中在利用中国余数定理时，通 过设置合理的参数，可以在保持水印获得最优鲁棒性的基础上，使得图像质量最优。对一般 的自然图像，加水印后图像的峰值信噪比一般在40dB(分贝)以上。(3)水印嵌入过程，在用含水印图像块替换原始图像块时不会造成水印信息的丢 失。由于实际实现时水印嵌入和提取使用的区域是以圆形区域为内接圆，边界部分补零的 方形图像块，如果采用常规的离散余弦变换或离散傅里叶变换嵌入水印，水印信息会散布 到边缘区域，在图像块替换的时候丢失水印信息。本发明中利用具有时频局部化特性的非 下采样轮廓波变换，通过修改低频子带内接圆内的系数进行水印嵌入，有效避免了水印信 息丢失的问题。(4)水印安全性好。本发明利用中国余数定理嵌入水印，嵌入时设置的互质的两个 整数起到密钥的作用。水印提取时，如果没有该密钥则无法进行水印提取，提高了方法的安 全性。
图1为本发明水印嵌入框图。图2为本发明水印提取框图。图3-A为本发明中水印嵌入/提取区域的生成过程初始特征点图。图3-B为本发明中水印嵌入/提取区域的生成过程局部最稳定特征点图。图3-C为本发明中水印嵌入/提取区域的生成过程最终采用的特征点图。图4为本发明水印嵌入过程中图像块的处理。图5为本发明具体水印嵌入过程。图6-A1为本发明水印嵌入效果的原始图像图。图6-A2为本发明水印嵌入效果的原始图像图。
图6-B1为本发明水印嵌入效果的加水印后的图像图。图6-B2为本发明水印嵌入效果的加水印后的图像图。图6-C1为本发明水印嵌入效果的原始图像与加水印后图像的差值图像图。图6-C2为本发明水印嵌入效果的原始图像与加水印后图像的差值图像图。图7-1为本发明水印提取80% JPEG压缩效果图。图7-12为本发明水印提取80% JPEG压缩图。图7-2为本发明水印提取锐化效果图。图7-22为本发明水印提取锐化图。图7-3为本发明水印提取缩放0. 9倍效果图。图7-32为本发明水印提取缩放0. 9倍图。图7-4为本发明水印提取缩放1. 2倍效果图。图4-42为本发明水印提取缩放1. 2倍图。图7-5为本发明水印提取逆时针旋转30度效果图。图7-52为本发明水印提取逆时针旋转30度图。图7-6为本发明水印提取顺时针旋转45度效果图。图7-62为本发明水印提取顺时针旋转45度图。
具体实施例方式实施例1 本发明中的非下采样轮廓波其英文表达为Contourlet ；尺度不变特征 变换其英文缩写表达为SIFT ；高斯差函数的英文表达为DoG ；篡改评价函数其英文缩写表 达为TAF ；描述符的英文表达为descriptor ；离散余弦变换的英文缩写表达为DCT ；离散傅 里叶变换的英文缩写表达为DFT。该大容量数字水印方法包括水印嵌入和水印提取二个过程；一、水印嵌入过程包括以下步骤(1)选取圆形的二值黑白图像作为水印图像；当需要嵌入的水印图像为其他形状 时，选其外接圆图像作为水印图像，记水印信号为W。(2)对原始图像，利用尺度不变特征变换算法提取多尺度特征点，并以特征点为中 心，特征尺度的整数倍作为半径确定圆形区域。(3)以尺度不变特征变换特征点检测过程中对应的高斯差函数值作为特征点稳定 性的定量度量，结合步骤(2)中确定的区域半径，从图像中筛选出局部最稳定的特征点，并 确定出其所对应的圆形区域。(4)以高斯差函数值作为特征点稳定性强弱的定量度量，对步骤(3)中相互之间 存在重叠的区域进行优化，确定互不重叠的圆形区域。(5)对提取的每一个区域，进行一层非下采样轮廓波变换，并取低频子带内接圆内 的系数作为水印嵌入所要修改的系数。(6)将水印图像W进行缩放，缩放后的二值水印图像与非下采样轮廓波变换的低 频子带大小相同。(7)水印嵌入时利用中国余数定理，根据水印比特对与其具有相同坐标处的低频 子带系数进行修改，对修改后的系数进行非下采样轮廓波重构，得到含水印的图像块。
(8)用含水印的图像块替换原始图像中的对应区域，得到整幅含水印的图像。二、所述 的水印提取过程包括如下步骤(1)对待检测图像，利用尺度不变特征变换算法提取多尺度特征点，并以特征点为 中心，特征尺度的整数倍作为半径确定圆形区域，特征尺度的放大倍数与水印嵌入时相同。(2)以尺度不变特征变换特征点检测过程中对应的高斯差函数值作为特征点稳定 性的定量度量，结合步骤(1)中确定的区域半径，从图像中筛选出局部最稳定的特征点，并 确定出其所对应的圆形区域。(3)以步骤⑵中特征点的稳定性程度作为筛选依据，对步骤(2)中相互间存在重 叠的区域进行优化，得到互不重叠的圆形区域。(4)对每一个局部区域，进行一层非下采样轮廓波变换，并取低频子带内接圆内的 系数作为水印提取所采用的系数。(5)利用中国余数定理进行水印提取，并将提取的水印图像变换到与原始水印图 像相同的大小。(6)计算提取的水印图像与原始水印图像的篡改评价函数，作为水印的相似度度 量。也可以根据提取的水印图像对水印存在与否做出判断。所述的水印嵌入过程的步骤(3)中所述筛选出局部最稳定的特征点为初步确定 水印嵌入的圆形区域，按如下步骤进行设从原始图像中提取的特征点中特征尺度处于 [4,8]之间的特征点集合为Q1 = {ρι，ρ2，ρ3，…，pm}，m为特征点的个数，对Q1中的特征点 P1，以其坐标作为圆心，特征尺度的整数倍为半径生成一个圆形区域R1，如果R1中包含有其 它特征点，则将Pji应的高斯差函数值与R1中其它特征点对应的高斯差函数值进行比较， 若为局部极大值，则保留P1，并从Q1中去掉存在于R1中的其它特征点，对Ωι中的所有特 征点进行上述操作，最后便得到局部最稳定的特征点，记为Ω2 = [qi, q2，q3，…，qn}，η为 局部最稳定特征点个数。所述的水印嵌入过程的步骤(4)中所述的以高斯差函数值作为特征点稳定性强 弱的定量度量，为对相互之间存在重叠的区域进行优化，筛选出互不重叠的区域，按如下 步骤进行(1)对局部最稳定特征点集合Ω2中的特征点qi，q2，q3，…，qn，首先筛选出稳定性 最好的点，设为由Qi与其特征尺度共同决定的圆形区域记为Ri,并将Qi从Ω2中去除， 记入集合Ω 3 ；(2)判断Ω2中剩余特征点对应的圆形区域与Ri是否重叠，若重叠则从Ω2中去掉 相应的特征点，否则保留；(3)对集合Ω2中的特征点，重复步骤⑴和步骤(2)，直到处理完所有的特征点， 得到新的特征点集合ω3;(4)对Ω3中的每一个特征点，结合其特征尺度，最终确定互不重叠的圆形区域，用 于水印嵌入。所述的水印嵌入过程的步骤(5)中所述的对每一个区域进行一层非下采样轮廓 波变换，并选取低频子带内接圆内的系数是用于水印嵌入，水印嵌入前，首先将水印图像变 换到与低频子带相同的大小。所述的水印嵌入过程的步骤(7)中所述的水印嵌入时利用中国余数定理，且在水印嵌入时，选取互质的两个整数分别为M1 = 13，M2 = 20，具体嵌入时根据所要嵌入的水印 比特对非下采样轮廓波变换低频子带内接圆内的系数进行修改，强制使得系数对礼、M2的 余数Rei、Re2满足特定的大小关系，即待嵌入水印位为1时，Re1 ^ Re2 ；待嵌入水印位为0 时，Re1 < Re20在图1、图3-A、图3-B、图3-C、图4和图5中，本发明充分利用多尺度特征点对图 像旋转、缩放、平移等几何变换的不变性，从图像中提取出局部最稳定的特征点，结合其特 征尺度确定圆形的水印嵌入区域，将二值水印图像经过适当的调整后在非下采样轮廓波域 结合中国余数定理进行嵌入。水印嵌入步骤如下第一步确定水印嵌入的区域。对一幅原始图像I(x，y)，首先利用尺度不变特征变换算法检测多尺度特征点。对 I (X，y)，SIFT特征点的检测步骤具体如下(1)构造原始图像的尺度空间L(x，y，σ )L (x，y，σ ) =I (χ, y) *G (χ, y, σ )其中，*表示卷积，G(x，y，σ)表示可变方差的高斯函数，定义如下G(x,y,a) =
2πσ(2)利用方差相差k倍的高斯差函数与原始图像卷积得到高斯差(Difference of Gaussian, DoG)图像 D(x，y，σ)，方法如下D (x，y，σ ) =I (χ, y) Φ [G (χ, y，k σ ) -G (χ，y，σ )]检测SIFT特征点时，将每一级高斯差图像中的每个像素与同尺度下的8个邻域像 素和上下尺度下对应位置的各9个像素点，共计26个点进行比较。当所比较的点为局部 极小值或极大值时，确定该点为特征点，否则舍弃。(3)在确定初始特征点之后，利用Hessian矩阵去除处于低对比度区域和 图像边缘的特征点。然后根据每个特征点的邻域信息确定特征点的主方向和描述符 (descriptor)0从图像中提取的尺度不变特征变换特征点信息包括位置、特征尺度、主方向、描述 符。本发明中利用其中的坐标信息和特征尺度确定水印嵌入和提取的位置。具体方法如 下设所采用的特征点坐标为(s，t)，特征尺度为σ，本发明中所确定的圆形图像块 为(x-s)2+ (y-t)2 = (k σ )2其中，k表示特征尺度放大的倍数，决定区域的半径。在实际实现中，特征尺度过 大或过小的特征点稳定性较差，尤其是当图像发生缩放变换的时候，很容易丢失。因此，本 发明中对尺度不变特征变换特征点进行进一步筛选，采用特征尺度处于区间[4，8]之间的 特征点。图3-A为利用尺度不变特征变换特征点在原始图像中确定的圆形区域。图3-A中的图像块之间存在重叠，如果直接嵌入水印则相互之间影响，致使水印 检测困难。为了获得互不重叠的区域，需对图3-A中的重叠区域进行优化，选择其中稳定性 最好的特征点对应的区域，去除稳定性差的区域。区域的稳定性由圆心处尺度不变特征变 换特征点对应的高斯差函数值决定，其值越大稳定性越好，否则稳定性差。本发明将首先从尺度不变特征变换特征点中筛选出局部最稳定特征点，具体方法如下设图3-A中所对应的特征点集合为Ω1= {ρι，ρ2，ρ3，…，pm}，m为特征点的个数。对Q1中的特征点P1，利用前面介绍的方法生成一个圆形区域R1。如果R1中包含有其它特征 点，则将P1对应的高斯差函数值与R1中其它特征点对应的高斯差函数值进行比较，若为局 部极大值，则保留P1，并从Q1中去掉存在于R1中的其它特征点。对Q1中的所有特征点进 行上述操作，最后便得到局部最稳定尺度不变特征变换特征点，记为Ω2= {qi，q2，q3，…， qn}，η为局部最稳定特征点个数。图3-Β所示为局部最稳定SIFT特征点对应的圆形区域。 在图3-Β中，每一个圆形区域内部仅存在一个尺度不变特征变换特征点，即局部最稳定的 点ο图3-Β中区域的重叠情况虽有改善，但仍未完全消除。为此，本发明进一步进行处 理。具体方法如下(1)对Ω2中的特征点Q1, q2, q3,…，qn，首先筛选出稳定性最好的点，设为q"由 Qi与其特征尺度共同决定的圆形区域记为Ri,并将Qi从Ω 2中去除，记入集合Ω3。(2)判断Ω 2中剩余特征点对应的圆形区域与Ri是否重叠，若重叠则从Ω2中去掉 相应的特征点；否则保留。(3)对集合Ω2中的特征点，重复步骤(1)和步骤(2)，直到处理完所有的特征点， 得到新的特征点集合ω3。(4)对Ω3中的每一个特征点，结合其特征尺度，最终确定互不重叠的圆形区域，用 于水印嵌入与提取。图3-C为本发明中所采用的用于水印嵌入和提取的区域。第二步对圆形区域进行预处理本发明中所采用的水印嵌入区域为圆形的图像区域，但由于数字图像的离散性， 实际实现中所采用的图像块必须为矩形。因此，必须首先对圆形图像块进行预处理。在图 4中，本发明中，在水印嵌入之前首先将圆形区域的边界部分补零，然后在变换域嵌入水印 信号。对含水印图像块，首先将边界部分去掉再替代原始图像中的对应部分得到加水印图像。第三步对图像块进行非下采样轮廓波变换为了获得较好的鲁棒性，本发明将水印信号嵌入到变换域。即首先对图像块进行 分解得到一系列系数，然后根据水印信号对系数进行修改，最后通过重构得到含水印图像。 本发明中，我们将大容量水印信号嵌入到非下采样轮廓波变换NSCT域。NSCT是多尺度几 何分析工具中的一种，是近年来出现的一种性能优良的时频分析工具。NSCT分解有更大的 灵活性，可以根据需要确定分解的方向数，并且分解得到的低频子带和高频子带大小与原 始图像相同。由于本发明中将水印嵌入到图像的局部区域中，因此相比于传统的小波变换， NSCT能够产生更多的系数，有利于实现大容量水印信号嵌入。第四步基于中国余数定理的大容量水印嵌入本发明中所采用的水印信号为圆形的二值图像，在图5中。在嵌入之前，首先对水 印图像进行一定的调整，将其变换到与非下采样轮廓波变换低频子带相同的大小。然后将 水印比特嵌入到与之具有相同坐标的NSCT系数中。为了获得更好的鲁棒性，选取低频子带 作为水印嵌入所要修改的系数。设所要嵌入的水印比特为w(x，y)，首先选取坐标(x，y)处的低频系数，记为C (x，y)。然后将w(x，y)结合中国余数定理进行嵌入。中国余数定理，也称孙子剩余定理，源自于“韩信点兵”的故事。韩信在点兵的时 候，为了保住军事机密，不让敌人知道自己部队的实力，先令士兵从1至3报数，然后记下最 后一个士兵所报之数；再令士兵从1至5报数，也记下最后一个士兵所报之数；最后令士兵 从1至7报数，又记下最后一个士兵所报之数；这样，他很快就算出了自己部队士兵的总人 数，而敌人则始终无法弄清他的部队究竟有多少名士兵。这个故事中所说的韩信点兵的计 算方法，就是现在被称为“中国余数定理”的一次同余式解法。数学描述如下设μ表示一个由r个整数组成的集合μ = (M1, M2,…，Mj，其中任意两个整数 Mi与Mj之间相对互质，即两者之间除了整数1之外没有公约数。另有一联立全等式为
其中Rei称为余数。整数Z可以通过下式计算
权利要求
1.一种结合多尺度特征与中国余数定理的大容量数字水印方法，其特征是该大容量 数字水印方法包括水印嵌入和水印提取二个过程；一、水印嵌入过程包括以下步骤(1)选取圆形的二值黑白图像作为水印图像；当需要嵌入的水印图像为其他形状时， 选其外接圆图像作为水印图像，记水印信号为W ；(2)对原始图像，利用尺度不变特征变换算法提取多尺度特征点，并以特征点为中心， 特征尺度的整数倍作为半径确定圆形区域；(3)以尺度不变特征变换特征点检测过程中对应的高斯差函数值作为特征点稳定性的 定量度量，结合步骤O)中确定的区域半径，从图像中筛选出局部最稳定的特征点，并确定 出其所对应的圆形区域；(4)以高斯差函数值作为特征点稳定性强弱的定量度量，对步骤(3)中相互之间存在 重叠的区域进行优化，确定互不重叠的圆形区域；(5)对提取的每一个区域，进行一层非下采样轮廓波变换，并取低频子带内接圆内的系 数作为水印嵌入所要修改的系数；(6)将水印图像W进行缩放，缩放后的二值水印图像与非下采样轮廓波变换的低频子 带大小相同；(7)水印嵌入时利用中国余数定理，根据水印比特对与其具有相同坐标处的低频子带 系数进行修改，对修改后的系数进行非下采样轮廓波重构，得到含水印的图像块；(8)用含水印的图像块替换原始图像中的对应区域，得到整幅含水印的图像；二、所述的水印提取过程包括如下步骤(1)对待检测图像，利用尺度不变特征变换算法提取多尺度特征点，并以特征点为中 心，特征尺度的整数倍作为半径确定圆形区域，特征尺度的放大倍数与水印嵌入时相同；(2)以尺度不变特征变换特征点检测过程中对应的高斯差函数值作为特征点稳定性的 定量度量，结合步骤(1)中确定的区域半径，从图像中筛选出局部最稳定的特征点，并确定 出其所对应的圆形区域；(3)以步骤( 中特征点的稳定性程度作为筛选依据，对步骤( 中相互间存在重叠的 区域进行优化，得到互不重叠的圆形区域；(4)对每一个局部区域，进行一层非下采样轮廓波变换，并取低频子带内接圆内的系数 作为水印提取所采用的系数；(5)利用中国余数定理进行水印提取，并将提取的水印图像变换到与原始水印图像相 同的大小；(6)计算提取的水印图像与原始水印图像的篡改评价函数，作为水印的相似度度量； 或者根据提取的水印图像对水印存在与否做出判断。
2.根据权利要求1所述的结合多尺度特征与中国余数定理的大容量数字水印方法，其 特征是所述的水印嵌入过程的步骤(3)中所述筛选出局部最稳定的特征点为初步确定水 印嵌入的圆形区域，按如下步骤进行设从原始图像中提取的特征点中特征尺度处于[4， 8]之间的特征点集合为Ω1= {Pl, p2, p3,…，pm}，m为特征点的个数，对Q1中的特征点 P1，以其坐标作为圆心，特征尺度的整数倍为半径生成一个圆形区域札，如果R1中包含有其 它特征点，则将Pji应的高斯差函数值与R1中其它特征点对应的高斯差函数值进行比较，若为局部极大值，则保留P1，并从Q1中去掉存在于R1中的其它特征点，对Ωι中的所有特 征点进行上述操作，最后便得到局部最稳定的特征点，记为Ω2 = {qi, q2，q3，…，qn}，η为 局部最稳定特征点个数。
3.根据权利要求1所述的结合多尺度特征与中国余数定理的大容量数字水印方法，其 特征是所述的水印嵌入过程的步骤(4)中所述的以高斯差函数值作为特征点稳定性强弱 的定量度量，为对相互之间存在重叠的区域进行优化，筛选出互不重叠的区域，按如下步骤 进行(1)对局部最稳定特征点集合Ω2中的特征点qi，q2,q3,…，qn，首先筛选出稳定性最 好的点，设为由Qi与其特征尺度共同决定的圆形区域记为Ri,并将Qi从Ω2中去除，记 入集合ω3;(2)判断Ω2中剩余特征点对应的圆形区域与Ri是否重叠，若重叠则从Ω 2中去掉相应 的特征点，否则保留；(3)对集合Ω2中的特征点，重复步骤(1)和步骤(2)，直到处理完所有的特征点，得到 新的特征点集合ω3;(4)对03中的每一个特征点，结合其特征尺度，最终确定互不重叠的圆形区域，用于水 印嵌入。
4.根据权利要求1所述的结合多尺度特征与中国余数定理的大容量数字水印方法，其 特征是所述的水印嵌入过程的步骤(5)中所述的对每一个区域进行一层非下采样轮廓波 变换，并选取低频子带内接圆内的系数是用于水印嵌入，水印嵌入前，首先将水印图像变换 到与低频子带相同的大小。
5.根据权利要求1所述的结合多尺度特征与中国余数定理的大容量数字水印方法，其 特征是所述的水印嵌入过程的步骤(7)中所述的水印嵌入时利用中国余数定理，且在水 印嵌入时，选取互质的两个整数分别为M1 = 13，M2 = 20，具体嵌入时根据所要嵌入的水印 比特对非下采样轮廓波变换低频子带内接圆内的系数进行修改，强制使得系数对礼、M2的 余数Rei、Re2满足特定的大小关系，即待嵌入水印位为1时，Re1 ^ Re2 ；待嵌入水印位为0 时，Re1 < Re20
全文摘要
一种结合多尺度特征与中国余数定理的大容量数字水印方法，属于用于隐秘通信和数字图像的大容量数字水印方法。该方法利用具有多重不变性的多尺度特征点确定水印嵌入和提取的位置，并将二值水印图像利用中国余数定理在非下采样轮廓波域进行嵌入。嵌入前根据嵌入区域对水印图像进行相应调整，使水印信号自适应于嵌入区域。提取时首先用相同的方法对含水印区域进行定位，然后在非下采样轮廓波域结合中国余数定理提取水印。本发明在嵌入大容量水印信号的同时，含水印图像可以获得很高的质量；提取的水印能抵抗常规的图像处理操作和旋转、缩放、平移等几何变换，且提取的水印与图像所受到的几何变换保持一致，用于隐秘通信和数字图像的版权保护。
文档编号G06T1/00GK102044055SQ20101028235
公开日2011年5月4日 申请日期2010年9月9日 优先权日2010年9月9日
发明者李雷达, 程德强, 钱建生 申请人:中国矿业大学