CN104537692B - 基于时空上下文信息辅助的关键点稳定跟踪方法 - Google Patents

基于时空上下文信息辅助的关键点稳定跟踪方法 Download PDF

Info

Publication number
CN104537692B
CN104537692B CN201410841442.9A CN201410841442A CN104537692B CN 104537692 B CN104537692 B CN 104537692B CN 201410841442 A CN201410841442 A CN 201410841442A CN 104537692 B CN104537692 B CN 104537692B
Authority
CN
China
Prior art keywords
target
context
key point
target location
model
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
CN201410841442.9A
Other languages
English (en)
Other versions
CN104537692A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201410841442.9A priority Critical patent/CN104537692B/zh
Publication of CN104537692A publication Critical patent/CN104537692A/zh
Application granted granted Critical
Publication of CN104537692B publication Critical patent/CN104537692B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)

Abstract

一种基于时空上下文信息辅助的关键点稳定跟踪方法,本发明包括以下步骤:首先通过人工方式选择或自动检测手段提取目标上的关键点,然后构建关键点与目标周围特征点的时空上下文模型,并通过在线学习方式对该模型进行更新,待获取新一帧图像,求取目标置信度响应,响应度最大的位置即为待跟踪的关键点在新一帧图像中的位置;当目标物体由于遮挡或表观变化严重而发生信息丢失时,利用上下文信息来预测目标位置,同时可以抑制由于错误表观信息导致的漂移。用本发明方法跟踪目标时,精度高,鲁棒性、适应性强,能实现遮挡、光照改变、快速运动等各种复杂情形下关键点的稳定跟踪。

Description

基于时空上下文信息辅助的关键点稳定跟踪方法
技术领域
本发明涉及一种计算机视觉领域的视觉目标跟踪方法。
背景技术
关键点是一类典型的目标特征点,对关键点的稳定跟踪是实现目标精确捕获的前提。在许多应用领域,都需要对目标进行精确定位,但由于噪声、遮挡、杂乱背景、光照变化、姿态变化及视角变化等干扰因素的影响,关键点位置也随着目标运动而发生改变,它可能被遮蔽成目标的内点,也可能外界遮档,因此即使跟踪到了目标,准确地跟踪关键点的位置也非常困难。目前对关键点的跟踪方法主要由两类,一是基于光流的跟踪,将特征点跟踪看做成是光流计算问题,通常在参考帧图像中选择一组特征点,并且假设特征点纹理在帧间保持不变,然后通过局部匹配实现对特征点的跟踪,然而光流场计算量大、实时性不强,还有孔径问题;另一类是基于特征的跟踪,首先提取目标稳定的特征点,然后对此特征点进行稳定跟踪。然而,当图像中目标发生遮挡时,可能出现特征隐藏或虚假特征及隐藏特征再现等情形,影响关键点跟踪的稳定性。
此外,图像序列中的运动目标物体一般都存在于一个上下文环境中,目标与上下文环境中的其他物体之间或自身不同时刻状态之间存在着紧密的联系。使用目标与上下文之间的关系,可以预测目标在下一时刻可能出现的位置。基于上下文的运动目标跟踪是指在跟踪过程中使用了上下文信息来确定目标物体位置的跟踪方法。常用的上下文信息包括与目标物体运动相关的近邻物体、物体自身的局部子区域、物体的运动轨迹以及轨迹之间的联系、物体运动的速度和加速度以及一些合理的推测等。当目标物体处于一个复杂的环境中,特别是当发生遮挡甚至完全遮蔽时,可利用目标周围上下文信息实现对目标的精确定位。使用上下文信息进行跟踪不仅可以利用目标物体自身的表观信息,还可以利用自身以外的多种信息,从而扩大了跟踪时可利用的信息来源。
发明内容
本发明目的是解决遮挡、光照改变、快速运动等各种复杂情形下关键点的稳定跟踪问题。
本发明利用时空上下文辅助信息实现对关键点的稳定跟踪,首先通过人工选择方式或自动检测手段提取目标上的关键点,然后构建关键点与目标周围特征点的时空上下文模型,并通过在线学习方式对该模型进行更新。待获取新一帧图像后,求取目标位置置信度响应,响应度最大的位置即为待跟踪的关键点在新一帧图像中的位置。当待跟踪目标由于遮挡或表观严重变化而丢失时,可利用上下文信息来预测目标位置。此外,还可用上下文辅助以抑制由于错误表观引起的跟踪漂移。本发明的具体设计方案包括步骤:
A建立目标位置估计模型
将关键点跟踪看成是求解目标位置估计似然函数的最大值问题,即求解目标位置似然的置信度:
z t = arg max z ∈ S ( z t ) l t ( z ) - - - ( 1 )
其中,lt(z)=P(zt|st),lt(z)为目标位置的似然函数,zt∈R2是当前帧中目标位置,st代表目标当前状态。
引入上下文信息对目标位置进行估计,令FS={f(z)=(I(y),y)|y∈S(zt)}表示与目标相关的上下文信息集合。I(y)是位于y的像素点的灰度值。S(zt)是当前帧目标位置zt周围上下文区域。将目标位置估计似然转换为目标周围上下文信息对目标位置的投影映射:
l ( z ) = P ( z | c ) = Σ f ( z ) ∈ F S P ( z , f ( y ) | c ) = Σ f ( z ) ∈ F S P ( z | f ( y ) , c ) P ( f ( y ) | c ) - - - ( 2 )
使置信度值最大的z*就是目标的位置。从公式中可以看到,似然函数可以分解为两个概率部分:一个是目标与周围上下文信息的空间关系的条件概率P(z|f(y),c),需要通过在线学习得出;一个是局部上下文位置点z的上下文先验概率P(f(y)|c)。
B建立目标与上下文的空间关系模型
空间上下文模型描述的是条件概率函数,将其建模为:
P(z|f(y),c)=h(z-y)(3)
其中,h(z-y)是一个关于目标z与上下文位置y的距离方向函数,它表征目标与上下文的空间关系,该函数径向不对称,有助于分辨二义性,防止跟踪漂移到相似目标上。
C建立上下文先验概率模型
根据生物视觉注意机理,处于某个视觉注意区域内的目标,离其越近的点关注得越多,反之关注得越小。该注意机理反映在物理描述上就是离目标越近的点被赋予的权重越大。据此,将上下文先验概率P(f(y)|c)建模为:
P(f(y)|c)=I(y)ωc(y-z*)(4)
其中I(y)是点y的灰度,ωc(·)为一个加权函数,y离z*越近,权值越大。
D时空上下文模型的在线学习更新
将式(3)、式(4)代入式(2),可得目标位置的似然估计值为
l ( z ) = Σ y ∈ F S h ( z - y ) I ( y ) ω c ( y - z * ) = h ( z ) ⊗ ( I ( z ) ω c ( z - z * ) ) - - - ( 5 )
其中,表示卷积,对式(5)两边同时进行快速傅里叶变换(FFT),可得:
其中,F(·)表示FFT,FFT表示快速傅里叶变换,表示点乘。空间上下文模型为:
h ( z ) = F - 1 ( F ( l ( z ) ) F ( I ( z ) ω c ( z - z * ) ) ) - - - ( 7 )
F-1(·)表示逆FFT,由(6)式及(7)式可知,只需要利用两个FFT和一个逆FFT,就能快速求取空间上下文条件概率。
将时空上下文模型更新看成是一个无限长滤波过程,即t+1时刻新的空间上下文与t时刻空间上下文满足以下关系:
Ht+1(z)=(1-ρ)Ht(z)+ρht(z)(8)
在获取t+1时刻新的时空上下文关系后,就可以两个正傅里叶变换和一个逆傅里叶变换来加速计算新一帧中目标位置的置信度响应:
目标在下一帧的位置对应置信度响应最大的位置
z t + 1 = arg max z ∈ S ( z t ) l t + 1 ( z ) - - - ( 10 )
用本发明方法跟踪目标时,精度高、鲁棒性、适应性强,能实现遮挡、光照改变、快速运动等各种复杂情形下关键点的稳定跟踪。
附图说明
图1是人工选定关键点及其周围区域置信度分布图;
图2为目标置信度投影分布图;
图3为具体跟踪效果图;
图4为关键点的中心定位误差曲线。
具体实施方式
本发明可以跟踪包括旋转、尺度缩放、光照变化等复杂场景中人工任意选定的关键点,并具有稳定、鲁棒的跟踪性能。下面结合附图,对本发明的具体原理做进一步的详细阐述。
实施例1短时间内目标与周围上下文的关系不会发生大的变化,因此只要保持这种关系并对其进行实时更新,当遮挡出现时,就能利用未遮挡条件下的上下文信息对目标位置进行预测。即使目标不可见或消失,只要这种关系存在且稳定,就能有效利用这些上下文信息预测出目标位置。
A建立目标位置估计模型
将关键点跟踪看成是求解目标位置估计似然函数的最大值问题,即求解目标位置似然的置信度:
z t = arg max z ∈ S ( z t ) l t ( z ) - - - ( 1 )
lt(z)=P(zt|st),lt(z)为目标位置的似然函数,zt∈R2是当前帧中目标位置,st代表目标当前状态。
引入上下文信息对目标位置进行估计,令FS={f(z)=(I(y),y)|y∈S(zt)}表示与目标相关的上下文信息集合。I(y)是位于y的像素点的灰度值。S(zt)是当前帧目标位置zt周围上下文区域。将目标位置估计似然转换为目标周围上下文信息对目标位置的投影映射:
l ( z ) = P ( z | c ) = Σ f ( z ) ∈ F S P ( z , f ( y ) | c ) = Σ f ( z ) ∈ F S P ( z | f ( y ) , c ) P ( f ( y ) | c ) - - - ( 2 )
使置信度值最大的z*就是目标的位置。从公式上可以看到,似然函数可以分解为两个概率部分:一个是目标与周围上下文信息的空间关系的条件概率P(z|f(y),c),需要通过在线学习得出。一个是局部上下文位置点z的上下文先验概率P(f(y)|c)。
B建立目标与上下文的空间关系模型
上下文的空间模型描述的是条件概率函数,将其建模为:
P(z|f(y),c)=h(z-y)(3)
h(z-f(y))是一个关于目标z与上下文位置y的距离方向函数,它表征目标与上下文的空间关系,该函数径向不对称,有助于分辨二义性,防止跟踪漂移到相似目标上。
C建立上下文先验概率模型
根据生物视觉注意机理,处于某个视觉注意区域内的目标,离其越近的点关注得越多,反之关注得越小。该注意机理反映在物理描述上就是离目标越近的点被赋予的权重越大。据此,将上下文先验概率P(f(y)|c)建模为:
P(f(y)|c)=I(y)ωc(y-z*)(4)
其中I(y)是点y的灰度,ωc(·)为一个加权函数,y离z*越近,权值越大。
D时空上下文模型的在线学习更新
将式(3)、式(4)代入式(2),得目标位置的似然估计值为
l ( z ) = Σ y ∈ F S h ( z - y ) I ( y ) ω c ( y - z * ) = h ( z ) ⊗ ( I ( z ) ω c ( z - z * ) ) - - - ( 5 )
表示卷积,对式(5)两边同时进行快速傅里叶变换,可得:
表示点乘。空间上下文模型为:
h ( z ) = F - 1 ( F ( l ( z ) ) F ( I ( z ) ω c ( z - z * ) ) ) - - - ( 7 )
只需要利用两个正傅里叶变换和一个逆傅里叶变换,就能快速求取空间上下文条件概率。
将时空上下文模型更新看成是一个IIR滤波过程,即t+1时刻新的空间上下文与t时刻空间上下文满足以下关系:
Ht+1(z)=(1-ρ)Ht(z)+ρht(z)(8)
在获取t+1时刻新的时空上下文关系后,就可以两个正傅里叶变换和一个逆傅里叶变换来加速计算新一帧中目标位置的置信度响应:
目标在下一帧的位置对应置信度响应最大的位置:
z t + 1 = arg max z ∈ S ( z t ) l t + 1 ( z ) - - - ( 10 )
图1给出了人工选定关键点及其周围区域置信度分布图,在目标位置置信图中,关键点所在位置具有最大的响应值。由于选定的关键点位于人脸,若采用传统基于特征点跟踪的方法对该关键点进行跟踪,由于人脸肤色相近,该关键点很容易漂移到与其颜色相近的特征点上,从而导致跟踪点在人脸上任意漂移。而本发明利用白框区域建立与关键点的时空上下文关系,同时通过在线学习方式更新该上下文关系,从而达到约束关键点运动与准确定位的目的。
图2为目标周围用于预测其位置的上下文信息对位置投影的分布平面图。为了直观区分不同区域上下文投影的不同,将观测视角旋转(-40,-7),由图2可见,待跟踪目标区域投影响应值最大,周围区域响应较小,故在利用上下文定位目标中心时,前者贡献大,后者贡献小,最终能准确捕获到新一帧中的目标位置。
图3给出了2组利用上下文信息辅助进行任意点跟踪的示例,由图3可见,虽然关键点周围存在与大量纹理相近的背景干扰区域,但本发明始终能准确追踪到初始选定待跟踪区域而不会漂移到相似背景区域上,避免了传统模板匹配、meanshift跟踪中出现的模板漂移问题,从而实现稳健、鲁棒的目标跟踪。
图4给出了利用本发明与CXT、l1等跟踪算法跟踪woman序列得到的中心定位误差曲线对比,该序列中的目标被严重遮挡。中心定位误差为目标真实位置与各算法跟踪结果的均方差,其中目标真实位置通过人工方式标定。由图可知,本发明得到的中心定位误差最小,从而说明本发明得到的跟踪位置与真实位置最为接近,跟踪精度最高,鲁棒性、适应性最强。

Claims (1)

1.一种基于时空上下文信息辅助的关键点稳定跟踪方法,其特征在于包括以下步骤:
首先通过人工选择方式或自动检测手段提取目标上的关键点,然后构建关键点与目标周围特征点的时空上下文模型,并通过在线学习方式对该模型进行更新,待获取新一帧图像,求取目标置信度响应,响应最大的位置即为待跟踪的关键点在新一帧图像中的位置;当跟踪目标由于遮挡或表观严重变化而丢失时,利用上下文信息来预测目标位置,此外,还可用上下文辅助以抑制由于错误表观引起的跟踪漂移,所述的目标置信度响应计算过程为:
A建立目标位置估计模型
将关键点跟踪看成是求解目标位置估计似然函数的最大值问题,即求解目标位置似然的置信度:
z t = a r g max z ∈ S ( z t ) l t ( z ) - - - ( 1 )
lt(z)=P(zt|st),lt(z)为目标位置的似然函数,zt∈R2是当前帧中目标位置,st代表目标当前状态;
引入上下文信息对目标位置进行估计,令FS={f(z)=(I(y),y)|y∈S(zt)}表示与目标相关的上下文信息集合;I(y)是位于y的像素点的灰度值;S(zt)是当前帧目标位置zt周围上下文区域,将目标位置估计似然转换为目标周围上下文信息对目标位置的投影映射:
l ( z ) = P ( z | c ) = Σ f ( z ) ∈ F S P ( z , f ( y ) | c ) = Σ f ( z ) ∈ F S P ( z | f ( y ) , c ) P ( f ( y ) | c ) - - - ( 2 )
使置信度值最大的z*就是目标的位置;从公式(2)知,似然函数可分解为两个概率部分:一个是目标与周围上下文信息的空间关系的条件概率P(z|f(y),c),通过在线学习得出;一个是局部上下文位置点z的上下文先验概率P(f(y)|c);
B建立目标与上下文的空间关系模型
空间上下文模型描述的是条件概率函数,将其建模为:
P(z|f(y),c)=h(z-y)(3)
h(z-y)是一个关于目标z与上下文位置y的距离方向函数,它表征目标与上下文的空间关系,该函数径向不对称,有助于分辨二义性,防止跟踪漂移到相似目标上;
C建立上下文先验概率模型
根据生物视觉注意机理,处于某个视觉注意区域内的目标,离其越近的点关注得越多,反之关注得越小;该注意机理反映在物理描述上就是离目标越近的点被赋予的权重越大;据此,将上下文先验概率P(f(y)|c)建模为:
P(f(y)|c)=I(y)ωc(y-z*)(4)
其中I(y)是点y的灰度,ωc(·)为一个加权函数,y离z*越近,权值越大;
D时空上下文模型的在线学习更新
将式(3)、式(4)代入式(2),得目标位置的似然估计值为
l ( z ) = Σ y ∈ F S h ( z - y ) I ( y ) ω c ( y - z * ) = h ( z ) ⊗ ( I ( z ) ω c ( z - z * ) ) - - - ( 5 )
表示卷积,对式(5)两边同时进行快速傅里叶变换,得到:
F(l(z))=F(h(z))⊙F(I(z)ωc(z-z*))(6)
⊙表示点乘,F(·)表示快速傅里叶变换,空间上下文模型为:
h ( z ) = F - 1 ( F ( l ( z ) ) F ( I ( z ) ω c ( z - z * ) ) ) - - - ( 7 )
F-1(·)表示逆FFT,FFT表示快速傅里叶变换,利用两个正傅里叶变换和一个逆傅里叶变换,求取空间上下文条件概率;
将时空上下文模型更新看成是一个IIR滤波过程,即t+1时刻新的空间上下文与t时刻空间上下文满足以下关系:
Ht+1(z)=(1-ρ)Ht(z)+ρht(z)(8)
在获取t+1时刻新的时空上下文关系后,用两个正傅里叶变换和一个逆傅里叶变换来加速计算新一帧中目标位置的置信度响应:
目标在下一帧的位置对应置信度响应最大的位置:
z t + 1 = arg max z ∈ S ( z t ) l t + z ( z ) - - - ( 10 ) .
CN201410841442.9A 2014-12-30 2014-12-30 基于时空上下文信息辅助的关键点稳定跟踪方法 Expired - Fee Related CN104537692B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410841442.9A CN104537692B (zh) 2014-12-30 2014-12-30 基于时空上下文信息辅助的关键点稳定跟踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410841442.9A CN104537692B (zh) 2014-12-30 2014-12-30 基于时空上下文信息辅助的关键点稳定跟踪方法

Publications (2)

Publication Number Publication Date
CN104537692A CN104537692A (zh) 2015-04-22
CN104537692B true CN104537692B (zh) 2016-02-17

Family

ID=52853211

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410841442.9A Expired - Fee Related CN104537692B (zh) 2014-12-30 2014-12-30 基于时空上下文信息辅助的关键点稳定跟踪方法

Country Status (1)

Country Link
CN (1) CN104537692B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104933735A (zh) * 2015-06-30 2015-09-23 中国电子科技集团公司第二十九研究所 一种基于时空上下文学习的实时人脸跟踪方法及系统
CN106023246A (zh) * 2016-05-05 2016-10-12 江南大学 一种基于局部敏感直方图的时空上下文跟踪方法
JP6755713B2 (ja) * 2016-05-25 2020-09-16 キヤノン株式会社 追尾装置、追尾方法及びプログラム
CN106127798B (zh) * 2016-06-13 2019-02-22 重庆大学 基于自适应模型的稠密时空上下文目标跟踪方法
CN106127766B (zh) * 2016-06-24 2018-12-25 国家新闻出版广电总局广播科学研究院 基于空间耦合关系和历史模型的目标跟踪方法
CN106203423B (zh) * 2016-06-26 2021-04-16 广东外语外贸大学 一种融合上下文探测的弱结构感知视觉目标跟踪方法
CN107045726B (zh) * 2017-04-27 2019-12-24 南京航空航天大学 一种基于时空上下文模型的目标跟踪方法
CN106960203B (zh) * 2017-04-28 2021-04-20 北京搜狐新媒体信息技术有限公司 一种面部特征点跟踪方法及系统
CN107451601A (zh) * 2017-07-04 2017-12-08 昆明理工大学 基于时空上下文全卷积网络的运动工件识别方法
CN109448023B (zh) * 2018-10-23 2021-05-18 武汉大学 一种卫星视频小目标实时跟踪方法
CN111107278B (zh) * 2018-10-26 2022-03-01 北京微播视界科技有限公司 图像处理方法、装置、电子设备及可读存储介质
CN111626082A (zh) * 2019-02-28 2020-09-04 佳能株式会社 检测装置和方法及图像处理装置和系统
CN110738685B (zh) * 2019-09-09 2023-05-05 桂林理工大学 一种融合颜色直方图响应的时空上下文跟踪方法
CN111815679B (zh) * 2020-07-27 2022-07-26 西北工业大学 一种基于双目相机的空间目标特征点丢失期间轨迹预测方法
CN113945207B (zh) * 2021-09-08 2023-11-24 杭州电子科技大学 一种基于目标三状态建模的遮挡目标多帧跟踪方法
CN116385497B (zh) * 2023-05-29 2023-08-22 成都与睿创新科技有限公司 用于体腔内的自定义目标追踪方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7751599B2 (en) * 2006-08-09 2010-07-06 Arcsoft, Inc. Method for driving virtual facial expressions by automatically detecting facial expressions of a face image
CN101763636B (zh) * 2009-09-23 2012-07-04 中国科学院自动化研究所 视频序列中的三维人脸位置和姿态跟踪的方法
CN103955950B (zh) * 2014-04-21 2017-02-08 中国科学院半导体研究所 一种利用关键点特征匹配的图像跟踪方法

Also Published As

Publication number Publication date
CN104537692A (zh) 2015-04-22

Similar Documents

Publication Publication Date Title
CN104537692B (zh) 基于时空上下文信息辅助的关键点稳定跟踪方法
CN105760811B (zh) 全局地图闭环匹配方法及装置
CN102222346B (zh) 一种车辆检测和跟踪方法
CN105856243A (zh) 一种移动智能机器人
CN105335986A (zh) 基于特征匹配和MeanShift算法的目标跟踪方法
CN104318782A (zh) 一种面向区域重叠的高速公路视频测速方法及系统
CN106952294B (zh) 一种基于rgb-d数据的视频跟踪方法
CN115451977A (zh) 车道标注数据的获取方法、计算机设备及存储介质
CN105758408A (zh) 局部地图构建方法及装置
AU2021255130B2 (en) Artificial intelligence and computer vision powered driving-performance assessment
Konrad et al. Localization in digital maps for road course estimation using grid maps
CN112017188B (zh) 一种空间非合作目标语义识别与重构方法
CN110246155A (zh) 一种基于模型交替的抗遮挡目标跟踪方法及系统
CN103577833A (zh) 基于运动模板的异常闯入检测方法
CN115407364A (zh) 点云地图处理方法、车道标注数据获取方法、设备及介质
Jia et al. Front-view vehicle detection by Markov chain Monte Carlo method
Xia et al. Extended object tracking with automotive radar using learned structural measurement model
KR20100041172A (ko) 영상 추적 장치의 이동표적 움직임 추적 방법
Zhu et al. Deepego: Deep instantaneous ego-motion estimation using automotive radar
KR20090093119A (ko) 움직이는 객체 추적을 위한 다중 정보의 융합 방법
Arnaud et al. Partial linear gaussian models for tracking in image sequences using sequential monte carlo methods
Xiu et al. Tracking algorithm based on the improved template matching
Stellet et al. Post processing of laser scanner measurements for testing advanced driver assistance systems
CN110673134A (zh) 一种雷达远近交替测距模式下航迹跟踪优化方法
CN115511970A (zh) 一种面向自主泊车的视觉定位方法

Legal Events

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

Granted publication date: 20160217

Termination date: 20181230