CN106448686B - 一种自适应音频空域隐写方法 - Google Patents

一种自适应音频空域隐写方法 Download PDF

Info

Publication number
CN106448686B
CN106448686B CN201610651588.6A CN201610651588A CN106448686B CN 106448686 B CN106448686 B CN 106448686B CN 201610651588 A CN201610651588 A CN 201610651588A CN 106448686 B CN106448686 B CN 106448686B
Authority
CN
China
Prior art keywords
audio
sample
wav audio
wav
steganography
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
CN201610651588.6A
Other languages
English (en)
Other versions
CN106448686A (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.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
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 Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN201610651588.6A priority Critical patent/CN106448686B/zh
Publication of CN106448686A publication Critical patent/CN106448686A/zh
Application granted granted Critical
Publication of CN106448686B publication Critical patent/CN106448686B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00Speech 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/0018Speech coding using phonetic or linguistical decoding of the source; Reconstruction using text-to-speech synthesis

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)
  • Signal Processing For Digital Recording And Reproducing (AREA)

Abstract

本发明涉及多媒体信息安全领域,具体涉及一种自适应音频空域隐写方法,包括秘密信息的嵌入和秘密信息的提取两个步骤。本发明方法以WAV音频信号的单个样本作为一个嵌入单元,根据原始WAV音频和经过一次Advanced Audio Coding(AAC)高比特率压缩解压后的WAV音频的信号差值设计每个嵌入单元修改的代价,并利用信息编码工具Syndrome Trellis Codes(STC)自适应嵌入信息并实现正确提取。经大量实验证明,本隐写方法与现有的基于LSB匹配的空域隐写方法相比,可使隐写后的音频具有更高的音频质量及安全性。因此,本发明方法对隐秘通信等应用领域具有重要意义。

Description

一种自适应音频空域隐写方法
技术领域
本发明涉及多媒体信息安全领域,具体涉及一种自适应音频空域隐写方法。
背景技术
音频隐写是信息隐藏技术的一个重要分支。音频信号作为日常生活中不可缺少的传播媒介,其信号中包含的信息冗余使音频成为了信息隐藏的良好载体。目前的音频隐写方法以嵌入信息时所采用的域划分,可分为:空域音频隐藏、压缩域音频隐藏等。现有的空域信息隐藏方法采用的大多是最低有效位修改法(LSB),且其嵌入位置的选择是由随机序列决定的。由于随机序列的不确定性,音频信号中人耳听觉上敏感的样本也有可能会被改动,从而导致音频听觉质量及安全性的降低。
发明内容
为了克服上述现有技术的不足,本发明提出一种自适应音频空域隐写方法,该隐写方法为空域音频隐藏,主要载体为WAV音频。本发明提出的WAV音频空域自适应隐写方法可有效地实现在音频信号的信息冗余位进行信息的嵌入,从而避开音频信号的听觉敏感位置,有效减小隐写时损失的音频质量,安全性也会有所提升。
为了实现上述发明目的,采用的技术方案如下:
一种自适应音频空域隐写方法,包括秘密信息的嵌入和提取两个步骤,所述秘密信息的嵌入包括如下子步骤:(11)高比特率AAC压缩原始WAV音频;(12)原始WAV音频与经过一次高比特率AAC压缩解压的WAV音频样本值相减得残差;(13)选定信息编码工具STC的嵌入载体;(14)根据残差设计STC修改嵌入载体样本值的代价;(15)设置嵌入率并计算需嵌入的秘密信息长度;(16)利用STC进行秘密信息的嵌入。
所述步骤(11)对原始WAV音频进行高比特率AAC压缩;
所述步骤(12)获得原始WAV音频和经过一次高比特率AAC压缩解压后的WAV音频样本之间的残差,具体步骤如下:分别用十进制读取原始WAV音 频和经过一次AAC压缩的WAV音频的所有样本,可得序列X={x1,x2,x3,…,xn}及X’={x′1,x′2,x′3,…,x′n},两种WAV音频的样本总数相同,统一用n表示,其中对xi及x′i分别表示两种WAV音频第i个样本的十进制整数值。序列X与序列X’相减可得残差序列Residual,可表示为{r1,r2,r3,…,rn},n表示WAV音频的样本总数,其中对ri=xi-x’i
所述步骤(13)选取原始WAV音频文件的所有样本值作为信息编码工具STC的嵌入载体,即序列X={x1,x2,x3,…,xn},n为WAV音频样本总数;
所述步骤(14)根据所得的残差序列Residual设计STC修改嵌入载体样本值的代价,设计原则如下:
(Ⅰ)对于原始WAV音频第i个样本值xi,修改的代价包含三种代价:在原值上加1的代价、保持不变的代价、在原值上减1的代价;
(Ⅱ)对于原始WAV音频第i个样本值xi,根据其相应的残差序列Residual中第i个残差值Ri,上述三种代价的设定如下:
情况1:样本值xi在原值上加1的代价为:
情况2:样本值xi保持不变的代价均为0
情况3:样本值xi在原值上减1的代价为:
所述步骤(15)设置秘密信息的嵌入率ratio,其中ratio表示原始WAV音频平均每个样本修改的比特数,单位为bps(bit per sample)。计算需嵌入的秘密信息的长度,记为m,m由WAV音频样本总数n和嵌入率ratio决定,即:m=[n×ratio]。
所述步骤(16)依据设定好的代价及应嵌入的m比特秘密信息对选定的嵌入载体做STC操作,实现最小化所有样本的修改代价之和,从而自适应选取嵌入载体中的若干位进行修改以达到秘密信息嵌入的目的。
所述秘密信息的提取包括如下子步骤:(21)获取隐写WAV音频的提取载体;(22)获取隐写WAV音频所嵌入的秘密信息的长度;(23)利用信息编码工具STC对所述提取载体进行信息提取。
所述步骤(21)读取隐写后的WAV音频,得到序列Y={y1,y2,y3,…,yn},其中n表示隐写WAV音频的样本总数,可知隐写后WAV音频的样本总数与原始WAV音频的样本总数保持一致;将所有的样本值作为STC操作的提取载体;
所述步骤(22)获取嵌入率ratio,其中ratio单位为bps(bit per sample),表示原始WAV音频平均每个样本修改的比特数;计算隐写WAV音频嵌入的秘密信息的长度m,其中m=[n×ratio]。
所述步骤(23)对提取载体进行STC解码操作,其中STC工具可根据提取载体的值与嵌入的秘密信息长度还原嵌入的m比特秘密信息。
与现有技术相比,本发明的有益效果是:本发明基于信息编码工具STC及高质量压缩编码AAC,选取WAV音频的单个样本作为嵌入单元,根据原始WAV音频与经过一次AAC高比特率压缩并解压后的WAV音频样本之间的差值设计STC修改嵌入载体的代价,从而使STC工具可自适应地修改WAV音频样本的若干位,避免了基于随机序列修改的LSB方法的随机性,可有效地实现在音频信号的信息冗余位进行信息的嵌入,从而保护了音频信号听觉敏感位置,有效减小隐写时损失的音频质量,安全性也会有所提升。
本发明的原理是经过一次高质量AAC压缩并解压的WAV音频保留下来的信号为原始WAV音频中与人耳感知系统相关的信息位,原始WAV音频信号与其的差值可判定为感知无关部分也即冗余部分。根据残差设定WAV音频各样本修改的代价可使STC工具自适应地选择嵌入载体中的冗余部分进行修改,从而获得更好的人耳听觉质量。其中信息编码工具STC的特点是,设定的样本修改代价越大,该样本被修改的几率就越低。
附图说明
图1所示为采用本发明进行秘密信息嵌入、提取的算法流程示意图。
图2(a)所示为原始WAV音频的信号波形。
图2(b)所示为经过一次高比特率AAC压缩解压后的WAV音频的信号波形。
图2(c)所示为图2 (a)和图2 (b)信号之间的残差。
图2(d)所示为嵌入率为0.3bps时原始WAV音频各信号段的修改率。由图可知修改的样本点集中在样本值变化较为复杂的信号段,其间包含更多的冗余信息,例如[6000,8000]范围内的样本;而对于接近静音听觉敏感的样本值,例如处于[0,2000]及[12000,14000]范围内的样本,修改的样本比例低于0.7%。
具体实施方式
下面结合附图对本发明做进一步的描述,但本发明的实施方式并不限于此。
图1所示为采用本发明进行秘密信息嵌入、提取的算法流程示意图。
一种自适应音频空域隐写方法,包括秘密信息的嵌入和提取两个步骤,其中秘密信息的嵌入包括如下子步骤:(11)高比特率AAC压缩原始WAV音频;(12)原始WAV音频与经过一次高比特率AAC压缩的WAV音频样本值相减得残差;(13)选定信息编码工具STC的嵌入载体;(14)根据残差设计STC修改嵌入载体样本值的代价;(15)设置嵌入率并计算需嵌入的秘密信息长度;(16)利用STC进行秘密信息的嵌入。
每个步骤具体说明如下:
(11)利用neroAACEnc软件以最高比特率400,000bps对原始WAV音频进行AAC音频压缩,并用neroAACDec软件对压缩后的AAC音频进行解码,得到经过一次高比特率AAC压缩并解压的WAV音频;
(12)分别用十进制读取原始WAV音频和经过一次高比特率AAC压缩并解压的WAV音频的所有样本值,可得序列X={x1,x2,x3,…,xn}及X’={x′1,x′2,x′3,…,x′n},其中n表示WAV音频的样本总数,两种WAV音频的样本总数相同,且对xi及x′i皆为整数。序列X与序列X’的信号波形如图2(a)和图2(b)所示。序列X与序列X’相减可得残差序列Residual,表示为{r1,r2,r3,…,rn},n为WAV音频样本总数,其中ri=xi-x’i。图2(c)为所述残差的波形示意图;
(13)选定原始WAV音频的所有样本值作为信息编码工具STC的嵌入载体,即序列X={x1,x2,x3,…,xn}表示,n为WAV音频样本总数;
(14)根据残差Residual序列设计嵌入载体各样本值被信息编码工具STC 修改的代价。具体设计原则如下:
(Ⅰ)对于原始WAV音频第i个样本值xi,修改的代价包含三种代价:在原值上加1的代价、保持不变的代价、在原值上减1的代价;
(Ⅱ)对于原始WAV音频第i个样本值xi,依据相应的残差序列Residual中第i个残差值ri,上述三种代价的设定如下:
情况1:样本值xi在原值上加1的代价为:
情况2:样本值xi保持不变的代价均为0
情况3:样本值xi在原值上减1的代价为:
(15)设置秘密信息的嵌入率ratio,其中ratio表示原始WAV音频平均每个样本修改的比特数,单位为bps(bit per sample)。计算原始WAV音频需嵌入的秘密信息长度,记为m,其中m由嵌入率ratio和原始WAV音频样本总数n决定,即:m=[n×ratio]。
(16)依据设定好的代价和当前所需嵌入的m比特秘密信息对选定的嵌入载体做STC操作,其中STC工具可选取嵌入载体的若干位进行修改以达到秘密信息嵌入的目的。
本发明秘密信息的提取包括步骤:(21)获取隐写WAV音频的提取载体;(22)获取隐写WAV音频所嵌入的秘密信息的长度;(23)利用信息编码工具STC对所述提取载体进行信息提取。
每个步骤具体说明如下:
(21)读取隐写后的WAV音频,得到序列Y={y1,y2,y3,…,yn},其中n表示隐写后WAV音频的样本总数,可知隐写后WAV音频的样本总数与原始WAV音频的样本总数保持一致;将所有的样本值作为STC操作的提取载体;
(22)获取秘密信息嵌入率ratio,计算隐写WAV音频嵌入的秘密信息长度m,其中m=[n×ratio];
(23)对提取载体进行STC操作,其中STC工具可根据提取载体的值与嵌入的秘密信息长度还原嵌入的m比特秘密信息。下面将结合具体的实验数据对本发明做更进一步的分析说明:
1、音频集合:
本实施例中,利用了10000个单声道长度为4秒的音乐及10000个单声道长度为4秒的语音做测试。
2、隐写工具及方法
最低有效位修改法(LSB)
本发明提出的自适应音频空域隐写方法
3、通用型隐写分析特征
“Q.Z.Liu,A.H.Sung,and M.Y.Qiao.Temporary derivative-based spectrumand mel-cepstrum audio steganalysis.Information IEEE Transactions onInformation Forensics and Security,vol.4,no.3,pp.369-368,Sept.2009.”论文提出了[D-MC]隐写分析特征。
“Q.Z.Liu,A.H.Sung,and M.Y.Qiao.Derivative-based audio stegnalysis.ACMTransactions on Multimedia Computing Communications and Applications,vol.7,no.3,pp.18:1-18:19,Sept.2011.”论文提出了[2D-Markov]隐写分析特征。
“W.Luo,H.Li,Q.Yan,R.Yang,and J.Huang.Improved audiosteganalyticfeature and its applications in audio forensics.Tech.Rep.,Sun Yat-SenUniversity,Guangzhou,P.R.China,510006,2016.”技术报告提出了[Combined time&frequency]隐写分析特征。
4、安全性能分析:
利用LSB最低有效位、本发明提出的自适应音频空域隐写方法生成不同嵌入容量的隐写音频,嵌入率分别设为0.1bps、0.3bps、0.5bps。针对每一种情况,采用三种通用的隐写分析特征,一半样本作为训练集,另一半样本作为测试集,采用集成分类器进行训练测试,下表1列出了各种情况下隐写音频的检测率。可以看出采用本发明描述的隐写方法隐写的音频检测率与基于LSB匹配的空域隐写方法相比均有降低,安全性更好。
表1
5、音频质量分析:
利用LSB最低有效位、本发明提出的自适应音频空域隐写方法生成不同嵌入容量的隐写音频,嵌入率分别为0.1bps、0.3bps、0.5bps。针对每一种情况,计算各隐写WAV音频与原始WAV音频的信噪比(单位:dB)及样本修改率。下表2列出了各种情况下的信噪比(单位:dB)和修改率。可以看出采用本发明描述的隐写方法后,隐写WAV音频的信噪比增加,修改率降低,音频质量有所提升。
表2
以上所述的本发明的实施方式,并不构成对本发明保护范围的限定。任何在本发明的精神原则之内所作出的修改、等同替换和改进等,均应包含在本发明的权利要求保护范围之内。

Claims (3)

1.一种自适应音频空域隐写方法,包括秘密信息的嵌入和秘密信息的提取两个步骤,其特征在于:
所述秘密信息的嵌入包括如下子步骤:
(11)对原始WAV音频数据进行高比特率AAC音频压缩,并对压缩后的AAC音频进行解码,得到经过一次高比特率AAC压缩并解压后的WAV音频;
(12)读取原始WAV音频和经过一次高比特率AAC压缩解压后的WAV音频的所有样本,分别用序列X={x1,x2,x3,...,xn}及X’={x’1,x′2,x′3,...,x′n}表示,其中两种WAV音频样本总数相同,用n表示;序列X与序列X’相减可得残差序列Residual,表示为{r1,r2,r3,...,rn},其中对ri=xi-x′i,n为WAV音频样本总数;
(13)选定原始WAV音频的所有样本值作为信息编码工具STC的嵌入载体,即序列X={x1,x2,x3,...,xn},n为WAV音频样本总数;
(14)根据残差序列Residual设计原始WAV音频各样本值被信息编码工具STC修改的代价;具体设计原则如下:
(I)对于原始WAV音频第i个样本值xi,修改的代价包含三种代价:在原值上加1的代价、保持不变的代价、在原值上减1的代价;
(II)对于原始WAV音频第i个样本值xi,根据其对应的残差序列Residual中第i个残差值ri,上述三种代价的设定如下:
情况1:样本值xi在原值上加1的代价为:
情况2:样本值xi保持不变的代价均为0
情况3:样本值xi在原值上减1的代价为:
(15)设置秘密信息的嵌入率ratio,ratio表示原始WAV音频平均每个样本修改的比特数,单位为bps(bit per sample);计算需嵌入的秘密信息长度,记为m,其中m由设定的嵌入率ratio与原始WAV音频样本总数n决定,即:m=[n×ratio];
(16)依据设定好的代价及需嵌入的m比特秘密信息对选定的嵌入载体做STC操作,实现最小化所有样本的修改代价之和,自适应地选取嵌入载体的若干位进行修改以达到信息嵌入的目的;
所述秘密信息的提取包含如下子步骤:
(21)读取隐写后的WAV音频,得到序列Y={y1,y2,y3,...,yn},其中n表示隐写后WAV音频的样本总数,可知隐写后WAV音频的样本总数与原始WAV音频的样本总数保持一致;将所有的样本值作为STC操作的提取载体;
(22)获取嵌入率ratio,计算隐写后的WAV音频中秘密信息的嵌入长度m,其中m=[n×ratio];
(23)对提取载体进行STC解码操作,其中STC工具可以根据提取载体的所有样本值和嵌入长度m还原嵌入的m比特秘密信息。
2.根据权利要求1所述的方法,其特征在于,步骤(12)采用十进制读取原始WAV音频和经过一次高比特率AAC压缩解压后的WAV音频的所有样本,样本值为整数。
3.根据权利要求2所述的方法,其特征在于,步骤(21)采用十进制读取隐写后的WAV音频的所有样本,样本值为整数。
CN201610651588.6A 2016-08-10 2016-08-10 一种自适应音频空域隐写方法 Expired - Fee Related CN106448686B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610651588.6A CN106448686B (zh) 2016-08-10 2016-08-10 一种自适应音频空域隐写方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610651588.6A CN106448686B (zh) 2016-08-10 2016-08-10 一种自适应音频空域隐写方法

Publications (2)

Publication Number Publication Date
CN106448686A CN106448686A (zh) 2017-02-22
CN106448686B true CN106448686B (zh) 2019-07-19

Family

ID=58184333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610651588.6A Expired - Fee Related CN106448686B (zh) 2016-08-10 2016-08-10 一种自适应音频空域隐写方法

Country Status (1)

Country Link
CN (1) CN106448686B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108271027B (zh) * 2018-01-10 2020-06-12 中国人民解放军战略支援部队信息工程大学 一种图像自适应隐秘信息的提取方法
CN109346091B (zh) * 2018-09-13 2021-12-14 中国科学院信息工程研究所 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统
CN109448744B (zh) * 2018-12-14 2022-02-01 中国科学院信息工程研究所 一种基于符号位自适应嵌入的mp3音频信息隐藏方法和系统
CN110619883B (zh) * 2019-08-23 2023-05-23 平安科技(深圳)有限公司 音乐的信息嵌入方法、提取方法、装置、终端及存储介质
CN111640443B (zh) * 2020-04-14 2023-09-26 宁波大学 一种音频可逆隐写方法和秘密信息提取方法
CN111816195B (zh) * 2020-07-03 2023-09-15 杭州秀秀科技有限公司 一种音频可逆隐写方法、秘密信息提取及载体音频恢复方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104485109A (zh) * 2013-01-22 2015-04-01 宁波大学 一种基于索引值的压缩域音频数字提取方法
CN104852799A (zh) * 2015-05-12 2015-08-19 陕西师范大学 基于分段序列的数字音频伪装及重构方法
EP2928154A2 (en) * 2014-04-03 2015-10-07 Barclays Bank PLC Mobile user authentication applying a call identifier
CN105118512A (zh) * 2015-07-23 2015-12-02 武汉大学 一种面向aac数字音频的通用隐写分析方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7076659B2 (en) * 2002-02-25 2006-07-11 Matsushita Electric Industrial Co., Ltd. Enhanced method for digital data hiding

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104485109A (zh) * 2013-01-22 2015-04-01 宁波大学 一种基于索引值的压缩域音频数字提取方法
EP2928154A2 (en) * 2014-04-03 2015-10-07 Barclays Bank PLC Mobile user authentication applying a call identifier
CN104852799A (zh) * 2015-05-12 2015-08-19 陕西师范大学 基于分段序列的数字音频伪装及重构方法
CN105118512A (zh) * 2015-07-23 2015-12-02 武汉大学 一种面向aac数字音频的通用隐写分析方法

Also Published As

Publication number Publication date
CN106448686A (zh) 2017-02-22

Similar Documents

Publication Publication Date Title
CN106448686B (zh) 一种自适应音频空域隐写方法
CN104937844B (zh) 在不同回放设备之间优化响度和动态范围
KR102155491B1 (ko) 다수의 미디어 처리 노드들을 갖는 적응적 처리
HUE027046T2 (en) Devices, processes and computer program products for adaptive encoding and decoding of a watermarked signal
CN101290773B (zh) 自适应的mp3数字水印嵌入和提取方法
CN101211562A (zh) 一种数字音乐作品的无损数字水印嵌入与提取方法
Li et al. Audio perceptual hashing based on nmf and mdct coefficients
Lu et al. Watermarking combined with CELP speech coding for authentication
Wu et al. Robust and blind audio watermarking algorithm in dual domain for overcoming synchronization attacks
Mandal et al. An approach for enhancing message security in audio steganography
CN1971734A (zh) 一种数字音乐作品认证信息的嵌入与提取方法
CN107071455A (zh) 基于数据流的jpeg图像信息隐藏方法
US20190189135A1 (en) Method and System for Data-Hiding Within Audio Transmissions
Sun et al. Data hiding in audio based on audio-to-image wavelet transform and vector quantization
Qian et al. A fragile watermarking algorithm for speech authentication by modifying least significant digits
Wei et al. Controlling bitrate steganography on AAC audio
Bazyar et al. A robust data embedding method for MPEG layer III audio steganography
Li et al. Audio steganography algorithm based on genetic algorithm for MDCT coefficient adjustment for AAC
Atoum New mp3 steganography data set
CN104184697B (zh) 一种音频指纹的提取方法及系统
CN104867496A (zh) 一种基于mdct量化系数的小值区的aac音频隐写和提取方法
Quan et al. Data hiding in MPEG compressed audio using wet paper codes
CN109346091B (zh) 一种基于等长熵码字替换的压缩音频自适应隐写方法和系统
CN109448744B (zh) 一种基于符号位自适应嵌入的mp3音频信息隐藏方法和系统
Vu et al. Data embedding in audio signal by a novel bit marking method

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190719