CN102123282A - 基于Wyner-Ziv视频编码系统的GOP层编码方法 - Google Patents

基于Wyner-Ziv视频编码系统的GOP层编码方法 Download PDF

Info

Publication number
CN102123282A
CN102123282A CN 201110057717 CN201110057717A CN102123282A CN 102123282 A CN102123282 A CN 102123282A CN 201110057717 CN201110057717 CN 201110057717 CN 201110057717 A CN201110057717 A CN 201110057717A CN 102123282 A CN102123282 A CN 102123282A
Authority
CN
China
Prior art keywords
gop
frame
coding
code check
key frame
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
CN 201110057717
Other languages
English (en)
Other versions
CN102123282B (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN 201110057717 priority Critical patent/CN102123282B/zh
Publication of CN102123282A publication Critical patent/CN102123282A/zh
Application granted granted Critical
Publication of CN102123282B publication Critical patent/CN102123282B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Compression Or Coding Systems Of Tv Signals (AREA)

Abstract

本发明公开了一种Wyner-Ziv视频编码系统中编码端基于率失真优化的图像组GOP层编码方法,主要解决现有码率分配技术不能灵活分配图像组GOP层码率分配比例以及不能应用于实际系统的问题。其编码过程是:在编码端估算码率分配比例因子k和失真D的关系;利用编码端重构的关键帧和估计的边信息估算满足用户对分布式视频编码DVC系统整体性能最优的最佳比例因子k*;用GOP内最佳比例因子k*分配关键帧码率Rk和WZ帧码率Rwz;用关键帧码率Rk和WZ帧码率Rwz编码当前GOP。本发明能灵活确定达到特定失真要求的GOP内码率分配的最佳比例,避免了速率分配的盲目性,可在编码端进行失真预估和码率分配,可用于实际的DVC系统。

Description

基于Wyner-Ziv视频编码系统的GOP层编码方法
技术领域
本发明属于视频编码技术领域,主要涉及分布式视频编码中基于率失真优化的无反馈编码速率控制中GOP(Group-of-Picture,图像组)层编码方法,可用于实用分布式视频编码系统和终端编码设备资源受限的视频通信系统。
背景技术
随着信息技术和互联网的飞速发展,多媒体已成为人们获取信息最主要的载体,多媒体业务将是未来无线宽带网络的关键业务。有数据显示,未来宽带网络上数据流量的80%将是以视频为主的多媒体数据流量。为了适应无线网络日益发展,越来越多的移动视频终端被用于多媒体通信中,如无线视频探测头、便携式摄像机等,以及广泛应用于传感器网络中的视频传感器。这些编码设备的处理能力和功耗等资源都十分有限,因此要求视频编码器简单易实现,且具有良好的抗误码性能和压缩效率。
传统的视频编码标准如MPEG、H.263和H.264等都采用非对称编码方式,编码器复杂度远大于解码器。另外,传统视频编码方法在编码端采用预测编码模式,网络传输导致的丢包或比特错误会造成编码端和解码端重构缓存的不匹配,从而出现误码扩散现象,容错能力较差,显然无法很好的适应无线网络下要求误码率高、终端设备资源受限等应用。
为了解决传统视频编码在视频通信中遇到的编码器复杂、抗误码能力差等问题,国内外研究一种新的视频压缩技术——分布式视频编码方法DVC,其中Wyner-Ziv视频编码是目前国内外分布式视频编码的研究热点。分布式视频编码系统采用帧内编码-帧间解码的结构,将编码视频序列分为关键帧和Wyner-Ziv帧,对关键帧进行传统的H.263或H.264帧内编码,而对Wyner-Ziv帧进行Wyner-Ziv视频编码。其中Wyner-Ziv帧的编码码率由它和边信息之间的统计相关性决定,由于原始Wyner-Ziv帧和边信息分别位于系统的编码侧和解码侧,这就使得在编码侧准确估计成功解码Wyner-Ziv帧所需要的比特数较为困难。为了使解码器能够正确解码,编码器需要传输足够多的校验比特,而传输过多的校验比特又会降低系统的压缩性能,所以码率控制在一定程度上决定了分布式编码系统的编码效率和应用环境。
现有的分布式视频编码系统绝大多数都采用有反馈信道的码率控制方法,编码得到的校验位根据解码侧的反馈请求分多次发送,直到解码成功,这样可以在一定程度上提高编码效率,但是大大增加了解码的复杂度和时延,不利于实时视频传输,同时使用反馈信道也使得分布式视频编码系统的应用环境受到限制;另外,现有的有反馈码率控制方法只是在编码侧简单设定关键帧的量化因子和Wyner-Ziv帧的量化矩阵,无法在带宽受限信道下有效进行编码速率控制。因此,无反馈的速率控制方法被广泛研究,以克服有反馈速率控制方法的缺点。现有的无反馈速率控制方法都是自顶向下设计,一般包括GOP层码率控制、帧层码率控制,系数带级码率控制和比特面级码率控制。使用无反馈速率控制方法法虽然能在基本上不增加或少增加编码端复杂度的基础上,取消有反馈速率控制中需要的反馈信道,增强了DVC系统的实用性,但在传统的无反馈速率控制方法中,由于其顶层GOP的码率控制只是根据一些简单的准则,如相邻图像相关性,以固定比例分配GOP内关键帧和WZ帧码率,而没有考虑这个分配比例对解码恢复序列质量的影响,因而导致码率分配具有盲目性,最终解码质量也无法达到整体最优。
发明内容
本发明的目的在于克服上述无反馈速率控制方法中已有码率分配方法的缺点,提出一种基于Wyner-Ziv视频编码系统的GOP层编码方法,以充分考虑GOP内关键帧和WZ帧的码率分配比例对解码恢复视频序列的质量影响,根据用户需求确定GOP内码率分配比例,避免码率分配的盲目性,提高解码图像整体质量。
实现本发明目的技术思路是在编码端通过估计GOP内关键帧和WZ帧的码率分配比例对解码恢复视频序列的质量影响,运用率失真优化理论,计算得到每一个GOP内最佳的码率分配比列,从而使得解码图像整体质量达到最优。具体实现步骤包括如下:
(1)确定关键帧码率在图像组GOP码率中所占比例因子k的定义域:k∈[ε,1]及比例因子k的变化步长step,ε取值为[0.01,0.3],step取值范围为[0.02,0.4],k的初始值k′=ε;
(2)根据当前GOP的总码率RGOP和当前比例因子k,计算关键帧码率Rk和WZ帧码率Rwz
Rk=RGOP·k
Rwz=RGOP-Rk
(3)根据关键帧码率Rk、WZ帧码率Rwz和率失真模型估计当前GOP内关键帧和WZ帧所需的量化因子Qp和量化矩阵Qm,并用量化因子Qp和量化矩阵Qm预编码当前GOP,得到当前GOP的相邻重构关键帧,对该相邻重构关键帧内插,得到边信息SI;
(4)根据相邻重构关键帧和边信息SI计算当前GOP内关键帧和边信息SI峰值信噪比PSNR的均值D1(k)=mean(PSNRkey,PSNRSI)和方差D2(k)=var(PSNRkey,PSNRSI),得到编码端估计的整体失真D(k)=αD1(k)+βD2(k),其中,mean和var分别为求均值和求方差运算,PSNRkey和PSNRSI分别代表关键帧和边信息的峰值信噪比,α为D1(k)的权重因子,β为D2(k)的权重因子,α和β取值范围均为[0,1],且α+β=1;
(5)判断比例因子k值是否超过定义域[ε,1]上界,如果是,则选择整体失真D(k)最小时对应的k作为最佳比例因子k*,执行步骤(6),否则,令k=k+step,返回步骤(2);
(6)根据量化因子Qp和量化矩阵Qm编码当前GOP,并判断当前GOP是否是最后一个,若是,则结束本次编码;否则执行步骤(7);
(7)重复执行步骤(1)到步骤(6),对剩余的GOP进行编码,直至所有的GOP编码完毕。
本发明与现有技术相比具有如下优点:
1)本发明由于在失真估计中灵活选取不同的权重因子α和β,进而计算得到不同的最佳分配比例因子k*,使DVC系统的总体性能满足用户对DVC系统中WZ帧、关键帧以及总体图像质量的不同要求,从而灵活确定达到这种特定的总体失真的最佳分配比例,避免了速率分配的盲目性。
2)本发明由于在编码端进行失真预估,在不增加或少增加编码端复杂度的情况下取消了DVC系统中的反馈信道,有利于DVC的实际应用;
3)本发明由于利用率失真优化理论分别计算每个GOP内速率分配最佳比例,确保了DVC整体的率失真性能达到最优。
附图说明
图1是本发明基于率失真优化的GOP层编码方法流程图;
图2是本发明Hall Monitor序列的比例因子与失真关系图;
图3是本发明Foreman序列的比例因子与失真关系图;
图4是本发明的Hall Monitor序列性能仿真图;
图5是本发明的Foreman序列性能仿真图。
具体实施方式
参照图1,本发明GOP层编码的具体实施步骤如下:
步骤1,确定比例因子k的定义域、变化步长及初值。
经过大量的实验发现,比例因子关键帧码率Rk在GOP总码率RGOP中所占比例
Figure BDA0000049654460000041
在DVC系统整体性能中起着关键作用,随着Rk在RGOP中的比例增大,即随着比例因子k的提高,DVC整体峰值信噪比PSNR提高,同时以方差表示的解码序列质量波动增大,反之亦然,如图2和图3所示,其中图2(a)是Hall序列比例因子k和PSNR的关系,图2(b)Hall序列比例因子k和方差的关系;图3(a)是Foreman序列比例因子k和PSNR的关系,图3(b)Hall序列比例因子k和方差的关系。
本发明根据图2和图3,确定比例因子k的定义域为:k∈[ε,1],k的变化步长为step,ε取值为[0.01,0.3],step取值范围为[0.02,0.4],k的初始值k′=ε。
步骤2,计算关键帧码率Rk和WZ帧码率Rwz
根据当前GOP的总码率RGOP和当前比例因子k,利用下式计算关键帧码率Rk和WZ帧码率Rwz
Rk=RGOP·k
Rwz=RGOP-Rk
步骤3,预编码当前GOP,并生成边信息SI。
根据关键帧码率Rk、WZ帧码率Rwz和率失真模型估计当前GOP内关键帧和WZ帧所需的量化因子Qp和量化矩阵Qm,并用量化因子Qp和量化矩阵Qm预编码当前GOP,得到当前GOP的相邻重构关键帧,对该相邻重构关键帧内插,得到边信息SI,具体步骤如下:
(3a)根据关键帧码率Rk,由率失真公式
Figure BDA0000049654460000042
计算得到量化因子Qp,
其中,m和n为率失真公式中两个不同的参量,取值范围为实数;
(3b)根据8个量化矩阵Qm和Wyner-Ziv视频编码方法,对WZ帧进行编码,得到8个预估计的WZ帧码率
Figure BDA0000049654460000051
其中i为量化矩阵的序号,取值为自然数1~8;
(3c)从(3b)中得到8个预估计的WZ帧码率
Figure BDA0000049654460000052
并从中选择与WZ帧码率Rwz差距最小的
Figure BDA0000049654460000053
所对应的量化矩阵Qm,作为当前WZ帧的量化矩阵;
(3d)根据步骤(3a)中计算得到的量化因子Qp,利用H.264/AVC INTRA编码方法第当前关键帧进行编码,得到重构关键帧;
(3e)根据量化矩阵Qm,利用Wyner-Ziv视频编码方法对当前WZ帧进行编码;
(3d)根据重构关键帧,利用均值内插算法或快速运动补偿内插算法,但不限于此类算法,计算得到边信息。
步骤4,在编码端进行失真估计。
根据相邻的重构关键帧和边信息SI,计算当前GOP内关键帧和边信息SI峰值信噪比PSNR的均值D1(k)和方差D2(k),从而计算得到编码端估计的整体失真D(k),如下式:
D1(k)=mean(PSNRkey,PSNRSI)
D2(k)=var(PSNRkey,PSNRSI)
D(k)=αD1(k)+βD2(k)
其中,mean和var分别为求均值和求方差运算,PSNRkey和PSNRSI分别代表关键帧和边信息的峰值信噪比,α为D1(k)的权重因子,β为D2(k)的权重因子,α和β取值范围均为[0,1],且α+β=1。
步骤5,判断当前比例因子k值是否越界。
判断比例因子k值是否超过定义域[ε,1]的上界,如果是,则选择整体失真D(k)最小时对应的k作为最佳比例因子k*,执行步骤(6);否则,令k=k+step,返回步骤(2);
步骤(6),编码图像组GOP。
(6a)用比例因子k中的最佳比例因子k*,计算关键帧码率Rk和WZ帧码率Rwz
Rk=RGOP·k*
Rwz=RGOP-Rk
其中,RGOP为当前GOP的总码率;
(6b)根据关键帧码率Rk,利用H.264/AVC INTRA编码方法,对当前关键帧进行编码;
(6c)根据WZ帧码率Rwz,利用Wyner-Ziv视频编码方法,对当前WZ帧进行编码。
(6d)判断当前GOP是否是最后一个,若是,则结束本次编码;否则执行步骤(7);
步骤7,重复执行步骤(1)到步骤(6),对剩余的GOP进行编码,直至所有的GOP编码完毕。
本发明的效果通过以下实验进一步说明:
1)实验条件
硬件环境:CPU INTEL Pentium(R),3.00GHZ,1.00GB内存;
软件环境:Windows XP(SP2),Microsoft Visual Studio.NET 2003;
GOP结构:关键帧,WZ帧,关键帧,WZ帧,......;
参考序列:Hall Monitor(测试帧数150帧,帧率15Hz),Foreman(测试帧数300帧,帧率15Hz);
图像大小:QCIF(176*144);
2)实验内容
本实验给出了在两种失真度量D1(α=0.2,β=0.8)和D2(α=0.8,β=0.2)下,采用本发明提出的基于率失真优化理论的GOP层编码方法得到的实验结果,如图4和图5。
在图4中,Hall Monitor图像在D1失真下图像平均质量为38.6dB,而此时的图像质量波动高达7dB;在D2失真下,图像质量波动仅2dB左右,平均质量相应为35.8dB。图5在D1和D2失真亦是如此。据此可以得出结论:采用本发明提供的率失真优化方法,可以根据用户对解码WZ帧、关键帧和总体图像质量的不同要求进行优化码率分配,得到符合要求的解码视频质量。

Claims (4)

1.一种基于Wyner-Ziv视频编码系统的GOP层编码方法,包括如下步骤:
(1)确定关键帧码率在图像组GOP码率中所占比例因子k的定义域:k∈[ε,1]及比例因子k的变化步长step,ε取值为[0.01,0.3],step取值范围为[0.02,0.4],k的初始值k′=ε;
(2)根据当前GOP的总码率RGOP和当前比例因子k,计算关键帧码率Rk和WZ帧码率Rwz
Rk=RGOP·k
Rwz=RGOP-Rk
(3)根据关键帧码率Rk、WZ帧码率Rwz和率失真模型估计当前GOP内关键帧和WZ帧所需的量化因子Qp和量化矩阵Qm,并用量化因子Qp和量化矩阵Qm预编码当前GOP,得到当前GOP的相邻重构关键帧,对该相邻重构关键帧内插,得到边信息SI;
(4)根据相邻重构关键帧和边信息SI计算当前GOP内关键帧和边信息SI峰值信噪比PSNR的均值D1(k)=mean(PSNRkey,PSNRSI)和方差D2(k)=var(PSNRkey,PSNRSI),得到编码端估计的整体失真D(k)=αD1(k)+βD2(k),其中,mean和var分别为求均值和求方差运算,PSNRkey和PSNRSI分别代表关键帧和边信息的峰值信噪比,α为D1(k)的权重因子,β为D2(k)的权重因子,α和β取值范围均为[0,1],且α+β=1;
(5)判断比例因子k值是否超过定义域[ε,1]上界,如果是,则选择整体失真D(k)最小时对应的k作为最佳比例因子k*,执行步骤(6),否则,令k=k+step,返回步骤(2);
(6)根据量化因子Qp和量化矩阵Qm编码当前GOP,并判断当前GOP是否是最后一个,若是,则结束本次编码;否则执行步骤(7);
(7)重复执行步骤(1)到步骤(6),对剩余的GOP进行编码,直至所有的GOP编码完毕。
2.根据权利要求1所述的GOP层编码方法,其中步骤(3)所述的根据关键帧码率Rk、WZ帧码率Rwz和率失真模型,估计当前GOP内关键帧和WZ帧所需的量化因子Qp和量化矩阵Qm,按如下步骤进行:
(2a)根据关键帧码率Rk,由公式计算得到量化因子Qp,其中,m和n为率失真模型中两个不同的参量,取值范围为实数;
(2b)根据8个量化矩阵Qm和Wyner-Ziv视频编码方法编码WZ帧,得到8个预估计的WZ帧码率
Figure FDA0000049654450000022
其中i为量化矩阵的序号,取值为自然数1~8;选择与WZ帧码率Rwz差距最小的预估计的WZ帧码率
Figure FDA0000049654450000023
所对应的量化矩阵,作为当前WZ帧的量化矩阵;
(2c)根据步骤(2a)中计算得到的量化因子Qp,利用H.264/AVC INTRA编码方法编码当前关键帧;根据步骤(2b)中计算得到的量化矩阵Qm,利用Wyner-Ziv视频编码方法编码当前WZ帧。
3.根据权利要求1所述的GOP层编码方法,其中步骤(3)所述的对相邻重构关键帧内插,是指用均值内插算法或快速运动补偿内插算法对相邻的两帧重构关键帧进行求均值或对相邻的两帧重构关键帧进行前后向低复杂度的运动搜索得到边信息。
4.根据权利要求1所述的GOP层编码方法,其中步骤(6)所述根据量化因子Qp和量化矩阵Qm编码当前GOP,按如下步骤进行:
(4a)用比例因子k中的最佳比例因子k*,计算关键帧码率Rk和WZ帧码率Rwz
Rk=RGOP·k*
Rwz=RGOP-Rk
其中,RGOP为当前GOP的总码率。
(4b)根据关键帧码率Rk,利用H.264/AVC INTRA编码方法,对当前关键帧进行编码;
(4c)根据WZ帧码率Rk,利用Wyner-Ziv视频编码方法,对当前WZ帧进行编码。
CN 201110057717 2011-03-10 2011-03-10 基于Wyner-Ziv视频编码系统的GOP层编码方法 Expired - Fee Related CN102123282B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110057717 CN102123282B (zh) 2011-03-10 2011-03-10 基于Wyner-Ziv视频编码系统的GOP层编码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110057717 CN102123282B (zh) 2011-03-10 2011-03-10 基于Wyner-Ziv视频编码系统的GOP层编码方法

Publications (2)

Publication Number Publication Date
CN102123282A true CN102123282A (zh) 2011-07-13
CN102123282B CN102123282B (zh) 2013-02-27

Family

ID=44251708

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110057717 Expired - Fee Related CN102123282B (zh) 2011-03-10 2011-03-10 基于Wyner-Ziv视频编码系统的GOP层编码方法

Country Status (1)

Country Link
CN (1) CN102123282B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105554500A (zh) * 2015-12-15 2016-05-04 华中科技大学 基于rcm译码符号不等分配的视频传输方法及系统
CN108401159A (zh) * 2017-02-07 2018-08-14 腾讯科技(深圳)有限公司 一种目标码率控制的方法以及电子设备
CN110800298A (zh) * 2018-08-21 2020-02-14 深圳市大疆创新科技有限公司 码率分配方法、码率控制方法、编码器和记录介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1589764A2 (en) * 2004-04-13 2005-10-26 Samsung Electronics Co., Ltd. Method and apparatus for supporting motion scalability
CN101146228A (zh) * 2007-10-24 2008-03-19 北京大学 一种基于率失真函数分析模型的svc平滑重建方法
CN101674472A (zh) * 2009-09-25 2010-03-17 上海大学 多视点视频编码的多级码率控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1589764A2 (en) * 2004-04-13 2005-10-26 Samsung Electronics Co., Ltd. Method and apparatus for supporting motion scalability
CN101146228A (zh) * 2007-10-24 2008-03-19 北京大学 一种基于率失真函数分析模型的svc平滑重建方法
CN101674472A (zh) * 2009-09-25 2010-03-17 上海大学 多视点视频编码的多级码率控制方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105554500A (zh) * 2015-12-15 2016-05-04 华中科技大学 基于rcm译码符号不等分配的视频传输方法及系统
CN105554500B (zh) * 2015-12-15 2018-08-21 华中科技大学 基于rcm译码符号不等分配的视频传输方法及系统
CN108401159A (zh) * 2017-02-07 2018-08-14 腾讯科技(深圳)有限公司 一种目标码率控制的方法以及电子设备
CN108401159B (zh) * 2017-02-07 2022-07-08 腾讯科技(深圳)有限公司 一种目标码率控制的方法以及电子设备
CN110800298A (zh) * 2018-08-21 2020-02-14 深圳市大疆创新科技有限公司 码率分配方法、码率控制方法、编码器和记录介质

Also Published As

Publication number Publication date
CN102123282B (zh) 2013-02-27

Similar Documents

Publication Publication Date Title
CN101835042B (zh) 基于无反馈速率控制的Wyner-Ziv视频编码系统及方法
CN103957389B (zh) 基于压缩感知的3g视频传输方法及系统
CN101232619B (zh) 嵌入帧内编码块的视频编码方法
CN100558170C (zh) 具有主动式缓存管理和复杂度控制功能的视频解码方法
CN103533359B (zh) 一种h.264码率控制方法
CN100574441C (zh) 一种感兴趣区域的率失真优化帧内刷新与码率分配方法
CN102186077B (zh) 基于Wyner-Ziv视频编码的Wyner-Ziv帧码率控制系统和方法
CN103634606B (zh) 视频编码方法和装置
CN106231320B (zh) 一种支持多机并行编码的联合码率控制方法及系统
CN105681797B (zh) 一种基于预测残差的dvc-hevc视频转码方法
CN106210721B (zh) 一种hevc快速码率转码方法
CN110933430B (zh) 二次编码优化方法
CN102572428B (zh) 面向多媒体传感网分布式编解码的边信息估计方法
CN103442228A (zh) 从h.264/avc标准到hevc标准的快速帧内转码方法及其转码器
CN102123282B (zh) 基于Wyner-Ziv视频编码系统的GOP层编码方法
CN102420987A (zh) 基于分层b帧结构的码率控制的自适应比特分配方法
CN102510496B (zh) 基于感兴趣区域的快速降尺寸转码方法
CN101127919A (zh) 一种视频序列的编码方法
CN101888547B (zh) 基于直接模式选择的h.264/avc快速转码方法及装置
CN103391439B (zh) 一种基于主动宏块隐藏的h.264/avc码率控制方法
CN101146229B (zh) 一种svc视频fgs优先级调度方法
CN102088608B (zh) 一种基于部分重建的可伸缩视频编码质量优化方法
CN103475878A (zh) 一种视频编码方法和编码器
Yang et al. Multimedia communication and scalable video coding
CN104796704B (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
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130227

Termination date: 20190310

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