发明名称 基于极谐变换的鲁棒数字音频水印嵌入及检测方法
摘要 本发明公开一种基于极谐变换的鲁棒数字音频水印嵌入及检测方法,是利用非下采样小波变换良好的频域局部化特征,多尺度、冗余性和平移不变性,可较好地刻画信号的非平稳性、很好地保持信号形状的特征,并且分解后各子带音频与原始音频具有相同长度的优点;加上极谐变换良好的数值稳定性,将水印嵌入在稳定的矩值上去抵抗常规信号攻击;再结合同步码技术。实验结果表明:本发明既能够对抗常规的信号处理,又可有效抵抗去同步攻击,具有良好的不可感知性和鲁棒性。此外,还具有设计简单、易于实现等优点。
申请公布号 CN105632506A 申请公布日期 2016.06.01
申请号 CN201510995516.9 申请日期 2015.12.28
申请人 辽宁师范大学 发明人 杨红颖;王思淼;王向阳;牛盼盼
分类号 G10L19/018(2013.01)I 主分类号 G10L19/018(2013.01)I
代理机构 大连非凡专利事务所 21220 代理人 闪红霞
主权项 一种基于极谐变换的鲁棒数字音频水印嵌入方法,其特征在于按照如下步骤进行:步骤11:水印预处理:将二值水印图像<img file="357798dest_path_image001.GIF" wi="17" he="18" />置乱加密为安全水印矩阵<img file="511699dest_path_image002.GIF" wi="21" he="19" />,其中<img file="331888dest_path_image003.GIF" wi="282" he="28" />再对其进行降维处理,即将二维水印图像转换为一维二进制水印序列:<img file="684372dest_path_image004.GIF" wi="525" he="28" />;步骤12:根据数字水印与同步码大小,对原始音频载体进行分段,然后利用时间域音频样本统计特性嵌入同步码,具体如下:步骤121:将<img file="564603dest_path_image005.GIF" wi="25" he="26" />按同步码长度<img file="877904dest_path_image006.GIF" wi="42" he="23" />分成<img file="298521dest_path_image006.GIF" wi="42" he="23" />段,每一段<img file="711703dest_path_image007.GIF" wi="63" he="24" />含有<img file="497257dest_path_image008.GIF" wi="11" he="15" />个音频样本,即<img file="360170dest_path_image009.GIF" wi="443" he="32" />步骤122:计算<img file="522161dest_path_image007.GIF" wi="63" he="24" />的平均值,即<img file="786921dest_path_image010.GIF" wi="342" he="52" />;步骤123:采用量化方法嵌入同步码,即对每一段<img file="8955dest_path_image007.GIF" wi="63" he="24" />,修改其均值<img file="421481dest_path_image011.GIF" wi="58" he="23" />,以嵌入一位同步码,修改策略为:<img file="121584dest_path_image012.GIF" wi="391" he="40" />其中,<img file="113286dest_path_image013.GIF" wi="209" he="29" />为修改前音频样本值,<img file="506222dest_path_image014.GIF" wi="195" he="33" />为修改后的音频样本值,且有:<img file="406045dest_path_image015.GIF" wi="433" he="118" />其中,mod为取模运算,S<sub>1</sub>为量化步长;步骤13:水印信号的嵌入,具体如下:步骤 131:将音频后半部分映射成二维矩阵的形式,并进行非下采样小波变换得到低频子带;步骤132:根据要嵌入水印位数的多少,在选取稳定矩的规则之下,选取相应数量的稳定的PHT矩进行水印嵌入,公式如下:Quantize the chosen<img file="909838dest_path_image016.GIF" wi="43" he="21" />by a quantizer q using:<img file="883610dest_path_image017.GIF" wi="108" he="25" />Modify the chosen<img file="447447dest_path_image016.GIF" wi="43" he="21" />value by:<img file="772249dest_path_image018.GIF" wi="388" he="66" />其中,<img file="814154dest_path_image019.GIF" wi="92" he="27" />为待嵌入水印的幅值,水印图像<img file="922661dest_path_image020.GIF" wi="63" he="23" />,mod为模运算函数,q为量化步长;步骤133:对选取的PHT矩嵌入水印信息前后分别进行重构,得到两个重构的二维形式的音频,对这两个重构的二维形式音频做差值运算,将得到的差值叠加到原载体二维形式的音频中,得到修改后的低频系数,再根据非下采样逆小波变换,即得到了含水印的二维形式音频,最后还原为一维形式的音频,再重构回原来的音频数据段当中得到含水印的数字音频;步骤14:重复步骤12~步骤13对其它音频数据段进行同步码与水印信息的嵌入。
地址 116029 辽宁省大连市沙河口区黄河路850号