CN109489812A - 扩大火焰探测器探测视角的方法 - Google Patents
扩大火焰探测器探测视角的方法 Download PDFInfo
- Publication number
- CN109489812A CN109489812A CN201811384073.XA CN201811384073A CN109489812A CN 109489812 A CN109489812 A CN 109489812A CN 201811384073 A CN201811384073 A CN 201811384073A CN 109489812 A CN109489812 A CN 109489812A
- Authority
- CN
- China
- Prior art keywords
- degree
- detection
- detection angle
- angle
- eigenvalue
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 143
- 230000000007 visual effect Effects 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000012512 characterization method Methods 0.000 claims abstract description 9
- 239000000284 extract Substances 0.000 claims abstract description 7
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 235000019504 cigarettes Nutrition 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000003053 immunization Effects 0.000 description 1
- 238000002649 immunization Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Fire-Detection Mechanisms (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明提供一种扩大火焰探测器探测视角的方法,包括:提取表征探测角度变化的第一特征值,建立多个不同探测角度关于第一特征值的隶属度函数;通过多个不同探测角度关于第一特征值的隶属度函数,计算预测火焰探测器的实际探测角度;提取表征火焰能量的第二特征值;建立阈值关于探测角度的隶属度函数,所述阈值在探测角度为0度时最大,随探测角度增大而减小;通过阈值关于探测角度的隶属度函数计算当前预测的实际探测角度对应的阈值,若第二特征值大于等于该对应阈值,则该实际探测角度的火焰能够被检测。本发明在不增加成本的情况下,有效地扩大火焰探测器的探测视角。
Description
技术领域
本发明涉及火焰探测领域,具体地说,涉及一种扩大火焰探测器探测视角的方法。
背景技术
光学火焰探测器是根据物质在燃烧过程中发出的光谱来探测火灾是否发生。根据光谱波长的不同,大致可以分为紫外火焰探测器,紫红外火焰探测器,红外火焰探测器以及视频火焰探测器等。相比于烟感和温感火灾探测器,利用光谱检测的方法灵敏度高,探测距离远,能够对风、雨、高温、高湿以及环境光和人工光源等干扰信号具有很好的免疫作用,并且使用寿命长。
通常点型火焰探测器,需要多个探测器组合安装在一个空间区域内,如果每个火焰探测器的探测视角不稳定或者探测视角较小,那么就会带来视角盲区,即使在发生火灾的时候,也会带来漏报的风险。
发明内容
本发明的目的在于克服现有技术中存在的不足,提供一种扩大火焰探测器探测视角的方法,采用自适应阈值,在不改变产品硬件的条件下,可以有效地扩大火焰探测器的探测视角。本发明采用的技术方案是:
一种扩大火焰探测器探测视角的方法,包括:
步骤S1,提取表征探测角度变化的第一特征值,建立多个不同探测角度关于第一特征值的隶属度函数;
步骤S2,通过多个不同探测角度关于第一特征值的隶属度函数,计算预测火焰探测器的实际探测角度;
步骤S3,提取表征火焰能量的第二特征值;建立阈值关于探测角度的隶属度函数,所述阈值在探测角度为0度时最大,随探测角度增大而减小;通过阈值关于探测角度的隶属度函数计算当前预测的实际探测角度对应的阈值,若第二特征值大于等于该对应阈值,则该实际探测角度的火焰能够被检测。
进一步地,第一特征值采用两个不同探测通道的探测信号时域值经过傅里叶变换,取设定的火焰闪烁频率范围内的频谱最大值之比,比值记为λ,作为第一特征值。
进一步地,所述建立多个不同探测角度关于第一特征值的隶属度函数,具体包括:
从中间向两侧划分不同的探测角度,在每个探测角度同等距离采集多组探测信号的数据样本,计算出每个探测角度对应第一特征值λ的最值和均值;
探测角度0度时λ的范围为[k1,k2],均值为s;±α1度处的λ的范围为[k3,k4],-α1度处的λ均值为s1,+α1度处的λ均值为s2;±α2度处的λ的范围为[k5,k6],-α2度处的λ均值为s3,+α2度处的λ均值为s4;±α3度处的λ的范围为[k7,k8],-α3度处的λ均值为s5,+α3度处的λ均值为s6;α3>α2>α1>0;;假定s1小于s2,s3小于s4;
多个不同探测角度关于第一特征值的隶属度函数如下式表示:
其中,P0(λ)为探测角度0度关于λ的隶属度函数,Pα1(λ)表示探测角度α1度关于λ的隶属度函数,Pα2(λ)表示探测角度α2度关于λ的隶属度函数,Pα3(λ)表示探测角度α3度关于λ的隶属度函数。
进一步地,所述计算预测火焰探测器的实际探测角度,具体包括:
根据实际采集的探测信号样本根据式(1)、(2)、(3)、(4)计算0、α1、α2、α3的隶属度值得到一个一维数组Pα,记作:
Pα=[P0,Pα1,Pα2,Pα3] (5)
根据一维数组Pα,计算预测出该探测信号样本的实际探测角度α′,记为:
进一步地,步骤S3具体包括:
根据采集的各探测角度处的各n组探测信号样本的数据,提取第二特征值为采集的该探测角度下火焰信号时域的峰峰值;
分别计算探测角度在0度、±α1度、±α2度、±α3度处时域峰峰值的均值,记为M1,M2,M3,M4;令M1,M2,M3,M4分别为0度、±α1度、±α2度、±α3度处的阈值;M1>M2>M3>M4;
各阈值关于探测角度α的隶属度函数,用公式表示为:
其中,PM1(α),PM2(α),PM3(α),PM4(α)为分别为阈值M1,M2,M3,M4关于探测角度α的隶属度函数;
然后,通过计算预测的实际探测角度α′代入式(7)、(8)、(9)、(10),计算阈值M1、M2、M3、M4的隶属度值得到一个一维数组PM,记作:
PM=[PM1,PM2,PM3,PM4] (11)
根据得到的一维数组PM计算出在该预测的实际探测角度下对应的阈值M,计算公式如下:
若实际提取的第二特征值大于等于该对应阈值M,则该预测的实际探测角度的火焰能够被检测。
本发明的优点在于:
1)火焰探测器的成本无需增加。
2)解决了多通道火焰探测器的探测视角较小或者探测视角不稳定的问题,有效地增加了探测视角。
附图说明
图1为本发明的探测视角划分示意图。
图2为本发明的探测角度隶属度函数示意图示意图。
图3为本发明的不同阈值的隶属度函数示意图。
图4为本发明的扩大探测视角的算法流程图。
具体实施方式
下面结合具体附图和实施例对本发明作进一步说明。
本发明提出的扩大火焰探测器探测视角的方法,是在火焰探测器具备探测火焰能力的基础上,进一步扩大其探测视角,避免探测视角集中于火焰中心部位;
通常,为了能够探测火焰,火焰探测器包含四个探测通道的红外传感器,其中通道1为4.4~4.6μm,通道2为4.1~4.3μm,通道3为2.0~2.4μm,通道4为4.8~5.4μm;在其它的实施例中,也可以采用紫外传感器等其它的探测光传感器;
传感器信号经过滤波放大,模数转换,再通过MCU进行数据的采集,再将数据传送给上位机进行分析;
(一)提取表征探测角度变化的第一特征值,建立多个不同探测角度关于第一特征值的隶属度函数;
如图1所示,首先根据与中心轴的夹角划分多个探测角度,中心轴对应的探测角度为0度;选用两个不同的探测通道,利用这两个探测通道对于探测角度增大能量衰减不同的特点,找到表征两个不同探测通道对于同一火焰在相同探测角度下探测信号差异的特征值,且该特征值的变化与探测角度的变化存在相关性;
本发明选用的第一特征值是两个不同探测通道的探测信号时域值经过傅里叶变换,取设定的火焰闪烁频率范围(通常为1~20Hz)内的频谱最大值之比,比值记为λ,作为第一特征值;
然后,基于探测信号的多组数据样本,建立多个不同探测角度关于第一特征值的隶属度函数;具体如下:
从中间向两侧划分不同的探测角度,在每个探测角度同等距离采集多组探测信号的数据样本,计算出每个探测角度对应第一特征值的最值和均值;
如图1所示,根据与中心轴的夹角划分探测角度,左侧为负,右侧为正,在0度、±α1度、±α2度、±α3度处采集同样大小火焰的探测信号数据样本;每处采集n组样本;本例中每处采集10组样本;0度角为中心角;中心轴是探测光传感器与火焰中心的轴线;α3>α2>α1>0;
分别计算各探测角度处的各n组样本的λ值,并计算出每处n组样本的λ的最大值,最小值以及均值;
0度时λ的范围为[k1,k2],均值为s;±α1度处的λ的范围为[k3,k4],-α1度处的λ均值为s1,+α1度处的λ均值为s2;±α2度处的λ的范围为[k5,k6],-α2度处的λ均值为s3,+α2度处的λ均值为s4;±α3度处的λ的范围为[k7,k8],-α3度处的λ均值为s5,+α3度处的λ均值为s6;假定s1小于s2,s3小于s4;
如图2所示,多个不同探测角度关于第一特征值的隶属度函数如下式表示:
其中,P0(λ)为探测角度0度关于λ的隶属度函数,Pα1(λ)表示探测角度α1度关于λ的隶属度函数,Pα2(λ)表示探测角度α2度关于λ的隶属度函数,Pα3(λ)表示探测角度α3度关于λ的隶属度函数。探测角度的隶属度函数采用梯形隶属度函数。
(二)利用建立的多个不同探测角度关于第一特征值的隶属度函数,计算预测火焰探测器的实际探测角度;
根据实际采集的探测信号样本根据式(1)、(2)、(3)、(4)计算0、α1、α2、α3的隶属度值得到一个一维数组Pα,记作:
Pα=[P0,Pα1,Pα2,Pα3] (5)
根据一维数组Pα,计算预测出该探测信号样本的实际探测角度α′,记为:
(三)根据表征火焰能量的第二特征值设计阈值,建立阈值关于探测角度的隶属度函数;该阈值是自适应阈值,能够随探测角度增大而减小;
本发明需要根据探测角度的值来调整阈值M的大小;具体如下:
根据采集的各探测角度处的各n组探测信号样本的数据,提取第二特征值为采集的该探测角度下火焰信号时域的峰峰值;
分别计算探测角度在0度、±α1度、±α2度、±α3度处时域峰峰值的均值,记为M1,M2,M3,M4;令M1,M2,M3,M4分别为0度、±α1度、±α2度、±α3度处的阈值;M1>M2>M3>M4;
如图3所示,各阈值关于探测角度α的隶属度函数,用公式表示为:
其中,PM1(α),PM2(α),PM3(α),PM4(α)为分别为阈值M1,M2,M3,M4关于探测角度α的隶属度函数。各阈值关于探测角度α的隶属度函数采用三角形隶属度函数。
然后,通过计算预测的实际探测角度α′代入式(7)、(8)、(9)、(10),计算阈值M1、M2、M3、M4的隶属度值得到一个一维数组PM,记作:
PM=[PM1,PM2,PM3,PM4] (11)
根据得到的一维数组PM计算出在该预测的实际探测角度下对应的阈值M,计算公式如下:
若实际提取的第二特征值大于等于该对应阈值M,则认为在不改变原始火焰大小的情况下,该预测的实际探测角度的火焰能够被检测。
扩大火焰探测器探测视角的算法如图4所示,
步骤一,判断检测到的信号是否为火焰信号;
步骤二,提取第一特征值,计算预测当前的实际探测角度;
步骤三,计算当前实际探测角度对应的阈值M;
步骤四,提取第二特征值,判断第二特征值是否大于M;
步骤五,如果大于M,则认为该实际探测角度下的火焰能够被识别;
中心角0度的阈值M最大,随着探测角度α的增大,对应的阈值M减小;如果提取样本的时域峰峰值大于阈值M,则判断该实际探测角度下的火焰能够被探测。相比于传统的采用固定阈值的方法,这种方法能够有效地增大边沿角度识别的距离,从而扩大火焰探测器的探测视角,解决传统火焰探测器中心角时识别距离远,而两侧识别距离近的问题。
本发明在不改变产品硬件的条件下,可以有效地扩大火焰探测器的探测视野,减小视野盲区。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种扩大火焰探测器探测视角的方法,其特征在于,包括:
步骤S1,提取表征探测角度变化的第一特征值,建立多个不同探测角度关于第一特征值的隶属度函数;
步骤S2,通过多个不同探测角度关于第一特征值的隶属度函数,计算预测火焰探测器的实际探测角度;
步骤S3,提取表征火焰能量的第二特征值;建立阈值关于探测角度的隶属度函数,所述阈值在探测角度为0度时最大,随探测角度增大而减小;通过阈值关于探测角度的隶属度函数计算当前预测的实际探测角度对应的阈值,若第二特征值大于等于该对应阈值,则该实际探测角度的火焰能够被检测。
2.如权利要求1所述的扩大火焰探测器探测视角的方法,其特征在于,
第一特征值采用两个不同探测通道的探测信号时域值经过傅里叶变换,取设定的火焰闪烁频率范围内的频谱最大值之比,比值记为λ,作为第一特征值。
3.如权利要求2所述的扩大火焰探测器探测视角的方法,其特征在于,
所述建立多个不同探测角度关于第一特征值的隶属度函数,具体包括:
从中间向两侧划分不同的探测角度,在每个探测角度同等距离采集多组探测信号的数据样本,计算出每个探测角度对应第一特征值λ的最值和均值;
探测角度0度时λ的范围为[k1,k2],均值为s;±α1度处的λ的范围为[k3,k4],-α1度处的λ均值为s1,+α1度处的λ均值为s2;±α2度处的λ的范围为[k5,k6],-α2度处的λ均值为s3,+α2度处的λ均值为s4;±α3度处的λ的范围为[k7,k8],-α3度处的λ均值为s5,+α3度处的λ均值为s6;α3>α2>α1>0;假定s1小于s2,s3小于s4;
多个不同探测角度关于第一特征值的隶属度函数如下式表示:
其中,P0(λ)为探测角度0度关于λ的隶属度函数,Pα1(λ)表示探测角度α1度关于λ的隶属度函数,Pα2(λ)表示探测角度α2度关于λ的隶属度函数,Pα3(λ)表示探测角度α3度关于λ的隶属度函数。
4.如权利要求3所述的扩大火焰探测器探测视角的方法,其特征在于,
所述计算预测火焰探测器的实际探测角度,具体包括:
根据实际采集的探测信号样本通过式(1)、(2)、(3)、(4)计算0、α1、α2、α3的隶属度值得到一个一维数组Pα,记作:
Pα=[P0,Pα1,Pα2,Pα3] (5)
根据一维数组Pα,计算预测出该探测信号样本的实际探测角度α′,记为:
5.如权利要求4所述的扩大火焰探测器探测视角的方法,其特征在于,
步骤S3具体包括:
根据采集的各探测角度处的各n组探测信号样本的数据,提取第二特征值为采集的该探测角度下火焰信号时域的峰峰值;
分别计算探测角度在0度、±α1度、±α2度、±α3度处时域峰峰值的均值,记为M1,M2,M3,M4;令M1,M2,M3,M4分别为0度、±α1度、±α2度、±α3度处的阈值;M1>M2>M3>M4;
各阈值关于探测角度α的隶属度函数,用公式表示为:
其中,PM1(α),PM2(α),PM3(α),PM4(α)为分别为阈值M1,M2,M3,M4关于探测角度α的隶属度函数;
然后,通过计算预测的实际探测角度α′代入式(7)、(8)、(9)、(10),计算阈值M1、M2、M3、M4的隶属度值得到一个一维数组PM,记作:
PM=[PM1,PM2,PM3,PM4] (11)
根据得到的一维数组PM计算出在该预测的实际探测角度下对应的阈值M,计算公式如下:
若实际提取的第二特征值大于等于该对应阈值M,则该预测的实际探测角度的火焰能够被检测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811384073.XA CN109489812B (zh) | 2018-11-20 | 2018-11-20 | 扩大火焰探测器探测视角的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811384073.XA CN109489812B (zh) | 2018-11-20 | 2018-11-20 | 扩大火焰探测器探测视角的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109489812A true CN109489812A (zh) | 2019-03-19 |
CN109489812B CN109489812B (zh) | 2021-03-16 |
Family
ID=65696331
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811384073.XA Active CN109489812B (zh) | 2018-11-20 | 2018-11-20 | 扩大火焰探测器探测视角的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109489812B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111664935A (zh) * | 2020-06-23 | 2020-09-15 | 安东仪器仪表检测有限公司 | 一种固定安装火焰探测器的原位校准方法及光源 |
CN115100811A (zh) * | 2022-06-22 | 2022-09-23 | 招商局重庆公路工程检测中心有限公司 | 一种公路隧道火焰探测器的探测空间调试方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4585966B2 (ja) * | 2005-12-27 | 2010-11-24 | ホーチキ株式会社 | 減光式煙感知器 |
CN105185027A (zh) * | 2015-09-11 | 2015-12-23 | 武汉理工光科股份有限公司 | 火焰探测报警方法、系统及装置 |
CN107025428A (zh) * | 2015-11-25 | 2017-08-08 | A.M.总建筑股份公司 | 用于有限空间的具有复合功能的红外辐射火灾检测器 |
CN206695909U (zh) * | 2017-04-26 | 2017-12-01 | 福建天广消防有限公司 | 一体式红外火焰探测装置 |
CN206788036U (zh) * | 2016-12-19 | 2017-12-22 | 中国科学技术大学 | 一种火灾烟雾散射和消光特性测量装置 |
-
2018
- 2018-11-20 CN CN201811384073.XA patent/CN109489812B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4585966B2 (ja) * | 2005-12-27 | 2010-11-24 | ホーチキ株式会社 | 減光式煙感知器 |
CN105185027A (zh) * | 2015-09-11 | 2015-12-23 | 武汉理工光科股份有限公司 | 火焰探测报警方法、系统及装置 |
CN107025428A (zh) * | 2015-11-25 | 2017-08-08 | A.M.总建筑股份公司 | 用于有限空间的具有复合功能的红外辐射火灾检测器 |
CN206788036U (zh) * | 2016-12-19 | 2017-12-22 | 中国科学技术大学 | 一种火灾烟雾散射和消光特性测量装置 |
CN206695909U (zh) * | 2017-04-26 | 2017-12-01 | 福建天广消防有限公司 | 一体式红外火焰探测装置 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111664935A (zh) * | 2020-06-23 | 2020-09-15 | 安东仪器仪表检测有限公司 | 一种固定安装火焰探测器的原位校准方法及光源 |
CN111664935B (zh) * | 2020-06-23 | 2022-09-09 | 安东仪器仪表检测有限公司 | 一种固定安装火焰探测器的原位校准方法及光源 |
CN115100811A (zh) * | 2022-06-22 | 2022-09-23 | 招商局重庆公路工程检测中心有限公司 | 一种公路隧道火焰探测器的探测空间调试方法及装置 |
CN115100811B (zh) * | 2022-06-22 | 2024-01-30 | 招商局重庆公路工程检测中心有限公司 | 一种公路隧道火焰探测器的探测空间调试方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109489812B (zh) | 2021-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105046868B (zh) | 一种狭长环境中基于红外热像仪的火灾预警方法 | |
Kashani et al. | Cluster-based roof covering damage detection in ground-based lidar data | |
CN113160513A (zh) | 基于多传感器的火焰探测装置 | |
US9224278B2 (en) | Automated method and system for detecting the presence of a lit cigarette | |
CN109489812A (zh) | 扩大火焰探测器探测视角的方法 | |
CN101770644A (zh) | 森林火灾远程视频监控烟火识别方法 | |
JP5308898B2 (ja) | 人検知センサ | |
WO2012125731A2 (en) | Mwir sensor for flame detection | |
CN108335454B (zh) | 一种火情检测方法及装置 | |
CN112542020A (zh) | 一种基于双波段的红外火焰探测方法及系统 | |
CN111951250A (zh) | 一种基于图像的火灾检测方法 | |
CN115097526A (zh) | 地震采集资料质量评价方法 | |
CN113218508A (zh) | 一种基于Himawari-8卫星数据的林火判别方法 | |
CN106097630A (zh) | 一种双红外通道火焰探测器的火焰识别方法 | |
CN209433517U (zh) | 一种基于多火焰图像与组合判据的火灾识别报警装置 | |
US5838242A (en) | Fire detection system using modulation ratiometrics | |
CN110031112A (zh) | 一种五波段火焰探测器及其探测方法 | |
CN110285755B (zh) | 一种用于筛选关联成像采样数据的装置和方法 | |
CN115778374B (zh) | 一种跌倒检测方法及系统 | |
CN111437556A (zh) | 一种火灾探测器、火灾探测方法和自动灭火系统 | |
CN115620102A (zh) | 一种基于多模态信息融合技术的火焰检测系统 | |
CN115841609A (zh) | 基于红外热图像分析的石化区火源检测方法、系统及应用 | |
CN201382895Y (zh) | 一种目标微生物的检测设备 | |
KR102171066B1 (ko) | 딥 러닝을 적용한 지진파형 탐지 방법 | |
CN104850826B (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 |