[发明专利]基于听觉特性与整型提升小波的数字音频水印嵌入与检测方法

摘要

本发明所述的基于听觉特性与整型提升小波的数字音频水印嵌入与检测方法是一种将灰度图像(或彩色图像)嵌入到原始数字音频作品的新数字水印方法，属于信息安全、多媒体信息处理领域。本发明首先利用静态图像压缩编码技术，将水印图像(灰度图像或彩色图像)编码成一维二进制序列并进行置乱加密；然后对数字音频信号进行分段处理，并依据人类听觉系统的时域掩蔽特性自适应选择音频段做整型提升小波变换；再依据人类听觉系统的频域掩蔽特性自适应确定小波系数，并通过量化处理过程将水印信息嵌入到小波变换后的听觉重要系数中；最后通过逆整型提升小波变换及音频数据段重新组合得到嵌入水印信息的数字音频作品；并利用嵌入过程的逆过程进行数字水印检测。

主权项

1、一种基于听觉特性与整型提升小波的数字音频水印嵌入与检测方法，其特征在于包括如下具体步骤：设原始数字音频信号为A＝{a(i)，0≤i＜N}，其中，N为数据个数，a(i)∈{0，1，2，Λ，(2^p-1)}是第i个音频数据的幅度值，p为表示每个数据所使用的比特数；  设数字水印为256级灰度图像(或24位真彩色图像)W＝{w(i，j)，0≤i＜P，0 ≤j＜Q}，其中，w(i，j)代表数字水印图像的第i行、第j列象素的灰度值(或颜色值)；数字水印的嵌入过程为：(1)数字水印图像的压缩编码与置乱加密当数字水印为256级灰度图像时，数字水印图像的压缩编码与置乱加密过程为：首先采用DCT(离散余弦变换)静态图像压缩技术，将数字水印图像W编码为一维的二进制序列V，工作步骤为：(a)将数字水印图像(灰度图像)W划分成8×8的图像子块；(b)对每个图像子块进行DCT变换，并量化所得到的DCT系数；于是，可得到压缩编码后的一维二进制序列V={v(k),0≤k<Nw,v(k)∈{0,1},Nw=15(P×Q)16};]]>其中：Nw=P×Q8×8((1+4)×8+5×4)=15(P×Q)16;]]>然后对上面得到的一维二进制序列V进行置乱加密，得到新的一维二进制序列X＝{x(k)，0≤k＜N_w，x(k)∈{0，1}}；注：当数字水印为24位真彩色图像时，数字水印图像的压缩编码与置乱加密过程为：(a)将彩色水印图像W分解为R、G、B三个彩色分量；(b)将彩色分量R划分成8×8的子块；(c)对每个子块进行DCT变换，并量化所得到的DCT系数；于是，可得到压缩编码后的一维二进制序列VR={vR(k),0≤k<Nw′,vR(k)∈{0,1},Nw′=15(P×Q)16};]]>这里，Nw′=P×Q8×8×((1+4)×8+5×4)=15(P×Q)16;]]>(d)对于彩色分量G、B实施步骤(b)和(c)的相同操作，即可得到彩色水印图像W的压缩编码后一维二进制序列V={v(k),0≤k<Nw,Nw=3Nw′,v(k)∈{0,1}};]]>对上面得到的一维二进制序列V进行置乱加密，得到新的一维二进制序列X＝{x(k)，0 ≤k＜N_w，x(k)∈{0，1}}；(2)音频信号的分段处理及重要音频段的自适应选取采用对数字音频信号进行分段的方法嵌入水印信息，即首先将整个音频信号划分为若干个音频数据段，进而依据人类听觉系统(HAS)自适应选取重要音频数据段用于嵌入水印，每个音频段仅嵌入1个比特的水印信息；设每个音频段含有L个数据，为了保证数字水印的透明性与鲁棒性，音频段长度L的取值以大于等于8为宜，则可将整个音频信号划分成个音频段，即A={A(k),0≤k<N/L},]]>A(k)＝{α(kL+i)，0≤i＜L}；根据听觉时域掩蔽特性，结合音频信号内容选取能量较大的音频段作为水印嵌入的候选段；下面给出基于内容的听觉重要音频段自适应选取方法；首先，计算每个音频段的样本能量E(k)(0≤k＜N/L)：E(k)=Σl=0L-1a2(kL+i)]]>然后，对所求取的N/L个音频段样本能量值进行降序排序，并选择能量较大的前N_w个音频段A_e＝{A_e(k)，0≤k＜N_w}作为水印嵌入的候选段；(3)整型提升小波变换及小波系数的选取对所选取的每个音频数据段A_e(k)(0≤k＜N_w)分别进行3级整型提升小波变换(ILWT)       CD_e(k)＝ILWT₃(A_e(k))＝{C_e(k)(3)，D_e(k)(3)，D_e(k)(2)，D_e(k)(1)}式中，C_e(k)(3)为第k音频段A_e(k)3级小波变换后的逼近信号，而{D_e(k)(3)，D_e(k)(2)，D_e(k)(1)}为第k音频段A_e(k)3级小波变换后的细节信号；根据人类听觉系统(HAS)的频域掩蔽特性，本算法从第k音频段A_e(k)3级小波变换细节信号D_e(k)(3)中选取绝对值大的小波系数(听觉重要系数)C(k)用于水印序列元素x(k)的嵌入；(4)数字水印的嵌入水印嵌入是通过对所选择的小波系数C(k)进行特殊量化完成的；设C(k)表示待量化小波系数，C′(k)表示量化处理后的小波系数，x(k)为待嵌入水印比特，Δ表示量化间隔；将量化间隔选取为基于人类听觉系统而计算出的可感知噪声阈值(JND)，即Δ＝JND(C(k))；则量化修改系数嵌入数字水印过程(步骤)如下：首先定义量化规则，将量化函数定义为用量化间隔JND(C(k))对待量化小波系数C(k)进行二次取模求余运算，即            Q(C(k))＝[C(k)mod(JND(C(k))]mod(2)其量化规则为：如果C(k)mod(JND(C(k))为奇数，则Q(C(k))取1，即C(k)归于1类；如果C(k)mod(JND(C(k))为偶数，则Q(C(k))取0，即C(k)归于0类；然后量化系数嵌入水印信息，即对所选小波系数C(k)进行适当修改并赋予二值信息的过程，所采用的量化嵌入方案为：如果Q(C(k))不等于x(k)，则如果x(k)＝0，则否则否则不对C(k)做任何处理；其中，(5)整型提升小波逆变换及音频段的重新组合为了得到含水印的数字音频信号，要对含有水印信息的小波系数进行3级整型提升小波逆变换，即用含有水印信息的小波系数C′(k)代替C(k)并对小波系数CD_e′进行3级整型提升小波逆变换(Inverse Integer Lifting Wavelet Transform，IILWT)，即Ae′=IILWT3(CDe′)]]>将A_e′代替A_e，再结合未被选用(嵌入水印)的音频数据段，即可最终得到含有水印的数字音频信号A′；数字水印的检测过程如下：(1)将待检测数字音频信号A′划分成个音频段(每个音频段含有L个数据)，根据水印嵌入过程所生成的候选段位置信息选择听觉重要音频数据段Ae′={Ae′(k),0≤k<Nw}]]>并对其分别进行3级整型提升小波变换(采纳嵌入过程所使用的小波基)；(2)根据嵌入水印过程中所生成的用于嵌入水印的小波系数位置信息，确定出含有水印序列元素x′(k)信息的小波系数C′(k)(k＝0，1，Λ，N_w-1)，并利用前面所定义的量化规则提取水印，即(3)按照数字水印嵌入算法的相反操作，对数字水印信息X′进行逆置乱加密，以得到原始水印的一维二进制序列V′＝{v′(k)，0≤k＜N_w，v′(k)∈{0，1}}；(4)对所得到的一维二进制序列V′进行解压缩，便可恢复嵌入的数字水印图像W′＝{w′(i，j)，0≤i＜P，0≤j＜Q}。