CN105548983A - 一种雷达组网识别目标方法 - Google Patents

一种雷达组网识别目标方法 Download PDF

Info

Publication number
CN105548983A
CN105548983A CN201510871256.4A CN201510871256A CN105548983A CN 105548983 A CN105548983 A CN 105548983A CN 201510871256 A CN201510871256 A CN 201510871256A CN 105548983 A CN105548983 A CN 105548983A
Authority
CN
China
Prior art keywords
target
attribute
classification
detection
radar
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
CN201510871256.4A
Other languages
English (en)
Other versions
CN105548983B (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.)
Ordnance Engineering College of PLA
Original Assignee
Ordnance Engineering College of PLA
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 Ordnance Engineering College of PLA filed Critical Ordnance Engineering College of PLA
Priority to CN201510871256.4A priority Critical patent/CN105548983B/zh
Priority claimed from CN201510871256.4A external-priority patent/CN105548983B/zh
Publication of CN105548983A publication Critical patent/CN105548983A/zh
Application granted granted Critical
Publication of CN105548983B publication Critical patent/CN105548983B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种雷达组网识别目标方法,属于雷达技术领域,该方法由s个雷达分别对探测目标的n个属性进行测量,形成s组雷达测量数据,并将探测目标归于m种可能目标中的一种,其中每种可能目标的每一个属性均具有特定的取值范围。本方法提高了识别系统的抗干扰能力,降低了识别难度,提高了识别准确度,该方法的识别分类过程是自上而下进行的,其相对于现有技术来说识别过程更加简单可靠,具有极强的推广应用价值。

Description

一种雷达组网识别目标方法
技术领域
本发明涉及雷达技术领域,尤其涉及一种雷达组网识别目标方法。
背景技术
现代战场环境日益复杂,在各种干扰存在的情况下,雷达回波得到的目标量测信息不确定性大,很难满足作战系统目标识别的需求。为此,现有技术中出现了基于雷达组网的目标融合识别方法。这种方法使用多个雷达同时对一个目标进行探测,并将探测结果进行融合,从而得出精确度优于单个雷达的测量结果。这种雷达组网有助于增强系统的抗干扰能力和环境适应能力,提高了目标识别的可信度及准确性,为作战指挥辅助决策提供了重要依据。
雷达组网探测目标方法的核心技术在于目标融合识别,即指利用多个传感器得到关于目标属性的数据并将这些数据进行融合,从而进行目标身份属性判别的过程。根据信息的抽象层次,目标融合识别可分为数据层融合、特征层融合和决策层融合。数据层融合直接融合来自多个传感器的原始观测信息,在融合的数据中进行特征提取并进行分类,从而得到目标的类别,具有较高的信息利用率。但数据层融合处理的信息量大,抗干扰能力差,且其仅对产生同类观测的传感器是适用的。决策层融合是较高层次的融合,因高层次融合的信息量小,且实现难度相对较小,所以目前目标识别大都采用此方法。但决策层融合损失了大量的特征信息,在量测数据受到干扰时,识别可信度降低。特征层融合属于中间层次,既保留了足够数量的重要信息,又实现了客观的信息压缩。
目前,在现有技术中,利用特征层融合进行目标识别的有关技术还处于空白状态。
发明内容
本发明要解决的技术问题是:针对上述现有技术的不足,提供一种雷达组网识别目标方法,该方法将特征层融合技术应用于雷达组网的目标识别过程中,提高了识别系统的抗干扰能力,降低了识别难度,提高了识别准确度,识别过程更加简单可靠。
为解决上述技术问题,本发明所采取的技术方案是:
一种雷达组网识别目标方法,该方法使用s个雷达分别对探测目标的n个属性进行测量,形成s组雷达测量数据,并将探测目标归于m种可能目标中的一种,其中每种可能目标的每一个属性均具有特定的取值范围;具体步骤包括:
(a)获得测量数据;
(b)根据测量数据计算每一个属性的信息增益;
(c)将具有最大信息增益的属性作为分类属性,并根据分类属性计算每组雷达测量数据对每个可能目标的隶属度,进而根据隶属度判定探测目标的所属类别;
(d)忽略已选取的分类属性,对剩余属性重复步骤(b)到步骤(c)的操作,直至将探测目标归结为唯一可能,该唯一可能即为对探测目标的识别结果。
具体地,信息增益按照式(1)计算:
(1),
其中,Gain(ck)表示第k个属性的信息增益,I为对探测目标进行识别的信息熵,E(ck)为利用ck划分当前数据集合所需要的信息熵;
式(1)中的I按照式(2)计算:
(2);
式(1)中的E(ck)按照式(3)计算:
(3),
其中,
(4),
(5)。
具体地,隶属度按照式(6)计算:
(6);
其中,
(7),
(8),
(9)。
进一步地,根据隶属度判定探测目标的所属类别的依据为:
时,判定探测目标为类别;当时,判定探测目标不是类别,而是其它类别。
本方法的原理可以简述如下:
一组雷达的量测数据可以表示为一个多维数组,其中每个数据都是某个属性变量下的相应取值,那么数组的维数由属性变量的个数决定。对同一目标观测,设雷达组网由部雷达构成,要求每部雷达产生的量测数据维数相同,都为维量测数据。表示成矩阵形式,可以得到维量测数据
假设有种可能的目标类别,对于量测数据集合,完全未知其类别,因此认为其属于任何一类的可能性均等,则对数据集合分类所需的信息熵为
依据目标数据库中已积累的数据,给出各类目标关于不同属性变量的量值及波动范围。设一个属性变量为,目标类别关于属性的标准值和标准偏差分别为,取。则目标关于属性的取值范围是
假设表示全体目标类别,这里取。则全体目标关于属性的取值范围是
若属性变量被选为分类属性,那么利用划分当前数据集合所需要的信息熵为
式中,是在分类特征下,对第组量测数据与目标类别之间距离的描述。
越小,表示第组量测数据属于目标类别的可能性越大;越大,表示第组量测数据属于目标类别的可能性越小。每组量测数据以不同的可能性属于不同的类别。
若属性变量为分类属性,此时状态的混乱程度用信息熵来度量。对归一化后,计算其熵值,即为利用划分第组量测数据的信息熵,该熵值表示分类后第组量测数据的不确定性。对组量测数据的信息熵求和,得到利用划分当前量测数据集合的信息熵。若分类后,数值较小,则表示利用划分当前量测数据集合后,该状态下的混乱程度较小,量测数据相对确定的属于某个类别。这正是分类识别要达到的目的。
利用属性变量作为分类属性,对当前分支节点数据集合进行分类,所得的信息增益为即为分类后信息熵的减少量。通过计算每一个属性变量的信息增益,将具有最大信息增益的属性作为分类属性,对当前节点数据集合进行分类。
若属性变量为分类属性,则第组量测数据对目标类别的隶属度为
其中,为第组量测数据与目标类别的关联系数,为第组量测数据与全体目标类别的关联系数,为目标类别的宽度。
表示了在当前分类属性下,第组量测数据所观测目标属于目标类别的程度。当时,判定所观测目标为目标类别;当时,判定所观测目标不是目标类别,而是其它类别。
为第组量测数据的权重,它表示对第组量测数据的信任程度。则判定第组量测数据是目标类别的可信度就可以定义为
对量测数据集合进行分类识别,判定识别结果为的可信度为。这样,正确识别结果的可信度越高,就表示准确性越高。若无事先说明,可以认为各雷达的量测数据权重相同,即取
基于以上原理,本方法的分类规则可以描述如下:
(1)由目标数据库给出各类别目标关于不同属性的量值及波动范围,各类别直接构成类决策树的叶节点;
(2)以全体量测数据集合创建第一个节点,即根节点;
(3)选取具有最大信息增益的属性作为分类属性,一旦一个属性变量在某节点处被选为分类属性,在该节点的任何后代上,选择分类属性时就不考虑该属性变量;
(4)根据对当前分支节点数据集合进行分类,产生类决策树的分支;
(5)当不同类别目标在同一属性变量上特性相近时,可能产生一组数据同时属于多类别目标的情况;由于雷达量测数据存在干扰,可能产生一组数据不属于任何类别的情况;以上情况可归纳为分类结果没有唯一确定;也就是说,经过一级分类不能得到识别结果;此时,将分类结果不唯一的量测数据组再次进行下一级分类,重复(3)、(4)过程,直到满足停止分类的条件;
(6)若某组量测数据中某个属性对应的值空缺,在选择分类属性时,只计算非空属性值分类所需的信息熵,忽略空缺值对分类属性的熵值的影响;对含有空缺值的数据进行分类时,若该分类属性下,量测数据空缺,则该组数据跳过本级分类,进入下一级分类,直到满足停止分类的条件;
(7)由于雷达量测数据不确定性大,会形成异类分支;该算法采用边构造边剪枝的方式,即当前节点产生相应的分支,暂时保留每个分支,进行下一级分类,若连续两次某分支相对其余分支为异类分支,则该分支视为无效,其他分支继续进行分类,直到确定所有量测数据类别;
(8)当且仅当以下条件之一成立时,递归步骤停止:
①所有有效的量测数据产生唯一、相同的分类结果,即到达类决策树的某叶节点;
②没有剩余属性变量可以用来对数据进行分类。
对于①,将分类结果判定为目标类别。对于②不能再进行分类时,选取最后一级分类产生的关联函数值中的最大值,将其相对应的类别判定为该级的分类结果。使用多数表决,将数据中绝大部分所属的类判定为目标类别。在分类完成后,计算Bi和B值即可得到识别结果的可信度,从而对识别的准确性进行评估。
采用上述技术方案所产生的有益效果在于:本发明方法将特征层融合技术应用于雷达组网的目标识别过程中,提高了识别系统的抗干扰能力,具有识别难度小、识别准确度高的优点;该方法的识别分类过程是自上而下进行的,其相对于现有技术,识别过程更加简单可靠,因而具有极强的推广应用价值。
附图说明
图1是本发明实施例的流程图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一:
假设现有雷达组网中三部雷达对一空中目标进行探测识别。已知三种目标类别。选取四项用于识别的属性变量,标准值分别为,标准偏差分别为,由先验知识给出三种目标类别的参数模板,如表1所示。
表1目标类别的标准值及标准偏差
每部雷达经时空配准和关联处理后的量测数据记录如表2所示。试确定该空中目标的类别。
表2雷达量测数据
首先,由表1先验信息,构成类决策树的叶节点。
其次,创建第一个节点。由式(2)得到,对量测数据集合分类所需的信息熵为。依据式(1)计算每一个属性变量的信息增益。将具有最大信息增益的属性作为分类属性,对当前节点数据集合依据式(6)分类。得到隶属度
由隶属度可以看出,第一组和第二组数据得到唯一的分类结果,即第一类目标。而第三组数据通过第一级分类未能做出有效识别。
再次,对第三组数据进行下一级分类,重复上述过程。计算除外每一个属性变量的信息增益,将具有最大信息增益的属性作为分类属性,对当前节点数据集合分类。得到隶属度
由隶属度可以看出,第三组数据得到唯一的分类结果,即第一类目标。至此,所有的量测数据产生唯一、相同的分类结果,即到达类决策树的叶节点,停止分类。识别结果为第一类目标。整个识别流程如图1所示。
最后,取,计算得到识别结果的可信度
实施例二:
三种目标类别的先验知识如表1所示。假设雷达组网的量测数据存在部分空缺值,将原量测数据任意设置空缺值得到表3。试确定该空中目标的类别。
表3雷达量测数据(含空缺值)
识别方法与实施例一略同,相同部分不再赘述。在第一个节点处,计算每一个属性变量的信息增益。由于存在空缺值,只计算非空属性值分类所需的信息熵。将具有最大信息增益属性作为分类属性,对当前节点数据集合分类。得到隶属度
由隶属度可以看出,第二组与第三组数据得到唯一的分类结果,即第一类目标。而第一组数据含有属性下的空缺值,则不参与本级分类,进行下一级分类,重复上述过程。计算除外每一个属性变量的信息增益,将具有最大信息增益的属性作为分类属性,对当前节点数据集合分类。得到隶属度
可以看出,第一组数据得到唯一的分类结果,即第一类目标。至此,所有的量测数据产生唯一、相同的分类结果,停止分类。识别结果为第一类目标,可信度为0.6506。
需要指出的是,以上具体实施方式只是本专利实现方案的具体个例,没有也不可能覆盖本专利的所有实现方式,因此不能视作对本专利保护范围的限定;凡是与以上案例属于相同构思的实现方案,均在本专利的保护范围之内。
通过以上实施例可见,本发明方法在特征层进行融合识别,相比决策层融合提高了信息利用率,能够对含有空缺值的量测数据进行正确分类,提高了目标识别结果的可信度。相比数据层融合,克服了数据量过大、抗干扰能力差的不足。并且,采用信息增益最大的属性作为分类属性进行分类,使得根节点到达叶节点的路径较短,提高了分类速度。
总之,本方法属于一种类决策树分类算法,它无需训练样本,是一种简单有效的特征层融合识别算法,普遍适用于战时雷达组网目标融合识别,特别是针对量测数据不准确、不完备的情况,能够对其进行合理的分类识别。此外,该方法不仅适用于雷达组网的目标识别,也可推广到其他多传感器的决策问题。

Claims (4)

1.一种雷达组网识别目标方法,其特征在于:由s个雷达分别对探测目标的n个属性进行测量,形成s组雷达测量数据,并将探测目标归于m种可能目标中的一种,其中每种可能目标的每一个属性均具有特定的取值范围;具体步骤包括:
(a)获得测量数据;
(b)根据测量数据计算每一个属性的信息增益;
(c)将具有最大信息增益的属性作为分类属性,并根据分类属性计算每组雷达测量数据对每个可能目标的隶属度,进而根据隶属度判定探测目标的所属类别;
(d)忽略已选取的分类属性,对剩余属性重复步骤(b)到步骤(c)的操作,直至将探测目标归结为唯一可能,该唯一可能即为对探测目标的识别结果。
2.根据权利要求1所述的雷达组网识别目标方法,其特征在于:所述信息增益按照式(1)计算:
(1),
其中,Gain(ck)表示第k个属性的信息增益,I为对探测目标进行识别的信息熵,E(ck)为利用ck划分当前数据集合所需要的信息熵;
式(1)中的I按照式(2)计算:
(2);
式(1)中的E(ck)按照式(3)计算:
(3),
其中,
(4),
(5)。
3.根据权利要求1所述的雷达组网识别目标方法,其特征在于:所述隶属度按照式(6)计算:
(6);
其中,
(7),
(8),
(9)。
4.根据权利要求3所述的雷达组网识别目标方法,其特征在于:所述根据隶属度判定探测目标的所属类别的依据为:
时,判定探测目标为类别
时,判定探测目标不是类别,而是其它类别。
CN201510871256.4A 2015-12-02 一种雷达组网目标识别方法 Expired - Fee Related CN105548983B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510871256.4A CN105548983B (zh) 2015-12-02 一种雷达组网目标识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510871256.4A CN105548983B (zh) 2015-12-02 一种雷达组网目标识别方法

Publications (2)

Publication Number Publication Date
CN105548983A true CN105548983A (zh) 2016-05-04
CN105548983B CN105548983B (zh) 2018-08-31

Family

ID=

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106054171A (zh) * 2016-05-27 2016-10-26 中国电子科技集团公司信息科学研究院 一种基于信息熵的多雷达节点自适应选择跟踪方法
CN107436424A (zh) * 2017-09-08 2017-12-05 中国电子科技集团公司信息科学研究院 一种基于信息增益的多雷达动态调节方法及装置
CN107561530A (zh) * 2017-07-27 2018-01-09 中国船舶重工集团公司第七二四研究所 一种基于协同被动探测多维信息的目标识别方法
CN107728122A (zh) * 2017-09-08 2018-02-23 中国电子科技集团公司信息科学研究院 一种基于信息熵的多雷达信息量的度量方法及装置
CN108169721A (zh) * 2017-11-24 2018-06-15 北京遥感设备研究所 一种雷达信号多目标检测方法
CN111308437A (zh) * 2020-02-27 2020-06-19 南京慧尔视智能科技有限公司 一种毫米波mimo交通雷达求熵解速度模糊的方法
CN112632463A (zh) * 2020-12-22 2021-04-09 中国航空工业集团公司沈阳飞机设计研究所 一种基于多属性的目标数据关联方法及装置
CN113687328A (zh) * 2021-09-14 2021-11-23 上海无线电设备研究所 一种弹载武器地面目标高分辨一维距离像识别方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070195680A1 (en) * 2006-02-14 2007-08-23 Intelliscience Corporation Methods and system for aggregating and using physical samples and data in a virtual environment
CN102663394A (zh) * 2012-03-02 2012-09-12 北京航空航天大学 基于多源遥感图像融合的大中型目标识别方法
CN103366184A (zh) * 2013-07-23 2013-10-23 武汉大学 基于混合分类器的极化sar数据分类方法及系统
CN104239712A (zh) * 2014-09-05 2014-12-24 电子科技大学 一种雷达抗干扰效能实时评估方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070195680A1 (en) * 2006-02-14 2007-08-23 Intelliscience Corporation Methods and system for aggregating and using physical samples and data in a virtual environment
CN102663394A (zh) * 2012-03-02 2012-09-12 北京航空航天大学 基于多源遥感图像融合的大中型目标识别方法
CN103366184A (zh) * 2013-07-23 2013-10-23 武汉大学 基于混合分类器的极化sar数据分类方法及系统
CN104239712A (zh) * 2014-09-05 2014-12-24 电子科技大学 一种雷达抗干扰效能实时评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘惠等: "基于模糊集理论和信息增益分析技术的分类算法", 《信号处理》 *
牛鹏飞等: "基于隶属度和D- S 理论的雷达目标识别方法", 《微计算机信息》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106054171A (zh) * 2016-05-27 2016-10-26 中国电子科技集团公司信息科学研究院 一种基于信息熵的多雷达节点自适应选择跟踪方法
CN106054171B (zh) * 2016-05-27 2021-07-13 中国电子科技集团公司信息科学研究院 一种基于信息熵的多雷达节点自适应选择跟踪方法
CN107561530A (zh) * 2017-07-27 2018-01-09 中国船舶重工集团公司第七二四研究所 一种基于协同被动探测多维信息的目标识别方法
CN107436424A (zh) * 2017-09-08 2017-12-05 中国电子科技集团公司信息科学研究院 一种基于信息增益的多雷达动态调节方法及装置
CN107728122A (zh) * 2017-09-08 2018-02-23 中国电子科技集团公司信息科学研究院 一种基于信息熵的多雷达信息量的度量方法及装置
CN108169721A (zh) * 2017-11-24 2018-06-15 北京遥感设备研究所 一种雷达信号多目标检测方法
CN111308437A (zh) * 2020-02-27 2020-06-19 南京慧尔视智能科技有限公司 一种毫米波mimo交通雷达求熵解速度模糊的方法
CN111308437B (zh) * 2020-02-27 2022-03-22 南京慧尔视智能科技有限公司 一种毫米波mimo交通雷达求熵解速度模糊的方法
CN112632463A (zh) * 2020-12-22 2021-04-09 中国航空工业集团公司沈阳飞机设计研究所 一种基于多属性的目标数据关联方法及装置
CN113687328A (zh) * 2021-09-14 2021-11-23 上海无线电设备研究所 一种弹载武器地面目标高分辨一维距离像识别方法

Similar Documents

Publication Publication Date Title
CN104699755B (zh) 一种基于数据挖掘的智能化多目标综合识别方法
CN105069093B (zh) 一种基于嵌入式索引的水文时间序列相似性搜索方法
CN109508000A (zh) 异构多传感器多目标跟踪方法
CN109299812B (zh) 一种基于深度学习模型和knn实时校正的洪水预测方法
CN104835178A (zh) 一种低信噪比运动小目标的跟踪与识别的方法
CN104239712B (zh) 一种雷达抗干扰效能实时评估方法
CN105488770A (zh) 一种面向对象的机载激光雷达点云滤波方法
CN116166960B (zh) 用于神经网络训练的大数据特征清洗方法及系统
CN102393912A (zh) 一种基于不确定推理的目标综合识别方法
CN103888541A (zh) 一种融合拓扑势和谱聚类的社区发现方法及系统
CN104268532B (zh) 网络环境下大规模图像与视频数据高效处理的方法
CN105678047A (zh) 融合经验模态分解降噪和复杂网络分析的风场表征方法
CN111222575B (zh) 一种基于hrrp目标识别的klxs多模型融合方法及系统
CN109765332A (zh) 基于隔离森林的变压器异常值实时检测和故障诊断方法
CN103942325A (zh) 一种融合气候分区思想的海陆气候事件关联规则挖掘方法
CN109150830B (zh) 一种基于支持向量机和概率神经网络的层次入侵检测方法
Mishra et al. Rainfall-runoff modeling using clustering and regression analysis for the river brahmaputra basin
CN114355298B (zh) 一种雷达复合调制脉冲信号识别方法
CN104199006B (zh) 一种基于多假设联合分布滤波器的随机集跟踪方法
CN103778342A (zh) 一种基于区间数和证据理论的辐射源识别算法
CN103413027B (zh) 一种社会网络重叠社区发现方法的评价方法
Bouaita et al. Rainfall estimation from MSG images using fuzzy association rules
CN105548983A (zh) 一种雷达组网识别目标方法
CN111768101B (zh) 一种顾及物候特征的遥感耕地变化检测方法及系统
CN105548983B (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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Shang Chaoxuan

Inventor after: Han Zhuangzhi

Inventor after: Wang Pin

Inventor after: Wang Xuefei

Inventor after: Hu Wenhua

Inventor after: Xie Hui

Inventor after: Liu Limin

Inventor before: Han Zhuangzhi

Inventor before: Shang Chaoxuan

Inventor before: Wang Pin

Inventor before: Wang Xuefei

Inventor before: Hu Wenhua

Inventor before: Xie Hui

Inventor before: Liu Limin

GR01 Patent grant
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: 20180831

Termination date: 20211202