CN105304091B - 一种基于dct的语音篡改恢复方法 - Google Patents
一种基于dct的语音篡改恢复方法 Download PDFInfo
- Publication number
- CN105304091B CN105304091B CN201510608654.7A CN201510608654A CN105304091B CN 105304091 B CN105304091 B CN 105304091B CN 201510608654 A CN201510608654 A CN 201510608654A CN 105304091 B CN105304091 B CN 105304091B
- Authority
- CN
- China
- Prior art keywords
- frame
- signal
- denoted
- content
- compressed signal
- 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
Links
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明公开了一种基于DCT的语音篡改恢复方法,首先将原始语音信号A分为P帧,每帧分为4段。采用基于DCT的语音信号压缩方法,得到压缩信号;同时对压缩的信号置乱。然后,采用量化信号样本值的方法,将各帧的帧号嵌入到该帧前两段中,并将该帧对应的置乱后的压缩信号嵌入到后两段中。本发明能够精确地定位被攻击的内容,同时对被攻击的信号能够篡改恢复。
Description
技术领域
本发明在给出了一种对语音内容真实性和完整性认证方法的同时,给出了对被攻击信号进行篡改恢复的方法,保障了数字语音内容的安全性,对被攻击内容的恢复使语音的内容更真实可靠。
背景技术
数字语音信号作为一种信息载体广泛应用在电话通讯、新闻报导等方面,已经成为人们日常生活不可分割的一部分。由于自身特点等原因,语音信号容易引起攻击者的兴趣而被攻击,被攻击的信号表达的意义和原始信号会有很大区别,甚至完全相反。这些给人们的生活带来了不便,增加了不和谐因素。于是,需要一种技术来对语音信号的真实性和完整性进行认证。同时,对被攻击信号的篡改恢复,将给人们的生活带来诸多便利。因此,可恢复的语音内容认证算法更具有研究意义和实用价值。
文献“Authenticity examination of compressed audio recordings usingdetection of multiple compression and encoders’identification”(Korycki R,Forensic Science International.Vol.238,pp.33-46,2014)针对压缩的音频信号提出了一种内容认证算法,基于MDCT系数的统计特征以及其它一些从压缩信号中提取的参数来对音频内容进行取证。该方法提高了水印系统的鲁邦性。然而,由于认证过程需要大量的样本进行训练,限制了其在实际生活中的应用。文献“Content-dependent watermarkingscheme in compressed speech with identifying manner and location of attacks”(Chen O T-C,Chia-Yi,Chia-Hsiung Liu,IEEE Transactions on Audio, Speech,andLanguage Processing.Vol.15,no.5,pp.1605-1616,2007)基于压缩编码技术提出了一种语音内容取证方法。一方面,该方法的水印信息在语音压缩过程中生成,对于采用不同方法压缩的语音信号,或者未被压缩的信号,该方法没有取证能力。另一方面,水印嵌入基于LSB的方法。由于LSB的水印嵌入方法是脆弱的,信号处理操作会被认为是恶意攻击,于是该方法在应用中存在局限性。以上方法对于检测到的恶意攻击,都不具有篡改恢复的能力。为此,研究具有篡改恢复能力的数字语音取证技术,不仅具有重要研究意义和实用价值,对提高数字语音内容的可靠性也具有积极作用。
发明内容
本发明的目的在于提供一种基于DCT的语音内容篡改恢复算法,该算法首先能够对恶意攻击进行有效监测并定位,然后对监测到的被攻击内容进行篡改恢复,以此来实现数字语音内容的真实性和完整性认证,并提高数字语音认证的可靠性。
为实现这样的目的,本发明给出了基于DCT的数字语音压缩方法,利用压缩信号来恢复被攻击的内容,设计了一种具有篡改恢复能力的数字语音取证方法。
一种基于DCT的数字语音篡改恢复方法,能够对恶意攻击进行有效监测并定位,以此来实现对数字语音内容的真实性和完整性取证,并通过对被攻击内容的篡改恢复使取证的结果更加真实可信,包括如下具体步骤:
(1)信号压缩:首先将原始语音信号A分为不重叠的P帧,第i帧记为Ai;同时对A进行重采样操作,重采样采用的采样频率F′小于原始信号采样频率F,采样后的信号记为A′;接着将A′分为不重叠的P帧,并对每帧进行DCT;取前 M个低频系数作为该帧的压缩信号;将各帧的压缩信号置乱,置乱后第i帧Ai对应的压缩信号记为Ci;
(2)嵌入帧号和压缩信号:将Ai分为4段,记为A1i、A2i、A3i和A4i;Ai对应的帧号i作为第i帧的标识嵌入在A1i和A2i中;同时,第i帧Ai对应的压缩信号Ci嵌入到A3i和A4i中;嵌入后的信号记为WA;
(3)内容认证:首先将待检测语音信号WA分为不重叠的P帧,第i帧记为 WAi;并将每帧分为4段,分别记为WA1i、WA2i、WA3i和WA4i;从WA1i和WA2i中提取帧号,同时从WA3i和WA4i中提取压缩信号。若从WA1i和WA2i中提取的帧号相同,同时,从WA3i和WA4i中提取压缩信号也相同,则表明该帧的内容是真实的,并跳入第i+1帧WAi+1进行内容认证操作;否则,则说明该帧被恶意攻击;
(4)篡改恢复:在检测到被攻击的内容之后,搜索找到下一个能够通过认证的帧,并提取该帧的帧号和上一个能通过认证的帧的帧号,两个帧号之差就是被攻击的信号;根据压缩信号置乱的方法,找到被攻击帧的压缩信号嵌入的位置;提取压缩信号并恢复被攻击的内容,来实现篡改恢复。
与现有的语音内容取证算法相比,本发明首先通过各帧的帧号来监测被攻击的信号,提高了篡改定位的精度;通过提取与被攻击信号对应的压缩信号,来恢复被攻击的内容,提高了算法的篡改恢复能力,也确保了被攻击的信号的可读性和认可度。既能够篡改监测有能够篡改恢复,有利于本发明在日常生活中的推广应用。
附图说明
图1为语音信号压缩过程框图。
图2为帧号和水印嵌入过程框图。
图3为语音内容取证过程框图。
图4为篡改定位和篡改恢复方法。
图5为本发明选取的含水印语音信号。
图6为删除攻击后的含水印语音信号。
图7为替换攻击后的含水印语音信号。
图8为删除攻击后的篡改定位结果。
图9为替换攻击后的篡改定位结果。
图10为删除攻击后的篡改恢复结果。
图11为替换攻击后的篡改恢复结果。
图12为样本值的分块方法图。
图13为样本值分块方法举例图。
图14为不同类型含水印语音信号的ODG值和SDG值图。
图15为常规信号处理后水印提取的BER值图。
具体实施方式
以下结合附图和实施例对本发明的技术方案作进一步描述。
1、信号压缩:
(1)将原始语音信号A={al,1≤l≤L}分为P帧,第i帧记为Ai, Ai={bi,t,1≤t≤L/P},其中bi,t表示为Ai的第t个样本。
(2)对A进行重采样操作,重采样的采样频率F′小于原始信号采样频率 F,采样后的信号记为A′,其长度记为L′,L′=L·F′/F。将A′分为不重叠的P 帧,第i帧记为A′i,
(3)对A′i进行DCT,所得系数记为Di={di,j,1≤j≤L′/P};取前M个系数并记为Gi,Gi={gi,j|gi,j=di,j1≤j≤M},M<<L′/P。Gi既是第i帧Ai压缩后的信号。信号压缩过程如图1所示。
2、嵌入帧号和压缩信号:
(1)将Ai分为4段,分别记为A1i、A2i、A3i和A4i;A1i和A2i的长为N,A3i和A4i的长为6M。将A3i和A4i等分为M个子段,每个子段含有6个样本点,第 j个子段记为和1≤j≤M。
(2)第i帧帧号i表示为Yi={y1,y2,…,yN},Yi作为第i帧的标识。Yi中的各元素可以由下式得到:
i=y1·10N-1+y210N-2+…+yN
记A1i的前N个样本点为a11,a12,…,a1N。用y1,y2,…,yN依次代替这N个样本点的倒数第二位的整数值,以此来完成Yi的嵌入。用同样的方法将Yi嵌入到 A2i中。
(3)对压缩信号Gi(1≤i≤P)采用混沌地址索引的方法进行置乱。混沌序列由下式Logistic混沌映射生成,其中x0表示混沌序列的初值。
xl+1=μxl(1-xl),3.5699≤μ≤4
记X={xl|l=1,2,…,P},将xl(1≤l≤P)按照下式从大到小排列,其中c(l)表示升序排列后混沌序列的地址索引。
xc(l)=ascend(xl)
置乱后的信号记为Ci={ci,j|1≤j≤M},1≤i≤P。以Ci的第一个系数ci,1为例来介绍嵌入方法,将ci,1嵌入到A3i的第一个子段中(6个样本点),过程如下:
①记将中的6个样本值分为6块,分别记为B1,B2, …,B6。B1由组成,B2由 组成,其中表示向下取整,其它块的组成详见图12,该分块方法的实例由图13给出。
②将ci,1的符号(“+”或“-”)嵌入到B1中,记B1中的三个数之和为T。如果0≤ci,1,并且Tmod2=1,量化B1中的某个数,使Tmod2=0;如果ci,1<0,并且 Tmod2=0,量化B1中的某个数,使Tmod2=1;对于其它情况,不对B1中的值做任何修改。
③记将嵌入到第二个块B2中。记B2中的三个数分别为z1,z2和z3,即由下式计算U=f(z1,z2,z3)。
f(z1,z2,z3)=(z1×1+z2×2+z3×3)mod10
如果z1,z2和z3保持不变;如果量化z1,z2或z3。量化方法为(在z1,z2和z3改变幅度最小的前提下)z1±1,z2±1或z3±1。值得注意的是,如果z1,z2或z3中有为0的情况,在计算U之前先将它们赋值为1。
采用相同的方法,将ci,1其它位的值( )嵌入到块B3,B4,B5,B6中。
④依据第①步到第③步方法,将Ci嵌入到A3i和A4i中。帧号和水印嵌入过程框图如图2所示。
3、内容认证:
(1)将待检测语音信号WA分为不重叠的P帧,第i帧记为WAi;并将每帧分为4段,分别记为WA1i、WA2i、WA3i和WA4i,其中WA3i和WA4i的长为6M。
(2)从WA1i和WA2i中提取帧号。取WA1i的前N个样本点,并记这些样本点倒数第二位的整数值分别为y′N,y′N-1,……y′1;同样的,取WA2i的前N个样本点,这些样本点倒数第二位的整数值记为
(3)从WA3i和WA4i中提取压缩信号。将WA3i和WA4i等分为M个子段,每个子段含有6个样本点,第j个子段记为和1≤j≤M。以从中提取系数c′i,1为例来介绍提取方法,记提取过程为:
①采用图12的方法将中的6个样本点分为6块,分别记为B1,B2, …,B6。
②提取c′i,1的符号。记B1中三个值的和为T,T=w1+w2+w3。如果Tmod2=0,则c′i,1的符号为“+”;如果Tmod2=1,则c′i,1的符号为“-”。
③c′i,1的最高位记为c′i,1的其它位整数,和可采用同样的方法从B3,B4,B5和B6中提取得到。
④根据提取的符号和各位的整数,由下式构造系数c′i,1,
依照上面的方法,分别从WA3i和WA4i中提取压缩信号,分别记为c′i,j和 1≤j≤M。
⑤内容认证。如果且1≤l≤N,1≤j≤M,则表明第i帧的内容是真实;否则,表明第i帧的内容被恶意攻击。内容认证过程框图如图3 所示。
4、篡改恢复:
在检测到被攻击的内容之后,移动样本,搜索找到下一个能够通过认证的语音帧。提取该帧的帧号和上一个通过认证的帧号,两个帧号之差就是被攻击的信号;根据压缩信号的置乱方法,找到被攻击帧对应压缩信号的嵌入位置;提取压缩信号并恢复被攻击的内容,以此实现篡改恢复。篡改定位和篡改恢复的方法如图4所示。
本发明方法的效果可以通过以下的性能分析验证:
1、不可听性
选取采样率为44.1kHz,样本长度为102400,16比特量化的单声道WAVE 格式语音信号作为测试样本。分为4类,分别记为T1,T2,T3和T4。T1为在安静的会议中录制的信号,T2为在讨论会中录制的信号,T3为在吵闹的车站录制的信号,T4为在野外录制的信号。图14给出了4类不同信号的ODG 值和SDG值。ODG值由PEAQ系统测试所得,SDG值由10位听众现场打分所得。由图14所给结果可以看出,本方法具有较好的不可听性。
2、鲁棒性
用误码率BER(bit error rate)来测试本方法对信号处理的鲁棒性,BER的定义为
其中,E为提取水印错误比特数,T为语音信号所嵌水印总比特数。BER值越小说明算法对信号处理的鲁棒性越强。
对4类在不同环境下录制的语音信号采用本发明方法嵌入水印,然后对含水印信号进行不同类型的信号处理,分别是低通滤波(截止频率16kHz),重采样(先下采样至22.05kHz,再上采样至44.1kHz),MP3压缩(压缩率为 64kbps)。信号处理后的BER值见图15。由图15数据可得,本发明对常见信号处理具有鲁棒性。
3、篡改恢复
对如图5所示的含水印语音信号进行删除攻击和替换攻击。攻击后的语音信号如图6和图7所示,篡改检测结果分别如图8和图9所示。对应的篡改恢复结果如图10和图11所示。为了更清晰地显示被攻击的帧,篡改检测结果中仅显示了部分能够正确提取的帧号,TL(i)=1表明第i帧的语音内容是真实的。由篡改恢复结果可以看出,本算法对恶意攻击具有篡改恢复能力。
Claims (1)
1.一种基于DCT的语音篡改恢复方法,对恶意攻击检测并精确定位被攻击的内容后对被攻击的内容进行篡改恢复,首先将原始语音信号A分为P帧,每帧分为4段,采用基于DCT的语音信号压缩方法,得到压缩信号;同时对压缩的信号置乱;然后,采用量化信号样本值的方法,将各帧的帧号嵌入到该帧前两段中,并将该帧对应的置乱后的压缩信号嵌入到后两段中,包括如下具体步骤:
(1)信号压缩:首先将原始语音信号A分为不重叠的P帧,第i帧记为Ai;同时对A进行重采样操作,重采样采用的采样频率F′小于原始信号采样频率F,采样后的信号记为A′;接着将A′分为不重叠的P帧,并对每帧进行DCT;取前M个低频系数作为该帧的压缩信号;将各帧的压缩信号置乱,置乱后第i帧Ai对应的压缩信号记为Ci;
(2)嵌入帧号和压缩信号:将Ai分为4段,记为A1i、A2i、A3i和A4i;Ai对应的帧号i作为第i帧的标识嵌入在A1i和A2i中;同时,第i帧Ai对应的压缩信号Ci嵌入到A3i和A4i中;嵌入后的信号记为WA;
(3)内容认证:首先将待检测语音信号WA分为不重叠的P帧,第i帧记为WAi;并将每帧分为4段,分别记为WA1i、WA2i、WA3i和WA4i;从WA1i和WA2i中提取帧号,同时从WA3i和WA4i中提取压缩信号;若从WA1i和WA2i中提取的帧号相同,同时,从WA3i和WA4i中提取压缩信号也相同,则表明该帧的内容是真实的,并跳入第i+1帧WAi+1进行内容认证操作;否则,则说明该帧被恶意攻击;
(4)篡改恢复:在检测到被攻击的内容之后,搜索找到下一个能够通过认证的帧,并提取该帧的帧号和上一个能通过认证的帧的帧号,两个帧号之差就是被攻击的信号;根据压缩信号置乱的方法,找到被攻击帧的压缩信号嵌入的位置;提取压缩信号并恢复被攻击的内容,来实现篡改恢复。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510608654.7A CN105304091B (zh) | 2015-06-26 | 2015-09-23 | 一种基于dct的语音篡改恢复方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2015103645597 | 2015-06-26 | ||
CN201510364559 | 2015-06-26 | ||
CN201510608654.7A CN105304091B (zh) | 2015-06-26 | 2015-09-23 | 一种基于dct的语音篡改恢复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105304091A CN105304091A (zh) | 2016-02-03 |
CN105304091B true CN105304091B (zh) | 2018-10-26 |
Family
ID=55201262
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510608654.7A Expired - Fee Related CN105304091B (zh) | 2015-06-26 | 2015-09-23 | 一种基于dct的语音篡改恢复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105304091B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105895109B (zh) * | 2016-05-10 | 2019-02-26 | 信阳师范学院 | 一种基于dwt和dct的数字语音取证和篡改恢复方法 |
CN106531176B (zh) * | 2016-10-27 | 2019-09-24 | 天津大学 | 音频信号篡改检测与恢复的数字水印算法 |
CN109754817A (zh) * | 2017-11-02 | 2019-05-14 | 北京三星通信技术研究有限公司 | 信号处理方法及终端设备 |
CN108877819B (zh) * | 2018-07-06 | 2023-05-26 | 信阳师范学院 | 一种基于系数自相关度的语音内容取证方法 |
CN109102799B (zh) * | 2018-08-17 | 2023-01-24 | 信阳师范学院 | 一种基于频域系数对数和的语音端点检测方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2301218A1 (en) * | 1997-08-26 | 1999-03-04 | Digimarc Corporation | Method and apparatus for watermarking video images |
US6442283B1 (en) * | 1999-01-11 | 2002-08-27 | Digimarc Corporation | Multimedia data embedding |
CN101247520A (zh) * | 2008-02-27 | 2008-08-20 | 北京航空航天大学 | 一种视频数据加密/解密方法 |
CN101290772A (zh) * | 2008-03-27 | 2008-10-22 | 上海交通大学 | 基于混合域系数矢量量化的音频零水印嵌入和提取方法 |
CN102097098A (zh) * | 2011-02-25 | 2011-06-15 | 宁波大学 | 一种以压缩音频为掩蔽载体的数字隐写及提取方法 |
CN102157154A (zh) * | 2011-01-28 | 2011-08-17 | 桂林电子科技大学 | 基于音频内容的非均匀离散余弦变换音频可靠性认证方法 |
CN102867513A (zh) * | 2012-08-07 | 2013-01-09 | 西南交通大学 | 一种基于伪Zernike矩的语音内容认证方法 |
CN103208288A (zh) * | 2013-03-13 | 2013-07-17 | 漳州职业技术学院 | 基于双重加密的dwt—dct域音频公开水印算法 |
CN103456308A (zh) * | 2013-08-05 | 2013-12-18 | 西南交通大学 | 一种可恢复的密文域语音内容认证方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20070119032A (ko) * | 2005-04-13 | 2007-12-18 | 코닌클리케 필립스 일렉트로닉스 엔.브이. | 위상 변조 이전에 워터마킹을 통한 인코딩 |
-
2015
- 2015-09-23 CN CN201510608654.7A patent/CN105304091B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2301218A1 (en) * | 1997-08-26 | 1999-03-04 | Digimarc Corporation | Method and apparatus for watermarking video images |
US6442283B1 (en) * | 1999-01-11 | 2002-08-27 | Digimarc Corporation | Multimedia data embedding |
CN101247520A (zh) * | 2008-02-27 | 2008-08-20 | 北京航空航天大学 | 一种视频数据加密/解密方法 |
CN101290772A (zh) * | 2008-03-27 | 2008-10-22 | 上海交通大学 | 基于混合域系数矢量量化的音频零水印嵌入和提取方法 |
CN102157154A (zh) * | 2011-01-28 | 2011-08-17 | 桂林电子科技大学 | 基于音频内容的非均匀离散余弦变换音频可靠性认证方法 |
CN102097098A (zh) * | 2011-02-25 | 2011-06-15 | 宁波大学 | 一种以压缩音频为掩蔽载体的数字隐写及提取方法 |
CN102867513A (zh) * | 2012-08-07 | 2013-01-09 | 西南交通大学 | 一种基于伪Zernike矩的语音内容认证方法 |
CN103208288A (zh) * | 2013-03-13 | 2013-07-17 | 漳州职业技术学院 | 基于双重加密的dwt—dct域音频公开水印算法 |
CN103456308A (zh) * | 2013-08-05 | 2013-12-18 | 西南交通大学 | 一种可恢复的密文域语音内容认证方法 |
Non-Patent Citations (5)
Title |
---|
A novel speech content authentication algorithm based on Bessel–Fourier moments;Zhenghui Liu et al.;《Digital Signal Processing》;20130919;第197–208页 * |
An Efficient Speech Content Authentication Algorithm Based on Coefficients Self-correlation Degree;Zhenghui Liu et al.;《Springer-Verlag Berlin Heidelberg》;20131231;第311–326页 * |
一种基于DCT变换域的数字图像水印技术的实现方法;冯岩 等;《信阳师范学院学报:自然科学版》;20081031;第21卷(第4期);第591-593页 * |
一种新的用于音频认证数字水印算法;申军伟 等;《计算机应用研究》;20121031;第29卷(第10期);第3850-3853页 * |
基于基音周期的语音内容认证算法;钱清 等;《铁道学报》;20140930;第36卷(第9期);第60-67页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105304091A (zh) | 2016-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105304091B (zh) | 一种基于dct的语音篡改恢复方法 | |
JP6535706B2 (ja) | データセットの3値ビットマップを作成するための方法 | |
CN101345054B (zh) | 用于声频文件的数字水印制作及识别方法 | |
CN101290772B (zh) | 基于混合域系数矢量量化的音频零水印嵌入和提取方法 | |
CN105550257A (zh) | 一种音视频指纹识别方法及一种基于音视频指纹流媒体的防篡改系统 | |
TW201447871A (zh) | 聲音訊號之特徵指紋擷取與匹配 | |
CN108198563A (zh) | 一种数字版权保护与内容认证的多功能音频保护方法 | |
CN106683687B (zh) | 异常声音的分类方法和装置 | |
CN114387977B (zh) | 基于双域深度特征和注意力机制的语音裁剪痕迹定位方法 | |
Chen et al. | Audio hash function based on non-negative matrix factorisation of mel-frequency cepstral coefficients | |
CN108682425B (zh) | 一种基于恒定水印的鲁棒数字音频水印嵌入系统 | |
CN107507626A (zh) | 一种基于语音频谱融合特征的手机来源识别方法 | |
CN110647656A (zh) | 一种利用变换域稀疏化和压缩降维的音频检索方法 | |
CN103594083A (zh) | 通过电视伴音自动识别电视节目的技术 | |
CN107993669A (zh) | 基于修改最低有效位数权重的语音内容认证和篡改恢复方法 | |
CN107274912B (zh) | 一种手机录音的设备来源辨识方法 | |
CN105895109B (zh) | 一种基于dwt和dct的数字语音取证和篡改恢复方法 | |
CN102867513B (zh) | 一种基于伪Zernike矩的语音内容认证方法 | |
CN103456308A (zh) | 一种可恢复的密文域语音内容认证方法 | |
CN101350198B (zh) | 基于骨导的语音压缩水印方法 | |
CN101609675A (zh) | 一种基于质心的脆弱音频水印方法 | |
Chen et al. | A robust feature extraction algorithm for audio fingerprinting | |
CN108877816B (zh) | 基于qmdct系数的aac音频重压缩检测方法 | |
CN113762110B (zh) | 一种执法即时证据固定方法及执法仪 | |
CN104134443A (zh) | 对称三进制串表示的语音感知哈希序列构造和认证方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20181026 Termination date: 20190923 |
|
CF01 | Termination of patent right due to non-payment of annual fee |