算法步骤—基于奇异值分解的小波域水印算法.docVIP

0 引言 数字水印是随着信息技术和计算机网络飞速发展而兴起的一种重要的版权保护手段，它通过在数字载体 (图像、音频、视频、文本等 )中加入隐秘信息来达到这个目的，在需要时这种隐秘信息可以通过特定方法检测出来，以验证作品是否合法或受到篡改。目前数字水印已成为多媒体信息安全研究领域的一个热点，也是信息隐藏技术研究领域的重要分支。 数字图像水印技术从实现过程上分空域算法和变换域算法。空域算法是在图像的空间域中嵌入水印的方法[1~4 ] ，空域算法通常具有较快的速度，但一般鲁棒性较差。变换域算法是指在图像的变换域中入水印的方法[5~7 ]。变换域包括DCT 域、DF T域、DWT 域等,由于图像的小波变换能够很好地匹配人类视觉系统( HVS) 的特性[8 ]，还可以与J PEG2000 标准相兼容，因此小波域数字水印技术具有很好的应用前景。 本文针对二值水印图像，结合奇异值分解与小波变换, 提出了基于奇异值分解的小波域水印算法。实验结果表明, 本文提供的方法具有较好的鲁棒性和不可见性。 1 小波变换与奇异值分解 作为一种数学工具，小波变换是对人们熟知的傅里叶变换和窗口傅里叶变换的一个重大突破，为信号分析、图像处理、量子物理及其他非线性科学的研究领域带来革命性的影响。小波变换具有许多良好的特性，这些性质奠定了小波域水印技术的基础。小波分解的空间－频率特性与HVS某些视觉特性有相似性。该特性是小波变换区别于DFT和DCT的一个重要方面，根据该特性可以将高强度的水印嵌入到 HVS 不太敏感的区域，这样在保证不影响图像视觉质量的前提下可以最大限度地增加嵌入水印的强度。 小波变换可以将图像分解成低频子带和高频子带。其中低频带表示由小波变换分解级数决定的最大尺度、最小分辨率下对原始图像的最佳逼近, 图像的大部分能量集中在此。高频带系列则分别是图像在不同尺度、不同分辨率下的细节信息，二维图像一级小波分解所得子图像按其重要性排序为LL1、HL1、LH1、HH1。 LL1 HL1 LH1 HH1 图1 二维静态图像的一级分解 奇异值分解是数值线性代数的有效算法之一，它在统计分析、信号与图像处理、系统理论与控制中被广泛应用。在图像处理中应用的优势在于图像奇异值的稳定性非常好，即当图像被施加小的扰动时，图像的奇异值不会有大的变化，并且奇异值所表现的是图像的内蕴特性而非视觉特性。 设图像矩阵为非负矩阵为非负矩阵，，其中表表示实数域，，则的奇异值分解定义如下： (1) 其中和均为酉矩阵，它们的列矢量分别为、，为对角阵，上标T表示矩阵转置，奇异置满足： (2) 非零奇异值的个数等于矩阵的秩。 2 水印算法2.1 水印嵌入算法 的值。 原始载体图像为灰度图像时水印嵌入过程如下： 步骤一：将原始的二值水印图像取反并进行Arnold变换得到，并将迭代次数作为密钥保存。 步骤二：将原始图像I进行一级小波分解，将得到中频子带HL1。 步骤三：将中频子带HL1进行SVD，得到正交矩阵、及一个对角矩阵。 步骤四：将水印叠加到矩阵上，对新产生的矩阵进行奇异值分解，得到、和。 步骤五：然后将矩阵、和相乘，得到处理后的包含水印的图像D。 步骤六：用D代替步骤二中的HL1，进行小波反变换，得到潜入水印的图像。 水印提取算法 水印提取算法恰好与水印的嵌入过程完全相反，过程如下： 步骤一：将含水印图像进行一级小波分解，可得到中频子带。 步骤二：对中频子带进行SVD，得到正交矩阵、和对角矩阵。 步骤三：用代替水印嵌入算法步骤三中的对角矩阵，然后与正交矩阵、相乘，即可得到矩阵图像D。 步骤四：利用水印嵌入算法步骤三、四中的及进行运算：，则可得到置乱后的二值水印图像。 步骤五：对进行次（表示Arnold变换的周期）Arnold变换并进行取反，即可得到最终提取的二值水印图像。