CN104834966A - 一种基于蚁群分类算法的渔场预报方法 - Google Patents
一种基于蚁群分类算法的渔场预报方法 Download PDFInfo
- Publication number
- CN104834966A CN104834966A CN201510177021.5A CN201510177021A CN104834966A CN 104834966 A CN104834966 A CN 104834966A CN 201510177021 A CN201510177021 A CN 201510177021A CN 104834966 A CN104834966 A CN 104834966A
- Authority
- CN
- China
- Prior art keywords
- ant
- data
- fishing ground
- fishing
- cpue
- 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
Links
- 238000013277 forecasting method Methods 0.000 title abstract 3
- 238000000034 method Methods 0.000 claims abstract description 38
- 230000007613 environmental effect Effects 0.000 claims abstract description 20
- 238000003306 harvesting Methods 0.000 claims abstract description 11
- 238000005516 engineering process Methods 0.000 claims abstract description 5
- 230000008569 process Effects 0.000 claims description 19
- 238000004519 manufacturing process Methods 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 9
- 238000012804 iterative process Methods 0.000 claims description 6
- 230000002159 abnormal effect Effects 0.000 claims description 5
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims description 5
- 229930002868 chlorophyll a Natural products 0.000 claims description 5
- 241000257303 Hymenoptera Species 0.000 claims description 3
- 238000000546 chi-square test Methods 0.000 claims description 3
- 235000013305 food Nutrition 0.000 claims description 3
- 230000000750 progressive effect Effects 0.000 claims description 3
- 239000013535 sea water Substances 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 abstract 1
- 241000238366 Cephalopoda Species 0.000 description 9
- 241000251468 Actinopterygii Species 0.000 description 4
- 238000002790 cross-validation Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 description 2
- 101100077702 Escherichia coli (strain K12) mog gene Proteins 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000003016 pheromone Substances 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于蚁群分类算法的渔场预报方法,包括以下步骤:利用近实时的卫星遥感技术得到的环境数据,并将环境数据和渔获数据进行预处理;对环境数据进行有监督的离散化处理;通过蚁群算法中的AntMiner方法,经训练后得到一系列的渔场判别规则,利用得到的渔场判别规则;将实时或现场环境数据代入渔场判别规则后即可得到实时的渔场预报数据和分布图。
Description
技术领域
本发明涉及渔场预报技术领域,特别是涉及一种基于蚁群分类算法的渔场预报方法。
背景技术
西北太平洋是中国远洋渔船从事大规模商业性捕捞的重要作业海域,柔鱼是该海域的主要捕捞对象之一。据FAO统计,2006-2010年日本在北太平洋捕捞柔鱼年捕捞量达0.8~7.4万t,平均为3.6万t。中国自20世纪90年代以来,鱿钓生产规模呈扩大趋势,近年来在西北太平洋海域的作业渔船数达200~280艘,年捕捞量在3.6~11.3万t,在中国远洋渔业生产中占重要地位。北太平洋柔鱼通常分为4个群体,分别为中部秋生群、东部秋生群、西部冬春生群和中东部冬春生群,其中在170°E以西为冬春生群,是中国鱿钓船的传统捕捞对象。但柔鱼CPUE分布及其渔场的形成和变化受海洋环境因子的影响较大,且渔场分布具有区域性、时空易变性等特点,传统渔业捕捞中,对渔场判断的方法较为落后,难以快速、准确确定渔场位置,束缚了柔鱼在该海域生产效率和捕捞产量。
发明内容
本发明所要解决的技术问题是提供一种基于蚁群分类算法的渔场预报方法,能够快速准确确定渔场位置,降低渔场寻鱼时间、提高渔获产量。
本发明解决其技术问题所采用的技术方案是:提供一种基于蚁群分类算法的渔场预报方法,包括以下步骤:
(1)利用近实时的卫星遥感技术得到的环境数据,并将环境数据和渔获数据进行预处理;
(2)对环境数据进行有监督的离散化处理;
(3)通过蚁群算法中的AntMiner方法,经训练后得到一系列的渔场判别规则,利用得到的渔场判别规则;
(4)将实时或现场环境数据代入渔场判别规则后即可得到实时的渔场预报数据和分布图。
所述步骤(1)中所述渔获数据包括日期、经度、纬度、产量、作业次数和平均日产量,时间分辨率为天,计算CPUE值:其中ΣCatch为1°×1°网格中一个月的总产量,ΣFishingdays为当月总的船天数,并将渔区按CPUE值的分位数划分为2类,分别定义为高CPUE渔区和低CPUE渔区;所述环境数据包括海域的海水表层温度、表层叶绿素a浓度、海面高度异常、SST梯度强度、chl梯度强度和SST异常数据;将渔获数据和环境数据按1°×1°的空间分辨率进行网格化处理。
所述步骤(2)具体包括:选择显著性水平,根据自由度得到对应的卡方值阈值,根据卡方分析方法将具有最小卡方值的相邻区间连续迭代合并,直到满足停止准则。
所述步骤(3)具体包括:利用Ant-miner算法模仿蚂蚁觅食行为进行分类规则的构造,每只蚂蚁的寻路过程相当于规则的生成过程,每次蚂蚁选择一个方向前进的过程相当于蚂蚁选择一个属性节点加入规则前件的过程,以上是一个迭代过程,每次迭代后将规则所正确覆盖的事例从训练集中剔除,在下一个迭代过程中是在更小的训练集中寻找新规则,这一过程重复进行直至找到的规则能够覆盖绝大部分训练集,得到一系列的渔场判别规则。有益效果
由于采用了上述的技术方案,本发明与现有技术相比,具有以下的优点和积极效果:本发明利用近实时的卫星遥感技术得到的环境数据,经过有监督的离散化后,通过蚁群算法中的AntMiner方法,经训练后得到一系列的渔场判别规则。利用这些规则,经代入实时或现场环境数据后即可实时地得到实时的渔场预报数据和分布图,从而能够快速准确确定渔场位置,降低渔场寻鱼时间、提高渔获产量。
附图说明
图1是不同离散化时显著性水平与CPUE分位点下交叉验证的精度分布图;
图2是不同离散化时显著性水平与CPUE分位点下交叉验证的精度提升值分布图;
图3是实施例中北太平洋柔鱼渔场预报与实际生产叠加图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
本发明的实施方式涉及一种基于蚁群分类算法的渔场预报方法,包括以下步骤:
(一)训练数据的准备和处理方法
渔获数据时间范围为2002~2011年,海域的空间范围为144°E-170°E、34°N-50°N,数据内容包括日期、经度、纬度、产量、作业次数和平均日产量。时间分辨率为天,计算CPUE值,即
其中,ΣCatch为1°×1°网格中一个月的总产量,ΣFishingdays为当月总的船天数,CPUE的单位为t/(船·天)。将渔区按CPUE的分位数划分为2类,分别定义为高CPUE渔区(High CPUE)和低CPUE渔区(Low CPUE)。
环境数据包括该海域的海水表层温度(SST)、表层叶绿素a浓度(chl)、海面高度异常(SSHA)、SST梯度强度(SSTG)、chl梯度强度(chlG)、SST异常(SSTA)数据。前两者均取自http://oceancolor.gsfc.nasa.gov/的月平均数据,海面高度异常数据来自于http://www.aviso.oceanobs.com/的月平均数据。为了使环境数据与渔获数据匹配,在此将生产数据和环境数据按1°×1°的空间分辨率进行网格化处理。其中,梯度强度可应用求梯度幅值(GradientMagnitude,简称GM)的公式来进行。设当前渔区的上下左右四个网格点的海表温度为SSTi±1,j和SSTi,j±1,以差分的格式表示表温梯度为
其中,i和j分别为当前渔区格点的列和行的索引值,Δy表示i-1行和i+1行之间的距离,Δx表示j-1列和j+1列之间的距离,单位为km。
同理可得到叶绿素a浓度的梯度为
其中,chli±1,j和chli,j±1为当前渔区的上下左右四个网格点的叶绿素a浓度。
SST异常通过公式(4)得到
其中,SSTi,j为第i行j列的SST值,分别为第i行j列多年SST平均值。
(二)解释变量的离散化
ChiMerge是有监督、自底向上的数据离散化,根据卡方分析方法将具有最小卡方值的相邻区间连续迭代合并,直到满足停止准则。要离散化某一属性,首先选择显著性水平α,根据自由度(输出变量的类别数-1)得到对应的卡方值阈值。当类别和属性独立时,有1-α的可能性使计算得到的卡方值会小于阈值,而大于阈值的卡方值说明属性和类别之间不是相互独立的,因此不能合并。离散化输出变量CPUE较简单,是通过确定一个0到1之间的分位数完成。
环境属性变量在不同显著水平α下离散程度不同,训练后得到的规则也不同。同样CPUE取不同的分位数q,也会对训练结果和渔场精度的验证产生不同的影响。因此,通过以下方法筛选出一对最佳(α,q)组合:设渔区总个数为N,划分高、低CPUE的分位数为q,意味着在测试数据集中,高、低CPUE渔区出现的概率分为q和1-q,则高CPUE渔场个数为NH=N(1-q),低CPUE渔场个数NL=Nq。在不经任何训练、按此高、低CPUE固有的概率随机为各个渔区分配渔场类型(High或Low)的情况下,仍然可以得到渔场精度
假设经训练后的预报精度为pα,q,显然,单纯以pα,q作为衡量训练的效果,是不合理的。而应从pα,q中去除pq'的影响,得到计算精度提升值
Δpα,q=pα,q-pq' (6)
当Δpα,q取最大值时,即可确定出(α,q)的最佳组合。
(三)机器学习过程
利用Ant-miner算法模仿蚂蚁觅食行为进行分类规则的构造,每只蚂蚁的寻路过程相当于规则的生成过程,每次蚂蚁选择一个方向前进的过程相当于蚂蚁选择一个属性节点加入规则前件的过程。以上是一个迭代过程,每次迭代后将规则所正确覆盖的事例从训练集中剔除,因此在下一个迭代过程中是在更小的训练集中寻找新规则,这一过程重复进行直至找到的规则能够覆盖绝大部分训练集,得到一系列的IF-THEN规则。其伪代码如下所示:
表1 Ant-Miner实现的步骤
在渔场环境中,假设环境属性节点i具有j值,记为(i|j),则(i|j)被蚂蚁选中的概率为
其中,a是所有环境属性值,bi是属性i的所有离散值,I为没有被蚂蚁选择过的属性集。
τij(t)为节点(i|j)在t时刻的信息素。当一条规则Rule完成,之后与规则相关的属性节点的信息素都将被更新,即
τij(t+1)=τij+τij×Q,i|j∈Rule,(8)
其中,Q为规则质量标准,为规则应用到训练集实例时的敏感度与特异度的乘积表示,具体如下:
其中TruePos是被规则覆盖,所属类与规则类相同的记录数目;FalsePos是被规则覆盖,所属类与规则类不同的记录数目;FalseNeg是未被规则覆盖,所属类别与规则类一致的记录数目;TrueNeg为不被规则覆盖,所属类与规则类不同的记录数目。
ηi,j为属性节点(i|j)的启发式函数,公式表示为
其中,infoTij为节点的信息熵,由下列公式表示
其中,k是数据集中类别数,|Ti,j|是数据集中所属性i的值为j的实例总数,为第i个属性值为j且属于w类的实例数。
现结合具体实施例进一步说明本发明:
(1)预报模型训练
设蚂蚁数量60只,每条规则实例数(No_Rules_Converg)为5,未覆盖的实例最大数(Max_Uncovered_Cases)为10,计算最大迭代次数1500。选取2002-2011年共10年的北太平洋柔鱼产量数据共计1024条记录作为训练集,设定不同的CPUE分位点qj=j×10-1(j=1,2,…,9)和不同的ChiMerger离散化显著性水平αi=i×10-3(i=1,2,…,10),对环境数据进行离散化并与分级的CPUE数据匹配,构成训练数据集,利用10折交叉验证测试模型的分类精度,得到分布图(如图1所示)。再利用公式(6)计算得到各个点的精度提升值其分布如图2所示。
从图2可以看出,当CPUE分位点为0.8和0.3时、均出现了精度提升值的高峰带,但为了避免预报时高CPUE渔场出现的范围过大,所以CPUE分位点取0.8,且当ChiMerger离散化显著性水平取0.006时,精度提升出现了一个高峰值13.82%,其10折交叉验证的精度为81.82%,标准差为±0.23%。因此将(0.006,0.8)作为离散化显著性水平和CPUE分位点,对训练数据集进行离散化并依Ant-Miner原理重新进行训练,得到表2的IF-THEN规则集。
表2 训练得到的渔场预报规则
(2)预报结果
将2012年6-12月份的环境数据代入表2中的规则,可得到各个月份每个渔区的预报渔场,经与实际生产渔场叠加后,如图3所示。
Claims (4)
1.一种基于蚁群分类算法的渔场预报方法,其特征在于,包括以下步骤:
(1)利用近实时的卫星遥感技术得到的环境数据,并将环境数据和渔获数据进行预处理;
(2)对环境数据进行有监督的离散化处理;
(3)通过蚁群算法中的AntMiner方法,经训练后得到一系列的渔场判别规则,利用得到的渔场判别规则;
(4)将实时或现场环境数据代入渔场判别规则后即可得到实时的渔场预报数据和分布图。
2.根据权利要求1所述的基于蚁群分类算法的渔场预报方法,其特征在于,所述步骤(1)中所述渔获数据包括日期、经度、纬度、产量、作业次数和平均日产量,时间分辨率为天,计算CPUE值:其中ΣCatch为1°×1°网格中一个月的总产量,ΣFishingdays为当月总的船天数,并将渔区按CPUE值的分位数划分为2类,分别定义为高CPUE渔区和低CPUE渔区;所述环境数据包括海域的海水表层温度、表层叶绿素a浓度、海面高度异常、SST梯度强度、chl梯度强度和SST异常数据;将渔获数据和环境数据按1°×1°的空间分辨率进行网格化处理。
3.根据权利要求1所述的基于蚁群分类算法的渔场预报方法,其特征在于,所述步骤(2)具体包括:选择显著性水平,根据自由度得到对应的卡方值阈值,根据卡方分析方法将具有最小卡方值的相邻区间连续迭代合并,直到满足停止准则。
4.根据权利要求1所述的基于蚁群分类算法的渔场预报方法,其特征在于,所述步骤(3)具体包括:利用Ant-miner算法模仿蚂蚁觅食行为进行分类规则的构造,每只蚂蚁的寻路过程相当于规则的生成过程,每次蚂蚁选择一个方向前进的过程相当于蚂蚁选择一个属性节点加入规则前件的过程,以上是一个迭代过程,每次迭代后将规则所正确覆盖的事例从训练集中剔除,在下一个迭代过程中是在更小的训练集中寻找新规则,这一过程重复进行直至找到的规则能够覆盖绝大部分训练集,得到一系列的渔场判别规则。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510177021.5A CN104834966B (zh) | 2015-04-15 | 2015-04-15 | 一种基于蚁群分类算法的渔场预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510177021.5A CN104834966B (zh) | 2015-04-15 | 2015-04-15 | 一种基于蚁群分类算法的渔场预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104834966A true CN104834966A (zh) | 2015-08-12 |
CN104834966B CN104834966B (zh) | 2018-07-24 |
Family
ID=53812840
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510177021.5A Expired - Fee Related CN104834966B (zh) | 2015-04-15 | 2015-04-15 | 一种基于蚁群分类算法的渔场预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104834966B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105654210A (zh) * | 2016-02-26 | 2016-06-08 | 中国水产科学研究院东海水产研究所 | 一种利用海洋遥感多环境要素的集成学习渔场预报方法 |
CN105787591A (zh) * | 2016-02-26 | 2016-07-20 | 中国水产科学研究院东海水产研究所 | 一种利用多尺度环境特征的渔场预报方法 |
CN106844777A (zh) * | 2017-03-05 | 2017-06-13 | 何钊荣 | 一种出海钓鱼信息分享与大数据挖掘系统 |
CN111240306A (zh) * | 2020-04-26 | 2020-06-05 | 南京市产品质量监督检验院 | 一种自适应配电变压器故障诊断系统及其诊断方法 |
CN112330029A (zh) * | 2020-11-08 | 2021-02-05 | 上海海洋大学 | 一种基于多层convLSTM的渔场预测计算方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103235982A (zh) * | 2013-04-16 | 2013-08-07 | 中国水产科学研究院东海水产研究所 | 一种基于贝叶斯网络模型的渔场预报方法 |
-
2015
- 2015-04-15 CN CN201510177021.5A patent/CN104834966B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103235982A (zh) * | 2013-04-16 | 2013-08-07 | 中国水产科学研究院东海水产研究所 | 一种基于贝叶斯网络模型的渔场预报方法 |
Non-Patent Citations (4)
Title |
---|
HU SOLANKI 等: "Fishery forecast using OCM chlorophyll concentration and AVHRR SST: validation results off Gujarat coast, India", 《INTERNATIONAL JOURNAL OF REMOTE SENSING》 * |
刘红岩: "《商务智能方法与应用》", 31 May 2013, 清华大学出版社 * |
崔雪森 等: "基于朴素贝叶斯的西北太平洋柔鱼渔场预报模型的建立", 《中国海洋大学学报》 * |
戴芹 等: "基于蚁群优化分类规则挖掘的遥感图像分类研究", 《计算机工程与应用》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105654210A (zh) * | 2016-02-26 | 2016-06-08 | 中国水产科学研究院东海水产研究所 | 一种利用海洋遥感多环境要素的集成学习渔场预报方法 |
CN105787591A (zh) * | 2016-02-26 | 2016-07-20 | 中国水产科学研究院东海水产研究所 | 一种利用多尺度环境特征的渔场预报方法 |
CN105787591B (zh) * | 2016-02-26 | 2019-08-20 | 中国水产科学研究院东海水产研究所 | 一种利用多尺度环境特征的渔场预报方法 |
CN106844777A (zh) * | 2017-03-05 | 2017-06-13 | 何钊荣 | 一种出海钓鱼信息分享与大数据挖掘系统 |
CN111240306A (zh) * | 2020-04-26 | 2020-06-05 | 南京市产品质量监督检验院 | 一种自适应配电变压器故障诊断系统及其诊断方法 |
CN112330029A (zh) * | 2020-11-08 | 2021-02-05 | 上海海洋大学 | 一种基于多层convLSTM的渔场预测计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104834966B (zh) | 2018-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
He et al. | Mining transition rules of cellular automata for simulating urban expansion by using the deep learning techniques | |
Valavanis et al. | Modelling of essential fish habitat based on remote sensing, spatial analysis and GIS | |
Tugores et al. | Habitat suitability modelling for sardine Sardina pilchardus in a highly diverse ecosystem: the Mediterranean Sea | |
CN104834966A (zh) | 一种基于蚁群分类算法的渔场预报方法 | |
Visser et al. | A computer model to forecast wetland vegetation changes resulting from restoration and protection in coastal Louisiana | |
CN106651616B (zh) | 一种中西太平洋鲣鱼围网渔业入渔的预测方法 | |
Conquest et al. | Sampling methods for stock assessment for small-scale fisheries in developing countries | |
CN110929939B (zh) | 一种基于聚类-信息量耦合模型下的滑坡灾害易发性空间预测方法 | |
CN106339416A (zh) | 基于网格快速搜寻密度峰值的数据聚类方法 | |
CN104914467A (zh) | 提取分类模型道的地震相聚类分析方法 | |
CN104809479A (zh) | 基于支持向量机的鱼类栖息地适宜性指数建模方法 | |
Hua et al. | Habitat suitability of Pacific saury (Cololabis saira) based on a yield-density model and weighted analysis | |
CN108665163A (zh) | 服务于主体功能区的生态功能综合评估方法 | |
CN111882244A (zh) | 基于层次框架的多源国土开发风险评估体系的构建方法 | |
Ayhan et al. | Analysis of image classification methods for remote sensing | |
CN108388683A (zh) | 一种基于生境因子的植被类型空间模拟方法 | |
CN117290684B (zh) | 基于Transformer的高温干旱天气预警方法、电子设备 | |
Chen et al. | The application of the genetic adaptive neural network in landslide disaster assessment | |
CN117787081A (zh) | 一种基于Morris和Sobol法的水文模型参数不确定性分析方法 | |
Noor et al. | Prediction map of rainfall classification using random forest and inverse distance weighted (IDW) | |
CN106204314A (zh) | 一种东南太平洋茎柔鱼中心渔场预测方法 | |
CN114997666A (zh) | 一种区域泥石流易发性评估方法 | |
Li et al. | Evaluating the relationship between spatial heterogeneity and temporal variability of larval fish assemblages in a coastal marine ecosystem (Haizhou Bay, China) | |
Amad et al. | Mapping the Mangrove Vulnerability Index Using Geographical Information System | |
Hennen et al. | Global-Detector-GIS-and Knowledge-based Tool for a Global Detection of the Potential for Production, Supply and Demand. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180724 |