CN106908218A - 河道表面流场的处理方法 - Google Patents

河道表面流场的处理方法 Download PDF

Info

Publication number
CN106908218A
CN106908218A CN201710179000.6A CN201710179000A CN106908218A CN 106908218 A CN106908218 A CN 106908218A CN 201710179000 A CN201710179000 A CN 201710179000A CN 106908218 A CN106908218 A CN 106908218A
Authority
CN
China
Prior art keywords
camera
flow field
surface flow
world coordinate
image
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
CN201710179000.6A
Other languages
English (en)
Other versions
CN106908218B (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.)
BEIJING SINFOTEK TECHNOLOGY Co Ltd
Original Assignee
BEIJING SINFOTEK TECHNOLOGY Co Ltd
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 BEIJING SINFOTEK TECHNOLOGY Co Ltd filed Critical BEIJING SINFOTEK TECHNOLOGY Co Ltd
Priority to CN201710179000.6A priority Critical patent/CN106908218B/zh
Publication of CN106908218A publication Critical patent/CN106908218A/zh
Application granted granted Critical
Publication of CN106908218B publication Critical patent/CN106908218B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Multimedia (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Closed-Circuit Television Systems (AREA)

Abstract

本发明河道表面流场的处理方法，它包括：在待测量区域的表面流场内撒若干个示踪粒子，相机对准上述示踪粒子的目标区域进行连续拍摄，标注出各个图像的序列，通过所拍摄的图像坐标，分别计算出相机所拍摄图像中每个点所对应的新世界坐标下的坐标(X”,Y”)，得到相机所拍摄的校正图像；根据上述相机校正的图像，在相邻时间的两个校正图像中，以最大相似的原则，对每个示踪粒子进行匹配，根据相邻时间的两个校正图像中所有粒子的匹配结果，得到每个示踪粒子的空间位移，得到所有相邻两个校正图像中所有示踪粒子的空间位移，从而获得表面流场，本发明河道表面流场的处理方法适用于粒子浓度不高的时候，解决了现有额粒子跟踪测速测量误差较大的问题，它可以进行大范围表面流场测量。

Description

河道表面流场的处理方法
技术领域
本发明涉及了一种计算流场的方法，特别是用于流场中目标粒子浓度比较少的情况下，进行目标粒子跟踪,从而获取大范围表面流场的测量方法。
背景技术
本发明涉及了一种计算流场的方法，特别是用于流场中目标粒子浓度比较少的情况下，进行目标粒子跟踪,从而获取大范围表面流场的测量方法。
发明内容
本申请的发明目的在于提供一种流场中目标粒子浓度比较少的情况下，大范围表面流场测量方法，利用目标粒子提取，并进行匹配跟踪的方法，获取大范围表面流场。
为了完成本申请的发明目的，本申请采用以下技术方案：
本发明的表面流场的计算方法，其中：它包括以下步骤：
1)在待测量区域的表面流场内，选取至少3个固定点,用测量设备获得上述至少3个固定点的世界坐标；
2)安置一台相机，使相机视野范围覆盖待测量区域的表面流场，相机对准待测量区域的表面流场进行拍摄，提取在世界坐标系下的上述至少3个固定点的坐标，和相机所拍摄图像中的上述至少3个固定点的图像坐标，
3)选出所有固定点的世界坐标X，Y的最小值Xmin，Ymin，以（Xmin，Ymin）点为新世界坐标系的原点位置，新世界坐标系X”、Y”轴方向与世界坐标系X、Y轴方向一致，建立新世界坐标系，并且在此基础上，对世界坐标系进行缩小倍数B，则可以得到新世界坐标系下的所有固定点坐标为
X”=(X-Xmin)/B
Y”=(Y-Ymin)/B；
4)根据相机固定点在新世界坐标系下的坐标和相机固定点在拍摄图像上的坐标，构建出相机固定点的新世界坐标和固定点在拍摄图像坐标的映射关系；固定点的新世界坐标为： …，固定点的拍摄图像坐标为： …
建立下列关系式：
计算得出式中的K值，
以新世界坐标系中的至少3个固定点为基点，根据下列公式，分别计算出相机拍摄图像中每个点所对应的新世界坐标系下的坐标(X”,Y”)
即通过计算得出X”,Y”坐标；
计算在世界坐标下任意两个固定点之间的距离L1，和该相机的在新世界坐标系下的上述两个固定点之间的距离L1”,得到像素的空间分辨率= L1/ L1”；
5)在待测量区域的表面流场内撒若干个示踪粒子，使示踪粒子漂浮在水面，每个示踪粒子在图像中成像的大小为2*2至3*3个像素，并且在5*5像素的区域范围内只有1个示踪粒子，上述示踪粒子为相机后续连续拍摄的目标区域；
6)在待测量区域的表面流场内，相机对准上述示踪粒子的目标区域进行连续拍摄，标注出各个图像的序列，通过所拍摄的图像坐标，根据下列公式，分别计算出相机所拍摄图像中每个点所对应的新世界坐标下的坐标(X”,Y”)
即通过计算得出新世界坐标系下的X”,Y”坐标，得到相机所拍摄的校正图像；
7)根据上述相机校正的图像，在相邻时间的两个校正图像中，以最大相似的原则，对每个示踪粒子进行匹配，根据相邻时间的两个校正图像中所有粒子的匹配结果，得到每个示踪粒子的像素的二维位移；
8)根据得到的像素空间分辨率，每个示踪粒子的空间速度=每个示踪粒子的像素的二维位移*像素的空间分辨率/两帧图像的时间间隔；
9)重复步骤7）和步骤8），得到所有相邻两个校正图像中所有示踪粒子的空间速度，从而获得表面流场。
本发明的表面流场的计算方法，其中：所述固定点为3-8个点；
本发明的表面流场的计算方法，其中：所述的B为>0的数。
本发明的表面流场的计算方法适用于粒子浓度不高的时候，解决了现有额粒子跟踪测速测量误差较大的问题，它可以进行大范围表面流场测量。
附图说明
图1为测量区域和相机示意图；
图2为相机所拍摄的固定点的示意图；
图3为由图2经过校正后的在新世界坐标系下的示意图；
图4是相机所拍摄的示踪粒子图像，经转换后，成为新世界坐标系下的相邻两个图像的叠加；
图5为图4匹配后的状态。在图4和图5中黑色示踪粒子为相邻图像中的前一桢的示踪粒子，灰色示踪粒子为相邻图像中的后一桢的示踪粒子
在图1至图4中，标号1为相机；标号2为固定点；标号3为示踪粒子。
具体实施方式
本发明的一种表面流场的计算方法包括以下步骤：
1)如图1所示，在待测量区域的表面流场内，选取4个固定点2,用测量设备获得4个固定点的世界坐标，图1至图3中，将固定点2画得很大是为了清楚起见，实际上只是一个坐标点；
2)安置一台相机1，使相机1视野范围覆盖待测量区域的表面流场，相机1对准待测量区域的表面流场进行拍摄，提取在世界坐标系下的4个固定点2的坐标，和图2所示的相机所拍摄图像中的4个固定点2的图像坐标，
3)选出所有固定点2的世界坐标X，Y的最小值Xmin，Ymin，以（Xmin，Ymin）点为新世界坐标系的原点位置，新世界坐标系X”、Y”轴方向与世界坐标系X、Y轴方向一致，建立新世界坐标系，并且在此基础上，对世界坐标系进行缩小倍数B，则可以得到新世界坐标系下的所有固定点坐标为
X”=(X-Xmin)/B
Y”=(Y-Ymin)/B，B为>0的数；
4)根据相机固定点2在新世界坐标系下的坐标和相机固定点在拍摄图像上的坐标，构建出相机固定点的新世界坐标和固定点在拍摄图像坐标的映射关系；固定点的新世界坐标为： …，固定点的拍摄图像坐标为： …
建立下列关系式：
计算得出式中的K值，
以新世界坐标系中的4个固定点为基点，根据下列公式，计算出相机拍摄图像中每个点所对应的新世界坐标系下的坐标(X”,Y”)
即通过计算得出X”,Y”坐标，得到图3,；
计算在世界坐标下任意两个固定点之间的距离L1，和该相机的在新世界坐标系下的上述两个固定点之间的距离L1”,得到像素的空间分辨率= L1/ L1”；
5)在待测量区域的表面流场内撒若干个示踪粒子3，使示踪粒子3漂浮在水面，每个示踪粒子在图像中成像的大小为2*2至3*3个像素，并且在5*5像素的区域范围内只有1个示踪粒子3，上述示踪粒子3为相机1后续连续拍摄的目标区域；
6)在待测量区域的表面流场内，相机1对准上述示踪粒子3的目标区域进行连续拍摄，标注出各个图像的序列，通过所拍摄的图像坐标，根据下列公式，分别计算出相机所拍摄图像中每个点所对应的新世界坐标下的坐标(X”,Y”)
即通过计算得出新世界坐标系下的X”,Y”坐标，得到相机所拍摄的校正图像；
7)根据上述相机校正的图像，如图4所示，在相邻时间的两个校正图像中，以最大相似的原则，对每个示踪粒子进行匹配，根据相邻时间的两个校正图像中所有粒子的匹配结果，得到每个示踪粒子的像素的二维位移，图5所示；
8)根据得到的像素空间分辨率，每个示踪粒子的空间速度=每个示踪粒子的像素的二维位移*像素的空间分辨率/两帧图像的时间间隔；
9)重复步骤7）和步骤8），得到所有相邻两个校正图像中所有示踪粒子的空间速度，从而获得表面流场。
以上描述是对本发明的解释，不是对发明的限定，本发明所限定的范围参见权利要求，在不违背本发明的精神的情况下，本发明可以作任何形式的修改。

Claims (3)

1.一种河道表面流场的处理方法，其特征在于：它包括以下步骤：
在待测量区域的表面流场内，选取至少3个固定点,用测量设备获得上述至少3个固定点的世界坐标；
安置一台相机，使相机视野范围覆盖待测量区域的表面流场，相机对准待测量区域的表面流场进行拍摄，提取在世界坐标系下的上述至少3个固定点的坐标，和相机所拍摄图像中的上述至少3个固定点的图像坐标，
选出所有固定点的世界坐标X，Y的最小值Xmin，Ymin，以（Xmin，Ymin）点为新世界坐标系的原点位置，新世界坐标系X”、Y”轴方向与世界坐标系X、Y轴方向一致，建立新世界坐标系，并且在此基础上，对世界坐标系进行缩小倍数B，则可以得到新世界坐标系下的所有固定点坐标为
X”=(X-Xmin)/B
Y”=(Y-Ymin)/B
根据相机固定点在新世界坐标系下的坐标和相机固定点在拍摄图像上的坐标，分别构建出相机固定点的新世界坐标和固定点在拍摄图像坐标的映射关系；固定点的新世界坐标为： …，固定点的拍摄图像坐标为：…
建立下列关系式：
计算得出式中的K值，
以新世界坐标系中的至少3个固定点为基点，根据下列公式，分别计算出相机拍摄图像中每个点所对应的新世界坐标系下的坐标(X”,Y”)
即通过计算得出X”,Y”坐标；
计算在世界坐标下任意两个固定点之间的距离L1，和该相机的在新世界坐标系下的上述两个固定点之间的距离L1”,得到像素的空间分辨率= L1/ L1”；
在待测量区域的表面流场内撒若干个示踪粒子，使示踪粒子漂浮在水面，每个示踪粒子在图像中成像的大小为2*2至3*3个像素，并且在5*5像素的区域范围内只有1个示踪粒子，上述示踪粒子为相机后续连续拍摄的目标区域；
在待测量区域的表面流场内，相机对准上述示踪粒子的目标区域进行连续拍摄，标注出各个图像的序列，通过所拍摄的图像坐标，根据下列公式，分别计算出相机所拍摄图像中每个点所对应的新世界坐标下的坐标(X”,Y”)
即通过计算得出新世界坐标系下的X”,Y”坐标，得到相机所拍摄的校正图像；
根据上述相机校正的图像，在相邻时间的两个校正图像中，以最大相似的原则，对每个示踪粒子进行匹配，根据相邻时间的两个校正图像中所有粒子的匹配结果，得到每个示踪粒子的像素的二维位移；
根据上述的像素空间分辨率，每个示踪粒子的空间速度=每个示踪粒子的像素的二维位移*像素的空间分辨率/两帧图像的时间间隔；
重复步骤7）和步骤8），得到所有相邻两个校正图像中所有示踪粒子的空间速度，从而获得表面流场。
2.如权利要求1所述的表面流场的计算方法，其特征在于：所述固定点为3-8个点。
3.如权利要求1所述的表面流场的计算方法，其特征在于：所述的B为>0的数。
CN201710179000.6A 2017-03-23 2017-03-23 河道表面流场的处理方法 Expired - Fee Related CN106908218B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710179000.6A CN106908218B (zh) 2017-03-23 2017-03-23 河道表面流场的处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710179000.6A CN106908218B (zh) 2017-03-23 2017-03-23 河道表面流场的处理方法

Publications (2)

Publication Number Publication Date
CN106908218A true CN106908218A (zh) 2017-06-30
CN106908218B CN106908218B (zh) 2018-11-30

Family

ID=59194697

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710179000.6A Expired - Fee Related CN106908218B (zh) 2017-03-23 2017-03-23 河道表面流场的处理方法

Country Status (1)

Country Link
CN (1) CN106908218B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110187142A (zh) * 2019-06-13 2019-08-30 上海彩虹鱼海洋科技股份有限公司 流速监测方法和系统
CN112067839A (zh) * 2020-09-18 2020-12-11 重庆交通大学 无控制点条件下的无人机超宽河流表面流场快速测量方法
CN112698053A (zh) * 2020-12-02 2021-04-23 武汉理工大学 粒子图像测速误差优化方法、系统和存储介质
CN114885097A (zh) * 2022-03-31 2022-08-09 浙江大华技术股份有限公司 一种智能相机监测水流的方法、装置及存储介质

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2485725A1 (fr) * 1980-06-28 1981-12-31 Kv Assoc Inc Procede et dispositif de mesure de direction et de debit de fluide, notamment eaux souterraines, par determination d'un champs thermique
US5177607A (en) * 1990-08-27 1993-01-05 Zexel Corporation Method and apparatus for measuring velocity of fluid
CN1116704A (zh) * 1994-12-23 1996-02-14 清华大学 河工模型试验流场量测中的图像摄取系统及图象处理方法
KR100345918B1 (ko) * 1999-03-25 2002-07-27 학교법인조선대학교 유동장의 화상처리 방법
CN102313684A (zh) * 2010-07-08 2012-01-11 中国科学院过程工程研究所 气固两相流流场实时测量系统及方法
CN102564508A (zh) * 2011-12-14 2012-07-11 河海大学 基于视频图像的河流流量在线测验实现方法
CN102866260A (zh) * 2012-09-18 2013-01-09 河海大学 非接触式河流表面流场成像量测方法
CN103308714A (zh) * 2013-06-06 2013-09-18 水利部交通运输部国家能源局南京水利科学研究院 一种水流中彩色粒子图像测速系统及其测量方法
CN103336144A (zh) * 2013-05-27 2013-10-02 中国船舶重工集团公司第七○二研究所 一种水面微速流场的测试方法
CN103558409A (zh) * 2013-10-24 2014-02-05 水利部交通运输部国家能源局南京水利科学研究院 一种水流中分布式ptv流场测量系统及其测量方法
CN103645341A (zh) * 2013-11-19 2014-03-19 北京信息科技大学 全流场3d可视化测速方法
CN203504692U (zh) * 2013-10-25 2014-03-26 北京尚水信息技术股份有限公司 河道表面流速的测量装置
CN104215795A (zh) * 2014-09-29 2014-12-17 南京昊控软件技术有限公司 基于昊控表面流场计算算法的大尺度粒子图像测速方法
CN106023097A (zh) * 2016-05-11 2016-10-12 华中科技大学 一种基于迭代法的流场图像预处理算法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2485725A1 (fr) * 1980-06-28 1981-12-31 Kv Assoc Inc Procede et dispositif de mesure de direction et de debit de fluide, notamment eaux souterraines, par determination d'un champs thermique
US5177607A (en) * 1990-08-27 1993-01-05 Zexel Corporation Method and apparatus for measuring velocity of fluid
CN1116704A (zh) * 1994-12-23 1996-02-14 清华大学 河工模型试验流场量测中的图像摄取系统及图象处理方法
KR100345918B1 (ko) * 1999-03-25 2002-07-27 학교법인조선대학교 유동장의 화상처리 방법
CN102313684A (zh) * 2010-07-08 2012-01-11 中国科学院过程工程研究所 气固两相流流场实时测量系统及方法
CN102564508A (zh) * 2011-12-14 2012-07-11 河海大学 基于视频图像的河流流量在线测验实现方法
CN102866260A (zh) * 2012-09-18 2013-01-09 河海大学 非接触式河流表面流场成像量测方法
CN103336144A (zh) * 2013-05-27 2013-10-02 中国船舶重工集团公司第七○二研究所 一种水面微速流场的测试方法
CN103308714A (zh) * 2013-06-06 2013-09-18 水利部交通运输部国家能源局南京水利科学研究院 一种水流中彩色粒子图像测速系统及其测量方法
CN103558409A (zh) * 2013-10-24 2014-02-05 水利部交通运输部国家能源局南京水利科学研究院 一种水流中分布式ptv流场测量系统及其测量方法
CN203504692U (zh) * 2013-10-25 2014-03-26 北京尚水信息技术股份有限公司 河道表面流速的测量装置
CN103645341A (zh) * 2013-11-19 2014-03-19 北京信息科技大学 全流场3d可视化测速方法
CN104215795A (zh) * 2014-09-29 2014-12-17 南京昊控软件技术有限公司 基于昊控表面流场计算算法的大尺度粒子图像测速方法
CN106023097A (zh) * 2016-05-11 2016-10-12 华中科技大学 一种基于迭代法的流场图像预处理算法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
唐洪武等: "实体模型表面流场、河势测量中图像技术应用研究进展", 《河海大学学报（自然科学版）》 *
李蔚: "基于图像处理的山区河道表面流场测算研究", 《人民长江》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110187142A (zh) * 2019-06-13 2019-08-30 上海彩虹鱼海洋科技股份有限公司 流速监测方法和系统
CN112067839A (zh) * 2020-09-18 2020-12-11 重庆交通大学 无控制点条件下的无人机超宽河流表面流场快速测量方法
CN112698053A (zh) * 2020-12-02 2021-04-23 武汉理工大学 粒子图像测速误差优化方法、系统和存储介质
CN112698053B (zh) * 2020-12-02 2022-02-22 武汉理工大学 粒子图像测速误差优化方法、系统和存储介质
CN114885097A (zh) * 2022-03-31 2022-08-09 浙江大华技术股份有限公司 一种智能相机监测水流的方法、装置及存储介质
CN114885097B (zh) * 2022-03-31 2024-02-06 浙江大华技术股份有限公司 一种智能相机监测水流的方法、装置及存储介质

Also Published As

Publication number Publication date
CN106908218B (zh) 2018-11-30

Similar Documents

Publication Publication Date Title
CN106908218A (zh) 河道表面流场的处理方法
CN105447888B (zh) 一种基于有效目标判断的无人机机动目标检测方法
US8755624B2 (en) Image registration device and method thereof
KR102523451B1 (ko) 드론 영상 기반의 하천 유속 측정을 위한 장치 및 방법
CN107886547B (zh) 一种鱼眼相机标定方法及系统
CN101521740A (zh) 基于多尺度不变特征的实时运动估计方法
CN104361603B (zh) 枪机图像目标标定方法及其系统
CN112801870B (zh) 一种基于网格优化的图像拼接方法，拼接系统及可读存储介质
CN106033614B (zh) 一种强视差下的移动相机运动目标检测方法
CN107392951A (zh) 遥感图像高精度快速配准方法
CN109214254B (zh) 一种确定机器人位移的方法及装置
CN107966137A (zh) 一种基于tdiccd拼接区图像的卫星平台颤振探测方法
CN113012234B (zh) 基于平面变换的高精度相机标定方法
CN112270719A (zh) 相机标定方法、装置及系统
CN106152978B (zh) 一种基于图像分析的非接触性面积测量方法
Yang et al. Precise measurement of position and attitude based on convolutional neural network and visual correspondence relationship
US11410459B2 (en) Face detection and recognition method using light field camera system
CN115201883A (zh) 一种运动目标视频定位测速系统及方法
CN107146242A (zh) 一种用于成像系统模糊核估计的高精度图像配准方法
CN104937608B (zh) 道路区域检测
CN105869108B (zh) 一种用于动平台移动目标侦测中的图像配准方法
Feng et al. Mountainous remote sensing images registration based on improved optical flow estimation
CN107424194A (zh) 键盘轮廓度的检测方法
Feng et al. Registration of multitemporal GF-1 remote sensing images with weighting perspective transformation model
CN113670268B (zh) 基于双目视觉的无人机和电力杆塔距离测量方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Treatment method of river surface flow field

Effective date of registration: 20210715

Granted publication date: 20181130

Pledgee: Beijing first financing Company limited by guarantee

Pledgor: BEIJING SINFOTEK SCIENCE AND TECHNOLOGY Co.,Ltd.

Registration number: Y2021990000605

PE01 Entry into force of the registration of the contract for pledge of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181130

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