CN111855608B - 基于融合特征波长选择算法的苹果酸度近红外无损检测方法 - Google Patents

基于融合特征波长选择算法的苹果酸度近红外无损检测方法 Download PDF

Info

Publication number
CN111855608B
CN111855608B CN202010745439.2A CN202010745439A CN111855608B CN 111855608 B CN111855608 B CN 111855608B CN 202010745439 A CN202010745439 A CN 202010745439A CN 111855608 B CN111855608 B CN 111855608B
Authority
CN
China
Prior art keywords
wavelength
apple
acidity
spectrum
characteristic
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.)
Active
Application number
CN202010745439.2A
Other languages
English (en)
Other versions
CN111855608A (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.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
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 Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN202010745439.2A priority Critical patent/CN111855608B/zh
Publication of CN111855608A publication Critical patent/CN111855608A/zh
Application granted granted Critical
Publication of CN111855608B publication Critical patent/CN111855608B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/253Fusion techniques of extracted features
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/68Food, e.g. fruit or vegetables

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于融合特征波长选择算法的苹果酸度近红外无损检测方法,其步骤包括:1采集苹果样本标记点区域的光谱信息,测量苹果样本标记点区域的酸度数据;2对采集到的光谱进行预处理;3分别利用连续投影算法SPA与竞争自适应重加权采样算法CARS进行特征波长选择,并将二者所选择的特征波长进行融合;4根据融合后的特征波长对应的光谱与酸度数据,在校正集上建立苹果酸度的偏最小二乘PLS预测模型,并在预测集上对模型结果进行评估。本发明能兼顾选取的波长变量数目与建立的模型准确性,模型简单、检测效率高、实用性强,可为快速无损检测苹果的酸度提供重要手段,降低特定用途便携式近红外仪器的开发成本。

Description

基于融合特征波长选择算法的苹果酸度近红外无损检测方法
技术领域
本发明实施例涉及苹果品质检测技术领域,具体涉及一种基于融合特征波长选择算法的苹果酸度近红外无损检测方法。
背景技术
酸度是衡量苹果内部品质的重要指标之一,影响着消费者的购买意愿。传统的苹果酸度测量方法是破坏式的,且费时费力。近红外光谱技术以其快速无损的优势在苹果内部品质检测方面得到了广泛应用,但较少将其运用于对酸度的检测。
在对基于近红外光谱的无损检测模型的建立过程中,波长变量选择是重要的一环。通过特定方法筛选特征波长或波长区间,一方面可以简化模型,另一方面由于不相关或非线性变量的剔除,可以得到预测能力强、稳健性好的校正模型。
在特征波长选择领域,最常使用的是连续投影算法SPA与竞争自适应重加权采样算法CARS,这两种方法各有优缺点。通过SPA算法选出的波长变量数目较少,但最终建立的模型准确性欠佳;通过CARS算法进行特征波长选择后最终建立的模型准确性较高,然而选出的波长变量数目也很多;若将二者联合使用,组成CARS-SPA算法,波长变量数目可以有效减少,但模型准确性也随之大幅下降。在苹果酸度无损检测领域,目前单一的或联合的特征波长选择算法均无法兼顾模型准确性与波长变量数目,存在或预测精度低,或模型复杂、检测效率低、实用性弱等问题,相应的便携式近红外仪器的开发成本相对较高。
发明内容
本发明是为了解决上述现有技术存在的不足之处,提出一种基于融合特征波长选择算法的苹果酸度近红外无损检测方法,以期能降低模型复杂度、提高模型检测效率与实用性,从而能实现对苹果酸度的快速无损检测,降低特定用途便携式近红外仪器的开发成本。
本发明为达到上述发明目的,采用如下技术方案:
本发明一种基于融合特征波长选择算法的苹果酸度近红外无损检测方法的特点也在于,步骤如下:
步骤一、采集苹果样本标记点区域的原始近红外反射率光谱;
步骤二、测量所述苹果样本标记点区域的酸度数据;
步骤三、对所述原始近红外反射率光谱进行预处理,得到预处理后的光谱;所述预处理包括:选取波长范围、将反射率光谱转换为吸光度光谱、对吸光度光谱进行归一化、平滑、散射校正的处理;
步骤四、利用连续投影算法与竞争自适应重加权采样算法分别对预处理后的光谱进行特征波长选择,并得到SPA特征波长集合S和CARS特征波长集合C;
步骤五、将SPA特征波长集合S和CARS特征波长集合C进行特征波长的融合,得到融合后的特征波长;
步骤六、以所述融合后的特征波长所对应的预处理后的光谱作为模型的输入,以所述酸度数据作为模型的输出,建立苹果酸度预测模型,从而用于苹果酸度的无损检测。
本发明所述的苹果酸度近红外无损检测方法的特点也在于:所述步骤五是按如下过程进行:
(1)令建模波长集合A=S,令剩余波长集合B=C–S,即集合C与集合S的差集;定义数组r,定义变量i,定义阈值ΔR;
(2)用所述建模波长集合A中波长对应的预处理后的光谱建立校正模型,并计算所述校正模型的决定系数R2
(3)令k为所述剩余波长集合B中元素的个数;初始化变量i=1;
(4)从所述剩余波长集合B中取出第i个波长加入所述建模波长集合A中,从而得到加入波长后的集合Ai
(5)用所述加入波长后的集合Ai中波长对应的预处理后的光谱建立第i个校正模型,并计算所述第i个校正模型的决定系数Ri 2后加入所述数组r中;
(6)将i+1复制给i,并判断i>k是否成立,若成立,则执行步骤(7),否则,返回步骤(4);
(7)取出所述数组r中的最大值,记为rmax,并判断rmax>R2是否成立,若成立,则执行步骤(8);否则,算法终止,从而得到建模波长集合A即为融合后的特征波长;
(8)将rmax所对应于剩余波长集合B中的第max个波长加入建模波长集合A,从而更新建模波长集合A,并将第max个波长从剩余波长集合B中剔除,从而更新剩余波长集合B;令R2=rmax,清空数组r;
(9)判断R2>ΔR是否成立,若成立,算法终止,从而得到建模波长集合A即为融合后的特征波长;否则,重复步骤(3)顺序执行,直至剩余波长集合B为空。
与现有技术相比,本发明的有益效果在于:
1、本发明提出的基于融合特征波长选择算法的苹果酸度无损检测方法,相比于传统的破坏式的苹果酸度检测方法,具有快速、无损、经济环保等优点;基于六个步骤建立预测模型,兼顾了选取的波长变量数目与建立的模型准确性,即在尽可能少的波长数目下得到准确性尽可能高的模型,具有极强的可行性。在预测模型的建立过程中,使用融合特征波长选择算法对光谱进行波长变量的筛选,克服了单一的或联合的波长选择算法均无法兼顾模型准确性与波长变量数目这一不足,较其它波长选择方法下建立的模型简单、预测精度高、检测速度快、实用性强。该方法能实现对苹果酸度的快速无损检测,压缩特定用途便携式近红外仪器的开发成本。
2、本发明融合特征波长选择算法考虑到了设计原理不同的SPA算法与CARS算法各自存在的优缺点,有效地将二者的优点结合了起来。基于该算法选取的波长变量建立的预测模型,与SPA算法下建立的模型相比,准确性有了大幅提升;与CARS算法下建立的模型相比,准确性虽略有损失,但选择的波长变量数目大幅减少,满足在尽可能少的波长数目下得到准确性尽可能高的模型这一思想。该算法在对苹果酸度的快速无损检测的过程中发挥着重要作用。
附图说明
图1为本发明苹果酸度近红外无损检测方法流程图;
图2为本发明实施例中涉及的黄金帅苹果的原始近红外反射率光谱;
图3为本发明实施例中涉及的黄金帅苹果经预处理后的近红外吸光度光谱;
图4为本发明实施例中涉及的融合特征波长选择算法选取的波长变量;
图5为本发明实施例中基于PLS模型得出的黄金帅苹果酸度的实测值与预测值的散点图。
具体实施方式
本实施例中,如图1所示,一种基于融合特征波长选择算法的苹果酸度近红外无损检测方法,步骤如下:
步骤一、采集苹果样本标记点区域的原始近红外反射率光谱。选取从超市采购的“黄金帅”苹果31个,在每个苹果样本赤道上均匀分布的四个标记点区域分别进行原始近红外反射率光谱采集,得到124条光谱数据。所用光谱仪为美国海洋光学公司生产的USB4000光谱仪,光谱采集范围为346~1046nm,分辨率2nm。所用采集软件为光谱仪配套软件OceanView,软件中光谱采集的积分时间设置为30ms,平均次数设置为5,平滑度设置为5。所得光谱如图2所示。
步骤二、测量苹果样本标记点区域的酸度数据。具体实施中,所用仪器为日本ATAGO公司生产的PAL-BX/ACID5型苹果糖酸一体机。取样本标记点区域的一小块果肉,挤出质量为1g的果汁进烧杯中,并加入去离子水至烧杯中的混合液质量为51g为止,用搅拌棒充分搅拌混合液,用滴管取出混合液滴3-4滴于仪器上,读出仪器读数,即为该点的酸度。共采集到124条酸度数据。
步骤三、对原始近红外反射率光谱进行预处理,得到预处理后的光谱。在此之前,首先剔除明显异常数据,保留样本数据110条,并用SPXY方法进行样本划分,划分校正集与预测集样本的比例为4:1。酸度数据如表1所示。
表1校正集和预测集样本参数统计
组别 样本数 最小值/% 最大值/% 平均值/% 标准差/%
校正集 88 0.10 0.68 0.30 0.12
预测集 22 0.10 0.44 0.24 0.08
由于原始光谱两端存在较大噪声,会影响建模效果,所以进行波长范围的选取,选取用于建模波长范围为475~925nm,该区域内共有波长变量2335个。将反射率光谱转换为吸光度光谱,转换公式为:
Figure BDA0002608188500000041
其中R为反射率,A为吸光度。为了校正由微小光程差异引起的光谱变化,对吸光度光谱进行归一化操作;为了消除光谱信号中叠加的随机误差,提高信噪比,对光谱进行窗口大小为5的SG卷积平滑;再进行多元散射校正MSC预处理以消除苹果表面散射及光程变化对光谱的影响。得到norm-sg-msc预处理后光谱,如图3所示。
步骤四、利用连续投影算法与竞争自适应重加权采样算法分别对预处理后的光谱进行特征波长选择,并得到SPA特征波长集合S和CARS特征波长集合C。
步骤五、将SPA特征波长集合S和CARS特征波长集合C进行特征波长的融合,得到融合后的特征波长。具体是按如下过程进行:
(1)令建模波长集合A=S,令剩余波长集合B=C–S,即集合C与集合S的差集;定义数组r,定义变量i,定义阈值ΔR;
(2)用建模波长集合A中波长对应的预处理后的光谱建立校正模型,并计算校正模型的决定系数R2
(3)令k为剩余波长集合B中元素的个数;初始化变量i=1;
(4)从剩余波长集合B中取出第i个波长加入建模波长集合A中,从而得到加入波长后的集合Ai
(5)用加入波长后的集合Ai中波长对应的预处理后的光谱建立第i个校正模型,并计算第i个校正模型的决定系数Ri 2后加入数组r中;
(6)将i+1复制给i,并判断i>k是否成立,若成立,则执行步骤(7),否则,返回步骤(4);
(7)取出数组r中的最大值,记为rmax,并判断rmax>R2是否成立,若成立,则执行步骤(8);否则,算法终止,经如上过程得到的建模波长集合A即为融合后的特征波长;
(8)将rmax所对应于剩余波长集合B中的第max个波长加入建模波长集合A,从而更新建模波长集合A,并将第max个波长从剩余波长集合B中剔除,从而更新剩余波长集合B;令R2=rmax,清空数组r;
(9)判断R2>ΔR是否成立,若成立,算法终止,经如上过程得到的建模波长集合A即为融合后的特征波长;否则,重复步骤(3)顺序执行,直至剩余波长集合B为空。
如图4所示,基于融合特征波长选择方法提取的波长变量共27个,分别为484.806,488.580,490.047,535.957,575.845,576.048,580.917,600.102,600.304,651.909,652.107,660.379,708.872,711.378,715.228,743.535,760.966,772.829,796.935,818.215,862.894,868.124,873.521,873.700,893.730,908.099,921.143,单位为nm。
步骤六、以融合后的特征波长所对应的预处理后的光谱作为模型的输入,以酸度数据作为模型的输出,建立苹果酸度预测模型,从而用于苹果酸度的无损检测。具体实施中,该步骤分为模型的建立与评估两部分。
在校正集上建立偏最小二乘PLS预测模型;将预测集中特征波长所对应的预处理后的光谱输入该模型,得到预测集样本酸度的预测值,与测定值进行比较。评估方法为四个指标:预测集的平均绝对误差MAEP、均方根误差RMSEP、相对分析误差RPD和决定系数
Figure BDA0002608188500000055
其计算公式分别为:
Figure BDA0002608188500000051
Figure BDA0002608188500000052
其中,yi
Figure BDA0002608188500000053
分别表示第i个样本点酸度的测定值和预测值,
Figure BDA0002608188500000054
表示所有样本点酸度的平均值。一个性能良好的预测模型,MAEP及RMSEP应尽可能低,RPD尽可能高,
Figure BDA0002608188500000061
尽可能接近于1。
分别建立不同波长选择方法下的PLS预测模型,将模型结果与融合算法的结果进行比较,如表2所示。
表2不同波长选择方法下模型的预测结果
Figure BDA0002608188500000062
从表2中可以看出,基于提出的融合特征波长选择算法建立的苹果酸度预测模型,其MAEP为0.0169,RMSEP为0.0188,RPD为4.2248,
Figure BDA0002608188500000063
为0.9440,四个指标均仅次于CARS方法,较SPA方法有极大的提高;选取特征波长变量数目为27,较CARS方法大幅减少。基于融合特征波长选择算法建立的苹果酸度预测模型结果如图5所示。

Claims (1)

1.一种基于融合特征波长选择算法的苹果酸度近红外无损检测方法,其特征在于,步骤如下:
步骤一、采集苹果样本标记点区域的原始近红外反射率光谱;
步骤二、测量所述苹果样本标记点区域的酸度数据;
步骤三、对所述原始近红外反射率光谱进行预处理,得到预处理后的光谱;所述预处理包括:选取波长范围、将反射率光谱转换为吸光度光谱、对吸光度光谱进行归一化、平滑、散射校正的处理;
步骤四、利用连续投影算法与竞争自适应重加权采样算法分别对预处理后的光谱进行特征波长选择,并得到SPA特征波长集合S和CARS特征波长集合C;
步骤五、将SPA特征波长集合S和CARS特征波长集合C进行特征波长的融合,得到融合后的特征波长;
(1)令建模波长集合A=S,令剩余波长集合B=C–S,即集合C与集合S的差集;定义数组r,定义变量i,定义阈值ΔR;
(2)用所述建模波长集合A中波长对应的预处理后的光谱建立校正模型,并计算所述校正模型的决定系数R2
(3)令k为所述剩余波长集合B中元素的个数;初始化变量i=1;
(4)从所述剩余波长集合B中取出第i个波长加入所述建模波长集合A中,从而得到加入波长后的集合Ai
(5)用所述加入波长后的集合Ai中波长对应的预处理后的光谱建立第i个校正模型,并计算所述第i个校正模型的决定系数Ri 2后加入所述数组r中;
(6)将i+1复制给i,并判断i>k是否成立,若成立,则执行步骤(7),否则,返回步骤(4);
(7)取出所述数组r中的最大值,记为rmax,并判断rmax>R2是否成立,若成立,则执行步骤(8);否则,算法终止,从而得到建模波长集合A即为融合后的特征波长;
(8)将rmax所对应于剩余波长集合B中的第max个波长加入建模波长集合A,从而更新建模波长集合A,并将第max个波长从剩余波长集合B中剔除,从而更新剩余波长集合B;令R2=rmax,清空数组r;
(9)判断R2>ΔR是否成立,若成立,算法终止,从而得到建模波长集合A即为融合后的特征波长;否则,重复步骤(3)顺序执行,直至剩余波长集合B为空;
步骤六、以所述融合后的特征波长所对应的预处理后的光谱作为模型的输入,以所述酸度数据作为模型的输出,建立苹果酸度预测模型,从而用于苹果酸度的无损检测。
CN202010745439.2A 2020-07-29 2020-07-29 基于融合特征波长选择算法的苹果酸度近红外无损检测方法 Active CN111855608B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010745439.2A CN111855608B (zh) 2020-07-29 2020-07-29 基于融合特征波长选择算法的苹果酸度近红外无损检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010745439.2A CN111855608B (zh) 2020-07-29 2020-07-29 基于融合特征波长选择算法的苹果酸度近红外无损检测方法

Publications (2)

Publication Number Publication Date
CN111855608A CN111855608A (zh) 2020-10-30
CN111855608B true CN111855608B (zh) 2023-04-07

Family

ID=72946463

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010745439.2A Active CN111855608B (zh) 2020-07-29 2020-07-29 基于融合特征波长选择算法的苹果酸度近红外无损检测方法

Country Status (1)

Country Link
CN (1) CN111855608B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113533236B (zh) 2021-09-13 2021-11-30 广东省农业科学院农业质量标准与监测技术研究所 水果光谱的特征提取方法、模型训练方法、检测方法
CN114166764A (zh) * 2021-11-09 2022-03-11 中国农业科学院农产品加工研究所 基于特征波长筛选的光谱特征模型的构建方法及装置
CN114136887B (zh) * 2021-12-07 2023-10-20 广东省农业科学院蔬菜研究所 一种基于近红外光谱技术快速检测冬瓜口感决定因子苹果酸的方法
CN114441468B (zh) * 2022-02-10 2023-06-27 合肥工业大学 一种基于近红外-中红外融合算法的窖泥品质检测方法
CN116026780B (zh) * 2023-03-28 2023-07-14 江西中医药大学 基于串联策略波长选择的包衣吸湿率在线检测方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002022657A (ja) * 2000-07-05 2002-01-23 Takara Keiki Seisakusho:Kk 果実の糖酸度非破壊計測法および装置
CN106769970A (zh) * 2016-12-07 2017-05-31 甘肃萃英大农科技有限公司 一种快速无损检测苹果的品种、糖度和酸度的方法
CN108956545B (zh) * 2018-06-15 2020-10-27 北京农业智能装备技术研究中心 一种水果内部品质无损检测模型建立方法及系统
CN110264050B (zh) * 2019-06-06 2021-05-25 北京农业质量标准与检测技术研究中心 农产品品质分析方法及分析仪

Also Published As

Publication number Publication date
CN111855608A (zh) 2020-10-30

Similar Documents

Publication Publication Date Title
CN111855608B (zh) 基于融合特征波长选择算法的苹果酸度近红外无损检测方法
CN107703097B (zh) 利用近红外光谱仪构建快速预测原油性质的模型的方法
CN105486655B (zh) 基于红外光谱智能鉴定模型的土壤有机质快速检测方法
CN107748146A (zh) 一种基于近红外光谱检测的原油属性快速预测方法
CN111968080A (zh) 一种基于高光谱和深度学习的肥城桃内外品质检测方法
CN104897607A (zh) 便携式近红外光谱食品建模与快速检测一体化方法和系统
CN104965973B (zh) 一种苹果霉心病多因子无损检测判别模型及其建立方法
CN107817223A (zh) 快速无损实时预测原油性质模型的构建方法及其应用
CN113267466B (zh) 基于光谱波长优选的水果糖度和酸度无损检测方法
CN109409350B (zh) 一种基于pca建模反馈式载荷加权的波长选择方法
CN110779875B (zh) 一种基于高光谱技术检测冬小麦麦穗水分含量的方法
CN112098357B (zh) 一种基于近红外光谱的草莓感官品质等级评价方法
CN111537469A (zh) 一种基于近红外技术的苹果品质快速无损检测方法
CN116578851A (zh) 一种高光谱土壤有效硼含量预测方法
WO2023108514A1 (zh) 茶叶近红外光谱分析中谱峰自动检测与重构方法及系统
CN111999258A (zh) 一种面向光谱基线校正的加权建模局部优化方法
CN113655027A (zh) 一种近红外快速检测植物中单宁含量的方法
CN111141809B (zh) 一种基于非接触式电导信号的土壤养分离子含量检测方法
CN117538287A (zh) 一种无损检测黄冠梨果磷含量的方法及装置
CN116718613A (zh) 基于太赫兹谱的环氧玻璃纤维复合材料内部缺陷检测方法
CN109459408A (zh) 一种基于稀疏回归lar算法的近红外定量分析方法
CN115392305A (zh) 基于改进时间卷积网络的土壤有机质含量高光谱建模方法
CN114414524A (zh) 快速检测航空煤油性质的方法
CN112861411A (zh) 基于近红外光谱主成分和神经网络的生物质钠含量测量与建模方法
CN112229817A (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