CN106971366A - 一种在音频信号中加入及提取水印的方法 - Google Patents
一种在音频信号中加入及提取水印的方法 Download PDFInfo
- Publication number
- CN106971366A CN106971366A CN201710068630.6A CN201710068630A CN106971366A CN 106971366 A CN106971366 A CN 106971366A CN 201710068630 A CN201710068630 A CN 201710068630A CN 106971366 A CN106971366 A CN 106971366A
- Authority
- CN
- China
- Prior art keywords
- lens
- light
- audio signal
- watermark information
- dfd
- 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
- 230000005236 sound signal Effects 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000006243 chemical reaction Methods 0.000 claims abstract description 22
- 238000005315 distribution function Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 238000000844 transformation Methods 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 abstract description 19
- 238000003384 imaging method Methods 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 6
- 239000000203 mixture Substances 0.000 description 4
- 230000004224 protection Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/0021—Image watermarking
-
- 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
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/06—Transformation of speech into a non-audible representation, e.g. speech visualisation or speech processing for tactile aids
- G10L21/10—Transforming into visible information
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Computational Linguistics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Quality & Reliability (AREA)
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
Abstract
本发明公开一种在音频信号中加入及提取水印的方法。通过将一维的音频信号转换成二维的音频图像,生成随机模板,基于光学衍射成像理论对音频图像、水印信息和随机模板进行衍射光学变换,然后进行加权叠加,实现了水印信息在音频图像中的信息隐藏。利用光学衍射成像具有的多维、大容量等优点,实现了大容量水印信息在音频信息中的隐藏,解决了现有方法水印信息嵌入量很小,信息隐藏受到很大限制的问题;利用光学衍射成像可撕毁性特点,实现了利用部分残缺音频图像提取完整水印信息的功能;通过提取水印信息,实现了音频信息来源可靠性的认证与证明;在提取水印时不需要原始载体图像,实现了水印信息的盲提取。
Description
技术领域
本发明属于信息安全领域,具体涉及一种在音频信号中加入及提取水印的方法。
背景技术
随着数字化多媒体信息应用及现代设备的普及,语音(音频信号)获得已变得越来越容易,并逐渐被普通百姓所掌握。语音还有其重要的真实和难以仿造特性,因此语音的应用范围在逐步扩大,例如,录音资料已成为司法认证的重要手段,也逐渐成为司法判断的重要依据。2002年出台的《最高人民法院关于民事诉讼证据的若干规定》的第七十条规定:“一方当事人提出的下列证据,对方当事人提出异议但没有足以反驳的相反证据的,人民法院应当确认其证明力:……(三)有其他证据佐证并以合法手段取得的、无疑点的视听资料或者与视听资料核对无误的复制件;……”但该规定第六十八条同时也明确了“以侵害他人合法权益或者违反法律禁止性规定的方法取得的证据,不能作为认定案件事实的依据。”这个规定说明,录音资料可以作为有效证据使用。这就从法律方面确定了录音资料有效性。然而,如何佐证视听(录音)资料是以合法手段取得的、无疑点的视听资料呢?比如,无签字能力的人的口头遗嘱,其合法性不成问题,但如何保证其未来不被篡改,成为一份“无疑点的视听资料”,这就需要进行音频保护,确保其真实性与完整性。
信息隐藏与认证技术是目前版权保护的有效手段,可以有效遏制数字化多媒体的非法使用。但在音频的传输保密与认证方面还相对落后,音频保护集中于音频的传输保密,有序列密码算法、分组密码算法及混沌密码算法,这些算法都是基于密码学实现的,仅能在数据从发送到接收的过程中进行数据保密,因此在此基础上又出现了音频信息隐藏技术,它将机密信息秘密地隐藏于一般的音频信息中。这是一种借助于现有图像信息隐藏技术发展而来的,这些方法隐藏的信息往往较少,非常有限,但目前这些算法都没有考虑信源的认证与公证问题。
发明内容
为了解决现有技术中存在的上述问题,本发明提出一种在音频信号中加入及提取水印的方法,将一维音频信号转换为二维音频图像,基于光学衍射成像理论经光学变换将水印信息嵌入二维音频图像中,第三方通过提取音频信号和水印信息实现音频信号的安全传输与信源认证。
为实现上述目的,本发明采用如下技术方案:
本发明提供一种在音频信号中加入水印的方法,包括以下步骤:
从一维的音频信号中截取一段,将截取的音频信号转换成二维的音频图像S,并保存为图像格式;
生成随机模板R,并保存为图像格式;
分别计算将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt;
对St、Wt和Rt进行加权叠加,将水印信息嵌入音频图像中,得到嵌入水印信息的音频图像Hw,用公式表示为:
Hw=w1St+w2Wt+w3Rt
式中,w1、w2、w3分别为St、Wt和Rt的加权系数;
定义DFD变换:当物平面x0o0y0的物光分布U0(x0,y0),物平面x0o0y0和观测平面xoy之间(没有透镜)的距离为d时观测平面的物光波场分布DFD[U0(x0,y0)],称DFD[U0(x0,y0)]为DFD变换;通过计算在d=z′1处得到音频图像的复振幅分布,z′1满足f为透镜的焦距;
将Sw、R以及用作解密密钥的参数传送给用户或第三方,所述参数包括光的波长λ、透镜的焦距f以及z1、z2和z3。
进一步地,所述将截取的音频信号转换成二维的音频图像S,包括:
将截取的音频信号归一化:
F(i)=f(i)/am
am=max[abs(f(i))]
式中,f(i)为截取的音频信号,i=1,2,…N,N为数据点数,F(i)为归一化后的音频信号。
将F(i)转换成不小于0的信号G(i):
G(i)=[F(i)+1]/2
将G(i)转化成二维矩阵,并保存为图像格式。
进一步地,所述分别计算将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt,包括:
计算透镜前d0处的物平面x0o0y0的物光分布为U0(x0,y0)时,透镜后d处的观测平面xoy的物光波场分布U(x,y):
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程;
将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt分别为:
式中,S0、W0、R0分别为音频图像S、水印图像W和随机模板R的像素强度分布函数。
进一步地,所述DFD变换的表达式为:
式中,j为虚数单位,k=2π/λ,π为圆周率,L为光波沿光轴的光程。
本发明还提供一种从嵌入水印信息的音频图像中提取水印信息的方法,包括以下步骤:
用户或第三方接收音频图像的复振幅分布Sw、随机模板R以及用作解密密钥的参数,所述参数包括光的波长λ、透镜的焦距f以及音频图像S、水印图像W和随机模板R置于透镜前的距离z1、z2和z3;
对Sw取模平方得到音频图像的强度分布|Sw|2;
将二维的音频图像的强度分布|Sw|2转化为一维的音频信号;
计算当收到的随机模板R位于透镜前z3时,在透镜后z′1处的物光场分布R′,并计算S′w=SW-R′;
定义DFD变换:当物平面x0o0y0的物光分布U0(x0,y0),物平面x0o0y0和观测平面xoy之间的距离为d时观测平面的物光波场分布DFD[U0(x0,y0)],称DFD[U0(x0,y0)]为DFD变换;
对S′w进行DFD逆变换得:
对Q进行DFD变换得到水印信息W′:
进一步地,所述将二维的音频图像的强度分布|Sw|2转化为一维的音频信号,包括:
将二维的|Sw|2转化为一维数据f′(i),i=1,2,…N,N为数据点数;
对f′(i)进行归一化:
F′(i)=f′(i)/a′m
a′m=max[f′(i)]
恢复负值:
G′(i)=2F′(i)-1.0
乘一增益系数km得到音频信号:
f″(i)=kmG′(i)。
进一步地,所述计算当收到的随机模板R位于透镜前z3时在透镜后z′1处的物光场分布R′,包括:
计算透镜前d0处的物平面x0o0y0的物光分布为U0(x0,y0)时,透镜后d处的观测平面xoy的物光波场分布U(x,y):
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程;
当收到的随机模板R位于透镜前z3时在透镜后z′1处的物光场分布R′为:
式中,R0为随机模板R的像素强度分布函数。
进一步地,所述DFD变换的表达式为:
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程。
与现有技术相比,本发明具有以下有益效果:
本发明通过将一维的音频信号转换成二维的音频图像,选择语音遗嘱人照片、录音音频的语言文字形式、公证处的完整文字或签章等作为水印信息,基于光学衍射成像理论对音频图像和水印信息进行衍射光学变换,然后进行加权叠加,实现了水印信息在音频图像中的信息隐藏。利用光学衍射成像具有的多维、大容量等优点,实现了大容量水印信息在音频信息中的隐藏,解决了现有方法水印信息嵌入量很小,信息隐藏受到很大限制的问题;利用光学衍射成像可撕毁性特点,实现了利用部分残缺音频图像提取完整水印信息的功能(利用语音图像残片可恢复完整的语音文字);通过提取水印信息,实现了音频信息来源可靠性的认证与证明,如水印中的口授人照片、公证处的完整文字、签章等信息可以充分显示与证明信息的来源及其真伪;在提取水印时不需要原始载体图像,实现了水印信息的盲提取。
附图说明
图1为由物平面、透镜、观测平面组成的光学衍射成像系统示意图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明实施例一种在音频信号中加入水印的方法,包括以下步骤:
步骤101,从一维的音频信号中截取一段,将截取的音频信号转换成二维的音频图像S,并保存为图像格式;
在本步骤中,可以利用格式工厂或MP3截取助手等音频截取软件,对输入的一维的音频信号进行截取。
步骤102,生成随机模板R,并保存为图像格式;
在本步骤中,可以利用MATLAB内置的rand(n,m)生成二维n×m大小矩阵,从而得到随机模板R;利用MATLAB内置函数Imwrite将二维矩阵保存为图像格式,即JPEG格式。
步骤103,如图1所示,分别计算将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt;
在本步骤中,随机模板R同音频图像S和水印图像W一起进行衍射光场分布计算,并进行后面步骤的加权叠加,可以防止当音频图像S和水印图像W距离较近时产生的相互干扰;另外,还能提高解密(被破解)的难度。
步骤104,对St、Wt和Rt进行加权叠加,将水印信息嵌入音频图像中,得到嵌入水印信息的音频图像Hw,用公式表示为:
Hw=w1St+w2Wt+w3Rt
式中,w1、w2、w3分别为St、Wt和Rt的加权系数。
在本步骤中,为了不影响音频信号的质量,水印图像和随机模板的加权系数w2和w3应较低,音频图像的加权系数w1应较高,具体数值通过实验确定。
步骤105,定义DFD变换:当物平面x0o0y0的物光分布U0(x0,y0),物平面x0o0y0和观测平面xoy之间(没有透镜)的距离为d时观测平面xoy的物光波场分布DFD[U0(x0,y0)],称DFD[U0(x0,y0)]为DFD变换;通过计算在d=z′1处得到音频图像的复振幅分布,z′1满足f为透镜的焦距;
在本步骤中,首先根据定义推导一般情况下观测平面xoy的物光波场分布DFD[U0(x0,y0)],然后通过计算在d=z′1处得到音频图像的复振幅分布。由于且z′1满足因此,只有音频图像能在d=z′1处成清晰的像,而水印图像和随机模板都只能成模糊的像,且趋于均匀分布,这相当于在近似均匀的背景上叠加了清晰的音频图像,所以在d=z′1处获得的光场分布Sw中,音频图像起决定作用,从Sw中可以有效恢复出一维的音频信号。因此,通过计算在d=z′1处得到了音频图像的复振幅分布。
步骤106,将Sw、R以及用作解密密钥的参数传送给用户或第三方,所述参数包括光的波长λ、透镜的焦距f以及z1、z2和z3。
在本步骤中,可以通过网络、无线电、拷贝等方式将Sw、R以及用作解密密钥的参数传送给用户或第三方。R及解密密钥用于从Sw中恢复音频信号和水印信息。
在本实施例中,通过对St、Wt和Rt进行加权叠加,实现了水印信息的隐藏。由于光学衍射成像具有的多维、大容量等优点,因此,可以将大容量水印信息隐藏在音频信息中,解决了现有方法水印信息嵌入量很小,信息隐藏受到很大限制的问题。
作为一种可选实施例,所述将截取的音频信号转换成二维的音频图像S,包括:
步骤1011,将截取的音频信号归一化:
F(i)=f(i)/am
am=max[abs(f(i))]
式中,f(i)为截取的音频信号,i=1,2,…N,N为数据点数,F(i)为归一化后的音频信号。
步骤1012,将F(i)转换成不小于0的信号G(i):
G(i)=[F(i)+1]/2
步骤1013,将G(i)转化成二维矩阵,并存储为图像格式。
在本实施例中,可利用MATLAB内置的单下标与全下标转变函数ind2sub将G(i)转化成二维矩阵;利用MATLAB内置函数Imwrite将二维矩阵保存为图像格式,即JPEG格式。
作为一种可选实施例,所述分别计算将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt,包括:
计算透镜前d0处的物平面x0o0y0的物光分布为U0(x0,y0)时,透镜后d处的观测平面xoy的物光波场分布U(x,y):
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程;
将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt分别为:
式中,S0、W0、R0分别为音频图像S、水印图像W和随机模板R的像素强度分布函数。
本发明还提供一种实施例从嵌入水印信息的音频图像中提取水印信息的方法,包括以下步骤:
步骤201,用户或第三方接收音频图像的复振幅分布Sw、随机模板R以及用作解密密钥的参数,所述参数包括光的波长λ、透镜的焦距f以及音频图像S、水印图像W和随机模板R置于透镜前的距离z1、z2和z3;
在本步骤中,用户或第三方可以通过网络、无线电、拷贝等方式接收音频图像的复振幅分布Sw、随机模板R以及用作解密密钥的参数。
步骤202,对Sw取模平方得到音频图像的强度分布|Sw|2;
在本步骤中,对复振幅分布Sw取模平方得到音频图像的强度分布|Sw|2,|Sw|2为二维数据。
步骤203,将二维的音频图像的强度分布|Sw|2转化为一维的音频信号;
在本步骤中,音频图像的强度分布|Sw|2为二维数据,将二维数据转换为一维数据后得到一维的音频信号。
本实施例得到的音频信号可以在音频播放装置中直接进行播放。由于水印图像和随机模板在Sw中可以看作是音频图像的近似均匀分布的背景场,而均匀分布的背景场在恢复出的音频信号中近似为直流成分,音频播放装置会自动滤除这些直流成分(比如通过隔直电容阻断直流成分),因此不会影响音频信号的正常播放。
步骤204,计算当收到的随机模板R位于透镜前z3时,在透镜后z′1处的物光场分布R′,并计算S′w=SW-R′;
在本步骤中,先计算随机模板R位于透镜前z3时,在d=z′1处所成像的物光场分布R′,然后通过计算S′w=SW-R′消除随机模板R在Sw中形成的近似均匀分布的背景成份。经这样处理后,S′w只包含音频图像和水印图像成份,消除随机模板成份后,便于提高水印信息提取的精度。
步骤205,定义DFD变换:当物平面x0o0y0的物光分布U0(x0,y0),物平面x0o0y0和观测平面xoy之间的距离为d时观测平面的物光波场分布DFD[U0(x0,y0)],称DFD[U0(x0,y0)]为DFD变换;
在本步骤中,定义DFD变换是为了后面步骤计算的方便,因为后面的步骤中多次用到DFD变换、逆变换,定义DFD变换后可以直接引用。DFD变换属于现有技术,只要给出DFD变换的含义就能得到DFD变换的表达式,而且采用不同的方法可以得到形式完全不同的表达式,但计算得到的结果数据相同。所以这里不给出具体的表达式。
步骤206,对S′w进行DFD逆变换得:
在本步骤中,对S′w进行逆DFD变换可通过求得到。对S′w进行逆DFD变换得到透镜后表面的物光场分布Q。
步骤207,对Q进行DFD变换得到水印信息W′:
在本步骤中,由于z′2满足因此,只有水印图像在d=z′2处能成清晰的像,大大抑制了Q中的音频图像成份,从而得到水印信息W′。
在本实施例中,通过提取水印信息,可以实现音频信息来源可靠性的认证与证明,如水印中的口授人照片、公证处的完整文字、签章等信息可以充分显示与证明信息的来源及其真伪;由于在提取水印时不需要原始载体图像,只需获得光的波长λ、透镜的焦距f以及音频图像、水印图像和随机模板置于透镜前的距离z1、z2和z3,就能提取水印信息,实现了水印信息的盲提取。
作为一种可选实施例,所述将二维的音频图像的强度分布|Sw|2转化为一维的音频信号,包括:
步骤2031,将二维的|Sw|2转化为一维数据f′(i),i=1,2,…N,N为数据点数;
在本步骤中,可利用全下标与单下标转变函数,将二维的|Sw|2转化为一维数据f′(i),如MATLAB内置的函数ind2sub。
步骤2032,对f′(i)进行归一化:
F′(i)=f′(i)/a′m
a′m=max[f′(i)]
步骤2033,恢复负值:
G′(i)=2F′(i)-1.0
步骤2034,乘一增益系数km得到音频信号:
f″(i)=kmG′(i)。
作为一种可选实施例,所述计算当收到的随机模板R位于透镜前z3时,在透镜后z′1处的物光场分布R′,包括:
计算透镜前d0处的物平面x0o0y0的物光分布为U0(x0,y0)时,透镜后d处的观测平面xoy的物光波场分布U(x,y):
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程;
当收到的随机模板R位于透镜前z3时,在透镜后z′1处的物光场分布R′为:
式中,R0为随机模板R的像素强度分布函数。
作为一种可选实施例,所述DFD变换的表达式为:
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程。
本实施例给出了DFD变换的一种表达式。如果采用不同的衍射场求解方法,DFD变换的表达式也不同,但最终的求解结果相同。
上述仅对本发明中的几种具体实施例加以说明,但并不能作为本发明的保护范围,凡是依据本发明中的设计精神所做出的等效变化或修饰或等比例放大或缩小等,均应认为落入本发明的保护范围。
Claims (7)
1.一种在音频信号中加入水印的方法,其特征在于,包括以下步骤:
从一维的音频信号中截取一段,将截取的音频信号转换成二维的音频图像S,并保存为图像格式;
生成随机模板R,并保存为图像格式;
分别计算将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt;
对St、Wt和Rt进行加权叠加,将水印信息嵌入音频图像中,得到嵌入水印信息的音频图像Hw,用公式表示为:
Hw=w1St+w2Wt+w3Rt
式中,w1、w2、w3分别为St、Wt和Rt的加权系数;
定义DFD变换:当物平面x0o0y0的物光分布U0(x0,y0),物平面x0o0y0和观测平面xoy之间的距离为d时观测平面的物光波场分布DFD[U0(x0,y0)],称DFD[U0(x0,y0)]为DFD变换;通过计算在d=z′1处得到音频图像的复振幅分布,z′1满足f为透镜的焦距;
将Sw、R以及用作解密密钥的参数传送给用户或第三方,所述参数包括光的波长λ、透镜的焦距f以及z1、z2和z3。
2.根据权利要求1所述的在音频信号中加入水印的方法,其特征在于,所述将截取的音频信号转换成二维的音频图像S,包括:
将截取的音频信号归一化:
F(i)=f(i)/am
am=max[abs(f(i))]
式中,f(i)为截取的音频信号,i=1,2,…N,N为数据点数,F(i)为归一化后的音频信号;
将F(i)转换成不小于0的信号G(i):
G(i)=[F(i)+1]/2
将G(i)转化成二维矩阵,并保存为图像格式。
3.根据权利要求1所述的在音频信号中加入水印的方法,其特征在于,所述分别计算将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt,包括:
计算透镜前d0处的物平面x0o0y0的物光分布为U0(x0,y0)时,透镜后d处的观测平面xoy的物光波场分布U(x,y):
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程;
将音频图像S、水印图像W和随机模板R置于透镜前z1、z2和z3时,在透镜后表面的衍射光场分布St、Wt和Rt分别为:
式中,S0、W0、R0分别为音频图像S、水印图像W和随机模板R的像素强度分布函数。
4.一种从嵌入水印信息的音频图像中提取水印信息的方法,其特征在于,包括以下步骤:
用户或第三方接收音频图像的复振幅分布Sw、随机模板R以及用作解密密钥的参数,所述参数包括光的波长λ、透镜的焦距f以及音频图像S、水印图像W和随机模板R置于透镜前的距离z1、z2和z3;
对Sw取模平方得到音频图像的强度分布|Sw|2;
将二维的音频图像的强度分布|Sw|2转化为一维的音频信号;
计算当收到的随机模板R位于透镜前z3时,在透镜后z′1处的物光场分布R′,并计算S′w=SW-R′;
定义DFD变换:当物平面x0o0y0的物光分布U0(x0,y0),物平面x0o0y0和观测平面xoy之间的距离为d时观测平面的物光波场分布DFD[U0(x0,y0)],称DFD[U0(x0,y0)]为DFD变换;
对S′w进行DFD逆变换得:
对Q进行DFD变换得到水印信息W′:
5.根据权利要求4所述的从嵌入水印信息的音频图像中提取水印信息的方法,其特征在于,所述将二维的音频图像的强度分布|Sw|2转化为一维的音频信号,包括:
将二维的|Sw|2转化为一维数据f′(i),i=1,2,…N,N为数据点数;
对f′(i)进行归一化:
F′(i)=f′(i)/a′m
a′m=max[f′(i)]
恢复负值:
G′(i)=2F′(i)-1.0
乘一增益系数km得到音频信号:
f″(i)=kmG′(i)。
6.根据权利要求4所述的从嵌入水印信息的音频图像中提取水印信息的方法,其特征在于,所述计算当收到的随机模板R位于透镜前z3时在透镜后z′1处的物光场分布R′,包括:
计算透镜前d0处的物平面x0o0y0的物光分布为U0(x0,y0)时,透镜后d处的观测平面xoy的物光波场分布U(x,y):
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程;
当收到的随机模板R位于透镜前z3时在透镜后z′1处的物光场分布R′为:
式中,R0为随机模板R的像素强度分布函数。
7.根据权利要求4所述的从嵌入水印信息的音频图像中提取水印信息的方法,其特征在于,所述DFD变换的表达式为:
式中,j为虚数单位,k=2π/λ,π为圆周率,λ为光的波长、f为透镜的焦距,L为光波沿光轴的光程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710068630.6A CN106971366B (zh) | 2017-02-08 | 2017-02-08 | 一种在音频信号中加入及提取水印的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710068630.6A CN106971366B (zh) | 2017-02-08 | 2017-02-08 | 一种在音频信号中加入及提取水印的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106971366A true CN106971366A (zh) | 2017-07-21 |
CN106971366B CN106971366B (zh) | 2020-09-11 |
Family
ID=59334774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710068630.6A Active CN106971366B (zh) | 2017-02-08 | 2017-02-08 | 一种在音频信号中加入及提取水印的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106971366B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108508772A (zh) * | 2018-04-26 | 2018-09-07 | 王琰 | 一种信号安全输出技术和识别方法 |
CN110135134A (zh) * | 2019-04-11 | 2019-08-16 | 北京印刷学院 | 一种基于视频的数字全息信息隐藏和处理方法 |
CN113515753A (zh) * | 2020-04-11 | 2021-10-19 | 阿里巴巴集团控股有限公司 | 一种信息隐藏方法及装置 |
WO2022179264A1 (zh) * | 2021-02-27 | 2022-09-01 | 腾讯音乐娱乐科技(深圳)有限公司 | 一种音频生成方法及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1487421A (zh) * | 2003-08-11 | 2004-04-07 | 深圳大学 | 基于虚拟光学的三维数字水印嵌入/检测方法及装置 |
CN1269046C (zh) * | 2003-08-04 | 2006-08-09 | 深圳大学 | 基于虚拟成像的信息隐藏加/解密方法 |
CN102136856A (zh) * | 2010-01-21 | 2011-07-27 | 中兴通讯股份有限公司 | 一种确定立方量度和最大回退功率的方法和装置 |
US20150242981A1 (en) * | 2010-03-05 | 2015-08-27 | Digimarc Corporation | Extending Watermark Detection Distortion Tolerances |
-
2017
- 2017-02-08 CN CN201710068630.6A patent/CN106971366B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1269046C (zh) * | 2003-08-04 | 2006-08-09 | 深圳大学 | 基于虚拟成像的信息隐藏加/解密方法 |
CN1487421A (zh) * | 2003-08-11 | 2004-04-07 | 深圳大学 | 基于虚拟光学的三维数字水印嵌入/检测方法及装置 |
CN102136856A (zh) * | 2010-01-21 | 2011-07-27 | 中兴通讯股份有限公司 | 一种确定立方量度和最大回退功率的方法和装置 |
US20150242981A1 (en) * | 2010-03-05 | 2015-08-27 | Digimarc Corporation | Extending Watermark Detection Distortion Tolerances |
Non-Patent Citations (3)
Title |
---|
尚春杰 等: "《低信噪比下的信号检测与参数估计算法》", 《信息传输与接入技术》 * |
张鹏: "《基于信息光学的多维数据加密及数字水印》", 《中国优秀博硕士学位论文全文数据库 (博士) 信息科技辑 (月刊) 计算机软件及计算机应用》 * |
钟金刚,陈家楠,赵雅静: "全数字全息术在音频信息加密中的应用", 《信息安全与通信保密》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108508772A (zh) * | 2018-04-26 | 2018-09-07 | 王琰 | 一种信号安全输出技术和识别方法 |
CN110135134A (zh) * | 2019-04-11 | 2019-08-16 | 北京印刷学院 | 一种基于视频的数字全息信息隐藏和处理方法 |
CN113515753A (zh) * | 2020-04-11 | 2021-10-19 | 阿里巴巴集团控股有限公司 | 一种信息隐藏方法及装置 |
CN113515753B (zh) * | 2020-04-11 | 2022-05-31 | 阿里巴巴集团控股有限公司 | 一种信息隐藏方法及装置 |
WO2022179264A1 (zh) * | 2021-02-27 | 2022-09-01 | 腾讯音乐娱乐科技(深圳)有限公司 | 一种音频生成方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN106971366B (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hsu et al. | Copyright protection scheme for digital images using visual cryptography and sampling methods | |
Arnold et al. | Techniques and applications of digital watermarking and content protection | |
CN106971366A (zh) | 一种在音频信号中加入及提取水印的方法 | |
Chen et al. | A novel digital watermarking based on general non-negative matrix factorization | |
CN106251278B (zh) | 一种基于直方图特征的图像加密域可逆信息隐藏方法 | |
Cox et al. | Digital watermarking and steganography | |
CN109189727A (zh) | 一种基于属性代理重加密的区块链密文云存储共享方法 | |
US8959357B2 (en) | Biometric encryption and key generation | |
CN107103577A (zh) | 用于完整性保护的半脆弱唯一性图像水印嵌入及提取方法 | |
Abood et al. | Audio steganography with enhanced LSB method for securing encrypted text with bit cycling | |
CN109151157A (zh) | 一种多媒体数字水印取证手机 | |
CN107426172A (zh) | 一种身份信息的匹配推送方法及装置 | |
CN109344578A (zh) | 基于混沌和小波变换的音频水印嵌入、提取方法 | |
CN101562681B (zh) | 基于状态参数的含水印数字照片的生成方法 | |
CN108614961A (zh) | 一种数据库数字水印嵌入和保护方法 | |
CN106653053B (zh) | 一种基于全息图像的音频加密解密方法 | |
CN1889428B (zh) | 网络环境下生物特征数据的安全隐匿传输方法 | |
Kshirsagar et al. | Anatomized study of security solutions for multimedia: deep learning-enabled authentication, cryptography and information hiding | |
Saraswat et al. | A review of digital image steganography | |
CN105337742A (zh) | 基于人脸图像特征和gps信息的lfsr文件加密及解密方法 | |
CN104463772A (zh) | 一种用于版权保护的零水印方法 | |
Umar et al. | I-Marks: An iris code embedding system for ownership identification of multimedia content | |
CN110580405B (zh) | 一种基于社交网络的图片版权保护系统及方法 | |
Mandal et al. | A mathematical model for secret message passing using Steganography | |
Mendua | A new approach of colour image encryption based on Henon like chaotic map |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |