信息隐蔽 试验十一 LSB信息隐蔽的RS分析

本文格式为Word版，下载可任意编辑——信息隐蔽试验十一LSB信息隐蔽的RS分析LSB信息隐蔽的RS分析

一，

试验目的：

了解RS隐写分析的原理，把握一种图像LSB隐写算法的分析方法，设计并实现一种基于图像的LSB隐写的RS隐写分析算法。二，

试验环境

（1）Windows7操作系统；（2）MATLABR2023b版本软件；（3）图像文件lena.bmp；（4）S-Tools工具；三，

试验原理

1.RS隐写分析原理

RS主要是针对采用伪随机LSB嵌入算法进行攻击的一种方法。RS方法不但能检测出图像是否隐蔽信息，而且还能比较确凿地估算出隐蔽的信息长度。

RS隐写分析算法考虑图像各个位平面之间具有一定的非线性相关性，当利用LSB隐写算法隐蔽机要信息后，这种相关性就会破坏。只要能找出衡量这一相关性的方法，并对隐蔽机要信息前后的状况加以对比，就有可能设计出隐写分析方法。

RS隐写分析方法的理论核心是：任何经过LSB隐写的图像，其最低比特位分布满足随机性，即0、1的取值概率均为1/2，而未经过隐写的图像不存在此特性。对于一个M?N像素的图片，设各个像素的值取自集合P，例如一个8bit的灰度图像，P?{0,1,2,?,255}。将这些像素分为有着n个相邻像素的子集，例如n可以取值为4，记为G?(x1,x2,x3,x4)。进一步利用如下函数表示图像块的空间相关性，即

f(x1,x2,x3,?,xn)??xi?1?xii?1n?1（1）

函数f使得每一个集合G都对应一个实数。G中的噪声越大，函数f的值越大。f的值越小，说明图像相邻像素之间的起伏越小，而图像块的空间相关性越强。然后，定义集合P上的3个函数：

（1）交换函数F1：2n?2n?1,0?1,2?3,?,254?255。（2）偏移函数F?1：2n?2n?1,?1?0,1?2,?,255?256。（3）恒等变换F0：F0(x)?x,x?P。

繁杂度分辩函数f和三种变换函数F定义以下三种类型的元素组R,S,U。（1）RegularG：G?R?f(F(G))?f(G)，也就是说对G中的元素进行变换之后增大了元素之间的区别程度，R表示正则组。

（2）SingularG：G?S?f(F(G))?f(G)，也就是说对G中的元素进行变换之后减小了元素之间的区别程度，S表示奇异组。

（3）UnusableG：G?U?f(F(G))?f(G)，也就是对G中的元素进行变

换之后元素之间的区别程度不变，G表示无用组。

引入伪装M，它的取值范围是-1、0、1。对应着F?1、F0、F1三个变换函数。记对于伪装M的R组的数量为RM，类似的S组记为SM。例如G=（27,28,26,25），M=（1,0,1,0），则-M=（-1,0，-1,0）。于是就得到FM(G)=（26,28,27,25），F?M(G)=（28,28,25,25）。

对于一个原始不包含隐蔽信息的载体来说具有下面的规律：RM?R?M?SM?S?M?S?M?RM?SM（2）

但是当在载体中嵌入了机要信息，就会有下面的式子成立：R?M（3）

也就是说对于一个待检测的载体，只需要得到上述的统计量，然后通过对比这四个统计参数之间的大小关系就能够判断出载体中是否含有机要信息。

RS算法不但能够判断载体是否含有机要信息，RS隐写分析算法还能够估计出隐写容量。上面算出的一组RM、SM、R?M、S?M是在比例为p/2（p为隐写率）的样点值已经作了一次F1变换状况下得到的。将该隐写体的所有样点值应用

F1变换，相对于原始载体则有比例为（1-p/2）的样点值作了F1变换，计算此时的RM、SM、R?M、S?M值。如下图1所示，得到的是两组统计量的数据。

图1RS曲线与像素更改比率间的关系示意图

具体步骤如下：

将待检测图像分成若干大小相等的图像块，寻常定义n=4，对每个小图像块定义一个掩膜算法函数，寻常为M（0,1,1,0），即判断此图像小块属于Regular或Singular，FM(G)?(F0(x1),F1(x2),F1(x3),F0(x4))，

遍历整幅图像，得到RM(p/2)，SM(p/2)。RM(p/2)，SM(p/2)定义为：RM(p/2)?RM（4）

bufsize/4SMSM(p/2)?（5）

bufsize/4改变掩膜算法函数，令M为（0，-1，-1,0）重复上一步操作得到R?M(p/2)和S?M(p/2)。

将图像数据的lSB位取反，用获得的数据重复上述操作得到RM(1?p/2)，

SM(1?p/2)以及R?M(1?p/2)，S?M(1?p/2)。

结合如下方程：

2(d1?d0)x其中：

2?d?0?d?1?3d0x?d0?d?1?0（6）

d0?RM(p/2)?SM(p/2),d1?RM(1?p/2)?SM(1?p/2)d?0?R?M(p/2)?S?M(p/2),d?1?R?M(1?p/2)?S?M(1?p/2)

解方程式（6），取绝对值较小的x，计算嵌入概率p为：p?x/(x?1/2)。2.S-Tools工具简介

S-Tools是一个时（空）域数字水印软件，支持WAV格式的音频文件、GIF和BMP格式的图像文件，下面介绍基于音频文件的隐蔽原理和使用方法。

WAV格式的音频文件在Windows中存储为8位或16位值，对于8位样本，这意味着取值范围介于0和255之间。16位的样本取值范围介于0和65535之间。S-Tools在WAV格式的音频文件中隐蔽信息时，首先用密码加密，然后根据密码决定要嵌入信息的位置，采用LSB（leastsignificantbits）方法，把机要信息代替载体音频文件的最不重要比特位。

例如，假设一个音频文件有以下8个字节的信息，分别为1321341371411211017438二进制表示为

10000100100001101000100110001101011001010100101000100110

假使要隐蔽二进制字节11010101（213）在这个序列里。简单地取代了音频文件的LSB（最低有效位）进行信息隐蔽。因此上述顺序将变更为

1331351361411201017439用二进制表示为

10000101100001111000100010001101011110000100101000100111

这样，机要信息就隐蔽在载体音频文件中，音频文件在听觉效果上和原始文件几乎没有区别。

四，试验分析1.隐蔽信息并分析

源代码RS.m如下：

clc,clearall

org=imread('lena.bmp');org=double(org);org=org(:)';

[rowsample_len]=size(org);rs=zeros(2,4);cor=zeros(1,3);m=floor(sample_len/4);M=randsrc(4,1,[01]);tmp=zeros(4,1);

%使用LSB，准备待检测图像，隐写率为40%msg_len=floor(sample_len*0.4);msg=randsrc(msg_len,1,[01;0.50.5]);stg=org;fori=1:msg_len

stg(i)=bitset(stg(i),1,msg(i));end

%对待检测图像应用非负和非正翻转，统计翻转前后的空间相关性i=1;forj=1:m

tmp=stg((j-1)*4+1:j*4);%把图像分为四个像素一组cor(1)=SpaceCor(tmp);cor(2)=SpaceCor(fpos(tmp,M));cor(3)=SpaceCor(fneg(tmp