CN106228981A - 一种基于压缩域的mp3自适应隐写方法 - Google Patents
一种基于压缩域的mp3自适应隐写方法 Download PDFInfo
- Publication number
- CN106228981A CN106228981A CN201610634457.7A CN201610634457A CN106228981A CN 106228981 A CN106228981 A CN 106228981A CN 201610634457 A CN201610634457 A CN 201610634457A CN 106228981 A CN106228981 A CN 106228981A
- Authority
- CN
- China
- Prior art keywords
- granule
- steganography
- information
- embedding
- secret information
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000006978 adaptation Effects 0.000 title claims abstract description 15
- 230000006835 compression Effects 0.000 title claims abstract description 14
- 238000007906 compression Methods 0.000 title claims abstract description 14
- 238000013139 quantization Methods 0.000 claims abstract description 38
- 238000000605 extraction Methods 0.000 claims abstract description 19
- 239000008187 granular material Substances 0.000 claims description 156
- 239000002245 particle Substances 0.000 claims description 14
- 238000013461 design Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 abstract description 4
- 230000005236 sound signal Effects 0.000 abstract description 3
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 208000011580 syndromic disease Diseases 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000006837 decompression Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
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
-
- 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
- G10L15/00—Speech recognition
- G10L15/06—Creation of reference templates; Training of speech recognition systems, e.g. adaptation to the characteristics of the speaker's voice
- G10L15/063—Training
-
- 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
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Image Processing (AREA)
Abstract
本发明涉及多媒体信息安全领域,具体涉及一种基于压缩域的MP3自适应隐写方法,包括秘密信息的嵌入和秘密信息的提取两个步骤。本发明方法以原始音频信号经过子带滤波器转换及量化编码后的频域系数为隐写的基本单位,利用信息编码工具Syndrome Trellis Codes(STC)进行自适应嵌入信息并实现正确信息提取。经大量实验证明,本隐写方法与现有的基于频域系数修改的音频隐写方法相比,具有更高的嵌入容量及安全性。因此,本发明方法对隐秘通信等应用领域具有重要意义。
Description
技术领域
本发明涉及多媒体信息安全领域,更具体地,涉及一种基于压缩域的MP3自适应隐写方法。
背景技术
音频隐写是信息隐藏技术的一个重要分支。目前的基于压缩域的MP3音频隐写算法依据隐写操作发生的位置分为三种类型:量化编码前隐写,边量化边隐写及哈夫曼熵编码隐写。量化编码前隐写的方法面对的挑战是如何正确提取秘密信息,哈夫曼熵编码隐写方法会直接降低隐写音频的听觉质量,而边量化边隐写可以在保证隐写音频听觉质量、隐蔽性的同时实现完整秘密信息的正确提取。经典的MP3stego音频隐写工具采用的隐写方法即为边量化边隐写,但其每个颗粒仅可以嵌入1个比特,嵌入容量低且安全性不高,因而不足以满足实际应用中的信息嵌入要求。
发明内容
本发明为克服上述现有技术所述的不足,提供一种基于压缩域的MP3自适应隐写方法。该方法基于信息编码工具STC,选取原始音频信号经过MP3量化编码后的频域系数作为嵌入载体,通过各频域系数修改代价的设定自适应地修改嵌入载体的若干位以达到秘密信息嵌入的目的,并可以实现秘密信息的正确提取。其中,信息编码工具STC可以根据嵌入单元设定的修改代价,自适应选取合适的频域系数进行修改,以实现最小化所有频域系数的修改代价之和。嵌入单元设置的代价值越大,该频域系数被期望修改的几率就越低。本发明可嵌入的秘密信息容量大,隐蔽性高,且经实验证明,本发明方法可有效抵抗现有的基于压缩域的MP3音频隐写分析方法,安全性更高。
为了实现上述发明目的,本发明的技术方案如下:
一种基于压缩域的MP3自适应隐写方法,包括秘密信息的嵌入和秘密信息的提取两个步骤;秘密信息的嵌入操作发生在MPEG-1Layer3(MP3)音频编码中频域系数的量化编码过程中;
所述秘密信息的嵌入为边量化边嵌入,包括如下子步骤:
(11)预编码获取原始音频量化编码时可嵌入秘密信息的颗粒总数N;
(12)依据可隐写颗粒总数N和秘密信息长度L设定各项参数;
(13)判断当前颗粒是否进行隐写并选定嵌入载体;
(14)设置信息编码工具STC修改频域系数的代价;
(15)利用STC进行秘密信息的嵌入;
所述秘密信息的提取包括如下子步骤:
(21)解码下一颗粒参数,并判断该颗粒是否为可隐写颗粒;
(22)获取信息嵌入的各项参数;
(23)判断当前颗粒是否为最后一个隐写颗粒;
(24)判断当前颗粒是否有嵌入信息,有则提取信息。
优选地,Layer3每个编码帧包含两个颗粒,每个颗粒包含576个频域系数,被分为大值区、小值区及零值区,且有块类型参数包含长块、短块、起始块及终止块等值,其中所述步骤(11)中可嵌入秘密信息的颗粒需同时满足如下条件:①颗粒块类型不是短块;②颗粒大值区频域系数总数≥50。
优选地,所述步骤(12)依据可隐写颗粒总数N和秘密信息长度L设定各项参数,包括如下子步骤:
(121)依据可隐写颗粒数N及秘密信息文件的比特数L,得到每个颗粒需嵌入的平均比特数
(122)由于使用STC工具嵌入秘密信息时,针对每组频域系数,最低需嵌入的秘密信息比特数为10;若平均比特数AvgBits<10,设间隔参数 将除第1、2个可隐写颗粒外所有可隐写颗粒划分为若干个间隔,在每个间隔中随机选取一个颗粒进行秘密信息的嵌入,嵌入容量m=10;根据秘密信息的容量,会产生个间隔,可标记为{Interv1,Interv2,Interv3,…,Intervn}。由初始密钥key生成n个范围在[1,Interval]内的随机数{k1,k2,k3,…,kn},其中对随机数ki的随机种子为key+i。针对第i(1≤i≤n)个间隔Intervi,选择该间隔中第ki个颗粒进行秘密信息的嵌入;
(123)平均比特数AvgBits≥10,则设定间隔参数Interval为1,即每个可嵌入秘密信息的颗粒均需嵌入AvgBits个比特的信息,嵌入容量m=AvgBits;嵌入容量的选择用如下公式直观说明:
(124)第1、2个可隐写颗粒用于嵌入参数信息,最后一个嵌入秘密信息的颗粒从最后一个间隔Intervn中选择,标记为N’,其嵌入秘密信息的长度为finalBits,计算公式如下:
N’=2+(n-1)×Interval+kn
优选地,步骤(13)对颗粒逐个进行量化编码,直到所有颗粒量化结束,并判断当前颗粒是否进行隐写,包含如下子步骤:
(131)选定颗粒量化后的频域系数,获取频域系数编码的大值区及块类型参数,若块类型不为短块且大值区频域系数总数≥50,则当前颗粒为可隐写颗粒,可隐写颗粒计数count加1,进行下一步骤;否则对下一个颗粒进行量化,重复本步骤;其中,可隐写颗粒计数count的初始值为0;
(132)可隐写颗粒计数count=1、2时分别嵌入步骤(12)中的参数信息,可隐写颗粒N’嵌入秘密信息文件最后finalBits比特信息;
(133)计算当前颗粒所处的间隔序号i,若该颗粒在其所处的间隔Intervi中处于第ki个位置,其中随机数ki的随机种子为初始密钥key与序号i的和,即key+i,则进行下一个步骤,否则量化下一个颗粒并返回步骤(131);
(134)选定当前颗粒大值区的所有频域系数作为STC操作的嵌入载体,记为Cover。
优选地,步骤(14)中设置信息编码工具STC修改嵌入载体中各频域系数的代价,其具体过程为:量化编码操作对每个颗粒的576个频域系数进行整体量化编码;每组颗粒量化前的系数记为序列{XR1,XR2,XR3,…,XR576},量化后的系数记为序列{IX1,IX2,IX3,…,IX576}。为了能让STC可以根据频域系数组的值自适应地修改若干位,需要对每个频域系数设定修改的代价,记为costs;具体设计原则如下:
①对于第i位量化后的频域系数IXi,修改的代价包含三种代价:在原值上减1的代价、保持不变的代价、在原值上加1的代价;
②考虑到在频域系数IXi上的修改会影响到IXi所在的比例因子频带的实际失真,在设定修改代价之前需要评估频域系数IXi加减1对实际失真的影响,即将IXi加减1后的失真与IXi所在的比例因子频带sb的平均失真DT(sb)对比,其中量化后的每个频域系数IXi与量化前的系数XRi的失真计算公式以及比例因子频带sb的平均失真计算公式如下:
其中qstep表示量化编码步骤的量化步长,BW(sb)表示频带sb的宽度即该频带包含频域系数的个数,low(sb)表示频带sb下边界,high(sb)表示频带sb的上边界。
③此外,为了避免修改后的频域系数组大值区与小值区边界混淆,将大值区最后4位频域系数的修改代价值设为无穷大,即不对其进行修改。对于第i位量化后的频域系数IXi,按如下两种情况设定其修改代价costs:
情况1:IXi≠0and i≤length(Cover)–4
规定函数
情况2:IXi=0or i>length(Cover)–4
其中,length(Cover)表示嵌入载体Cover的长度,也即当前量化颗粒大值区频域系数的总数。
优选地,步骤(15)利用信息编码工具STC进行秘密信息的嵌入,具体是:依据设定好的各频域系数的修改代价和当前颗粒需嵌入的秘密信息,对选定的嵌入载体做STC操作,自适应地选取嵌入载体中的若干位进行修改以达到信息嵌入的目的,并实现最小化嵌入载体所有频域系数的修改代价之和;若成功嵌入则重新计算修改后的频域系数的各项参数,需保证对修改后的频域系数进行哈夫曼编码所需要的比特位数要小于当前颗粒可用于编码的最大位数,满足条件则对下一个颗粒进行秘密信息的嵌入,否则需要对当前颗粒重新进行信息嵌入;返回(13)步骤。
优选地,所述秘密信息的提取操作发生在MPEG-1Layer 3(MP3)音频解码中哈夫曼解码还原频域系数的过程中,具体子步骤为:
(210)解码还原下一个颗粒的块类型参数、大值区的参数及该颗粒的576个频域系数;若块类型不为短块且大值区频域系数总数≥50,则当前颗粒为可隐写颗粒,可隐写颗粒计数count加1,进行下一步骤;否则重复本步骤;其中可隐写颗粒计数count的初始值为0;
(220)若当前可隐写颗粒计数count=1、2,提取嵌入到第1、2个隐写颗粒中的参数:每个颗粒嵌入的秘密信息比特数m,初始密钥key,间隔Interval,终止颗粒N’,终止颗粒嵌入比特数finalBits;
(230)若当前可隐写颗粒计数count=N’,则将大值区内的所有频域系数当作提取载体,提取的秘密信息长度为finalBits,并利用STC工具对提取载体做解码操作得到该颗粒嵌入的finalBits比特秘密信息;整合解码所得的所有二进制流秘密信息即为最终的秘密信息文件,终止解码操作,否则进行下一步骤;
(240)从密钥key生成随机数组{k1,k2,k3,…,kn},对ki∈[1,Interval],且随机数ki的随机种子为key+i。若当前颗粒在其所处的第i个间隔Intervi中处于第ki个位置,则选取当前颗粒大值区内的所有频域系数作为提取载体,提取的秘密信息长度为m,并利用STC工具对提取载体进行解码得到该颗粒嵌入的m比特秘密信息。返回步骤(210)。
与现有技术相比,本发明技术方案的有益效果是:
1)利用MP3Stego、本发明提出的基于压缩域的MP3自适应隐写方法生成不同嵌入容量的隐写音频,其中MP3Stego采用最高嵌入容量即每个颗粒嵌入1bit,本发明描述的自适应隐写方法分别嵌入1bit、5bits、10bits。针对每一种情况,采用三种通用的隐写分析特征,一半样本作为训练集,另一半样本作为测试集,并利用优化参数的SVM分类器进行训练测试,本发明描述的隐写方法检测率较MP3Stego有大幅降低(接近50%),安全性更好。
2)利用MP3Stego、本发明提出的基于压缩域的MP3自适应隐写方法生成不同嵌入容量的隐写音频,其中MP3Stego采用最高嵌入容量即每个颗粒嵌入1bit,本发明描述的自适应隐写方法分别嵌入1bit、5bits、10bits。针对每一种情况,采用本发明描述的方法生成的隐写音频信噪比更高,音频质量更好。
附图说明
图1所示为采用本发明进行秘密信息嵌入、提取的算法流程示意图。
图2所示为频域系数的参数示意图。
图3所示为隐写颗粒选取示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。下面结合附图和实施例对本发明的技术方案做进一步的说明。
图1所示为采用本发明进行秘密信息嵌入、提取的算法流程示意图。
其中秘密信息的嵌入过程如图1(a)、1(b)所示,包括如下子步骤:(11)预编码获取原始音频量化编码时可嵌入秘密信息的颗粒总数N;(12)依据可隐写颗粒总数N和秘密信息长度L设定各项参数;(13)判断当前颗粒是否进行隐写并选定嵌入载体;(14)设置信息编码工具STC修改频域系数的代价;(15)利用STC进行秘密信息的嵌入。
每个步骤具体说明如下:
(11)预编码获取原始音频可嵌入秘密信息的颗粒总数N。Layer3每个编码帧包含两个颗粒,量化编码后每个颗粒包含576个频域系数,被分为大值区、小值区及零值区,且有块类型参数取值分别为长块、短块、起始块和终止块,如图2所示。其中可嵌入秘密信息的颗粒需同时满足如下条件:①颗粒块类型不是短块;②颗粒大值区频域系数总数≥50。
(12)设定嵌入秘密信息时所需参数,包括如下子步骤:
①依据可隐写颗粒数N及秘密信息文件的比特数L,可得每个颗粒需嵌入的平均比特数
②由于使用STC工具嵌入秘密信息时,针对每组频域系数,最低需嵌入的秘密信息比特数为10,因此需对平均比特数AvgBits<10的情况设间隔参数将除第1、2个可隐写颗粒外所有可隐写颗粒划分为若干个间隔,每个间隔中随机选取一个颗粒嵌入信息,嵌入容量m=10;
根据秘密信息的容量,会产生个间隔,如图3所示,表示为{Interv1,Interv2,Interv3,…,Intervn}。由初始密钥key生成n个范围在[1,Interval]内的随机数{k1,k2,k3,…,kn},每个随机数对应相应间隔内随机颗粒的选择。生成随机数方法如图3所示,每个随机数的随机种子设定为key+i,因此每个初始key都对应一组特定的随机数序列。针对第i(1≤i≤n)个间隔,选择该间隔中第ki个颗粒进行秘密信息的嵌入,过程如图3所示。
③若平均比特数AvgBits≥10,则Interval为1,即每个可嵌入秘密信息的颗粒均需嵌入AvgBits个比特的信息,嵌入容量m=AvgBits。嵌入容量的选择可用如下公式直观说明:
④第1、2个可隐写颗粒用于嵌入参数信息,最后一个嵌入秘密信息的颗粒从最后一个间隔Intervn中选择,标记为N’,其嵌入秘密信息的长度为finalBits,计算公式如下:
N’=2+(n-1)×Interval+kn
(13)逐个对颗粒进行量化编码,直到所有颗粒量化结束,并判断当前颗粒是否进行隐写,包含如下子步骤:
①选定颗粒量化后的频域系数,获取频域系数编码的大值区及块类型参数,若大值区频域系数总数≥50且块类型不为短块,则当前颗粒为可隐写颗粒,可隐写颗粒计数count加1,进行下一步骤;否则对下一个颗粒进行量化,重复本步骤;其中可隐写颗粒计数count的初始值为0;
②可隐写颗粒计数count=1、2时分别嵌入步骤(12)中的参数信息,可隐写颗粒N’嵌入秘密信息最后finalBits比特信息;
③计算当前颗粒所处的间隔序号i,若该颗粒在其所处的第i个间隔Intervi中处于第ki个位置,其中ki是以key+i为随机种子生成的随机数,则进行下一个步骤,否则量化下一个颗粒并返回步骤①。
④选定当前颗粒的所有大值区的频域系数作为STC操作的嵌入载体,记为Cover。
(14)设置STC修改嵌入载体中各频域系数的代价:量化编码操作对每个颗粒的576个频域系数进行整体量化编码。每组颗粒量化前的系数由序列{XR1,XR2,XR3,…,XR576}表示,量化后的系数记为{IX1,IX2,IX3,…,IX576}。为了能让STC可以根据频域系数组的值自适应地修改若干位,需要对每个频域系数设定修改的代价,记为costs。具体设计原则如下:
①对于第i位量化后的频域系数IXi,修改的代价包含三种代价:在原值上减1的代价、保持不变的代价、在原值上加1的代价;
②考虑到在频域系数IXi上的修改会影响到IXi所在的比例因子频带的实际失真,在设定修改代价之前需要评估频域系数IXi加减1对实际失真的影响,即将IXi加减1后的失真与IXi所在的比例因子频带sb的平均失真DT(sb)对比,其中量化后的每个频域系数IXi与量化前的系数XRi的失真计算公式以及比例因子频带sb的平均失真计算公式如下:
其中qstep表示量化编码步骤的量化步长,BW(sb)表示频带sb的宽度即该频带包含频域系数的个数,low(sb)表示频带sb下边界,high(sb)表示频带sb的上边界。
③此外,为了避免修改后的频域系数大值区与小值区边界混淆,将大值区最后4位频域系数的修改代价值设为无穷大,即不对其进行修改。对于第i位量化后的频域系数IXi,按如下两种情况设定其修改代价costs:
情况1:IXi≠0and i≤length(Cover)–4
规定函数
情况2:IXi=0or i>length(Cover)–4
其中,length(Cover)表示嵌入载体Cover的长度,也即当前量化颗粒大值区频域系数的总数。
(15)依据设定好的各频域系数的修改代价和当前颗粒需嵌入的秘密信息,对嵌入载体做STC操作,选取嵌入载体中的若干位进行修改以达到信息嵌入的目的,并实现最小化嵌入载体各频域系数的修改代价之和;若成功嵌入则重新计算修改后的频域系数组的各项参数,需保证对修改后的频域系数进行哈夫曼编码所需要的比特位数要小于当前颗粒可用于编码的最大位数,满足条件则对下一个颗粒进行秘密信息的嵌入,否则需要对当前颗粒重新进行信息嵌入;返回(13)步骤。
本发明秘密信息的提取操作发生在MPEG-1Layer 3(MP3)音频解码中哈夫曼解码还原频域系数的过程中,如图1(c)所示,包括步骤:(21)解码下一颗粒参数,判断该颗粒是否为可隐写颗粒;(22)获取信息嵌入的各项参数;(23)判断当前颗粒是否为最后一个隐写颗粒;(24)判断当前颗粒是否有嵌入信息,有则提取信息。
(21)解码还原下一个颗粒的块类型参数,大值区的参数及该颗粒的576个频域系数;若大值区频域系数总数≥50且块类型不为短块,则当前颗粒为可隐写颗粒,可隐写颗粒计数count加1,进行下一步骤;否则重复本步骤;其中可隐写颗粒计数count的初始值为0;
(22)若当前可隐写颗粒计数count=1、2,提取嵌入到第1、2个可隐写颗粒中的参数:每个颗粒嵌入的秘密信息比特数m,密钥key,间隔Interval,终止颗粒N’,终止颗粒嵌入比特数finalBits;
(23)若当前可隐写颗粒计数count=N’,将大值区内的所有频域系数当作提取载体,提取的秘密信息长度为finalBits,并利用STC工具对提取载体做解码操作得到该颗粒嵌入的finalBits比特秘密信息;整合解码所得的所有二进制流秘密信息即为最终的秘密信息文件,终止解码操作,否则进行下一步骤;
(24)从密钥key生成随机数组{k1,k2,k3,…,kn},对ki∈[1,Interval|,且随机数ki的随机种子为key+i。若当前颗粒在其所处的第i个间隔Intervi中处于第ki个位置,则选取当前颗粒大值区内的所有频域系数作为提取载体,提取的秘密信息长度为m,并利用STC工具对提取载体进行解码得到该颗粒嵌入的m比特秘密信息。返回步骤(21)。
下面将结合具体的实验数据对本发明做更进一步的分析说明:
1、音频集合:
本实施例中,原始音频信号采用无损的WAV格式音频;利用5000个单声道长度为20秒的无损WAV音频做测试,采用8Hz编码器对原始WAV音频进行MP3压缩编码,并在编码过程中进行隐写操作。
2、隐写工具及方法
MP3Stego工具
本发明提出的基于压缩域的MP3自适应隐写方法
3、通用型隐写分析特征
“M.Y.Qiao,A.H.Sung,and Q.Z.Liu.MP3audio steganalysis.InformationScience.231(2013),123-134.”论文提出Qiao294-D隐写分析特征;
“D.Q.Yan,R.D.Wang,X.M.Yu,and J.zhu.Steganalysis for MP3stego usingdifferential statistics of quantization step.Digitall Signal Processing.23(2013)1181-1185.”论文提出Yan2-D隐写分析特征;
“C.Jin,R.D.Wang,D.Q.Yan.Stegnalysis of MP3Stego with embedding-rateusing Markov feature.Multimed.Tools.Appl.DOI,10,1007/s 11042-016-3264-y.”论文提出Jin338-D隐写分析特征;
4、安全性能分析:
利用MP3Stego、本发明提出的基于压缩域的MP3自适应隐写方法生成不同嵌入容量的隐写音频,其中MP3Stego采用最高嵌入容量即每个颗粒嵌入1bit,本发明描述的自适应隐写方法分别嵌入1bit、5bits、10bits。针对每一种情况,采用三种通用的隐写分析特征,一半的样本用作训练集,另一半的样本用作测试集,采用优化参数的SVM分类器进行训练测试,下表1列出了各种情况下的检测率。可以看出本发明描述的自适应隐写方法检测率较MP3Stego有大幅降低(接近50%),安全性更好。
表1
5、音频质量分析:
利用MP3Stego、本发明提出的基于压缩域的MP3自适应隐写方法生成不同嵌入容量的隐写音频,其中MP3Stego采用最高嵌入容量即每个颗粒嵌入1bit,本发明描述的自适应隐写方法分别嵌入1bit、5bits、10bits。针对每一种情况,计算隐写MP3音频和未进行秘密信息嵌入的MP3音频解压回空域后与原始WAV音频的信噪比,单位dB。
下表2列出了各种情况下的信噪比。可以看出采用本发明描述的自适应隐写方法隐写的音频信噪比数值有显著的提升,音频质量更好。
表2
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (7)
1.一种基于压缩域的MP3自适应隐写方法,包括秘密信息的嵌入和秘密信息的提取两个步骤;秘密信息的嵌入操作发生在MPEG-1 Layer 3(MP3)音频编码中频域系数的量化编码过程中,其特征在于,
所述秘密信息的嵌入为边量化边嵌入,包括如下子步骤:
(11)预编码获取原始音频量化编码时可嵌入秘密信息的颗粒总数N;
(12)依据可隐写颗粒总数N和秘密信息长度L设定各项参数;
(13)判断当前颗粒是否进行隐写并选定嵌入载体;
(14)设置信息编码工具STC修改频域系数的代价;
(15)利用STC进行秘密信息的嵌入;
所述秘密信息的提取包括如下子步骤:
(21)解码下一颗粒参数,判断当前颗粒是否为可隐写颗粒;
(22)获取信息嵌入的各项参数;
(23)判断当前颗粒是否为最后一个隐写颗粒;
(24)判断当前颗粒是否有嵌入信息,有则提取信息。
2.根据权利要求1所述的方法,其特征在于,Layer3每个编码帧包含两个颗粒,每个颗粒包含576个频域系数,被分为大值区、小值区及零值区,且有块类型参数取值分别为长块、短块、起始块和终止块,其中所述步骤(11)中可嵌入秘密信息的颗粒需同时满足如下条件:①颗粒块类型不是短块;②颗粒大值区频域系数总数≥50。
3.根据权利要求2所述的方法,其特征在于,所述步骤(12)依据可隐写颗粒总数N和秘密信息长度L设定各项参数,包括如下子步骤:
(121)依据可隐写颗粒数N及秘密信息文件的比特数L,得到每个颗粒需嵌入的平均比特数
(122)平均比特数AvgBits<10时,设间隔参数将音频除第1、2个可隐写颗粒外的所有颗粒划分为若干个间隔,每个间隔包含Interval个颗粒,从中随机选取一个颗粒进行秘密信息的嵌入,嵌入容量m=10;根据秘密信息的容量,会产生个间隔,表示为{Interv1,Interv2,Interv3,…,Intervn};由初始密钥key生成n个范围在[1,Interval]内的随机数{k1,k2,k3,…,kn},其中对随机数ki的随机种子为key+i;每个随机数对应一个间隔内随机颗粒的选择;针对第i个间隔Intervi,1≤i≤n;选择该间隔中第ki个颗粒进行秘密信息的嵌入;
(123)平均比特数AvgBits≥10,则设定间隔参数Interval为1,即每个可嵌入秘密信息的颗粒均需嵌入AvgBits个比特的信息,嵌入容量m=AvgBits;嵌入容量的选择用如下公式直观说明:
(124)第1、2个可隐写颗粒用于嵌入参数信息;最后一个嵌入秘密信息的颗粒从最后一个间隔Intervn中选择,标记为N’,其嵌入秘密信息的长度为finalBits,计算公式如下:
N’=2+(n-1)×Interval+kn
4.根据权利要求3所述的方法,其特征在于,步骤(13)是依次对颗粒进行量化编码,直到所有颗粒量化结束,并判断当前颗粒是否进行隐写,包含如下子步骤:
(131)选定颗粒量化后的频域系数,获取频域系数编码的大值区及块类型参数,若块类型不为短块且大值区频域系数总数≥50,则当前颗粒为可隐写颗粒,可隐写颗粒计数count加1,进行下一步骤;否则对下一个颗粒进行量化并重复本步骤;其中,可隐写颗粒计数count的初始值为0;
(132)可隐写颗粒计数count=1、2时分别嵌入步骤(12)中的参数信息,可隐写颗粒N’嵌入秘密信息文件最后finalBits比特信息;
(133)计算当前颗粒所处的间隔序号i,若该颗粒在其所处的间隔Intervi中处于第ki个位置,其中随机数ki的随机种子为初始密钥key与序号i的和,即key+i,则进行下一个步骤,否则量化下一个颗粒并返回步骤(131);
(134)选定当前颗粒大值区内的所有频域系数作为STC操作的嵌入载体,记为Cover。
5.根据权利要求4所述的方法,其特征在于,步骤(14)中设置信息编码工具STC修改嵌入载体中各频域系数的代价,其具体过程为:量化编码操作对每个颗粒的576个频域系数进行整体量化编码;每组颗粒量化前的系数记为{XR1,XR2,XR3,…,XR576},量化后的系数记为{IX1,IX2,IX3,…,IX576};为了能让STC可以根据频域系数组的值自适应地修改若干位,需要对每个频域系数设定修改的代价,记为costs;具体设计原则如下:
①对于第i位量化后的频域系数IXi,修改的代价包含三种代价:在原值上减1的代价、保持不变的代价、在原值上加1的代价;
②考虑到在频域系数IXi上的修改会影响到IXi所在的比例因子频带的实际失真,在设定各频域系数的修改代价之前需要评估频域系数IXi加减1对实际失真的影响,即将IXi加减1后的失真与IXi所在的比例因子频带sb的平均失真DT(sb)对比,其中量化后的每个频域系数IXi与量化前的系数XRi的失真计算公式以及比例因子频带sb的平均失真计算公式如下:
其中qstep指量化编码步骤的量化步长,BW(sb)表示频带sb的宽度即该频带包含频域系数的个数,low(sb)表示频带sb下边界,high(sb)表示频带sb的上边界;
③对于第i位量化后的频域系数IXi,按如下两种情况设定其修改代价costs:
情况1:IXi≠0and i≤length(Cover)–4
规定函数
情况2:IXi=0or i>length(Cover)–4
其中,length(Cover)表示嵌入载体Cover的长度,也即当前量化颗粒大值区频域系数的总数。
6.根据权利要求5所述的方法,其特征在于,步骤(15)利用信息编码工具STC进行秘密信息的嵌入,具体是:依据设定好的各频域系数的修改代价和当前颗粒需嵌入的秘密信息,对选定的嵌入载体做STC操作,选取嵌入载体中的若干位进行修改以达到信息嵌入的目的,并实现最小化嵌入载体所有频域系数的修改代价之和;若成功嵌入则计算修改后的频域系数的各项参数,需保证对修改后的频域系数进行哈夫曼编码所需要的比特位数要小于当前颗粒可用于编码的最大位数,满足条件则对下一个颗粒进行秘密信息的嵌入,否则需要对当前颗粒重新进行信息嵌入;返回(13)步骤。
7.根据权利要求6所述的方法,其特征在于,所述秘密信息的提取操作发生在MPEG-1Layer 3(MP3)音频解码中哈夫曼解码还原频域系数的过程中,具体子步骤为:
(210)解码还原下一个颗粒的块类型参数、大值区参数及该颗粒的576个频域系数;若块类型不为短块且大值区频域系数总数≥50,则当前颗粒为可隐写颗粒,可隐写颗粒计数count加1,进行下一步骤;否则重复本步骤;其中可隐写颗粒计数count的初始值为0;
(220)若当前可隐写颗粒计数count=1、2,提取嵌入到第1、2个可隐写颗粒中的参数:每个颗粒嵌入的秘密信息比特数m,密钥key,间隔Interval,终止颗粒N’,终止颗粒嵌入的比特数finalBits;
(230)若当前可隐写颗粒计数count=N’,将当前颗粒大值区内的所有频域系数当作提取载体,提取的秘密信息长度为finalBits,并利用STC工具对提取载体进行解码得到嵌入的finalBits比特秘密信息;整合解码所得的所有二进制流秘密信息即为最终的秘密信息文件,终止解码操作,否则进行下一步骤;
(240)从密钥key生成随机数组{k1,k2,k3,…,kn},对ki∈[1,Interval],且随机数ki的随机种子为key+i;计算当前颗粒所处的间隔序号i,若该颗粒在其所处的间隔Intervi中处于第ki个位置,则将当前颗粒大值区内所有频域系数当作提取载体,提取信息长度为m,并利用STC工具对提取载体进行解码操作得到该颗粒嵌入的m比特秘密信息;返回步骤(210)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610634457.7A CN106228981B (zh) | 2016-08-03 | 2016-08-03 | 一种基于压缩域的mp3自适应隐写方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610634457.7A CN106228981B (zh) | 2016-08-03 | 2016-08-03 | 一种基于压缩域的mp3自适应隐写方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106228981A true CN106228981A (zh) | 2016-12-14 |
CN106228981B CN106228981B (zh) | 2019-07-19 |
Family
ID=57547098
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610634457.7A Expired - Fee Related CN106228981B (zh) | 2016-08-03 | 2016-08-03 | 一种基于压缩域的mp3自适应隐写方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106228981B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107886962A (zh) * | 2017-11-17 | 2018-04-06 | 南京理工大学 | 一种面向ip语音的高安全性隐写方法 |
CN108447492A (zh) * | 2018-03-13 | 2018-08-24 | 中山大学 | 一种基于分段的自适应直方图平移可逆信息隐藏方法 |
CN109346091A (zh) * | 2018-09-13 | 2019-02-15 | 中国科学院信息工程研究所 | 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统 |
CN109448744A (zh) * | 2018-12-14 | 2019-03-08 | 中国科学院信息工程研究所 | 一种基于符号位自适应嵌入的mp3音频信息隐藏方法和系统 |
CN110120228A (zh) * | 2019-04-28 | 2019-08-13 | 武汉大学 | 基于声谱图及深度残差网络的音频通用隐写分析方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1394409A (zh) * | 2000-11-06 | 2003-01-29 | 索尼公司 | 加密设备和方法、解密设备和方法、及存储介质 |
CN101527139A (zh) * | 2009-02-16 | 2009-09-09 | 成都九洲电子信息系统有限责任公司 | 一种音频编码解码方法及其装置 |
CN103886864A (zh) * | 2014-03-03 | 2014-06-25 | 南京邮电大学 | 一种对dwt域失真补偿量化嵌入的隐秘信息检测方法 |
CN103914839A (zh) * | 2014-03-27 | 2014-07-09 | 中山大学 | 一种基于隐写分析的图像拼接篡改检测方法及装置 |
US20150039320A1 (en) * | 2013-07-31 | 2015-02-05 | The Nielsen Company (Us), Llc | Systems Apparatus and Methods for Encoding/Decoding Persistent Universal Media Codes to Encoded Audio |
WO2015017679A1 (en) * | 2013-07-31 | 2015-02-05 | The Nielsen Company (Us), Llc | Apparatus, system and method for merging code layers for audio encoding and decoding |
CN104821169A (zh) * | 2013-04-08 | 2015-08-05 | 宁波大学 | 一种针对MP3Stegz的隐写检测方法 |
-
2016
- 2016-08-03 CN CN201610634457.7A patent/CN106228981B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1394409A (zh) * | 2000-11-06 | 2003-01-29 | 索尼公司 | 加密设备和方法、解密设备和方法、及存储介质 |
CN101527139A (zh) * | 2009-02-16 | 2009-09-09 | 成都九洲电子信息系统有限责任公司 | 一种音频编码解码方法及其装置 |
CN104821169A (zh) * | 2013-04-08 | 2015-08-05 | 宁波大学 | 一种针对MP3Stegz的隐写检测方法 |
US20150039320A1 (en) * | 2013-07-31 | 2015-02-05 | The Nielsen Company (Us), Llc | Systems Apparatus and Methods for Encoding/Decoding Persistent Universal Media Codes to Encoded Audio |
WO2015017679A1 (en) * | 2013-07-31 | 2015-02-05 | The Nielsen Company (Us), Llc | Apparatus, system and method for merging code layers for audio encoding and decoding |
CN103886864A (zh) * | 2014-03-03 | 2014-06-25 | 南京邮电大学 | 一种对dwt域失真补偿量化嵌入的隐秘信息检测方法 |
CN103914839A (zh) * | 2014-03-27 | 2014-07-09 | 中山大学 | 一种基于隐写分析的图像拼接篡改检测方法及装置 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107886962A (zh) * | 2017-11-17 | 2018-04-06 | 南京理工大学 | 一种面向ip语音的高安全性隐写方法 |
CN108447492A (zh) * | 2018-03-13 | 2018-08-24 | 中山大学 | 一种基于分段的自适应直方图平移可逆信息隐藏方法 |
CN108447492B (zh) * | 2018-03-13 | 2021-08-20 | 中山大学 | 一种基于分段的自适应直方图平移可逆信息隐藏方法 |
CN109346091A (zh) * | 2018-09-13 | 2019-02-15 | 中国科学院信息工程研究所 | 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统 |
CN109346091B (zh) * | 2018-09-13 | 2021-12-14 | 中国科学院信息工程研究所 | 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统 |
CN109448744A (zh) * | 2018-12-14 | 2019-03-08 | 中国科学院信息工程研究所 | 一种基于符号位自适应嵌入的mp3音频信息隐藏方法和系统 |
CN109448744B (zh) * | 2018-12-14 | 2022-02-01 | 中国科学院信息工程研究所 | 一种基于符号位自适应嵌入的mp3音频信息隐藏方法和系统 |
CN110120228A (zh) * | 2019-04-28 | 2019-08-13 | 武汉大学 | 基于声谱图及深度残差网络的音频通用隐写分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106228981B (zh) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106228981A (zh) | 一种基于压缩域的mp3自适应隐写方法 | |
CN1659785B (zh) | 信号多速率点阵矢量量化的方法和系统 | |
Tao et al. | Significantly improving lossy compression for scientific data sets based on multidimensional prediction and error-controlled quantization | |
EP2054883B1 (en) | Audio encoding system | |
CN105071815B (zh) | 算术编码设备和算术解码设备 | |
CN102132494B (zh) | 通信方法和通信装置 | |
CN103858433A (zh) | 分层熵编码及解码 | |
CN102687403A (zh) | 使用算术阶段压缩未充分利用的代码空间的编码器和解码器 | |
CN107578452B (zh) | 一种具有兼容格式和不变尺寸的jpeg图像加密方法 | |
CN103985389B (zh) | 一种针对amr音频文件的隐写分析方法 | |
CN104064193A (zh) | 一种线性预测语音编码的信息隐藏方法与提取方法 | |
CN103929210A (zh) | 一种基于遗传算法与神经网络的硬判决译码方法 | |
CN100489965C (zh) | 音频编码系统 | |
CN102057425A (zh) | 用于数字信号的压缩或解压缩的方法、系统和设备 | |
US20080292028A1 (en) | Method and Apparatus for Signal Processing and Encoding and Decoding Method, and Apparatus Therefor | |
CN101290774B (zh) | 音频编码和解码系统 | |
CN114328898A (zh) | 文本摘要生成方法及其装置、设备、介质、产品 | |
US20160111104A1 (en) | Signal encoding and decoding methods and devices | |
CN101908943B (zh) | 一种应用于td-scdma系统的信道解码方法 | |
CN105981100A (zh) | 用于改善对声场的高阶高保真度立体声响复制表示进行编码所需的边信息的编码的方法和装置 | |
CN109448744B (zh) | 一种基于符号位自适应嵌入的mp3音频信息隐藏方法和系统 | |
KR20120009837A (ko) | 컨텍스트 기반의 무손실 부호화 장치 및 방법, 그리고 복호화 장치 및 방법 | |
Kostadinov et al. | On digital watermarking for audio signals | |
CN105632505B (zh) | 主成分分析pca映射模型的编解码方法及装置 | |
KR102420763B1 (ko) | 뉴럴 네트워크 시스템 및 뉴럴 네트워크의 필터 데이터의 처리 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |
Granted publication date: 20190719 |