基于MATLAB的数字水印系统

素颜马尾好姑娘i 2022-09-11 10:10 490阅读 0赞

## 基于MATLAB的数字水印系统 ##

## 一、课题背景和研究意义 ##

从二十世纪九十年代开始，网络信息技术在全世界范围内取得了迅猛发展，它极大方便了人们之间的通信和交流。借助计算机网络所提供的强大的多媒体通信功能，人们可以方便、快速地将数字信息（数字音乐、图像、影视等方面作品）传到世界各地，一份电子邮件可以瞬息间传遍全球\[1-7\]。但在同时计算机网络也成为犯罪集团、非法组织和有恶意的个人利用的工具。从恶意传播计算机病毒，到非法入侵重要部门信息系统，窃取重要机密甚至可使系统瘫痪；从用计算机进行金融犯罪，到利用表面无害的多媒体资料传递隐蔽的有害信息，对计算机信息系统进行恶意攻击的手段可谓是层出不穷。

信息隐藏技术在保密通信、版权保护等领域中都具有相当广泛的应用价值，根据不同应用背景，信息隐藏技术可以分为两个重要分支隐写术（Stegano Graphy）和数字水印\[8\]（Digital Watermarking）。数字水印主要是为了保护知识产权的，通过在原始媒体数据中嵌入信息来证实该媒体的所有权归属问题。数字水印的主要目的不是限制对媒体访问，而是确保媒体中水印不被篡改或消除。因此稳健性是数字水印最基本要求之一。数字水印中的稳健性是指水印图像经过一些常见改变后，水印仍具有较好可检测性。这些改变包括常见的图像处理手段（如数据压缩、低通滤波、图像增强、一次抽样、二次量化、A/D和D/A转换等等）、几何变换和几何失真、噪声干扰、多重水印（Multiple Watermarking）的重叠等。对不同应用场合，要求具有不同的稳健性。但需要指出的是，存在另一种与稳健水印性质相反的水印，被称为易损水印（Fragile Watermarks），它们被用来证实原始媒体是否被改动过。稳健性在整个水印系统设计中具有很重要的分量，这也是将隐写术和数字水印分别对待的原因之一。隐写术主要考虑的是安全性（即统计特性上无法检测隐密信息的存在性）和嵌入容量，不可见性等。

因此，信息隐藏技术的研究在信息安全领域中具有很重要的地位，它对于军事、情报、国家安全方面的重要意义是不言而喻。它包括数字隐写与隐写分析两个方面。一方面要以尽可能隐蔽的方式将信息深藏于浩如烟海的数字多媒体信号之中，毫不引起对方怀疑而达到隐蔽通信的目的；另一方则以各种手段检测可疑信息的存在，搜寻敌对隐蔽通信的信源，隔断隐蔽通信的信道。设计一种高度安全的隐写方法是一项富于挑战性的课题，而对隐写的准确性判定往往比隐写本身更加困难\[12\]。数字隐写与隐写分析的交互发展正是方兴未艾，成为互联网时代信息战技术的一个新的课题。信息网络上的攻防技术水平将反映的是一个国家的科技水平和防范意识。

## 二、LSB算法介绍 ##

最低有效位（Least Significant Bits，LSB）\[18\]方法是最早提出来的最基础的空域图像信息隐藏方法，其它的许多空域算法都是从它的基本原理进行改进扩展得到的，使得LSB方法成为使用最广泛的隐藏技术之一。现在有一些简单信息隐藏软件大多是运用LSB和调色板调整等相关技术将信息隐藏在24 bit图像或256色图像之中，如Hide and Seek，Stego - Dos，White Noise Storm，S-tools等经典信息隐藏软件。

1、隐藏原理

LSB方法是通过调整载体图像像素值的最低若干有效位来实现数据的嵌入方法，使所隐藏信息在视觉上很难被他人发觉，而且只有知道了秘密信息嵌入的位置才能正确提取出秘密的信息。显然，LSB隐藏算法最低位被改变的概率大概是50%，它在原始图像里面引入了极小噪声，在视觉上是不可见的。实际上，对于24 bit真彩色的图像，我们在其最低的两位甚至三位来隐藏信息在视觉上仍然是不可见的，对于灰度图像，改变其最低的两位也能取得较好的效果。

另外，在LSB方法之中，也可不采用直接嵌入的方法，根据异或的可逆准则性，采用替换的准则来实现信息的隐藏效果。异或的简单原理如下：

![v2-9ce5ca879497f5f6dec502c2ff617b71\_b.jpg][v2-9ce5ca879497f5f6dec502c2ff617b71_b.jpg]

。因此在嵌入数据位之时，嵌入的是数据位与1或0的异或值。基于异或的运算也有许多可改进的算法，在嵌入的过程之中，首先计算出每个像素灰度值的每一位的异或值，并把所得到结果与要嵌入的信息进行异或运算，然后把像素灰度值的最低位全部清零或置为1，再根据异或运算结果的值来改变最低位信息，实际上，这相当于对信息进行了一层加密的处理，嵌入的不再是原始信息了，而是原始信息的另外一种表达形式，不知道密钥的攻击者很难从中提取出有用信息。

2、鲁棒性分析

LSB算法具有相当弱的鲁棒性。对于许多变换，即使是有益的，也是很脆弱的。有损压缩典型如JPEG压缩，就有可能彻底破坏隐藏的信息。因为LSB算法试图利用人类的视觉系统上的漏洞，而有损压缩算法所依赖的是对附加噪声的不敏感性，正是利用它来减少数据量的。几何变换的移动像素尤其是改变像素在原栅格中的位置都有可能破坏嵌入的信息。任何其它的图像的变换如模糊、滤波等，通常都也会破坏隐藏的数据。

3、隐蔽效果

图2.1为512×512的Lena的灰度图，图2.2为LSB直接替换的最低有效位后的Lenna图，我们看不出两个图片之间的区别，隐蔽性比较好。

![v2-b3ede55b59d8e11c7f668aa97c42e722\_b.jpg][v2-b3ede55b59d8e11c7f668aa97c42e722_b.jpg]

图2.1 原来的Lenna灰度图 图2.2 Lenna的隐秘灰度图

## 三、 水印攻击性实验 ##

由于数字水印在实际应用中可能会遭到各种各样的攻击，因此对算法进行攻击测试是衡量一个水印算法优劣的重要手段。

1、JPEG 压缩实验

首先对嵌入水印后的图像进行 JPEG压缩(如图3.3)，而后从压缩的图像中提取出水印(如图3.4)，对融合图像进行不同JPEG压缩比压缩，得出数据如表3.1所示。程序代码如下：

imwrite(Y,'watermarked1.jpg','jpeg','quality',45);

Y=imread('watermarked1.jpg');

![v2-5157db68a32a641230b98782aef1cd96\_b.jpg][v2-5157db68a32a641230b98782aef1cd96_b.jpg]

图3.3 压缩后的融合图形

![v2-e6215742d43da9cb5d88d39636e18aaa\_b.jpg][v2-e6215742d43da9cb5d88d39636e18aaa_b.jpg]

![v2-b86a6717c42659ba5829bdd7b2687b33\_b.jpg][v2-b86a6717c42659ba5829bdd7b2687b33_b.jpg]

![v2-9da3e224ca87d84d3e4cb8cfd4bb5305\_b.jpg][v2-9da3e224ca87d84d3e4cb8cfd4bb5305_b.jpg]

(a) 未经压缩的提取水印 ( b) 压缩90%后的提取水印 (c) 压缩80%后的提取水印

![v2-b347c8a5b9cb018ad8cfc99ab08b913c\_b.jpg][v2-b347c8a5b9cb018ad8cfc99ab08b913c_b.jpg]

![v2-ca9161bc1d8fb469e4d3dbbd2d55f614\_b.jpg][v2-ca9161bc1d8fb469e4d3dbbd2d55f614_b.jpg]

![v2-80735e8773e0e34ad9cbc36860045d91\_b.jpg][v2-80735e8773e0e34ad9cbc36860045d91_b.jpg]

(d) 压缩70%后的提取水印 (e) 压缩60%后的提取水印 (f) 压缩50%后的提取水印

图3.4 压缩实验后的水印效果图

表3.1 JPEG压缩实验提取与嵌入水印的峰值信噪比和相似度

<table> 
 <tbody> 
  <tr> 
   <td>JPEG压缩比(quality)</td> 
   <td>信噪比(PSNR)</td> 
   <td>相似度(NC)</td> 
  </tr> 
  <tr> 
   <td>50</td> 
   <td>9.4534 db</td> 
   <td>0.9270</td> 
  </tr> 
  <tr> 
   <td>60</td> 
   <td>12.0882 db</td> 
   <td>0.9350</td> 
  </tr> 
  <tr> 
   <td>70</td> 
   <td>14.4313 db</td> 
   <td>0.9588</td> 
  </tr> 
  <tr> 
   <td>80</td> 
   <td>17.6216 db</td> 
   <td>0.9887</td> 
  </tr> 
  <tr> 
   <td>90</td> 
   <td>21.1069 db</td> 
   <td>0.9938</td> 
  </tr> 
 </tbody> 
</table>

由实验可知，小波域水印法对JEPG压缩实验有一定的抗性，压缩比直接影响*PSNR*和*NC*的取值，当Quality取值越大时，压缩质量越好，提取出的水印效果也就越好。

2、噪声实验

加入噪声是对水印鲁棒性考验的一种常见的攻击 。本实验选取两种常用的噪声，高斯噪声和椒盐噪声作为攻击手段。先对融合图像加入高斯噪声(如图6.5)，通过对高斯噪声的方差进行变化，得到不同效果图和数据，如下图6.6和表6.2所示。程序代码如下：

sigma=1;

Y=Y+sigma\*randn(size(Y));

![v2-6b67078da06c8858e7e3d24982b440ad\_b.jpg][v2-6b67078da06c8858e7e3d24982b440ad_b.jpg]

图3.5 加入高斯噪声后的融合图像

![v2-e6215742d43da9cb5d88d39636e18aaa\_b.jpg][v2-e6215742d43da9cb5d88d39636e18aaa_b.jpg]

![v2-ccfa4ac13b1c3e1272b750fba5482311\_b.jpg][v2-ccfa4ac13b1c3e1272b750fba5482311_b.jpg]

![v2-97a517632b3daceb077b8e5c27ff314d\_b.jpg][v2-97a517632b3daceb077b8e5c27ff314d_b.jpg]

(a) 未经处理的提取水印 (b) sigma=1时的提取水印 (c) sigma= 2时的提取水印

![v2-09d5530a1836a1df407e4d02fa86c4e4\_b.jpg][v2-09d5530a1836a1df407e4d02fa86c4e4_b.jpg]

![v2-362108ef6836a3514a974ac0f337d6af\_b.jpg][v2-362108ef6836a3514a974ac0f337d6af_b.jpg]

![v2-9b89f07ae4fb588bf9ce2609ee17f357\_b.jpg][v2-9b89f07ae4fb588bf9ce2609ee17f357_b.jpg]

(d) sigma=3时的提取水印 (e) sigma= 4时的提取水印 (f) sigma= 5时的提取水印

图3.6 高斯噪声实验后的水印效果图

表3.2 高斯噪声实验提取与嵌入水印的峰值信噪比和相似度

<table> 
 <tbody> 
  <tr> 
   <td>高斯噪声的方差(sigma)</td> 
   <td>信噪比(PSNR)</td> 
   <td>相似度(NC)</td> 
  </tr> 
  <tr> 
   <td>1</td> 
   <td>22.030 db</td> 
   <td>0.9963</td> 
  </tr> 
  <tr> 
   <td>2</td> 
   <td>16.1241 db</td> 
   <td>0.9859</td> 
  </tr> 
  <tr> 
   <td>3</td> 
   <td>12.7570 db</td> 
   <td>0.9706</td> 
  </tr> 
  <tr> 
   <td>4</td> 
   <td>9.9961 db</td> 
   <td>0.9453</td> 
  </tr> 
  <tr> 
   <td>5</td> 
   <td>8.1260 db</td> 
   <td>0.9203</td> 
  </tr> 
 </tbody> 
</table>

由上述实验可知，小波域水印法对高斯噪声攻击性实验具有一定的抗性，但效果不太好，高斯方差直接影响*PSNR*和*NC*的取值。

再对融合图像加入椒盐噪声(如图3.7)，通过对椒盐噪声的密度取值进行变化，得到不同的图像和数据，如下图3.8和表3.3所示。程序代码如下：

U=0.5\*ones(size(Y));

J=imnoise(U,'salt & pepper',0.0001);

Y(J==0)=0;

Y(J==1)=255；

![v2-4c0bb83517a170a577b6fbb7ccd36e1c\_b.jpg][v2-4c0bb83517a170a577b6fbb7ccd36e1c_b.jpg]

图3.7 加入椒盐噪声后的融合图像

![v2-e6215742d43da9cb5d88d39636e18aaa\_b.jpg][v2-e6215742d43da9cb5d88d39636e18aaa_b.jpg]

![v2-30bdae09ec68bcc5e37ac23874598eb0\_b.jpg][v2-30bdae09ec68bcc5e37ac23874598eb0_b.jpg]

![v2-b54972b91f337006f602009d9c761a75\_b.jpg][v2-b54972b91f337006f602009d9c761a75_b.jpg]

(a) 未经处理的提取水印 (b)密度为0.0002时的提取水印 (c) 密度为0.0004时的提取水印

![v2-d2f23639845fd8ff6933bed951a00ad9\_b.jpg][v2-d2f23639845fd8ff6933bed951a00ad9_b.jpg]

![v2-e8bb5071722313f21aac83d17ae365e9\_b.jpg][v2-e8bb5071722313f21aac83d17ae365e9_b.jpg]

![v2-30fa61b869fea8f75b34603ddec1e641\_b.jpg][v2-30fa61b869fea8f75b34603ddec1e641_b.jpg]

(d) 密度为0.006时的提取水印 (e) 密度为0.0008时的提取水印 (f)密度为0.001时的提取水印

图3.8 椒盐噪声实验后的水印效果图

表3.3 椒盐噪声实验提取与嵌入水印的峰值信噪比和相似度

<table> 
 <tbody> 
  <tr> 
   <td>椒盐噪声的密度</td> 
   <td>信噪比(PSNR)</td> 
   <td>相似度(NC)</td> 
  </tr> 
  <tr> 
   <td>0.0002</td> 
   <td>16.3944 db</td> 
   <td>0.9939</td> 
  </tr> 
  <tr> 
   <td>0.0004</td> 
   <td>14.7578 db</td> 
   <td>0.9807</td> 
  </tr> 
  <tr> 
   <td>0.0006</td> 
   <td>12.8045 db</td> 
   <td>0.9694</td> 
  </tr> 
  <tr> 
   <td>0.0008</td> 
   <td>11.4493 db</td> 
   <td>0.9588</td> 
  </tr> 
  <tr> 
   <td>0.001</td> 
   <td>10.4967 db</td> 
   <td>0.9451</td> 
  </tr> 
 </tbody> 
</table>

由上述实验可知，小波域水印法对椒盐噪声攻击性实验也具有一定的抗性，但效果不是太好，改变椒盐噪声的密度可直接影响*PSNR*和*NC*的取值。

综上所述，小波域水印法在盲提取和检测中具有良好的鲁棒性和不可见性。在进行抗攻击性实验室，也表现出一定的抗性。但抗压缩性的能力明显优于抗噪声性的能力，以后的改进需向抗噪声干扰这方面进行。

[v2-9ce5ca879497f5f6dec502c2ff617b71_b.jpg]: https://ss.csdn.net/p?https://pic2.zhimg.com/v2-9ce5ca879497f5f6dec502c2ff617b71_b.jpg
[v2-b3ede55b59d8e11c7f668aa97c42e722_b.jpg]: /images/20220828/f0d0b41a91ce40029b3de01b26052426.png
[v2-5157db68a32a641230b98782aef1cd96_b.jpg]: /images/20220828/ab57f07c8e6f421a92206ba3b695e9f8.png
[v2-e6215742d43da9cb5d88d39636e18aaa_b.jpg]: /images/20220828/16701a217b5e4cdda15bcd7abd93f384.png
[v2-b86a6717c42659ba5829bdd7b2687b33_b.jpg]: /images/20220828/95d6760de4c94fdc874432c29290c427.png
[v2-9da3e224ca87d84d3e4cb8cfd4bb5305_b.jpg]: /images/20220828/aa3132b0aa9e4a938b168e67ca3ec947.png
[v2-b347c8a5b9cb018ad8cfc99ab08b913c_b.jpg]: /images/20220828/9b5ebd50431d4a3598ef4da6ab9331ba.png
[v2-ca9161bc1d8fb469e4d3dbbd2d55f614_b.jpg]: /images/20220828/bf21ea16098c4107b43dcb59d8008d24.png
[v2-80735e8773e0e34ad9cbc36860045d91_b.jpg]: /images/20220828/ae60745e35fe48b9ac525c9d07001394.png
[v2-6b67078da06c8858e7e3d24982b440ad_b.jpg]: /images/20220828/4cfbce78b49d4790911ce49788e7280d.png
[v2-ccfa4ac13b1c3e1272b750fba5482311_b.jpg]: /images/20220828/3af86b09e148442a9f656966a0fa721d.png
[v2-97a517632b3daceb077b8e5c27ff314d_b.jpg]: /images/20220828/58be0b465e214a02879026f3dbcf03f9.png
[v2-09d5530a1836a1df407e4d02fa86c4e4_b.jpg]: /images/20220828/9e3b7ce67b954573965d90ef497b9eaa.png
[v2-362108ef6836a3514a974ac0f337d6af_b.jpg]: /images/20220828/0cbdc7c0289144d09b582200026ce5ea.png
[v2-9b89f07ae4fb588bf9ce2609ee17f357_b.jpg]: /images/20220828/640c06cd26204a4ba57d525430cbbc4a.png
[v2-4c0bb83517a170a577b6fbb7ccd36e1c_b.jpg]: /images/20220828/e784f93d86024411a693c380d5604dfd.png
[v2-30bdae09ec68bcc5e37ac23874598eb0_b.jpg]: /images/20220828/d5e453a8569b461c8d6bc115f9f452f5.png
[v2-b54972b91f337006f602009d9c761a75_b.jpg]: /images/20220828/65523c8f912b43f3b5b5ff47908d34c5.png
[v2-d2f23639845fd8ff6933bed951a00ad9_b.jpg]: /images/20220828/3655edc14260473598b4c3a02230db18.png
[v2-e8bb5071722313f21aac83d17ae365e9_b.jpg]: /images/20220828/6bf27d762ecb4e989e7b0c5af08f65ce.png
[v2-30fa61b869fea8f75b34603ddec1e641_b.jpg]: /images/20220828/c9396594668a4f0fa8525f0c1b4f389c.png