CN102360526A - 一种快速道路路段状态的实时监测方法 - Google Patents

一种快速道路路段状态的实时监测方法 Download PDF

Info

Publication number
CN102360526A
CN102360526A CN2011103023868A CN201110302386A CN102360526A CN 102360526 A CN102360526 A CN 102360526A CN 2011103023868 A CN2011103023868 A CN 2011103023868A CN 201110302386 A CN201110302386 A CN 201110302386A CN 102360526 A CN102360526 A CN 102360526A
Authority
CN
China
Prior art keywords
traffic flow
traffic
accident
flow modes
check point
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
CN2011103023868A
Other languages
English (en)
Other versions
CN102360526B (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN 201110302386 priority Critical patent/CN102360526B/zh
Publication of CN102360526A publication Critical patent/CN102360526A/zh
Application granted granted Critical
Publication of CN102360526B publication Critical patent/CN102360526B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)

Abstract

一种快速道路路段状态的实时监测方法,在待检查的快速道路上设置交通流检测点,各个检测点之间布设的间距为800米,每个检测点覆盖上游400米到下游400米范围的路段,收集历史交通流数据和事故数据,通过聚类分析和回归分析提取危险交通流运行状态,并通过决策树方法给出交通流运行状态的监测判别方法。本发明方法处理简单易行,便于实时鉴别快速道路事故“黑点”,使用本发明方法实时监测快速道路事故危险路段具有实际的工程运用价值。

Description

一种快速道路路段状态的实时监测方法
技术领域
本发明属于交通智能管理和控制技术领域,为一种快速道路路段状态的实时监测方法,可实时鉴别快速道路上有发生交通事故风险的路段。
背景技术
随着改革的深化和城市化、机动化进程的不断深入,自20世纪80年代中期开始,我国交通运输业迅速发展,路网结构得到不断完善,交通基础设施条件不断改善,机动车保有量和交通需求迅速增加。但是,车辆与道路增长速度的不协调以及管理手段的不当,也带来了严重的交通安全问题。在这样的背景下,道路交通安全问题成为越来越受到关注的热点问题,提高我国交通安全管理工作水平,对于减少人员伤亡和财产损失,发展新时期和谐交通具有重要的现实意义。
为了预防和降低交通事故数量,如何鉴别事故高发路段是首要解决的问题。目前,交通管理相关部门主要通过一些传统的方法来实现对事故黑点的鉴别,这些方法往往需要采集较长一段时间内的事故数据、交通数据、天气数据和道路参数,一般国际上推荐采用3年的事故数据进行分析。虽然这些方法在鉴别快速道路事故黑点方面是有效的,但是这些方法难以在短时间内对事故易发路段进行鉴别。
而另一方面,在交通安全前沿研究中,学者们开始关注实时交通流信息与交通安全之间的关系。目前,交通流的不同运行状态是否与交通安全相关,它们之间关联度究竟多大,事故发生前是否存在某种可作为事故征兆的“危险”交通流状态,一直没有得到很好的解答。本发明提出一种提取快速道路“危险”交通流状态的方法,并将该方法应用于快速道路事故危险路段的实时鉴别中,通过检测快速道路上各个检测点的交通流状态,实时鉴别事故危险路段,从而解决已有事故黑点鉴别方法难以在短时间内鉴别危险路段的缺陷。
发明内容
本发明要解决的问题是:克服已有事故黑点鉴别方法的不足,利用现有智能交通***检测的实时交通流数据,提出一种能够实时判别快速道路危险路段的方法。
本发明的技术方案为:快速道路事故路段状态的实时监测方法,包括以下步骤:
1)在快速道路上每隔800米设置一个交通流检测点,收集待检查的快速道路的事故数据,确立各起交通事故的发生地点,从而确定与交通事故发生地点距离最近的一个交通流检测点,提取该检测点在事故发生前半小时内的交通数据,包括交通流率和速度,其采样精度为每5分钟采集一次;
2)对于每起交通事故,采用病例-对照研究法按1∶4比例选取事发路段正常交通流状态下,即没有发生事故时的交通占有率数据,所述1∶4指,对应于每个事故范例,选取4个相应的正常交通流状态范例;
3)将事故发生前的交通数据和步骤2)选取的对应的正常交通流状态数据以及周围发生事故最少的一个检测点2个月内交通数据组合成交通流数据样本,采用K-means动态聚类分析方法,选用交通流率和速度2个交通流参数对交通流数据样本进行聚类分析;
4)对由K-means动态聚类分析方法划分得到的各交通流状态进行Logistic回归,自变量为交通流状态,应变量为是否发生交通事故,根据Logistic回归中各交通流状态前面系数的显著水平P-value是否显著来确定交通流状态是否为危险交通流运行状态,当显著水平P-value小于0.05时,代表该交通流状态为危险交通流状态,反之,当显著水平P-value大于0.05时,代表该交通流状态不是危险交通流状态;
5)采用QUEST决策树模型确立各种交通流状态的边界用于判别实时交通流状态:以交通密度作为QUEST决策树模型中的解释变量,以交通流状态为目标变量,确立各种交通流状态的边界,即各种交通流状态对应的交通密度范围,其中交通密度=交通流率/速度,交通流率与速度由检测点获得;
6)每隔5分钟采集各个检测点的交通流率和速度,计算得到各个点的交通密度,根据QUEST决策树模型确立的各种交通流状态对应的交通密度范围,判断整条快速道路是否有检测点出现危险交通流状态,如果某个检测点得到的交通密度落在危险交通流状态对应的交通密度范围,说明该检测点附近出现危险交通流状态,则表明该检测点覆盖的800米路段为危险路段,有发生交通事故的危险性,进行预警提示。
本发明利用快速道路上各个检测点采集的交通数据,实时检测各个检测点附近是否出现危险交通流状态,从而达到实时鉴别事故易发路段的目的。相对于传统快速道路事故黑点鉴别方法而言,本方法能够在较短时间内发现事故危险路段,并且处理简单易行。因而,使用本方法鉴别快速道路危险交通流运行状态具有实际的工程运用价值。
附图说明
图1为本发明实施例的检测点示意图。
图2为本发明监测方法的流程图。
图3为本发明综合运用K-means聚类、Logistic回归和QUEST决策树模型的建模程序框图。
图4为本发明实施例的聚类后交通流状态划分图。
具体实施方式
本发明将K-means聚类、Logistic回归和QUEST决策树模型用于实时监测事故易发路段,一种快速道路事故危险路段的实时监测方法,包括以下步骤:
1)在快速道路上每隔800米设置一个交通流检测点,收集待检查的快速道路的事故数据,确立各起交通事故的发生地点,从而确定与交通事故发生地点距离最近的一个交通流检测点,提取该检测点在事故发生前半小时内的交通数据,包括交通流率和速度,其采样精度为每5分钟采集一次;
2)对于每起交通事故,采用病例-对照研究法按1∶4比例选取事发路段正常交通流状态下,即没有发生事故时的交通占有率数据,所述1∶4指,对应于每个事故范例,选取4个相应的正常交通流状态范例;
3)将事故发生前的交通数据和步骤2)选取的对应的正常交通流状态数据以及周围发生事故最少的一个检测点2个月内交通数据组合成交通流数据样本,采用K-means动态聚类分析方法,选用交通流率和速度2个交通流参数对交通流数据样本进行聚类分析;
4)对由K-means动态聚类分析方法划分得到的各交通流状态进行Logistic回归,自变量为交通流状态,应变量为是否发生交通事故,根据Logistic回归中各交通流状态前面系数的显著水平P-value是否显著来确定交通流状态是否为危险交通流运行状态,当显著水平P-value小于0.05时,代表该交通流状态为危险交通流状态,反之,当显著水平P-value大于0.05时,代表该交通流状态不是危险交通流状态;
5)由于聚类分析方法是对收集的历史数据进行分类,只能给出各个状态的聚类中心,无法判别当前交通流究竟处于何种状态下,因而本发明采用QUEST决策树模型确立各种交通流状态的边界,用于判别实时交通流状态:以交通密度作为QUEST决策树模型中的解释变量,以交通流状态为目标变量,确立各种交通流状态的边界,即各种交通流状态对应的交通密度范围,其中交通密度=交通流率/速度,交通流率与速度由检测点获得;
6)每隔5分钟采集各个检测点的交通流率和速度,计算得到各个点的交通密度,根据QUEST决策树模型确立的各种交通流状态对应的交通密度范围,判断整条快速道路是否有检测点出现危险交通流状态,如果某个检测点得到的交通密度落在危险交通流状态对应的交通密度范围,说明该检测点附近出现危险交通流状态,则表明该检测点覆盖的800米路段为危险路段,有发生交通事故的危险性,进行预警提示。
步骤1)-4)将交通流划分不同状态,并利用P-value从中确立危险交通流状态,完成对交通流状态的划分;步骤5)确立危险交通流状态对应的交通密度范围,当当前交通密度落在危险交通流状态对应的交通密度范围时,即判定当前交通流为危险状态。
下面用结合附图说明本发明。在快速道路待检测区域上下游设置检测点,在检测点布置检测设备,例如电磁感应线圈或者其他交通流检测设备,各检测点之间相隔800米,如图1所示。
检测点以5分钟为一个采样单位采集交通流率和速度。令X为上述5分钟内交通流参数交通流率和速度,y为是否发生交通事故,发生交通事故时y取值为1,没有发生交通事故时y取值为-1,则,
Figure BDA0000096592430000041
Figure BDA0000096592430000042
Y = y 1 y 2 . . . y n = 1 / - 1 1 / - 1 . . . 1 / - 1
实际运用过程分为模型标定和模型使用两个过程。
模型标定:收集或采集上述各检测点一段时间内的交通数据,含事故样本和非事故样本。为了保证模型的准确性,保证模型能够反映交通流检测参数与交通事故发生之间的关系,样本应该足够大,根据已有研究事故范例不应该小于50个。根据前述步骤1)~5),计算得到交通流状态划分,并在其中确立危险交通流状态,由分类树最终得到危险交通流状态的边界值。
模型使用:实时采集各个检测点5分钟交通流率和速度,将各个点检测到的交通流率和速度与模型标定中确立的危险交通流状态边界值进行比较,当发现某点检测到的交通流率和速度处于危险交通流状态所在范围时,表明该点覆盖的800米路段有发生交通事故的危险性。此时应该通过快速道路可变信息显示板向驾驶员发布该点覆盖的路段有发生交通事故的危险性,提醒驾驶员谨慎驾驶车辆。
本实施例用加利福尼亚旧金山海湾地区的I-880高速公路上采集的真实交通数据和事故数据(I-880数据)来展示本方法的工作流程。数据包括待检查点的交通流率和速度,每隔5分钟采集一次。
在I-880数据库中总共提取了80个交通事故范例,选取了80起交通事故发生的交通流数据(包括交通流率和速度),同时利用病例-对照研究法按1∶4比例选取对应于各起交通事故的正常交通流状态下(没有发生交通事故)的交通数据。将事故组数据和正常交通流数据以及一个检测点2个月内的交通流数据组合成交通流数据样本。
利用运用K-means动态聚类法首先将交通流状态分为如图4中的6个状态,并利用Logistic回归发现状态5和状态6为交通事故发生前的危险交通流状态。
在确立危险交通流状态后,利用QUEST决策树模型确立6类交通流状态的边界值如表1:
表1 6类交通流状态的边界值
 交通流状态   交通密度(交通流率除以速度)边界值
 交通流状态1   0<交通密度≤17.6
 交通流状态2   17.6<交通密度≤27.2
 交通流状态3   27.2<交通密度≤36.8
 交通流状态4   36.8<交通密度≤56.0
 交通流状态5(危险交通流状态)   56.0<交通密度≤76.8
 交通流状态6(危险交通流状态)   交通密度≥76.8
在确立各交通流状态的边界值之后,每隔5分钟采集各个检测点的交通流率和速度,将各个点检测到的实时交通流数据与危险交通流状态的边界值进行比较,当发现某点的检测到的实时交通流数据处于交通流状态5或者6(危险交通流状态)所在范围时,表明该点覆盖的800米路段有发生交通事故的危险性。此时应该通过高速公路可变显示板向驾驶员发布该点有发生交通事故的危险性,提醒驾驶员谨慎驾驶车辆。
由于本方法处理简单易行,不受人为因素影响,便于实时鉴别快速道路事故危险路段;因此,使用本方法监测快速道路事故危险路段具有实际的工程运用价值。

Claims (1)

1.一种快速道路路段状态的实时监测方法,其特征是包括以下步骤:
1)在快速道路上每隔800米设置一个交通流检测点,收集待检查的快速道路的事故数据,确立各起交通事故的发生地点,从而确定与交通事故发生地点距离最近的一个交通流检测点,提取该检测点在事故发生前半小时内的交通数据,包括交通流率和速度,其采样精度为每5分钟采集一次;
2)对于每起交通事故,采用病例-对照研究法按1∶4比例选取事发路段正常交通流状态下,即没有发生事故时的交通占有率数据,所述1∶4指,对应于每个事故范例,选取4个相应的正常交通流状态范例;
3)将事故发生前的交通数据和步骤2)选取的对应的正常交通流状态数据,以及周围发生事故最少的一个检测点2个月内交通数据组合成交通流数据样本,采用K-means动态聚类分析方法,选用交通流率和速度2个交通流参数对交通流数据样本进行聚类分析;
4)对由K-means动态聚类分析方法划分得到的各交通流状态进行Logistic回归,自变量为交通流状态,应变量为是否发生交通事故,根据Logistic回归中各交通流状态前面系数的显著水平P-value是否显著来确定交通流状态是否为危险交通流运行状态,当显著水平P-value小于0.05时,代表该交通流状态为危险交通流状态,反之,当显著水平P-value大于0.05时,代表该交通流状态不是危险交通流状态;
5)采用QUEST决策树模型确立各种交通流状态的边界用于判别实时交通流状态:以交通密度作为QUEST决策树模型中的解释变量,以交通流状态为目标变量,确立各种交通流状态的边界,即各种交通流状态对应的交通密度范围,其中交通密度=交通流率/速度,交通流率与速度由检测点获得;
6)每隔5分钟采集各个检测点的交通流率和速度,计算得到各个点的交通密度,根据QUEST决策树模型确立的各种交通流状态对应的交通密度范围,判断整条快速道路是否有检测点出现危险交通流状态,如果某个检测点得到的交通密度落在危险交通流状态对应的交通密度范围,说明该检测点附近出现危险交通流状态,则表明该检测点覆盖的800米路段为危险路段,有发生交通事故的危险性,进行预警提示。
CN 201110302386 2011-09-28 2011-09-28 一种快速道路路段状态的实时监测方法 Active CN102360526B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110302386 CN102360526B (zh) 2011-09-28 2011-09-28 一种快速道路路段状态的实时监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110302386 CN102360526B (zh) 2011-09-28 2011-09-28 一种快速道路路段状态的实时监测方法

Publications (2)

Publication Number Publication Date
CN102360526A true CN102360526A (zh) 2012-02-22
CN102360526B CN102360526B (zh) 2013-09-18

Family

ID=45585849

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110302386 Active CN102360526B (zh) 2011-09-28 2011-09-28 一种快速道路路段状态的实时监测方法

Country Status (1)

Country Link
CN (1) CN102360526B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102610098A (zh) * 2012-03-09 2012-07-25 西安费斯达自动化工程有限公司 一种基于离散模型的路段交通流检测器设置方法
CN103116979A (zh) * 2013-01-17 2013-05-22 东南大学 基于***安全指标分布法的道路事故黑点鉴别***
CN103198713A (zh) * 2013-03-21 2013-07-10 东南大学 基于交通数据和天气数据的减少交通事故的车辆调控方法
CN103198707A (zh) * 2013-03-21 2013-07-10 东南大学 一种基于晴天状况下危险交通流特征的车辆调控方法
CN104182453A (zh) * 2014-06-20 2014-12-03 银江股份有限公司 一种面向海量历史浮动车数据的分布式地图匹配方法
CN104408918A (zh) * 2014-11-20 2015-03-11 江苏物联网研究发展中心 基于MapReduce计算框架的浮动车数据与道路匹配的方法
CN105046952A (zh) * 2015-06-17 2015-11-11 天津职业技术师范大学 稀疏道路交通事件地空联合检测方法
CN105869398A (zh) * 2016-05-10 2016-08-17 浙江工业大学 一种基于K-means聚类的道路交通畅通度判断方法
CN106485919A (zh) * 2016-10-11 2017-03-08 东南大学 一种判定快速道路定点测速仪对交通事故数量影响的方法
CN106781474A (zh) * 2016-12-21 2017-05-31 东南大学 一种基于交通视频判断交通事故发生地自愈能力的方法
CN106935052A (zh) * 2015-12-30 2017-07-07 沈阳美行科技有限公司 一种基于行驶数据的变速提示方法及装置
CN107590766A (zh) * 2017-09-28 2018-01-16 东南大学 一种与道路交通事故风险相关的用地组合形态的判别方法
CN108417091A (zh) * 2018-05-10 2018-08-17 武汉理工大学 基于网联车的行车风险路段辨识与预警***与方法
CN108682149A (zh) * 2018-05-21 2018-10-19 东南大学 基于二项Logistic回归的公路事故黑点路段线形致因分析方法
CN115188158A (zh) * 2022-06-22 2022-10-14 中国电子科技集团公司电子科学研究院 一种基于历史案例空间分布的智能安全预警方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6177885B1 (en) * 1998-11-03 2001-01-23 Esco Electronics, Inc. System and method for detecting traffic anomalies
CN101706996A (zh) * 2009-11-12 2010-05-12 北京交通大学 一种基于信息融合的高速公路交通状态识别方法
CN101739815A (zh) * 2009-11-06 2010-06-16 吉林大学 道路交通拥挤状态在线识别方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6177885B1 (en) * 1998-11-03 2001-01-23 Esco Electronics, Inc. System and method for detecting traffic anomalies
CN101739815A (zh) * 2009-11-06 2010-06-16 吉林大学 道路交通拥挤状态在线识别方法
CN101706996A (zh) * 2009-11-12 2010-05-12 北京交通大学 一种基于信息融合的高速公路交通状态识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张鹏: "交通事故黑点智能排查***的研究与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102610098B (zh) * 2012-03-09 2014-05-07 西安费斯达自动化工程有限公司 一种基于离散模型的路段交通流检测器设置方法
CN102610098A (zh) * 2012-03-09 2012-07-25 西安费斯达自动化工程有限公司 一种基于离散模型的路段交通流检测器设置方法
CN103116979A (zh) * 2013-01-17 2013-05-22 东南大学 基于***安全指标分布法的道路事故黑点鉴别***
CN103116979B (zh) * 2013-01-17 2015-01-28 东南大学 基于***安全指标分布法的道路事故黑点鉴别***
CN103198713B (zh) * 2013-03-21 2016-01-06 东南大学 基于交通数据和天气数据的减少交通事故的车辆调控方法
CN103198713A (zh) * 2013-03-21 2013-07-10 东南大学 基于交通数据和天气数据的减少交通事故的车辆调控方法
CN103198707A (zh) * 2013-03-21 2013-07-10 东南大学 一种基于晴天状况下危险交通流特征的车辆调控方法
CN103198707B (zh) * 2013-03-21 2016-03-02 东南大学 一种基于晴天状况下危险交通流特征的车辆调控方法
CN104182453A (zh) * 2014-06-20 2014-12-03 银江股份有限公司 一种面向海量历史浮动车数据的分布式地图匹配方法
CN104408918A (zh) * 2014-11-20 2015-03-11 江苏物联网研究发展中心 基于MapReduce计算框架的浮动车数据与道路匹配的方法
CN105046952A (zh) * 2015-06-17 2015-11-11 天津职业技术师范大学 稀疏道路交通事件地空联合检测方法
CN105046952B (zh) * 2015-06-17 2017-12-08 天津职业技术师范大学 稀疏道路交通事件地空联合检测方法
CN106935052A (zh) * 2015-12-30 2017-07-07 沈阳美行科技有限公司 一种基于行驶数据的变速提示方法及装置
CN105869398A (zh) * 2016-05-10 2016-08-17 浙江工业大学 一种基于K-means聚类的道路交通畅通度判断方法
CN106485919A (zh) * 2016-10-11 2017-03-08 东南大学 一种判定快速道路定点测速仪对交通事故数量影响的方法
CN106781474A (zh) * 2016-12-21 2017-05-31 东南大学 一种基于交通视频判断交通事故发生地自愈能力的方法
CN106781474B (zh) * 2016-12-21 2019-03-12 东南大学 一种基于交通视频判断交通事故发生地自愈能力的方法
CN107590766A (zh) * 2017-09-28 2018-01-16 东南大学 一种与道路交通事故风险相关的用地组合形态的判别方法
CN108417091A (zh) * 2018-05-10 2018-08-17 武汉理工大学 基于网联车的行车风险路段辨识与预警***与方法
CN108682149A (zh) * 2018-05-21 2018-10-19 东南大学 基于二项Logistic回归的公路事故黑点路段线形致因分析方法
CN115188158A (zh) * 2022-06-22 2022-10-14 中国电子科技集团公司电子科学研究院 一种基于历史案例空间分布的智能安全预警方法及装置

Also Published As

Publication number Publication date
CN102360526B (zh) 2013-09-18

Similar Documents

Publication Publication Date Title
CN102360526B (zh) 一种快速道路路段状态的实时监测方法
CN102360525B (zh) 基于判别分析的快速道路实时交通事故风险预测方法
CN103646534B (zh) 一种道路实时交通事故风险控制方法
CN103021176B (zh) 基于断面检测器的城市道路交通状态判别方法
CN104157156B (zh) 一种高速公路危险路段车速动态管理预警方法
CN104240499B (zh) 一种基于微波数据的异常拥堵点判别方法
CN104021671B (zh) 一种svm与模糊判断相结合的道路实时路况的判断方法
CN103150900A (zh) 一种基于视频的交通拥堵事件自动检测方法
CN110120153A (zh) 一种公交驾驶事故风险评估***及其方法
CN104484996A (zh) 一种基于多源数据的路段交通状态判别方法
CN103198711B (zh) 一种降低不同严重程度交通事故概率的车辆调控方法
CN102968901A (zh) 获取区域拥堵信息的方法及区域拥堵分析装置
CN101783074A (zh) 一种城市道路交通流状态实时判别方法及***
CN101739825A (zh) 基于gps浮动车的交通数据故障识别与修复方法
CN108922168B (zh) 一种中观层面事故多发道路判别的方法
CN102881162A (zh) 大规模交通信息的数据处理及融合方法
CN101882374A (zh) 一种车辆行程分布综合信息分析***及方法
CN102735320A (zh) 基于桥梁动应变识别车辆重量的方法
CN103198713A (zh) 基于交通数据和天气数据的减少交通事故的车辆调控方法
CN102360524B (zh) 快速道路危险交通流特征的自动检测确认方法
CN101465059A (zh) 城市道路交通安全态势鉴判预警***
CN100481153C (zh) 基于偏最小二乘原理的高速公路交通事件自动检测方法
CN103198707B (zh) 一种基于晴天状况下危险交通流特征的车辆调控方法
CN101551940A (zh) 城市快速道路交通状态判别与发布***及其方法
CN105869398A (zh) 一种基于K-means聚类的道路交通畅通度判断方法

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