CN104199821A - 一种基于概要结构的流数据立方体构建方法 - Google Patents

一种基于概要结构的流数据立方体构建方法 Download PDF

Info

Publication number
CN104199821A
CN104199821A CN201410323039.7A CN201410323039A CN104199821A CN 104199821 A CN104199821 A CN 104199821A CN 201410323039 A CN201410323039 A CN 201410323039A CN 104199821 A CN104199821 A CN 104199821A
Authority
CN
China
Prior art keywords
dimension
data
flow data
value
data cube
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
CN201410323039.7A
Other languages
English (en)
Other versions
CN104199821B (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 University City College ZUCC
Original Assignee
Zhejiang University City College ZUCC
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 University City College ZUCC filed Critical Zhejiang University City College ZUCC
Priority to CN201410323039.7A priority Critical patent/CN104199821B/zh
Publication of CN104199821A publication Critical patent/CN104199821A/zh
Application granted granted Critical
Publication of CN104199821B publication Critical patent/CN104199821B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24Querying
    • G06F16/245Query processing
    • G06F16/2458Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2474Sequence data queries, e.g. querying versioned data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/22Indexing; Data structures therefor; Storage structures
    • G06F16/2228Indexing structures

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Computational Linguistics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Fuzzy Systems (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

本发明涉及一种基于概要结构的流数据立方体构建方法(Sketch Cube),包括以下步骤:把任意维度组合映射成唯一自然数;根据上下限单调原则对维度组合裁剪;在类线性空间中保存有效数据单元信息;构建时间序列索引。本发明的有益效果为:可以在类线性存储空间中满足实时分析要求,并且可以有效地控制正确率。

Description

一种基于概要结构的流数据立方体构建方法
技术领域
本发明属于计算机数据统计与分析领域,尤其涉及一种基于概要结构的流数据立方体构建方法。 
背景技术
移动互联网的发展带来了越来越多的流数据,如路由器中的IP地址包、微博微信等社交媒体中的短文本内容、用户商品浏览行为日志信息等,流数据是一组顺序、大量、快速且连续到达的数据序列,一般情况下,其可被视为随时间延续而无限增长的动态数据的集合。 
目前,常见的流数据分析方法有基于数据抽样和基于数据压缩两种方式。这些方法把数据看成整体而进行全局的分析,而并没有把流数据看成是一串连续的多维度数据而在不同层次进行分析挖掘。对不同层次和大小的数据单元进行聚合操作分析,可以挖掘出应用在特定时间下的使用场景,因而具有重要意义。传统数据仓库联机分析OLAP处理工具把数据仓库中的记录投射到不同空间中的数据立方进行分析,从而提供多角度综合分析能力。OLAP结合概要结构可以提供在实时流数据中的数据立方方法。 
发明内容
针对现有技术中存在的不足,本发明的目的在于提供一种基于概要结构的流数据立方体构建方法,以解决现有海量流数据下的数据立方体在存储空间、处理速度等方面的缺点。 
本发明的技术方案如下: 
一种基于概要结构的流数据立方体构建方法,包括以下步骤: 
步骤S100:对流数据的任意维度组合通过配对函数映射成一唯一标识,该标识为一自然数; 
步骤S100进一步包括有: 
步骤S110:将所需构建的流数据立方体中的相同属性值按照频率的高低进行排序,并按照自然数从小向大的次序进行映射; 
步骤S120:对于所需构建的流数据立方体中的不同属性值,按其数值的基数进行倒排; 
步骤S200:提供一概要统计模型,其根据所需构建的流数据立方体的数据分布特性和存储模型特点,依据上下限单调原则对各维度组合裁剪; 
步骤S200进一步包括有: 
步骤S210:提供一裁剪模型,对单个维度进行统计,由所需构建的流数据立方体聚合的单调性获得裁剪公式; 
步骤S220:根据最终效果,判断使用单维度裁剪模型或者组合维度的裁剪模型; 
步骤S300:以固定时间片为单位进行索引,使用的主键为流数据所在的ID和时间值,对应的值为概要统计模型中的多维度数组,并将有效的数据单元信息保存在类线性空间中; 
步骤S300进一步包括: 
步骤S310:不同的流数据在特定的时间段内,对所有维度组合产生的唯一标识进行最小计数概要统计; 
步骤S320:对于统计结果先放置在计算机的内存,然后定时存储于NoSql数据库中; 
步骤S400:实时查询时,给定流数据所在的ID、时间段以及相关属性 值,返回相应流数据立方体中的度量值。 
步骤S400进一步包括: 
步骤S410:不同的流数据在特定的时间段内,所有的维度都会在先映射后添加到Hash表中,给定流数据所在的ID、时间段以及相关属性值,生成含有这些属性的所有维度组合; 
步骤S420:将这些维度组合映射成为映射值后,找出各映射值所对应的在Hash表中的位置信息; 
步骤S430:根据步骤S420中获得的位置信息,查询并使用最小技术概要方法得到该组合的统计数据; 
步骤S440:将这些维度组合按照所要求的属性,进行组合或者抵消后形成具体的属性,然后将结果返回。 
本发明的有益效果是:对于互联网中出现的越来越复杂的流数据,将其看作是一连串的多维度数据并当作整体进行全局分析,将流数据的各维度组合进行一一映射生成一永不冲突且唯一的标识,并以连续的自然数依小到大的方式呈现,如此实现了在类线性存储空间中的实时分析,且可以有效地控制正确率。 
说明书附图
图1是本发明的流程框图; 
图2是本发明的操作框图; 
图3是小颗粒时间累加得到大颗粒时间度量值示意图。 
具体实施方式
以下描述和附图充分地示出本发明的具体实施方案,以使本领域的技术人员能够实践它们。实施例仅代表可能的变化。本发明的实施方案的范 围包括权利要求书的整个范围,以及权利要求书的所有可获得的等同物。在本文中,本发明的这些实施方案可以被单独地或总地用术语“发明”来表示,这仅仅是为了方便,并且如果事实上公开了超过一个的发明,不是要自动地限制该应用的范围为任何单个发明或发明构思。 
实施例1 
本发明提出一种基于概要结构的流数据立方体构建方法,其用于解决海量流数据下的数据立方体在存储空间、处理速度等方面的缺点。 
构建概要结构的流数据立方体包括如下步骤: 
步骤S100:提出一种可扩展的数据单元标识模型,对流数据的任意维度组合通过配对函数映射生成唯一标识,且永不冲突。 
该步骤S100进一步包括: 
步骤S110:将所需要的流数据立方体中的相同属性值按照其频率的高低排序,按自然数从小向大的次序进行映射; 
步骤S120:对于所需要的流数据立方体中的不同属性值,按其数值的基数进行倒排,使得属性值分布广泛的值放置靠后,可以有效减少配对函数产生的值。 
在步骤S100中运用一数据单元标识器进行映射,其功能是把不同的维度组合映射成一唯一整数,且该数据单元标识器支持维度的扩展,在新增维度或修改维度值都不会影响原有的映射值。 
数据单元标识器为一数据单元标识函数,具体包括: 
算法1:数据单元标识算法DCI 
输入:流数据记录r=(a1,a2,…,ai,m,t)。 
输出:r所有维度组合标识。 
①将记录r中的维度映射成(n1,n2,…,nn
②根据(n1,n2,…,nn)计算出其对应的所有组合set<Combination>。 
③将set<result>设置为空; 
④对于set<Combination>中的任何一个元素x,做以下操作: 
⑤对这个元素进行求Pairing Function(x)操作后,将结果存入到set<result>集中; 
⑥结束循环; 
⑦将set<result>集返回; 
算法1(Data Cell Identifier)把流数据中每条记录映射成由小到大的自然数,通过递归函数得到与之相关的所有维度组合,并使用康托尔配对函数获得每个组合唯一标识。数据单元标识算法输出所有相关数据单元标识,且配对函数保证标识符可扩展不冲突。 
具体地,如,在流数据R=(A1,A2,…,An,M,T)中,|Ai|表示第i维的基数。 
把记录R的维度Ai映射成连续的自然数。 
对于上述的记录R产生的n个自然数Ni,这些自然数形成了一个集合S,使用配对函数对S中的任意非空子集生成唯一自然数。 
继续进行一第二步骤S200,为一数据分析步骤,该步骤S200为提出一种改进的概要统计模型,根据数据分布特性和存储模型特点,有效裁剪结果无效的数据单元。不但可以提高计算能力和存储空间效率,而且可大幅度提高统计精确度。 
提出一种通过概要技术对流数据进行OLAP统计的方法,存储和计数由数据单元标识模型产生的维度组合信息。 
计数最小概要模型CM Sketch(Count-Min Sketch)是一个使用互相独立哈希函数族函数统计流数据元素出现频率的模型。 
该模型进一步包括: 
下面所示为w×d的二维数组。 
其中d表示互相独立哈希族函数的个数,w表示每个哈希函数的映射范围,如下式所示 
hk:{1...N}→{1...w},(1≤k≤d) 
同时,对于2n元素的集合只需n位就能表示。因此,设计包含d个函数的互相独立哈希函数族,可使用个不同元素两两组合表示。 
在该数据表示时,给定数据集合如下式所示 
SeedSet = { 1,2 , . . . , n | n d + 1 }
从SeedSet中随机取不同的元素a和b,设计哈希函数为式 
ha,b(Z)=a×z+b 
由ha,b(Z)产生的坐标信息是变长的,下式可把变长值规约到长度为w大小的数组中: 
ha,b(Z)=((a×z+b)%p)%w,    (1)其中p为大质数。 
上式(1)为Carter-Wegman模型,以保证元素的平均分布。 
该步骤S200进一步包括: 
步骤S210:提供一种裁剪模型,对单个维度进行统计,由数据立方体聚合的单调性给出裁剪公式; 
步骤S220:根据最终效果,判断使用单维度裁剪模型或者组合维度的裁剪模型。 
继续进行一第三步骤S300,以固定时间片为单位进行索引,使用的主键为流数据所在的ID和时间值,对应的值为概要统计模型中的多维度数 组,并将有效的数据单元信息保存在类线性空间中; 
步骤S300进一步包括: 
步骤S310:不同的流数据在特定的时间段内,对所有维度组合产生的唯一标识进行最小计数概要统计; 
步骤S320:对于统计结果先放置在计算机的内存,然后定时存储于NoSql数据库中。 
然后进行步骤S400,为在实时查询时,给定流数据所在的ID、时间段以及相关属性值,返回相应流数据立方体中的度量值。 
步骤S400进一步包括: 
步骤S410:不同的流数据在特定的时间段内,所有的维度都会在先映射后添加到Hash表中,给定流数据所在的ID、时间段以及相关属性值,生成含有这些属性的所有维度组合; 
步骤S420:将这些维度组合映射成为映射值后,找出各映射值所对应的在Hash表中的位置信息; 
步骤S430:根据步骤S420中获得的位置信息,查询并使用最小技术概要方法得到该维度组合的统计数据; 
步骤S440:将这些维度组合按照所要求的属性,进行组合或者抵消后形成具体的属性,然后将结果返回。 
具体如下,对于该概要统计模型的填充操作为:对于第t次到达的元素c,计数最小概要模型更新操作如下式所示。 
&ForAll; 1 &le; j &le; d : CM [ j , h j ( i t ) ] &LeftArrow; CM [ j , h j ( i t ) ] + c t - - - ( 2 )
该概要统计模型的更新时间复杂度为
该概要统计模型中统计元素ai在CM Sketch中的操作如下式: 
a ^ i = min 1 &le; j &le; d CM [ j , h j ( i ) ] - - - ( 3 )
即通过哈希函数族中的每个函数计算该元素在对应数组中的下标值,获得所有可能值中的最小值即为该元素的估计值。其查询时间复杂度为O(1)。 
该概要统计模型的压缩率为: 
对于(A1,A2,…,An)数据立方体,其所有数据单元个数T为: 
T = &Pi; 1 n ( | A i | + 1 ) - - - ( 4 )
则Sketch Cube模型的压缩率P为下式: 
P = &Pi; 1 n ( | A i | + 1 ) w &times; d &times; 100 % - - - ( 5 )
流数据具体有内在时序性,对流数据挖掘用倾斜时间窗口TTW(Tilted-Time Window)在不同时间粒度(Multiple Time Granularities)上进行分析,即对流数据挖掘用倾斜时间窗口在不同时间粒度上进行分析。Sketch Cube按时间片段对元素组合进行统计把结果放入计数最小概要模型。如此不但可以提高计算能力和存储空间效率,而且可大幅度提高统计精确度。 
SketchCube设计的存储结构可支持任意时间粒度的组合,其合并公式如下。 
&ForAll; 1 &le; j &le; d : CM [ j , h j ( i t a + t b ) ] &LeftArrow; CM [ j , h j ( i t a ) ] + CM [ j , h j ( i t b ) ] - - - ( 6 )
对于给定哈希函数族,相同维度组合在不同时间中的映射地址相等,可单次扫描小颗粒时间累加得到大颗粒时间度量值(参见图3)。 
上述实施例仅仅是为清楚地说明本发明创造所作的举例,而并非对本发明具体实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所引伸出的任何显 而易见的变化或变动仍处于本权利要求的保护范围之中。 

Claims (6)

1.一种基于概要结构的流数据立方体构建方法,其特征在于,包括以下步骤:
步骤S100:对流数据的任意维度组合通过配对函数映射成一唯一标识,该标识为一自然数;
步骤S200:提供一概要统计模型,其根据所需构建的流数据立方体的数据分布特性和存储模型特点,依据上下限单调原则对各维度组合裁剪;
步骤S300:以固定时间片为单位进行索引,使用的主键为流数据所在的ID和时间值,对应的值为概要统计模型中的多维度数组,并将有效的数据单元信息保存在类线性空间中;
步骤S400:实时查询时,给定流数据所在的ID、时间段以及相关属性值,返回相应流数据立方体中的度量值。
2.根据权利要求1所述的构建方法,其特征在于,其中的步骤S100进一步包括有:
步骤S110:将所需构建的流数据立方体中的相同属性值按照频率的高低进行排序,并按照自然数从小向大的次序进行映射;
步骤S120:对于所需构建的流数据立方体中的不同属性值,按其数值的基数进行倒排。
3.根据权利要求1所述的构建方法,其特征在于,其中的步骤S200进一步包括有:
步骤S210:提供一裁剪模型,对单个维度进行统计,由所需构建的流数据立方体聚合的单调性获得裁剪公式;
步骤S220:根据最终效果,判断使用单维度裁剪模型或者组合维度的裁剪模型。
4.根据权利要求1所述的构建方法,其特征在于,其中的步骤S300进一步包括:
步骤S310:不同的流数据在特定的时间段内,对所有维度组合产生的唯一标识进行最小计数概要统计;
步骤S320:对于统计结果先放置在计算机的内存,然后定时存储于NoSql数据库中。
5.根据权利要求1所述的构建方法,其特征在于,其中的步骤S400进一步包括:
步骤S410:不同的流数据在特定的时间段内,所有的维度都会在先映射后添加到Hash表中,给定流数据所在的ID、时间段以及相关属性值,生成含有这些属性的所有维度组合;
步骤S420:将这些维度组合映射成为映射值后,找出各映射值所对应的在Hash表中的位置信息;
步骤S430:根据步骤S420中获得的位置信息,查询并使用最小技术概要方法得到该组合的统计数据;
步骤S440:将这些维度组合按照所要求的属性,进行组合或者抵消后形成具体的属性,然后将结果返回。
6.根据权利要求1所述的构建方法,其特征在于,所述步骤S100中,运用一数据单元标识器进行映射,该数据单元标识器支持维度扩展。
CN201410323039.7A 2014-07-08 2014-07-08 一种基于概要结构的流数据立方体构建方法 Expired - Fee Related CN104199821B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410323039.7A CN104199821B (zh) 2014-07-08 2014-07-08 一种基于概要结构的流数据立方体构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410323039.7A CN104199821B (zh) 2014-07-08 2014-07-08 一种基于概要结构的流数据立方体构建方法

Publications (2)

Publication Number Publication Date
CN104199821A true CN104199821A (zh) 2014-12-10
CN104199821B CN104199821B (zh) 2017-06-20

Family

ID=52085114

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410323039.7A Expired - Fee Related CN104199821B (zh) 2014-07-08 2014-07-08 一种基于概要结构的流数据立方体构建方法

Country Status (1)

Country Link
CN (1) CN104199821B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016134579A1 (zh) * 2015-02-28 2016-09-01 中兴通讯股份有限公司 一种查找流数据的方法和查询服务器
CN106844713A (zh) * 2017-02-07 2017-06-13 北京微影时代科技有限公司 一种数据立方体生成的方法及装置
CN107203635A (zh) * 2017-06-07 2017-09-26 南开大学 一种基于最小略图的流模式下有向标签图的略图构建方法
CN108241692A (zh) * 2016-12-26 2018-07-03 北京国双科技有限公司 数据的查询方法及装置
CN112328464A (zh) * 2021-01-05 2021-02-05 北京基调网络股份有限公司 指标数据存储、相关性分析方法及计算机可读存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5799300A (en) * 1996-12-12 1998-08-25 International Business Machines Corporations Method and system for performing range-sum queries on a data cube
US20030093424A1 (en) * 2001-09-10 2003-05-15 Seok-Ju Chun Dynamic update cube and hybrid query search method for range-sum queries
CN1564160A (zh) * 2004-04-22 2005-01-12 重庆市弘越科技有限公司 建立及查询多维数据立方体的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5799300A (en) * 1996-12-12 1998-08-25 International Business Machines Corporations Method and system for performing range-sum queries on a data cube
US20030093424A1 (en) * 2001-09-10 2003-05-15 Seok-Ju Chun Dynamic update cube and hybrid query search method for range-sum queries
CN1564160A (zh) * 2004-04-22 2005-01-12 重庆市弘越科技有限公司 建立及查询多维数据立方体的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LIN C X, EL AT.: "Text cube: Computing IR measures for multidimensional text database analysis", 《EIGHTH IEEE INTERNATIONAL CONFERENCE ON DATA MINING》 *
ROSENBERG A L, ET AL.: "Efficient pairing functions--and why you should care", 《INTERNATIONAL PARALLEL AND DISTRIBUTED PROCESSING SYMPOSIUM》 *
李爱平等: "基于Dwarf的数据流立方体的研究与实现", 《计算机研究与发展》 *
韩矞: "基于网络安全数据流的数据立方体技术的研究与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016134579A1 (zh) * 2015-02-28 2016-09-01 中兴通讯股份有限公司 一种查找流数据的方法和查询服务器
CN105989141A (zh) * 2015-02-28 2016-10-05 中兴通讯股份有限公司 一种查找流数据的方法和查询服务器
CN105989141B (zh) * 2015-02-28 2021-02-23 中兴通讯股份有限公司 一种查找流数据的方法和查询服务器
CN108241692A (zh) * 2016-12-26 2018-07-03 北京国双科技有限公司 数据的查询方法及装置
CN106844713A (zh) * 2017-02-07 2017-06-13 北京微影时代科技有限公司 一种数据立方体生成的方法及装置
CN107203635A (zh) * 2017-06-07 2017-09-26 南开大学 一种基于最小略图的流模式下有向标签图的略图构建方法
CN107203635B (zh) * 2017-06-07 2020-08-11 南开大学 一种基于最小略图的流模式下有向标签图的略图构建方法
CN112328464A (zh) * 2021-01-05 2021-02-05 北京基调网络股份有限公司 指标数据存储、相关性分析方法及计算机可读存储介质
CN112328464B (zh) * 2021-01-05 2021-03-26 北京基调网络股份有限公司 指标数据存储、相关性分析方法及计算机可读存储介质

Also Published As

Publication number Publication date
CN104199821B (zh) 2017-06-20

Similar Documents

Publication Publication Date Title
Agarwal Range searching
CN104199821A (zh) 一种基于概要结构的流数据立方体构建方法
Fox et al. Spatio-temporal indexing in non-relational distributed databases
CN103064689B (zh) 一种动态生成报表页面的实现方法及系统
Kane et al. An optimal algorithm for the distinct elements problem
CN104166651B (zh) 基于对同类数据对象整合的数据搜索的方法和装置
CN107193967A (zh) 一种多源异构行业领域大数据处理全链路解决方案
US6980980B1 (en) Summary-detail cube architecture using horizontal partitioning of dimensions
CN105205146B (zh) 一种计算微博用户影响力的方法
CN102306176B (zh) 一种基于数据仓库内在特征的olap关键词查询方法
CN107066476A (zh) 一种基于物品相似度的实时推荐方法
CN104794221A (zh) 一种基于业务对象的多维数据分析系统
CA2361242A1 (en) A calculation engine for use in olap environments
US9153051B2 (en) Visualization of parallel co-ordinates
CN105488231A (zh) 一种基于自适应表维度划分的大数据处理方法
CN105930388B (zh) 一种基于函数依赖关系的olap分组聚集方法
CN103927337B (zh) 用于联机分析处理中关联关系的数据处理方法和装置
CN103123650A (zh) 一种基于整数映射的xml数据库全文索引方法
Cuzzocrea et al. Semantics-aware advanced OLAP visualization of multidimensional data cubes
CN110134879A (zh) 一种基于差分隐私保护的兴趣点推荐算法
Wang et al. Continuous visible k nearest neighbor query on moving objects
CN110096646A (zh) 品类关联信息的生成及其视频推送方法和相关设备
Xie et al. OLAP over probabilistic data cubes I: Aggregating, materializing, and querying
Wu et al. Efficient evaluation of object-centric exploration queries for visualization
CN110019252A (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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170620

Termination date: 20200708

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