专利名称：盲水印的三维网格分割方法
技术领域：
本发明涉及一种盲水印的三维网格分割方法。
背景技术：
近十年来，数字水印已经成为版权保护和所有权声明的一种重要方法。大部分可应用的 数字水印技术主要集中于传统的数据媒体类型如音频、图像和视频，他们广泛地分布在互连 网上。直到若干年前，学者们才开始对三维水印进行研究。这是三维对象越来越多地在建筑 设计、虚拟现实、娱乐、机械工程和文化遗产等领域广泛传播的结果。但是，对三维对象的 研究还远没有像对图像和视频的研究那样的深入。
任何水印算法的两个主要要求是水印的鲁棒性和安全性。鲁棒性要求水印能够被检测出
来即使水印模型经历各种有意或无意的变换和攻击。为什么大部份的三维图形对象水印算法 是非盲水印，在检测时需要载体的信息，或者一些算法需要复杂的注册的过程？ 一般而言， 一旦有意或无意的操作发生，如果没有载体的信息的话水印的嵌入和提取将会不同步，结果 导致水印检测失败。然而，除非多携带密钥或一些附加信息，否则在水印的检测过程中安全 需求是实现不了的。这样可以避免多传输或复制载体也能增加安全性。
在大部分的数字水印应用中存在抗剪切或分割的问题。在这种情况下，任何原始相关信 息的丢失都将导致在提取水印吋需要同步恢复丢失的信息才能提取嵌入的水印，在某些情况 下这会导致嵌入的数据会丢失。不仅三维水印需要处理这种问题，许多变换的空域水印算法 也需要处理这种问题。
本发明所述的网格分割算法是以改进三维模型盲水印抵抗剪切、变形等拓扑攻击的鲁棒 性为目的的。有两种方法来解决这种问题。 一种方法是把三维模型分割成若干块，每个水印 位重复的嵌入到那些水印块中，在剪切和简化攻击的情况下，可以从剩下的块中提取水印。 另一种方法是将每个水印位复制c次，也就是同一水印位被重复嵌入c次，这样可以增强水
印的抗攻击性。然而后者存在下述的缺点
1) 计算代价高。如特征值分解或重构水印模型时高阶的Lapkcian矩阵求逆将花费相当 多的时间来执行。
2) 缺少具有显著几何特征的替代网格作为水印目标。
当网格被分割成特殊个数的替代网格的时候，计算时间将大大的减少，并且水印技术可 以利用网格可被开发的完全分开的局部特征。
有效的网格分块是高质量三维水印算法的关键。但是，由于他们会遇到任意的攻击，三维网格分块不是件容易的事。在水印提取阶段，三维网格在没有载体网格的情况下能够被同 步的分割成块，这在三维水印架构里扮演重要的角色。
分割算法把网格分割成特定数量的块，分割将通过以下几种方式之一来实施。他们对水 印应用有着不同的理论贡献。
基于拓扑学的分割技术通过邻接顶点来分割网格。过去几年解决图像分割问题的算法还 很少。当给定网格连通性时，通过邻接关系来对顶点分组并不需要网格的几何信息。Ohbudli 的水印算法从多边形网格的拓扑结构来构造Kirchhoff矩阵。通过修改这个矩阵的光谱系数来 起嵌入水印。 一个简单的网格分割算法虽然不能抵抗网格简化攻击，但也可以在一定程度上 增强抗切割攻击的鲁棒性。
基于几何学的分隔技术是通过位置的相近性来分隔网格独立于网格连接关系。Lee的算 法就是基于几何学分割的例子，在没有连接信息的情况下也能正确计算相邻顶点。算法通过 一个重复的操作来得到若干个中心点并且网格顶点按中心和顶点之间距离最小的原则群集成 块。但是被剪切的网格的块中心点会偏离于原来载体网格的中心点。在遇到剪切攻击的情况 下这种网格分割不稳定。为了进行水印提取和重排列处理，这些块的中心点必须用一个列表 来传输。
以上提到的两种方法，块的数量是人为给定的。大部分基于语意学的网格分块算法，根 据网格区域的同质性自动决定分块的数量。他们利用三维多边形网格几何和拓扑的特性进行 分组的方法来把网格分割成有意义的子网格。到现在为止基于语意学的网格分块方法的水印 应用还很少。
虽然网格分割问题已经受到了很大的关注，但主要的研究集中在为了实现高效的压縮和 编码的三维模型分割方面。将现有的这些分割算法应用于三维水印，存在的主要问题是真实 的分割可能不满足水印的需要，因为三维压缩架构依赖于一些相关的变换参数和剩余的补偿 替代网格之间差距的误差。本发明的目标是只须很少的辅助数据来开发一种能够预测水印网 格分割边界的分割方法。
发明内容
本发明的目的在于提供一种盲水印的三维网格分割方法，该方法不仅在水印嵌入过程中 只需携带少量信息就能够实现对攻击后的模型实施同步分块，而且水印检测过程中不需要原 始模型就能提取水印，只需很少的辅助数据就能够预测水印网格分割边界。
为实现上述之目的，本发明的技术方案是 一种盲水印的三维网格分割方法，其特征在 于水印的嵌入过程包括以下步骤
(1)用NPCA算法对原模型进行处理，并将所有顶点分成复数个集合；(2) 进行网格分割；将所述复数个顶点集合根据平均距离映射成正视图，并从正视图中 选取适当的相素作为附加信息，从而得到附加信息抽样视(3) 嵌入水印信息； 水印的提取过程包括以下步骤-
(a) 用NPCA算法对受攻击的模型进行处理；
(b) 提取受攻击模型的附加信息抽样视(c) 将原模型的附加信息抽样视图与受攻击模型的附加信息抽样视图相减，进行网格对 齐然后估算分割网格；
(d) 提取水印。
本发明的有益效果是通过从载体网格的正视图中选取适当的相素作为附加信息来匹配三 维网格正视图，并不需要原载体网格就能决定载体网格的关键位置，然后划分出稳定的分块 边界。因此，水印嵌入过程中只需携带少量信息就能够实现对攻击后的模型实施同步分块， 水印检测过程中不需要原始模型就能提取水印，是一种只需很少的辅助数据就能够预测水印 网格分割边界的盲水印分割方法。
下面结合附图及具体实施例对本发明作进一步的详细说明。
图1是本发明实施例的流程图。
具体实施例方式
我们以将原始网格和攻击网格分割成跟水印嵌入和检测阶段类似的特殊块为目的。水印 检测过程中除了需要一些附加信息外不需要原始网格。如果水印网格遭受轻微的攻击，原始 网格和水印网格之间的相似部分的同步就比较容易。攻击会导致原始网格和水印网格的同步 分割信息严重地丢失。
本发明的盲水印的三维网格分割方法，水印的嵌入过程包括以下步骤
(1) 用NPCA算法对原模型进行处理，并将所有顶点分成复数个集合；
(2) 进行网格分割；将所述复数个顶点集合根据平均距离映射成正视图，并从正视图中 选取适当的相素作为附加信息，从而得到附加信息抽样视(3) 嵌入水印信息； 水印的提取过程包括以下步骤
(a) 用NPCA算法对受攻击的模型进行处理；
(b) 提取受攻击模型的附加信息抽样视(c) 将原模型的附加信息抽样视图与受攻击模型的附加信息抽样视图相减，进行网格对齐然后估算分割网格； (d)提取水印。
三维网格同步分块的主要步骤如图1所示。下面将详细叙述本发明的减少网格分割误差 问题的方法。
1. NPCA的作用
所有的三维模型包括水印模型存在许多不同的类型，使得很难找到一种对齐方法可以适 用于所以可能的分类。最著名的对齐方法是PCA方法和Karhunen-Loeve变换。PCA算法输 入一些顶点集合得到一些互不相关的特性。模型的顶点集合和模型三角形的质心作为PCA的 输入。CPCA和NPCA协方差计算如下
Ci = [f(Ai十玩+ Ci)+9争十Bi十Cj)/3)]
A '=1
其中，E是三维模型的总面积，Ai，Bi,Ci是三角形Ti的顶点,i^,2，…，N， N是模型三角形 的个数。
&,叫是三角形r,的表面积和单位法向量。f(V)=(V-m)(v-m)T 。这里把m尸丄^^,.Hi做
￡ ,=i
为质心m，而不是mfl^A.fAi+Bi+CO/G做为质心。因为后者对拓扑攻击比较敏感。表l
比较了两种m的选择。敏感性有欧氏距离《nM，m严),《m2,m严)来定义，其中m" , Z=l，2 ,是
被攻击模型的质心。 考虑如下攻击类型
1) 网格不完全切除，面的一边切除掉10%, 20%， 30%;
2) 网格简化，基于边删除的方法，简化30%和50%。
那么单位特征向量d或C2就被计算出来了，他代表模型的主要方向并形成正交的诉3 空间。这些向量按特征值的递增或递减进行分类并用来对齐模型坐标的主要的轴。也就是模 型被平移使得模型质心与坐标原点重合并旋转主要的轴和坐标轴重合。
然而，如果对模型实施拓扑攻击的话，模型的主方向会发生偏离。即使如此，我们也可
以发现PCA主方向的偏离度比NPCA主方向的偏离度大。记《。，/l^，f为d的按顺序递减 的特征值(《)2《、《)〉0),0(/),11(/)是相应的特征向量。同理才)2《^豸)X)是Q的特征值，ui、t4n)，ii^是相应的特征向量。当模型受到拓扑攻击时，CPCA和NPCA协方差矩 阵重新计算，记为C,W和C"。向量《，v『，v^对应于C^的特征值/4^^^/^，特征向 量vf)W)，v^对应于Ct)的W) 2/4"") 2/4n)>0特征值。记a"g(u"),v"))为向量u!s)和v")的夹 角，其中se(c,n)。
通过实验，质心m2和正交坐标C2被攻击后跟d相比受到了轻微的影响。这使得我们 用m2做为模型的质心和NPCA算法来分割模型，这将在下一部分详细叙述。
2. 分割网格
我们的网格分割算法比较适合于把网格分割成2， 4或8块。用方程来表示，这里我们记 诉s空间平面方程ax+by+cz+d-0为p-(a,b，c，d),平面p和向量v的关系记为pdc(p,v)。如果 pdc(p,v)X)，那么v就归到同一组，否则归到其他组。
计算Q的特征向量并且按照特征值递减进行分类。使模型对齐到公认的系统。定义两平 面Pl=(0， 1,0,0)和p2=(0, 0,1,0)。以分割成4部分为例，对齐后模型的所有的顶点Vj简记 为
Gb if pdc(PbVj)〉0 & pdc(p2，Vi)>0 G2， if pdc(ppVi)〉0 & pdc(p2,Vj)《0 G3, if pdc(Pi，Vj)SO & pdc(p2,Vi)>0 G4， if pdc(Pi,Vi)SO & pdc(p2,Vj)S0
其中N是模型顶点数。
3. 攻击模型的分割边界评估
总的来说，盲水印算法比较实用，因为水印检测过程不需要原始模型的信息。然而众所 周知的是盲水印算法的鲁棒性比非盲水印的鲁棒性差，在鲁棒性方面圩多优化算法取得了很 大的进步。在这一部分将详细的叙述在水印检测阶段如何同步地分割在水印嵌入阶段只嵌入 少量附加信息的网格。
实际上，对于网格分割提出了一个新的图像匹配分割算法，简称EM。首先，得到三维被 攻击的模型的公认位置的视图。这个视图由三维模型面到质心迈距离的二维灰度图像来表示。 那么方块内平均错误的方法被用來评估带有辅助信息的原始模型的公认位置参数，也可以用 来评估被攻击的模型公认位置参数。
最后，执行一个三维模型的逆相似变换，这样分割边界就画出来了。详细过程如下 首先将被攻击的模型变换到公认位置，它的表面被绘制成视图。开始，我们设置 《=(27' + 1); /(45)，化-2^/(2S)，/,hO，l，…,25-l其中5是抽样率。所有顶点vi=(xi,yi,zi),^=1,2，...^)被转换到球面坐标
《=Kv,) = VV+x2+z,2 ，
《=e(v,) = arctan(^,2 +z,2 )/乂, "=p(v,.) = arctan(z, /;c,). 然后我们定义O(y，A)和craA:)如下 Ca) = {v, i & S<S < .
0ca)，其中i-i是模操作.如果0(M)为空，那么crca"o.
然后，范围在
的ol/，W和二维灰度图像也叫做视图，从xoy-平面前的一个视点得
到。也就是，这个视图把平均距离信息ol/,W做为xoy平面上的投影。因为xoy-平面前的信
息和xoy-平面后的信息重叠在xoy-平面上，除前者之外后者被拼凑成一块。为了方便起见， 我们把攻击模型的视图叫做攻击视图，也由矩阵M,d 来表示。
我们使用EM算法来决定分块边界。设^…^为选择的像素集合，他们分散在载体网格视 图里，将被用来识别攻击模型的分割边界在分割阶段。有很多方法来确定这些像素，例如选 择满足/mod8-4条件的像素。按照经验我们设置4^/6""4^/32为附加信息量。
EM算法搜索攻击视图矩阵和^…^的平移距离和旋转角度使得方块平均错误最小
,!X+W6- ,]2 (2)
其中",厶)是^在载体网格视图中的位置并且(^， )是沿水平和垂直方向偏离",/J
的值。很显然，攻击视图可以通过平移和旋转载体网格得到。所以方程(2)等价于-
廳S[°>i+ i",r),p2+ 2(w) 一&]2
，ipos(&)+jposfe).tan(力 力,pos(^)-ipos(^)-tan(力 其中化，化是/W'i的偏离象素，函数ipos(仏)，jpos(&)分别代表《,的水平和垂直位置， 这里一个格子表示一个象素。我们得到平移和旋转的度数是轻微的这个结论。整数外，P2的搜索范围限定在[^/8，B/8]之间并且rm/(2^),选择象-^/4^^5/4这样小的整数，平均方块误 差的最小值的计算时间也随着减少。
而且，我们需要通过攻击模型的质心m来估算原始网格的质心。质心修正的过程将改进 分割的精度，因为我们的网格分割算法在模型对齐到公认的相关框架之后执行的。在平移距 离^,P2和旋转角度/W基础上，我们定义增量m那么
最后，实施旋转和平移变换并且所有的顶点根据方程l的方法分成若干块。平移使得m' 与坐标原点一致。那么模型里的所有顶点通过下列式子得到
其中y^/^兀/5,WziY(ca)是绕y轴或z轴弧度为co的旋转矩阵。更多不规则的块可以通过指定的相 关路径比如Dijkstra路径来获得。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超 出本发明技术方案的范围时，均属于本发明的保护范围。
权利要求
1、一种盲水印的三维网格分割方法，其特征在于水印的嵌入过程包括以下步骤(1)用NPCA算法对原模型进行处理，并将所有顶点分成复数个集合；(2)进行网格分割；将所述复数个顶点集合映射成映射视图，并从映射视图中选取适当的相素作为附加信息，从而得到附加信息抽样视图；(3)嵌入水印信息；水印的提取过程包括以下步骤(a)用NPCA算法对受攻击的模型进行处理；(b)提取受攻击模型的附加信息抽样视图；(c)将原模型的附加信息抽样视图与受攻击模型的附加信息抽样视图相减，进行网格对齐然后估算分割网格；(d)提取水印。
2、根据权利要求1所述的盲水印的三维网格分割方法，其特征在于在步骤(2)中，将所 述复数个顶点集合根据平均距离映射成正视图，并从正视图中选取适当的相素作为附加信息。
全文摘要
本发明涉及一种盲水印的三维网格分割方法，水印的嵌入过程包括以下步骤1)用NPCA算法对原模型进行处理，并将所有顶点分成复数个集合；2)进行网格分割；将顶点集合映射成正视图，并从正视图中选取适当的相素得到附加信息抽样视图；3)嵌入水印信息；水印的提取过程包括以下步骤a)用NPCA算法对受攻击的模型进行处理；b)提取受攻击模型的附加信息抽样视图；c)将原模型的附加信息抽样视图与受攻击模型的附加信息抽样视图相减，进行网格对齐然后估算分割网格；d)提取水印。该方法只需携带少量信息就能够实现对攻击后的模型实施同步分块，且不需要原始模型就能提取水印，只需很少的辅助数据就能够预测水印网格分割边界。
文档编号G06T1/00GK101533509SQ20091011135
公开日2009年9月16日 申请日期2009年3月23日 优先权日2009年3月23日
发明者姚志强, 孔祥增, 潘日晶, 陈丽萍 申请人:福建师范大学