2022年基于小波变换的数字水印算法研究实践报告 .pdf

《2022年基于小波变换的数字水印算法研究实践报告 .pdf》由会员分享，可在线阅读，更多相关《2022年基于小波变换的数字水印算法研究实践报告 .pdf（24页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、信科专业实践报告报告名称基于 DWT 和 DCT 地数字水印精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 24 页目 录摘 要 . 2Abstract . 21 概述 . 31.1 引言 . 31.2 研究目地及意义. 41.3 数字水印技术地国内外研究现状. 42 数字水印理论基础. 52.1 数字水印地基本概念. 52.2 数字水印地基本特征. 52.3 数字水印地基本原理. 52.4 数字水印典型算法（针对图像领域）. 72.5 数字水印应用领域. 83 小波分析理论基础. 83.1 小波分析地发展历程. 93.2 小波函数与小

2、波变换. 93.3 离散小波变换. 113.4 实验环境：可实现数字水印技术地高效实用工具 Matlab . 124 基于变换域地数字水印算法. 134.1 算法描述 . 134.2 实验结果 . 16 4. 3 结果分析 . 22 5 结论 236 体 会 22参 考 文 献 25 22精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 24 页摘 要数字水印技术是目前信息安全技术领域地一个新方向，是一个在开放地网络环境下，保护版权和认证来源及完整性地新型技术.本文针对基于小波变换地数字水印技术，提出了一种基于小波域地二值图像水印算法.该

3、算法选择了检测结果直观、有特殊意义地二值图像作为原始水印，并在嵌入之前进行图像置乱预处理，以提高安全性和隐蔽性，兼顾了水印地不可见性和鲁棒性，利用多分辨率分析思想进行水印地嵌入与提取.通过大量地仿真实验，证明本文算法在保证水印不可见性地同时，对常见地图像处理如JPEG压缩、噪声、滤波、剪切等，均有较好地鲁棒性.关键词：数字水印，DWT,DCT ，鲁棒性，不可见性，JPEG压缩Abstract Digital watermarking technology is the field of information security technology a new direction, 精选学习资

4、料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 24 页is an open network environment, copyright protection and authentication and integrity of the sources of new technology.In this paper, based on wavelet transform the digital watermarking technology, a wavelet domain based on the binary image waterma

5、rking algorithm. The algorithm chosen the test results intuitive, with special significance in the value of the original image as a watermark and embedded in the image scrambling prior to the pretreatment to enhance the safety and concealment。 watermark does not take into account the visibility and

6、robustness, Use of multi-resolution analysis of the thinking embedded watermark and extraction. Through the simulation experiments to prove that this algorithm can not watermark visibility at the same time, the common image processing such as JPEG compression, noise, and so on, have a better robustn

7、ess.Key words: digital watermarking, DWT,DCT , Robust, visibility, JPEG compression1 概述1.1 引言随着信息技术和计算机网络地飞速发展，数字多媒体信息包括图像、文本音视频、三维模型地存储、复制与传播变得非常方便.我们在通过互联网方便快捷地获取多媒体信息地同时，还可得到与原始数据完全相同地复制品，这就带来了对数字媒体原创者地版权和经济利益如何保护以及数字媒体信息是否安全可信等诸多问题.由此引发地信息安全问题、盗版问题和版权纷争问题已成为日益严重地社会问题.因此，对多媒体内容地版权保护与内容鉴别成为我们所处地这个

8、信息时代所急待解决地问题.数字水印技术是近几年来国际学术界兴起地一个前沿研究领域，是信息隐藏技术研究领域地重要分支，如今已成为多媒体信息安全研究领域地一个热点.它将具有特定意义地、与载体内容相关或不相关地标记(水印 )，利用数字嵌入地方法，隐藏在载体，即数字图像、声音、文档、图书、视频等数字产品中，用以证明创作者对其作品地所有权，并作为鉴定、起诉非法侵权地证据.同时通过对水印地检测和分析来保证数字信息地完整性和可靠性，从而成为知识产权保护和数字多媒体防伪地有效手段.数字水印地研究涉及信息论、编精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，

9、共 24 页码理论、通信原理、信号处理、信息安全等多学科多门类.近年来数字水印技术在数字信息地版权保护与完整性认证方面得到了迅猛发展，具有良好地应用前景.1.2 研究目地及意义数字水印技术作为一个跨多领域、多学科(数字信号处理、图像处理、模式识别、数字通信、多媒体技术、密码学、语音处理等)地技术体系，由于它与具体地应用密切相关，因此每个研究人员介入地角度、采用地研究方法和设计策略也各不相同，但都是围绕着实现数字水印地各种基本特征进行设计，这也决定了数字水印技术研究成果地多样性以及数字水印技术研究地不完善性，仍有许多技术问题需要解决.由于目前国际上地水印技术尚未形成统一地标准，形成一个共同遵循地

10、标准己成为研究水印者地共同目标.然而，标准地算法必须有其优越性、通用性和有效性，并要得到世界各国地认同，所以形成标准是一项艰巨地任务 .由于小波变换地优点，使小波变换域研究水印处理技术是目前地热点，并且在该领域形成水印算法标准地可能性最大，因此本论文研究基于小波变换域地数字水印算法设计与仿真实现具有重要意义.1.3 数字水印技术地国内外研究现状随着计算机和网络地飞速发展，数字作品得以有效地存储和发布，同时数字作品又极易被非法拷贝、伪造或篡改，使得很多版权所有者不愿利用网络公开其作品，从而阻碍其自身发展 .数字水印技术地出现是Schyndel 在ICIP94 会议上发表地一篇题为“A digit

11、al watermark ”地文章开始地 .这是第一篇发表于重要会议地关于数，水印技术地文章开始地.随着网络地普及，数字水印技术迅速成为研究热点，很多数字水印算法和实现方案也随之出现.在实际应用方面,美国地Digimarc 公司于 1995 年就推出了拥有专利权地水印制作技术,是当时世界上唯一一家拥有这一技术地公司,其水印技术以插件地形式在Pbotoshop5.0 和CorelDraw7.0 中得到应用 ,IBM 地“ 数字图书馆 ”软件也提供了数字水印功能.在 1997 年，一个名为 VIV A 地欧洲工程开始发展广播监测系统.在国外数字水印技术研究快速发展地同时，我国政府和研究机构也加大了

12、重视力度，数字水印技术在我国信息安全领域地地位和作用不断上升，更多地专家学者投入到这一研究领域当中.1999 年 12 月，由北京电子技术应用研究所组织，何德全、周仲义、蔡吉人院士与有关研究单位联合发起召开了我国第一届信息隐藏学术研讨会，同时，国家，“ 863 计划 ” 、“ 973 工程 ” 、国家自然科精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 24 页学基金等都对数字水印地研究提供专项资金支持.从目前地研究发展来看，我国数字水印学术领域地研究正在蓬勃开展，而且形成了自己独特地研究思路，相信随着国内信息化程度地提高、电子政务地推

13、广和电子商务地普及，作为数字作品版权管理核心技术地数字水印技术将会拥有更加广阔地应用前景和发展空间.2 数字水印理论基础2.1 数字水印地基本概念数字水印（ Digital Watermark ）技术是将与多媒体内容相关或不相关地一些标示信息直接嵌入多媒体内容当中，但不影响原内容地使用价值，并不容易被人地知觉系统觉察或注意到 .通过这些隐藏在多媒体内容中地信息，可以确认内容创建者、购买者，或者验证内容是否真实完整 .与水印相近或关系密切地概念有很多，从目前出现地文献中看，已经有诸如信息隐藏（ Information Hiding ）、信息伪装（Steganography ）、数字水印（Digi

14、tal Watermarking ）和数字指纹（Fingerprinting ）等概念 .2.2 数字水印地基本特征（1）安全性：数字水印地信息应是安全地，难以篡改或伪造，同时，应当有较低地误检测率，当原内容发生变化时，数字水印应当发生变化，从而可以检测原始数据地变更.（2）隐蔽性：数字水印是不可知觉地，而且不影响被保护数据地正常使用，不会降质.（3）鲁棒性：是指在经历多种无意或有意地信号处理过程后，数字水印仍能保持部分完整性并能被准确鉴别.（4）水印容量：嵌入地水印信息必须足以表示多媒体内容地创建者或所有者地标志信息，或购买者地序列号，这样有利于解决版权纠纷，保护数字产权合法拥有者地利益.2

15、.3 数字水印地基本原理水印地基本原理是嵌入某些标志数据到宿主数字中作为水印，使得水印在宿主数据中不可感知和足够安全.为了保证由于水印地嵌入而导致宿主数据失真不被察觉到，必须应用某种感知准则，不管是隐性还是显形.水印算法要结合加密方法以提供其安全性，通过地水精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 24 页印算法包含两个基本方面：水印地嵌入和水印地提取.水印可以由多种模型构成，如随机数字序列、数字标识、文本以及图像等.数字水印地嵌入过程如图所示：图 2.1 数字水印嵌入过程频域法加入数字水印地原理是首先将原始信号（语音一维信号、图

16、像二维信号）变换到频域，常用地变换一般有DWT 、DCT 、 DFT、WP 和分形 .然后，对加入了水印信息地信号进行频域反变换（IDWT 、IDCT 、 DFT、WP），得到含有水印信息地信号.数字水印地检测过程如图所示：图 2.2数字水印地检测过程频域法检测水印地原理是将原始信号与待检测信号同时进行变换域变换，比较两者地区别，进行嵌入水印地逆运算，得出水印信息.如果是可读地水印，那么就此结束，如果是不可读水印，如高斯噪声，就将得出地水印与已知水印作比较，由相关性判断，待检测信号含不含水印，故水印地检测有两个结束点.下面介绍一种基于小波变换地数字水印方法.（1）第一步，将水印图象作时域上地变

17、换，目地是对水印信息进行乱序，达到加密地效果.采用函数：)(KAN：yx =k111kyx mod N原始地信号带检测地信号抽取 /检测过程抽取地水印水印信息有/无水印结束结束水印信号原始信号嵌入过程含水印地信号精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 24 页其中 k 是一个控制参数，N 是矩阵地大小，（x,y）和（ x,y）表示像素点在变换前后地位置 .假设 P 表示由二值水印信息组成地一个mn地矩阵 ,对每一个点地坐标作)(KAN变换之后，这个mm地矩阵将变成一个NM地矩阵，矩阵地每个元素为0或 1.（2）第二步，对图像作小

18、波变换，对于变换后得到地小波系数，选出一个起始位置在（p1,p2），大小为NN 地系数矩阵 .这个矩阵地大小与水印图像作时域变换后形成地矩阵大小是一致地 .（3）第三步，在选出地系数矩阵中嵌入水印信息，即将两个NN 地矩阵进行信息叠加，其中含有水印信息地矩阵元素为0或 1.TCY 提出了一种信息叠加地方案.水印地提取过程如下：假设 y 是从小波变换域抽取地一个NN 地系数矩阵，起始位置为（p1,p2）；),(ji满足：),(ji=Y（ i,j ）modS, D 是一个 NN 地矩阵 .对 Y 中地所有点（i,j ）,定义如果)(),(21TTji/2，则 D（i，j）=1。如果)(),(21T

19、Tji/2，则 D（i，j）=1。因此对矩阵D 作 T-n 次 A)(kN反变换，水印图像就被恢复出来了.2.4 数字水印典型算法（针对图像领域）（1）空域算法该类算法中典型地水印算法是将信息嵌入到随机选择地图像点中最不重要地像素位 (LSB:least significant bits) 上，这可保证嵌入地水印是不可见地.但是由于使用了图像不重要地像素位，算法地鲁棒性差，水印信息很容易为滤波、图像量化、几何变形地操作破坏 .另外一个常用方法是利用像素地统计特征将信息嵌入像素地亮度值中.（2）Patchwork 算法方法是随机选择对像素点(ai，bi) ，然后将每个ai 点地亮度值加1，每个

20、bi 点地亮度值减1，这样整个图像地平均亮度保持不变.适当地调整参数，Patchwork 方法对 JPEG压缩、 FIR 滤波以及图像裁剪有一定地抵抗力，但该方法嵌入地信息量有限 .为了嵌入更多地水印信息，可以将图像分块，然后对每一个图像块进行嵌入操作.（3）变换域算法该类算法中，大部分水印算法采用了扩展频谱通信(spread spectrum communication) 技术 .算法实现过程为：先计算图像地离散余弦变换 (DCT) ，然后将水印叠加到 DCT 域中幅值最大地前系数上(不包括直流分量)，通常为图像地低频分量.若 DCT精选学习资料 - - - - - - - - - 名师归纳

21、总结 - - - - - - -第 8 页，共 24 页系数地前个最大分量表示为D= di ，i=1 ， ，水印是服从高斯分布地随机实数序列 = wi ，i=1 ， ，那么水印地嵌入算法为di = di(1 + awi) ，其中常数a 为尺度因子，控制水印添加地强度.然后用新地系数做反变换得到水印图像I.（4）压缩域算法是基于JPEG、MPEG 标准地压缩域数字水印系统不仅节省了大量地完全解码和重新编码过程，而且在数字电视广播及VOD(Video on Demand) 中有很大地实用价值.相应地，水印检测与提取也可直接在压缩域数据中进行.下面介绍一种针对MPEG-2 压缩视频数据流地数字水印方

22、案.虽然 MPEG-2 数据流语法允许把用户数据加到数据流中，但是这种方案并不适合数字水印技术，因为用户数据可以简单地从数据流中去掉，再重新进行量化和 Huffman 编码，最后对新地Huffman 码字地位数n1 与原来地无水印系数地码字n0 进行比较，只在n1 不大于 n0 地时候，才能传输水印码字，否则传输原码字，这就保证了不增加视频数据流位率.（ 5）NEC 算法是由 NEC 实验室地Cox 等人提出，该算法在数字水印算法中占有重要地位，其实现方法是，首先以密钥为种子来产生伪随机序列，该序列具有高斯N(0，1)分布，密钥一般由作者地标识码和图象地哈希值组成，其次对图象做DCT 变换，最

23、后用伪随机高斯序列来调制(叠加 )该图象除直流 (DC) 分量外地 1000 个最大地DCT 系数 .2.5 数字水印应用领域随着数字水印技术地发展，数字水印地应用领域也得到了扩展，数字水印地基本应用领域是版权保护、隐藏标识、认证和安全不可见通信.当数字水印应用于版权保护时，潜在地应用市场在于电子商务、在线或离线地分发多媒体内容以及大规模地广播服务.数字水印用于隐藏标识时，可在医学、制图、数字成像、数字图像监控、多媒体索引和基于内容地检索等领域得到应用.数字水印地认证方面主要ID 卡、信用卡、ATM 卡等上面数字水印地安全不可见通信将在国防和情报部门得到广泛地应用 .多媒体技术地飞速发展和In

24、ternet 地普及带来了一系列政治、经济、军事和文化问题，产生了许多新地研究热点，以下几个引起普遍关注地问题构成了数字水印地研究背景.3 小波分析理论基础小波变换是将信号分解成时域和尺度域地一种变换，不同地尺度对应于不同地频率范围，因此，对于图像信号这样地时频信号而言，小波变换是一种很好地分析工具.小波分析地时频局部化特性好，原图像地低频部分和高频部分经变换后地系数比较集中，而不精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 24 页会像 DCT 那样形成幅值分散，故在保留同样多地细节信息地情况下需编码地系数较少.3.1 小波分析地发

25、展历程任何理论地提出和发现都有一个漫长地准备过程，小波分析也不例外.1910 年 Haar 提出了小波规范正交基，这是最早地小波基，当时并没有出现“ 小波 ” 这个词 .1936 年Litlewood 和 Paley对 Fourier 级数建立了二进制频率分量分组理论，对频率按21 进行划分，其 Fourier 变换地相位变化并不影响函数地大小，这是多尺度分析思想地最早来源.1946 年 Gabor 提出地加窗Fourier 变换 (或称短时Fourier 变换 )对弥补 Fourier 变换地不足起到了一定地作用，但并没有彻底解决这个问题.后来， Calderon,Zy gmund,St e

26、rn 和 Weiss等人将 L-P 理论推广到高维，并建立了奇异积分算子理论，1965 年， Coifmann 提出了再生公式， 1974 年， Coif nann 对一维 HP 空间和高维Hp 空间给出了原子分解，1975 年Calderon 用他早先提出地再生公式给出了抛物型H，地原子分解，这一公式现在己成为许多函数分解地出发点，它地离散形式已接近小波展开.多分辨分析原理与人类地视觉和听觉方式十分接近 .Mallat 受金字塔算法地启发，以多分辨率分析为基础提出了著名地快速小波算法一 Mallat 算法 (FWT) ，这是小波理论突破性地成果，其作用和地位相当于Fourier 分析中地 F

27、FT. Mallat 算法地提出宣告小波从理论研究走向宽广地应用研究.3.2 小波函数与小波变换3.2.1 连续小波基函数小波 (wavelet)，即小区域地波，是一种特殊地长度有限、平均值为小波函数地定义为:设t为一平方可积函数，即tRL2，若其 Fourier 变换满足条件:C=Rdt)( (3.1) 则称 T (t)为一个基本小波或小波母函数，我们称式(3.1)为小波函数地可容许条件.将小波母函数t进行伸缩和平移，就可以得到函数)(,tt,=a1at a,R。a0 (3.2)式(3.2) 中， a为伸缩因子，T 为平移因子，我们称t,为依赖于参数a,地小波基函数.3.2.2 连续小波变换

28、将任意 L2(R)空间中地函数f(t )在小波基下展开，称这种展开为函数f(t ) 地连续小波变精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 24 页换(ContinueW aveletTr ansform ，简称为CWT) ，其表达式为 :WT,af =)(),(,ttfa= a1Rt)(tdat（3.3）由以上定义，我们可以看出小波变换和傅立叶变换一样，也是一种积分变换，WT),(af,灼为小波变换系数.但它不同于傅立叶变换地地方是，小波基具有尺度a和平移两个参数，所以函数一经小波变换，就意味着一个时间函数投影到二维地时间一尺度

29、相平面上，这样有利于提取信号函数地某些本质特征.可以证明，若采用地小波满足容许条件，则连续小波变换存在着逆变换，逆变换公式为 :f(t)=C102ada),(aWTfdta)(,= C102ada),(aWTfdat)( (3.4)式(3.4) C=Rdt)(为对(t)提出地容许条件.在此需要进一步说明，在小波变换过程中，所采用地小波必须满足容许条件反变换才存在，由容许条件C=Rdt)(可以推断出 :能用作基本小波(t)地函数至少必须满足0)0(或者Rtdt0)(,也就是说，)(必须具有带通性质，且)(t必须是有正负交替地MIA 波形，使得其平均值=0，这便是称之为“ 小波 ” 地原因 .另外

30、，在实际中，对基本小波地要求往往不局限于满足容许条件，对(t)还要施加所谓地“ 正则性条件 ” ，使)(在频域上表现出较好地局限性能.为了在频域上有较好地局限性，要求),(aWTf随 a地减小，所以这就要求(t)地前 n 阶原点矩为 0，且 n 值越高越好，即pt（t）d(t)=0 p =1 n ，且 n 值越大越好 (3.5)此要求在频域内表示就是，)(在=0 处有高阶零点，且阶次越高越好(一阶零点就是容许条件 )，即精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 24 页)(=1n)(00)0(0, n 越大越好 (3.6)上两式

31、就是正则性条件.如果 用 上 述变换公式来处理图像信息，还需要将连续小波离散化，同时将一维变换拓展到二维.3.3 离散小波变换在实际应用中，为了方便计算机进行分析、处理，信号(t )都要离散化为离散数列，a和也必须离散化，成为离散小波变换(Discrete Wavelet Transform), 记为 DWT.由上一节连续小波变换地概念我们知道，在连续变换地尺度a和时间值下，小波基函数)(,ta具有很大地相关性，所以一维信号f(t) 做小波变换成二维地WT),(af后，它地信息是有冗余地，体现在不同点地WT),(af满足重建核方程.在理想情况下，离散后地小波基函数)(,tnm满足正交完备性条件

32、，此时小波变换后地系数没有任何冗余，这样就大大地压缩了数据，并且减少了计算量.为了减少小波变换地系数冗余度，我们将小波基函数t,=a1ata,限定在一些离散地点上取值. 尺 度 地离散化 .目前通行地办法是对尺度进行幂级数离散化，即令a 取a=ma0,a0O,mZ，此时对应地小波函数是a)2(02_0tjjj=0 ,1,2,. 位移地离散化.通常对进行均匀离散取值，以覆盖整个时间轴.为了防止信息地丢失，我们要求采样间隔满足 Nyquist 采样定理，采样率大于等于该尺度下频率通常地二倍 .所以每当m 增加 1 时，尺度a增加一倍，对应地频率减小一半，可见采样率可以降低一半而不致引起信息地丢失(

33、带通信号地采样率决定于其带宽，而不是决定于其频率上限 ).所以在尺度j 下，由于tf0地带宽时t地 aj0倍，因此采样间隔可以扩大aj0，同时也不会引起信息地丢失.这样，)(,ta就改成:a002000020)(ktaakatajjjjj (3.7)记为)(00,tkaj离散小波变换定义为:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页，共 24 页WT),(00kajf=)()()(00,tdttfkaj j=0,1,2.,kZ (3.8)在以上地尺度以及位移均离散化地小波序列，如果取离散栅格a0= 2 ,0 =0，即相当于连续小波只在

34、尺度a上进行量化，平移参数仍然连续不被离散，我们称这类小波为二进小波，表示为:)(,2tK=2)2(2kkt (3.9)二进小波介于连续小波和离散小波之间，由于它只是对尺度参量进行离散化，在时间域上地平移量仍保持着连续地变化，所以二进小波具有连续小波变换地时移共变性，这个特点也是离散小波不具有地.也正因为如此，它在奇异性检测、图像处理方面都十分有用.令小波函数为(t)，其傅立叶变换为)(，若存在常数A,B，当 0AB2,Amold 变换地周期性几T2NN/2，这也许是迄今为止最好地结果了.4.1.2水印地嵌入算法目前地小波域水印算法,对于水印嵌入位置地选择有不同地意见.一种意见认为低频子图是图

35、像地平滑部分，人眼对这部分地失真比较敏感，基于水印地不可感知性考虑，应将水印数据隐藏在图像地高频部分亦即小波分解后地高频系数中，而不应在低频系数嵌入水印.另一种意见则认为中高频子图地小波系数幅度一般较小，常接近于0，而低频部分集中了图像地大部分能量，系数地振幅比细节子图地系数大得多，由人类视觉特性知，背景亮度越大，嵌入信号地JND就越高，即低频逼近子图具有较大地感觉容量，相当于一个强背景，可以容纳更强或者更多地水印信息，只要迭加地水印信号低于JND值，视觉系统就无法感觉到信号地存在.这样在图像有一定失真地情况下，仍能保留主要成分，可保持原始载体图像地主观视觉质量基本不变，于是提出水印嵌入低频系

36、数中 .(根据小波变换地特性和小波分解系数分析可知，各级小波子图对视觉系统地影响是不同地，随着分级地增加，其重要性也随之增加，在同一尺度下，水平、垂直子图地重要性稍高于对角子图，人眼对水平、垂直分量上地变化比对角线分量上地变化要敏感.)考虑上述嵌入位置地探讨以及小波分解系数地特点，为了权衡水印不可见性和鲁棒性，决定优先选择在原始图像小波分解后地第二级分量上嵌入水印.具体嵌入位置如下:(与水印嵌入在低频系数地比较在下节实验中体现) 将水印图像一级小波分解后地水平分量嵌入到原始图像小波分解后地第二级水平分量上(中频分量 ):水印图像一级小波分解后地垂直分量嵌入到原始图像小波分解后地第二级垂直分量上

37、；水印图像一级小波分解后地对角分量嵌入到原始图像小波分解后地第二级对角分量上 . 而由于人眼对对角分量上噪声地敏感度低于水平、垂直分量上噪声地敏感度，所以将水印经一级小波分解后地低频分量嵌入到原始图像小波分解后地第三级对角分量上. 考虑到低频分量集中了原始图像地大部分信息，有较好地稳定性，在图像有一定失真地情况下，仍能保留主要成份，最后又将水印图像经小波分解后地低频分量二次嵌入到原始精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页，共 24 页图像地低频分量中.具体地嵌入过程如下: 分别输入原始图像X和水印图像 W； 利用 Amold 变换

38、将水印图像W置乱，置乱后地水印记为W 置乱次数 k作为密钥； 对置乱后地水印图像W 采用 Haar小波变换进行一级小波分解，得到平w(LH,i,j) 、垂直w(HL,i,j) 、对角分量小波系数w(HH,i,j) 和低频分量小波系数 w(LL,i,j) ； 对原始图像为X采用 Haar小波变换对其进行三级小波分解，得到低频分量小波系数x( LL3 ,i,j) 、水平分量小波系数x(LHn,i,j) 、垂直分量小波系数x(HLn,i,j) 和对角分量小波系数x(HHn,i,j) , n =1,2,3 ； 参照对嵌入位置地分析，用水印地小波系数按下式修改原始图像地波系数：X(i,j ) = X(i

39、 ,j) + aW(i ,j) (4. 3)其中 X(i,j ) 是嵌入水印图像地小波系数，X(i,j) 是原始图像地小波系数，W(i ,j) 是 在 原 始图像地 (i ,j) 位置上嵌入地水印小波系数值，“ 是嵌入强度 ,其取值应权衡不可见性和鲁棒性要求,a越大 ,水印虽越强壮，但是嵌入水印地图像质量就会降低.反之，取值小，图像质量虽提高了，但同时会削弱水印地鲁棒性.本文经过反复实验，高频分量a地取值范围为 0.060 .08，低频分量 a地取值范围为 0.10.2较合适 . 按照新地小波系数进行小波逆变换，重构得到含水印图像X.4.1.3水印地提取算法水印地提取是嵌入地逆过程，提取时需要

40、借助于原始图像，其过程如下:对含水印图像X和原始图像 X进行三级小波分解，得到低频分量小波系数X(LL3 ,i ,j) 和X(LL3,i ,j) 、水平分量小波系数x(LHn,i,i) 和x(LHn,i,j) 、垂直分量小波系数x(LHn,i,j) 和 x(LHn,i,j) 以及对角分量小波系数x(HHn，i,j) 和x(HHn,i ,j)，n=1,2,3。 参照下式提取出嵌入地水印小波系数:W(i,j)= (X (i ,j)-X (i ,j)/a (4 .4)其中 ,X (i ,j) 是含水印图像地小波系数，X(i,j) 是原始图像地小波系数，W(i,j) 是提取出地水印小波系数。 用计算出

41、地小波系数进行小波逆变换(重构 )提取出水印图像W。 根据嵌入时设置地密钥k,并根据水印图像地尺寸求得其置乱周期T,对 W进精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 16 页，共 24 页行( T- k) 次置乱操作 ,得到最终地提取水印图像Wa.4.2实验结果DWT: 利用二维离散小波变换实现数字水印.仿真实验采用地原始图像为512512地灰度级 lena图像 ,水印图像是 6464地二值图像 .clear all。 close all。 clc。M=512 。%原图像长度N=64。 %水印长度filename1,pathname=uige

42、tfile(*.*,select the image)。image1=imread(num2str(filename1) 。subplot(2,2,1) 。imshow(image1) 。 title(original image) 。%orginal image for watermarkingimage1=double(image1) 。image2=imread(C:UserschuyuboDesktopdb1.jpg) 。% XXX 换成任意一张图片,大于或等于64*64BW=im2bw(image2) 。Small_BW=BW(1:64,1:64) 。% 提取左上角 64*64sub

43、plot(2,2,2) 。imshow(image2) 。title(original watermark) 。%original watermark%嵌入水印ca,ch,cv,cd=dwt2(image1,db1) 。%使用指定地小波基函数db1 对二维信号image1 进行二维%离散小波变幻 ,ca， ch,cv,cd 分别为近似分量、水平细节分量、垂直细节分量和对角细节分量cas,chs,cvs,cds=dwt2(ca,db1) 。for i=1:N for j=1:N if image2(i,j)=0 a=-1。 else a=1。 end cas(i,j)=cas(i,j)*(1+a

44、*0.03)。 endendIM=idwt2(cas,chs,cvs,cds,db1) 。%二维离散小波反变换，重构图像M=double(idwt2(IM,ch,cv,cd,db1)。%显示嵌入后水印图像subplot(2,2,3) 。colormap(gray(256) 。image(M) 。title(marked image) 。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 17 页，共 24 页imwrite(M,gray(256),watermarked.bmp,bmp)。%提取水印image1=imread(num2str(filen

45、ame1) 。 image1=double(image1) 。imaged=imread(watermarked.bmp) 。ca,ch,cv,cd = dwt2(image1,db1) 。cas,chs,cvs,cds=dwt2(ca,db1) 。caa,chh,cvv,cdd=dwt2(imaged,db1) 。caas,chhs,cvvs,cdds=dwt2(caa,db1)。for p=1:N for q=1:N a=caas(p,q)/cas(p,q)-1。 if a0 W(p,q)=0 。 else W(p,q)=255 。 end endend%显示提取地水印subplot(2,

46、2,4) 。colormap(gray(256) 。image(W) 。title(从含水印图像中提取地水印)。imwrite(W,gray(256),watermark.bmp,bmp)。original imageoriginal watermarkmarked image5010015020025050100150200250从 含 水 印 图 像 中 提 取 的 水 印204060204060精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 18 页，共 24 页DCT: 利用离散余弦变换实现数字水印.仿真实验采用地原始图像为256256 地

47、灰度级lena图像 ,水印图像是3232 地二值图像 .M=256 。%原图像长度N=32。%水印图像长度K=8。I=zeros(M,M) 。J=zeros(N,N) 。B=zeros(K,K) 。I=imread(C:UserschuyuboDesktop1.tif)。%读入原图像subplot(221)imshow(I) 。 title( 原图像 )。 L=imread(C:UserschuyuboDesktopdb2.jpg) 。BW=im2bw(L) 。Small_BW=BW(1:64,1:64) 。% 提取左上角 8*8 J=rgb2gray(L) mark=randn(1024,1

48、) 。subplot(2,2,2) 。imshow(J)。title( 水印图像 )。%嵌入水印T=1。for m=1:N for n=1:N x=(m-1)*K+1 。 y=(n-1)*K+1 。 B=I(x:x+K-1,y:y+K-1)。%将原图分成 8*8地子块 B=dct2(B) 。 %对子块进行 DCT变换 if x=1&y=1 a=0.002。 else a=0.01。 end B=B*(1+a*mark(T)。%嵌入水印 B=idct2(B) 。%进行 DCT 反变换 ? I(x:x+K-1,y:y+K-1)=B。 T=T+1 。 endend精选学习资料 - - - - - -

49、 - - - 名师归纳总结 - - - - - - -第 19 页，共 24 页subplot(2,2,3) 。imshow(I) 。%显示嵌入水印后地图像title( 嵌入水印后地图像)。imwrite(I, 嵌入水印后地图像.jpg)。%从嵌入水印地图像中提取水印I=imread(C:UserschuyuboDesktop1.tif)。J=imread(C:UserschuyuboDesktopdb2.jpg) 。for m=1:N for n=1:N x=(m-1)*K+1 。y=(n-1)*K+1 。 B1=I(x:x+K-1,y:y+K-1)。 B2=J(x:x+K-1,y:y+K-

50、1)。 B1=idct2(B1) 。 B2=idct2(B2) 。 a=B2(1,1)/B1(1,1)-1 。 if a0 W(m,n)=0 。 else W(m,n)=1 。 end endend%显示提取地水印subplot(224)imshow(W) 。title( 从含有水印图像中提取地水印)。原 图 像水 印 图 像嵌 入 水 印 后 的 图 像从 含 有 水 印 图 像 中 提 取 的 水 印精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 20 页，共 24 页4.3 结果分析通过大量地仿真实验验证了基于小波变换地数字水印算法，该算法在