基于LPM和双谱分析的抗旋转攻击零水印算法

何 冰

(1.渭南师范学院 数理学院，陕西 渭南 714099;2.陕西省X射线检测与应用研究开发中心,陕西 渭南 714099)



【现代应用技术研究】

基于LPM和双谱分析的抗旋转攻击零水印算法

何 冰1,2

(1.渭南师范学院 数理学院，陕西 渭南 714099;2.陕西省X射线检测与应用研究开发中心,陕西 渭南 714099)

零数字水印算法因其具有一定的鲁棒性和良好的不可见性，在版权保护领域被广泛应用。然而已有的零数字水印方案抵抗几何变换攻击的性能不强，假如将宿主图像旋转一定小的角度就容易导致数字水印在检测时失败。为了提升零数字水印抵抗旋转几何攻击的能力，提出了一种LPM(对数极坐标)和双谱分析相结合的抵抗旋转几何变换的零数字水印方法。先将二维载体宿主图像由笛卡尔坐标系转换到对数极坐标系中，再对转换后的图像进行双谱变换；然后计算双谱分析后的绝对值矩阵的代表特征：均方差、均值、主对角线元素之和及能量，并将其组合后得到特征向量；最后利用这些特征向量来设计和构建零数字水印。仿真实验表明，该算法不但能够获得较好的视觉效果，而且对于旋转几何变换攻击具有较强的稳健性，同时对于非几何攻击信号处理，如均值和中值滤波、JPEG有损压缩、裁切攻击等均具有一定的稳健性。

对数极坐标变换;零数字水印;稳健性;双谱分析

0 引言

网络技术和人工智能技术的快速发展，在便利民众的同时，也带来了一些负面问题：版权认证的归属问题、交易过程中防篡改复制问题、保护信息媒体数据原创性问题等都成为工业4.0时代下急需解决的问题。数字水印技术作为解决这些问题的一种有效方法，已成为信息安全领域的研究热点。国内数字水印领域专家谭铁牛院士、长江学者杨义先教授等都为数字水印技术在版权保护、防篡改性上作出了重要贡献。但现有的数字水印技术因其自身的局限性，并未在电商和微商时代的背景下取得广泛的实际应用。

目前，零水印技术仍处在探索及发展阶段，2016年12月我们在国外多个大型学术数据库以“zero-watermarking”为主题词进行检索，在2011—2016年间共检索出50篇相关论文，其中Science Direct(elsevier)数据库检索到1篇，Springerlink数据库检索到8篇，SPIE数据库检索到1篇，IEEE Xplore Digital Library 数据库检索到3篇，会议论文37篇；国内在中国知网数据库(CNKI)以“零水印”为主题进行检索，共检索出513篇相关学术论文。统计和分析以上研究成果，说明关于零水印的研究在国内处于研究发展阶段，而在国外尚属起步，但该方向的研究已引起相关领域专家、学者的关注。零水印算法的核心是提取能够表征图像本质特性的图像特征信息，然后和版权水印图像绑定生成零水印，最后将零水印在版权认证数据库中(IPR)注册保存。因此有关零水印的研究可以转换为如何有效提取原始宿主图像特征问题。目前研究者在零水印图像特征的构造问题上已提供了很多思路和解决方法，这些方法可以归结为两大类：(1)基于空间域的零水印图像特征构造方法[1-2]；(2)基于图像变换(变换域)的零水印图像特征构造方法[3-10]。在变换域中，针对旋转几何攻击，本文提出一种基于LPM和双谱分析的抗旋转攻击零数字水印算法，可以有效抵抗旋转造成的几何变换。

1 对数极坐标变换(Log Polar Mapping，LPM)

令二维函数I(x,y)表示笛卡儿坐标(直角坐标系)(x,y)上的灰度图像，该灰度图像顺时针旋转φ角度后可以表示为：

Ir(x,y)=I(xcosφ+ysinφ,ycosφ-xsinφ)。

(1)

在笛卡儿坐标系(直角坐标系)下，令：

x=eρcosθ，y=eρsinθ，0≤θ<2π，ρ∈R2，

(2)

将(2)式代入(1)式得：

Ir(eρcosθ,eρsinθ) =I(eρcosθcosφ+eρsinθsinφ,eρsinθcosφ-eρcosθsinφ)

=I(eρcos(θ-φ),eρsin(θ-φ))。

(3)

简记为：

Ir(ρ,θ)=I(ρ,θ-φ)。

(4)

其中：ρ、θ分别表示极轴和极角的采样点的个数(对数极坐标系下的二维坐标)。

由(4)式可知灰度图像在笛卡儿坐标系(直角坐标系)下若将图像旋转角度φ，则该图像会在对数极坐标系下沿着极角θ方向平移φ个单位。

2 双谱分析

定义1 某平稳随机时间序列S(k)的三阶累积量可以表示为：

LSSS(m,n)=E[S(k)·S(k+n)·S(k+m)]，

(5)

则其双谱密度为：

(6)

其中：双谱密度Ds(ω1,ω2)一般为复数表示，即它不但有幅度信息，而且包含相应的相位信息，(6)式又可表示为：

Ds(ω1,ω2)=|Ds(ω1,ω2)|exp[jφD(ω1,ω2)]。

(7)

定义2 假设S(k)是确定性实数随机序列，且时长为有限值，则其傅里叶变换S(ω)存在，则有：

Ds(ω1,ω2) =S(ω1)·S(ω2)·S*(ω1+ω2)

=S(ω1)·S(ω2)·S*(-ω1-ω2)。

(8)

其中：*表示相应复数的共轭。

如果{S(k)}为一零均值高斯随机平稳序列，则对于所有的m和n，有：

Rs(m,n)=E[S(k)S(k+m)S(k+n)]=0,

(9)

所以其高斯过程的双谱必然恒定等于0。

3 算法的实现过程

3.1旋转不变性实现

Ir(ρ,θ)=I(ρ,θ-φ) 。

(10)

由傅里叶变换的性质可知，Ir(ρ,θ)对应的傅里叶变换的结果可以表示为：

Zr(r,k)=Z(r,k)exp(iφk)。

(11)

其中：Z(r,k)表示I(ρ,θ)的傅里叶变换的结果，Zr(r,k)表示Ir(ρ,θ)对应的傅里叶变换的结果，将(11)式代入(8)式，可得：

=Z(r,k1)·Z(r,k2)·Z(r,-k1-k2)

=Ds(r,ω1,ω2)。

(12)

可见，图像f(x,y)旋转前后其LPM后的双谱具有不变性。

为了降低计算的复杂度，对于LPM后的二维图像矩阵，取其n个列向量，然后对每一列向量组成的双谱绝对值矩阵进行分析，则有：

对n个ri对应的双谱绝对值矩阵分别计算其均值、方差、能量、主对角线元素之和并组合得到对应的特征向量。

3.2零数字水印注册

具体步骤如下：

步骤1：将笛卡尔坐标系的二维灰度图像通过对数极坐标变换(LPM)，得到对数极坐标系下的二维灰度图像矩阵；

步骤2：对对数极坐标系下的二维灰度图像矩阵，提取n列局部矩阵进行双谱变换，得到双谱绝对值矩阵，计算其矩阵的均值、方差、能量、主对角线元素之和并组合得到对应的特征向量(|E1||V2||S3||D4|)，将其在版权认证中心进行注册；

步骤3：对特征向量|E1||V2||S3||D4|加盖时间戳，并且连同用户的相关签名信息一并在版权保护中心进行实名注册，此时则宣布原始载体图像已在版权保护之下。

图1 零水印注册流程图

图2 零数字水印检测流程图

3.3零数字水印检测

具体步骤如下：

步骤1：首先对待检测图像分别进行对数极坐标变换和双谱变换操作，然后计算双谱绝对值矩阵的均值、方差、能量及主对角线上的和，并将最后计算的值作为数字向量B。

步骤2：计算出数字向量B与零数字水印在注册时版权认证中心的数字向量A的绝对差之和：

(13)

当d≥ε，同时时间戳记录的时间与版权保护中心提供的信息不相符时(ε为经验系数，取值为0.05)则验证结束，说明待检测二维图像中不包含版权水印信息，否则证明版权水印存在。

4 Matlab仿真结果及分析

利用Matlab 2013a软件进行仿真实验，二维宿主图像采用256×256的“Lena”灰度图像。图像的视觉质量的客观定量描述采用峰值信噪比PSNR。

(14)

PSNR=10·lg(2552/MSE)。

(15)

(a)Lena宿主灰度图像

(b)Pepper宿主灰度图像

(c)Plane宿主灰度图像

(d)Baboon宿主灰度图像

本文进行了以下相关的验证性实验：实验1是二维宿主图像遭受旋转几何变换攻击后的实验，结果如表1所示；实验2是对Lena二维宿主灰度图像进行JPEG有损压缩操作，再进行水印检测，结果如表2所示；实验3是Lena二维宿主灰度图像感染高斯噪声及滤波后进行的水印检测，结果如表3所示；最后对于其他宿主载体图像是否会出现误判(错误检测)的情况，实验4给出了零数字水印误检率结果。分别选取Plane、Pepper和Baboon 3幅二维宿主灰度图像，结果如图4所示。实验结果表明，本方案对旋转、有损压缩、平滑滤波具有一定的鲁棒性，同时通过误检率实验也可以看出本方案具有很好的检测精度。

表1 原始宿主图像遭受几何攻击后的实验结果

表2 有损压缩(JPEG压缩)的实验结果

表3 椒盐噪声、高斯噪声及平滑滤波后的分析实验结果

图4 不同载体图像的旋转实验误检曲线图

5 结语

本文提出了一种对数极坐标和双谱分析的零数字水印算法，结合Matlab仿真实验结果，说明该方法对不同旋转角度下的水印攻击有很强的稳健性，同时该方法对于感染了椒盐噪声、高斯白噪声、均值和中值滤波、JPEG有损压缩的待检测图像也具有一定的稳健性；最后误检率实验结果也显示了本方法对于没有嵌入数字水印的待检测图像其错误检测的概率也非常小，说明该算法具有很好的检测精度。因此本文所提出的零数字水印是一种较为实用的版权保护方法。

[1] 温泉,孙锬锋,王树勋.零水印的概念与应用[J].电子学报，2003,31(2):214-216.

[2] 吴伟民,丁冉,林志毅,等.基于混沌的医学图像篡改定位零水印算法[J].计算机应用研究，2014,31(12):3685-3688.

[3] 曲长波,杨晓陶,袁铎宁.小波域视觉密码零水印算法[J].中国图像图形学报，2014,19(3):367-371.

[4] 谢勇,张金龙,张雯.一种基于奇异值的抗打印/扫描的彩色图像零水印方案[J].包装学报，2014,16(3):11-13.

[5] 马建糊,何甲兴.基于小波变换的零水印算法[J].中国图像图形学报,2007,12(4):582-585.

[6] 陈伟琦,李倩.基于DWT-SVD的图像双零水印算法[J].计算机工程与应用,2014,36(10):1992-1996.

[7] 张春凯,杨德志,丛佩丽.基于Krawtchouk矩和NSCT变换的鲁棒性零水印算法[J].计算机应用与软件,2014,31(3):282-285.

[8] 冯银波,陈善学.基于四元数域的彩色图像双重零水印算法[J].计算机应用与软件,2014,31(9):264-267.

[9] 吴伟民,丁冉,林志毅,等.基于细胞自动机与奇异值分解的零水印算法[J].计算机应用，2014,34(6):1690-1693.

[10] 赵玉霞.基于混沌系统与提升小波的抗剪切攻击的彩色图像盲水印算法[J].工程图像学报,2010,24(4):216-219.

【责任编辑牛怀岗】

AZeroImageWatermarkingResistingtoRotationAttackBasedonLPMandBispectrumAnalysis

HE Bing1, 2
(1. School of Mathematics and Physics, Weinan Normal University, Weinan 714099, China; 2. Center of X-Ray Detection and Application of Shaanxi, Weinan 714099, China)

Zero image watermarking has been widely used in copyright protection, due to its robustness and better invisibility. The existing image watermarking methods cannot resist geometric transformation, for example, the original image is rotated by little angles and the watermark could not be detected. In order to improve the ability to resist the geometric rotation transformation, a zero watermarking algorithm resisting to rotation attack is proposed based on log polar mapping and bispectrum analysis. Firstly, the original image is converted to log polar coordinates from Cartesian coordinates by LPM, and then the converted image is transformed by bispectrum analysis. And calculating the numbers of simple feature about matrix of bispectrum, such as mean, variance, energy, etc. Finally, the zero watermark information is designed and constructed by using these feature vectors. The experiment results show that this proposed method can obtain better visual effect and it is robust enough to rotation attack. Meanwhile, it is also robust enough to some image degradation process such as filtering, JPEG compression and cropping attack.

LPM; Zero watermarking; robustness; bispectrum analysis

TP309.7

：A

：1009-5128(2017)16-0029-06

2017-03-19

陕西省教育厅自然科学专项计划项目：校企联合开发的第二代水印技术在版权保护中的应用(16JK1280)；渭南师范学院自然科学一般研究项目：X光数字图像知识版权防伪研究(17YKS12)

何冰(1982—)，男，陕西合阳人，渭南师范学院数理学院副教授，西安电子科技大学博士研究生，主要从事特殊函数分析、图像处理研究。