CN111898440A - 一种基于三维卷积神经网络的山火检测方法 - Google Patents

一种基于三维卷积神经网络的山火检测方法 Download PDF

Info

Publication number
CN111898440A
CN111898440A CN202010607252.6A CN202010607252A CN111898440A CN 111898440 A CN111898440 A CN 111898440A CN 202010607252 A CN202010607252 A CN 202010607252A CN 111898440 A CN111898440 A CN 111898440A
Authority
CN
China
Prior art keywords
mountain fire
neural network
videos
dimensional
size
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
CN202010607252.6A
Other languages
English (en)
Other versions
CN111898440B (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.)
Chengdu Sihan Technology Co ltd
Original Assignee
Chengdu Sihan 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 Chengdu Sihan Technology Co ltd filed Critical Chengdu Sihan Technology Co ltd
Priority to CN202010607252.6A priority Critical patent/CN111898440B/zh
Publication of CN111898440A publication Critical patent/CN111898440A/zh
Application granted granted Critical
Publication of CN111898440B publication Critical patent/CN111898440B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/40Scenes; Scene-specific elements in video content
    • G06V20/41Higher-level, semantic clustering, classification or understanding of video scenes, e.g. detection, labelling or Markovian modelling of sport events or news items

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Linguistics (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Multimedia (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Remote Sensing (AREA)
  • Astronomy & Astrophysics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Fire-Detection Mechanisms (AREA)
  • Image Analysis (AREA)

Abstract

本发明属于深度学习技术领域,具体涉及一种基于三维卷积神经网络的山火检测方法。本发明以深度学习取代特征工程,自动提取图像序列的时空特征,大大提升特征描述子的开发效率;以三维卷积取代传统的二维卷积,不仅能学习图像的空间模式,同时也能学习图像序列之间的运动模式,大大提升了视频特征描述子的表达和判别能力;以3DCNN取代传统的二维卷积网络搭建目标检测器,大大提升了山火的检测精度。

Description

一种基于三维卷积神经网络的山火检测方法
技术领域
本发明属于深度学习技术领域,具体涉及一种基于三维卷积神经网络的山火检测方法。
背景技术
森林是维持地球生态平衡的关键。森林火灾和山火不仅造成巨大的经济损失,对生态环境也有严重的破坏。人们曾采用瞭望塔上目视的方式,发现火情和发出报警信息。这种方式不仅成本高昂、效率低下,也容易受人为疏忽的影响。因此,设计一套具有一定自主决策能力的无人化智能山火检测技术,对于维护生态环境正常运转和保持国泰民安至关重要。
为了实现山火的远距离探测,出现了的利用遥感卫星来检测大面积山火的技术。不幸的是,遥感图像的空间分辨率受限于时间分辨率和遥感距离,易受气象条件的影响,同时这类方法无法做到针对关注区域的实时监控,也无法检测早期小面积山火。而地基山火检测技术弥补了天基检测技术的不足,利用可见光或红外摄像机采集图像,通过图像处理或计算机视觉技术检测山火。这类方法早期依赖特征工程技术,提取图像的颜色、纹理、形态等空间特征,或者通过频谱分析提取图像的频谱特征。通过监督学习训练分类器,结合图像金字塔和滑动窗口技术实现山火的分类和定位。近年来随着深度学习的崛起,出现了将以卷积神经网络(ConvolutionNeural Network,CNN)为基础的分类或者目标检测技术引入山火检测领域的尝试。然而早期的方法过于依赖特征工程,难以获得适合山火检测的性能优异的特征描述子;而随后的二维CNN类方法只关注于空间特征的建模,虽然适用于图像级别的分类或检测,但应用在有显著运动特征的山火检测上仍然难以获得较高的精度。
发明内容
为了克服传统的山火检测技术时空特征建模能力低下和误报率过高的不足,本发明提出了一种以三维卷积神经网络(3-Dimensional Convolution Neural Network,3DCNN)为基础的山火检测技术,通过深度学习对山火的时空特征进行建模,获得性能优异的视频特征描述子,并在此基础上实现山火的高精度检测。
为实现上述目的,本发明采用如下技术方案:
一种基于三维卷积神经网络的山火检测方法,包括以下步骤:
S1、构建训练数据集:获取多个相同区域发生山火与未发生山火的视频,其中未发生山火的视频数量大于等于发生山火的视频数量,并分别用标签1和0标注;将所有视频裁剪为长度2秒的视频片段,对发生山火的视频,对视频片段中第一帧中山火区域的坐标进行标注获得训练数据集;
S2、采用10个三维卷积层、5个池化层和Loss层构建三维卷积神经网络:将三维卷积层定义为Convn,池化层定义为Pooln,n是指层数,依次为:Conv1、Pool1、Conv2、Pool2、Conv3、Conv4、Pool3、Conv5、Conv6、Pool4、Conv7、Conv8、Pool5、Conv9、Conv10、Softmax层,其中,Conv1-Conv8的卷积核大小均为3x3x3,步长为[1,1,1],所有池化层的池化核大小为2x2x2,Pool1的步长为[1,2,2],Pool2-Pool5的步长为[2,2,2],Conv9的卷积核大小为1x3x3,步长为[1,1,1],Conv10的卷积核大小为1x1x1,步长为[1,1,1];采用的损失函数包括物体存在与否的损失(Lobj)、类别概率损失(Lclass)和位置大小(Lxywh)的损失,其中,Lobj和Lclass均采用二值交叉熵,其计算公式为:
li=-(yilogxi+(1-yi)log(1-xi)),i∈0,1,...,N-1
其中的x表示预测值,y表示目标值,N表示批量大小;
Lxywh采用均方误差,其计算公式为:
li=(xi-yi)2,i∈0,1,...,N-1
最后总的损失是Lobj、Lclass和Lxywh之和;
S3、对训练数据集,采用随机时空扰动、随机水平翻转及等比例缩放的方式将其规范化为16x224x224的视频片段后,对构建的三维卷积神经网络进行训练:采用以批量大小为30的随机梯度下降法进行训练,采用学习率预热方式,最开始的1000个批量逐渐将学习率从0提升至0.0005,然后保持,直到总迭代次数的90%和95%的时候分别将学习率降低5倍,总的训练周期为200;获得训练好的三维卷积神经网络;
S4、对实时获取的视频,将其分为视频片段后输入训练好的三维卷积神经网络3DNetDet,获得山火检测结果。
本发明的有益效果为:(1)以深度学习取代特征工程,自动提取图像序列的时空特征,大大提升特征描述子的开发效率;(2)以三维卷积取代传统的二维卷积,不仅能学习图像的空间模式,同时也能学习图像序列之间的运动模式,大大提升了视频特征描述子的表达和判别能力;(3)以3DCNN取代传统的二维卷积网络搭建目标检测器,大大提升了山火的检测精度。
附图说明
图1为山火检测流程图;
图2为以3DCNN为基础的分类器示意图;
图3为山火检测子网络示意图;
图4为矩形交集面积和并集面积示意图。
具体实施方式
下面结合附图和实施例对本发明进行进一步说明。
实施例
本例的处理流程如图1所示,在本例中,收集山火视频数据集的方式包括:从互联网新闻网站、视频网站或者大型的森林防火项目数据发布官网获取山火视频,还可以安排专人在保证安全的情况下模拟山火,并通过摄像机录像获取山火视频。同时收集至少同等数量的没有发生山火的关注类型区域的视频。将视频裁剪成长度为2秒的视频片段,分别用标签“0”和标签“1”标注没有山火和有山火的视频片段的类别。如果视频片段有山火,标注视频片段第一帧中山火区域的坐标,包括最上、最左、最下和最右边的坐标。汇集所有标注好的视频片段作为山火视频数据集,在本例中,通过测试集对所提出的方案进行验证,因此随机获取其80%的视频片段作为训练集,剩下的视频片段构成测试集。
搭建以3DCNN为基础的分类网络,将其命名为3DNet。如表1所示是3DNet的网络配置,表中通道顺序均遵循时间、高度、宽度的先后顺序。所有的三维卷积都采用3x3x3大小的卷积核,总共有8个三维卷积层。在卷积层Conv1、Conv2、Conv4、Conv6和Conv8之后分别接最大池化层,5个池化核的大小均为2x2x2,除了第一个池化层以外,其余池化层的池化步长均为2,而第一个池化层在时间维度的步长为1,两个空间维度的步长为2。在最后一个池化层之后接两个包含4096个神经元的全连接层,最后接Softmax层执行分类。
表1分类网络3DNet网络架构
Figure BDA0002561291220000031
Figure BDA0002561291220000041
在大型视频分类数据集Sports-1M上训练3DNet。对于每个视频,从中随机提取5个2秒钟长的视频片段,并缩放至127x171。训练的时候,采用时空随机扰动的方法从视频片段中截取16x112x112的视频片段用于训练,并以50%的概率对视频片段进行水平翻转。以批量大小为30的随机梯度下降(SGD)训练,初始学习率设为0.003,每经过150000次迭代将学习率减半,总共训练20个周期。
搭建以3DNet为骨干网的山火检测网络3DNetDet。目标检测与图像分类相比,需要更多的细节信息,所以将网络的输入大小扩大2倍至16x224x224。将最后一个最大池化层之后的所有层去除,添加用于执行检测任务的模块。表2是3DNetDet网络检测模块配置。检测模块需要预测山火的类别和位置。Loss/Decoder层用于在训练的时候计算损失,或者在推理的时候解码山火的类别概率、位置和大小。3DNetDet的损失包含物体存在与否的损失(Lobj)、类别概率损失(Lclass)和位置大小(Lxywh)等三个部分的损失。其中,Lobj和Lclass均采用二值交叉熵,其计算公式为:
li=-(yi logxi+(1-yi)log(1-xi)),i∈0,1,...,N-1 (公式1)
其中的x表示预测值,y表示目标值,N表示批量大小。Lxywh采用均方误差,其计算公式为:
li=(xi-yi)2,i∈0,1,...,N-1 (公式2)
其中变量的含义同公式1。
表2 3DNetDet检测模块配置
Figure BDA0002561291220000051
在山火视频数据集的训练集上训练3DNetDet。以3DCNN分类网络的骨干网参数初始化3DNetDet的骨干网。对于训练集中的每个视频片段,采用随机时空扰动、随机水平翻转及等比例缩放的方式将其规范化为16x224x224的视频片段。与此同时,采用随机的色度、饱和度和亮度调整,以及对比度增强等方法处理规范化之后的视频片段,作为训练时3DNetDet的输入。以批量大小为30的随机梯度下降(SGD)训练。采用学习率预热技术,最开始的1000个批量逐渐将学习率从0提升至0.0005,然后保持,直到总迭代次数的90%和95%的时候分别将学习率降低5倍。总的训练周期为200。
在山火视频数据集的测试集上测试3DNetDet的检测性能。对于3DNetDet的每个预测值,将置信度低于0.005的预测值移除,同时以非最大值抑制合并同一物体的检测框,然后和山火的标注数据比较,计算准确率和召回率。对于每一个保留的预测值,如果预测框和真实框的最大的交并面积比(Intersection Over Union,IOU)大于0.5,就认为这个预测值是有效的,对应类别的“真正”(True Positive,TP)计数器累加一;反之,则对应的“假正”(False Positive,FP)计数器累加一。测试集中“假负”(False Negative,FN)的数量和TP之和就是测试集中类别标签为“1”的视频片段的总数。分别以公式3和公式4计算准确率和召回率。
Figure BDA0002561291220000061
Figure BDA0002561291220000062
在部署3DNetDet检测山火时,从摄像机获取实时视频流,并将其切分成帧长为16的视频片段,同时要确保相邻的视频片段有8帧的重合。对于3DNetDet输出的预测值,如果其概率大于设定的阈值(阈值的典型值为0.5),就认为测试的视频片段描述了一例山火事件,于是结合对应的预测框位置、防山火设备安装位置信息(经纬度、海拔)、摄像头的内外参数以及地理信息***推算着火点的空间坐标。防山火***将3DNetDet的预测值和着火点空间坐标作为报警信息的一部分发送到相关责任单位。

Claims (1)

1.一种基于三维卷积神经网络的山火检测方法,其特征在于,包括以下步骤:
S1、构建训练数据集:获取多个相同区域发生山火与未发生山火的视频,其中未发生山火的视频数量大于等于发生山火的视频数量,并分别用标签1和0标注;将所有视频裁剪为长度2秒的视频片段,对发生山火的视频,对视频片段中第一帧中山火区域的坐标进行标注获得训练数据集;
S2、采用10个三维卷积层、5个池化层和Loss层构建三维卷积神经网络:将三维卷积层定义为Convn,池化层定义为Pooln,n是指层数,依次为:Conv1、Pool1、Conv2、Pool2、Conv3、Conv4、Pool3、Conv5、Conv6、Pool4、Conv7、Conv8、Pool5、Conv9、Conv10、Softmax层,其中,Conv1-Conv8的卷积核大小均为3x3x3,步长为[1,1,1],所有池化层的池化核大小为2x2x2,Pool1的步长为[1,2,2],Pool2-Pool5的步长为[2,2,2],Conv9的卷积核大小为1x3x3,步长为[1,1,1],Conv10的卷积核大小为1x1x1,步长为[1,1,1];采用的损失函数包括物体存在与否的损失(Lobj)、类别概率损失(Lclass)和位置大小(Lxywh)的损失,其中,Lobj和Lclass均采用二值交叉熵,其计算公式为:
li=-(yilogxi+(1-yi)log(1-xi)),i∈0,1,...,N-1
其中的x表示预测值,y表示目标值,N表示批量大小;
Lxywh采用均方误差,其计算公式为:
li=(xi-yi)2,i∈0,1,...,N-1
总的损失是Lobj、Lclass和Lxywh之和;
S3、对训练数据集,采用随机时空扰动、随机水平翻转及等比例缩放的方式将其规范化为16x224x224的视频片段后,对构建的三维卷积神经网络进行训练:采用以批量大小为30的随机梯度下降法进行训练,采用学习率预热方式,最开始的1000个批量逐渐将学习率从0提升至0.0005,然后保持,直到总迭代次数的90%和95%的时候分别将学习率降低5倍,总的训练周期为200;获得训练好的三维卷积神经网络;
S4、对实时获取的视频,将其分为视频片段后输入训练好的三围卷积神经网络,获得山火检测结果。
CN202010607252.6A 2020-06-30 2020-06-30 一种基于三维卷积神经网络的山火检测方法 Active CN111898440B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010607252.6A CN111898440B (zh) 2020-06-30 2020-06-30 一种基于三维卷积神经网络的山火检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010607252.6A CN111898440B (zh) 2020-06-30 2020-06-30 一种基于三维卷积神经网络的山火检测方法

Publications (2)

Publication Number Publication Date
CN111898440A true CN111898440A (zh) 2020-11-06
CN111898440B CN111898440B (zh) 2023-12-01

Family

ID=73207241

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010607252.6A Active CN111898440B (zh) 2020-06-30 2020-06-30 一种基于三维卷积神经网络的山火检测方法

Country Status (1)

Country Link
CN (1) CN111898440B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113190031A (zh) * 2021-04-30 2021-07-30 成都思晗科技股份有限公司 基于无人机的森林火灾自动拍照跟踪方法、装置及***
CN114973584A (zh) * 2022-05-10 2022-08-30 云南电网有限责任公司电力科学研究院 山火告警方法及装置、计算机设备及存储介质

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106897714A (zh) * 2017-03-23 2017-06-27 北京大学深圳研究生院 一种基于卷积神经网络的视频动作检测方法
CN107480729A (zh) * 2017-09-05 2017-12-15 江苏电力信息技术有限公司 一种基于深度时空域特征的输电线路山火检测方法
CN108764142A (zh) * 2018-05-25 2018-11-06 北京工业大学 基于3dcnn的无人机影像森林烟雾检测和分类方法
CN109389185A (zh) * 2018-11-15 2019-02-26 中国科学技术大学 使用三维卷积神经网络的视频烟雾识别方法
CN109522819A (zh) * 2018-10-29 2019-03-26 西安交通大学 一种基于深度学习的火灾图像识别方法
CN109829583A (zh) * 2019-01-31 2019-05-31 成都思晗科技股份有限公司 基于概率编程技术的山火风险预测方法
CN109919993A (zh) * 2019-03-12 2019-06-21 腾讯科技(深圳)有限公司 视差图获取方法、装置和设备及控制***
CN110570615A (zh) * 2019-09-04 2019-12-13 云南电网有限责任公司带电作业分公司 天空地联合的输电线路通道山火趋势预警方法、装置、***及存储介质
US20200110966A1 (en) * 2018-10-09 2020-04-09 Naver Corporation Systems and methods for detecting a point of interest change using a convolutional neural network

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106897714A (zh) * 2017-03-23 2017-06-27 北京大学深圳研究生院 一种基于卷积神经网络的视频动作检测方法
CN107480729A (zh) * 2017-09-05 2017-12-15 江苏电力信息技术有限公司 一种基于深度时空域特征的输电线路山火检测方法
CN108764142A (zh) * 2018-05-25 2018-11-06 北京工业大学 基于3dcnn的无人机影像森林烟雾检测和分类方法
US20200110966A1 (en) * 2018-10-09 2020-04-09 Naver Corporation Systems and methods for detecting a point of interest change using a convolutional neural network
CN109522819A (zh) * 2018-10-29 2019-03-26 西安交通大学 一种基于深度学习的火灾图像识别方法
CN109389185A (zh) * 2018-11-15 2019-02-26 中国科学技术大学 使用三维卷积神经网络的视频烟雾识别方法
CN109829583A (zh) * 2019-01-31 2019-05-31 成都思晗科技股份有限公司 基于概率编程技术的山火风险预测方法
CN109919993A (zh) * 2019-03-12 2019-06-21 腾讯科技(深圳)有限公司 视差图获取方法、装置和设备及控制***
CN110570615A (zh) * 2019-09-04 2019-12-13 云南电网有限责任公司带电作业分公司 天空地联合的输电线路通道山火趋势预警方法、装置、***及存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
XIUQING LI: "3D parallel fully convolutional networks for real-time video wildfire smoke detection", 《IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY》, vol. 30, no. 1, pages 89 - 103, XP011766409, DOI: 10.1109/TCSVT.2018.2889193 *
付光远;辜弘炀;汪洪桥;: "基于卷积神经网络的高光谱图像谱-空联合分类", 科学技术与工程, no. 21, pages 273 - 279 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113190031A (zh) * 2021-04-30 2021-07-30 成都思晗科技股份有限公司 基于无人机的森林火灾自动拍照跟踪方法、装置及***
CN114973584A (zh) * 2022-05-10 2022-08-30 云南电网有限责任公司电力科学研究院 山火告警方法及装置、计算机设备及存储介质

Also Published As

Publication number Publication date
CN111898440B (zh) 2023-12-01

Similar Documents

Publication Publication Date Title
CN109522819B (zh) 一种基于深度学习的火灾图像识别方法
CN112287816B (zh) 一种基于深度学习的危险工作区事故自动检测报警方法
Zhu et al. Msnet: A multilevel instance segmentation network for natural disaster damage assessment in aerial videos
Zhan et al. A high-precision forest fire smoke detection approach based on ARGNet
CN110827505A (zh) 一种基于深度学习的烟雾分割方法
CN111126293A (zh) 一种火焰和烟雾异常情况检测方法及***
CN113963301A (zh) 一种时空特征融合的视频火灾烟雾探测方法及***
CN111898440B (zh) 一种基于三维卷积神经网络的山火检测方法
Qiang et al. Forest fire smoke detection under complex backgrounds using TRPCA and TSVB
CN112541403A (zh) 一种利用红外摄像头的室内人员跌倒检测方法
CN115690564A (zh) 基于Recursive BIFPN网络的室外火灾烟雾图像检测方法
Luo Research on fire detection based on YOLOv5
CN117612249A (zh) 一种基于改进OpenPose算法的井下矿工危险行为识别方法和装置
CN116994209A (zh) 基于人工智能的图像数据处理***及其方法
Munawar et al. Fire detection through Image Processing; A brief overview
KR20240028290A (ko) 산불 검출 장치 및 방법
Supangkat et al. Moving Image Interpretation Models to Support City Analysis
CN116188442A (zh) 一种适用于任意场景下的高精度森林烟火检测方法
CN116704360A (zh) 一种基于时空特征融合的卫星遥感火点实时监测方法
CN116469029A (zh) 一种基于多层感知机网络的油田现场违章行为识别方法
Xue et al. A Fire Detection and Assessment Method based on YOLOv8
CN110852174A (zh) 一种基于视频监控的早期烟雾检测方法
Song et al. PDD: Post-Disaster Dataset for Human Detection and Performance Evaluation
Rajan et al. Forest fire detection using machine learning
CN115995051A (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