CN109346091B - 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统 - Google Patents

一种基于等长熵码字替换的压缩音频自适应隐写方法和系统 Download PDF

Info

Publication number
CN109346091B
CN109346091B CN201811067226.8A CN201811067226A CN109346091B CN 109346091 B CN109346091 B CN 109346091B CN 201811067226 A CN201811067226 A CN 201811067226A CN 109346091 B CN109346091 B CN 109346091B
Authority
CN
China
Prior art keywords
code word
bit string
secret
huffman
huffman code
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201811067226.8A
Other languages
English (en)
Other versions
CN109346091A (zh
Inventor
易小伟
杨坤
赵险峰
于海波
刘长军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Information Engineering of CAS
Original Assignee
Institute of Information Engineering of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Information Engineering of CAS filed Critical Institute of Information Engineering of CAS
Priority to CN201811067226.8A priority Critical patent/CN109346091B/zh
Publication of CN109346091A publication Critical patent/CN109346091A/zh
Application granted granted Critical
Publication of CN109346091B publication Critical patent/CN109346091B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/018—Audio watermarking, i.e. embedding inaudible data in the audio signal

Landscapes

  • Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明公开了一种基于等长熵码字替换的压缩音频自适应隐写方法和系统。本方法的信息嵌入部分为：建立哈夫曼码字映射表，生成可相互替换的哈夫曼码字对，根据映射表将音频帧中大值区的可用哈夫曼码字映射成比特串，并计算每个比特的修改代价；利用STC编码将秘密信息嵌入到哈夫曼码字映射成的比特串中，得到载密二进制串；根据哈夫曼码字映射表将载密比特串逆映射成原始的哈夫曼码字，得到载密哈夫曼码字序列。信息提取部分根据哈夫曼码字映射表将音频帧大值区中可隐写码字映射成二进制比特串，利用STC解码从二进制比特串中提取秘密信息。本发明的隐写方法与现有方法相比具有更高的安全性、更大的隐写容量和更高的嵌入速率。

Description

一种基于等长熵码字替换的压缩音频自适应隐写方法和系统
技术领域
本发明涉及一种信息隐藏技术，尤其是涉及一种以MP3等压缩音频为载体的信息自适应嵌入及提取的方法和系统。
背景技术
随着互联网信息技术的迅速发展，信息安全问题日益严重，已经成为当前社会的热点，并引起学术界的重视。隐写术是确保信息安全传输、实现隐蔽通信的重要方法，是信息安全领域的重要分支。隐写术主要是通过隐藏通信行为实现秘密信息的安全传输，它将秘密信息嵌入到数字媒体文件，如音视频、图像和文档，并且不会引起原始载体的视觉和听觉上的感知失真，使秘密信息在不引起第三方注意的情况下完成传递。
音频作为最常见的多媒体形式之一，是一种良好的信息隐藏载体。目前，已有多种音频隐写方法适用于时域音频和压缩域音频。由于时域音频需要较大的存储空间和传输带宽，给存储和传播带来挑战，因此互联网上的音频多数以MP3和AAC等压缩格式存在。鉴于压缩格式音频的广泛传播和使用，以压缩格式音频为载体的隐写算法具有较大的实用价值。压缩域音频隐写算法已经逐渐成为音频隐写的研究重点。英国剑桥大学的Petitcolas开发的MP3Stego成为最经典、使用最广泛的MP3隐写软件之一。MP3Stego算法通过调整量化编码块长度的奇偶性来实现信息嵌入。严迪群等人为解决MP3Stego算法在低码率嵌入下的死循环问题，提出了改进的基于量化步长奇偶性的隐写方法(参考文献：严迪群.压缩域音频隐写与隐写分析中若干问题的研究[D].宁波大学,2012)。此外，严迪群等人还提出了基于码字替换的大容量隐写算法(参考文献：Yan D Q,Wang R D,Zhang L G.A high capacityMP3steganography based on Huffman coding[J].Journal of Sichuan University,2011,48(6):1281-1286)。张垚等人提出了直接修改小值区哈夫曼码字的隐写算法(参考文献：张垚,潘峰,申军伟.基于MP3的后置式自适应隐写算法[J].计算机科学,2016,43(8):114-117)和基于哈夫曼码表索引的隐写算法(参考文献：张垚,潘峰,申军伟.基于MP3的内嵌型自适应隐写算法[J].计算机工程与设计,2016,37(6):1537-1542)。董亚坤提出了基于符号位修改的MP3隐写算法和基于linbits位修改的MP3隐写算法(参考文献：董亚坤.基于MP3的信息隐藏技术研究[D].北京邮电大学,2015)。相关的发明专利主要包括：中国发明专利CN103106901B提出了基于音频颗粒索引值的MP3隐写方法；中国发明专利CN102097098A提出了以窗口类型切换规则为依据的压缩域音频隐写方法；中国发明专利CN106228981A提出了一种直接修改频域量化系数的MP3自适应隐写算法。
经分析，已有针对MP3的隐写算法存在安全负载率低和抗检性差等问题。例如，针对MP3Stego等算法的检测正确率已经能达到90％以上(参考文献：陈益如,王让定,严迪群.基于Huffman码表索引的MP3Stego隐写分析方法[J].计算机工程与应用,2012,48(9):124-126；王让定,羊开云,严迪群,等.一种针对MP3Stego隐写后的音频的隐写检测方法,CN104282310A[P].2015)。
发明内容
针对现有MP3等压缩音频隐写方法容量低和安全性差等问题，本发明提出了一种基于等长熵码字替换的压缩音频自适应隐写方法和系统。
本发明所述“等长熵”是指哈夫曼编码表中码字长度相等的码字集合所具备的一种属性，它反映了码字在编码过程中出现的统计频次相等或很接近。该方法通过构造特定的哈夫曼码字空间，生成可相互替换哈夫曼码字对，并对所有可隐写码字进行奇偶分配，以码字奇偶分配后的二进制串为载体进行信息嵌入。信息嵌入之前，以心理声学模型中人耳静音阈值为依据构造合理的代价函数，计算每个可嵌入点的修改代价。信息嵌入以STC(Syndrome-Trellis Codes)编码实现，结合每个嵌入点的修改代价，实现最优路径选择，使得嵌入总失真最小。本发明能够实现大容量隐写，能够满足实际应用的需求，并且实验证明本发明设计的自适应隐写方法在较高负载率下也能有效抵抗现有的多种MP3等压缩音频隐写分析方法，具有较好的抗检性。
本发明是一种自适应隐蔽通信方法，能在大容量隐写的前提下达到较高的安全性，它主要包括秘密信息嵌入和信息提取两部分。
秘密信息嵌入部分发生在MP3等压缩音频的音频编码过程中量化编码之后和哈夫曼编码之前，在音频编码的过程中实现信息的嵌入操作，包括如下子步骤：
(S1)构造哈夫曼码字映射表，生成可相互替换的哈夫曼码字对。在任意哈夫曼编码表内，两个哈夫曼码字相互替换必须满足三个条件，即码字长度相等、码字对应系数符号位个数相等、码字linbits位一致。满足以上三个条件的两个哈夫曼码字相互替换不会影响码流的结构和长度，并且不会破坏音频的正常编解码。每一张码表中的哈夫曼码字可以根据此码表内是否存在可与其相互替换的码字，划分为可隐写码字和不可隐写码字。对可隐写码字空间中的哈夫曼码字进行奇偶分配，进一步划分为两个子空间，其中一个子空间中的码字表示比特‘0’，另一个空间中的码字表示比特‘1’。
(S2)MP3编码音频帧包含两个颗粒，每个颗粒包含576个频域系数。编码过程中系数按照顺序被分为大值区、小值区和零值区(MP3编码标准中的术语)，本发明中仅使用大值区中的哈夫曼码字进行信息嵌入。针对每个帧音频，根据步骤(S1)构造的哈夫曼码字映射表将两个颗粒中大值区所有的可隐写码字映射为比特串C。在映射过程中，依据心理声学模型与哈夫曼码字对应系数的修改幅度计算所有可隐写码字的修改代价ρ。
(S3)将比特串C和修改代价ρ使用相同的算法及密钥进行置乱。
(S4)使用STC编码进行信息嵌入。STC编码器根据嵌入点的修改代价，自适应的选取最优嵌入路径将秘密信息嵌入到比特串C，使修改总失真最小，得到载密的比特串S。
(S5)使用与步骤(S3)相同的置乱算法与密钥将S进行逆置乱。
(S6)将载密比特串S逆映射成哈夫曼码字。如果当前码字是可隐写码字，则根据步骤(S1)中构造的哈夫曼码字映射表和当前哈夫曼码字，将S中的每个比特逆映射成原始的哈夫曼码字。
(S7)重复步骤(S2)～步骤(S6)，直到所有的秘密信息嵌入完成或音频编码结束，生成新的隐藏秘密信息的MP3音频文件。
秘密信息提取部分发生在MP3等压缩音频的音频解码中哈夫曼解码过程和频域系数反量化过程之间，在音频解码的过程中实现信息的提取操作，包括如下子步骤：
(P1)根据嵌入过程步骤(S1)中哈夫曼码字映射表的构造方法构造哈夫曼码字映射表。
(P2)对音频进行部分解码，得到哈夫曼码字。根据哈夫曼码字映射表将每个可隐写码字映射成比特串S(S是载密比特串)。
(P3)使用步骤(S3)中的置乱算法与密钥对S进行置乱。
(P4)使用STC编码从S中提取秘密信息，提取过程使用的参数如生成矩阵的宽度和高度与嵌入过程中使用的参数保持一致。
(P5)重复上述步骤(P2)～步骤(P4)，直到音频解码结束或信息完整提取，得到还原的秘密信息。
与上面方法对应地，本发明还提供一种基于等长熵码字替换的压缩音频自适应隐写系统，其包括码字-二进制映射模块、失真函数模块、置乱模块、STC编码模块、逆置乱模块和二进制-码字映射模块；所述码字-二进制映射模块将原始哈夫曼码流转换为二进制比特串；所述失真函数模块根据心理声学模型和系数修改幅度计算每个可嵌入位置的修改代价；所述置乱模块将二进制比特串和修改代价使用相同的置乱算法和密钥进行置乱；所述STC编码模块利用STC编码结合可嵌入位置的修改代价，将待嵌信息嵌入到置乱后的比特串中，得到载密比特串；所述逆置乱模块将载密比特串采用相同的置乱算法和密钥进行逆置乱；所述二进制-码字映射模块将逆置乱后的载密比特串逆映射成原始的哈夫曼码。
进一步地，在提取秘密信息时，所述码字-二进制映射模块根据哈夫曼码字映射表将哈夫曼码字串映射成比特串；所述置乱模块将比特串进行置乱；所述STC编码模块利用STC解码使用STC编码从置乱后的比特串中提取秘密信息。
本发明与现有技术相比的有益效果在于：
(1)本发明中，隐写算法具有更高的安全性。与传统的非自适应MP3隐写算法相比，本发明计算每个可嵌入点的修改代价，并且结合了STC编码实现自适应隐写，每次嵌入选择最佳嵌入路径使修改总失真最小，提高隐写方法的抗隐写分析能力。
(2)本发明与多数传统的MP3隐写算法相比，具有更大的隐写容量。多数MP3隐写方法隐写容量相对较小，难以满足实际应用的需要。本发明中的方法，320kbps码率音频最大隐写容量可达到11kbps，128kbps码率时最大隐写容量也能达到5kbps。
(3)本发明与专利CN106228981A中的音频隐写方法相比，具有更高的嵌入速率。CN106228981A中的自适应隐写方法不能确保一次就嵌入成功，可能需要多次重复嵌入；本发明中的隐写方法不需要迭代嵌入的方式，只需一次嵌入就能成功，降低了嵌入的时间开销。
附图说明
图1是本发明信息隐写方法的实现流程图；
图2是与图1所述方法对应的信息隐写系统的模块构成及工作流程图；
图3是本发明秘密信息隐写过程示意图；
图4是本发明码表空间搜索示意图；
图5是本发明信息提取方法实现流程图；
图6是本发明秘密信息提取过程示意图。
具体实施方式
本发明的基于等长熵码字替换的压缩音频自适应隐写方法包括信息嵌入和信息提取，具体实施详细说明如下。
图1是本发明信息隐写方法的实现流程图，图2是与该方法对应的隐写系统的模块构成及工作流程图。主要可以分为六个模块：码字-二进制映射模块、失真函数模块、置乱模块、STC编码模块、逆置乱模块和二进制-码字映射模块。首先，原始哈夫曼码流Hc经过码字-二进制映射模块会被转换为载体比特串C。其次，失真函数模块根据心理声学模型和系数修改幅度计算每个可嵌入位置的修改代价ρ。置乱模块将C和ρ使用相同的置乱算法和密钥进行置乱，得到置乱后的比特串C1和置乱后的修改代价ρ1。接着，STC编码模块利用STC编码结合可嵌入位置的修改代价，将待嵌信息M嵌入到置乱后的比特串C1中，得到载密比特串S1。逆置乱模块将S1采用相同的置乱算法和密钥进行逆置乱，得到比特串S。最后，经过二进制-码字映射模块将S逆映射成原始的哈夫曼码Hs。嵌入算法的数学关系式可表示为：
Figure BDA0001798603990000051
图3是本发明秘密信息隐写过程示意图，显示了本发明在信息嵌入过程中数据流的变化。嵌入操作完成后，最终的载密码流与原始码流相比，有两个码字Hk和Hk+2发生了改变。信息嵌入方法的具体实现过程如下：
步骤1：基于等长熵码字替换的隐写方法首先要构造哈夫曼码字映射表，在此之前要确定哈夫曼码字可替换的条件。大值区两个频域系数编码后的码流结构，称为一个编码单元，信息嵌入操作不能改变此编码单元的长度和结构，否则会造成编解码错误。假定
Figure BDA0001798603990000052
表示第k个码表的第i个码字，
Figure BDA0001798603990000053
表示
Figure BDA0001798603990000054
对应的QMDCT(Quantized ModifiedDiscrete Cosine Transform)系数。因此，对于
Figure BDA0001798603990000055
可相互替换的哈夫曼码字对
Figure BDA0001798603990000056
和
Figure BDA0001798603990000057
必须满足如下三个条件：
(C1)
Figure BDA0001798603990000058
的长度和
Figure BDA0001798603990000059
的长度相等，L·表示计算哈夫曼码字长度，即：
Figure BDA00017986039900000510
(C2)
Figure BDA00017986039900000511
的符号位个数和
Figure BDA00017986039900000512
符号位个数相等，Θ(·)表示计算QMDCT系数对非零系数的个数，即：
Figure BDA00017986039900000513
(C3)
Figure BDA00017986039900000514
和
Figure BDA00017986039900000515
的linbits存在性一致，
Figure BDA00017986039900000516
和
Figure BDA00017986039900000517
的linbits存在性一致，G(·)表示计算QMDCT系数的linbits位数，即：
Figure BDA00017986039900000518
根据上述三个条件，哈夫曼码字映射表的详细构造方法如下：
(1)记集合Π(k)表示包含第k张码表中所有的哈夫曼码字。
(2)将Π(k)分为两个不重叠的子集：
Figure BDA00017986039900000519
和
Figure BDA00017986039900000520
其中
Figure BDA00017986039900000521
表示不可隐写码字集合，
Figure BDA00017986039900000522
表示可隐写码字集合。首先，将
Figure BDA00017986039900000523
和
Figure BDA00017986039900000524
都初始化为
Figure BDA00017986039900000525
对于
Figure BDA00017986039900000526
Figure BDA0001798603990000061
如果
Figure BDA0001798603990000062
满足上述三个相互替换条件，就将
Figure BDA0001798603990000063
和
Figure BDA0001798603990000064
移入
Figure BDA0001798603990000065
中，否则放入
Figure BDA0001798603990000066
重复上述过程，直到
Figure BDA0001798603990000067
(3)划分
Figure BDA0001798603990000068
为子空间
Figure BDA0001798603990000069
和
Figure BDA00017986039900000610
依据(2)假定找到
Figure BDA00017986039900000611
为一对可替换码字,并按照zig-zag顺序遍历假定为
Figure BDA00017986039900000612
则将先序
Figure BDA00017986039900000613
分配到集合
Figure BDA00017986039900000614
将后序
Figure BDA00017986039900000615
分配到集合
Figure BDA00017986039900000616
其中，在隐写嵌入时使用
Figure BDA00017986039900000617
中的码字表示比特‘0’，而
Figure BDA00017986039900000618
中的码字表示比特‘1’。
根据哈夫曼码字统计分布特征发现，一般地，两个可相互替换的码字在以zig-zag顺序搜索时距离最近。如图4所示，是第7号哈夫曼码表进行码字空间搜索原理示意图。根据可隐写码字条件(C2)，大值区的QMDCT系数对被划分为3个子区域R0、R1和R2，每个区域都独立进行搜索，搜索顺序如图4中箭头所示的方向。以长度为8的哈夫曼码字为例，存在h<2,3>、h<3,2>、h<5,1>、h<4,2>、h<2,4>和h＜1,5＞共六个码字都能同时符合条件(C1)-(C3)。按照zig-zag顺序遍历时，6个码字组成了3个码字对：
Figure BDA00017986039900000619
并按照码字搜索顺序将h<2,3>、h<5,1>、h<2,4>放入
Figure BDA00017986039900000620
将h<3,2>、h<4,2>、h<1,5>放入
Figure BDA00017986039900000621
其中h的码字-二进制映射过程fctb(h)可以表示为：
Figure BDA00017986039900000622
相反地，其逆映射二进制-码字映射过程fbtc(m,g)可以表示为：
Figure BDA00017986039900000623
其中，m表示经过STC编码后的比特值，g表示m对应的哈夫曼码字，
Figure BDA00017986039900000624
表示g的替换码字。
步骤2：计算每一个可嵌入位置的修改代价。信息嵌入操作是基于可隐写码字替换实现的，码字替换操作最终会造成QMDCT系数的变化。嵌入点的修改失真主要受两个方面影响，一方面是QMDCT系数的修改幅度，另一方面是人耳对QMDCT系数频线的听觉灵敏度。
(1)使用哈夫曼码字hi替换为hj时，即从<xi,yi>变成<xj,yj>，QMDCT系数对之间的变化可用平面点的曼哈顿距离dij来计算，
dij＝|xi-xj|+|yi-yj|
其中，dij越大表示使用这对哈夫曼码字进行隐写嵌入带来的失真越大。
(2)人耳对不同频率的听觉灵敏度可用绝对静音阈值衡量。绝对静音阈值Tf是通过大量实验得到的，其数学式可近似表示为，
Figure BDA0001798603990000071
其中，f表示频率值，Tf越小表示人耳对当前频率越敏感。
综合(1)和(2)可得，本发明中隐写方法的代价函数定义为，
Figure BDA0001798603990000072
其中，i是哈夫曼码字在颗粒中的索引值，ρi表示在第i个哈夫曼码字进行隐写嵌入的失真代价，Tf2i表示第2个频线的绝对听觉阈值，σ是定值常量，保证
Figure BDA0001798603990000073
值大于1。
依据码字映射公式，当前码字h在集合
Figure BDA0001798603990000074
中时将其映射为比特‘0’，在集合
Figure BDA0001798603990000075
中时将其映射为比特‘1’，h在集合
Figure BDA0001798603990000076
中时表示h为不可隐写码字，在码字映射过程中根据代价函数计算当前嵌入位置的修改代价。当前帧遍历结束时，得到比特串C和对应的修改代价ρ。
步骤3：将比特串C和对应的修改代价ρ进行置乱，得到C1和ρ1。本发明中使用Logistic混沌映射算法对C和ρ进行置乱，对C和ρ进行置乱时使用相同的代价置乱密钥。置乱的目的是使得嵌入点的代价能够均匀分布，达到更优的嵌入路径选择。
步骤4：利用STC编码，将秘密信息M嵌入到载体C1中，得到载密的比特串S1。嵌入过程可表示为，
S1＝STC(C1,M,ρ,height,width)
其中，height和width分别是STC编码器的子矩阵的高度和宽度。STC编码能够根据修改代价选择最优嵌入路径，使得每次嵌入总失真最小。
步骤5：使用步骤3中的混沌置乱算法，并使用相同的密钥将载密串S1进行逆映射得到S，使其次序与原始次序相同。
步骤6：遍历帧，如果大值区哈夫曼码字hi在集合
Figure BDA0001798603990000077
中，则根据码字逆映射公式将hi替换成hj，否则跳过不作处理，继续遍历下一个码字。逆映射之后得到载密的哈夫曼码字串，继续后续的MP3编码流程。
步骤7：重复上述步骤2～步骤6，直到遍历所有音频帧或秘密消息全部嵌入完成，得到载密MP3音频。
图5是本发明信息提取方法实现流程图。信息提取与信息嵌入是相反的过程，主要分为三个部分，首先在码字-二进制映射模块根据映射表将哈夫曼码字串Hs映射成比特串S；其次，将S进行置乱得到S1；最后，利用STC解码从S1中提取信息M。图6是本发明秘密信息提取过程实例示意图，信息提取方法的具体实现过程：
步骤1：使用信息嵌入过程中哈夫曼码字映射表的构造方法生成映射表，具体过程如下：
(1)记集合Π(k)表示包含第k张码表中所有的哈夫曼码字。
(2)将Π(k)分为两个不重叠的子集：
Figure BDA0001798603990000081
和
Figure BDA0001798603990000082
其中
Figure BDA0001798603990000083
表示不可隐写码字集合，
Figure BDA0001798603990000084
表示可隐写码字集合。首先，将
Figure BDA0001798603990000085
和
Figure BDA0001798603990000086
都初始化为
Figure BDA0001798603990000087
对于
Figure BDA0001798603990000088
Figure BDA0001798603990000089
如果
Figure BDA00017986039900000810
满足上述三个相互替换条件，就将
Figure BDA00017986039900000811
和
Figure BDA00017986039900000812
移入
Figure BDA00017986039900000813
中，否则放入
Figure BDA00017986039900000814
重复上述过程，直到
Figure BDA00017986039900000815
(3)划分
Figure BDA00017986039900000816
为子空间
Figure BDA00017986039900000817
和
Figure BDA00017986039900000818
依据(2)假定找到
Figure BDA00017986039900000819
为一对可替换码字,并按照zig-zag顺序遍历假定为
Figure BDA00017986039900000820
则将先序
Figure BDA00017986039900000821
分配到集合
Figure BDA00017986039900000822
将后序
Figure BDA00017986039900000823
分配到集合
Figure BDA00017986039900000824
其中，在隐写嵌入时使用
Figure BDA00017986039900000825
中的码字表示比特‘0’，而
Figure BDA00017986039900000826
中的码字表示比特‘1’。
步骤2：对于每个音频帧，当大值区中的哈夫曼码字在集合
Figure BDA00017986039900000827
中时，根据码字映射公式
Figure BDA00017986039900000828
将哈夫曼码字序列映射为比特串S，如果h在集合
Figure BDA00017986039900000829
中则将其映射成比特‘0’；如果h在集合
Figure BDA00017986039900000830
中则将其映射成比特‘1’；否则，跳过当前码字。
步骤3：使用混沌置乱算法对S进行置乱操作，置乱算法密钥设置与嵌入过程置乱算法密钥相同，得到新的比特串S1。
步骤4：使用STC解码器从载密串S1中提取信息，即，
M＝STC-1(S1,height,width)
其中，height是STC解码器中子矩阵高度，width是STC解码器中子矩阵宽度。
步骤5：重复上述步骤2～步骤4，直到所有音频帧结束或者信息完全提取。
本发明中，置乱算法和代价函数可以采用其它不同的算法替代，但不影响系统结构和流程。另外，本发明方法同样适用于AAC(MPEG-4标准)压缩音频。
以上实施例仅用以说明本发明的技术方案而非对其进行限制，本领域的普通技术人员可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明的精神和范围，本发明的保护范围应以权利要求书所述为准。

Claims (8)

1.一种基于等长熵码字替换的压缩音频自适应隐写方法，其特征在于，采用以下步骤嵌入秘密信息：
(S1)构造哈夫曼码字映射表，生成可相互替换的哈夫曼码字对；
(S2)将大值区中的可隐写哈夫曼码字映射为二进制比特串，并计算可隐写哈夫曼码字的修改代价；
(S3)将二进制比特串和修改代价使用相同的置乱算法及密钥进行置乱；
(S4)使用STC编码将秘密信息嵌入到二进制比特串中，得到载密二进制比特串；
(S5)将载密二进制比特串使用与步骤(S3)相同的置乱算法与密钥进行逆置乱；
(S6)将逆置乱后的载密二进制比特串逆映射成哈夫曼码字；
(S7)重复步骤(S2)～步骤(S6)，直到所有的秘密信息嵌入完成或音频编码结束，生成新的隐藏秘密信息的音频文件；
步骤(S1)所述哈夫曼码字映射表中，两个哈夫曼码字相互替换必须满足三个条件，即码字长度相等、码字对应系数符号位个数相等、码字linbits位一致；
所述哈夫曼码字映射表的构造方法如下：
(1)记集合Π(k)表示包含第k张码表中所有的哈夫曼码字；
(2)将∏(k)分为两个不重叠的子集：
Figure FDA0003263959980000011
和
Figure FDA0003263959980000012
其中
Figure FDA0003263959980000013
表示不可隐写码字集合，
Figure FDA0003263959980000014
表示可隐写码字集合：首先将
Figure FDA0003263959980000015
和
Figure FDA0003263959980000016
都初始化为
Figure FDA0003263959980000017
对于
Figure FDA0003263959980000018
Figure FDA0003263959980000019
其中i≠j，如果
Figure FDA00032639599800000110
满足三个相互替换条件，就将
Figure FDA00032639599800000111
和
Figure FDA00032639599800000112
移入
Figure FDA00032639599800000113
中，否则放入
Figure FDA00032639599800000114
重复上述过程，直到
Figure FDA00032639599800000115
(3)划分
Figure FDA00032639599800000116
为子空间
Figure FDA00032639599800000117
和
Figure FDA00032639599800000118
依据步骤(2)假定找到
Figure FDA00032639599800000119
为一对可替换码字，并按照zig-zag顺序遍历，将先序
Figure FDA00032639599800000120
分配到集合
Figure FDA00032639599800000121
将后序
Figure FDA00032639599800000122
分配到集合
Figure FDA00032639599800000123
其中，在隐写嵌入时使用
Figure FDA00032639599800000124
中的码字表示比特‘0’，而
Figure FDA00032639599800000125
中的码字表示比特‘1’。
2.根据权利要求1所述的方法，其特征在于，步骤(S2)依据心理声学模型与哈夫曼码字对应系数的修改幅度计算所有可隐写哈夫曼码字的修改代价。
3.根据权利要求1所述的方法，其特征在于，步骤(S2)包括：
(1)采用平面点的曼哈顿距离dij来计算QMDCT系数对之间的变化；
(2)采用绝对静音阈值衡量人耳对不同频率的听觉灵敏度，绝对静音阈值Tf表示为，
Figure FDA0003263959980000021
其中，f表示频率值，Tf越小表示人耳对当前频率越敏感；
(3)综合步骤(1)和(2)得到失真函数：
Figure FDA0003263959980000022
其中，i是哈夫曼码字在颗粒中的索引值，ρi表示在第i个哈夫曼码字进行隐写嵌入时的失真代价，
Figure FDA0003263959980000023
表示第2i个频线的绝对听觉阈值，σ是定值常量，保证
Figure FDA0003263959980000024
值大于1。
4.根据权利要求1所述的方法，其特征在于，步骤(S3)使用Logistic混沌映射算法进行置乱。
5.根据权利要求1所述的方法，其特征在于，步骤(S4)中，STC编码器根据嵌入点的修改代价，自适应的选取最优嵌入路径将秘密信息嵌入到比特串，使修改总失真最小，得到载密的比特串。
6.根据权利要求1所述的方法，其特征在于，采用以下步骤提取秘密信息：
(P1)根据步骤(S1)中哈夫曼码字映射表的构造方法构造哈夫曼码字映射表；
(P2)对音频进行部分解码，得到哈夫曼码字，根据哈夫曼码字映射表将每个可隐写码字映射成载密比特串；
(P3)使用步骤(S3)中的置乱算法与密钥对载密比特串进行置乱；
(P4)使用STC编码从载密比特串中提取秘密信息，提取过程使用的参数与嵌入过程中使用的参数保持一致；
(P5)重复步骤(P2)～步骤(P4)，直到音频解码结束或信息完整提取，得到还原的秘密信息。
7.一种用于权利要求1所述方法的基于等长熵码字替换的压缩音频自适应隐写系统，其特征在于，包括码字-二进制映射模块、失真函数模块、置乱模块、STC编码模块、逆置乱模块和二进制-码字映射模块；所述码字-二进制映射模块将原始哈夫曼码流转换为二进制比特串；所述失真函数模块根据心理声学模型和系数修改幅度计算每个可嵌入位置的修改代价；所述置乱模块将二进制比特串和修改代价使用相同的置乱算法和密钥进行置乱；所述STC编码模块利用STC编码结合可嵌入位置的修改代价，将待嵌信息嵌入到置乱后的比特串中，得到载密比特串；所述逆置乱模块将载密比特串采用相同的置乱算法和密钥进行逆置乱；所述二进制-码字映射模块将逆置乱后的载密比特串逆映射成原始的哈夫曼码。
8.根据权利要求7所述的系统，其特征在于，在提取秘密信息时，所述码字-二进制映射模块根据哈夫曼码字映射表将哈夫曼码字串映射成比特串S；所述置乱模块将比特串S进行置乱；所述STC编码模块利用STC解码从置乱后的比特串中提取秘密信息。
CN201811067226.8A 2018-09-13 2018-09-13 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统 Expired - Fee Related CN109346091B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811067226.8A CN109346091B (zh) 2018-09-13 2018-09-13 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811067226.8A CN109346091B (zh) 2018-09-13 2018-09-13 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统

Publications (2)

Publication Number Publication Date
CN109346091A CN109346091A (zh) 2019-02-15
CN109346091B true CN109346091B (zh) 2021-12-14

Family

ID=65305497

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811067226.8A Expired - Fee Related CN109346091B (zh) 2018-09-13 2018-09-13 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统

Country Status (1)

Country Link
CN (1) CN109346091B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6768980B1 (en) * 1999-09-03 2004-07-27 Thomas W. Meyer Method of and apparatus for high-bandwidth steganographic embedding of data in a series of digital signals or measurements such as taken from analog data streams or subsampled and/or transformed digital data
CN101546559A (zh) * 2009-04-13 2009-09-30 宁波大学 一种用于压缩数字音频的秘密信息密写及提取方法
CN105118512A (zh) * 2015-07-23 2015-12-02 武汉大学 一种面向aac数字音频的通用隐写分析方法
CN106228981A (zh) * 2016-08-03 2016-12-14 中山大学 一种基于压缩域的mp3自适应隐写方法
CN106448686A (zh) * 2016-08-10 2017-02-22 中山大学 一种自适应音频空域隐写方法
CN107623855A (zh) * 2016-07-13 2018-01-23 谭心瑶 一种基于压缩编码的高嵌入率数据隐写装置及隐写方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007007627A1 (de) * 2006-09-15 2008-03-27 Rwth Aachen Steganographie in digitalen Signal-Codierern

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6768980B1 (en) * 1999-09-03 2004-07-27 Thomas W. Meyer Method of and apparatus for high-bandwidth steganographic embedding of data in a series of digital signals or measurements such as taken from analog data streams or subsampled and/or transformed digital data
CN101546559A (zh) * 2009-04-13 2009-09-30 宁波大学 一种用于压缩数字音频的秘密信息密写及提取方法
CN105118512A (zh) * 2015-07-23 2015-12-02 武汉大学 一种面向aac数字音频的通用隐写分析方法
CN107623855A (zh) * 2016-07-13 2018-01-23 谭心瑶 一种基于压缩编码的高嵌入率数据隐写装置及隐写方法
CN106228981A (zh) * 2016-08-03 2016-12-14 中山大学 一种基于压缩域的mp3自适应隐写方法
CN106448686A (zh) * 2016-08-10 2017-02-22 中山大学 一种自适应音频空域隐写方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Steganalysis for MP3Stego using differential statistics of quantization step;Diqun Yan等;《Digital Signal Processing》;20130218;第23卷;全文 *
基于Huffman编码的大容量MP3隐写算法;严迪群等;《四川大学学报(自然科学版)》;20111130;第48卷(第6期);全文 *
基于码表分布特征和重编码的MP3Stego隐写分析;万威 等;《中国科学院研究生院学报》;20120131;第29卷(第1期);全文 *

Also Published As

Publication number Publication date
CN109346091A (zh) 2019-02-15

Similar Documents

Publication Publication Date Title
Ali et al. High capacity, transparent and secure audio steganography model based on fractal coding and chaotic map in temporal domain
He et al. Reversible data hiding in JPEG images based on negative influence models
Zhang et al. Dither modulation based adaptive steganography resisting JPEG compression and statistic detection
Kumar et al. Reversible data hiding: A contemporary survey of state-of-the-art, opportunities and challenges
Lee et al. Lossless data hiding for VQ indices based on neighboring correlation
Yang et al. Reversible watermarking in the VQ-compressed domain
CN111755018B (zh) 基于小波变换及量化嵌入密钥的音频隐藏方法及装置
Yin et al. Robust adaptive steganography based on dither modulation and modification with re-compression
Kumar et al. Performance comparison of robust steganography based on multiple transformation techniques
Shie et al. Data hiding based on compressed VQ indices of images
CN102664014B (zh) 一种基于对数量化索引调制的盲音频水印实现方法
Lee et al. Information hiding based on block match coding for vector quantization-compressed images
Wei et al. Lightweight AAC Audio Steganalysis Model Based on ResNeXt
Pan et al. Robust VQ-based digital watermarking for memoryless binary symmetric channel
CN109346091B (zh) 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统
Sun et al. Reversible data hiding for 3D meshes in the PVQ-compressed domain
Ramkumar et al. Signaling methods for multimedia steganography
Qin et al. Reversible data embedding for vector quantization compressed images using search‐order coding and index parity matching
Jiang et al. A high capacity steganographic method based on quantization table modification
Wu et al. Robust and blind audio watermarking scheme based on genetic algorithm in dual transform domain
Zheng et al. An adjustable RDH method for AMBTC-compressed images using unilateral pixel value modification strategy
CN104867496A (zh) 一种基于mdct量化系数的小值区的aac音频隐写和提取方法
Hemalatha A Robust MP3 Audio Steganography with Improved Capacity
Liu et al. A new SMVQ-based reversible data hiding scheme using state-codebook sorting
CN109448744B (zh) 一种基于符号位自适应嵌入的mp3音频信息隐藏方法和系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20211214