抗任意剪切和翻录的音频水印方法及系统与流程

1.本发明涉及音频水印技术领域，具体涉及一种抗任意剪切和翻录的音频水印方法及系统。


背景技术：

2.随着互联网覆盖率程度的不断加深以及人们阅读习惯的逐渐改变，作为传播知识重要载体的有声读物也越来越受到市场的关注，并被越来越多的读者所接受。然而随着市场规模的不断扩大，不法分子开始使用盗版音频来牟取非法利益，损害了有声读物、音乐等音频所有者的利益。因此，如何保护有声读物、音乐等音频所有者的版权已成为亟待解决的问题。
3.目前，在音频水印技术领域，有基于稳健特征点的平稳小波域数字水印算法，也有基于离散余弦系数均值特征的水印方案,等等，这些音频水印技术能够在音频完整、遭受较小攻击的情况下实现水印的嵌入和提取，一定程度上满足了技术人员的需要。
4.然而，这些现有的音频水印技术难以抵抗对音频任意剪切的攻击，特别是剪切长度过多的情况，同时也难以抵抗翻录软件和录音笔 line
‑
in模式的翻录攻击。由此可见，现有技术无法在音频被极大破坏的情况下完整恢复水印。


技术实现要素：

5.(一)解决的技术问题
6.针对现有技术的不足，本发明提供了一种抗任意剪切和翻录的音频水印方法及系统，解决了现有技术无法在音频被极大破坏的情况下完整恢复水印的问题。
7.(二)技术方案
8.为实现以上目的，本发明通过以下技术方案予以实现：
9.第一方面，本发明提出了一种抗任意剪切和翻录的音频水印方法，所述方法包括：
10.提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改以获取特征值序列，并基于所述特征值序列生成对应的量化特征序列；
11.基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频；
12.基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频。
13.优选的，所述提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改以获取特征值序列包括：
14.对待嵌入水印的音频信号进行分段，并对每一段依次进行三级离散小波变换和离散余弦变换，并将获得系数进行奇异值分解，从而获得音频的变换域特征；
15.对音频的变换域特征进行修改以获取特征值序列。
16.优选的，所述基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频包括：
17.选择本原多项式并生成对应的m序列，然后将m序列与音频的所述量化特征序列进行异或生成密钥1，并将一维的水印序列与m序列进行异或生成密钥2。
18.优选的，所述基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频时，基于启发式搜索算法对参数进行寻优，以不断调整水印嵌入时的嵌入强度。
19.优选的，所述方法还包括：采用同嵌入水印时相同的量化规则对水印进行提取，其具体过程包括：
20.s41、采用改进的谱熵法对已嵌入水印的音频进行时域上的恢复；
21.s42、对已嵌入水印的音频进行分段，并对每一段依次进行三级离散小波变换，离散余弦变换，奇异值分解，然后采用嵌入水印时相同的规则对变换域特征进行量化以获取量化特征序列；
22.s43、将现有音频的特征与所述量化特征序列进行对比，然后结合密钥1和密钥2提取水印图像。
23.优选的，所述方法还包括：
24.s44、当利用s41
‑
s43步骤提取水印失败后，采用模式匹配
‑
回滚法对已嵌入水印的音频进行时域上的恢复，然后继续进行s42
‑
s43步骤对水印进行提取。
25.第二方面，本发明还提出了一种抗任意剪切和翻录的音频水印系统，所述系统包括：
26.嵌入预处理模块，用于提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改以获取特征值序列，并基于所述特征值序列生成对应的量化特征序列；
27.水印嵌入模块，用于基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频；
28.水印绑定模块，用于基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频。
29.优选的，所述嵌入预处理模块提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改以获取特征值序列包括：
30.对待嵌入水印的音频信号进行分段，并对每一段依次进行三级离散小波变换和离散余弦变换，并将获得系数进行奇异值分解，从而获得音频的变换域特征；
31.对音频的变换域特征进行修改以获取特征值序列。
32.优选的，所述水印绑定模块基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频包括：
33.选择本原多项式并生成对应的m序列，然后将m序列与音频的所述量化特征序列进行异或生成密钥1，并将一维的水印序列与m序列进行异或生成密钥2。
34.优选的，所述系统还包括水印提取模块：采用同嵌入水印时相同的量化规则对水印进行提取，其具体过程包括：
35.s41、采用改进的谱熵法对已嵌入水印的音频进行时域上的恢复；
36.s42、对已嵌入水印的音频进行分段，并对每一段依次进行三级离散小波变换，离散余弦变换，奇异值分解，然后采用嵌入水印时相同的规则对变换域特征进行量化以获取量化特征序列；
37.s43、将现有音频的特征与所述量化特征序列进行对比，然后结合密钥1和密钥2提取水印图像；
38.s44、当利用s41
‑
s43步骤提取水印失败后，采用模式匹配
‑
回滚法对已嵌入水印的音频进行时域上的恢复，然后进行s42
‑
s43步骤以提取水印。
39.(三)有益效果
40.本发明提供了一种抗任意剪切和翻录的音频水印方法及系统。与现有技术相比，具备以下有益效果：
41.本发明的技术方案针对待嵌入水印的音频，首先提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改和量化以获取特征值序列和量化特征序列，基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频，然后基于m序列将水印与所述量化特征序列进行绑定。该技术方案将m序列与dwt
‑
dct
‑
svd域的水印技术相结合，通过该法嵌入至音频的水印可在音频遭受极大破坏的情况下仍能完整性恢复。
附图说明
42.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
43.图1为本发明抗任意剪切和翻录的音频水印嵌入过程的流程图；
44.图2为本发明中音频特征提取与量化过程流程图；
45.图3为本发明中水印与音频特征绑定流程图；
46.图4为本发明中水印嵌入流程图；
47.图5为本发明中参数λ优化流程图；
48.图6为本发明中水印提取流程图；
49.图7为本发明中对剪切后音频进行时域恢复流程图；
50.图8为使用本法的水印鲁棒性测试结果图；
51.图9为使用方法一的水印鲁棒性测试结果图；
52.图10为使用方法二的水印鲁棒性测试结果图；
53.图11为本发明实施例检测水印不可感知性时的音频信号波形图。
具体实施方式
54.为使本发明实施例的目的、技术方案和优点更加清楚，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
55.本技术实施例通过提供一种抗任意剪切和翻录的音频水印方法及系统，解决了现有技术无法在音频被极大破坏的情况下完整恢复水印的问题，实现了在音频经过任意剪切和翻录后，仍然可以正确的提取出水印，以保护音频作品版权的目的。
56.本技术实施例中的技术方案为解决上述技术问题，总体思路如下：
57.本发明的技术方案为了实现音频在被极大破坏的情况下仍能完整恢复水印的目的，在音频嵌入水印时，首先，提取待嵌入水印的音频特征并进行修改和量化以获取特征值序列和量化特征序列，基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频，然后基于m序列将水印与所述量化特征序列进行绑定。而在提取音频水印时，根据音频水印被攻击的具体情况，选用短时改进子带谱熵法与双门限法或者使用模式匹配
‑
回滚法，以实现对遭受大面积攻击后的音频的时域特征进行恢复，从而最终实现水印的提取，可提高水印的鲁棒性。
58.为了更好的理解上述技术方案，下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
59.实施例1：
60.本发明首先提出了一种抗任意剪切和翻录的音频水印方法，所述方法包括：
61.s1、提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改以获取特征值序列，并基于所述特征值序列生成对应的量化特征序列；
62.s2、基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频；
63.s3、基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频。
64.由此可见，本发明的技术方案针对待嵌入水印的音频，首先提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改和量化以获取特征值序列和量化特征序列，基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频，然后基于m序列将水印与所述量化特征序列进行绑定。该技术方案将m序列与 dwt
‑
dct
‑
svd域的水印技术相结合，通过该法嵌入至音频的水印可在音频遭受极大破坏的情况下仍能完整性恢复。
65.在本发明实施例的实现过程当中，为了提取待嵌入水印的音频在 dwt
‑
dct
‑
svd这个变换域的特征，同时对相应的特征进行修改与量化，使它能够嵌入水印信息，一种较佳的处理方式是，提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改以获取特征值序列时包括：
66.对待嵌入水印的音频信号进行分段，并对每一段依次进行三级离散小波变换和离散余弦变换，并将获得系数进行奇异值分解，从而获得音频的变换域特征；
67.对音频的变换域特征进行修改以获取特征值序列。
68.实际上，为了根据很少的连续音频特征完整地恢复出整个水印，从而有效地保护音频与视频发行人的版权，本发明实施例使用了m序列与本原多项式，而在基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频包括：
69.选择本原多项式并生成对应的m序列，然后将m序列与音频的所述量化特征序列进行异或生成密钥1，并将一维的水印序列与m序列进行异或生成密钥2。
70.另外，为了在嵌入水印时选择一个合适的水印嵌入强度，一种较佳的处理方式是，在基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频时，基于启发式搜索算法对参数进行寻优，以不断调整水印嵌入时的嵌入强度。
71.在本发明实施例的实现过程当中，为了在音频遭受大规模剪切和翻录时仍能很方便的提取出水印，一种较佳的处理方式是，本发明的抗任意剪切和翻录的音频水印方法还
包括：采用同嵌入水印时相同的量化规则对水印进行提取，其具体过程包括：
72.s41、采用改进的谱熵法对已嵌入水印的音频进行时域上的恢复；
73.s42、对已嵌入水印的音频进行分段，并对每一段依次进行三级离散小波变换，离散余弦变换，奇异值分解，然后采用嵌入水印时相同的规则对变换域特征进行修改与量化以获取量化特征序列；
74.s43、将现有音频的特征与所述量化特征序列进行对比，然后结合密钥1和密钥2提取水印图像。
75.另外，当利用s41
‑
s43步骤提取水印失败后，执行s44、采用模式匹配
‑
回滚法对已嵌入水印的音频进行时域上的恢复，然后继续进行 s42
‑
s43步骤对水印进行提取。在提取音频水印时利用模式匹配
‑
回滚法，可实现对遭受大面积剪切攻击后的音频的时域特征的恢复，从而间接地提高了水印的鲁棒性。
76.下面结合对s1
‑
s3具体步骤的解释，来详细说明本发明的具体实现过程。
77.s1、提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改以获取特征值序列，并基于所述特征值序列生成对应的量化特征序列。
78.提取待嵌入水印的音频在dwt
‑
dct
‑
svd这个变换域的特征，同时对相应的特征进行修改与量化以分别获取特征值序列和量化特征序列，使音频能够嵌入水印信息。参见图2，音频特征提取与量化的过程具体为：
79.首先，对待嵌入水印的音频信号进行分段，并对每一段音频信号依次进行离散小波变换(discrete wavelet transformation，dwt)和离散余弦变换(discrete cosine transform，dct)，并将获得的系数进行奇异值分解(singular value decomposition，svd)，从而获得相应的变换域的特征。例如，以每10000个样本点为一段对音频进行分段，然后采用db1小波将每一段音频进行三级dwt变换，再进行dct变换，并对其进行奇异值分解后，然后将奇异值保存在列表feature中。
80.将feature中的每n个特征分为一组，将分得的每一组分别记为f (1)，f(2)，
…
，f(x)，最后会得到x组【x＝length(feature)
‑ꢀ
mod(length(feature),2n)】，其中，mod表示对x的值取余数的意思。而在对变换域的特征进行量化与修改时，不会对feature中序号超出n*x 的元素进行量化与修改。此时f(1)，f(2)，
…
，f(x)每一组中都含有n个元素，将每组号为i中的第j个元素记为f(i)(j)，其中，i 为组号，j为组内的元素的序号，且1<＝x<＝i，1<＝j<＝n。对于上述组号为f(2*t
‑
1)与f(2*t)中的每个元素(1≤t≤x/2)，将f(2*t
‑
1)(k)与f(2*t) (k)进行比较。将f(2*t
‑
1)(k)与f(2*t)(k)分别记作a，b，并计算a/(a b)，若abs(a/(a b)
‑
0.5)<0.1，则判断a/(a b)是否大于0.5，若大于0.5，则令a＝a λ(a b),b＝b
‑
λ(a b)，反之，则令a＝a
‑
λ(a b),b＝b λ(a b)，将a，b分别保存在feature((2*t
‑
2) *n k)，feature((2*t
‑
1)*n k)中，并且如果a>b则令存储音频特征的列表feature_watermark中的feature_watermark((2*t
‑
2)*n k)＝1， feature_watermark((2*t
‑
1)*n k)＝0，反之若a<b,则令feature_watermark ((2*t
‑
2)*n k)＝0，feature_watermark((2*t
‑
1)*n k)＝1。其中，k 代表组内元素的序号，(2*t
‑
2)*n k与(2*t
‑
1)*n k的结果所代表的含义是此元素在feature，feature_watermark中的位置。在一般情况下令 n＝2。
81.s2、基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频。
82.采用扩频技术将修改过的上述特征值序列写回到音频的 dwt
‑
dct
‑
svd域中获取经过特征强化后的音频，同时为了保障水印的鲁棒性与不可感知性，在嵌入水印时采用基于启发式搜索算法的参数调整方案，不断地调整嵌入强度，以满足出版商对水印的鲁棒性与不可感知性的高标准要求。参见图4，水印嵌入过程具体为：
83.首先，将分段后的音频的每一段依次进行三级离散小波变换，离散余弦变换，奇异值分解，然后采用扩频技术嵌入水印，将feature(i) 替换原音频对应的奇异值v，之后对嵌入后含水印的音频依次进行 svd反变换，idct变换，以及idwt变换，从而得到了含水印信息的音频。
84.然后，进行音频品质与水印鲁棒性的判断。如果在遭受了模拟攻击后，仍然能够提取出水印，且音频的信噪比与水印的不可感知性符合要求，则证明所选择的嵌入强度符合要求，最终，便以此时的强度嵌入水印，从而得到包含版权水印信息并且具有极高品质的音频。此处，以λ表示水印嵌入音频的嵌入强度，并采用上述的嵌入原理来嵌入水印。在此过程中，为了得到最佳的水印嵌入强度，可利用启发式搜索方法调整参数λ，具体过程为：
85.(1)初始化参数，设置参数λ的上界max，下界min、信噪比snr 的阈值thred，此处，取thred＝30db；
86.(2)计算δ＝max
‑
min,min＝(max min)/2,λ＝roundn(min,
‑
2)，其中，λ＝roundn(min,
‑
2)表示精确到小数点后第二位；
87.(3)对含水印的音频计算其信噪比snr，若snr<thred，则 max＝λ，否则min＝λ；计算δ＝max
‑
min。
88.(4)重复步骤(2)和步骤(3)直到δ＞0.01。
89.至此获得参数λ的一个大致范围，此范围保证了水印的不可感知性，接着在取值范围[min，max]内通过基于目标函数的线性规划对参数λ进行寻优。参见图5，为参数λ的优化过程，图中，目标函数计算方法如下：
[0090]
object＝αsnr β(1
‑
mber)
[0091]
α＝(1
‑
mber)/(snr 1
‑
mber)
[0092]
β＝(1
‑
mber)/(snr 1
‑
mber)
[0093]
其中，mber表示八种基本信号处理攻击下误码率ber的均值， best_object表示最优目标函数值，best_λ表示最优参数，step为线性规划采用的步长，这里设为step＝0.001。
[0094]
按照上述步骤获取最优参数，并以此强度进行最终的嵌入。
[0095]
s3、基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频。
[0096]
在获得量化特征序列后，基于m序列将水印与量化特征序列进行绑定，从而使音频特征真正的成为水印的载体，参见图3，水印与音频特征的绑定过程为：
[0097]
首先，根据本原多项式表选择的12阶本原多项式生成对应的m序列，然后将m序列与音频量化后的特征进行异或，从而生成密钥1，再将一维的水印序列与m序列进行异或，从而生成密钥2。具体地，选定本原多项式d^12 d^11 d^10 d^7 d^6 d^5 d^4 d^1 1，根据此本原多项式生成对应的m序列m_sequence，m序列的长度为4095。读取音频水印图像，大小为64*64，并将此图像转成一维的序列 watermark。将m序列与音频量化后的特征进行异或，从而生成密钥1，再将一维的水印序列与m序列进行异或，从而生成密钥2。其中，对于feature_watermark中的每个元素，与m序列中序号与之对应的元素进行异或，将结果保存在secret1
中，用代码形式可表示为：
[0098]
secret1＝[]
[0099]
fori＝1:length(feature_watermark)
[0100]
secret1(i)＝mod(feature_watermark(i) m_sequence(i),2)；
[0101]
对于m序列中的每个元素，与水印序列watermark中序号与之对应的元素进行异或，将结果保存在secret2中，用代码形式可表示为：
[0102]
secert2＝[]
[0103]
for i＝1:length(m_sequence)
[0104]
secret2(i)＝mod(watermark(i) m_sequence(i),2)。
[0105]
此外，在实际应用中，为了实现版权的认定，上述抗任意剪切和翻录的音频水印方法还包括对嵌入音频的水印进行提取过程。采用同嵌入水印时相同的量化规则，在音频的dwt
‑
dct
‑
svd这个变换域中对水印进行提取，从而实现版权的认定。参见图6，水印的提取过程为：
[0106]
首先，采用改进的谱熵法对音频进行动态帧检测，对音频进行时域上的恢复，再对音频进行分段，对每一段分别进行三级离散小波变换，离散余弦变换，奇异值分解；然后根据嵌入时选定的规则对相应的特征进行量化；在此之后，将现有音频的特征与原始特征进行对比，如果满足条件，则将此信息送入对应的处理函数中，进行水印的恢复；最后根据两级密钥，解密出相应的水印图像。具体地：
[0107]
首先，采用改进的谱熵法(《matlab在语音信号分析与合成中的应用2013版》，宋知用)对音频进行动态帧检测，对音频进行时域上的恢复，然后将音频进行分段，每10000个样本点为一段，将每一段音频采用db1小波进行三级dwt变换，再采用dct变换，再将其进行奇异值分解后，并将奇异值保存在列表feature中，再根据上述嵌入时的原理，得出它的量化特征ex_feature_watermark，将此时的特征与原音频的量化特征feature_watermark进行对比，找到连续的12个与原音频特征完全相同的量化特征段symbol，记录下这一特征组的第一个值对应于原音频特征ex_feature_watermark的第几个值，将此位置记为 i，同时将ex_feature_watermark中从i开始的12个元素存入特征列表 symbol中，并将位置i与特征段symbol作为参数传入函数 watermark_recovery中，水印恢复watermark_recovery的工作流程为：
[0108]
1)水印恢复函数watermark_recovery会根据参数i(i表示是这段与对应的原始特征相比，是正确的特征的起始位置)与symbol,利用 secret1对特征进行解密，得到此时的m序列的12个初值；
[0109]
2)根据参数i与m序列的12个初值恢复出对应的m序列(当i 不等于1时，恢复出的m序列与原m序列相比可能有移位)，(根据m 序列的周期性)；
[0110]
3)根据i的值，将恢复出m序列的循环右移i位，得到正确嵌入时的顺序，将其命名为result；
[0111]
4)将result与secret2的每一项进行模2加法，即异或操作，得出解密后的水印序列watermark，此水印序列长度为4095，令watermark (4096)＝0，至此，便可完成全部水印的提取。
[0112]
若音频遭受大规模的减切，通过上述过程无法提取出水印时，则可以采用模式匹
配
‑
回滚法对音频进行时域上的恢复，再利用上述水印提取流程的后续步骤提取出水印。参见图7，利用模式匹配
‑
回滚法对剪切后音频进行时域恢复的大致流程为：
[0113]
1)将原音频a分段，每10000个样本点为1段，每一段分别记为 a(i)(1≤i≤n)；
[0114]
2)在剪切后的音频b中，设置一个滑动窗口，使窗口内的值为 b(j,j 9999)，此滑动窗口大小为10000，滑动值为20，其中j为窗口起始值的样本点编号；
[0115]
3)将窗口内的样本点进行5级dwt，并将结果与原音频a(i)的5 级dwt后的系数进行对照，即求相关系数；
[0116]
4)当corr>＝0.4时，认为这一段未经剪切，将它放入恢复结果列表list中，并将滑动计数器的值置为0，并且使参照对象变为原始音频 a(i)的下一段a(i 1)，令窗口起始值j＝j 10000；
[0117]
5)若corr<0.4,则令窗口继续向前滑动，滑动计数器counter的值增加相应的滑动样本点数,即令j＝j 20；
[0118]
6)反之滑动计数器的值若超过了规定的滑动次数11000个样本点，还未满足corr>＝0.4的条件，则认定这一帧音频被剪掉了，将这一帧的所有样本点补0，然后存入恢复结果列表list，并且使参照对象变为原始音频a(i)的下一段a(i 1),将滑动计数器的值置为0，音频b窗口的起始值变为j
‑
counter，窗口范围不变；
[0119]
7)重复上述步骤3)
‑
6)，直至i>n，或窗口滑至音频b的最后一个样本点。
[0120]
然后再按照上述水印提取的相同步骤继续提取水印，此处不再赘述。
[0121]
至此，则完成了本发明一种抗任意剪切和翻录的音频水印方法的全部过程。
[0122]
为了验证本发明实施例的有效性，我们通过实验来验证经过本发明的方法嵌入至音频的水印的鲁棒性和不可感知性。具体地：
[0123]
1、水印鲁棒性的测试。分别对使用本法、使用基于稳健特征点的平稳小波域数字水印算法(方法一)，以及使用基于dct
‑
clm的抗翻录攻击的鲁棒音频水印算法(方法二)嵌入水印的音频文件test1(茉莉花)和test2(blues)(音频文件均为wav格式，采样率44100，量化精度16bit)进行水印鲁棒性测试，使它们遭受常见的8种攻击与录音笔翻录，其中，翻录条件1：使用电脑内部自带的系统混音进行翻录(相当于低品质的da
‑
ad转换)；翻录条件2：使用录音笔进行翻录(不经过空气信道，利用录音笔的line
‑
in模式)，在翻录的同时，降低采样频率，(从44100降到7100)；剪切1：音频每隔400个样本点减去5 个样本点；剪切2：一次性减去音频的75％的样本点。最终取出的水印与平均误码率如图8
‑
10所示。其中，图8为使用本法嵌入水印经过攻击后的误码率以及水印完整性情况，图9为使用方法一嵌入水印的水印鲁棒性测试结果，图10为使用方法二嵌入水印的水印鲁棒性测试结果。通过对比可知，利用本方案嵌入的水印在遭受mp3压缩、噪声攻击、低通滤波、中值滤波等8种常见的信号处理攻击后误码率为零，而且在遭受任意剪切攻击和录音笔翻录后的误码率也为0。因此，本方案在音频中嵌入的水印具有极高的鲁棒性。
[0124]
2.水印的不可感知性测试。
[0125]
itu
‑
r的bs.1387标准[20]是目前国际上采用最多的评价音频不可听性的客观标准。该标准弥补了snr没有考虑到人耳听觉系统(has) 特性的缺陷，其测试工具peaq计算出的听觉质量客观区分度odg是不可听性的衡量指标。主观标准我们采用itu
‑
r bs.1116中定义的听觉质量主观区分度sdg[11]。该方法具体做法是提供给听众原始音频和含水印音频，
让他们通过打分评价二者的区分度，然后取他们打分的平均值。odg和sdg的评分标准如下表4所示。从表中可以看出， odg值和sdg值越接近0，意味着嵌入水印的不可听性越好。
[0126]
表4听觉质量客观区分度odg评价标准
[0127][0128]
对音频文件test 1(茉莉花)和test2(blues)(wav格式，采样率 44100，量化精度16bit)进行水印不可听性测试。
[0129]
原始音频信号和采用本法嵌入水印的音频信号的波形图分别如图11(a)和(b)所示，从图中可以看出，二者在波形图上几乎没有差别。
[0130]
将我们的水印方案与采用方法一和方法二的水印方案进行对比，其结果如下表5所示。
[0131]
表5 snr、odg和sdg测试结果
[0132][0133]
从表5可以看出，本系统嵌入水印后的音频具有很好的不可感知性，其snr，odg，sdg指标均优于其他方案。
[0134]
实施例2：
[0135]
第二方面，本发明还公开了一种抗任意剪切和翻录的音频水印系统，该系统包括：
[0136]
嵌入预处理模块，用于提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改以获取特征值序列，并基于所述特征值序列生成对应的量化特征序列；
[0137]
水印嵌入模块，用于基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频；
[0138]
水印绑定模块，用于基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频。
[0139]
优选的，所述嵌入预处理模块提取待嵌入水印的音频的音频特征，并对所述音频
特征进行修改以获取特征值序列包括：
[0140]
对待嵌入水印的音频信号进行分段，并对每一段依次进行三级离散小波变换和离散余弦变换，并将获得系数进行奇异值分解，从而获得音频的变换域特征；
[0141]
对音频的变换域特征进行修改以获取特征值序列。
[0142]
优选的，所述水印绑定模块基于m序列将水印与所述量化特征序列进行绑定获取包含水印信息的音频包括：
[0143]
选择本原多项式并生成对应的m序列，然后将m序列与音频的所述量化特征序列进行异或生成密钥1，并将一维的水印序列与m序列进行异或生成密钥2。
[0144]
优选的，所述系统还包括水印提取模块：采用同嵌入水印时相同的量化规则对水印进行提取，其具体过程包括：
[0145]
s41、采用改进的谱熵法对已嵌入水印的音频进行时域上的恢复；
[0146]
s42、对已嵌入水印的音频进行分段，并对每一段依次进行三级离散小波变换，离散余弦变换，奇异值分解，然后采用嵌入水印时相同的规则对变换域特征进行量化以获取量化特征序列；
[0147]
s43、将现有音频的特征与所述量化特征序列进行对比，然后结合密钥1和密钥2提取水印图像；
[0148]
s44、当利用s41
‑
s43步骤提取水印失败后，采用模式匹配
‑
回滚法对已嵌入水印的音频进行时域上的恢复，然后进行s42
‑
s43步骤以提取水印。
[0149]
可理解的是，本发明实施例提供的一种抗任意剪切和翻录的音频水印系统与上述抗任意剪切和翻录的音频水印方法相对应，其有关内容的解释、举例、有益效果等部分可以参考抗任意剪切和翻录的音频水印方法中的相应内容，此处不再赘述。
[0150]
综上所述，与现有技术相比，具备以下有益效果：
[0151]
1、本发明的技术方案针对待嵌入水印的音频，首先提取待嵌入水印的音频的音频特征，并对所述音频特征进行修改和量化以获取特征值序列和量化特征序列，基于扩频技术将所述特征值序列写回待嵌入水印的音频中获取经过特征强化后的音频，然后基于m序列将水印与所述量化特征序列进行绑定。该技术方案将m序列与dwt
‑
dct
‑
svd 域的水印技术相结合，通过该法嵌入至音频的水印可在音频遭受极大破坏的情况下仍能完整性恢复；
[0152]
2、本发明使用m序列与本原多项式，可以根据很少的连续音频特征完整地恢复出整个水印，从而有效地保护音频与视频发行人的版权；
[0153]
3、本发明在提取音频水印时，利用短时改进子带谱熵法与双门限法，对翻录后的音频在时域特征上进行恢复，以便于提取出水印；
[0154]
4、本发明在提取音频水印时，利用模式匹配
‑
回滚法，实现对遭受大面积剪切攻击后的音频的时域特征的恢复，从而间接地提高了水印的鲁棒性。
[0155]
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在
包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0156]
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。