CN1890983A - 改进的内插像素值的计算 - Google Patents

改进的内插像素值的计算 Download PDF

Info

Publication number
CN1890983A
CN1890983A CNA2004800356241A CN200480035624A CN1890983A CN 1890983 A CN1890983 A CN 1890983A CN A2004800356241 A CNA2004800356241 A CN A2004800356241A CN 200480035624 A CN200480035624 A CN 200480035624A CN 1890983 A CN1890983 A CN 1890983A
Authority
CN
China
Prior art keywords
pixels
pixel
location
integer
value
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
CNA2004800356241A
Other languages
English (en)
Other versions
CN100589580C (zh
Inventor
G·比约恩特盖得
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.)
Cisco Systems International SARL
Original Assignee
Tandberg Telecom AS
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 Tandberg Telecom AS filed Critical Tandberg Telecom AS
Publication of CN1890983A publication Critical patent/CN1890983A/zh
Application granted granted Critical
Publication of CN100589580C publication Critical patent/CN100589580C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/50Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using predictive coding
    • H04N19/503Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using predictive coding involving temporal prediction
    • H04N19/51Motion estimation or motion compensation
    • H04N19/523Motion estimation or motion compensation with sub-pixel accuracy

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Color Television Image Signal Generators (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明涉及一种视频压缩系统,并且尤其涉及数字视频系统中的压缩/解压。本发明公开了一种计算视频图像中内插像素位置的值的可替换方法。

Description

改进的内插像素值的计算
技术领域
本发明涉及一种视频压缩系统,并且尤其涉及数字视频系统的压缩/解压。
背景技术
在一些应用中采用了实时动态图像发送,例如视频会议、网络会议、TV广播以及视频电话等等。
然而,由于数字视频通常通过以8比特(1字节)表示图像中的每个像素来描述,所以表示动态图像需要大量信息。这些未压缩的视频数据导致了很大的比特量,并且由于带宽有限而不能在传统的通信网络和传输线路上实时传输。
因此,实时视频传输需要很大程度的数据压缩。然而数据压缩可能与图像质量平衡。因此,为了开发允许在带宽有限的数据连接上进行高质量视频的实时传输的压缩技术而进行了大量研究。
在视频压缩系统中,主要目标是以尽可能少的容量表示视频信息。容量定义为以比特为单位,表示为常数值或者比特/时间单位。在这两种情况下,主要目标均是减少比特数量。
在MPEG*和H.26*标准中描述了最常见的视频编码方法,这些方法均使用基于先前编码和解码图像的预测的块。视频数据在发送前经历四个主要过程:预测、变换、量化和熵编码。
预测过程很大程度上减少了待发送的视频序列中各个图像所需的比特量。它利用了序列中一部分与序列中其他部分的相似性。由于预测器部分对于编码器和解码器都是已知的,因此仅需要传输差值。这种差值的表示通常需要少得多的容量。预测主要基于来自先前重建的图像的图像内容,其中内容位置由运动向量定义。
在典型的视频序列中,当前块M的内容可能与先前解码的图像中对应块的内容相似。如果从上次解码图像以后没有发生改变,M的内容将与先前解码图像中相同位置的块相同。在其他情况下,图像中的目标可能移动,从而M的内容与先前解码图像中不同位置的块更加相似。这种移动可以通过运动向量(V)表示。例如,运动向量(3;4)表示M的内容在上次解码图像之后向左移动了3个像素并且向上移动了4个像素。
通过执行运动搜索来确定与块相关的运动向量。通过将块的内容与不同空间偏移的先前图像中的块连续比较来进行所述搜索。将与具有和当前块最佳匹配的比较块相关的相对于当前块的偏移确定为关联运动向量。在H.262、H.263、MPEG1、MPEG2中,对相同概念进行扩展使得运动向量也可以附带1/2像素值。向量分量5.5则表示运动为5个像素和6个像素的中值。更具体的,通过计算表示运动5的像素和表示运动6的像素之间的平均值来获得预测值。由于对两个像素操作以获得其间像素的预测值,这被称为2阶滤波器。这种运动向量通常称为具有分数像素分辨率或者分数运动向量。所有滤波器操作可以通过脉冲响应来定义。对两个像素进行平均的操作可以表示为脉冲响应(1/2,1/2)。类似地,对4个像素进行平均表示脉冲响应(1/4,1/4,1/4,1/4)。
发明内容
本发明的一个目标是提供一种在视频编码和解码中用于计算视频图像中位于整数像素位置之间的内插像素位置的值的方法,可以克服或者减少现有技术的缺陷。
本发明的一个特定目标是提供一种需要更小的计算容量同时保持或者改善性能和质量的方法。
为了实现上述和进一步的目标,本发明提供了一种在视频编码和解码中用于计算视频图像中位于整数像素位置之间的内插像素位置的像素值的方法,其中与整数像素位置水平或者垂直排列的一组内插像素位置的像素值是通过在相关的水平或者垂直排列的内插像素位置上应用第一离散滤波器而计算的,并且其中未与整数像素位置水平或者垂直排列的内插的第一像素位置的像素值是通过在第二和第三像素位置上应用第二离散滤波器而计算的,所述第二和第三像素位置为整数像素位置或者与所述内插的第一像素位置排列在一起的内插像素位置。
在一个优选实施例中,所述第二离散滤波器为2阶滤波器。
在一个优选实施例中,所述第一离散滤波器为对称6阶滤波器。
第二离散2阶滤波器优选地为对所述第二和第三像素位置的像素值进行平均的对称2阶滤波器。
第二和第三像素位置优选地分别位于第一内插像素位置的两侧并且相对于第一内插像素位置具有相同距离。
根据H.264/AVC编码标准,所述方法的上述实施例中的任何一者可以优选地在像素运动补偿过程中使用。
附图说明
通过参考下面的详细描述和附图可以更好地理解本发明,其中:
图1显示了根据H.264/AVC标准的整数和分数像素位置的示例;以及
图2显示了通过根据本发明的方法计算内插像素值时整数和分数像素位置的示例。
具体实施方式
下面通过描述优选实施例来描述本发明。然而,本领域技术人员应当理解,本发明范围内的其他应用和修改在所附权利要求书中定义。
近来ITU和ISO/IEC联合开发了一种新的视频压缩标准。在两个标准化实体中公共标准的正式名称为“ITU-T Recommendation H.264”和“ISO/IECMPEG-4(部分10)Advanced Video Coding”。所述公共标准在下面将被称为H.264/AVC。
图1显示了根据H.264/AVC标准的整数和分数像素位置的示例。
在H.264/AVC中,在运动分辨率和各个内插使用的像素数量方面对编码方法进行了改进。该方法使用高达1/4像素精度的运动补偿预测。甚至定义了1/8像素精度,但是没有包含在任何文件中。图1中显示了整数和分数像素位置的示例。为了简化起见,仅显示位置A、E、U和Y之间的内插。
位置A、E、U和Y表示整数像素位置。位置A”、E’、A’和E”表示A-E直线上其他的整数位置。
位置c、k、m、o和w表示半像素位置。这些位置中的内插值可以通过使用在整数像素值上操作的脉冲响应为(1/32,-5/32,20/32,20/32,-5/32,1/32)的6阶滤波器而获得。例如,c可以通过下式计算:
c=1/32·A″-5/32·E′+20/32·A+20/32·E-5/32·A′+1/32·E″
该滤波器适当地在水平或者垂直方向上操作。而且为了获得m的值,该滤波器不在整数值上操作,而是在其他方向上的已经内插的值上操作。以上描绘的方块中的剩下位置通过对各个整数和相邻半像素位置进行平均(2阶滤波器)而获得:
b=(A+c)/2,d=(c+E)/2,f=(A+k)/2,g=(c+k)/2,h=(c+m)/2,i=(c+o)/2,j=(E+o)/2
l=(k+m)/2,n=(m+o)/2,p=(U+k)/2,q=(k+w)/2,r=(m+w)/2,s=(w+o)/2,t=(Y+o)/2
v=(U+w)/2,x=(w+Y)/2
如上所述,因为位置m的计算必须在两个方向上进行计算,所以它比其他1/2像素位置需要更多计算容量。当如上所述使用6阶滤波器计算1/2位置时,首先在水平方向上应用6阶滤波器操作以得到m的垂直直线中的1/2像素位置,然后这些位置经过垂直方向上的6阶滤波器操作以确定m的值(反之亦然)。因为m是未与任何整数像素位置水平或者垂直排列的唯一的1/2像素位置,所以这种两次6阶操作对于2x2块中的位置m是仅有的。计算该位置所需的资源通常大于计算其他1/2像素位置所需资源的两倍。这种增加的复杂度进一步引入到现有技术中依赖于m的值的h、l、n、r的计算中。
本发明通过提出对m值以及至少某些1/4像素位置的值的可替换推导而避免了这种增加的复杂度。下面通过一个示例实施例描述本发明的原理。该示例使用与先前描述中相同的字母以表示2×2块中的整数位置和相关的1/2和1/4像素内插。
假定已经计算了与整数像素位置水平或者垂直排列的1/2像素位置。这些位置特别地表示为c、k、o和w。它们通过如上所述的6阶滤波器计算,但是也可以应用其他滤波器和计算方法。
m的值是通过在两个邻近的整数像素位置上应用2阶滤波器操作而计算出的,而不是通过在与m水平或者垂直排列的1/2像素位置上应用完全的滤波器操作计算出来。在此示例中,2阶滤波器操作对所述两个位置的值进行平均,即所述2阶滤波器具有脉冲响应(1/2,1/2)。在m的各个相对侧可以找到所述两个邻近位置。由于没有整数像素位置相对于m水平或者垂直排列,最近的整数像素位置在对角方向上,m表示为下式:
                          m=(E+U)/2
根据本发明,所涉及的1/4像素位置即在现有技术中从m得到的1/4像素位置同样通过对应方式得到,从而所有像素位置通过A、E、U、Y、c、k、o和/或w直接得到。通过这种方式,确定1/4像素位置值和特定情况的1/2像素位置需要更少的中间计算。这会减少计算的总体复杂度。
根据本发明优选实施例的一组推荐的表达式如下:
b=(A+c)/2,d=(c+E)/2,f=(A+k)/2,g=(c+k)/2,h=(A+o)/2,i=(c+o)/2,j=(E+o)/2
l=(c+U)/2,n=(E+w)/2,p=(U+k)/2,q=(k+w)/2,r=(k+Y)/2,s=(w+o)/2,t=(Y+o)/2
v=(U+w)/2,x=(w+Y)/2,m=(E+U)/2
本发明并不限于上述示例实施例中使用的对称2阶滤波器。基本原理是通过对整数像素位置和与整数像素位置水平或者垂直排列的1/2像素位置进行倾斜的滤波器操作的方式,来计算所有的1/4像素位置和特定情况的1/2像素位置的值。通过选择最邻近的已经计算的1/2和/或整数像素位置可以实现与现有技术计算相同的精度,所有这些选择的像素与待计算的像素位置位于同一斜线上。
对上述的示例实施例的某些表达式在图2中显示了这条斜线。粗体字母表示在(a)-(e)中各个情况中计算哪个像素。
因此图2显示了通过根据本发明的方法计算内插像素值时整数和分数像素位置的示例。
在图2(a)中,g表示的1/4像素位置的像素值通过表达式g=(c+k)/2计算。斜线穿过1/4像素位置c和1/4像素位置k。
在图2(b)中,i表示的1/4像素位置的像素值通过表达式i=(c+o)/2计算。斜线穿过1/4像素位置c和1/4像素位置o。
在图2(c)中,r表示的1/4像素位置的像素值通过表达式r=(k+Y)/2计算。斜线穿过1/4像素位置k和1/4像素位置Y。
在图2(d)中,s表示的1/4像素位置的像素值通过表达式s=(w+o)/2计算。斜线穿过1/4像素位置w和1/4像素位置o。
在图2(e)中,m表示的1/4像素位置的像素值通过表达式m=(E+U)/2计算。斜线穿过1/4像素位置E和1/4像素位置U。
基于上述示例描述可以对本发明进行很多修改和变化。本发明的范围由所附权利要求书定义。

Claims (9)

1.一种在视频编码和解码中用于计算视频图像中位于整数像素位置之间的内插像素位置的像素值的方法,其中与整数像素位置水平或者垂直排列的一组内插像素位置的像素值是通过在相关的水平或者垂直排列的内插像素位置上应用第一离散滤波器而计算的,
其特征在于所述方法包括通过在第二和第三像素位置上应用第二离散滤波器来计算未与整数像素位置水平或者垂直排列的内插的第一像素位置的像素值,所述第二和第三像素位置为整数像素位置或者与所述内插的第一像素位置排列在一起的内插像素位置。
2.根据权利要求1所述的方法,其中,所述第二离散滤波器为2阶滤波器。
3.根据权利要求2所述的方法,其中,所述第一离散滤波器为对称6阶滤波器。
4.根据权利要求2或3所述的方法,其中,所述第二离散2阶滤波器为对所述第二和第三像素位置的像素值进行平均的对称2阶滤波器。
5.根据权利要求4所述的方法,其中,所述第二和第三像素位置分别位于所述第一内插像素位置的两侧并且相对于所述第一内插像素位置具有相同距离。
6.根据权利要求2或3所述的方法,在根据H.264/AVC编码标准的像素运动补偿过程中使用。
7.根据权利要求4所述的方法,在根据H.264/AVC编码标准的像素运动补偿过程中使用。
8.根据权利要求5所述的方法,在根据H.264/AVC编码标准的像素运动补偿过程中使用。
9.根据权利要求1至5中任何一项所述的方法,在根据H.264/AVC编码标准的像素运动补偿过程中使用。
CN200480035624A 2003-12-05 2004-10-01 改进的内插像素值的计算方法 Expired - Fee Related CN100589580C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
NO20035444 2003-12-05
NO20035444A NO320114B1 (no) 2003-12-05 2003-12-05 Forbedret utregning av interpolerte pixelverdier

Publications (2)

Publication Number Publication Date
CN1890983A true CN1890983A (zh) 2007-01-03
CN100589580C CN100589580C (zh) 2010-02-10

Family

ID=30439628

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200480035624A Expired - Fee Related CN100589580C (zh) 2003-12-05 2004-10-01 改进的内插像素值的计算方法

Country Status (6)

Country Link
US (1) US7616689B2 (zh)
EP (1) EP1695560B1 (zh)
JP (1) JP4761390B2 (zh)
CN (1) CN100589580C (zh)
NO (1) NO320114B1 (zh)
WO (1) WO2005055615A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102265619A (zh) * 2008-09-23 2011-11-30 高通股份有限公司 切换式内插滤波器中的偏移计算

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9330060B1 (en) 2003-04-15 2016-05-03 Nvidia Corporation Method and device for encoding and decoding video image data
US8660182B2 (en) 2003-06-09 2014-02-25 Nvidia Corporation MPEG motion estimation based on dual start points
US8731071B1 (en) 2005-12-15 2014-05-20 Nvidia Corporation System for performing finite input response (FIR) filtering in motion estimation
US8724702B1 (en) 2006-03-29 2014-05-13 Nvidia Corporation Methods and systems for motion estimation used in video coding
US8660380B2 (en) 2006-08-25 2014-02-25 Nvidia Corporation Method and system for performing two-dimensional transform on data value array with reduced power consumption
US8756482B2 (en) 2007-05-25 2014-06-17 Nvidia Corporation Efficient encoding/decoding of a sequence of data frames
US9118927B2 (en) * 2007-06-13 2015-08-25 Nvidia Corporation Sub-pixel interpolation and its application in motion compensated encoding of a video signal
US8873625B2 (en) 2007-07-18 2014-10-28 Nvidia Corporation Enhanced compression in representing non-frame-edge blocks of image frames
US20100296587A1 (en) * 2007-10-05 2010-11-25 Nokia Corporation Video coding with pixel-aligned directional adaptive interpolation filters
US8416861B2 (en) * 2007-10-14 2013-04-09 Nokia Corporation Fixed-point implementation of an adaptive image filter with high coding efficiency
WO2010001562A1 (ja) * 2008-07-02 2010-01-07 株式会社ルネサステクノロジ フィルタ処理装置及び半導体装置
WO2010063881A1 (en) * 2008-12-03 2010-06-10 Nokia Corporation Flexible interpolation filter structures for video coding
US8666181B2 (en) 2008-12-10 2014-03-04 Nvidia Corporation Adaptive multiple engine image motion detection system and method
NO330275B1 (no) * 2008-12-19 2011-03-21 Tandberg Telecom As Fremgangsmate i en videokodings-/-dekodingsprosess
EP2442569A1 (en) * 2009-06-09 2012-04-18 Sony Corporation Image processing device and method
KR101682147B1 (ko) * 2010-04-05 2016-12-05 삼성전자주식회사 변환 및 역변환에 기초한 보간 방법 및 장치
CN105069747B (zh) * 2015-07-13 2018-07-20 中国科学院遥感与数字地球研究所 图像插值方法及装置

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0294958B1 (en) 1987-06-09 1995-08-23 Sony Corporation Motion compensated interpolation of digital television images
US5870497A (en) * 1991-03-15 1999-02-09 C-Cube Microsystems Decoder for compressed video signals
US5638128A (en) * 1994-11-08 1997-06-10 General Instrument Corporation Of Delaware Pixel interpolation filters for video decompression processor
KR100226684B1 (ko) * 1996-03-22 1999-10-15 전주범 반화소 움직임 추정장치
US6288745B1 (en) * 1997-04-24 2001-09-11 Mitsubishi Denki Kabushiki Kaisha Scanner line interpolation device
JPH11122626A (ja) * 1997-10-17 1999-04-30 Nikon Corp 画像処理方法及び装置並びに画像処理プログラムを記録した記録媒体
US6259734B1 (en) * 1998-04-03 2001-07-10 Matsushita Electric Industrial Co., Ltd. Image processing method, image processing apparatus, and data storage media
US6628714B1 (en) * 1998-12-18 2003-09-30 Zenith Electronics Corporation Down converting MPEG encoded high definition sequences to lower resolution with reduced memory in decoder loop
US6658157B1 (en) * 1999-06-29 2003-12-02 Sony Corporation Method and apparatus for converting image information
JP4361991B2 (ja) 1999-08-20 2009-11-11 メディア・テック・ユーエスエイ・インコーポレーテッド 画像処理装置
US6510178B1 (en) * 1999-12-15 2003-01-21 Zenith Electronics Corporation Compensating for drift in the down conversion of high definition sequences to lower resolution sequences
US6950469B2 (en) * 2001-09-17 2005-09-27 Nokia Corporation Method for sub-pixel value interpolation
US7245326B2 (en) 2001-11-19 2007-07-17 Matsushita Electric Industrial Co. Ltd. Method of edge based interpolation
US8175159B2 (en) * 2002-01-24 2012-05-08 Hitachi, Ltd. Moving picture signal coding method, decoding method, coding apparatus, and decoding apparatus
JP4120301B2 (ja) * 2002-04-25 2008-07-16 ソニー株式会社 画像処理装置およびその方法
JP4318019B2 (ja) * 2002-05-28 2009-08-19 ソニー株式会社 画像処理装置および方法、記録媒体、並びにプログラム
US9137541B2 (en) * 2003-05-23 2015-09-15 Broadcom Corporation Video data cache

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102265619A (zh) * 2008-09-23 2011-11-30 高通股份有限公司 切换式内插滤波器中的偏移计算
US8750378B2 (en) 2008-09-23 2014-06-10 Qualcomm Incorporated Offset calculation in switched interpolation filters
CN102265619B (zh) * 2008-09-23 2014-06-18 高通股份有限公司 切换式内插滤波器中的偏移计算方法

Also Published As

Publication number Publication date
NO320114B1 (no) 2005-10-24
NO20035444L (no) 2005-06-06
NO20035444D0 (no) 2003-12-05
EP1695560B1 (en) 2019-06-19
JP2007515115A (ja) 2007-06-07
JP4761390B2 (ja) 2011-08-31
WO2005055615A1 (en) 2005-06-16
US7616689B2 (en) 2009-11-10
US20050123040A1 (en) 2005-06-09
EP1695560A1 (en) 2006-08-30
CN100589580C (zh) 2010-02-10

Similar Documents

Publication Publication Date Title
US11240498B2 (en) Independently coding frame areas
CN1264355C (zh) 视频编码方法
CN1890983A (zh) 改进的内插像素值的计算
US7499495B2 (en) Extended range motion vectors
CN1170436C (zh) 压缩图像比特流的变码方法
US20050063471A1 (en) Flexible range reduction
CN1679330A (zh) 视频错误隐藏的系统和方法
JP2011505781A (ja) 高解像度デジタル静止画像をビデオと並行して符号化するためのavc規格の拡張
WO2011086836A1 (ja) 符号化装置、復号装置、および、データ構造
KR20040077774A (ko) 화상 정보 부호화 장치 및 방법, 및 화상 정보 복호 장치및 방법
CN1875628A (zh) 校正内插像素值的方法
CN1767644A (zh) 用于视频编码的非整数像素共享
CN1806444A (zh) 自适应加权运动估计方法和采用该方法的帧频转换设备
CN1969556A (zh) 色度去块的方法
US7502415B2 (en) Range reduction
CN1813480A (zh) 变换域中零的早期检测
EP2532163B1 (en) Improved method and apparatus for sub-pixel interpolation
CN1232933C (zh) 图像序列的格式转换方法
CN1157954C (zh) 编码视频数据的转换
US20060256860A1 (en) Transcoding with look-ahead
US20050141608A1 (en) Pipeline-type operation method for a video processing apparatus and bit rate control method using the same
WO2008097104A1 (en) Method for pixel prediction with low complexity
JP4734301B2 (ja) 量子化誤差低減装置
CN1929617A (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
ASS Succession or assignment of patent right

Owner name: CISCO SYSTEMS INTERNATIONAL BV

Free format text: FORMER OWNER: TANDBERG TELECOM AS

Effective date: 20120406

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20120406

Address after: Swiss basil

Patentee after: Cisco Systems International Sarl

Address before: Norway lsaac

Patentee before: Tandberg Telecom AS

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: 20100210

Termination date: 20201001