CN104484851A - 自适应高容量的图像信息隐藏的方法 - Google Patents

自适应高容量的图像信息隐藏的方法 Download PDF

Info

Publication number
CN104484851A
CN104484851A CN201410812054.8A CN201410812054A CN104484851A CN 104484851 A CN104484851 A CN 104484851A CN 201410812054 A CN201410812054 A CN 201410812054A CN 104484851 A CN104484851 A CN 104484851A
Authority
CN
China
Prior art keywords
block
sub
image
information
self
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
Application number
CN201410812054.8A
Other languages
English (en)
Other versions
CN104484851B (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.)
Xihua University
Original Assignee
Xihua 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 Xihua University filed Critical Xihua University
Priority to CN201410812054.8A priority Critical patent/CN104484851B/zh
Publication of CN104484851A publication Critical patent/CN104484851A/zh
Application granted granted Critical
Publication of CN104484851B publication Critical patent/CN104484851B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Editing Of Facsimile Originals (AREA)
  • Facsimile Transmission Control (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种自适应高容量的图像信息隐藏的方法其步骤如下,输入载体图像,同时对载体图像进行分块处理,并将待嵌入的秘密信息转化为比特;对载体图像分块后的子块的利用情况进行判断,子块使用完时直接结束信息隐藏过程,相反则继续对剩余的子块进行信息隐藏操作;选取载体图像分块后的第i块子块的中心元素,定为Ci,并以Ci为圆心,将该子块旋参照Z×90°转到指定的角度值,其中Z为旋转参数;将转化后的比特秘密信息嵌入到该子块中,并将该子块以Ci为圆心,反方向旋转到起始位置;判断秘密信息嵌入的情况进行判断,直至完成载密图像的输出,信息隐藏方法结束。该算法不仅有效提高了隐藏的容量,并且其算法实现简单、计算量小。

Description

自适应高容量的图像信息隐藏的方法
技术领域
本发明涉及图像隐藏领域,特别涉及一种自适应高容量的图像信息隐藏的方法。
背景技术
高容量是信息隐藏研究的主要内容之一。近年来,对自适应嵌入的研究已成为信息隐藏领域研究的又一热点。图像信息隐藏技术是通过设计的隐藏算法把秘密信息嵌入在载体图像(Cover image)中,形成载密图像(Stego image)。基于数字图像的信息隐藏技术可分为空域和变换域信息隐藏技术。一个好的图像信息隐藏算法应该包含较高的容量(bpp)和好的图像质量。通常情况下,容量与图像质量是相互制约的。即,图像容量越大,其图像的质量就越差。如何既能够提高容量,又能够减小或者保持图像的质量,这也是信息隐藏设计者必须所面临的问题。针对高容量的信息隐藏,研究者提出一些算法,取得一定成果,但有待进一步提高。其中的空域信息隐藏技术通常采用LSB替换法,算法简单、复杂度低,是一种常用的信息隐藏技术,但其隐藏数据量非常有限,抗干扰性较差。随后学者提出的数字水印方法(LDT),进一步提高了容量,基于直方图的可逆信息隐藏方法(HMPD),通过修改直方图来隐藏信息,利用直方图和整数变化,提出的大容量可逆的水印方法(EGIT),以及自适应的信息隐藏方法(ARW),虽然进一步提高隐藏的容量,而获得较低的图像失真,但其容量和图像质量仍有待于提高;虽然LCIH方法提高了的数据嵌入容量和能够获得较好的隐藏图像质量。但是,LCIH算法的嵌入数据的容量和稳定性,还有待进一步提高。
发明内容
本发明要解决的是传统方式的多种图像隐藏方法在隐藏数据量、图象精确度和嵌入数据稳定性方面有缺陷的技术问题。
为解决上述技术问题,本发明所采取的技术方案是:提供一种自适应高容量的图像信息隐藏的方法,该方法的实现步骤如下:
步骤S1、输入载体图像,并对载体图像进行分块处理;将待嵌入的秘密信息转化为比特;
步骤S2、对载体图像分块后的子块的利用情况进行判断,子块使用完时直接结束信息隐藏过程,并在子块未使用完时进行步骤S3;
步骤S3、选取载体图像分块后的第i块子块的中心元素,定为Ci,并以Ci为圆心,将该子块旋参照Z×90°转到指定的角度值,其中Z为旋转参数;
步骤S4、将转化后的比特秘密信息嵌入到步骤S3得到的转动后的子块中,并将该子块以Ci为圆心,反方向旋转到起始位置;
步骤S5、判断秘密信息嵌入是否完成,若完成则输出载密图像,信息隐藏方法结束,否则返回步骤S2继续执行。
较优的,在上述技术方案中,对载体图像进行分块处理具体为:按照的比例分成m×n个子块,其中M×N为载体图像的像素大小。
较优的,在上述技术方案中,步骤S3中旋转参数Z的计算公式为k,其中表示取小于*的步骤S3中旋转参数Z为奇数时,子块以其中心元素Ci为圆心顺时针旋转Z×90°;Z为偶数时,子块以其中心元素Ci为圆心逆时针旋转Z×90°。
较优的,在上述技术方案中,将转化后的比特秘密信息嵌入到步骤S3得到的转动后的子块中,并将该子块以Ci为圆心,反方向旋转到起始位置具体为:
步骤S41、计算第i块中每个像素隐藏信息的量di,其计算公式为:
其中,|*|表示对数*取绝对值,C(i-1)和C(i)分别表示前后两个相邻子块的中心元素,表示取小于*的最大整数,l=1,2,3,…;
步骤S42、将嵌入第i块中的第j个像素的二进制秘密信息binj转换为十进制数mj,其计算公式为:
bj=(C(j,i))mod 2di
mj=bin2dec(binj)
其中,j=1,2,3,4,...8。
较优的,在上述技术方案中,步骤S41中di的确定具体为:在i=1时,di=l;在i>1时,根据步骤S41中的di得到最终di的修正式如下:
d i = 5 , d i ≥ 5 d i , otherwise .
较优的,在上述技术方案中,步骤S4中秘密信息嵌入的公式如下:
S ( j , i ) = C ( j , i ) - ( b j - m j ) , b j ≥ m j C ( j , i ) + ( m j - b j ) , otherwise
其中,j=1,2,3,4,...8。
较优的,在上述技术方案中,还包括秘密信息的提取算法,其步骤如下:
步骤S10、将载密图像按照步骤S1中的方式进行分块,并计算其旋转参数Z;
步骤S20、选取载密图像分块后的第i块像素子块,并选取其中心元素S(i),并根据Z×90°将该像素子块旋转到制定角度;
步骤S30、计算第j块中每个像素隐藏信息的量di,并根据S(j,i)、S(i)和di计算出mj,最终依据di将将mj转为成二进制数binj,其中,S(j,i)是第i块中的第j个像素;
步骤S40、重复上述步骤S30,直到全部恢复秘密信息比特流,最终从比特流中恢复秘密信息。
本发明提供的自适应高容量的图像信息隐藏的方法,结合了相邻子块的特性,以及自适应的和高容量的信息隐藏算法特点,通过子块的中心元以及子块随机旋转和自适应的嵌入信息,其相关参数作为密钥的一部分,提高了秘密信息隐藏的安全性。该方法具有较大的数据嵌入容量,同时隐藏图像具有较好的图像质量;并且该隐藏图像在不需要额外的信息的条件下,还能够从中正确提取嵌入数据,适用于多领域保密通讯信息的隐藏。
附图说明
下面结合附图对本发明作进一步说明:
图1是本发明提供的自适应高容量的图像信息隐藏的方法的流程。
具体实施方式
本发明提供的自适应高容量的图像信息隐藏的方法,该方法的实现步骤如下:
步骤S1、输入载体图像,并对载体图像进行分块处理;将待嵌入的秘密信息转化为比特;
步骤S2、对载体图像分块后的子块的利用情况进行判断,子块使用完时直接结束信息隐藏过程,并在子块未使用完时进行步骤S3;
步骤S3、选取载体图像分块后的第i块子块的中心元素,定为Ci,并以Ci为圆心,将该子块旋参照Z×90°转到指定的角度值,其中Z为旋转参数;
步骤S4、将转化后的比特秘密信息嵌入到步骤S3得到的转动后的子块中,并将该子块以Ci为圆心,反方向旋转到起始位置;
步骤S5、判断秘密信息嵌入是否完成,若完成则输出载密图像,信息隐藏方法结束,否则返回步骤S2继续执行。
作为一种可实施方式,对载体图像进行分块处理具体为:按照的比例分成m×n个子块,其中M×N为载体图像的像素大小。
作为一种可实施方式,步骤S3中旋转参数Z的计算公式为其中表示取小于*的步骤S3中旋转参数Z为奇数时,子块以其中心元素Ci为圆心顺时针旋转Z×90°;Z为偶数时,子块以其中心元素Ci为圆心逆时针旋转Z×90°。
作为一种可实施方式,将转化后的比特秘密信息嵌入到步骤S3得到的转动后的子块中,并将该子块以Ci为圆心,反方向旋转到起始位置具体为:
步骤S41、计算第i块中每个像素隐藏信息的量di,其计算公式为:
其中,|*|表示对数*取绝对值,C(i-1)和C(i)分别表示前后两个相邻子块的中心元素,表示取小于*的最大整数,l=1,2,3,…;
步骤S42、将嵌入第i块中的第j个像素的二进制秘密信息binj转换为十进制数mj,其计算公式为:
bj=(C(j,i))mod 2di
mj=bin2dec(binj)
其中,j=1,2,3,4,...8。
作为一种可实施方式,步骤S41中di的确定具体为:在i=1时,di=l;在i>1时,根据步骤S41中的di得到最终di的修正式如下:
d i = 5 , d i ≥ 5 d i , otherwise .
作为一种可实施方式,步骤S4中秘密信息嵌入的公式如下:
S ( j , i ) = C ( j , i ) - ( b j - m j ) , b j ≥ m j C ( j , i ) + ( m j - b j ) , otherwise
其中,j=1,2,3,4,...8。
作为一种可实施方式,还包括秘密信息的提取算法,其步骤如下:
步骤S10、将载密图像按照步骤S1中的方式进行分块,并计算其旋转参数Z;
步骤S20、选取载密图像分块后的第i块像素子块,并选取其中心元素S(i),并根据Z×90°将该像素子块旋转到制定角度;
步骤S30、计算第j块中每个像素隐藏信息的量di,并根据S(j,i)、S(i)和di计算出mj,最终依据di将将mj转为成二进制数binj,其中,S(j,i)是第i块中的第j个像素;
步骤S40、重复上述步骤S30,直到全部恢复秘密信息比特流,最终从比特流中恢复秘密信息。
针对传统方式的图像隐藏方法,以Lena图像为例,3×3分块为标准,对LDT、HMPD、EGIT、ARW、LCIH和本申请的自适应高容量的图像信息隐藏的方法,(此处缩写为AHCA),各方法的容量及隐藏图像的质量进行比较,其中bpp为容量,其对比结果如下:
从表中可以看出,本申请的自适应高容量的图像信息隐藏的方法AHCA,在图像隐藏的容量及隐藏图像的质量方面具有更高的优越性,很好的解决了传统方式的图像隐藏信息的方法在隐藏数据量、图像精确度和嵌入数据稳定性方面有缺陷的技术问题。
本发明提供的自适应高容量的图像信息隐藏的方法,结合了相邻子块的特性,以及自适应的和高容量的信息隐藏算法特点,通过子块的中心元以及子块随机旋转和自适应的嵌入信息,其相关参数作为密钥的一部分,提高了秘密信息隐藏的安全性。该方法具有较大的数据嵌入容量,同时隐藏图像具有较好的图像质量;并且该隐藏图像在不需要额外的信息的条件下,还能够从中正确提取嵌入数据,适用于多领域保密通讯信息的隐藏。
上述实施方式旨在举例说明本发明可为本领域专业技术人员实现或使用,对上述实施方式进行修改对本领域的专业技术人员来说将是显而易见的,故本发明包括但不限于上述实施方式,任何符合本权利要求书或说明书描述,符合与本文所公开的原理和新颖性、创造性特点的方法、工艺、产品,均落入本发明的保护范围之内。

Claims (8)

1.一种自适应高容量的图像信息隐藏的方法,其特征在于,该方法的实现步骤如下:
步骤S1、输入载体图像,并对所述载体图像进行分块处理;将待嵌入的秘密信息转化为比特;
步骤S2、对所述载体图像分块后的子块的利用情况进行判断,所述子块使用完时直接结束信息隐藏过程,并在所述子块未使用完时进行步骤S3;
步骤S3、选取所述载体图像分块后的第i块子块的中心元素,定为Ci,并以Ci为圆心,将该子块旋参照Z×90°转到指定的角度值,其中Z为旋转参数;
步骤S4、将转化后的比特秘密信息嵌入到步骤S3得到的转动后的子块中,并将该子块以Ci为圆心,反方向旋转到起始位置;
步骤S5、判断所述秘密信息嵌入是否完成,若完成则输出所述载密图像,信息隐藏方法结束,否则返回步骤S2继续执行。
2.根据权利要求1所述的自适应高容量的图像信息隐藏的方法,其特征在于,所述对所述载体图像进行分块处理具体为:按照的比例分成m×n个子块,其中M×N为所述载体图像的像素大小。
3.根据权利要求1所述的自适应高容量的图像信息隐藏的方法,其特征在于,步骤S3中旋转参数Z的计算公式为其中表示取小于*的最大整数,k=2,3,4。
4.根据权利要求1至3任一项所述的自适应高容量的图像信息隐藏的方法,其特征在于,步骤S3中旋转参数Z为奇数时,所述子块以其中心元素Ci为圆心顺时针旋转Z×90°;Z为偶数时,所述子块以其中心元素Ci为圆心逆时针旋转Z×90°。
5.根据权利要求1所述的自适应高容量的图像信息隐藏的方法,其特征在于,所述将转化后的比特秘密信息嵌入到步骤S3得到的转动后的子块中,并将该子块以Ci为圆心,反方向旋转到起始位置具体为:
步骤S41、计算第i块中每个像素隐藏信息的量di,其计算公式为:
其中,|*|表示对数*取绝对值,C(i-1)和C(i)分别表示前后两个相邻子块的中心元素,表示取小于*的最大整数,l=1,2,3,…;
步骤S42、将嵌入第i块中的第j个像素的二进制秘密信息binj转换为十进制数mj,其计算公式为:
bj=(C(j,i))mod 2di
mj=bin2dec(binj)
其中,bj为平衡因子,j=1,2,3,4,...8。
6.根据权利要求5所述的自适应高容量的图像信息隐藏的方法,其特征在于,步骤S41中di的确定具体为:在i=1时,di=l;在i>1时,根据步骤S41中的di得到最终di的修正式如下:
d i = 5 , d i ≥ 5 d i , otherwise .
7.根据权利要求5所述的自适应高容量的图像信息隐藏的方法,其特征在于,步骤S4中秘密信息嵌入的公式如下:
S ( j , i ) = C ( j , i ) - ( b j - m j ) , b j ≥ m j C ( j , i ) + ( m i - b j ) , otherwise
其中,j=1,2,3,4,...8。
8.根据权利要求6所述的自适应高容量的图像信息隐藏的方法,其特征在于,还包括秘密信息的提取算法,其步骤如下:
步骤S10、将所述载密图像按照步骤S1中的方式进行分块,并计算其旋转参数Z;
步骤S20、选取所述载密图像分块后的第i块像素子块,并选取其中心元素S(i),并根据Z×90°将该像素子块旋转到制定角度;
步骤S30、计算第j块中每个像素隐藏信息的量di,并根据S(j,i)、S(i)和di计算出mj,最终依据di将将mj转为成二进制数binj,其中,S(j,i)是第i块中的第j个像素;
步骤S40、重复上述步骤S30,直到全部恢复秘密信息比特流,最终从比特流中恢复秘密信息。
CN201410812054.8A 2014-12-23 2014-12-23 自适应高容量的图像信息隐藏的方法 Expired - Fee Related CN104484851B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410812054.8A CN104484851B (zh) 2014-12-23 2014-12-23 自适应高容量的图像信息隐藏的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410812054.8A CN104484851B (zh) 2014-12-23 2014-12-23 自适应高容量的图像信息隐藏的方法

Publications (2)

Publication Number Publication Date
CN104484851A true CN104484851A (zh) 2015-04-01
CN104484851B CN104484851B (zh) 2017-10-03

Family

ID=52759391

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410812054.8A Expired - Fee Related CN104484851B (zh) 2014-12-23 2014-12-23 自适应高容量的图像信息隐藏的方法

Country Status (1)

Country Link
CN (1) CN104484851B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105488822A (zh) * 2015-12-15 2016-04-13 西华大学 基于ambtc算法的可逆图像隐藏算法
CN108881276A (zh) * 2018-07-10 2018-11-23 北京贞宇科技有限公司 数据加密传输方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050053259A1 (en) * 2003-09-05 2005-03-10 Wataru Asano Digital watermark detection method and apparatus
CN101582157A (zh) * 2009-03-06 2009-11-18 中山大学 一种基于相邻像素差异的自适应空域隐写方法
US20100014705A1 (en) * 2003-11-19 2010-01-21 Gustafson Ammon E Optimized Digital Watermarking Functions for Streaming Data
CN103279917A (zh) * 2013-05-30 2013-09-04 上海理工大学 变换域抗几何攻击灰度图像数字水印技术
CN104143174A (zh) * 2014-08-15 2014-11-12 陕西师范大学 基于Tangram算法和2维双尺度矩形映射的图像伪装及重构方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050053259A1 (en) * 2003-09-05 2005-03-10 Wataru Asano Digital watermark detection method and apparatus
US20100014705A1 (en) * 2003-11-19 2010-01-21 Gustafson Ammon E Optimized Digital Watermarking Functions for Streaming Data
CN101582157A (zh) * 2009-03-06 2009-11-18 中山大学 一种基于相邻像素差异的自适应空域隐写方法
CN103279917A (zh) * 2013-05-30 2013-09-04 上海理工大学 变换域抗几何攻击灰度图像数字水印技术
CN104143174A (zh) * 2014-08-15 2014-11-12 陕西师范大学 基于Tangram算法和2维双尺度矩形映射的图像伪装及重构方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HENGFU YANG,ET AL.: "A high-capacity image data hiding scheme using adaptive LSB substitution", 《RADIOENGINEERING》 *
唐明伟等: "一种高效的随机分块图像隐写算法", 《计算机工程与应用》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105488822A (zh) * 2015-12-15 2016-04-13 西华大学 基于ambtc算法的可逆图像隐藏算法
CN105488822B (zh) * 2015-12-15 2019-06-21 西华大学 基于ambtc算法的可逆图像隐藏算法
CN108881276A (zh) * 2018-07-10 2018-11-23 北京贞宇科技有限公司 数据加密传输方法

Also Published As

Publication number Publication date
CN104484851B (zh) 2017-10-03

Similar Documents

Publication Publication Date Title
Duan et al. A new high capacity image steganography method combined with image elliptic curve cryptography and deep neural network
Yao et al. Improved dual-image reversible data hiding method using the selection strategy of shiftable pixels' coordinates with minimum distortion
Arham et al. Multiple layer data hiding scheme based on difference expansion of quad
Peng et al. Adaptive reversible data hiding scheme based on integer transform
CN110163801B (zh) 一种图像超分辨和着色方法、系统及电子设备
CN109685716B (zh) 一种基于高斯编码反馈的生成对抗网络的图像超分辨率重建方法
CN111292265A (zh) 一种基于生成式对抗神经网络的图像修复方法
CN103413269A (zh) 一种图像隐写方法和秘密信息提取方法
Sakshi et al. Least Significant Bit Steganography for Text and Image hiding
CN104899835A (zh) 基于盲模糊估计与锚定空间映射的图像超分辨处理方法
Lee et al. Reversible data hiding using Lagrange interpolation for prediction-error expansion embedding
Wang et al. HidingGAN: High capacity information hiding with generative adversarial network
CN106875324B (zh) 基于sbde的无损图像信息隐藏方法
Wang et al. A robust image watermarking algorithm using SVR detection
CN104484851A (zh) 自适应高容量的图像信息隐藏的方法
Li et al. Privacy protection method based on multidimensional feature fusion under 6G networks
CN105488822A (zh) 基于ambtc算法的可逆图像隐藏算法
CN104766263A (zh) 基于四元数Legendre矩校正的彩色图像水印嵌入及检测方法
Shi et al. A Blind Watermarking Technique for Color Image based on SVD with Circulation.
CN109285106A (zh) 一种基于slt-svd的医学影像的可逆还原方法
CN109920014B (zh) 3d网格模型隐写方法
Yi et al. Improved reduced difference expansion based reversible data hiding scheme for digital images
Zhu et al. Image sanitization in online social networks: A general framework for breaking robust information hiding
CN116188652A (zh) 一种基于双尺度循环生成对抗的人脸灰度图像着色方法
Wang et al. A novel (n, n) secret image sharing scheme based on Sudoku

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Tang Mingwei

Inventor after: Hu Jinrong

Inventor after: Cheng Xuwen

Inventor after: Guandengshui

Inventor after: Fang Yicheng

Inventor after: Rao Wei

Inventor after: Wang Chun

Inventor after: Zhang Yiwei

Inventor after: Wang Liang

Inventor after: Zeng Chengke

Inventor after: Chen Xiaoliang

Inventor after: Song Wen

Inventor after: Hu Jie

Inventor after: Du Yajun

Inventor after: Huang Zexia

Inventor after: Xia Meichen

Inventor after: Zhao Chengfang

Inventor before: Zeng Chengke

Inventor before: Cheng Xuwen

Inventor before: Guandengshui

Inventor before: Fang Yicheng

Inventor before: Rao Wei

Inventor before: Wang Chun

Inventor before: Zhang Yiwei

Inventor before: Wang Liang

Inventor before: Chen Xiaoliang

Inventor before: Song Wen

Inventor before: Hu Jie

Inventor before: Du Yajun

Inventor before: Huang Zexia

Inventor before: Xia Meichen

Inventor before: Zhao Chengfang

Inventor before: Hu Jinrong

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: CENG SHENGKE CHEN XIAOLIANG SONG WEN HU JIE DU YAJUN HUANG ZEXIA XIA MEICHEN ZHAO CHENGFANG HU JINRONG CHENG XUWEN GUAN DENGSHUI FANG YICHENG RAO WEI WANG CHUN ZHANG YIWEI WANG LIANG TO: TANG MINGWEI CENG SHENGKE CHEN XIAOLIANG SONG WEN HU JIE DU YAJUN HUANG ZEXIA XIA MEICHEN ZHAO CHENGFANG HU JINRONG CHENG XUWEN GUAN DENGSHUI FANG YICHENG RAO WEI WANG CHUN ZHANG YIWEI WANG LIANG

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171003

Termination date: 20211223

CF01 Termination of patent right due to non-payment of annual fee