CN103278225B - 基于图像技术检测的皮带称重装置 - Google Patents

基于图像技术检测的皮带称重装置 Download PDF

Info

Publication number
CN103278225B
CN103278225B CN201310219904.9A CN201310219904A CN103278225B CN 103278225 B CN103278225 B CN 103278225B CN 201310219904 A CN201310219904 A CN 201310219904A CN 103278225 B CN103278225 B CN 103278225B
Authority
CN
China
Prior art keywords
belt
height
image
weighing area
hot spot
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
CN201310219904.9A
Other languages
English (en)
Other versions
CN103278225A (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.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN201310219904.9A priority Critical patent/CN103278225B/zh
Publication of CN103278225A publication Critical patent/CN103278225A/zh
Application granted granted Critical
Publication of CN103278225B publication Critical patent/CN103278225B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于图像技术的皮带称重装置,该装置由皮带转速传感器、激光源、CCD成像系统、预装有软件的计算机及连接CCD成像系统和计算机的传输光缆组成;激光源含有多个激光器产生一组平行激光束将称重区域离散成连续小块;计算机控制CCD成像系统按图像采间隔连续获取称重区域光斑图像,每获取一幅图像对图像进行处理,提取光斑像点位置,根据光斑像点位置计算每个小块的重量相加成称重区域的重量,累加每次采样间隔内称重区域的重量即得到当前采样时刻皮带上总的物料量;该装置称重区域大小可灵活调整,依据每次存贮的称重区域物料重量,可查询任意时间段的运料量及生成运料量报表。

Description

基于图像技术检测的皮带称重装置
技术领域
本发明总体涉及传送带散装物料称重领域,尤其涉及利用图像技术检测皮带运送物料重量的皮带称重装置。
背景技术
在煤矿生产领域,皮带运输是煤矿原煤运输的重要方式之一,对皮带运输的原煤进行称重是遏制煤矿超能力生产,保障生产安全和加强税收管理等的有效措施。有多种皮带称重的装置,然而,目前没有已知皮带称重装置的设计是基于图像技术的。皮带运输原煤的称重装置有两种:即基于重力检测原理的电子皮带称和基于射线吸收原理检测的核子皮带称。
电子皮带称称重有一些内在的问题,电子皮带称称重时原煤在皮带之上,称重架在皮带之下,皮带夹在原煤与称重架之间,皮带上原煤的重量必须通过称重架和拉紧的皮带来称重,称重架受力状态和皮带张力直接影响着电子皮带称的测量精度,由于称重架受力状态和皮带张力随着皮带的运动在不断改变,电子皮带称要保持稳定的测量精度,必须经常性的维护使电子皮带称工作在皮带张力恒定、皮带不跑偏、称架不积灰、托辊校准好、称重传感器参数匹配等运行条件,因此,电子皮带称的维护工作很频繁,维护工作量很大。另外皮带在各处的张力并不相同,其变化也不同,电子皮带称往往要求在皮带张力小及张力变化小的地方安装,对安装位置有要求。精度稍高的电子皮带称称架都比较庞大和复杂,且称架多为钢结构件,因此电子皮带称往往比较笨重,安装不容易。
核子皮带称通过射线吸收原理进行称重,物料特性的变化,如品种、成分、含水量、在皮带上断面形状的变化等待都对测量精度有影响,如称重物料由煤改为矿石甚至煤矸石时,需要重新进行校准才能保证测量精度,而校准复杂,耗力大,耗时长,另外,放射源申办繁琐,需要配置辐射剂量专用测量设备,工作人员享受特殊保健并需要定期检查身体等。
两种称重装置价格都是比较昂贵,这些原因都导致了企业生产成本的增加,经济效益的缩小。
需要一种解决或至少改善现有技术中固有的一个或多个问题的皮带原煤称重装置。
发明内容
因此,本发明的目的是提供一种基于图像技术的皮带称重装置,该装置成本低,组件少,安装简单,维护工作量小,并且安装灵活,测量精度满足需求。
根据一种实例形式,提供一种基于图像技术的皮带称重装置,包括皮带转速传感器、激光源、CCD成像系统、预装有软件的计量管理计算机及连接CCD成像系统和所述计算机的传输线缆;皮带转速传感器提供皮带速度,激光源由至少2个以上在皮带上方沿皮带宽度线成直线排列的激光器组成,激光源向称重区域投射平行激光束,平行激光束与皮带面有大于0度且小于90度的夹角;称重区域大小为w×v×Δt,其中w为皮带宽度,v为皮带转动速度,Δt为图像采样间隔;优选地,称重区域位于皮带下方有托辊支撑处。CCD成像系统位于皮带中线位置的上方,光轴垂直于皮带面,CCD成像系统视场能覆盖大于或等于称重区域的范围,CCD成像系统为面阵成像系统,优选地,CCD成像系统为方阵CCD。计量管理计算机控制CCD成像系统按Δt间隔连续获取称重区域光斑图像,对图像进行处理,提取光斑像点位置,计算称重区域物料重量及皮带上总的物料重量。
在进一步特定的但非限制性的形式中,激光源在一个带有透明视窗的长形箱体内,长形箱体内有两条平行滑轨,其中一侧滑轨带有刻度标尺,以提供激光器间的间距,激光器可在平行滑轨上任意滑动并在任意位置固定,激光器的滑动范围至少在皮带宽度范围内,平行滑轨可旋转并带动所有激光器一起旋转而激光器自身保持相对静止。
在进一步特定的但非限制性的形式中,若装置应用于煤矿井下,传输线缆是防爆型的,CCD成像系统是防爆型的,激光源是防爆型的。
在进一步特定的但非限制性的形式中,装置称重方法包括如下步骤:
A.测量装置安装完毕后,首先对测量系统光斑像点位置与光斑相对于皮带面来说的高度之间的等量关系进行拟合;
B.计量开始后,计算机控制CCD成像系统每隔Δt时间摄取一幅称重区域光斑的图像,并将图像信号传送给计算机;
C.计算机对光斑图像进行处理,得到光斑像点位置,依据步骤A拟合的等量关系计算光斑处物料的高度;
D.计算由光斑将称重区域离散成的每个小区块的物料高度hi
E.由下式计算称重区域物料的重量ΔG,并将其存贮;
ΔG = Σ i = 1 N h i · w i · v · Δt · ρ
其中:其中wi为小区块的宽度,v为皮带转动速度,Δt为CCD成像系统摄取图像间隔,ρ为物料的体密度;
F.累加每次E的结果,即为从计量开始到当前采样时刻皮带上总的重量Gt,Gt=Gt-1+ΔG,Gt-1为t-1时刻皮带上物料的总重量。
在进一步特定的但非限制性的形式中,步骤A等量关系由下式描述:
h = ad d + b
其中,h为光斑所在处的物料相对于皮带面的高度,a,b为测量系统参数,d为光斑像点相对于参考像点的像素偏移量,其中,参考像点为光斑位于皮带面时的像点。
有利的是,测量系统参数a=H0其中,f为镜头焦距,L为光斑与光轴与皮带面交点连线的距离,H0为镜头光心到皮带面的垂直距离,θ为激光束与皮带面夹角。
优选地,测量系统参数a,b由最小二乘法拟合得到。
在进一步特定的但非限制性的形式中,步骤D小区块的物料高度为该区块相邻两光斑物料高度的平均值;步骤D小区块的物料高度为该区块相邻光斑在前一采样时刻和当前采样时刻物料高度的平均值,其中前一采样时刻为t-1时刻,当前时刻为t时刻,时间差为Δt即采样间隔。
在进一步特定的但非限制性的形式中,装置的称重区域大小是可调整的,当皮带速度改变时,可调整采样间隔时间来保持称重区域的不变,也可根据精度要求调整称重区域大小。
在进一步特定的但非限制性的形式中,存贮每次得到的称重区域物料的重量,可查询任意时间段的运料量及生成运料报表。
附图说明
通过以下说明,附图实施例变得显而已见,其仅以结合附图描述的至少一种优选但非限制性实施例的示例方式给出。
图1示出了一种基于图像技术的皮带称重装置结构示意图;
图2示出了一种基于图像技术的皮带称重装置光斑处物料高度测量原理示意图;
图3示出了一种基于图像技术的皮带称重装置小区块物料高度计算示意图;
图4示出了一种基于图像技术的皮带称重装置小区块物料高度计算示意图;
附图符号:
1-激光源2-摄像系统3-光缆4-计算机5-激光器6-透明窗7-箱体8-旋转装置9-平行滑轨10-激光束11-转速传感器12-皮带13-托辊
具体实施方式
图1为基于图像技术的皮带称重装置示意图。该装置由皮带传动传感器11、激光源1、CCD成像系统2、传输线缆3和预装有处理软件的计量管理计算机4组成。其中,考虑到CCD成像系统与计量管理计算机间的距离、信号的传输质量及信息量大小,传输线缆3优选光缆。激光源1至少含有2个激光器5产生两束相互平行的激光束10,激光器5在一个一面带有透明视窗6的长形箱体7内,长形箱体7位于皮带12上方,透明视窗6面朝向皮带,长形箱体7的长边与皮带宽度线平行,长形箱体7内有两条平行滑轨9,其中一侧滑轨带有刻度标尺,激光器5可在平行滑轨9上任意滑动并在任意位置固定,激光器5的滑动范围至少在皮带宽度范围内,平行滑轨9可旋转并带动激光器5旋转,激光器5采用红色点光源,圆型光斑模式。固定好激光器并记录下激光器间的间距,调节激光器光路相对于水平面的出射角使平行激光光线照射在称重区域,由于称重区域皮带在物料重力作用下会下沉,为减小皮带下沉给称重带来的误差,称重区域优选在下方有托辊13支撑的皮带面处,称重区域大小为w×v×Δt,其中w为皮带宽度,v为皮带转动速度,Δt为图像采样间隔。
CCD成像系统2在皮带上方位于皮带中线位置,其光轴垂直于皮带面,CCD视场能覆盖大于或等于称重区域,CCD成像系统优选方阵CCD成像系统。
若装置用在煤矿井下,传输线缆应是矿用防爆型的,激光源可是矿用本安型或矿用隔爆型,CCD成像系统可是矿用本安型或矿用隔爆型。
当皮带上物料高度不同时,物料上的光斑离CCD镜头中心的距离也不同,光斑在CCD上的成像点位置也不同,通过光斑在CCD上成像点位置可以确定光斑点处物料的高度。
图2示出了一种基于图像技术的皮带称重装置光斑点处物料高度测量原理示意图。如图2所示,当皮带为空时,激光器103发射的激光束在皮带面202形成光斑201,在CCD面阵上的像点为204,当皮带上有物料时,激光束发射单元103发射的激光束则在物料面上形成激光光斑207,在CCD面阵上的像点为205,由激光光斑位置和CCD成像系统的几何位置关系,根据成像原理,光斑点207处的物料高度h为:
h = ad d + b , a=H0 b = fctgθ + fL H 0
其中,d为像点204与像点205的距离,f为镜头焦距,L为光斑201到光轴与皮带面交点206的连线距离,H0为镜头光心203到皮带面的垂直距离,即203与206之间连线距离。
将每个激光器发射的激光束在空皮带时形成的光斑距光轴与皮带面交点206的连线距离代入式中即得该激光器所发射的激光束形成的光斑所在处的物料高度测量模型。
当摄像机以Δt时间间隔采样光斑图像时,采样间隔时间Δt内对应皮带上大小为w×v×Δt区域,此即为称重区域,此时光斑像点位置反映了称重区域物料高度信息,通过光斑图像计算出称重区域物料的高度即可计算出采样间隔Δt时间内称重区域物料重量,由于称重区域物料表面是起伏的,意味着物料表面高度是变化的,因此用平行光束将称重区域离散成多个小块来计算,每个小块的大小为Δt×v×wi,其中wi为相邻两激光器间的间距,连续采样并累加每个采样间隔时间Δt内物料的重量,即可得到从计时开始到当前时刻皮带上的物料总重量。
为提高计量精度并且避免精确定位光斑在空皮上的位置和CCD光轴的位置以及CCD镜头精确聚焦,测量装置安装后对每个激光器所产生的光斑的物料高度测量模型中的参数a和b用实际的物理位置进行拟合,拟合步骤如下:
A1.对于每个光斑,先摄像取光斑位于皮带面上时的光斑图像,再获取位于测量范围内至少3个不同高度位置处的光斑图像,优选地,为5个高度位置;
A2.基于步骤A1获取的图像,确定不同高度位置处光斑像点位置;
A3.以位于皮带面上的光斑像点位置为参考点,计算其它光斑像点与之的像素偏移量;
A4.用最小二乘法求解得到参数a和b;
测量装置标定后,即可进行称重。
具体称重过程包括如下步骤:
A.首先用最小二乘法拟合光斑像点位置与光斑相对于皮带面来说的高度之间的等量关系;
B.计量开始后,计算机控制CCD成像系统每隔Δt时间摄取一幅称重区域光斑的图像,并将图像信号传送给计算机;
C.计算机对光斑图像进行处理,得到光斑像点位置,依据步骤A拟合的等量关系计算光斑处物料的高度;
D.计算由光斑将称重区域离散成的每小区块的物料高度hi
E.由下式计算称重区域物料的重量ΔG,并将其加以存贮,
ΔG = Σ i = 1 N h i · w i · v · Δt · ρ
其中:ρ为物料的体密度;
F.累加每次E的结果,即为从计量开始到当前采样时刻皮带上总的运料量Gt,Gt=Gt-1+ΔG,Gt-1为t-1时刻时皮带上物料的总重量。
采样间隔Δt大小可根据测量精度来选取,在满足测量精度的前提下,应尽可能选大,以避免系统处理数据量过大,另外,当皮带速度改变时,可通过调整Δt来维持称重区域恒定。
步骤E数据保存的好处在于,可以利用其查询任一时间段内的物料重量并形成运送量报表等。
小区块的物料高度可为该区块当前采样时刻相邻两光斑处物料高度的平均值,图3示出了这种计算示意图。小区块1的煤位高度为光斑301的煤位高度和光斑302的煤位高度的平均,小区块2的煤位高度为光斑302的煤位高度和光斑303的煤位高度的平均,依此类推。
为进一步提高测量精度,小区块的物料高度可为该区块相邻光斑在前一采样时刻和当前采样时刻物料高度的平均值;其中前一采样时刻为t-1时刻,当前时刻为t时刻,时间差为Δt即采样间隔。图4示出了这种计算示意图。小区块1的煤位高度为光斑302、光斑304、光斑301和光斑3034个点的煤位高度的平均值,依此类推其后小区块的煤位高度。

Claims (4)

1.一种基于图像检测技术的皮带称重装置,包括激光源、CCD成像系统、预装有软件的计算机、连接CCD成像系统和所述计算机的传输线缆及提供皮带速度的转速传感器,其特征在于:
激光源由至少2个以上在皮带上方沿皮带宽度成直线排列的激光器组成,所述激光源在一个带有透明视窗的长形箱体内,长形箱体内有两条平行滑轨,其中一侧滑轨带有刻度标尺,以提供所述激光器间的间距,所述激光器可在平行滑轨上任意滑动并在任意位置固定,所述激光器的滑动范围至少在皮带宽度范围内,平行滑轨可旋转并带动所有所述激光器一起旋转而所述激光器自身保持相对静止;所述激光源向称重区域投射平行激光束,平行激光束与皮带面有大于0度且小于90度的夹角;所述CCD成像系统位于所述皮带中线位置的上方,光轴垂直于皮带面,视场能覆盖大于或等于称重区域的范围,通过光缆连接到计算机;所述计算机控制CCD成像系统按Δt间隔连续获取称重区域光斑图像,对图像进行处理,提取光斑像点位置,计算称重区域物料重量及至当前时刻皮带上总的物料重量;所述装置的测量系统存在如下等量关系:h为光斑所在处的物料相对于皮带面的高度,a,b为系统参数,用最小二乘法拟合得到;d为光斑像点相对于参考像点的像素偏移量,所述参考像点为该光斑位于皮带面时所成的像点。
2.根据权利要求1所述的基于图像检测技术的皮带称重装置,其特征在于所述装置通过关系式计算称重区域物料重量,其中wi为平行激光束将称重区域离散成的小区块的宽度,υ为皮带转动速度,Δt为CCD成像系统摄取图像间隔,ρ为物料的体密度,N为小区块数,hi为小区块内的物料高度;
由Gt=Gt-1+ΔG计算直至当前时刻皮带上总的物料重量,其中Gt-1为t-1时刻皮带上物料的总重量。
3.根据权利要求2所述的基于图像检测技术的皮带称重装置,其特征在于所述装置计算所述小区块的物料高度时取相邻两光斑处物料高度的平均值为所述小区块的物料高度。
4.根据权利要求2所述的基于图像检测技术的皮带称重装置,其特征在于所述装置计算所述小区块的物料高度时取相邻两光斑在前一采样时刻和当前采样时刻物料高度的平均值为所述小区块的物料高度。
CN201310219904.9A 2013-06-05 2013-06-05 基于图像技术检测的皮带称重装置 Expired - Fee Related CN103278225B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310219904.9A CN103278225B (zh) 2013-06-05 2013-06-05 基于图像技术检测的皮带称重装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310219904.9A CN103278225B (zh) 2013-06-05 2013-06-05 基于图像技术检测的皮带称重装置

Publications (2)

Publication Number Publication Date
CN103278225A CN103278225A (zh) 2013-09-04
CN103278225B true CN103278225B (zh) 2016-05-04

Family

ID=49060811

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310219904.9A Expired - Fee Related CN103278225B (zh) 2013-06-05 2013-06-05 基于图像技术检测的皮带称重装置

Country Status (1)

Country Link
CN (1) CN103278225B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103837222B (zh) * 2014-03-12 2015-11-18 北京工业职业技术学院 皮带秤图像校验方法和计量监测系统
CN104006869A (zh) * 2014-06-04 2014-08-27 上海诺尚信息技术有限公司 一种用于皮带输送的激光称量系统及其称量方法
CN105258777A (zh) * 2014-07-18 2016-01-20 北京斯凯尔工业自动化技术研究所 皮带秤空带监测方法、装置和系统
CN104534990A (zh) * 2015-01-13 2015-04-22 中国矿业大学(北京) 基于激光测距的皮带物料测量系统
CN106743339A (zh) * 2016-12-08 2017-05-31 中国神华能源股份有限公司 一种自动调节胶带运输机速度的方法和装置
CN108254281B (zh) * 2016-12-29 2021-07-20 湖南长天自控工程有限公司 一种皮带机物料水分率的在线检测方法与装置
CN106949854A (zh) * 2017-04-26 2017-07-14 苏州睿牛机器人技术有限公司 一种激光视觉传感器及其检测方法
CN107128661B (zh) * 2017-05-17 2020-01-14 安徽理工大学 非接触式皮带输送计重及偏离警示装置
CN107255498B (zh) * 2017-06-02 2020-06-26 西安理工大学 基于双目视觉的输送带物料重量测量系统及测量方法
CN107727207A (zh) * 2017-09-19 2018-02-23 安凯 一种光学皮带秤及其测量方法
CN108593062B (zh) * 2018-05-25 2020-10-30 太原科技大学 基于条纹投影的带式输送机在线称重设备的测试方法
CN108414064B (zh) * 2018-06-01 2023-06-09 成都睿畜电子科技有限公司 遥感式料塔称重设备及称重方法
CN110376598A (zh) * 2019-08-14 2019-10-25 深圳市镭神智能系统有限公司 一种视觉系统与激光雷达融合装置及融合系统
DE102019215614A1 (de) * 2019-10-11 2021-04-15 BSH Hausgeräte GmbH Bestimmen eines Objektgewichts eines Kühlobjekts
CN110902313A (zh) * 2019-12-05 2020-03-24 天津成科传动机电技术股份有限公司 皮带轮廓检测、输送带流量检测方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5331118A (en) * 1992-11-27 1994-07-19 Soren Jensen Package dimensional volume and weight determination system for conveyors
CN1397786A (zh) * 2002-08-16 2003-02-19 清华大学 在线测量物料重量的装置
CN201464023U (zh) * 2009-06-05 2010-05-12 北京斯凯尔工业科技有限公司 电子皮带秤
CN103090936A (zh) * 2013-01-21 2013-05-08 中国计量学院 一种新型煤炭运输量测量装置和测量方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080035390A1 (en) * 2006-08-09 2008-02-14 Wurz David A Dimensioning and weighing system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5331118A (en) * 1992-11-27 1994-07-19 Soren Jensen Package dimensional volume and weight determination system for conveyors
CN1397786A (zh) * 2002-08-16 2003-02-19 清华大学 在线测量物料重量的装置
CN201464023U (zh) * 2009-06-05 2010-05-12 北京斯凯尔工业科技有限公司 电子皮带秤
CN103090936A (zh) * 2013-01-21 2013-05-08 中国计量学院 一种新型煤炭运输量测量装置和测量方法

Also Published As

Publication number Publication date
CN103278225A (zh) 2013-09-04

Similar Documents

Publication Publication Date Title
CN103278225B (zh) 基于图像技术检测的皮带称重装置
CN204613057U (zh) 煤炭灰分在线测量系统
CN103308127B (zh) 基于图像技术的皮带称重方法
US7570372B2 (en) Optical device for measuring the thickness of an at least partially transparent medium
US7098409B2 (en) Apparatus for weighing materials online
CN103217252B (zh) 一种移动容器式高精度微压检测装置
EP2944983B1 (en) Scanning imaging systems
CN104949620A (zh) 用于光学测量装置的校正设备及校正方法
US9612213B2 (en) Automatic z-correction for basis weight sensors
CN104006827A (zh) 一种惯组标定用北向基准稳定性评估测试方法
JP2007263818A (ja) 厚さ計測装置の調整方法及びその装置
CN110207625A (zh) 一种阴极板平面度动态检测系统及方法
CN102650602B (zh) 放射线检查装置
Xu et al. Investigation on dynamic calibration for an optical-fiber solids concentration probe in gas-solid two-phase flows
CN110902313A (zh) 皮带轮廓检测、输送带流量检测方法、装置、设备及存储介质
Shimizu et al. A liquid-surface-based three-axis inclination sensor for measurement of stage tilt motions
CN105143914A (zh) 用于在近程区域和远程区域中光学间距测量的方法
CN106249427A (zh) 一种基于激光成像的光轴调整方法
Guedri et al. Bulk material flow measurement based only on a smart camera fixed above a moving belt conveyor
KR102114744B1 (ko) 디지털 선로 게이지 및 이를 이용한 데이터 측정 및 관리 방법
CN204228590U (zh) 用于计算物料灰分的参数检测装置及物料灰分的测量设备
Song et al. Drone photogrammetry for accurate and efficient rock joint roughness assessment on steep and inaccessible slopes
CN1940544B (zh) 高精度在线灰分测量装置
US6493418B1 (en) Method for precisely measuring the weight of mass material and nucleonic weigher employing this method
CN208765674U (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160504

Termination date: 20190605