1.一种抗几何攻击的立体图像零水印嵌入与提取方法，其特征在于，包括：在原始彩色立体图像中嵌入零水印和对待认证合法性的彩色立体图像进行零水印提取；

所述在原始彩色立体图像中嵌入零水印，包括：

采用基于矩阵变换的图像置乱方法对原始二值版权水印W进行置乱，得到一次置乱后的水印图像W1，采用基于完备拉丁方的图像置乱方法对水印图像W1进行再次置乱，得到二次置乱后的水印图像W2，利用二维Chebyshev‑Singer Map系统对二次置乱后的水印图像W2进行混沌加密，得到双重置乱及加密处理后的水印图像WE；

对原始立体图像I左、右视点各自的R、G、B通道进行分离，将分离后的六个分量进行联合并构造一个三阶张量ELR；对张量ELR进行Tucker分解，获得原始立体图像I的第一能量图EM；

对EM进行非重叠分块并计算每个子块的均值，由所有子块的均值构造EM的均值子图EMS；

计算均值子图EMS的分数阶雅克比‑傅里叶矩，并生成其混合低阶矩特征VM；

对混合低阶矩特征VM按照其均值进行二值化，得到二值序列H，再对二值序列H进行升维变换，得到二值过渡矩阵P1，对过渡矩阵P1冗余化和扩充，获得二值鲁棒特征矩阵T；

对二值鲁棒特征矩阵T进行混沌排序置乱，得到置乱后的特征矩阵TS；

将双重置乱和加密处理后的水印图像WE和排序置乱后的特征矩阵TS执行异或运算，得到最终的认证零水印WZ，将WZ保存到注册机构的水印数据库中，并将零水印嵌入过程中的密钥保存，即完成零水印的嵌入过程；

所述对待认证合法性的彩色立体图像进行零水印提取，包括：对待认证立体图像I′左、右视点各自的R、G、B通道进行分离，将分离后的六个分量进行联合并构造一个三阶张量E′LR；然后对张量E′LR进行Tucker分解，获得待认证立体图像I′的第一能量图E′M；

对E′M进行非重叠分块并计算每个子块的均值，然后由所有子块的均值构造E′M的均值子图E′MS；

计算均值子图E′MS的分数阶雅克比‑傅里叶矩，并生成其混合低阶矩特征V′M；

对混合低阶矩特征V′M按照其均值进行二值化，得到二值序列H′，再对二值序列H′进行升维变换，得到二值过渡矩阵P′1，对过渡矩阵P′1冗余化和扩充，获得二值鲁棒特征矩阵T′；

对二值鲁棒特征矩阵T′进行混沌排序置乱，得到置乱后的特征矩阵T′S；

取出保存在注册机构版权鉴别数据库中的认证零水印WZ，与置乱后的特征矩阵T′S进行异或操作，得到未解密的二值水印图像W′E；

对未解密水印图像W′E依次进行混沌解密和基于完备拉丁方、矩阵变换的先后两次反置乱操作，即可提取出最终的水印信息W′，最后根据W′显示的内容信息来鉴别待认证的彩色立体图像I′的版权归属；

所述获得原始立体图像I的第一能量图EM，包括：

将一组大小为M×N的双视点彩色立体图像I的左、右视点各自的R、G、B通道进行分离，将分离后的共六个分量联合并视作一个三阶张量ELR，对其进行Tucker分解，则立体图像I可表示如下，，

(1) (2) (3)

式中ELR表示张量形式下的立体图像I，表示核心张量，V 、V 、V 是3个大小分别为M×M、N×N和6×6的因子矩阵，设张量 为张量ELR的一个子张量，且表示子张量 的第i层，则第1层子图Z1即为立体图像Tucker分解后的第一能量图EM。

2.根据权利要求1所述的方法，其特征在于，所述利用二维Chebyshev‑Singer Map系统对二次置乱后的水印图像W2进行混沌加密，包括：所述二维Chebyshev‑Singer Map系统，采用如下公式表达：，

式中，a,b,c,d均为系统固定参数，τ和μ为系统控制参数，xi和yi为系统迭代过程中的中间变量，xi+1和yi+1分别为xi和yi的下一个状态；

采用二维Chebyshev‑Singer Map系统产生一个长度为原始版权水印W两倍大小的混沌随机序列Y1，截取随机序列Y1的后半段，并按照其序列均值进行二值化，对二值化后的序列进行升维变换，得到与原始版权水印W相同大小的二维混沌矩阵，将所述二维混沌矩阵与二次置乱后的水印图像W2执行异或运算，得到加密处理后的水印图像WE。

3.根据权利要求1所述的方法，其特征在于，所述计算均值子图EMS的分数阶雅克比‑傅里叶矩，并生成其混合低阶矩特征VM，包括：将均值子图EMS转换到极坐标下，获得极坐标图像 ，根据下式计算 的分数阶雅克比‑傅里叶矩Fnm，，

+

式中，分数参数α∈R ；阶数n∈N；重复度m∈Z；参数p,q∈R且满足p‑q>‑1，q>0；径向基函数Jn(α,p,q,r)可表示为 ，为权函数， 为归一化常

数， 为分数阶Jacobi多项式， 为伽

马函数；

由Fnm生成混合低阶矩特征VM，通过如下方式实现：

,

其中，n表示阶数，m表示重复度，NW表示水印图像尺寸，分数参数α∈{0.25,0.5,1,2,4}，符号⌈•⌉表示向上取整；

所述对过渡矩阵P1冗余化和扩充，获得二值鲁棒特征矩阵T的方法为：对过渡矩阵P1冗余化，得到矩阵 ，将矩阵P2均分成左上、右上、左下、右下4个区域块，分别记作B1、B2、B3和B4，分别用 替换B1、B2、B3、B4区域块中的0元素，分别用 替换B1、B2、B3、B4区域块中的1元素，得到扩充后的矩阵即为二值鲁棒特征矩阵T。

4.根据权利要求1所述的方法，其特征在于，所述采用基于矩阵变换的图像置乱方法对原始二值版权水印W进行置乱，得到一次置乱后的水印图像W1，包括：利用矩阵变换公式对原始版权水印W进行第一次置乱，获得第一次置乱后的水印图像W1，所述矩阵变换公式表达为：，

式中，A表示置乱矩阵，且满足det(A)=±1， ，NW为原始版权水印图像W的尺寸大小；(x,y)和(x′,y′)分别表示原始版权水印图像像素经过矩阵变换前后的坐标，e、f为控制参数。

5.根据权利要求1所述的方法，其特征在于，所述采用基于完备拉丁方的图像置乱方法对水印图像W1进行再次置乱，得到二次置乱后的水印图像W2，包括：构造NW阶完备拉丁方矩阵C1，再将矩阵C1与一个元素全为1的NW阶矩阵相加得到矩阵C2；

2

将矩阵C2进行列扩展，构造出一个含有NW 个互异的有序数对的置乱矩阵S，利用置乱矩阵S提供的位置信息，对一次置乱后的水印图像W1进行像素位置置乱，获得二次置乱后的水印图像W2。

6.根据权利要求1所述的方法，其特征在于，所述对二值鲁棒特征矩阵T进行混沌排序置乱，得到置乱后的特征矩阵TS，包括：将特征矩阵T通过Z字形扫描成行向量T1，利用二维Chebyshev‑Singer Map系统生成长度与原始水印W大小相同的随机序列Y2，然后对序列Y2进行升序排序，并将排序后的位置索引记为U，最后基于索引U对行向量T1进行重排序，再将重排序后的结果进行Z字形逆扫描得到置乱后的特征矩阵TS。

7.根据权利要求1所述的方法，其特征在于，所述对过渡矩阵P′1冗余化和扩充，获得二值鲁棒特征矩阵T′，包括：对过渡矩阵P′1冗余化，得到矩阵 ，将矩阵P′2均分成左上、右上、左下、右下4个区域块，分别记作B′1、B′2、B′3和B′4，分别用 替换B′1、B′2、B′3和B′4区域块中的0元素，分别用 替换B′1、B′2、B′3和B′4区域块中的1元素，得到填充后的矩阵即为二值鲁棒特征矩阵T′。

8.根据权利要求1所述的方法，其特征在于，所述密钥包括：将对原始二值图像水印图像W进行矩阵变换置乱时的控制参数e、f和置乱次数作为密钥Key1；

将对一次置乱后的水印图像W1进行完备拉丁方置乱时的置乱矩阵S和置乱次数作为密钥Key2；

将采用二维Chebyshev‑Singer Map系统对二次置乱后的水印图像W2进行混沌加密生成随机序列Y1时使用的系统固定参数a,b,c,d和控制参数τ、μ作为密钥Key3；

将采用Chebyshev‑Singer Map系统对二值鲁棒特征矩阵T和T′进行混沌排序置乱生成随机序列Y2时使用的系统固定参数a,b,c,d和控制参数τ、μ作为密钥Key4。