基于空域和混沌映射的数字水印嵌入算法(1)

王

伟，等：基于空域和混沌映射的数字水印嵌入算法

４期基于空域和混沌映射的数字水印嵌入算法

王

伟

朱理

毛旺农１

（湖南大学软件学院，长沙４１００８２；Ｅ－ｍａｉｌ：ｘｘｙｙｗｗｙｗｚ＠１６３．ｃｏｍ；

湖南大学北校区网络中心，长沙４１００７９）摘要提出一种ＬＳＢ改进算法，利用混沌序列对水印预处理，而后利用固定归一化相关嵌入图像数据低位，从而改进了传统

ＬＳＢ算法稳健性差的缺点。实验结果表明，采用该算法的水印在无原始图像的情况下，通过密钥能够清晰地检测到水印的存在，对常见的图像处理有较好的鲁棒性。关键词

空域

混沌鲁棒性

数字水印

图像处理

中图法分类号ＴＰ３０９；文献标识码

Ａ

２００５年１０月２５日收到

第６卷第４期２００６年２月

１６７１－１８１５（２００６）０４－０３８３－０４科学技术与工程

ＳｃｉｅｎｃｅＴｅｃｈｎｏｌｏｇｙａｎｄＥｎｇｉｎｅｅｒｉｎｇＶｏｌ．６Ｎｏ．４Ｆｅｂ．２００６

２００６Ｓｃｉ．Ｔｅｃｈ．Ｅｎｇｎｇ．

ｃ

近年来，图像数字水印技术的研究取得了很大

进展，陆续提出了诸如空域、变换域、压缩域、基于

统计学、基于生理模型等多种数字水印算法［１—３］

，有

一些水印算法［４］已经具备了较好的鲁棒性，但是这

些算法有的只能给出“有”“无”水印的回答，而无法对嵌入的所有水印比特进行准确提取。典型数字水印的空域算法是将信息嵌入到图像点中最不重要的像素位上，再通过记录提取这些信息来检测水印，该种方法具有算法简单、速度快、容易实现等特

性，已受到人们高度重视

［５，６］

。本文利用混沌随机性好的特性，借鉴了文献［５］中的一些思想，针对ＬＳＢ算法对于压缩、滤波及旋转等图像处理稳健性差的缺点，改进了传统的ＬＳＢ算法，

如图１所示。实验结果表明，该水印方案安全可靠、不可见性好、对常见的图像处理有较好的鲁棒性、易于实现。

１水印信息预处理

混沌区的数据有两个特性：迭代不重复性和初值敏感性。当选定适当系数使方程进入混沌状态时，方程将进行无限不循环迭代，因此不会出现重复的迭代值。任何人如果不知道其初值都无法预测下一个迭代值，这种迭代的结果可以用来产生随机序列。

考虑一维离散混沌映射ｆ：Ｕ→Ｕ，Ｕ"Ｒ，

ｘｎ＋１＝ｆ（ｘｎ，λ），ｘｎ∈Ｕ，λ∈Ｒ

（１）

其中ｎ＝１，２，

…表示迭代次数，λ是控制系统混沌行为的参数。先由二值水印图像获得水印信号Ｗ，然后利用行扫描将其映射为一维向量Ｗ１，Ｗ１＝｛Ｗ１（

ｉ）｜Ｗ１（ｉ）∈｛０，１｝，０≤ｉ＜ｍ×ｎ｝，ｍ×ｎ为水印图像大小。由式（１）生成的实数混沌序列二值化后可得｛０，１｝二值混沌序列，然后利用密钥Ｋ１可得十进制混沌伪随机序列Ｐ，Ｐ＝｛Ｐ（ｉ）｜Ｐ（ｉ）∈｛０，１｝，０≤ｉ＜ｍ×ｎ｝

。最后用生成的序列Ｐ对Ｗ１进行混沌调制，得到待嵌入水印信号Ｗｔ，Ｗｔ（ｉ）＝Ｗ１（ｉ）’ｂ（ｉ）（０≤ｉ＜ｍ×ｎ），’在这里表示异或运算。

２水印嵌入

度量鲁棒性最简单的方法就是假设具有高检

图１水印嵌入和提取框图

科学技术与工程６卷

测值的水印比低检测值的水印鲁棒性更好，这样，检测统计量本身就成了鲁棒性的度量尺度。这里的解决方案是用一个固定归一化相关算法来进行调整，而不是调整线性相关。对每一幅作品，通过直方图处理，将每个像素值由８比特压缩为７比特，将所有像素值的第８个比特作为嵌入区域，在嵌入区域内找到与水印归一化相关期望值最接近的点。嵌入区域是一个以参考向量Ｗｒ为中心的锥体表面部分，因此，只需要找出锥体上与给定点ｖ０距离最近的点就可以了。嵌入的操作说明如图２。

很明显这些点位于同时包含ｗｒ和ｖ０的平面内，通

过Ｇｒａｍ－Ｓｃｈｍｉｄｔ规格化正交，为这个平面寻找两个相互垂直的坐标轴，从而把这个问题转换成二维问题。

Ｘ＝Ｗｒ

Ｗｒ

（２）Ｙ＝ｖ０－Ｘ

（ｖ０Ｘ）ｖ０－Ｘ

（ｖ０Ｘ）（３）平面上每个点都能用坐标值ｘ和ｙ表示为ｘＸ＋ｙＹ，

Ｘ轴与参考时间向量重合。ｖ０的坐标是

（ｙｖ０

必须是正值）：

ｘｖ０

＝ｖ０Ｘ

（４）ｙｖ０

＝ｖ０Ｙ

（５）

因为ｙｖ０

为正值，所以最近的点应该在上半平面

的线上，这条线用向量（ｘｔ，ｙｔ）表示，其中ｘｔ＝ｃｏｓθ，ｙｔ＝ｓｉｎθ，而角θ是嵌入区域的锥形内角的半角，由于θ＝

ａｒｃｃｏｓ（ｔｎｃ），ｔｎｃ是期望的归一化相关性，

于是有：ｘｔ＝ｔｎｃ

（６）ｙｔ＝１－ｔ２ｎｃ

!（７）

算法步骤如下：

Ｓｔｅｐ１：用经过散列处理的混沌序列

（其混沌初值作为密钥Ｋ１）置乱水印Ｗ；Ｓｔｅｐ２：给定一个混沌初值（密钥Ｋ２）生成混沌序列｛Ｐｒ｜ｒ＝１，２，…｝其个数为水印大小；Ｓｔｅｐ３：在通过ｘｔ，ｙｔ的直线上，找到离ｖ０最近的点ｖｗ

（在Ｗｒ的检测区域内）：ｘｖｗ＝ｘｔ（ｘｔｘｖｗ＋ｙｔｙｖ０）（８）ｙｖｗ＝ｙｔ（ｘｔｘｖｗ＋ｙｔｙｖ０）

（９）Ｓｔｅｐ４：取Ｐｒ对计算出的新向量ｖｗ嵌入水印，得到水印图像Ｉｗ。

３水印提取

水印提取过程中所用密钥Ｋ１和Ｋ２与嵌入过程相同，具体步骤如下：Ｓｔｅｐ１：给定密钥Ｋ２在图像Ｉｗ中找到对应的向量ｖｗ，

得到水印信息Ｗｔ；Ｓｔｅｐ２：再由密钥Ｋ１生成二值混沌序列Ｐ对Ｗ１解调，即：

Ｗ１′（ｔ）＝Ｗｔ（ｔ）

"Ｐ（ｔ），ｔ＝０，１，２，…，ｍ×ｎ－１；Ｓｔｅｐ３：最后按行扫描顺序将Ｗ１′

映射成二维矩阵Ｗ′，Ｗ′＝｛Ｗ′（ｉ，ｊ），０≤ｉ＜ｍ，０≤ｊ＜ｎ｝，Ｗ′（ｉ，ｊ）∈｛０，１｝，Ｗ′即是恢复出的水印信号；Ｓｔｅｐ４：采用如下的标准衡量二值水印的失真率（其中Ｗ表示Ｗ的逻辑非运算）：

ｗｓ＝

∑ｉ＝１ｉ＝ｍ∑ｊ＝１ｊ＝ｎ

Ｗ＇（ｉ，ｊ）"Ｗ（ｉ，ｊ）ｍ×ｎ

４实验结果

本文用大小为２５６×２５６

的２５６灰度Ｌｅｎａ图像用

ＭＡＴＬＡＢ进行水印嵌入和提取的仿真实验，采用大小为６４×６４视觉可辨认的二值图像作为数字水印，嵌入水印后的图像如图３所示，结果表明该方法有很好的不可见性。

图２固定归一化的嵌入方法

３８４