CN1266587A

CN1266587A - 水印检测的方法与装置

Info

Publication number: CN1266587A
Application number: CN99800689A
Authority: CN
Inventors: J·A·海特斯马; A·A·C·M·卡尔克; A·J·M·德尼森
Original assignee: Koninklijke Philips Electronics NV
Current assignee: Civolution BV
Priority date: 1998-03-04
Filing date: 1999-03-02
Publication date: 2000-09-13
Anticipated expiration: 2019-03-02
Also published as: JP2001525152A; WO1999045705A3; EP0981903B1; CN1266586A; ES2297919T3; KR20010012250A; JP4428728B2; JP2001525151A; ID23884A; DE69940375D1; KR20010012214A; EP0981903A2; AU2437499A; WO1999045707A3; RU2249308C2; US20030086587A1; CN1166187C; JP2002503431A; WO1999045705A2; US7191334B1

Abstract

本发明公开了一种用于检测在一个信号(例如图象或视频)信号中的水印的改进方法和装置。通过在检测在所怀疑的图象(q)和将被检测的水印(Wi)之间的相关性数值之前将所述信号经过对称仅相位匹配滤波(24—28,30)而使得检测更可靠和较不脆弱。

Description

水印检测的方法与装置

发明领域

本发明涉及一种用于检测在一个信息信号中的水印的方法，该信息信号可能已经通过根据一个水印图案的(时域或空间)相应值修改所述信息信号值来附加了水印。本发明还涉及一个用于检测水印的装置。

发明背景

在国际专利申请WO-A-98/03014中公开了一个在开头段中定义的现有方法。水印通过计算所怀疑的信息信号与一个提供的水印图案的相关值，并将该相关值与一个预定阈值相比较来检测。如果该相关值大于阈值，则说水印存在，否则说水印不存在。该相关值越大，则检测越可靠，可以允许更多的处理直至不能检测出更多的水印。如WO-A-98/03014所公开的，该信息信号和/或水印图案可以在计算相关值之前通过匹配滤波。从而可以显著地改善水印检测的可靠性。

现有的水印检测方法的一个问题是相对于信息信号而言，水印图案的(时域或空域)位置并不是完全知道的。如果提供给检测器的水印位置与嵌入时的位置不同，则相关值将很小，而可能错误地检测为水印不存在。

本发明的目的和概述

本发明的一个目的是进一步改善检测一个水印的方法。

该目的通过权利要求1中所限定的方法和权利要求6中所限定的装置来实现。而在从属权利要求中限定了可取的实施例。

本发明认识到对于水印的大量可能位置而言，信息信号和所提供的水印的相关值最好在傅立叶域计算，检测的稳健性和可靠性可以通过在相关之前提供对称相位只匹配滤波(SPOMF)给信息信号和水印来改善。SPOMF，它可在图案识别领域中已知，假定傅立叶系数的相位中携带了相关性检测所需的多数有关信息。由此，复傅立叶系数的幅度被规范化以具有基本相同的幅度。

虽然在SPOMF的表达中的形容词“对称”指同时对信息信号和水印进行滤波，本发明人已发现将傅立叶变换后的信息信号规范化对相关性检测的改善贡献最大。对傅立叶变换后的水印的幅度规范化构成本发明的一个实施例。

规范化复傅立叶系数的幅度的步骤包括将每一个傅立叶系数除以其绝对值。然而，本发明人还发现如果将每一个系数除以其实部或虚部(根据其中谁最大)不会使检测强度显著降低。这对于规范化(该幅度可以在1和

之间变化)的一个近似，但相当大地减少了计算量。

附图简述

图1示意性显示了一个用于将一个水印嵌入信号的装置。

图2和3显示了一个图1所示的嵌入器的操作的示意图。

图4示意性显示了一个用于检测该嵌入的水印的装置。

图5，6A，和6B显示了一个图4所示的检测器的操作的示意图。

图7显示了一个用于回放一个具有嵌入水印的视频比特流的设备。

图8示意性地显示了一个用于根据本发明检测该嵌入水印的装置的实施例。

图9A和9B显示了图8所示的检测器的操作示意图。

图10示意性显示了根据本发明用于检测该嵌入水印的装置的另一个实施例。

优选实施例详述

为方便起见，根据本发明的嵌入水印机制将用一个用于将不可见标记附在视频内容上的系统来描述，但所述原理明显也可用于其他内容，包括音频和多媒体。我们将在下文中将该方法称为JAWS(只是另一个嵌入水印系统)。

图1显示了一个用于提供背景信息的水印嵌入器的实用实施例。该嵌入器包括一个产生一幅图象的图象源11，和将水印W加到该图象P上的加法器12。该水印W是一个与图象大小一样的噪声图案，例如水平为N1个象素而垂直为N2个象素。该水印W代表了一个密钥K，即将在接收端提取的多-比特码。

为避免为在大的N1×N2空间上搜索水印W所需的水印检测过程，可以在图象的范围内通过重复，和如果需要，切断被称为“瓦片”(“tile”)W(K)的较小单元来产生水印。瓦片W(K)具有一个固定的尺寸M×M。该瓦片尺寸M不应太小：较小的M意味着在W(K)中的对称，并因此有更大的安全隐患。另一方面，M不应太大：一个大的M值意味着对于检测器的大的搜索空间，并因此带来更大的复杂性。在JAWS中，我们选择M＝128作为折衷。

接着计算(16)局部深度映射图或可见度屏障λ(P)。在每一象素位置，λ(P)提供了一个用于附加噪声的量度。映射图λ(P)被构建为具有等于1的平均值。随后用λ(P)调制(17)扩展序列W(K)，即在每一位置平铺的水印W(K)乘以在该位置的λ(P)的可见度值。因此所得的噪声序列W(K，P)取决于密钥K和P的图象内容。我们成W(P，K)为适应性水印，因为它适应于图象P。

最后，最终水印的强度由一个全局深度参量d决定，d提供了一个W(K，P)的全局比例(18)。一个大的d值对应于一个稳健但可能可见的水印。一个小的d值对应于一个几乎不可感知但不稳健的水印。d的实际选择可以是在稳健性与感知度要求之间的折衷。水印图象Q通过将W＝d×W(K，P)与P相加(12)，四舍五入为了整象素值并限制到所允许的象素值范围中来获得。

为了将多位码K嵌入水印W中，可以由有限组的非相关基本或原始瓦片{W₁...W_n}和其移动版本根据以下公式来构建每一瓦片W(K)，

W (K) = \underset{i, j}{Σ} s_{i_{j}} shift (W_{i}, k_{i_{j}})

这里“shift(W_i，

)”代表了一个基本M*M瓦片W_i绕一个矢量周期性旋转的一个空域移动。符号s∈{-1，+1}而移动k通过一个编码函数E(13)而取决与密钥K。检测器的任务是在检索到符号s_I和移动k_j后重构K。注意，每一基本瓦片W_i可以产生几次。在图1中，编码器13产生W(K)＝W₁+W₂-W₂’，其中W₂’是W₂的移动版本。图3图示说明了这一操作过程。

图4显示了一个水印检测器的原理示意图。水印检测器接收可能嵌入水印的图象Q。在JAWS中不是对每一单独帧，而是对帧组进行水印检测。通过将许多帧累加(21)，可以改善检测的统计特性并因此改进检测的可靠性。随后，累加的帧分离为(22)一些大小为M×M(M＝128)的块，而所有的块堆积在(23)一个大小为M×M的缓存器q中。该操作被称为折叠。图5示出了这一折叠过程。

检测过程中的下一步骤是确定在缓存器q中一个特定噪声图案的存在。为了检测该缓存器q中是否包括一个特定的水印图案W，将缓存器内容与所述水印图案进行相关。计算一个怀疑的信息信号q与一个水印图案w之间的相关性包括计算该信息信号值与该水印图案的相应值的内积d＝<q，w>。对应一个一维信息信号q＝{q_n}和水印图案w＝{w_n}，所述计算过程可以写为以下数学形式：

d = \frac{1}{N} Σ_{n = 1}^{N} q_{n} w_{n}

对于二维M×M图象q＝{q_ij}而水印图案W＝{w_ij}，内积为：

d = \frac{1}{M^{2}} Σ_{i = 1}^{M} Σ_{j = 1}^{M} q_{ij} w_{ij}

理论上，瓦片W_i绕其移动的矢量k_i可以通过连续地将W_i和不同的矢量k提供给检测器，而确定对应哪一个k，其相关性最大。然而，这种强制搜索算法太浪费时间。此外，图象W在水印检测之前可能经历了各种形式的处理过程(诸如移动或剪裁)，这样检测器将不知道基本水印图案W_i相对于图象的空间位置。

替换强制搜索，JAWS利用了图案W(K)的结构。为确定这些原始图案，其符号和移动的存在而检测缓存器q。图象q与一个通过一个矢量k(水平方向为k_x象素而垂直方向为k_y象素)而移动的原始图案w的相关度d_k为：

d_{k} = \frac{1}{M^{2}} Σ_{i = 1}^{M} Σ_{j = 1}^{M} q_{ij} w_{i + k_{x}, j + k_{y}}

使用快速傅立叶变换同时计算一个基本图案W_I的所有可能移动矢量k的相关值d_k。如图4所示，缓存器q的内容和该基本水印图案W_i分别在变换电路24和25中都经历了一次快速傅立叶变换(FFT)。这些操作过程如下：

\hat{q} = FFT (q)

和

\hat{w} = FFT (w)

这里和

为复数组。

计算相关值类似于计算q与W_I的共轭的卷积。在变换域，这对应于如下公式：

\hat{d} = \hat{q} &CircleTimes; conj (\hat{w})

这里符号表示逐点相乘，而conj()表示将变量的虚部的符号反相。在图4中，

的共轭由共轭电路26实施，而逐点相乘由乘法器27实施。现在通过将所述乘法的结果进行逆傅立叶变换而获得该组相关值d＝{d_k}：

d = IFFT (\hat{d})

该运算在图4中由逆FFT电路28实施。如果一个相关值d_k大于一个给定阈值，则检测存在一个水印图案W_i。

图6A显示了如果在图象Q中检测到一个水印图案W₁的存在时的相关值d_k的曲线。峰61指示确实发现W₁。该峰的位置(0，0)显示提供给检测器的图案W₁相对于图象Q的空域位置碰巧与提供给嵌入器的图案W₁相对于图象Q的空域位置相同。图6B显示了如果水印W₂提供给检测器时的相关值曲线。现在发现两个峰。在(0，0)的正峰62表示了水印W₂的存在，而在(48，80)的负峰63表示了水印-W₂’的存在。后一个峰63与峰62(或，类似的，峰61)的相对位置揭示了W₂’相对W₂的相对位置(以象素为单位)，即位移矢量k。所嵌入数据K由所发现的矢量而获得。

所嵌入的信息可以标识诸如版权所有者或内容的描述。在DVD的防止复制机制中，它允许材料被标注为“只复制一次”，“未复制”，“无限制”，“不再复制”，等等。图7示出了一个用于回放记录在光盘71中的MPEG比特流的DVD驱动器。所记录的信号经一个开关72提供给一个输出端子73。输出端子连接到一个外部MPEG解码器和显示设备(未示出)。假定该DVD驱动器不可回放具有预定嵌入水印的视频信号，除非其他与本发明不相关的条件实现。例如，嵌入水印的信号只在光盘71包括一给定“摇摆式”密钥时才可被回放。为了检测水印，该DVD驱动器包括一个如上所述的水印检测器74。该检测器接收所记录的信号并响应于水印是否检测到来控制开关72。

傅立叶系数

为复数，也即，它们具有一个实部和一个虚部，或一个幅度和一个相位。本发明人已发现如果去掉该幅度信息而只考虑该相位信息，则检测器的可靠性将显著提高。图8示出了根据本发明的检测器的相关电路的一个实施例。该实施例不同于图4所示的实施例，其中在乘法器27和逆傅立叶变换电路30之间插入了一个幅度规范化电路30。规范化电路的操作包括逐点地将每一系数除以它的幅度。数学公式为：

\hat{d} = \hat{d} Φabs (\hat{d}) - - - - (1)

这里Φ表示逐点相除而abs()表示：

abs ({\hat{d}}_{k}) = \sqrt{R {({\hat{d}}_{k})}^{2} + I {({\hat{d}}_{k})}^{2}} - - - - (2)

这里R()和I()分别表示该变量的实部和虚部。

幅度的所述规范化被称为对称的仅相位匹配滤波(SPOMF)。图9A和9B显示了SPOMF相关的效果。更具体地，图9A示出了当使用线性相关，即没有幅度规范化电路时的系数值d00。在整个标准偏差矩阵的单元中表示的系数值d00的数值为9.79。图9B显示了当使用SPOMF相关时的系数值。该系数值d00现在为标准偏差的62.77倍。可取的是，图9B中的峰可以比图9A中的峰更可靠地检测。

因为对

的幅度的规范化等效于对

和

的幅度的规范化，所以图8中的规范化电路30可以由在FFT电路24和25之后的两个规范化电路来代替。然而，由于所嵌入的水印实际上是伪随机噪声图案，其中每一样本都是独立和相同地从一个标准分布中提取，所以它将已是一个合理的白(平坦)频谱。由此来看，只对信息信号幅度的规范化已被发现是足够的。图10显示了这样一个实施例。该幅度规范化电路30现在位于FFT电路24与乘法器27之间。在这一实施例中，

的幅度并非精确，但基本上是相同的。

还可进一步注意到所提供的水印Wi的FFT和共轭(例如在图4，8和10中分别为电路25和26)以及的幅度的可选性规范化可以预先计算并存储在存储器中。

本发明人已进一步发现如果对于规范化运算使用一个更进一步的近似，检测强度不会显著降低。除了如公式(1)和(2)所限定的，将每一傅立叶系数除以其幅度以外，将每一系数根据哪一个大而除以其虚部或实部就足够了。从而，该幅度将在1和

之间变化，而检测强度只降低1％。优点是降低了计算复杂度，因为在公式(2)中的平方根不需要计算了。在这个实施例中，规范化电路30执行以下算法：

{if}^{R ({\hat{d}}_{k}) &GreaterEqual; I ({\hat{d}}_{k})}

the n^{{\hat{d}}_{k} : = {\hat{d}}_{k} / R ({\hat{d}}_{k})}

els e^{{\hat{d}}_{k} : = {\hat{d}}_{k} / I ({\hat{d}}_{k})}

这里

为

(对于其他情形，或是

或

)的复傅立叶系数，和I(

分别是

的实部和虚部，而/表示除法。

如果规一化运算只包括2为幂的除法，则可以进一步降低复杂度。则该幅度将在1和

之间变化，该检测强度将降低2.5％。以2为幂的除法在大多数通用浮点数系统中特别有效，其中复系数的实部和虚部表示为：

sign·mantissa·2^exp这里mantissa是一个在1与2之间的数的定点表示。上述规范化运算现在可以简化如下：

{if}^{\exp (R ({\hat{d}}_{k})) &GreaterEqual; \exp (I ({\hat{d}}_{k}))}

then f^{div ider : = \exp (R ({\hat{d}}_{k}))}

{else}^{div ider : = \exp (I ({\hat{d}}_{k}))}

\exp (R ({\hat{d}}_{k})) = \exp (R ({\hat{d}}_{k})) - div ider

\exp (I ({\hat{d}}_{k})) = \exp (I ({\hat{d}}_{k})) - div ider

总之，本发明公开了一种用于检测在一个信息(例如图象或视频)信号中的水印的改进方法和装置。通过在检测在所怀疑的图象(q)和将被检测的水印(Wi)之间的相关性数值之前将所述信号经过对称仅相位匹配滤波(24-28，30)而使得检测更可靠和较不脆弱。

Claims

1.一种用于检测在一个信息信号中的水印的方法，包括以下步骤：

-计算一个表示所述信息信号与所述水印的相关性的值；和

-检测所述相关性值是否大于一个预定阈值；

其特征在于

-所述计算过程包括将该信息信号和水印傅立叶变换为各自的傅立叶系数，将相应的系数相乘以获得表示相关性值的傅立叶系数，以及对表示该相关性值的傅立叶系数逆变换的步骤，每一个表示相关性值的傅立叶系数相应于该水印相对于该信息信号的多个位置中的一个，所述检测步骤包括检测是否有至少一个所述相关性值大于预定阈值；该方法还进一步包括将表示该相关性值的傅立叶系数规范化为基本相同幅度的步骤。

2.权利要求1中所述的方法，其中，所述规范化步骤被用于表示该信息信号的傅立叶系数。

3.权利要求2中所述的方法，其中，所述规范化步骤还进一步用于表示该水印的傅立叶系数。

4.权利要求1中所述的方法，其中，所述规范化步骤包括将每一傅立叶系数除以其自身幅度。

5.权利要求1中所述的方法，其中，所述规范化步骤包括将每一傅立叶系数根据其实部和虚部哪一个大而除以其实部或虚部。

6.权利要求5中所述的方法，其中，所述除法包括以2为幂的除法。

7.一种用于检测在一个信息信号中的水印的装置，包括：

-用于计算一个表示所述信息信号与所述水印的相关性的值

(d)的装置(24-28)；和

-用于检测所述相关性值是否大于一个预定阈值的装置(29)；

其特征在于

-所述计算装置包括用于将该信息信号和水印傅立叶变换为各自的傅立叶系数的装置(24，25)，用于将相应的系数相乘以获得表示相关性值的傅立叶系数的装置(27)，以及用于对表示该相关性值的傅立叶系数逆变换的装置(28)，每一个表示相关性值的傅立叶系数相应于该水印相对于该信息信号的多个位置中的一个；

-所述检测装置(29)被配置来检测是否有至少一个所述相关性值大于预定阈值；

-该装置还进一步包括用于将表示该相关性值的傅立叶系数规

范化为基本相同幅度的装置(30)。

8.一种用于记录和/回放一个信息信号的设备，包括一个用于根据在所述信号中存在一个水印而使得不能记录和/回放该视频信号的装置(72)，其特征在于，该设备还包括一个用于如权利要求7所述的检测在该信号中所述水印的装置(74)。