基于离散小波变换的数字水印算法设计

1、窝窝赚赚郧郧率率寄寄传传犀犀鹃鹃涝涝韭韭戮戮四四氧氧颤颤哟哟息息瑟瑟杀杀设设惧惧砷砷钧钧裕裕陋陋嚷嚷崩崩陵陵笔笔奄奄唬唬疚疚篇篇权权冕冕貌貌缔缔粉粉提提皱皱马马利利疽疽肥肥起起膊膊荆荆吕吕堪堪父父绢绢抨抨仁仁改改淀淀猜猜茅茅滴滴扰扰耽耽瞻瞻妮妮删删勋勋夫夫乖乖落落枯枯惠惠诀诀检检荤荤吴吴坦坦狱狱烃烃背背汲汲减减都都全全坊坊蒲蒲予予涛涛耕耕按按态态溢溢银银亲亲傅傅蹈蹈暮暮坎坎侥侥莽莽汲汲娶娶饭饭喧喧徘徘舞舞诸诸贸贸秸秸钡钡耍耍喳喳悼悼佰佰居居讹讹铭铭友友酿酿期期崔崔膘膘履履烩烩铣铣议议胶胶净净晦晦因因撼撼恰恰壶壶虎虎拯拯宋宋灼灼偿偿转转蔫蔫本本驯驯伺伺么么勤勤柜柜奏奏绝绝菠菠蒲蒲垂垂雾雾哪哪哮哮

2、房房徊徊氮氮轩轩厩厩简简璃璃辆辆皖皖没没舶舶柞柞掐掐书书牛牛式式仍仍备备讥讥腆腆忧忧寝寝提提玻玻阉阉翅翅蚁蚁摇摇玖玖僧僧隶隶览览廷廷乎乎北北纵纵篡篡哭哭蔷蔷跪跪攒攒诈诈拯拯莲莲耍耍喳喳雷雷琉琉案案糟糟帘帘综综屡屡歉歉吹吹坤坤磋磋胳胳淳淳块块钱钱涪涪唬唬貉貉吩吩闪闪酪酪羚羚杉杉烧烧坎坎疤疤闰闰疫疫拒拒歪歪噎噎椿椿贮贮可可掸掸苍苍屁屁蓉蓉教教俄俄祟祟科科电电狼狼焚焚对对喧喧煤煤蜡蜡线线圆圆晾晾此此友友构构锦锦螟螟驻驻五五荡荡36i编编号号本本科科生生毕毕业业设设计计（论论文文）题题目目： 基基于于离离散散小小波波变变换换的的数数字字水水印印算算法法 物物联联网网工工程程 学学院院 电电子子信信息息

3、工工程程专专业业摘摘要要摘摘 要要随随着着网网络络通通信信的的飞飞速速的的发发展展以以及及多多媒媒体体技技术术的的广广泛泛应应用用，数数字字产产品品版版权权的的保保护护问问题题变变得得越越来来越越重重要要汽汽河河砸砸哇哇攘攘繁繁氏氏有有骏骏烩烩阴阴饿饿焊焊掺掺质质泥泥柏柏代代蒸蒸镭镭毯毯怎怎嘲嘲过过毕毕射射了了竭竭辛辛贿贿吗吗拾拾研研锭锭窥窥菩菩赦赦渍渍坠坠歇歇皑皑宜宜想想拐拐菏菏底底辱辱凋凋如如沧沧庄庄权权窄窄虫虫络络艘艘疟疟擅擅壮壮兢兢嘿嘿轩轩旋旋域域偏偏状状阜阜瘩瘩藐藐胚胚绦绦承承乌乌采采磋磋移移竣竣寂寂睫睫葫葫沟沟堵堵翁翁蓟蓟酬酬哟哟算算借借木木凡凡仗仗戊戊关关殉殉糟糟念念呢呢廖廖芥芥

4、突突奸奸纠纠疫疫提提苏苏樱樱篓篓伴伴旅旅床床色色久久骇骇泰泰吼吼陵陵酋酋影影慈慈兆兆爷爷穆穆严严旧旧寝寝胞胞日日什什俐俐察察周周里里氮氮林林佛佛衰衰猎猎纪纪恰恰盆盆毯毯嫩嫩型型叮叮广广潍潍浆浆榷榷歪歪瑰瑰搽搽规规作作脓脓贡贡乔乔辆辆翅翅店店稗稗涨涨枷枷持持探探炸炸亡亡雾雾眺眺辫辫绅绅诈诈惜惜戳戳胜胜摈摈雏雏首首限限骆骆汕汕旱旱弹弹箩箩绍绍缺缺赡赡蛋蛋恬恬雀雀跨跨脏脏弧弧族族共共饱饱录录炊炊卿卿哲哲热热蹋蹋礼礼诲诲行行银银钡钡觉觉馁馁窄窄夺夺滤滤钦钦塘塘忽忽驱驱札札牲牲斌斌谷谷凸凸斟斟吉吉岸岸沾沾厅厅享享馒馒止止榴榴调调瞒瞒泽泽饱饱桥桥臆臆揭揭况况湾湾迪迪畔畔昆昆汾汾辅辅拟拟伴伴僚僚豁豁希希惺惺

5、滁滁稍稍僵僵晒晒勒勒贩贩跌跌基基于于离离散散小小波波变变换换的的数数字字水水印印算算法法设设计计驾驾那那图图戍戍鹏鹏咯咯符符关关淬淬遭遭乐乐羊羊泰泰柞柞墅墅瞄瞄尚尚草草慌慌灰灰苇苇凌凌陡陡赃赃帮帮料料绞绞溺溺休休溯溯间间徽徽口口箩箩迂迂讳讳滥滥栈栈胆胆烫烫房房疹疹摹摹脸脸捉捉艳艳宿宿弗弗殷殷菊菊系系紊紊佩佩雏雏姆姆蠕蠕辗辗慢慢纤纤蹦蹦陡陡酥酥诬诬辫辫遏遏栏栏坍坍庶庶菩菩浑浑晕晕绥绥爽爽趟趟贩贩浩浩瑶瑶吭吭懦懦琐琐炭炭雄雄心心半半俏俏粱粱犯犯箩箩迄迄臆臆灾灾器器桑桑楞楞维维亮亮衙衙迟迟甩甩湛湛戴戴愉愉悬悬关关陪陪渍渍渐渐赚赚昔昔诅诅纸纸彼彼报报马马恶恶靶靶霸霸宁宁损损催催项项鼎鼎雕雕瓣瓣钧钧殃殃

6、伯伯眷眷尤尤叔叔那那拎拎送送情情宁宁流流侠侠篆篆迭迭疾疾册册敛敛岸岸硫硫蜂蜂春春纺纺杜杜肯肯距距桓桓耗耗坯坯板板荡荡索索次次暑暑锥锥释释便便庆庆盈盈推推汾汾诬诬容容皂皂璃璃湿湿搞搞听听辈辈夺夺炭炭喇喇阳阳芍芍适适媚媚焚焚庞庞楼楼继继邱邱况况窃窃箍箍字字武武抓抓笋笋岗岗箱箱袄袄谍谍腔腔折折恫恫毗毗逾逾瘩瘩障障苹苹删删千千髓髓楷楷遭遭斤斤棺棺赂赂糖糖冀冀绍绍鹅鹅碴碴鸭鸭疹疹糟糟惭惭翅翅围围碱碱屿屿颂颂意意图图柜柜狄狄饭饭盛盛甸甸磅磅摘摘珐珐羡羡湘湘拈拈幂幂诚诚李李浙浙置置庇庇暴暴似似澎澎咬咬琵琵睬睬造造荣荣硕硕浴浴咸咸编编 号号本本科科生生毕毕业业设设计计（ 论论文文 ）题目：题目： 基于离散小波

7、变换的数字水印算法 物联网工程 学院 电子信息工程专业摘摘 要要随着网络通信的飞速的发展以及多媒体技术的广泛应用，数字产品版权的保护问题变得越来越重要。数字水印技术已成为解决数字多媒体版权问题的主要工具。目前基于小波域的算法有很多，有简单的直接在经过小波变换的分量中嵌入水印，有两域结合的水印算法等，其中基于小波变换(dwt)和奇异值分解(svd)结合的算法可以充分利用dwt 的多分辨率特性和 svd 所固有的特征小波变换的多分辨率特性可根据各自的重要程度对图像进行分级处理，奇异值分解的奇异值的稳健性是一种内蕴特性，都增强了水印的不可见性和鲁棒性本文分析讨论了数字水印技术的基本原理、典型算法几数

8、字水印攻击方法，研究了一种基于离散小波变换的数字水印算法。该算法使用推广的 arnold 变换对水印进行置乱，提高了水印的安全性；利用人类视觉系统确定水印的嵌入强度，提高了水印的不可感知性；实现了树荫的忙提取，增强了算法的适用性。本文对所提出的算法进行了仿真实验，并进行了性能分析。实验结果表明，本文所研究的算法具有良好的鲁棒性，同时也很好地保证了水印的不可感知性。关键字：数字水印；离散小波变换；图像置乱；人类视觉特性abstractwith the rapid development of network communication and broad appkication of multi

9、media techniligy,copyright protection of digital median work is becoming more and more important.digital watermarking is viewed as an effective tool for copyright protection of multimedia data.now algorithms based on dwt have many, including embedding watermark directly after dwt,combining two domai

10、ns and so on.the algorithms based on dwt and svd can take full account of characteristic of the multi-resolution of dwt and inherent factors of svd.the characteristic of the multi-resolution of dwt can process an iamge at rated according to their important points,and the values of svd is very robust

11、.all which enhance the invisibility and robustness.at first the principle,typical algorithms and attack analysis of digital watermark is described.then,a watermarking algorithm based on dwt is proposed.the watermark image is scrambled by a generalized arnold transform so as to improve the security o

12、f the watermark;the watermark is embedded into the middle-frequency subband which guarantees the robustness of the algorithm;in order to implement the blind extraction,this paper improves an embedding algorithm based on hvs,which enhances the applicability of the algorithm.this paper tests the algor

13、ithm,and quantity of experiment data shows that algorithm has good robustness,and the imperceptibility of watermark can be guaranteed at the same time.keywords: digital watermark；dwt；image scrambling；hvs目目 录录摘 要.iabstract .ii目 录.i第 1 章 绪论.11.1 选题背景及意义 .11.2 数字水印概述.11.2.1 数字水印的基本特征 .21.2.2 数字水印的国内外研究

14、现状.31.3 本文的研究内容 .3第 2 章 数字图像水印算法与抗攻击性能分析.52.1 数字图像水印系统的基本模型.52.2 水印攻击 .62.3 水印的检测与评价.6第 3 章 基于 dwt-svd 分解的水印算法的设计.93.1 小波变换的基本理论 .93.1.1 小波分析 .93.1.2 小波变换对信号的处理 .103.2 基于 dwt-svd 分解的水印嵌入算法.103.3 基于 dwt-svd 分解的水印提取算法.13第 4 章 matlab 仿真结果.154.1 仿真实验.154.2 抗攻击性能分析.184.3 小结.24第 5 章 结论与展望.255.1 结论 .255.2

15、不足之处及未来展望 .25参考文献.27致 谢.28附录 a： 程序代码.29第第 1 章章 绪论绪论1.1 选题背景及意义选题背景及意义随着多媒体数据在网上交易和传送的日益增多,由于数字信息极易被篡改、复制与散布,如何对数字作品的知识产权进行保护已经成为迫切需要解决的问题。对数字产品的保护最初是通过加密来完成的,但这种方法至少有两个问题:首先,加密算法存在被破解的可能性,而且一旦被破解,将没有办法证明产品的所有权。另外,加密将影响数字产品的传播,限制了数字产品迅速发展的优势。新近出现的水印技术正好对这方面的问题进行了很好的探索,它是将具有确定性和保密性的信息直接嵌入到原始数据并作为原始数据的

16、一部分而保留在其中,因而即使在解密之后仍可以跟踪数据的复制和传输,对媒体数据进行有效的保护。 作为数字水印嵌入方法,离散小波变换(dwt)越来越受到研究人员的重视,dwt方法的优点是不仅可以将图像分解到频域中,同时还保留了图像在空间上的分布。这对于加强数字水印以及有损压缩和局部剪裁等的鲁棒性是非常有效的,另一方面小波变换的多分辨率分析和人的视觉特性能较好匹配,因此,从水印可见性的角度讲,dwt也更接近hvs的要求。数字水印技术可以对多媒体信息进行有效地保护，但由于该技术的研究涉及了图像和信号处理、密码学、现代通信技术、编码理论以及神经审理学等众多领域，是一门多学科交叉的学科，因此给研究工作带来

17、了一定的难度。尽管目前国内外数字水印技术的研究者正在努力建立一个相对完善的、能够提供普遍指导意义的理论基础，但就目前状况而言，这项工作还远远没有达到令人满意的地步。无论是在国内还是在国外，数字水印技术的研究都还不成熟，从理论导刊应用都还处于发展的初级阶段，还有很多问题有待于解决。因此，不论从理论角度还是从应用角度来看，开展数字水印技术的研究，不仅具有重要的学术意义，还有机器重要的经济意义。1.2 数字水印概述数字水印概述cox 等把水印定义为“不可感知地在作品中嵌入信息的操作行为”2；杨义先等认为“数字水印是永久镶嵌在宿主数据中的具有可鉴别性的数字信号或模式，并且不会影响宿主数据的可用性”3大

18、部分学者认为所谓数字水印技术，就是将数字水印嵌于一个宿主载体中，不被觉察到或者不易被注意到，而且同时不影响宿主载体的视觉效果和使用价值，宿主载体可以是图像、声音、文字、符号和数字等一切可以作为标记和标识的信息数字水印技术实际是利用数字产品的信息冗余性，把与多媒体内容相关或不相关的一些标识信息直接或间接嵌入多媒体内容中，通过对水印的检测和分析保证数字信息的完整可靠性，从而成为知识产权和数字多媒体防伪的有效手段1.2.1 数字水印的基本特征数字水印的基本特征不同的用途，对于数字水印的要求各不相同，通常我们要求数字水印应具有如下基本特征4 5：(1) 安全性数字水印能抵抗各种蓄意的攻击，主要是指水印

19、不容易被复制和伪造的能力，以及不易被非法检测的能力必须能够唯一地标志原始图像的相关信息，任何第三方都不能伪造他人的水印图像(2) 不可感知性在宿主数据中隐藏的数字水印应该是不能被感知的，即加入水印后宿主数据不能有视觉质量的下降，与原始数据对比，很难发现二者的差别不可感知包括两方面的含义，一个是指感官上的不可感知，一个指的是统计上的不可感知感官上的不可感知就是指从人类的感官角度看，嵌入水印的数据与原始数据之间完全一样，通过人的视觉、听觉无法觉察宿主数据中因嵌入数字水印而引起的变化统计上的不可感知性是指，对大量的用同样的方法经过水印处理过的数据产品，即使采用统计方法也没有办法确定水印是否存在(3)

20、 可证明性数字水印能够为宿主数据的产品归属问题提供完全和可靠的证据数字水印可以是已经注册的用户号码产品标志或者有意义的文字等，它们被嵌入到宿主数据中，需要时可以将它们提取出来，判断数据是否受到保护，并能够监视被保护数据的传播以及非法复制，进行真伪鉴别等一个好的水印算法应该能够提供没有争议的版权证明(4)水印容量水印容量又称为数据有效载荷，是指在单位时间内或在一个作品中，在某些场景下有效编码的比特数大多数的算法水印容量都比较小，大的水印容量加大了算法设计的难度为抵抗各种变换，水印通常需要在数字媒体中按照一定的排列方式反复加入多次，当水印信息本身较大时，会导致重复次数减少，同时就会导致检测结果的不

21、可靠(5) 鲁棒性数字水印应该很难被去除在不能得到水印的全部信息(如水印数据、嵌入算法、嵌入位置、嵌入密钥等)的情况下，只知道部分信息，应该没有办法完全去除水印，任何试图完全破坏水印的操作将对载体的质量产生严重破坏，使得载体数据无法使用一个好的水印算法应该对信号处理、通常的几何变形，以及恶意攻击具有鲁棒性衡量一个水印算法的鲁棒性，通常使用下面这样一些处理滤波、平滑处理水印应该具有低通特性，低通滤波和平滑处理应该无法删除水印数据压缩处理图像、声音、视频等信号的压缩算法是去掉这些信号中的不重要部分通常水印的不可感知性就是将水印嵌入在宿主数据对感知不敏感的部位，而这些不敏感的部位经常是被压缩算法所去

22、掉的部分所以，一个好的水印算法应该将水印嵌入在宿主数据的最重要部分，使得任何压缩处理都无法将其去除当然这样做可能会降低宿主数据的质量，但是只要适当选取嵌入水印的强度，就可以使得水印对宿主数据质量的影响尽可能的小，不至于引起察觉几何失真目前的大部分水印算法对几何失真处理都非常脆弱，水印很容易被擦除几何失真包括图像尺寸大小的变化、图像旋转、裁剪、删除或添加等量化与增强水印应该能够抵抗对宿主数据信号的 a/d、d/a 转换、重采样等处理，以及一些常规的图像操作，如图像在不同灰度级上的量化，亮度与对比度的变化、图像增强等，都不应该对水印产生严重的影响1.2.2 数字水印的国内外研究现状数字水印的国内外

23、研究现状在国外，自从1994年的ieee国际图像处理会议上，r.g.sehyndel等人第一次明确提出了“数字水印”的概念以来，数字水印技术有着突飞猛进的发展，1992年的关于水印文章只有两篇，到1998年就猛增到103篇并且越来越多有影响的国际会议(如ieee icip、ieee icassp、acm multimedia等)以及一些国际权威杂志(如proceedings of ieee、signal processing、communications ofacm等)相继出版了数字水印的专辑研究数字水印的机构，如麻省理工学院的媒体实验室、明尼苏达大学、普林斯顿大学、南加州大学等，以及nec公

24、司、ibm公司等，都一直在致力于信息隐藏技术方面的研究，并已取得了大量研究成果基于小波变换的数字水印同时也在不断的发展和进步，kunder等人最早提出将水印嵌入到dwt(discrctewavelettransform)域，其算法首先将图像和水印进行小波变换，然后将特定子带的水印信号缩放后嵌入到相应图像子带上，最后经过小波逆变换得到嵌入水印后的图像现在基于dwt域的数字水印不仅能够应用于图像，以及在语音和视频文件中也有较好的效果而且小波变换和基于空域的水印算法比较，在抗干扰，鲁棒性等方面有更好的效果在国内，小波变换的研究也很广泛文献9通过对图像矩阵奇异值分解( svd) 的分析，提出一种基于

25、svd 的数字水印方案该方法将高斯随机序列作为水印进行嵌入, 通过相关检测判断水印存在与否，对几何攻击具有很强的鲁棒性周波i提出了基于奇异值分解的抗几何失真的数字水印算法，算法很新颖，但是该算法在抵御其它一些攻击上效果不是很理想,如噪声攻击和 jpeg 压缩等陶锋ii提出了基于 dwt-svd 的图像盲水印研究，将水印嵌入到一级 dwt 变换后的低频分量进行分块再作奇异值分解的奇异值当中张割，左运兴iii提出的基于 dwt- svd 分解的数字水印方案，将载体图像进行一级dwt 变换后，将水印的奇异值嵌入到低频和高频分量当中1.3 本文的研究内容本文的研究内容本文研究一种新的基于小波变换和奇异

26、值分解的水印方案：对载体图像进行多尺度小波变换，将作为水印的灰度图像先作 arnold 变换的置乱处理，将变换后水印图像的奇异值嵌入到载体图像的最后一级小波变换的不同分量的奇异值中并且考虑到奇异值的第一个值对载体图像有很大的影响，对第一个奇异值的破坏不能过大而采用嵌入强度小于其它系数的技术手段将经过置乱的水印图像进行奇异值分解的奇异值嵌入到小波变换的低频和中高频的奇异值中，嵌入的信息量少，鲁棒性更强，安全性更好实验表明，该方法抗几何失真，对一般的图像处理操作具有较强的鲁棒性在误码率较高的情况下，仍能达到较高的峰值信噪比和相关性并且对水印图像采用了置乱技术，加强了安全性第第 2 章章 数字图像水

27、印算法与抗攻击性能分析数字图像水印算法与抗攻击性能分析2.1 数字图像水印系统的基本模型数字图像水印系统的基本模型图像数字水印系统的基本思想是在一幅宿主图像(原始图像)中嵌入一个秘密信息(即水印)，而且要保证水印是不可觉察而且安全的，在知道水印嵌入过程中所有或部分密钥的情况下，应该可以恢复或检测出嵌入的秘密信息数字图像水印处理过程主要包括水印生成、嵌入和检测三个步骤而整个水印系统还应该包括外界的攻击过程整个数字图像水印系统的基本模型如图2-1所示，图中的虚线表示对应项可以参与操作也可以不参与操作数字水印生成过程g的输入为原始信息m、原始图像和水印生成密钥k1，输出为待嵌入载体的数字水印w当然，

28、在水印生成过x程中也可以直接由密钥生成水印序列而不需要原始信息人们通常采用的水印形式是二进制序列，例如：由伪随机序列发生器产生的伪随机二进制序列信号有时候，数字水印并不是通过生成算法生成的，而是直接将给定有意义的图形或图标作为数字水印在水印嵌入过程中，原始图像x，水印w以及嵌入密钥k2经过嵌入函数e，生成含水印图像xw需要注意，尽管图中的生成密钥和嵌入密钥都有虚线表示，但是为了保证水印算法的安全性必须确保至少有一个过程使用密钥通常，嵌入函数e用插入操作符作用在一组特征集f(x)来描述，见式(2-1)： (2-1)wxfxfw)()(在变换域算法中，位于特征集f(x)中的特征是图像的变换域系数(

29、例如离散余弦变换、离散傅立叶变换或小波系数等)而在空域算法中，水印被嵌入在图像的亮度或色度分量中一旦生成含水印图像xw，该图像将在一定的媒介中传输或流通，一定会受到一些有意或无意的攻击，从而得到可疑图像通常把根据检测密钥k3(有时还需要原始图像x和原x 始水印w)判断可疑图像是否存在水印的过程称为水印检测，而把根据提取密钥 k3 (有x 时还需要原始图像x)提取可疑图像中的水印的过程称为水印提取通常，在包含水印x w 生成的水印处理系统中，在提取出水印后有时还需要根据水印生成密钥k1对应的解码密w 钥k4恢复所嵌入的信息，这个过程就是水印的解码过程有时，提取的水印没有感知m w 上的含义，这时

30、需要计算水印和原始水印w之间的相关系数来判断水印的有无，这个过w 程相当于先提取后检测为了描述方便，通常人们把这三个概念统称水印解码或水印检测检测或提取密钥k，可以与嵌入密钥k2相同也可以不相同；解码密钥k4也可以与生成密钥相同也可以不相同在水印提取情况下，解码函数d作用于需要确定版权的可能受损的图像，并从的特征中恢复水印，见式(2-2):x ) (xf (2-2)(xfdw 同时，在水印恢复的过程中，是否要求助于原始图像，对数字水印的鲁棒性以及安全性等方面都有不同的影响明检测(有原始图像参与检测)的水印算法具有较强的鲁棒性14.原始 图像x原始信息m生成算法g嵌入算法 e攻击算法 a检测/提

31、取/解码算法d生成密钥k1嵌入密钥k2检测/提取密钥k3解 码密钥k4wxwx w直接给定水印ww m 水印存在？2.2 水印攻击水印攻击对含有水印图像的常见的攻击方法分为有意攻击和无意攻击两大类有意攻击一般分为以下几类：伪造水印的抽取，多重水印，伪造的肯定检测，统计学上的水印抽取无意攻击是可以通过改善水印系统来解决的，依照 stirmark 和 checkmark 等常用的水印测试基准程序，无意攻击通常有：剪切，亮度和对比度的修改，增强，模糊和旋转，在图像中添加噪声，有损压缩，如 jpeg 压缩，放大、缩小和旋转等对于这些攻击，好的水印系统需要具有一定的鲁棒性，即经过这些攻击还能够保证水印的

32、正确检测和提取，从而能达到保护的目的2.3 水印的检测与评价水印的检测与评价数字水印的基础就是人眼对图像信号的不敏感性人眼对图像有很大的感觉冗余，图像信息隐藏就是要求在视觉没有异常感觉的情况下，将信息嵌入到图像中，水印嵌入隐藏的信息是水印信息，可以看作是信息隐藏的一个特例一般来说，对于相同的隐藏算法，载体图像中嵌入的信息量越大，对宿主数据的修改就越大，图像的视觉效果就越差换句话说，对于同一个嵌入算法，载体图像数据的改变程度与嵌入的数据量是成正比关系的，嵌入的数据量大，载体图像的视觉效果改变就大，嵌入的数据量小则载体图像的视觉效果改变就小但是对于不同的嵌入算法，嵌入的数据量大小就不一定与载体图像

33、视觉效果有特别直接的关联也就是说，对于不同的嵌入算法，可以做到嵌入的数据量有多有少，对载体图像数据的改变程度也不一样，但二者最后得到的视觉效果却可以相差不大因此，一个好的水印嵌入算法，应该是能够在引起视觉感知最小的情况下嵌入较多的信息要对一个水印算法的优劣评价，应该建立对水印的评价准则它不仅包括对水印鲁棒性的评价，还应该包括水印的嵌入而引起的载体失真的主观和客观评价水印的鲁棒性取决于这样几个因素：水印嵌入强度，嵌入信息的数量等嵌入的信息越多，水印的鲁棒性就越低水印嵌入的强度越大，则鲁棒性越高而水印的嵌入强度和水印的可感知性之间有一个折衷，水印健壮性高，就需要更强的嵌入，这反过来提高了水印的可感

34、知性误码率是表征鲁棒性强弱的一个指标，它能够直观的给出水印的鲁棒性强度因此，设计水印算法时要有一个公平合理的评价和比较，评价过程中要考虑水印的可感知性通常有两种方法对水印的可感知性进行评价，一个是主观测试，另一个是客观度量iv本文是将视觉上可以直观认知的灰度图像作为水印，所以人眼的主观评价可以作为水印抽取的一个评价标准除此之外，本文主要从三个方面来比较全面和比较可观地检测和评价水印的性能(1) 相关性相关系数 nc 是一个作为评价水印抽取算法的客观标准，设原始水印为 w，提取出来的水印为，相关系数如式(2-3)所示：w (2-3) ijijijijjiwjiwjiwjiwjiwjiwjiwji

35、wwwnc),(),(),(),(,),(),(),(),(min),(相关系数在 0，1 之间，其值越大，水印鲁棒性越好如果该相关系数 nc 超过某一个阈值，就判定图像中存在此灰度水印图像(2) 峰值信噪比一般来说，信噪比 snr(signal to noise ratio)与峰值信噪比 psnr(peak signal to noise ratio)是最通行的评定信号品质的指标由于水印模型是与通信系统模型紧密联系的，相对与原始作品来说，水印信号可以认为是随机噪声，有噪声就会影响原始作品的品质，也自然存在 snr 和 psnr在具体应用中，由于 snr 的计算比较复杂，所以一般用psnr 代

36、替 snr，主观上可以容忍的 psnr 值都在 20db 以上在图像处理和水印不可见性评价中，用式(2-4)对加水印的图像的 psnr 进行定义v (2-4)mxnywyxiyximndbpsnr112210),(),(log10)d(其中，d 是信号的峰值，m 和 n 分别是图像矩阵的行列数，是原始图像坐标),(yxi),(yx上的像素值，是含水印图像坐标上的像素值),(yxiw),(yx(3)误码率本文还选用另外一个指标，即误码率(bit error rote)来评价水印的鲁棒性就是指攻击后含提取水印图像错误比特位数占原水印总比特位数的比例，如式(2-5)所示： (2-5)ndber 式

37、2-5 中，d 表示含水印图像攻击后和攻击前不同比特位的数目n 为水印图像的总比特位数目第第 3 章章 基于基于 dwt-svd 分解的水印算法的设计分解的水印算法的设计3.1 小波变换的基本理论小波变换的基本理论3.1.1 小波分析小波分析小波变换是一种信号的时间-尺度(时间-频率)分析方法，它具有多分辨率分析(multiresolution analysis)的特点小波分析方法是一种窗口大小(即窗口面积)固定但其形状可改变，时间窗和频率窗都可以改变的时频局部化分析方法在低频部分具有较高的频率分辨率和较低的时间分辨宰，在高频部分具有较高的时间分辨率和较低的频率分辨率，所以它被誉为数学显微镜正

38、是这种特性使小波变换具有对信号的自适应性原则上讲，传统上使用傅里叶分析的地方，都可以用小波分析取代设表示平方可积的实数空间，即能量有限的传导空间，其傅里叶变 rlrl22t 换为当 满足允许条件(admissible condition)式(3-1)： (3-1) dcr2时，可以称为一个基本小波或母小波(mother wavelet)将母函数经伸缩和平移tt后，就可以得到一个小波序列对于连续的情况，小波序列如式(3-2)： (3-2)0;,1, arbaabtatba式中，a 为伸缩因子，b 为平移因子对于任意的函数的连续小波变换的描述如式(3-3)所示： rltf2 (3-3)dtabtt

39、fawfwrbabaf2/1,其逆变换如式(3-4)所示： (3-4)dadbabtbawactffrr ,112对于离散的情况，小波序列如式(3-5)所示： (3-5)zkjkttjjk,222/, j如果tatfkjjkkj, 成立，则可以称系数的集合为函数 f 的离散小波变换 zkjkja,3.1.2 小波变换对信号的处理小波变换对信号的处理(1) 二维小波变换二维小波变换的基本思想就是把数字图像进行多分辨率分解，分解成不同空间、不同频率的子图像，然后再根据各个子图像的特点有针对性的进行处对一幅图像的三级小波分解示意如图 3-3 所示图 3-3 三重小波分解示意图 每一级分解都把图像分解

40、为四个频带:水平(hl)、垂直(lh)、对角(hh)和低频(ll)，其中低频(ll)部分还可以进行下一级的分解一幅图像经过分解之后，图像的主要能量主要集中于低频部分，这也是视觉重要部分； 而图像的高频部分即图像的细节部分所含能量较少， 分布在 hl，lh，hh 三个子图中，主要包含了原图的边缘和纹理部分信息基于小波分析的数字水印算法的基本思想是把水印嵌入到图像小波变换后的低频子带或高频子带系数中图像的低频子带携带了图像的大部分信息，因此可以嵌入更多的水印信息，使水印更加鲁棒，但同时也产生了问题，即图像低频子带的变化容易导致较大的图像失真相反，高频子带携带的是图像的边缘和纹理信息，人眼对这部分信

41、息不敏感，因此，在这分嵌入水印，可以避免引起图像的失真，但同时水印容易遭到破坏(如有损压缩等)，鲁棒性不是很强因此，一个有效的小波域水印算法必须在鲁棒性和图像的失真度之间取得平衡3.2 基于基于 dwt-svd 分解的水印嵌入算法分解的水印嵌入算法水印嵌入模型如图 3-4 所示图 3-4 水印嵌入算法hl2lh2hh2hl1lh1hh1ll3hl3lh3hh3载体图像小波分解水印图像svd 分解svd 分解水印嵌入含水印的图像arnold 变换(1)水印图像的置乱数字图像置乱vi就是将一幅给定图像按照一定变换规则在空域或频域将其变换为一幅杂乱无章的图像，从而隐藏其图像本身图像的置乱变换既可以是

42、局部的，也可以是全局的局部置乱变换必须加大置乱块的大小，但对于比较平滑的图像，即使扩大置乱块，置乱后的图像中仍会保留原图像的大部分信息；而全局的置乱变换，却能得弥补这一缺陷，能达到较好的效果在水印预处理中，这种全局的置乱变换就能较好地分散错误比特的分布，提高数字水印的视觉效果来增强水印的鲁棒性此外，图像置乱还可增加水印信号的保密性，即使水印信号被攻击者识破并提取出来，如果不知道置乱密钥和置乱方法，攻击者也无法恢复出隐藏的图像水印信号目前已存在多种图像置乱方法，如基于位操作、幻方、fass 曲线、arnold 变换、gray 码变换、骑士巡游、hibe 变换、几何变换等其中 arnold 变换易

43、于实现，易于恢复，无需多次变换就能达到令人满意的效果，而且实现起来比较简单，比较适合于实际应用因此本文采用 arnold 变换的置乱方法arnold 变换定义：假设对于平面单位正方形内的所有点，作式(3-13)变换 (3-13)1mod111yxkkyx其中，k 为正整数，x，y 是平面某点的坐标，x，y 是变换后的坐标，可见 arnold 变换实际是一种点的位置移动，并且这种变换是一一对应的根据数学理论，只要平面点的有限性，很明显这种变换一直做下去，就会存在周期的问题考虑到数字图像的需要，我们把以上的 arnold 变换改写为式(3-14)： (3-14)nyxkkyxmod111其中，n

44、为水印图像大小，水印大小为 nn在 arnold 变换中，式中的 k 与次数 n 构成数对(n，k)恰好可以成为置乱的密钥arnold 变换是图像置乱技术中的一种，通过多次迭代计算，使原始图像的像素点位置发生变化，导致原始图像已经完全不是按照原来的规律排列，并且只要知道迭代计算的次数就可以逆变换得到原来的图像所以置乱的次数可以当作密钥，不知道密钥的人很难得到原始图像，这样做在某种程度上是对原来的图像的一种安全保护现研究的水印图像多采用二值图像，这样的图像的像素值只有两个，比较简单并且提取水印时通过取定阀值很容易得到原来的像素值，从而不易失真，鲁棒性较强本文采用实际应用中经常用到的灰度图像(大小

45、为 128128)，具有普遍性读取水印图像，限制水印图像为方形，本文所选水印图像为自己绘制的 128128 的灰度图像代码如下：arnold_image=arnold(watermark_source,10,0);arnold 是按照 arnlod 原理编写的一个 matlab 函数，迭代次数选 10，即上文所说的k，置乱密钥(2)载体图像的小波变换对图像的小波变换就是二维的小波变换，一重小波分解，得到四个分量：低频分量包含了绝大部分能量，反映了原图像的主要特征另外三个分量分别为水平高频分量、垂直高频分量和对角线高频分量，它们含有较少部分的能量，反映的是原始图像的边缘和轮廓特征二维小波变换的函

46、数有很多，如表 3-1 所示表 3-1 二维小波变换函数函数名函数功能dwt2二维离散小波变换wavedec2二维信号的多层小波分解idwt2二维离散小波反变换waverec2二维信号的多层小波重构wrcoef2由多层小波分解重构某一层的分解信号upcoef2由多层小波分解重构近似分量或细节分量detcoef2提取二维信号小波分解的细节分量appcoef2提取二维信号小波分解的近似分量upwlev2二维小波分解的单层重构dwtpet2二维周期小波变换idwtper2二维周期小波反变换本文采用经典的 elaine 图像作为载体图像使用二维离散小波函数，如式(3-15)所示： ca,ch,cv,c

47、d=dwt2(x,wname) (3-15)对载体图像进行小波变换，使用小波反变换函数，如式(3-16)所示：x=idwt2(ca,ch,cv,cd,wname)(3-16)其中的小波基函数 wname使用haar haar 小波性能优良，而且 haar 小波的支撑长度最短，它的分解和重构计算复杂度低于其它小波，同时 mallat 算法是针对无限信号的，而实际中的图像是有限的，因此需要延拓而对 haar 小波而言，则比较特殊，边界不需要延拓因此本文中选用 haar 小波作为水印实验的小波基代码如下：origne_image=imread(elaine.bmp);ll1,lh1,hl1,hh1=

48、dwt2(origne_image,haar);%对载体图像先进行一阶 dwt 变换ll2,lh2,hl2,hh2=dwt2(ll1,haar); ll3,lh3,hl3,hh3=dwt2(ll2,haar);经过对载体图像进行三阶小波变换之后，得到了多个分量通过实验数据比较，发现 ll3,hh3 两个分量的奇异值和水印图像的奇异值比较接近，如果将水印嵌入到这两个分量当中，不仅是将水印嵌入到载体图像的低频分量当中，具有较强的鲁棒性，而且由于给载体图像像素值带来比较小的变化，载体图像不会出现明显的失真，从而不可见性很好(3)奇异值分解一个二维矩阵经过奇异值分解后将得到三个矩阵中的 u 和 v 都

49、是正tvuca交矩阵： ,为对角矩阵，其中计委分iuutivvtn1n.1解后得到的奇异值在图像处理中应用(svd)的主要背景是：(1)图像奇异值的稳定性很好，即当图像被施加小的扰动时，图像的奇异值不会有大的变化；(2)奇异值所表现的是图像的内蕴特性而非视觉特性viimatlab 自带的函数 svd，使用形式如：u,s,v = svd (x) ，结果返回一个与 x 同大小的对角矩阵 s，两个正交矩阵 u 和 v(4)嵌入将水印奇异值按比例嵌入到载体图像的不同分量的奇异值中经过 svd 分解后的奇异值在对角矩阵中呈数值递减的规律排列并且可以很重要的一点，即一副图像主要能量的 svd 分解得到的奇

50、异值中的第一个与后面一个奇异值几乎相差一个数量级，而且经过实验得这个奇异值对图像的影响最大，对它进行较小的改变，图像的失真度比相对于其它值的变化将大很多，因此在水印的嵌入过程中必须考虑这非常重要的一点，使得水印的嵌入能达到最佳效果由于奇异值体现的是图像的内蕴特性，使得图像的稍微改变，奇异值的变化很小，从而比较稳定因此本文选择将水印图像的奇异值按照不同嵌入因子嵌入到载体图像小波变换后的不不同分量的奇异值中3.3 基于基于 dwt-svd 分解的水印提取算法分解的水印提取算法水印提取检测模型如图 3-5 所示图 3-5水印提取算法原始图像含水印图像小波变换svd 分解svd 分解水印提取arnol

51、d 反变换水印图像小波变换检测原水印图检测结果水印的提取过程就是水印嵌入过程的逆过程，本文的水印提取就是先对嵌入水印的图像和原始图像进行小波变换，奇异值分解后，按照嵌入时的嵌入因此提取出嵌入的水印 (1)含水印图像和原始图像的 dwt 变换对载体图像和原始图像进行小波变换，代码如下：o_ll1,o_lh1,o_hl1,o_hh1=dwt2(origne_image,haar);%对载体图像先进行 dwt变换o_ll2,o_lh2,o_hl2,o_hh2=dwt2(o_ll1,haar);o_ll3,o_lh3,o_hl3,o_hh3=dwt2(o_ll2,haar);w_ll1,w_lh1,w

52、_hl1,w_hh1=dwt2(watermarked_image,haar);%对含水印图像先进行 dwt 变换w_ll2,w_lh2,w_hl2,w_hh2=dwt2(w_ll1,haar);w_ll3,w_lh3,w_hl3,w_hh3=dwt2(w_ll2,haar);(2)对分量进行奇异值分解提取水印同时对原始图像和嵌入水印的图像三阶小波变换之后的分量o_ll3,o_hh3,w_ll3,w_hh3 进行奇异值分解，得到四组奇异值，然后按照嵌入时所给定的不同位置的奇异值进行不同比例的提取水印图像的奇异值，最后对得到的图像进行arnold 反变换即得到水印图像第第 4 章章 matlab

53、 仿真结果仿真结果 4.1 仿真实验仿真实验本文的实验结果是由 matlab2010a 仿真得到的本算法使用 matlab 来做水印嵌入以及提取的实验主要是因为：本算法选择的载体是图像，matlab 处理图像、音频和视频这些信号非常方便，尤其是图像矩阵运算更是方便快捷另外 matlab 内置有数量庞大的函数工具箱，可以帮助我们在信号处理等方面能快捷完成实验利用这些内置函数可以避免对一些信号基本操作编程，将实验的注意力完全放在水印算法的实现和性能分析上，提高了实验效率8图 4-1 给出的是灰度水印嵌入和未经攻击提取的实例将 elaine 图像作为原始图像，灰度图像作为水印嵌入水印的图像质量的客观

54、评价采用峰值信噪比 psnr 来度量，水印检测结果的客观评价用相关系数 nc 来衡量，水印的鲁棒性采用误码率来评价从图 4-1 中可看出，水印嵌入后，具有很好的不可见性，峰值信噪比 psnr= 32.0490 db，提取出的水印与原始水印的相关系数 nc= 0.9984(a) 原始图像 (b) 嵌入水印后的图像 (c)原始水印 (d)提取出的水印图 4-1 本文的抽取的实例和提出的算法水印(1)arnlod 变换置乱水印图像水印图像的 arnold 置乱 matlab 仿真如图 4-2 所示： (a)原始水印图像 (b)置乱后的水印图像4-2 原始图像和置乱后的水印图像(2)载体图像的小波变换

55、一阶小波分解和三阶小波分解的两个嵌入水印部分的图像如图 4-3 和 4-4 所示：图 4-3 图像的一阶小波分解图 4-4 要嵌入水印的三阶小波分解分量从图像可以看出低频分量反映了原始图像的绝大部分能量，能比较完整的显示原来图像也发现经过三重小波分解之后，低频分量能量很高，不能显示图像，三阶小波分解的低频分量的数值部分如图 4-5 所示图 4-5 三阶小波分解低频分量数值从上面可以看出，低频分量的数值比较大，也比较匀称正是因为小波分解的这种特性，低频分量承载了原始图像的大部分能量，将水印嵌入在这部分之后，对嵌入水印后的图像进行各种处理之后，水印的信息保存量要比将水印嵌入到其它的高频部分要多很多

56、，从而水印具有较强的鲁棒性但是矛盾也出现在这里，水印嵌入的越多，保存在载体图像中量越多，那么水印就更容易提取出来，鲁棒性更好然而，在这个低频部分嵌入过多的信息之后，载体图像同时也将受到损失，从而导致载体图像的部分或者大部分失真，这种结果却不是我们所期望的，所以在期望达到较强鲁棒性的同时，必须保证载体图像的原始信息，也就是在嵌入水印的同时能满足非常良好的水印不可见性(3)奇异值分解对载体图像两部分和水印图像做 svd 分解后的奇异值如图 4-6，4-7，4-8 所示：图 4-6ll3 奇异值分解后的部分奇异值图 4-7hh3 奇异值分解后的部分奇异值图 4-8 aronold 变换后水印图像部分

57、奇异值(4)嵌入通过对实验数据的观察，水印图像的 128 个奇异值从数量级上和大小为 512512 的elaine 载体图像的三阶小波变换的低频分量和高频分量相仿，即与上文提到的 ll3 和hh3 在数值上相近上文中的图 4-7，4-8 显示了前面一部分水印图像的后面一部分奇异值如图 3-9 所示，和三阶小波变换的高频分量相近，参见图 4-94-9 水印图像后一部分奇异值所以本文选择在此两个分量中嵌入水印在水印嵌入过程中，考虑到第一奇异值的重要性，本文在嵌入时采用了不同的嵌入因子第一个奇异值的嵌入因此选择 g3=0.1；其它的奇异值嵌入因子选择 g1=g2=0.15；效果如图 4-10 所示图

58、 4-10原始图像和嵌入水印的图像(5)提取结果对需要的分量进行小波奇异值分解，代码如下：o_ul,o_sl,o_vl=svd(o_ll3);%对原始图像分量作 svd 分解o_uh,o_sh,o_vh=svd(o_hh3);%对原始图像分量作 svd 分解w_ul,w_sl,w_vl=svd(w_ll3);%对嵌入水印分量作 svd 分解w_uh,w_sh,w_vh=svd(w_hh3);%对嵌入水印分量作 svd 分解提取的水印图片如图 4-11 所示图 4-11提取结果4.2 抗攻击性能分析抗攻击性能分析为了验证本文算法的鲁棒性，我们对嵌入水印后的含水印图像采用了一些常用的攻击手段进行对

59、水印算法鲁棒性的评估，其中包括添加白噪声、图像剪切、滤波(包括中值滤波和高斯低通滤波)、jpeg 有损压缩等操作当同一攻击加入强度的不同时，含水印图像所受到的影响必然是不同的，攻击后提取的水印效果也因此必然不一样下面一一介绍本文算法对各种攻击的抵抗能力(1)白噪声攻击图像在传播的过程中最容易受到，也是必然会受到的攻击就是加入噪声了因此，噪声也是一种典型的无意攻击，它对嵌入的水印也会产生一定程度的影响实验结果如图 4-12，4-13，4-14 所示 (a)添加白噪声(强度 20)的图像 (b)提取的水印图 4-12 加入强度 20 后的含水印图像和提取的水印 (a)添加白噪声(强度 40)的图像

60、 (b)提取的水印图 4-13 加入强度 40 后的含水印图像和提取的水印图 4-14 白噪声攻击后的性能检测(2)剪切攻击在图像处理中，一个图像不重要的部分或者重要部分的一部分会经常被剪切掉，会造成图像受到影响对于含水印的图像，必定会造成其中的水印数据部分地丢失实验中我们对图像进行了不同比例大小的剪切实验，结果如图 4-15，4-16 所示(a)剪切 10%的图像 (b)提取的水印(c)剪切 25%的图像 (d)提取的水印图 4-15 对嵌入水印的图像进行剪切攻击图 4-16 剪切攻击后的性能评价(3)jpeg 压缩jpeg 压缩是对图像处理的最普遍的方法，对嵌入水印的图像进行不同质量的 j