CN105526874A - 一种基于光谱特征参数的油膜厚度识别方法 - Google Patents
一种基于光谱特征参数的油膜厚度识别方法 Download PDFInfo
- Publication number
- CN105526874A CN105526874A CN201510882371.1A CN201510882371A CN105526874A CN 105526874 A CN105526874 A CN 105526874A CN 201510882371 A CN201510882371 A CN 201510882371A CN 105526874 A CN105526874 A CN 105526874A
- Authority
- CN
- China
- Prior art keywords
- spectral
- oil film
- film thickness
- characteristic
- characteristic parameters
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于光谱特征参数的油膜厚度识别方法,整理室外模拟实验测量的纯油种的ASD光谱数据库,进行油膜光谱的分析工作,再经过高光谱遥感图像数据预处理、光谱重建,利用ASD光谱数据库对重建的光谱特征进行评价,识别特征谱段及特征参数,建立波谱特征拟合模型。本发明考虑实际应用时大气透过率及传感器增益产生的背景噪声,优选可消除影响的特征参数,建立特征参数与厚度的查找表,解决油膜厚度反演的不稳定和精度不高的问题。
Description
技术领域
本发明涉及油膜测量技术领域,具体是一种基于光谱特征参数的油膜厚度识别方法。
背景技术
水上溢油厚度识别属目前待解决的难题。现有研究主要利用空基/天基传感器的光谱数值与油膜厚度建立查找表,如Jan等(2008)利用航空平台搭载四通道传感器DMSC-MKII,探讨了基于紫外-可见光-近红外(UV-visible-nearIR,下同)的油膜厚度识别算法,确定了厚度可探测范围和厚度查找表;OceanImaging公司(2010)使用机载可见光-近红外多光谱设备实现了应急下的相对厚度识别;刘丙新等(2014)基于地面实测光谱特征建立了决策树分类方法,用于机载多光谱遥感数据的油膜相对厚度信息提取。而高光谱为准确、有效区分各种海面目标提供了新的技术思路,Clark等(2010)基于地面实测光谱与AVIRIS航空数据反演了BP事故油膜厚度分布,估算了溢油量,指出基于航空、卫星平台一次性获取大型事故油膜厚度分布的技术可行性,但工作量巨大。孙鹏等(2013)从地面实测光谱与油膜厚度响应特性出发,找出油膜厚度光谱特征参数,为多光谱影像油膜厚度的定量评估提供了方法基础。
现有技术一的技术方案:基于空基实验建立光谱通道与厚度的查找表。缺点:(1)缺乏地面实验数据的支持,只适用于一种环境条件,结果不可靠。(2)光谱通道数据难以消除传感器噪声和大气透过率的影响。
现有技术二的技术方案:基于地面光谱实验建立厚度特征参数与油膜厚度的查找表。缺点:未考虑传感器噪声和大气透过率的影响,难于应用到空基/天基传感器,而空基/天基传感器是应急应用的关键。
上述研究很难投入实际应用,其主要原因在于油种特性与传感器光学特性对厚度反演信号的影响。
发明内容
本发明的目的在于提供一种油膜厚度反演的稳定和精度提高的基于光谱特征参数的油膜厚度识别方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
一种基于光谱特征参数的油膜厚度识别方法,步骤如下:
1)整理室外模拟实验测量的纯油种的ASD光谱数据库:ASD光谱数据库包含不同纯油种、不同油水比例的油膜厚度时光谱反射率;整理同步获取的高光谱遥感图像数据;
2)进行油膜光谱的分析工作:首先分析不同纯油种成分及分子结构变化的属性光谱表征,找出不同厚度尺度下C-H特征的变化规律、获取不同特征谱段对油膜厚度的敏感度及变化规律,进一步定量评价厚度识别范围和反演精度,并用于评价光谱重建后的光谱;
3)高光谱遥感图像数据预处理:对高光谱遥感图像数据进行几何校正与辐射校正,并对校正后的数据进行定量化评估;
4)光谱重建:光谱重建通过同步获取遥感影像和地面对应的ASD地物光谱仪同步实测光谱以及地面辐射控制点,建立线性方程,实现辐射校正对获取的高光谱遥感图像数据进行光谱重建;
5)利用步骤1)的室外模拟实验的ASD光谱数据库对重建的光谱特征进行评价:利用相似系数法对重建的光谱特征进行精度评价;并以光谱重建评价为基础,进一步对光谱重建方案进行优化;
6)识别特征谱段及特征参数:依据油膜厚度-光谱响应规律,有效划分并优选不同厚度尺度下厚度识别的特征谱段及特征参数,其中特征谱段包括可见光近红外谱段、中红外谱段、远红外谱段内的特征波谱曲线;特征参数包括Db、绿峰反射率Rg、Dy、Rg/Ro、(Rg-Ro)/(Rg+Ro);
7)建立波谱特征拟合模型:将特征参数与油膜厚度进行线性相关分析;所建波谱特征拟合模型如下:
Rg/Ro的拟合函数为4.02168*10-4+0.00168x;
(Rg-Ro)/(Rg+Ro)的拟合函数为0.39608+0.72537x。
作为本发明进一步的方案:辐射校正是指大气校正。
作为本发明进一步的方案:步骤4)中利用相似系数法对光谱重建后的数据进行精度评价的方法,包括以下步骤:
设实测时反射率光谱数据为
Xi=(xi1,…,xim),i=1,…,n,于是ASD光谱数据库数据矩阵L表示为
(式1)
设高光谱遥感图像数据中待评价像元a的重建光谱矢量为a=(a1,…,am),选择相似系数评价函数,相似系数值大者判断为优;相似系数计算见式2:
其中,
作为本发明进一步的方案:步骤7)中油膜的厚度上限为300μm。
与现有技术相比,本发明的有益效果是:
本发明以同步开展地-空实验,考虑实际应用时大气透过率及传感器增益产生的背景噪声,优选可消除影响的特征参数,建立特征参数与厚度的查找表,解决油膜厚度反演的不稳定和精度不高的问题。
附图说明
图1是本发明实施例1的流程图。
具体实施方式
下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
请参阅图1,本发明实施例中,本发明基于光谱特征参数建立油膜厚度识别方案,关键在于解决海洋环境下的溢油目标光谱重建问题和特征参数优选问题,设置地面-空基同步实验实现光谱重建,获取空基视反射光谱曲线-厚度关系,进一步提取特征参数,建立波谱特征拟合模型以优选特征参数。具体步骤如下:
1)整理室外模拟实验测量的纯油种的ASD光谱数据库:ASD光谱数据库包含不同纯油种、不同油水比例的油膜厚度时光谱反射率;整理同步获取的高光谱遥感图像数据;
2)进行油膜光谱的分析工作:首先分析不同纯油种成分及分子结构变化的属性光谱表征,找出不同厚度尺度下C-H特征的变化规律、获取不同谱段对油膜厚度的敏感度及变化规律,进一步定量评价厚度识别范围和反演精度,并用于评价光谱重建后的光谱;
3)高光谱遥感图像数据预处理:对高光谱遥感图像数据进行几何校正与辐射校正,并对校正后的数据进行定量化评估;其中辐射校正是指大气校正,针对同步获取的高光谱遥感图像数据采用建立的光谱重建方案进行大气校正;
4)光谱重建:光谱重建通过同步获取遥感影像和地面对应的ASD地物光谱仪同步实测光谱以及地面辐射控制点,建立线性方程,实现大气校正对获取的高光谱遥感图像数据进行光谱重建,真实反映不同油膜厚度的光谱特征;
5)利用步骤1)的室外模拟实验的ASD光谱数据库对重建的光谱特征进行评价:利用相似系数法进行精度评价;并以光谱重建评价为基础,进一步对重建方案进行优化;
其中利用相似系数法对光谱重建后的数据进行精度评价的方法,包括以下步骤:
设实测时反射率光谱数据为
Xi=(xi1,…,xim),i=1,…,n,
于是ASD光谱数据库数据矩阵L表示为
(式1)
设高光谱遥感图像数据中待评价像元a的重建光谱矢量为a=(a1,…,am),选择相似系数评价函数,相似系数值大者判断为优;相似系数计算见式2:
其中,
6)识别特征谱段及特征参数:依据油膜厚度-光谱响应规律,有效划分并优选不同厚度尺度下厚度识别特征谱段及特征参数,其中特征谱段包括可见光近红外谱段、中红外谱段、远红外谱段内的特征波谱曲线;特征参数包括Db、绿峰反射率Rg、Dy、Rg/Ro、(Rg-Ro)/(Rg+Ro);
7)建立波谱特征拟合模型:光学传感器探测原油厚度上限是200μm~300μm,实际应用中,光学传感器将只能区分厚度为300μm的油膜;故将特征参数与厚度上限300μm油膜进行线性相关分析;最终所建波谱特征拟合模型如表1:
表1
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (5)
1.一种基于光谱特征参数的油膜厚度识别方法,其特征在于,步骤如下:
1)整理室外模拟实验测量的纯油种的ASD光谱数据库:ASD光谱数据库包含不同纯油种、不同油水比例的油膜厚度时光谱反射率;整理同步获取的高光谱遥感图像数据;
2)进行油膜光谱的分析工作:首先分析不同纯油种成分及分子结构变化的属性光谱表征,找出不同厚度尺度下C-H特征的变化规律、获取不同特征谱段对油膜厚度的敏感度及变化规律,进一步定量评价厚度识别范围和反演精度,并用于评价光谱重建后的光谱;
3)高光谱遥感图像数据预处理:对高光谱遥感图像数据进行几何校正与辐射校正,并对校正后的数据进行定量化评估;
4)光谱重建:光谱重建通过同步获取遥感影像和地面对应的ASD地物光谱仪同步实测光谱以及地面辐射控制点,建立线性方程,实现辐射校正对获取的高光谱遥感图像数据进行光谱重建;
5)利用步骤1)的室外模拟实验的ASD光谱数据库对重建的光谱特征进行评价:利用相似系数法对重建的光谱特征进行精度评价;并以光谱重建评价为基础,进一步对光谱重建方案进行优化;
6)识别特征谱段及特征参数:依据油膜厚度-光谱响应规律,有效划分并优选不同厚度尺度下厚度识别的特征谱段及特征参数,其中特征谱段包括可见光近红外谱段、中红外谱段、远红外谱段内的特征波谱曲线;特征参数包括Db、绿峰反射率Rg、Dy、Rg/Ro、(Rg-Ro)/(Rg+Ro);
7)建立波谱特征拟合模型:将特征参数与油膜厚度进行线性相关分析;所建波谱特征拟合模型如下:
Rg/Ro的拟合函数为4.02168*10-4+0.00168x;
(Rg-Ro)/(Rg+Ro)的拟合函数为0.39608+0.72537x。
2.根据权利要求1所述的基于光谱特征参数的油膜厚度识别方法,其特征在于,辐射校正是指大气校正。
3.根据权利要求1所述的基于光谱特征参数的油膜厚度识别方法,其特征在于,步骤4)中利用相似系数法对光谱重建后的数据进行精度评价的方法,包括以下步骤:
设实测时反射率光谱数据为
Xi=(xi1,…,xim),i=1,…,n,
于是ASD光谱数据库数据矩阵L表示为
设高光谱遥感图像数据中待评价像元a的重建光谱矢量为a=(a1,…,am),选择相似系数评价函数,相似系数值大者判断为优;相似系数计算见式2:
其中,
4.根据权利要求1所述的基于光谱特征参数的油膜厚度识别方法,其特征在于,步骤7)所建识别模型,采用可消除大气透过率及传感器增益产生的背景噪声的比值特征参数。
5.根据权利要求1所述的基于光谱特征参数的油膜厚度识别方法,其特征在于,步骤7)中油膜的厚度上限为300μm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510882371.1A CN105526874B (zh) | 2015-12-03 | 2015-12-03 | 一种基于光谱特征参数的油膜厚度识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510882371.1A CN105526874B (zh) | 2015-12-03 | 2015-12-03 | 一种基于光谱特征参数的油膜厚度识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105526874A true CN105526874A (zh) | 2016-04-27 |
CN105526874B CN105526874B (zh) | 2019-09-20 |
Family
ID=55769251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510882371.1A Expired - Fee Related CN105526874B (zh) | 2015-12-03 | 2015-12-03 | 一种基于光谱特征参数的油膜厚度识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105526874B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106767457A (zh) * | 2016-12-19 | 2017-05-31 | 中国科学院烟台海岸带研究所 | 一种基于拉曼光谱测量的水面油膜厚度测量方法及装置 |
CN108956538A (zh) * | 2018-06-28 | 2018-12-07 | 中国石油天然气股份有限公司 | 河流溢油量遥感检测方法及装置 |
CN109509209A (zh) * | 2018-11-15 | 2019-03-22 | 上海卫星工程研究所 | 利用高光谱技术探测海空环境空中动目标的分析方法 |
CN111121653A (zh) * | 2019-12-31 | 2020-05-08 | 武汉颐光科技有限公司 | 一种单层薄膜临界厚度估值计算方法 |
CN111291665A (zh) * | 2020-01-22 | 2020-06-16 | 广东牵引信息科技有限公司 | 一种海洋溢油乳化预测方法、设备及存储介质 |
CN111406198A (zh) * | 2020-02-24 | 2020-07-10 | 长江存储科技有限责任公司 | 用于半导体芯片表面形貌计量的系统和方法 |
CN111597751A (zh) * | 2020-03-24 | 2020-08-28 | 自然资源部第一海洋研究所 | 基于自扩展深度置信网络的原油油膜绝对厚度反演方法 |
CN113267132A (zh) * | 2021-07-16 | 2021-08-17 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于双色成像的全局油膜厚度在线测量方法 |
CN113424017A (zh) * | 2019-03-15 | 2021-09-21 | 陶氏环球技术有限责任公司 | 通过近红外高光谱成像进行测量的薄膜厚度测量仪 |
US11243067B2 (en) | 2020-02-24 | 2022-02-08 | Yangtze Memory Technologies Co., Ltd. | Systems and methods for semiconductor chip surface topography metrology |
US11448499B2 (en) | 2020-02-24 | 2022-09-20 | Yangtze Memory Technologies Co., Ltd. | Systems and methods for semiconductor chip surface topography metrology |
US11454491B2 (en) | 2020-02-24 | 2022-09-27 | Yangtze Memory Technologies Co., Ltd. | Systems having light source with extended spectrum for semiconductor chip surface topography metrology |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070146725A1 (en) * | 2005-11-10 | 2007-06-28 | Belov Michael L | Method of and device for thickness measurements of thin films |
CN101532825A (zh) * | 2009-04-27 | 2009-09-16 | 天津大学 | 基于差分激光三角法的海面溢油油膜厚度测量方法 |
CN102507474A (zh) * | 2011-10-28 | 2012-06-20 | 大连海事大学 | 一种船舶溢油目标的识别方法及系统 |
-
2015
- 2015-12-03 CN CN201510882371.1A patent/CN105526874B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070146725A1 (en) * | 2005-11-10 | 2007-06-28 | Belov Michael L | Method of and device for thickness measurements of thin films |
CN101532825A (zh) * | 2009-04-27 | 2009-09-16 | 天津大学 | 基于差分激光三角法的海面溢油油膜厚度测量方法 |
CN102507474A (zh) * | 2011-10-28 | 2012-06-20 | 大连海事大学 | 一种船舶溢油目标的识别方法及系统 |
Non-Patent Citations (2)
Title |
---|
兰国新等: "海上溢油遥感光谱信息挖掘与应用研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
刘秀英等: "《森林树种高光谱遥感研究》", 31 December 2011 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106767457A (zh) * | 2016-12-19 | 2017-05-31 | 中国科学院烟台海岸带研究所 | 一种基于拉曼光谱测量的水面油膜厚度测量方法及装置 |
CN108956538B (zh) * | 2018-06-28 | 2021-03-09 | 中国石油天然气股份有限公司 | 河流溢油量遥感检测方法及装置 |
CN108956538A (zh) * | 2018-06-28 | 2018-12-07 | 中国石油天然气股份有限公司 | 河流溢油量遥感检测方法及装置 |
CN109509209A (zh) * | 2018-11-15 | 2019-03-22 | 上海卫星工程研究所 | 利用高光谱技术探测海空环境空中动目标的分析方法 |
CN109509209B (zh) * | 2018-11-15 | 2023-08-15 | 上海卫星工程研究所 | 利用高光谱技术探测海空环境空中动目标的分析方法 |
CN113424017A (zh) * | 2019-03-15 | 2021-09-21 | 陶氏环球技术有限责任公司 | 通过近红外高光谱成像进行测量的薄膜厚度测量仪 |
CN111121653A (zh) * | 2019-12-31 | 2020-05-08 | 武汉颐光科技有限公司 | 一种单层薄膜临界厚度估值计算方法 |
CN111291665A (zh) * | 2020-01-22 | 2020-06-16 | 广东牵引信息科技有限公司 | 一种海洋溢油乳化预测方法、设备及存储介质 |
CN111291665B (zh) * | 2020-01-22 | 2022-04-01 | 广东牵引信息科技有限公司 | 一种海洋溢油乳化预测方法、设备及存储介质 |
CN111406198B (zh) * | 2020-02-24 | 2021-02-19 | 长江存储科技有限责任公司 | 用于半导体芯片表面形貌计量的系统和方法 |
US11243067B2 (en) | 2020-02-24 | 2022-02-08 | Yangtze Memory Technologies Co., Ltd. | Systems and methods for semiconductor chip surface topography metrology |
US11448499B2 (en) | 2020-02-24 | 2022-09-20 | Yangtze Memory Technologies Co., Ltd. | Systems and methods for semiconductor chip surface topography metrology |
US11454491B2 (en) | 2020-02-24 | 2022-09-27 | Yangtze Memory Technologies Co., Ltd. | Systems having light source with extended spectrum for semiconductor chip surface topography metrology |
US11562919B2 (en) | 2020-02-24 | 2023-01-24 | Yangtze Memory Technologies Co., Ltd. | Systems and methods for semiconductor chip surface topography metrology |
CN111406198A (zh) * | 2020-02-24 | 2020-07-10 | 长江存储科技有限责任公司 | 用于半导体芯片表面形貌计量的系统和方法 |
US11796307B2 (en) | 2020-02-24 | 2023-10-24 | Yangtze Memory Technologies Co., Ltd. | Systems and methods for semiconductor chip surface topography metrology |
CN111597751A (zh) * | 2020-03-24 | 2020-08-28 | 自然资源部第一海洋研究所 | 基于自扩展深度置信网络的原油油膜绝对厚度反演方法 |
CN111597751B (zh) * | 2020-03-24 | 2023-10-24 | 自然资源部第一海洋研究所 | 基于自扩展深度置信网络的原油油膜绝对厚度反演方法 |
CN113267132A (zh) * | 2021-07-16 | 2021-08-17 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于双色成像的全局油膜厚度在线测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105526874B (zh) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105526874B (zh) | 一种基于光谱特征参数的油膜厚度识别方法 | |
Ma et al. | Coupling urbanization analyses for studying urban thermal environment and its interplay with biophysical parameters based on TM/ETM+ imagery | |
Li et al. | MODIS observations of water color of the largest 10 lakes in China between 2000 and 2012 | |
Galvao et al. | On intra-annual EVI variability in the dry season of tropical forest: A case study with MODIS and hyperspectral data | |
CN102708307B (zh) | 一种应用于城市的植被指数构造方法 | |
CN103456122B (zh) | 一种森林火灾烟雾识别方法及装置 | |
Hill et al. | Characterizing vegetation cover in global savannas with an annual foliage clumping index derived from the MODIS BRDF product | |
CN102540165B (zh) | Modis地表反射率数据的预处理方法及系统 | |
CN107688003B (zh) | 一种消除植被冠层结构和地表背景影响的叶片反射率卫星遥感提取方法 | |
CN110100262B (zh) | 用于从图像移除云的图像处理设备、方法和存储介质 | |
US20050114027A1 (en) | Cloud shadow detection: VNIR-SWIR | |
CN106845326B (zh) | 一种基于航空高光谱遥感数据的冰川识别方法 | |
CN107368617B (zh) | 基于Lowtran7大气软件的地空探测红外成像系统作用距离计算方法 | |
Zhai | Inversion of organic matter content in wetland soil based on Landsat 8 remote sensing image | |
CN109766909B (zh) | 基于谱图融合的海岸环境微塑料老化行为解析方法 | |
CN104820250A (zh) | 一种极轨气象卫星virr海洋上云检测的处理方法 | |
CN107421892A (zh) | 一种用于水体信息提取的高光谱数据处理方法 | |
CN116977873B (zh) | 一种基于高分卫星遥感的裸露沙地决策树提取方法 | |
CN107991243B (zh) | 一种基于航空高光谱遥感数据的高海拔地区水体提取方法 | |
Lee et al. | New approach for snow cover detection through spectral pattern recognition with MODIS data | |
CN105718924B (zh) | 基于多特征综合及机器学习的高分卫星影像云检测方法 | |
Wu et al. | A comparison of 4 shadow compensation techniques for land cover classification of shaded areas from high radiometric resolution aerial images | |
CN107657206B (zh) | 一种基于遥感技术估算绿地覆盖率的方法 | |
CN110009584B (zh) | 基于参考光谱匹配的多光谱遥感影像大气校正系统及方法 | |
Moreira et al. | Variation in spectral shape of urban materials |
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: 20190920 Termination date: 20201203 |
|
CF01 | Termination of patent right due to non-payment of annual fee |