CN103871077B - 一种道路车辆监控视频中的关键帧提取方法 - Google Patents

一种道路车辆监控视频中的关键帧提取方法 Download PDF

Info

Publication number
CN103871077B
CN103871077B CN201410080355.6A CN201410080355A CN103871077B CN 103871077 B CN103871077 B CN 103871077B CN 201410080355 A CN201410080355 A CN 201410080355A CN 103871077 B CN103871077 B CN 103871077B
Authority
CN
China
Prior art keywords
image
images
sequence
sample image
sample
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
CN201410080355.6A
Other languages
English (en)
Other versions
CN103871077A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201410080355.6A priority Critical patent/CN103871077B/zh
Publication of CN103871077A publication Critical patent/CN103871077A/zh
Application granted granted Critical
Publication of CN103871077B publication Critical patent/CN103871077B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)
  • Traffic Control Systems (AREA)
  • Closed-Circuit Television Systems (AREA)

Abstract

本发明公开了一种道路车辆监控视频中的关键帧提取方法,首先提取采集到的道路车辆监控视频的样本图像集,通过前景运动对象检测提取车辆最小外接矩形图像作为样本图像源,针对每个从进入到离开监控区域的车辆目标的外接矩形序列图像,判断符合车辆监控分析要求的作为正样本图像,不符合要求的作为负样本图像,得到用于训练AdaBoost分类器的训练样本图像集;提取前景车辆目标的面积特征和积分通道特征,训练得到AdaBoost分类器;利用训练好的分类器对从进入到离开监控区域的车辆目标的外接矩形图像序列进行分类,AdaBoost分类器对每帧图像打分,选择打分最高的图像作为关键帧。本发明可提高道路车辆监控视频关键帧提取的有效性和准确性。

Description

一种道路车辆监控视频中的关键帧提取方法
技术领域
本发明涉及数字图像处理领域,特别是一种道路车辆监控视频中的关键帧提取方法。
背景技术
随着数字视频处理技术的发展,各种道路视频监控系统广泛应用于道路安全管理。如何对海量道路监控视频数据进行有效处理是一个重要的问题。道路监控视频具有区别于一般视频的特点:摄像机固定于一个位置对道路场景进行拍摄,几乎不存在场景的变换;人们最关心的重要信息是从进入视频到离开监控区域的运动车辆的车牌、车型、车辆颜色等特征信息。
现有单纯使用MPEG4或H.264视频编码标准的道路视频监控系统存在监控视频数据量与视频检索效率的矛盾。视频数据量越大,检索效率越低。如何从车辆进入到离开监控范围的视频序列中自动得到最能清晰反映车辆车牌、车型等信息的显著关键帧是解决该矛盾的关键。道路车辆监控系统中一般都采用设定虚拟线圈的方法提取关键帧,即在道路监控视频中设定一系列的预定区域,模拟现实中的感应线圈,当进入虚拟线圈的运动车辆产生触发时提取相应的图像作为关键帧。该方法计算速度快,成本低,但虚拟线圈的位置、大小和数目都需要针对不同的道路监控视频具体而定,具有很大局限性。
针对这个问题,研究者们提出了各种改进的视频数据量压缩方法,目的是提取能够充分反映车牌、车型、车辆颜色等特征信息的关键帧。ZHAOShu-long在利用运动车辆检测对视频分段的基础上,提出一种基于车牌定位的关键帧提取方法。该方法选择车牌距离图像底部最近的图像帧作为关键帧,只适用于摄像机正对车道拍摄的情况。当监控摄像机位于路口侧面,与车道成一个角度时该方法并不适用。CongcongLi提出一种贪婪算法,先检测出运动目标,然后不断将检测出运动目标的图像帧进行融合重叠,得到虚拟的关键帧图像,该方法只适用于有固定速度运动目标的场景。YuanfengYang通过设置虚拟检测线,根据当前帧与背景的差异衡量当前帧的运动信息量,从而进行视频分段,并采用聚类的方法进行关键帧选取。YanYang等提取图像的边缘直方图和信息熵作为运动信息特征,选择对应的局部最大值作为关键帧,该方法简单直观,但没有考虑人对道路监控视频关注的重点,准确率较低。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种道路车辆监控视频中的关键帧提取方法,快速、准确地提取关键帧,从而实现道路车辆监控视频分析数据量的有效压缩。
为解决上述技术问题,本发明所采用的技术方案是:一种道路车辆监控视频中的关键帧提取方法,该方法为:
1)通过前景运动对象检测方法从采集到的道路车辆监控视频序列中提取前景车辆目标,解码所述采集到的道路车辆监控视频序列,得到多帧序列图像,从每一帧序列图像中截取前景车辆目标最小外接矩形图像,计算所有前景车辆目标最小外接矩形图像在其所对应的序列图像中的面积比例其中,分别为前景车辆目标最小外接矩形图像的像素高和宽;分别为序列图像的像素高和宽;得到面积特征向量;
2)将所有前景车辆目标最小外接矩形图像转换成像素大小一致的转换图像,对所述转换图像进行梯度幅值通道变换和梯度方向通道变换,得到积分通道特征向量;
3)扩充上述面积特征向量,使得面积特征向量与积分通道特征向量维数一致,合并维数一致的面积特征向量与积分通道特征向量,得到总的特征向量;
4)将上述所有前景车辆目标最小外接矩形图像作为样本图像源,从所述样本图像源中选取符合车辆监控分析要求的作为正样本图像,将样本图像源中不符合车辆监控分析要求的作为负样本图像;符合车辆监控分析要求是指样本图像在对应的序列图像中占的面积比例为0.4%以上,不符合车辆监控分析要求是指样本图像在对应的序列图像中占的面积比例小于0.4%;
5)标记所有的正样本图像和负样本图像:其中表示负样本图像,表示正样本图像,表示样本图像源中的样本图像;
6)初始化权重其中分别为负样本图像和正样本图像的个数;
7)训练分类器:令
8)利用下式归一化权重:为训练第个分类器的第个样本图像的权重;为归一化后的第个分类器的第个样本图像的权重;
9)即得到阈值和方向其中是总的特征向量中的元素最大值,弱分类器为:
h j = 1 p j f < p j &theta; j 0 else
其中决定不等式的方向,为1或者0;
则第个弱分类器的迭代误差和为:
10)选择迭代误差和最小的弱分类器作为候选分类器;
11)利用下式更新权重:被正确分类时,反之,
12)令作为训练第个分类器的第个样本图像的权重,重复上述步骤8)~步骤11),直到得到个候选分类器,利用所述个候选分类器确定强分类器
h ( x ) = 1 , &Sigma; t = 1 T &alpha; t h t ( x ) &GreaterEqual; 1 2 &Sigma; t = 1 T &alpha; t ; 0 , else
13)解码需要检测的道路车辆监控视频序列,得到多帧需检测的序列图像,利用上述强分类器对所述多帧需检测的序列图像打分,将得分最高的作为关键帧。
与现有技术相比,本发明所具有的有益效果为:本发明通过提取样本图像的面积特征和积分通道特征,利用AdaBoost方法训练得到一个分类器,最后利用AdaBoost分类器对测试样本进行分类,能快速、准确地提取出道路车辆监控视频中的关键帧,能够有效地提取出一个运动车辆从进入到离开监控区域的序列图像帧中最清晰显著的图像,从而实现道路车辆监控视频分析数据量的有效压缩。与现有的其他关键帧提取方法相比,在针对道路监控视频中特定车辆图像序列的关键帧提取问题上,本发明结合人对道路车辆监控视频关注重点,提取样本图像的面积特征和积分通道特征,表现出更高的准确率,更低的漏检率。同时,只要前期一次训练出有效的AdaBoost分类器,使用时提取关键帧速度快,方便快捷,具有较好的适用性。
附图说明
图1为本发明的方法流程图;
图2(a)正样本示例,车牌车型清晰,显著度高;
图2(b)负样本示例,因角度、面积等原因,车牌、车型等显著度低;
图3中间结果对比;
图4(a)测试视频示例;
图4(b)前景运动检测提取车辆目标最小外接矩形;
图4(c)分类器分类结果;
图4(d)选择分类结果中打分最高的图像作为关键帧。
具体实施方式
如图1,本发明首先对采集到的道路车辆监控视频提取训练样本图像集,通过前景运动对象检测提取车辆最小包围矩形图像作为样本图像源,针对每个从进入到离开监控区域的车辆目标的外接矩形序列图像,选取肉眼主观判断符合车辆监控分析要求的作为正样本图像,不符合要求的作为负样本图像,得到用于训练AdaBoost分类器的训练样本图像集,符合车辆监控分析要求是指可以清晰辨识车牌号码,且样本图像在对应的序列图像中占的面积比例为0.4%以上,不符合车辆监控分析要求是无法清晰辨识车牌号码,并且指样本图像在对应的序列图像中占的面积比例小于0.4%;如图2(a)和图2(b)所示为部分样本图像示意图;其次,先对所有样本提取面积特征,再将所有样本统一转换成64×64像素大小,提取积分通道特征,结合面积特征和积分通道特征,得到2560维的特征向量;再次,将特征向量作为AdaBoost的输入,设定训练轮数T=1024,训练得到分类器;最后,选取若干段监控视频对算法进行测试。本发明的具体方法步骤如下:
首先,对采集到的道路车辆监控视频提取训练样本图像集。在长沙三一大道、双拥路和车站北路等路口采集上午和下午真实路况下具有代表性的道路车辆视频,共采集48段视频,总时长52分钟,总数据量约为2.1G,包含208个运动车辆,最长的运动车辆图像序列为66帧,最短的为23帧。从视频序列中截取包含运动车辆的最小外接矩形图像共2677张,选取的正样本图像共683张,负样本图像共1994张,用于训练分类器,部分样本如图2(a)和图2(b)所示。
其次,对所有样本提取面积特征和积分通道特征。先提取面积特征,采用以下公式计算前景车辆目标的面积比例:
P = h &times; w H &times; W
h,w分别为前景车辆目标最小外接矩形的像素高和宽;
H,W分别为采集到的视频解码后得到序列图像的像素高和宽。
将所有样本统一转换成64×64像素大小,提取积分通道特征。为实现特征的快速计算,将特征定义为通道图上随机矩形区域的像素值和,可以通过通道积分图在常数复杂度下完成计算,即:
ck表示第K类积分通道。
随机选取通道图上两个像素点,像素坐标分别为(x1,y1),(x2,y2),得到四个随机矩形区域,矩形的四个顶点像素坐标分别为(0,0)(0,y1)(x1,0)(x1,y1),(0,0)(0,y1)(x2,0)(x2,y1),(0,0)(0,y2)(x1,0)(x1,y2),(0,0)(0,y2)(x2,0)(x2,y2)fk(x1,y1)、fk(x1,y2)、fk(x2,y1)和fk(x2,y2)分别表示相应随机矩形区域所有点的像素值和。
计算七类积分通道(一个梯度幅值通道、六个梯度方向通道)。
(1)梯度幅值通道:梯度幅值通道C1的计算基于图像I的灰度图Igrey,计算方法如下:
分别表示图像I的灰度图Igray对x,y所求偏导数。
(2)梯度方向通道:首先要基于灰度图Igray计算每个像素位置的梯度方向:
然后把θ(x,y)离散为六个梯度方向得到六个方向上的梯度幅值通道图:Ck(x,y)=C1(x,y)·1[θ(x,y)=θk-2],其中k=2,3,...,7。
结合面积特征和积分通道特征,得到2560维的特征向量。
再次,将特征向量作为AdaBoost的输入,设定训练轮数T=1024,训练得到强分类器。具体的AdaBoost算法过程如下:
步骤(1):标记n个样本图像(x1,y1),...,(xnyn),其中m个负正样本图像标记为yi=0,、l个正样本图像标记为yi=1。
步骤(2):初始化权重。每个负样本图像的初始权值为每个正样本图像的初始权值为
步骤(3):选择T个弱分类器(迭代T次)。对于t=1,...,T,
权值归一化,使ωt成概率分布。
&omega; t , i &prime; &LeftArrow; &omega; t , i &Sigma; i = 1 n &omega; t , i
对每一个特征j,训练一个弱分类器hj,则第j个弱分类器的迭代误差和 &epsiv; j = &Sigma; i &omega; t , i &prime; | h j ( x i ) - y i | ;
选择最小迭代误差和εt对应的弱分类器ht
更新权重:若样本xi被正确分类(即当标记样本xi为正样本,之后又将xi分类为正样本时,分类正确;反之,不正确),ei=0;否则,ei=1。
步骤(4):将T个弱分类器线性组合在一起,得到一个强分类器:
h ( x ) = 1 &Sigma; t - 1 &tau; &alpha; t h t ( x ) &GreaterEqual; 1 2 &Sigma; t - 1 &tau; &alpha; t 0 otherwise
其中
最后,选取若干段监控视频对算法进行测试。测试视频共包含67个运动车辆目标,从中截取了运动车辆最小外接矩形图像2000张,有67个从进入到离开监控区域的车辆目标的外接矩形图像序列。以人工目测法作为评价其他算法参考准则,通过比较各种算法提取出的图像帧区间分布及其与人工目测法提取的图像帧区间的重合程度来比较算法中间结果的优劣。人工目测方法从肉眼主观判断符合车辆监控分析要求的图像帧区间中任选一帧作为关键帧,K-Means聚类方法和本文算法提取的关键帧落在肉眼主观判断符合车辆监控分析要求的图像帧区间时,即认为与人工目测法提取的关键帧相同。
通过提取关键帧的准确率、误检率和漏检率指标来衡量关键提取算法的有效性。各指标定义如下:
以测试视频中提取的车牌号为湘AKL489的运动车辆图像序列为例,如图4(a)-图4(d)所示,在监控区域共记录53帧,肉眼主观判断符合车辆监控分析要求的图像帧区间是[47,53],如图3所示。从图3可以看出,本文算法提取第50、51、52、53帧,均在肉眼主观判断符合车辆监控分析要求的图像帧区间[47,53]内,且分布集中;K-Means聚类方法提取第11、12、44、45、46、47、48、53帧,其中第47、48和53帧在肉眼主观判断符合车辆监控分析要求的图像帧区间[47,53]内,但重合度小,且分布不集中,呈现跳跃性。本文算法在中间结果的提取上表现出更优越的性能。
关键帧提取的实验结果如表1所示,可以看出,与K-Means聚类方法相比,本文方法表现好出更高的准确率、更低的漏检率和误检率,所提取的结果与人工目测法提取的结果最为相近。
通过以上本发明所提出的针对道路监控视频中特定车辆图像序列,在前景运动对象检测的基础上的关键帧提取方法,能够有效地提取出一个运动车辆从进入到离开监控区域的序列图像帧中最清晰显著的图像,从而实现道路车辆监控视频分析数据量的有效压缩。与传统聚类的关键帧提取算法相比,本发明的方法可提高道路车辆监控视频关键帧提取的有效性和准确性。

Claims (1)

1.一种道路车辆监控视频中的关键帧提取方法,其特征在于,该方法为:
1)通过前景运动对象检测方法从采集到的道路车辆监控视频序列中提取前景车辆目标,解码所述采集到的道路车辆监控视频序列,得到多帧序列图像,从每一帧序列图像中截取前景车辆目标最小外接矩形图像,计算所有前景车辆目标最小外接矩形图像在其所对应的序列图像中的面积比例P:其中,分别为前景车辆目标最小外接矩形图像的像素高和宽;分别为序列图像的像素高和宽;得到面积特征向量;
2)将所有前景车辆目标最小外接矩形图像转换成像素大小一致的转换图像,对所述转换图像进行梯度幅值通道变换和梯度方向通道变换,得到积分通道特征向量;
3)扩充上述面积特征向量,使得面积特征向量与积分通道特征向量维数一致,合并维数一致的面积特征向量与积分通道特征向量,得到总的特征向量;
4)将上述所有前景车辆目标最小外接矩形图像作为样本图像源,从所述样本图像源中选取肉眼主观判断符合车辆监控分析要求的作为正样本图像,将样本图像源中不符合车辆监控分析要求的作为负样本图像;符合车辆监控分析要求是指能清晰辨识车牌号码、且样本图像在对应的序列图像中占的面积比例为0.4%以上,不符合车辆监控分析要求是指无法清晰辨识车牌号码、且样本图像在对应的序列图像中占的面积比例小于0.4%;
5)标记所有的正样本图像和负样本图像:其中表示负样本图像,表示正样本图像,表示样本图像源中的样本图像;
6)初始化权重 &omega; 1 , i : &omega; 1 , i = 1 2 m , y i = 0 1 2 l , y i = 1 ; 其中分别为负样本图像和正样本图像的个数; n = m + l ;
7)训练分类器:令t=1;
8)利用下式归一化权重: 为训练第个分类器的第个样本图像的权重;为归一化后的第个分类器的第个样本图像的权重;
9)对所述总的特征向量中的每个元素训练一个弱分类器即得到阈值和方向其中 是总的特征向量中的元素最大值,所述弱分类器为:
h j = 1 p j f < p j &theta; j 0 else
其中决定不等式的方向,为1或者0;
则第个弱分类器的迭代误差和为:
10)选择迭代误差和最小的弱分类器作为候选分类器;
11)利用下式更新权重: &omega; t + 1 , i = &omega; &prime; t , i &beta; t 1 - e i , 被正确分类时, e i = 0 , 反之, e i = 1 ; &beta; t = &epsiv; t 1 - &epsiv; t ; &alpha; t = log 1 &beta; t ; 为最小的迭代误差和;
12)令t=t+1,将作为训练第个分类器的第个样本图像的权重,重复上述步骤8)~步骤11),直到得到T个候选分类器,利用所述T个候选分类器确定强分类器
h ( x ) = 1 , &Sigma; t = 1 T &alpha; t h t ( x ) &GreaterEqual; 1 2 &Sigma; t = 1 T &alpha; t 0 , else ;
13)解码需要检测的道路车辆监控视频序列,得到多帧需检测的序列图像,利用上述强分类器对所述多帧需检测的序列图像打分,将得分最高的作为关键帧。
CN201410080355.6A 2014-03-06 2014-03-06 一种道路车辆监控视频中的关键帧提取方法 Active CN103871077B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410080355.6A CN103871077B (zh) 2014-03-06 2014-03-06 一种道路车辆监控视频中的关键帧提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410080355.6A CN103871077B (zh) 2014-03-06 2014-03-06 一种道路车辆监控视频中的关键帧提取方法

Publications (2)

Publication Number Publication Date
CN103871077A CN103871077A (zh) 2014-06-18
CN103871077B true CN103871077B (zh) 2016-06-15

Family

ID=50909582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410080355.6A Active CN103871077B (zh) 2014-03-06 2014-03-06 一种道路车辆监控视频中的关键帧提取方法

Country Status (1)

Country Link
CN (1) CN103871077B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104504732A (zh) * 2014-12-25 2015-04-08 合肥寰景信息技术有限公司 一种基于关键帧提取的视频内容检索方法
CN104537392B (zh) * 2014-12-26 2017-10-17 电子科技大学 一种基于判别性语义部件学习的对象检测方法
CN105404901B (zh) * 2015-12-24 2019-10-18 上海玮舟微电子科技有限公司 分类器的训练方法、图像检测方法及各自系统
CN106250912A (zh) * 2016-07-21 2016-12-21 成都之达科技有限公司 基于图像的车辆位置获取方法
CN106529522A (zh) * 2016-11-07 2017-03-22 湖南源信光电科技有限公司 一种基于积分通道特征的车牌定位方法
CN106845520B (zh) * 2016-12-23 2018-05-18 深圳云天励飞技术有限公司 一种图像处理方法及终端
CN106991668B (zh) * 2017-03-09 2020-08-18 南京邮电大学 一种天网摄像头拍摄画面的评价方法
CN109993035B (zh) * 2017-12-29 2021-06-29 深圳市优必选科技有限公司 基于嵌入式系统的人体检测的方法以及装置
CN109086682B (zh) * 2018-07-11 2021-07-27 东南大学 一种基于多特征融合的智能视频黑烟车检测方法
CN110992420B (zh) * 2019-12-18 2023-07-21 天地伟业技术有限公司 一种检测粮仓车辆运载情况的系统及方法
CN113053136A (zh) * 2019-12-26 2021-06-29 上海晋沙智能科技有限公司 一种道路智能安防监控系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102184388A (zh) * 2011-05-16 2011-09-14 苏州两江科技有限公司 人脸和车辆自适应快速检测系统及检测方法
CN102855500A (zh) * 2011-06-27 2013-01-02 东南大学 一种基于Haar和HoG特征的前车检测方法
CN103514460A (zh) * 2013-07-30 2014-01-15 深圳市智美达科技有限公司 视频监控多视角车辆检测方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102184388A (zh) * 2011-05-16 2011-09-14 苏州两江科技有限公司 人脸和车辆自适应快速检测系统及检测方法
CN102855500A (zh) * 2011-06-27 2013-01-02 东南大学 一种基于Haar和HoG特征的前车检测方法
CN103514460A (zh) * 2013-07-30 2014-01-15 深圳市智美达科技有限公司 视频监控多视角车辆检测方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
车牌识别系统中关键技术的研究;许利显;《<中国优秀硕士学位论文全文数据库>电子期刊》;20081231(第8期);第17-23,28-32页 *

Also Published As

Publication number Publication date
CN103871077A (zh) 2014-06-18

Similar Documents

Publication Publication Date Title
CN103871077B (zh) 一种道路车辆监控视频中的关键帧提取方法
CN106469299B (zh) 一种车辆搜索方法及装置
CN108921083B (zh) 基于深度学习目标检测的非法流动摊贩识别方法
CN104463196B (zh) 一种基于视频的天气现象识别方法
CN103530600B (zh) 复杂光照下的车牌识别方法及系统
CN100545867C (zh) 航拍交通视频车辆快速检测方法
CN103530638B (zh) 多摄像头下的行人匹配方法
Derpanis et al. Classification of traffic video based on a spatiotemporal orientation analysis
CN107273832B (zh) 基于积分通道特征与卷积神经网络的车牌识别方法及系统
CN103020985B (zh) 一种基于场量分析的视频图像显著性检测方法
CN102214309B (zh) 一种基于头肩模型的特定人体识别方法
CN104700099A (zh) 识别交通标志的方法和装置
CN105608456A (zh) 一种基于全卷积网络的多方向文本检测方法
CN111582339B (zh) 一种基于深度学习的车辆检测与识别的方法
CN106295532B (zh) 一种视频图像中的人体动作识别方法
CN103810707B (zh) 一种基于移动视觉焦点的图像视觉显著性检测方法
CN104978567A (zh) 基于场景分类的车辆检测方法
Juranek et al. Real-time pose estimation piggybacked on object detection
CN112016605A (zh) 一种基于边界框角点对齐和边界匹配的目标检测方法
CN114067444A (zh) 基于元伪标签和光照不变特征的人脸欺骗检测方法和系统
CN103853724A (zh) 多媒体数据分类方法及装置
CN102867183A (zh) 一种车辆遗撒物检测方法、装置及智能交通监控系统
CN111008574A (zh) 一种基于形体识别技术的重点人员轨迹分析方法
CN106557740A (zh) 一种遥感图像中油库目标的识别方法
CN104732534A (zh) 一种图像中显著目标的抠取方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant