CN102214366A - 三维点采样数据高性能渲染方法 - Google Patents

三维点采样数据高性能渲染方法 Download PDF

Info

Publication number
CN102214366A
CN102214366A CN 201110203177 CN201110203177A CN102214366A CN 102214366 A CN102214366 A CN 102214366A CN 201110203177 CN201110203177 CN 201110203177 CN 201110203177 A CN201110203177 A CN 201110203177A CN 102214366 A CN102214366 A CN 102214366A
Authority
CN
China
Prior art keywords
point
sampled
dimensional
sampled point
sampling
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 201110203177
Other languages
English (en)
Other versions
CN102214366B (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.)
Zhejiang Wanli College
Original Assignee
Zhejiang Wanli College
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 Zhejiang Wanli College filed Critical Zhejiang Wanli College
Priority to CN 201110203177 priority Critical patent/CN102214366B/zh
Publication of CN102214366A publication Critical patent/CN102214366A/zh
Application granted granted Critical
Publication of CN102214366B publication Critical patent/CN102214366B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Generation (AREA)

Abstract

本发明提供一种三维点采样数据高性能渲染方法,包括下述步骤:通过构建3D采样点最近邻域集的局部隐式曲面并结合微分几何理论,估算其法向信息;通过协方差分析邻域点集,确定3D采样点在其切平面上的椭圆表示;利用具有良好反走样效果的椭圆加权滤波对3D采样点数据进行绘制,以获得高质量的数据渲染效果;依据采样点在当前视角下深度值得到的排序,设计三维点采样数据的单遍绘制算法,以加快渲染速度。采用本方法可对3D点采样数据进行高性能即快速且高质量的渲染。

Description

三维点采样数据高性能渲染方法
技术领域
本发明涉及计算机图形处理领域,特别涉及一种三维点采样数据高性能渲染方法。
背景技术
3D扫描获取技术的快速发展,使3D点采样数据模型已成为继一维的声音数据、二维的图像数据与视频数据之后的一种新兴数字媒体。在逆向工程、工业产品创新设计、数字娱乐、影视动画、电子商务、通信、物理模拟、文物保护与修复等领域中,3D点采样数据模型有着广泛的应用,且产生了越来越深远的影响。近年来,已涌现出多种3D点采样数据的建模和造型方法,这使得高性能、快速高质量地渲染绘制3D点采样数据成为一项非常重要的工作。
自2000年Rusinkiewicz和Levoy开发了QSplat点绘制系统以来,已出现了多种基于点的绘制方法和技术。QSplat利用层次包围球数据结构,实现了3D点采样数据的快速显示,然而其不足是产生严重的图形走样。Zwicker等将信号处理的思想应用到点的绘制中,提出了椭圆加权Elliptical Weighted Average,EWA算法,将屏幕空间的重采样和滤波巧妙地结合在一起,很好的解决了点绘制的走样问题,得到了高质量的绘制结果;但该算法是纯软件实现的。Ren等推导出物体空间的EWA滤波算子,并基于GPU实现该算法,因此提高了绘制速度。Dachsbacher等将点模型的层次结构序列化为线性数组,基于此将绘制算法由层次树的遍历转化为线性数组的遍历,从而提高了绘制效率。
基于硬件的EWA算法尽管提高了绘制速度,然而也存在相应的不足。其一,通常需要进行多遍绘制,在第一遍绘制的时候,需要将物体向远离视角的方向移动一个微小的偏移量后再对深度缓存进行绘制,然而偏移量大小的确定不是件容易的事,偏移量太大或者太小都会影响到绘制效果。其二,多遍绘制算法只绘制位于视角最前端的点,因此无法逼真的绘制半透明体。
发明内容
本发明克服了上述现有技术中的不足,提供了一种三维点采样数据高性能渲染方法,通过设计单遍绘制方案提高渲染速度,通过EWA增强渲染的质量。
本发明的技术方案是这样实现的:
一种三维点采样数据高性能渲染方法,所述渲染方法包括下述步骤:(1)构建采样点最近邻域集的局部隐式曲面,基于此估算法向;(2)协方差分析邻域点集,确定采样点在其切平面上的椭圆表示;(3)利用椭圆加权滤波,绘制3D采样点;(4)按照采样点在当前视角下深度值进行的动态排序,设计三维点采样数据从远到近的单遍绘制方法。
为实现发明目的,在步骤(1)中,采样点的法向是通过计算隐式曲面的梯度得到,具体是:
(a)利用kD树快速搜索采样点的k最近邻域Nk(pi),依据三维采样点采样数据的规模n取k∈[9,30];
(b)基于径向基函数RBF,RBF即radial basis functions,并采用法矢约束的方式加入约束点,在采样点集Nk(pi)上重构采样点pi的局部隐式曲面。
给定R3空间的l个散乱点{p1,...,pl},插值它们的RBF隐式曲面一般形式为:
f ( r ) = Σ i = 1 l w i φ ( r - p i ) + P ( r ) = 0 式①
由l个散乱点所给定的方程约束条件为
f ( p i ) = Σ j = 1 l w j φ ( p i - p j ) + P ( p i ) = h i , i = 1,2 , . . . , l ,
以及使能量方程取得最小化所需满足的正交条件为
Σ j = 1 l w j = Σ j = 1 l w j p j x = Σ j = 1 l w j p j y = Σ j = 1 l w j p j z = 0 ,
令φij=φ(pi-pj),建立下面的线性方程组
φ 11 φ 12 . . . φ 1 l 1 p 1 x p 1 y p 1 z φ 21 φ 22 . . . φ 2 l 1 p 2 x p 2 y p 2 z M M O M M M M M φ l 1 φ l 2 . . . φ ll 1 p l x p l y p l z 1 1 . . . 1 0 0 . . . 0 p 1 x p 2 x . . . p l x 0 0 . . . 0 p 1 y p 2 y . . . p l y 0 0 . . . 0 p 1 z p 2 z . . . p l z 0 0 . . . 0 w 1 w 2 M w l a 0 a 1 a 2 a 3 = h 1 h 2 M h l 0 0 0 0 式②
为了使式②有唯一一组解(w1,w2,...,wn,a0,a1,a2,a3),采用法矢约束的方式加入约束点;将求得的结果代入式①,便可得到隐式曲面方程:
f ( x , y , z ) = Σ j = 1 n w j ( ( x - p j x ) 2 + ( y - p j y ) 2 + ( z - p j z ) 2 ) 3 + a 0 + a 1 x + a 2 y + a 3 z 式③
(c)计算隐式曲面在采样点pi处的梯度得其法向并采用最小生成树传播法对ni进行全局一致化处理。
为实现发明目的,所述的三维点采样数据高性能渲染方法,在步骤(2)中,定义采样点pi的k最近邻域Nk(pi)的协方差矩阵并进行分析,其3个特征值为非负值,由此确定Nk(pi)的最佳包围体椭圆球的3个轴长;其3个特征向量组成一正交基,由此确定采样点pi在其切平面上的椭圆表示。更为详细的是:
点集Nk(pi)的协方差矩阵为 C = p i 1 - p ‾ . . . p i k - p ‾ T · p i 1 - p ‾ . . . p i k - p ‾ , 其中
Figure BDA0000077048550000036
为Nk(pi)的质心,
Figure BDA0000077048550000037
矩阵C是对称半正定的,其三个特征值λi(i=0,1,2)为非负实值(设0≤λ0≤λ1≤λ2),所对应的三个特征向量vi(i=0,1,2)组成一个正交基;平面
Figure BDA0000077048550000041
使得
Figure BDA0000077048550000042
周围的点到此平面的距离和为最小,此平面可视为采样点pi的最近邻域Nk(pi)的切平面,v0则可作为局部曲面在采样点pi的近似法向ni。于是在局部标架
Figure BDA0000077048550000043
中,Nk(pi)的最佳包围体椭圆球为:
Figure BDA0000077048550000044
其中f为缩放因子;因此,在采样点pi切平面上的局部二维标架(pi,v2,v1)中,采样点pi的椭圆表示为:
Figure BDA0000077048550000045
为实现发明目的,在步骤(3)中,为了高质量地渲染三维点采样数据,利用了椭圆加权滤波的方法,具体是:在采样点pi切平面的局部标架中,定义pi的重构核,以便连续地重构点采样模型;由此,二维屏幕空间的椭圆加权滤波重采样函数便是景物空间点数据重构滤波与二维屏幕低通滤波的卷积。
为实现发明目的,所述的三维点采样数据高性能渲染方法,在步骤(4)中,为了快速地渲染三维点采样数据,设计了单遍绘制方法,具体是;首先,根据采样点的x、y和z坐标由小到大排序得到3组顺序表;然后,根据当前视点判定绘制时所用的顺序表,以避免每次绘制的排序问题,加快了绘制速度。对于不可见采样点的剔除,利用八叉树来加速该过程:遍历八叉树,当结点的包围球投影在屏幕之外或结点的法锥面完全背向于视点,则剔除该节点所包含的采样点。同时,该方法遵循从后往前的绘制原则,故直接用开放性图形库OpenGL绘制可得到透明的数据绘制效果。
采用了上述技术方案的本发明的原理及有益效果是:
本发明利用具有良好反走样效果的椭圆加权滤波对3D采样点数据进行绘制,从而能够获得高质量的数据渲染效果;本发明依据采样点在当前视角下深度值得到的排序,设计出三维点采样数据的单遍绘制算法,从而能够更加快速地进行渲染。因此,采用本发明能够对3D点采样数据进行高性能即快速且高质量的渲染。
附图说明
图1为本发明实施例中人手骨骼的三维点采样数据图;
图2为与图1对应的渲染绘制效果图;
图3为本发明实施例中雕花球的三维点采样数据图;
图4为与图3对应的渲染绘制效果图;
图5为本发明实施例中的摇臂模型的渲染效果图;
图6为本发明实施例中机械零件模型的渲染效果图;
图7为本发明实施例中马克斯普朗克头像模型的渲染效果图;
图8为本发明实施例中伊丽莎白头像模型的渲染效果图;
图9为本发明实施例中龙的渲染效果图;
图10为本发明实施例中兔子的渲染效果图;
图11是采用QSplat方法渲染的花结的效果图;
图12为图11的A部放大图;
图13是采用本发明方法渲染的花结的效果图;
图14为图13的A部放大图。
具体实施方式
本发明的具体实施方式如下:
实施例:参见图1~图4,一种三维点采样数据高性能渲染方法,所述方法包括:
(1)构建采样点最近邻域集的局部隐式曲面,基于此估算法向。具体如下:(a)利用kD树快速搜索采样点pi的k最近邻域Nk(pi);(b)基于径向基函数RBF,并采用法矢约束的方式向式② φ 11 φ 12 . . . φ 1 l 1 p 1 x p 1 y p 1 z φ 21 φ 22 . . . φ 2 l 1 p 2 x p 2 y p 2 z M M O M M M M M φ l 1 φ l 2 . . . φ ll 1 p l x p l y p l z 1 1 . . . 1 0 0 . . . 0 p 1 x p 2 x . . . p l x 0 0 . . . 0 p 1 y p 2 y . . . p l y 0 0 . . . 0 p 1 z p 2 z . . . p l z 0 0 . . . 0 w 1 w 2 M w l a 0 a 1 a 2 a 3 = h 1 h 2 M h l 0 0 0 0 加入约束点,在采样点集Nk(pi)上重构采样点pi的局部隐式曲面f(x,y,z),以上由式③ f ( x , y , z ) = Σ j = 1 n w j ( ( x - p j x ) 2 + ( y - p j y ) 2 + ( z - p j z ) 2 ) 3 + a 0 + a 1 x + a 2 y + a 3 z 得到;(c)计算采样点pi的隐式曲面在采样点pi处的梯度得其法向
Figure BDA0000077048550000063
并采用最小生成树传播法对所有的ni进行全局一致化处理。
(2)协方差分析邻域点集,确定采样点在其切平面上的椭圆表示。具体如下:点集Nk(pi)的协方差矩阵为 C = p i 1 - p ‾ . . . p i k - p ‾ T · p i 1 - p ‾ . . . p i k - p ‾ , 其中
Figure BDA0000077048550000065
为Nk(pi)的质心,
Figure BDA0000077048550000066
矩阵C是对称半正定的,其三个特征值λi(i=0,1,2)为非负实值(设0≤λ0≤λ1≤λ2),所对应的三个特征向量vi(i=0,1,2)组成一个正交基;平面
Figure BDA0000077048550000067
可视为Nk(pi)的切平面,v0可作为局部曲面在pi点的近似法向ni。在局部标架
Figure BDA0000077048550000068
中,Nk(pi)的最佳包围体椭圆球为:
Figure BDA0000077048550000069
因此,在pi切平面上的局部二维标架(pi,v2,v1)中,pi的椭圆表示为: x 2 λ 2 2 + y 2 λ 1 2 = f 2 .
(3)利用椭圆加权滤波,绘制3D采样点。具体如下:在采样点pi切平面的局部标架中,定义采样点pi的重构核,以连续地重构点采样模型;于是二维屏幕空间的椭圆加权滤波重采样函数便是景物空间点数据重构滤波与二维屏幕低通滤波的卷积。因此,实现三维点采样数据的高质量渲染绘制。
(4)按照采样点在当前视角下深度值进行的动态排序,三维点采样数据从远到近进行单遍绘制。具体如下:首先,根据采样点的x、y和z坐标由小到大排序得到3组顺序表;然后,根据当前视点判定绘制时所用的顺序表,以避免每次绘制的排序问题,从而加快了绘制速度。对于不可见采样点的剔除,利用八叉树来加速该过程:遍历八叉树,当结点的包围球投影在屏幕之外或结点的法锥面完全背向于视点,则剔除该节点所包含的采样点。
本发明中,“三维点采样数据”是由n个采样点构成,而各个采样点是独立的;于是对每个采样点(比如一采样点pi)进行的运算是在其最近邻域(可理解为邻居)上进行的,那么最近邻域通常取k个(可理解为取k个最近的邻居)。
图2、图4~图10、图13、图14显示了采用本发明方法得到的多个模型的渲染效果图。可以看出,图形质量较高,没有走样。
表1给出了三维点采样数据模型的渲染绘制实验统计数据,从中表明该技术的快速高效性,其中微机的运行环境为P4 2.8GHz、1GB内存、Radeon X300显卡。实验时,绘制窗口分辨率为512*512,由表1的帧率栏(前行是本发明的,后行是Ren的EWA法的)可以看出,本发明的绘制明显快于后者的,达到了4MPS(每秒百万个点元)的快速绘制速度。
表1
Figure BDA0000077048550000071
Figure BDA0000077048550000081
参照图11~图14给出了QSplat渲染效果和本发明渲染效果的对比,可以看出本发明的绘制效果消除了QSplat方法的图形走样现象,生成了高质量的绘制效果。

Claims (5)

1.一种三维点采样数据高性能渲染方法,其特征在于,所述渲染方法包括下述步骤:(1)构建采样点最近邻域集的局部隐式曲面,基于此估算法向;(2)协方差分析邻域点集,确定采样点在其切平面上的椭圆表示;(3)利用椭圆加权滤波,绘制3D采样点;(4)按照采样点在当前视角下深度值进行的动态排序,设计三维点采样数据从远到近的单遍绘制方法。
2.根据权利要求1所述的三维点采样数据高性能渲染方法,其特征在于:在步骤(1)中,采样点的法向是通过计算隐式曲面的梯度得到,具体是:(a)利用kD树快速搜索采样点的k最近邻域Nk(pi),依据三维点采样数据的规模n取k∈[9,30];(b)基于径向基函数RBF,并采用法矢约束的方式加入约束点,在采样点集Nk(pi)上重构采样点pi的局部隐式曲面;(c)计算采样点pi的隐式曲面在采样点pi处的梯度得其法向ni,并采用最小生成树传播法对ni进行全局一致化处理。
3.根据权利要求1所述的三维点采样数据高性能渲染方法,其特征在于:在步骤(2)中,定义采样点的k最近邻域Nk(pi)的协方差矩阵并进行分析,其3个特征值为非负值,由此确定采样点的k最近邻域Nk(pi)的最佳包围体椭圆球的3个轴长;其3个特征向量组成一正交基,由此确定采样点pi在其切平面上的椭圆表示。
4.根据权利要求1所述的三维点采样数据高性能渲染方法,其特征在于:在步骤(3)中,具体是:在采样点pi切平面的局部标架中,定义pi的重构核,以便连续地重构点采样模型;二维屏幕空间的椭圆加权滤波重采样函数便是景物空间点数据重构滤波与二维屏幕低通滤波的卷积。
5.根据权利要求1所述的三维点采样数据高性能渲染方法,其特征在于:在步骤(4)中,具体是;首先,根据采样点的x、y和z坐标由小到大排序得到3组顺序表;然后,根据当前视点判定绘制时所用的顺序表,以避免每次绘制的排序问题,加快了绘制速度;对于不可见采样点的剔除,利用八叉树来加速该过程:遍历八叉树,当结点的包围球投影在屏幕之外或结点的法锥面完全背向于视点,则剔除该节点所包含的采样点;同时,该方法遵循从后往前的绘制原则,故直接用开放性图形库OpenGL绘制可得到透明的数据绘制效果。
CN 201110203177 2011-07-20 2011-07-20 三维点采样数据高性能渲染方法 Expired - Fee Related CN102214366B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110203177 CN102214366B (zh) 2011-07-20 2011-07-20 三维点采样数据高性能渲染方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110203177 CN102214366B (zh) 2011-07-20 2011-07-20 三维点采样数据高性能渲染方法

Publications (2)

Publication Number Publication Date
CN102214366A true CN102214366A (zh) 2011-10-12
CN102214366B CN102214366B (zh) 2013-01-16

Family

ID=44745661

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110203177 Expired - Fee Related CN102214366B (zh) 2011-07-20 2011-07-20 三维点采样数据高性能渲染方法

Country Status (1)

Country Link
CN (1) CN102214366B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103177474A (zh) * 2013-03-14 2013-06-26 腾讯科技(深圳)有限公司 三维模型的邻域点坐标确定方法及装置、构建方法及装置
CN103279972A (zh) * 2013-05-08 2013-09-04 中国科学院软件研究所 一种基于kd树的并行多维自适应采样与重构方法
CN104751507A (zh) * 2013-12-31 2015-07-01 北界创想(北京)软件有限公司 图形内容渲染方法和装置
CN106780746A (zh) * 2016-11-28 2017-05-31 中国科学院自动化研究所 基于简单互斥操作的蓝噪声采样方法
CN108805261A (zh) * 2017-04-28 2018-11-13 微软技术许可有限责任公司 基于八叉树的卷积神经网络

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1430184A (zh) * 2001-12-29 2003-07-16 田捷 利用半边数据结构实现三维网格模型的简化方法
US20040012603A1 (en) * 2002-07-19 2004-01-22 Hanspeter Pfister Object space EWA splatting of point-based 3D models
US20060077204A1 (en) * 2004-10-08 2006-04-13 Hanspeter Pfister Sample rate adaptive filtering for volume rendering
CN101763649A (zh) * 2009-12-30 2010-06-30 北京航空航天大学 一种增强模型轮廓的表面点绘制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1430184A (zh) * 2001-12-29 2003-07-16 田捷 利用半边数据结构实现三维网格模型的简化方法
US20040012603A1 (en) * 2002-07-19 2004-01-22 Hanspeter Pfister Object space EWA splatting of point-based 3D models
US20060077204A1 (en) * 2004-10-08 2006-04-13 Hanspeter Pfister Sample rate adaptive filtering for volume rendering
CN101763649A (zh) * 2009-12-30 2010-06-30 北京航空航天大学 一种增强模型轮廓的表面点绘制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《EUROGRAPHICS》 20020930 Liu Ren et al. Object Space EWA Surface Splatting: A Hardware Accelated Approach to High Quality Point Rendering 全文 1-5 第21卷, 第3期 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103177474A (zh) * 2013-03-14 2013-06-26 腾讯科技(深圳)有限公司 三维模型的邻域点坐标确定方法及装置、构建方法及装置
CN103177474B (zh) * 2013-03-14 2016-04-20 腾讯科技(深圳)有限公司 三维模型的邻域点坐标确定方法及装置、构建方法及装置
CN103279972A (zh) * 2013-05-08 2013-09-04 中国科学院软件研究所 一种基于kd树的并行多维自适应采样与重构方法
CN104751507A (zh) * 2013-12-31 2015-07-01 北界创想(北京)软件有限公司 图形内容渲染方法和装置
CN104751507B (zh) * 2013-12-31 2018-01-19 北界创想(北京)软件有限公司 图形内容渲染方法和装置
CN106780746A (zh) * 2016-11-28 2017-05-31 中国科学院自动化研究所 基于简单互斥操作的蓝噪声采样方法
CN106780746B (zh) * 2016-11-28 2020-03-27 中国科学院自动化研究所 基于简单互斥操作的蓝噪声采样方法
CN108805261A (zh) * 2017-04-28 2018-11-13 微软技术许可有限责任公司 基于八叉树的卷积神经网络
CN108805261B (zh) * 2017-04-28 2021-11-12 微软技术许可有限责任公司 基于八叉树的卷积神经网络
US11704537B2 (en) 2017-04-28 2023-07-18 Microsoft Technology Licensing, Llc Octree-based convolutional neural network

Also Published As

Publication number Publication date
CN102214366B (zh) 2013-01-16

Similar Documents

Publication Publication Date Title
Lin et al. Magic3d: High-resolution text-to-3d content creation
CN102157008B (zh) 一种大规模虚拟人群实时绘制方法
Zeng et al. Octree-based fusion for realtime 3D reconstruction
Jalba et al. Surface and curve skeletonization of large 3D models on the GPU
CN103530907B (zh) 基于图像的复杂三维模型绘制方法
CN109147048A (zh) 一种利用单张彩色图的三维网格重建方法
CN103021017A (zh) 基于gpu加速的三维场景重建方法
CN110490917A (zh) 三维重建方法及装置
CN102214366B (zh) 三维点采样数据高性能渲染方法
CN1818977A (zh) 由一幅正面图像实现快速人脸模型重建的方法
CN103679808A (zh) 利用二维核磁图像重建头部三维模型的方法及系统
Rabby et al. BeyondPixels: A comprehensive review of the evolution of neural radiance fields
Zhu et al. Nonlocal low-rank point cloud denoising for 3-D measurement surfaces
Yang et al. A completely parallel surface reconstruction method for particle-based fluids
Vyatkin Method of binary search for image elements of functionally defined objects using graphics processing units
Meyer et al. PEGASUS: Physically Enhanced Gaussian Splatting Simulation System for 6DOF Object Pose Dataset Generation
Wu et al. DeferredGS: Decoupled and Editable Gaussian Splatting with Deferred Shading
CN110738719A (zh) 一种基于视距分层优化的Web3D模型渲染方法
Li et al. 3d shape reconstruction of furniture object from a single real indoor image
CN112002019A (zh) 一种基于mr混合现实的模拟人物阴影的方法
Hao et al. Image completion with perspective constraint based on a single image
Legrand et al. Morton integrals for high speed geometry simplification
Gao et al. Mani-GS: Gaussian Splatting Manipulation with Triangular Mesh
Jeong et al. Detail control in line drawings of 3D meshes
Yang et al. Spatial Display Model of Oil Painting Art Based on Digital Vision Design.

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

Granted publication date: 20130116

Termination date: 20130720