专利名称：一种针对时域音频lsb隐写的检测方法
技术领域：
本发明涉及信息隐藏检测领域，特别涉及一种针对时域音频LSB隐写的检 测方法。
背景技术：
近年来，隐写术和数字水印技术的研究得到较大的发展。许多隐秘术和数 字水印软件能够从网上直接下载，使普通人利用这些技术实现"隐蔽"通讯成 为可能。这种现象客观上要求隐写检测技术的进步，以便检测和阻止不法信息 隐藏在看似无异动普通载体(图像、音频和视频等)中传递。
在众多的音频信息隐藏方法中，空域LSB(最低比特位)隐藏以其隐蔽性较 好、信息隐藏量大且易于实现等优点而被广泛应用。因此，对LSB信息隐藏进 行有效、准确而又可靠的检测，对于保护信息安全、防止隐秘通信有着十分重 要的意义。
目前针对LSB隐写算法的检测方法均针对图像载体，典型的有Chi-Square 方法，Fridrich提出的RS检测方法，Dumitrescu等提出的样本对方法(SPA)。 其中Chi-Square方法是把隐秘对象的理论期望频率分布和从可能被修改的载体 中观察到样本分布进行比较，从而找出差异，检测是否有信息嵌入。对音频信 号的检测虚警率很高，原因是最后一个平面的值几乎是随机的，比如Windows 注销音。RS检测方法考虑到起点是图像各个位平面之间具有一定的非线性相关 性，而当LSB隐藏信息后，这种相关性将被破坏。但是，由于载体音频信号本身可能具有一定的随机性，通过RS分析可能会得出含有少量隐藏密文的结论， 使其虚警率偏高。音频信号含有的噪声程度也影响对隐密信息的检测，音频和
其LSB具有相当大的随机性，导致计算准确度下降。另外，消息位在音频中的 位置也是影响估计精度的重要因素。SPA方法基于抽样对及其有限状态机分析， 本质是基于随机性特征的高阶统计分析。在音频信号中，相邻点的样本值相差 很大，绝大多数载体在隐藏操作前后相邻点的差值大于3，使算法失效。可见， 本领域缺乏针对时域音频LSB隐写的高效检测技术手段。
发明内容
本发明目的在于解决现有技术不足，提供一种针对时域音频LSB隐写的检 测方法。
为实现上述目的，本发明提供一种针对时域音频LSB隐写的检测方法，具 体步骤如下，
歩骤一，读取待测音频文件，将文件内各样点的样点值整数化；
步骤二，通过考察待测音频文件中样本对样点值的变化，得出样本对统计变化
规律，
首先，将含有^个样点的待测音频文件分为^-1个样本对(1,^+,)，其中 0</<W， x,， x,+i为相邻样点；样点x,的幅值表示为2" + L其中* = 0，1;样点;c,+， 的幅值表示为2m + y，其中_/ = 0,1;
然后，考察上述所得样本对样点值变化，包括以下三部分内容， 考察待测音频文件附= 时样本对(、,、1)的四种情况，用6。表示"0,) = 0时样 本对(x,，;c,+,;)的集合，统计该集合中元素的个数，用IW表示；用6,表示b(U-l时 样本对(x,，x,J的集合，统计该集合中元素的个数，用hl表示；用62表示"1,) = 0时样本对b,x,J的集合，统计该集合中元素的个数，用|62|表示；用63表示 "U、l时样本对(x,，x,+,)的集合；统计该集合中元素的个数，用l^表示；
考察待测音频文件/^m + l时样本对(x,，;c,+,)的情况，用《表示"0,7 = 1时样本 对0c,，、+,)的集合，统计该集合中元素的个数，用l刈表示；
考察待测音频文件"=m-1时样本对(u,+,)的情况，用《表示"l，y = 0时样本 对(x,，x,+,)的集合，统计该集合中元素的个数，用^| 步骤三，通过样本对统计规律分析嵌入强度，
首先，将求得的ki，
62|， |63|， |《|， IW值分别带入下面32个等式;
& ,。l + l63l + 2l^卜41
t =4 6, -4k
4k
A:5 =|60| + |63| + 2
A:7 =|60| + |63| + 2|62|_4|《 & H6ol + H + 2l&l-41《
一
/i13 —1"0
A,， = |6r
^23 = 16。
&_ I/)
""29 — |"0
+ |63 +2
+ |
+ |63| + 2^2|-4|<i, + |63| + 2|62|-4|《 + |63| + 2|62|_4|《 + 1631 + 21^1-4^2 + 1631 + 2^1-4^2 + |&| + 2k|-4k
+ W + 2k -4
& =,1-41^
A。 =4H-4|<i2 A;12 =來|-4|《
夂6 "l^l-^ 夂s =來|-4|<^
A:20 =4p2
4 k
&2 =,|-4|<i2 =來|-4|<^
*26 =4|62|-4|(i2
A28 =4|^|_丰2
&。=來|_4|《
=|6o| + |63| + 2|62|-4|<i2|， A"=來|-4|c/2| 然后，将上述求得的32个值分别代入下面等式，求出16个嵌入强度估计值 ar，其中x二l…16;
= 1 _ 、/l -々.
2W 2_r-
最后，取16个嵌入强度估计值^中的最小值作为嵌入强度分析值，用S表
6不；
a = min, <z2, or3, a4, 5 ， "6, "7, 8 ， "9 ，0, an, 2, 3, a, 4, a, 5, 6) 步骤四，判断嵌入强度分析值S是否达到强度阈值，若嵌入强度分析值^大于或 等于强度阀值则检测得出待测音频文件中含有LSB隐秘信息，否则检测结果为 不含LSB隐秘信息。
而且，强度阀值设为0.05。
本发明提高了一种针对音频LSB隐写的主动分析检测方法，无需先验条件， 是一种"盲"的检测方案。本发明提供的方案不但能检测音频信号中是否含有 LSB隐写嵌入的隐秘信息，而且因为是从估计隐密信息的嵌入强度入手，对后续 破解信息具有实用意义。采用本发明进行隐写检测具有高效性，并且正确率高， 虚警率和漏检率低。本发明适用于对多种类型的音频文件测试，性能稳健；同
时算法复杂性低、计算开销小。可见这种隐写分析检测方法在我国信息安全建
设中将可广泛发挥重要作用。
图1是本发明流程图。
具体实施例方式
本发明提供一种针对音频LSB隐写的检测方法，通过考察LSB隐写带来的 音频特定样本对样点值的变化，得出样本对统计变化规律，最后推导出嵌入强 度的二次方程，预估嵌入强度，并通过强度阀值判断待测音频文件包含的音频 信息中是否存在LSB隐写嵌入的隐秘信息。具体步骤为 步骤一，读取待测音频文件，将文件内各样点的样点值整数化；
所述待测音频文件一般为wav格式，若是压縮格式，则应当解码出完整的音频信息后进行处理。从待测音频文件内选取样点的具体方式属于现有成熟技 术，本发明不予赘述。
步骤二，通过考察待测音频文件中样本对的变化，得出样本对统计变化规律，
首先，将含有W个样点的待测音频文件分为W-l个样本对(x,，;c,+，;)，其中 0d<AA， x,， x,+,为相邻样点；从待测音频文件中抽取样点的数目根据具体实施 时的精度要求而定，即7V为整数，W值大小可以根据精度要求选定。样点值整 数化后得到样点的幅值，考察样点的幅值即是考察该样点的样点值。整数可能 为奇数或偶数情况，因此本发明将样点x,的幅值表示为2" + h其中& = 0，1;样点 x,+,的幅值表示为+./ ，其中y' = 0， 1 。
然后，考察上述所得样本对样点值变化，包括以下三部分内容。这三种考 察统计能够精练而准确地反应LSB隐写带来的样本对统计变化规律，只需统计 出6个数值，因此有利于本发明技术方案的高效性，特别是具有实施时漏检率 低的优点。
(1)考察待测音频文件m = n时样本对(x,,)的四种情况，用6。表示"0,= 0
时样本对(X,，;C,J的集合，统计该集合中元素的个数，用l6。l表示；用6,表示"(^、1
时样本对(x,，x,J的集合，统计该集合中元素的个数，用IW表示；用62表示
"1，/=0时样本对(^《+1;)的集合，统计该集合中元素的个数，用iw表示；用63
表示^ = 1,/ = 1时样本对(1,^+1)的集合；统计该集合中元素的个数，用|63|表示； (2 )考察待测音频文件"=附+1时样本对(x,, ;c,+,)的情况，用《表示"0,7' = 1时
样本对(x,，x,+,)的集合，统计该集合中元素的个数，用l刈表示；
(3 )考察待测音频文件"=m _ 1时样本对(x,,x,+l)的情况，用《表示"W = 0时
样本对(x,，x,+,)的集合，统计该集合中元素的个数，用l叫表示；步骤三，通过统计规律分析嵌入强度， 首先，将求得的|6。
1，1~卜W ,|<5 ,| ，值分别带入下面
+ |63|
=+ 2W_4K|:4W-
&=+ 2W-來」:4|、|--4《|
= |-4|t/2|，&8 =--4W
=60|+ 1^1+丰|-4W，=4W--來
&u=|6。+丰|， 夂2=來|
二|~+丰_4W，=來|-4 t/2|
夂5=|6。l + lh+ 2 62，=,-4|《|
辟2-4《，^18=4 62|_4W
=+ 2 62-4《， 0=來
=+丰1-4|《，、2=剩-4|t/2|
、3=|6。+丰-4|《，、4=4 62-4W
允25=|6。+丰-4|<i2，&26-4W
克27|&l + lh+ 2 62_4《t 8_4|"2|
=|6。+ 2 62_4《， 0-4W
|60+丰2-4|《，*32=來-來l
然后，将上述求得的32个值分别代入下面等式，求出16个嵌入强度估计值
，简称为16个《值，其中^ = 1...16;
"v = 1 —、/l —:
以上算法实际上来自根据样本对统计变化规律推导得出的关于嵌入强度的
欠方程
(I 6。 | + |63 |+2|Z^—4|^|)a2—(2|6。 I+2163 |+41 6'卜81 J |) + 416'卜41 d' |=0 求解其根如下
' (x = 1...16)， A;2l—,=卜0| + |63| + 2|6'|-卞|， 、 =4|6卜卞|
a = 1 — 、/l - k
/
其中M为W或IW， l刈为W或KI。考虑全部组合，即共得到16组、—简 称为32个H直。这种嵌入强度的二次方程具有通用性，具体实施时可以直接按
9照二次方程的求解根式，将步骤三样本对统计得到的值带入本发明技术方案提 供的公式进行计算即可。因为算法简单，降低了计算开销，进一步保证了高速 高效，降低了虚警率。
最后，取16个嵌入强度估计值^中的最小值作为嵌入强度分析值，用S表
示。选取最小值是为了最大限度防止误警，具体实施时也可考虑采用平均值。 "=min(",，"2,a:3，"4,"5,"6，6ir7,"8,"9，6^o,"u，a:,2,"B,",4,c^5,",6)
步骤四，判断嵌入强度分析值^是否达到强度阈值，若嵌入强度分析值S大于或
等于强度阀值则检测得出待测音频文件内的音频信息中含有LSB隐秘信息，否 则检测结果为不含LSB隐秘信息。
强度阀值可以预先设定，将判断过程通过软件方式自动完成。或者具体实 施时根据检测灵敏度要求设定。本发明人广泛采用含有LSB隐秘信息的音频文 件进行实验，得出一般含有LSB隐秘信息的音频文件嵌入强度大于或等于0. 05 的结论。因此，本发明提出最优选实施方案应当将强度阔值设定在0.05。参见 图l，本发明的整体技术方案可以采用软件方式实现，自动进行检测的基本流程 可以简单总结为打开待测音频文件；将样点值整数化；统计得到hl， W， |62|，
|63|， |《|， |^| ;计算32个A值；计算16个"值；通过16个《值的最小值求得^;
根据判断条件^w.Q5，分析出是否不含LSB隐秘信息。
权利要求
1.一种针对时域音频LSB隐写的检测方法，其特征是具体步骤如下，步骤一，读取待测音频文件，将文件内各样点的样点值整数化；步骤二，通过考察待测音频文件中样本对样点值的变化，得出样本对统计变化规律，首先，将含有N个样点的待测音频文件分为N-1个样本对(xi，xi+1)，其中0＜i＜N，xi，xi+1为相邻样点；样点xi的幅值表示为2n+k，其中k＝0，1；样点xi+1的幅值表示为2m+j，其中j＝0，1；然后，考察上述所得样本对样点值变化，包括以下三部分内容，考察待测音频文件m＝n时样本对(xi，xi+1)的四种情况，用b0表示k＝0，j＝0时样本对(xi，xi-1)的集合，统计该集合中元素的个数，用|b0|表示；用b1表示k＝0，j＝1时样本对(xi，xi+1)的集合，统计该集合中元素的个数，用|b1|表示；用b2表示k＝1，j＝0时样本对(xi，xi+1)的集合，统计该集合中元素的个数，用|b2|表示；用b3表示k＝1，j＝1时样本对(xi，xi+1)的集合；统计该集合中元素的个数，用|b3|表示；考察待测音频文件n＝m+1时样本对(xi，xi+1)的情况，用d1表示k＝0，j＝1时样本对(xi，xi+1)的集合，统计该集合中元素的个数，用|d1|表示；考察待测音频文件n＝m-1时样本对(xi，xi+1)的情况，用d2表示k＝1，j＝0时样本对(xi，xi+1)的集合，统计该集合中元素的个数，用|d2|表示；步骤三，通过样本对统计规律分析嵌入强度，首先，将求得的|b0|，|b1|，|b2|，|b3|，|d1|，|d2|值分别带入下面32个等式；k1＝|b0|+|b3|+2|b1|-4|d1|，k2＝4|b1|-4|d1|k3＝|b0|+|b3|+2|b2|-4|d1|，k4＝4|b1|-4|d1|k5＝|b0|+|b3|+2|b1|-4|d2|，k6＝4|b1|-4|d1|k7＝|b0|+|b3|+2|b2|-4|d2|，k8＝4|b1|-4|d1|k9＝|b0|+|b3|+2|b1|-4|d1|，k10＝4|b1|-4|d2|k11＝|b0|+|b3|+2|b1|-4|d1|，k12＝4|b2|-4|d1|k13＝|b0|+|b3|+2|b1|-4|d1|，k14＝4|b2|-4|d2|k15＝|b0|+|b3|+2|b2|-4|d1|，k16＝4|b1|-4|d2|k17＝|b0|+|b3|+2|b2|-4|d1|，k18＝4|b2|-4|d1|k19＝|b0|+|b3|+2|b2|-4|d1|，k20＝4|b2|-4|d2|k21＝|b0|+|b3|+2|b1|-4|d2|，k22＝4|b1|-4|d2|k23＝|b0|+|b3|+2|b1|-4|d2|，k24＝4|b2|-4|d1|k25＝|b0|+|b3|+2|b1|-4|d2|，k26＝4|b2|-4|d2|k27＝|b0|+|b3|+2|b2|-4|d2|，k28＝4|b1|-4|d2|k29＝|b0|+|b3|+2|b2|-4|d2|，k30＝4|b2|-4|d1|k31＝|b0|+|b3|+2|b2|-4|d2|，k32＝4|b2|-4|d2|然后，将上述求得的32个值分别代入下面等式，求出16个嵌入强度估计值αx，其中x＝1...16；<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mi>α</mi> <mi>x</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><msqrt> <mn>1</mn> <mo>-</mo> <msub><mi>k</mi><mrow> <mn>2</mn> <mi>x</mi></mrow> </msub> <mo>/</mo> <msub><mi>k</mi><mrow> <mn>2</mn> <mi>x</mi> <mo>-</mo> <mn>1</mn></mrow> </msub></msqrt> </mrow>]]></math></maths>最后，取16个嵌入强度估计值αx中的最小值作为嵌入强度分析值，用 id="icf0002" file="A2008100467570003C2.tif" wi="2" he="4" top= "180" left = "180" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>表示；<maths id="math0002" num="0002" ><math><![CDATA[ <mrow><mover> <mi>α</mi> <mo>^</mo></mover><mo>=</mo><mi>min</mi><mrow> <mo>(</mo> <msub><mi>α</mi><mn>1</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>2</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>3</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>4</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>5</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>6</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>7</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>8</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>9</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>10</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>11</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>12</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>13</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>14</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>15</mn> </msub> <mo>,</mo> <msub><mi>α</mi><mn>16</mn> </msub> <mo>)</mo></mrow> </mrow>]]></math></maths>步骤四，判断嵌入强度分析值 id="icf0004" file="A2008100467570003C4.tif" wi="2" he="4" top= "212" left = "82" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>是否达到强度阈值，若嵌入强度分析值 id="icf0005" file="A2008100467570003C5.tif" wi="2" he="4" top= "213" left = "169" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>大于或等于强度阀值则检测得出待测音频文件中含有LSB隐秘信息，否则检测结果为不含LSB隐秘信息。
2.根据权利要求1所述的隐写检测方法，其特征是强度阀值设为0.05。
全文摘要
本发明公开了一种针对时域音频LSB隐写的检测方法，具体步骤如下读取待测音频文件，将文件内各样点的样点值整数化；通过考察待测音频文件中样本对样点值的变化，得出样本对统计变化规律；通过样本对统计规律分析嵌入强度；判断嵌入强度分析值a是否达到强度阈值，若嵌入强度分析值a大于或等于强度阀值则检测得出待测音频文件中含有LSB隐秘信息，否则检测结果为不含LSB隐秘信息。该方法不但能检测出LSB隐写信息的存在，还能估算出嵌入强度，保证了检测的高正确率以及低虚警率和低漏检率，而且该方法使用简单，时间复杂度低。本发明可应用于隐秘信息检测领域，对防止非法的隐蔽通信具有重要的意义。
文档编号G10L19/00GK101494051SQ20081004675
公开日2009年7月29日 申请日期2008年1月23日 优先权日2008年1月23日
发明者伟 曾, 胡瑞敏, 艾浩军, 尚 高 申请人:武汉大学