本發明涉及數字音頻水印技術領域，尤其涉及基于嵌入區域選擇的音頻水印方法及系統。

背景技術：

數字音頻水印是向音頻信號中添加某些數字信息以達到文件真偽鑒別、版權保護、信息隱藏等目的的信號處理操作。音頻水印嵌入區域的選擇技術是指在水印嵌入到音頻信號之前，選擇合適音頻區域嵌入水印。傳統音頻水印技術，沒有考慮到音頻信號的特征，對整個音頻文件都進行水印的嵌入，這樣會導致1)音頻信號幅度低的區域嵌入水印后，幅值超出了掩蔽閾值產生噪音，破壞了感知透明性；2)對于音頻信號中出現變化劇烈的瞬態信號，該區域的音頻信號的方差很大，嵌入水印后導致檢測水印時的水印誤碼率很高；3)在頻域嵌入水印，如果選擇人耳感知不顯著的區域嵌入水印，在經過信號處理或音頻有損壓縮后，水印將會丟失一部分，導致水印檢測誤碼率高。

技術實現要素：

本發明的目的是提供選擇區域嵌入的音頻水印技術，使水印能夠嵌入到合適的音頻區域中去，避免出現不必要的噪聲以及減少誤碼的發生。

為達到上述目的，本發明提供的技術方案提供一種基于嵌入區域選擇的音頻水印方法，包括嵌入過程和檢測過程，

所述嵌入過程包括以下步驟，

步驟A1，讀取音頻文件，得到采樣率fs1和分幀后第n幀時域音頻的信號x_n，幀長為N，

先對每幀信號x_n進行是否能夠作為嵌入區域的判斷，

然后針對能夠作為嵌入區域的各幀信號x_n，進行音頻水印的嵌入頻率段的選擇，進行音頻水印的嵌入頻率段的選擇，設根據人耳感知敏感的頻率部分預設的嵌入的開始頻率為FWMIN、結束頻率為FWMAX，一幀的開始嵌入點freqmin1和嵌入結束點freqmax1求取如下，

freqmin1＝floor((FWMIN×2.0/fs1)×N)

freqmax1＝floor((FWMAX×2.0/fs1)×N)

其中，floor為向下取整函數；

步驟A2，對能夠嵌入水印的各幀信號x_n，進行離散傅立葉變換得到頻域信號X_n；

步驟A3，利用密鑰key作為隨機數種子，生成長度為freqmax1-freqmin1+1的二進制偽隨機的擴頻序列u；

步驟A4，根據擴頻序列u、頻域信號X_n和水印比特b，進行水印的嵌入，得到嵌入水印后的頻域信號，計算如下，

|X′_n|＝|X_n|+bαu

其中，α為常數，控制水印的嵌入強度，|X_n|和|X′_n|分別表示嵌入水印前的頻域幅值和嵌入水印后的頻域幅度，然后通過歐拉公式得到嵌入水印后的頻域信號

其中，∠X_n表示頻域信號的相位，X′_n表示嵌入水印后的頻域信號，e為數學自然指數；

步驟A5，將嵌入水印后的頻域信號X′_n變換到時域，生成嵌入水印的音頻文件；

所述檢測過程包括以下步驟，

步驟B1，讀取待檢測的音頻文件，得到的時域分幀后的第n幀信號z_n和采樣率fs2，

先對每幀信號x_n進行是否能夠作為嵌入區域的判斷；

針對能夠作為嵌入區域的各幀信號x_n，作為待檢測的信號，計算檢測范圍的起始點freqmin2和頻域結束點freqmax2

freqmin2＝floor((FWMIN×2.0/fs2)×N)

freqmax2＝floor((FWMAX×2.0/fs2)×N)

步驟B2，進行離散傅立葉變換得到待檢測信號的頻域信號Z_n，相應頻域幅度值記為|Z_n|；

步驟B3，利用密鑰key作為隨機數種子，生成長度為freqmax2-freqmin2+1的二進制偽隨機的擴頻序列u；

步驟B4，根據擴頻序列u和待檢測信號的頻域幅度值|Z_n|，計算出檢測的充分統計量r_n如下，

如果充分統計量r_n≥0，那么檢測到的水印比特為b＝1；否則，檢測到的水印比特為b＝0。

而且，步驟A1和步驟B1中，對每幀信號x_n進行是否能夠作為嵌入區域的判斷，實現方式如下，

1)信號x_n的平均能量的大小超出預設的相應閾值τ₁，是則為靜音區，不允許嵌入水印；

2)如果信號x_n內包含瞬態信號，則不允許嵌入水印。

而且，信號x_n內是否包含瞬態信號，通過以下方式判斷，

設將一幀信號分解為S個塊，分別計算出S個塊的能量，比較最大能量的塊與最小能量塊的能量比rate和預設的相應閾值τ₂，如果rate大于τ₂則認為該幀信號包含瞬態信號。

本發明還相應提供一種基于嵌入區域選擇的音頻水印系統，包括音頻水印嵌入子系統和水印檢測子系統，

所述音頻水印嵌入子系統包括以下模塊，

選擇合適區域嵌入模塊，用于讀取音頻文件，得到采樣率fs1和分幀后第n幀時域音頻的信號x_n，幀長為N，

先對每幀信號x_n進行是否能夠作為嵌入區域的判斷，

然后針對能夠作為嵌入區域的各幀信號x_n，進行音頻水印的嵌入頻率段的選擇，進行音頻水印的嵌入頻率段的選擇，設根據人耳感知敏感的頻率部分預設的嵌入的開始頻率為FWMIN、結束頻率為FWMAX，一幀的開始嵌入點freqmin1和嵌入結束點freqmax1求取如下，

freqmin1＝floor((FWMIN×2.0/fs1)×N)

freqmax1＝floor((FWMAX×2.0/fs1)×N)

其中，floor為向下取整函數；

第一時頻轉換模塊，用于對能夠嵌入水印的各幀信號x_n，進行離散傅立葉變換得到頻域信號X_n；

第一擴頻序列生成模塊，用于利用密鑰key作為隨機數種子，生成長度為freqmax1-freqmin1+1的二進制偽隨機的擴頻序列u；

水印嵌入模塊，用于根據擴頻序列u、頻域信號X_n和水印比特b，進行水印的嵌入，得 到嵌入水印后的頻域信號，計算如下，

|X′_n|＝|X_n|+bαu

其中，α為常數，控制水印的嵌入強度，|X_n|和|X′_n|分別表示嵌入水印前的頻域幅值和嵌入水印后的頻域幅度，然后通過歐拉公式得到嵌入水印后的頻域信號

其中，∠X_n表示頻域信號的相位，X′_n表示嵌入水印后的頻域信號，e為數學自然指數；

時頻逆變換模塊，用于將嵌入水印后的頻域信號X′_n變換到時域，生成嵌入水印的音頻文件；

所述水印檢測子系統包括以下模塊，

選擇合適區域檢測模塊，用于讀取待檢測的音頻文件，得到的時域分幀后的第n幀信號z_n和采樣率fs2，

先對每幀信號x_n進行是否能夠作為嵌入區域的判斷；

針對能夠作為嵌入區域的各幀信號x_n，作為待檢測的信號，計算檢測范圍的起始點freqmin2和頻域結束點freqmax2

freqmin2＝floor((FWMIN×2.0/fs2)×N)

freqmax2＝floor((FWMAX×2.0/fs2)×N)

第二時頻轉換模塊，用于進行離散傅立葉變換得到待檢測信號的頻域信號Z_n，相應頻域幅度值記為|Z_n|；

第二擴頻序列生成模塊，用于利用密鑰key作為隨機數種子，生成長度為freqmax2-freqmin2+1的二進制偽隨機的擴頻序列u；

相關檢測模塊，用于根據擴頻序列u和待檢測信號的頻域幅度值|Z_n|，計算出檢測的充分統計量r_n如下，

如果充分統計量r_n≥0，那么檢測到的水印比特為b＝1；否則，檢測到的水印比特為b＝0。

而且，選擇合適區域嵌入模塊和選擇合適區域檢測模塊中，對每幀信號x_n進行是否能夠 作為嵌入區域的判斷，實現方式如下，

1)信號x_n的平均能量的大小超出預設的相應閾值τ₁，是則為靜音區，不允許嵌入水印；

2)如果信號x_n內包含瞬態信號，則不允許嵌入水印。

而且，信號x_n內是否包含瞬態信號，通過以下方式判斷，

設將一幀信號分解為S個塊，分別計算出S個塊的能量，比較最大能量的塊與最小能量塊的能量比rate和預設的相應閾值τ₂，如果rate大于τ₂則認為該幀信號包含瞬態信號。

本發明提出了通過幀內最大能量與最小能量比來濾除瞬態信號提升水印檢測的準確率，通過將水印嵌入在人耳感知顯著的頻段來提升水印的魯棒性，進一步地，提出利用平均能量來濾除安靜區域提升感知透明性。本發明技術方案具有重要的市場價值。

附圖說明

圖1是本發明實施例的嵌入子系統結構框圖。

圖2是本發明實施例的檢測子系統結構框圖。

圖3是本發明實施例的嵌入過程流程圖

圖4是本發明實施例的檢測過程流程圖。

具體實施方式

下面以具體實施例結合附圖對本發明的技術方案作進一步說明。

本發明實施例提供一種基于嵌入區域選擇的音頻水印系統，包括音頻水印嵌入子系統和水印檢測子系統。

參見圖1，本發明實施例提供的嵌入區域選擇的音頻水印技術嵌入子系統，包括選擇合適區域嵌入模塊1、第一時頻轉換模塊2、第一擴頻序列生成模塊3、水印嵌入模塊4和時頻逆變換模塊5，具體實施時可以采用軟件固化技術實現各模塊。

所述選擇合適區域嵌入模塊1，對讀取的時域音頻信號幀進行判斷，具體實施時可以逐幀判斷是否能夠滿足嵌入水印的條件：不滿足就跳過此幀，繼續下一幀的判斷；如果滿足就將信號輸出給第一時頻變換模塊2，根據讀取到的時域音頻信號的采樣率和人耳較為敏感的頻率范圍計算此頻域信號嵌入水印的范圍，并將可嵌入范圍內的頻域信號輸出給水印嵌入模塊4，將該嵌入范圍的最大值和最小值輸出給第一擴頻序列生成模塊3；

所述第一時頻轉換模塊2，用于將讀取到的時域音頻信號轉換為頻域信號，輸出給水印嵌入模塊4；

所述第一擴頻序列生成模塊3，用于根據隨機數種子和選擇合適區域嵌入模塊1輸入的 嵌入范圍的最大值和最小值生成與嵌入范圍同長度的幅值為1或-1均勻分布的隨機序列，并將此隨機序列輸出給水印嵌入模塊4；

所述水印嵌入模塊4，對于頻域信號中的幅度譜，生成頻域的帶有水印信息的音頻信號輸出給時頻逆變換模塊5；

所述時頻逆變換模塊5，用于將水印嵌入模塊4輸入的頻域的帶有水印信息的音頻信號轉換為時域的帶有水印信息的音頻信號，并將此時域的帶有水印信息的音頻信號生成音頻文件，就得到帶有水印信息的音頻文件。

參見圖2，本發明實施例提供的基于相位編碼的自適應音頻水印檢測子系統，包括選擇合適區域檢測模塊6、第二時頻轉換模塊7、第二擴頻序列生成模塊8、相關檢測模塊9，具體實施時可以采用軟件固化技術實現各模塊。

所述選擇合適區域檢測模塊6與選擇合適區域嵌入模塊1的功能基本相同，不滿足水印嵌入條件的區域，一般也不含有水印，檢測時可以不用考慮：具體實施時可以逐幀判斷，對于不滿足檢測條件的幀，跳過不檢測，繼續下一幀的判斷；滿足檢測條件的音頻信號輸出給第二時頻變換模塊7，同樣將頻率檢測區域的最大值與最小值輸出給第二時頻轉換模塊7和第二擴頻序列生成模塊8；

所述第二時頻轉換模塊7，用于將讀取到的時域音頻信號轉換為頻域信號，輸出給相關檢測模塊9；

所述第二擴頻序列生成模塊8與第一擴頻序列生成模塊3的功能基本相同，將產生的結果輸出給相關檢測模塊9；

所述相關檢測模塊9，用于根據檢測范圍對輸入的待檢測的頻域幅值信號和擴頻序列生成模塊9輸入的擴頻序列，計算相關值，根據相關值的符號，判斷出水印。

各模塊具體實現參見方法相應步驟，本發明不予贅述。本發明實施例提供的基于嵌入區域選擇的音頻水印方法，包括嵌入過程和檢測過程。

參見圖3，本發明實施例提供的基于選擇區域的音頻水印嵌入過程可以采用計算機軟件技術手段自動進行流程，具體包括以下步驟：

步驟A1，讀取音頻文件，對時域的音頻信號x先分幀，得到采樣率fs1和分幀后的第n幀時域音頻信號x_n(幀長為N)，對每幀信號x_n進行是否能夠作為嵌入區域的判斷，判斷包含兩方面的判斷：

1)判斷x_n的平均能量的大小是否超出設定的閾值，來判斷當前幀x_n是否為靜音區，如 果是靜音區就不允許嵌入水印，否則超出閾值就不是靜音區，可能進行嵌入。通過下面的公式計算第n幀的平均能量

其中，N為幀長，即一幀內的樣本點數；i為一幀內的樣本點索引序號，取值在0到N-1之間；x_n²(i)表示第n幀時域信號x_n在幀內第i點的能量；τ₁為平均能量的判決閾值，具體實施時本領域技術人員可自行預設取值，例如根據經驗得到；如果超出閾值，則滿足條件1)，進行下面條件2)的判斷。

2)對于一幀內出現瞬態信號的情況，由于其頻率劇烈變化，會造成的較大的方差，在檢測時信號方差越大造成的水印檢測的錯誤概率越高，這種情況也不應該嵌入水印。通過將一幀分解為S個塊，分別計算出S個塊的能量，通過最大能量的塊與最小能量塊的能量比rate和閾值τ₂的比較，rate大于τ₂則認為是該幀信號包含瞬態信號不予嵌入水印，否則可以嵌入水印。具體實施時，本領域技術人員可自行預設S的取值。

具體實現方式如下：

首先將一幀信號x_n分成S個塊，則每個子塊內的樣本點數M為

M＝N/S (2)

每個塊的能量E_i計算如下

其中，i表示幀內塊的索引序號，j表示幀內樣本點的索引序號，x_n²(j)表示第n幀時域信號x_n在幀內第j點的能量。

找出塊能量中的最大能量E_Max和最小能量E_Min

E_Max＝MAX{E_i}，E_Min＝MIN{E_i}，i∈[0,S-1] (4)

其中，MAX，MIN分別表示求最大值函數和最小值函數。

最大能量和最小能量的比rate計算如下：

如果rate＞τ₂，就認為信號幀x_n內存在瞬態信號，該幀不嵌入水印；否則，可以嵌入水印。其中τ₂為閾值，具體實施時本領域技術人員可自行預設取值，例如τ₂為瞬態信號的檢測閾值，根據經驗得到。

然后針對能夠作為嵌入區域的各幀信號x_n，對于音頻水印的嵌入頻率段的選擇，應為人耳感知較為顯著的區域，本領域技術人員可根據人耳感知特性自行預先設定，例如1000-7000Hz。因為這些區域的信號在經過濾波、音頻壓縮等攻擊后，不會被去除。所以將水印嵌入到感知明顯的區域，在經受一些信號攻擊后不會被抹掉，能夠檢測出來。設設根據人耳感知敏感的頻率部分預設的嵌入的開始頻率為FWMIN、結束頻率為FWMAX，對應一幀的開始嵌入點freqmin1和嵌入結束點freqmax1求取如下，

freqmin1＝floor((FWMIN×2.0/fs1)×N) (6)

freqmax1＝floor((FWMAX×2.0/fs1)×N) (7)

其中，floor為向下取整函數。

根據開始嵌入點freqmin1和嵌入結束點freqmax1，選取此范圍內的頻域音頻信號。

具體實施時可以逐幀判斷，不滿足條件的跳過，進行下一幀的判斷。

步驟A2，對能夠嵌入水印的信號幀x_n，進行FFT變換(快速離散傅立葉變換)為頻域信號X_n。

步驟A3，利用密鑰key作為隨機數種子，生成長度為freqmax1-freqmin1+1的二進制偽隨機擴頻序列u。

在MATLAB中的實施例具體過程如下：

首先，利用密鑰key，調用RandStream函數(隨機種子函數)對rand函數(隨機數生成函數)進行初始化，然后調用rand函數生成隨機數，由于rand函數生成的隨機數是0～1之間的數，還需對這些數進行四舍五入變成0和1的二進制偽隨機序列，然后將此單極性的偽隨機序列，轉為雙極性只含有+1和-1的偽隨機序列u。

步驟A4，根據擴頻序列u、頻域信號X_n和水印比特b，利用下面的公式(8)進行水印的嵌入，得到嵌入水印后的頻域信號，計算實現如下

|X′_n|＝|X_n|+bαu (8)

其中，α為常數，控制水印的嵌入強度，具體實施時本領域技術人員可預設取值；|X_n|和|X′_n|分別表示嵌入水印前的頻域幅值和嵌入水印后的頻域幅度，然后通過歐拉公式得到嵌入 水印后的頻域信號。

其中，∠X_n表示頻域信號的相位，X′_n表示嵌入水印后的頻域信號，e為數學自然指數。

步驟A5，將嵌入水印后的頻域信號X′_n變換到時域，最后生成音頻文件，即得到嵌入水印的音頻文件。

參見圖4，本發明實施例提供的基于選擇區域嵌入的音頻水印檢測過程，可以采用計算機軟件技術手段自動進行流程，具體包括以下步驟：

步驟B1，讀取待檢測的音頻文件，得到的時域分幀后的第n幀信號z_n和采樣率fs2，對各時域信號z_n采取步驟A1中一樣的判決方法，

即考慮如下兩個條件，

1)信號x_n的平均能量的大小超出預設的相應閾值τ₁，是則為靜音區，不允許嵌入水印；

2)如果信號x_n內包含瞬態信號，則不允許嵌入水印。

則不為靜音區且不包含瞬態信號的幀信號，能夠嵌入水印并有待檢測。

具體實施時可以逐幀判斷，不滿足條件的跳過，進行下一幀的判斷。

針對能夠作為嵌入區域的各幀信號x_n，作為待檢測的信號，計算檢測范圍的頻域起始點freqmin2和頻域結束點freqmax2

freqmin2＝floor((FWMIN×2.0/fs2)×N) (10)

freqmax2＝floor((FWMAX×2.0/fs2)×N) (11)

步驟B2，對于滿足檢測條件的信號z_n，進行離散傅立葉變換得到待檢測信號的頻域信號Z_n，相應頻域幅度值記為|Z_n|。

步驟B3，利用密鑰key，生成二進制擴頻序列u(與上面嵌入方法得到的u方式相同)，即利用密鑰key作為隨機數種子，生成長度為freqmax2-freqmin2+1的二進制偽隨機擴頻序列u。

步驟B4，根據擴頻序列u和待檢測信號的頻域幅度值|Z_n|，通過計算擴頻序列u和待檢測信號的頻域幅度值|Z_n|的相關值，計算出檢測的充分統計量r_n

其中，<·>表示信號的內積計算。

如果充分統計量r_n≥0，那么檢測到的水印比特為b＝1；否則，檢測到的水印比特為b＝0。

本發明中所描述的具體實施例僅僅是對本發明精神作舉例說明。本發明所屬技術領域的技術人員可以對所描述的具體實施例做各種各樣的修改或補充或采用類似的方式替代，但并不會偏離本發明的精神或者超越所附權利要求書所定義的范圍。