CN108847022A - 一种微波交通数据采集设备的异常值检测方法 - Google Patents
一种微波交通数据采集设备的异常值检测方法 Download PDFInfo
- Publication number
- CN108847022A CN108847022A CN201810585465.6A CN201810585465A CN108847022A CN 108847022 A CN108847022 A CN 108847022A CN 201810585465 A CN201810585465 A CN 201810585465A CN 108847022 A CN108847022 A CN 108847022A
- Authority
- CN
- China
- Prior art keywords
- data
- cluster
- collection equipment
- classification
- valid
- 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
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/04—Detecting movement of traffic to be counted or controlled using optical or ultrasonic detectors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
- G08G1/0133—Traffic data processing for classifying traffic situation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Data Mining & Analysis (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Probability & Statistics with Applications (AREA)
- Traffic Control Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种微波交通数据采集设备的异常值检测方法,包括:(1)读取微波交通数据采集设备采集的数据,对同一维度的数据进行有效维度筛选,得到多维度有效数据;(2)对多维度有效数据进行标准化处理;(3)对标准化后的多维有效数据进行聚类,得到相应的聚类结果;(4)对聚类结果中各个类别中的数据的表现进行判定,确定哪些类别的数据为异常数据。本发明的有益效果:可以有效利用微波数据的多维度特征,在全局角度筛选出数据异常值,摒弃了之前只从单维度分析异常值的办法。提升异常值检测效率,在海量数据层面上也具有可操作性,节省时间和技术人员的精力。
Description
技术领域
本发明属于城市交通领域,涉及一种基于Kmeans聚类算法的微波交通数据采集设备的异常值检测方法。
背景技术
在交通领域存在多种数据采集方式和设备,其中微波设备由于其稳定性而被广泛安装在城市道路上。对交通数据采集设备收集的数据进行分析和利用在整治交通拥堵,规划道路等城市发展问题上具有指导意义。然而由于客观原因的存在,交通采集设备所收集的数据往往具有异常值,从而导致数据质量下降,影响了后续对数据的分析,所以需要对原始数据进行异常值检测,保证数据质量。根据以往经验,交通采集设备所收集的数据会包含统计性异常值和业务性异常值。统计性异常值一般都是由于设备故障而产生的,而业务性异常值属于不符合业务规律的。
现有的异常值检测方案一般还是以技术人员的经验判断为准人为筛选出异常值。这种方法效率不高而且难以在多维度数据上使用。与此同时,人为判断异常值也常常会因为数据量的增大忽略掉一些潜在的异常值。跟现有的异常值检测技术相比,聚类算法可以应对多维度数据的情况,通过把数据点按照之间在多维空间上的欧式距离远近归类找出“少且异”的异常值。
发明内容
本发明的目的是提供一种提升异常值检测效率,在海量数据层面上也具有可操作性,节省时间和技术人员的精力的微波交通数据采集设备的异常值检测方法。
本发明采用的技术方案是:
一种微波交通数据采集设备的异常值检测方法,其特征在于:包括:
(1)读取微波交通数据采集设备采集的数据,对同一维度的数据进行有效维度筛选,得到多维度有效数据;
(2)对多维度有效数据进行标准化处理;
(3)对标准化后的多维有效数据进行聚类,得到相应的聚类结果;
(4)对聚类结果中各个类别中的数据的表现进行判定,确定哪些类别的数据为异常数据。本发明针对交通采集设备微波所采集的数据特点,运用聚类算法对其多维度特征进行分类,达到异常值检测的效果。同时,在判断聚类类别是否为异常值时,需要运用交通理论的知识。
进一步,所述多维度有效数据包括车流量、车辆过车速度、车辆占有率。原始微波数据包含很多维度,但是有些维度如纪录自增ID对异常值检测起不了显著作用,而有些维度虽然有字段存在但是并没有采集到任何数据,这些维度都可以视为无效维度,剔除这些之后剩下的即有意义的维度,即车流量,车辆过车速度,车辆占有率。
进一步,步骤(2)中的标准化处理采用min-max标准化处理。筛选出有效维度之后,因为各个维度的取值范围不一,在进行聚类的时候会对结果造成不良影响,所以需要对数据进行标准化处理。在对数据进行min-max标准化处理之后,各个维度对结果都具有同等的解释作用。
进一步,步骤(3)聚类采用Kmeans聚类,具体步骤包括:
1)设定聚类类别数目k,需要聚类的类别数目需要事先确定,最后的结果会把所有原始数据点分为k个类别,而每个类别中的数据都具有相似性;
2)随机分配k个聚类中心,根据给定的数目k自生成k个聚类中心;
3)对每一数据点都计算其与所有聚类中心的欧式距离,然后把数据点分配给最近的聚类中心;
4)重复2)和3)的内容,直到聚类中心收敛,即所有数据点不会随着聚类中心的变化改变其所依附的聚类中心。本发明最后得到的数据是按照以k个聚类中心所代表的k个类别组织的;每个类别内的成员都具有相似性而与其他类别不相似。
或者,步骤(3)中的聚类是根据已有的异常数据的样本值进行训练,对获得异常数据的特征聚类并建立异常数据模型,根据异常数据模型对在新数据中找出类似的异常数据。
进一步,步骤(4)中的判定方法包括业务员的经验和相关维度之间的相关性关系。
进一步,业务人员的经验指的是聚类出来的一个类别之中将近99%的数据表现为该类别中的数据包括多个维度是无效数据或者是仅有一个维度是有效数据,则该类别的数据是异常数据。比如:聚类出来的一个类别之中将近99%的数据表现为流量为0,占有率为0,速度为240,这些属于有意义的有效数据但是这显然属于异常值。
进一步,相关维度之间的相关性关系指的是多维度有效数据中的各维度之间存在一定的线性关系,若某个类别中的数据在所述线性关系的置信区间内,则该类别的数据可以视作异常数据。比如车流量=车辆过车速度*密度;车辆占有率=车辆平均长度*密度,根据交通理论中这两个关于车流量,车辆过车速度,车辆占有率的公式,可以看到通过密度这个维度可以把三者直接放在一个公式下:车流量=车辆过车速度*车辆占有率/车辆平均长度,即三者存在一定的线性关系。当得到聚类结果的时候,可以通过对各个类别中数据的表现进行判定,若这个类别中的数据在这个线性关系的置信区间内,则这个类别的数据可以视作异常数据。
进一步,对步骤(4)中确定的异常数据打上数据异常的标签。
本发明的有益效果:
1、可以有效利用微波数据的多维度特征,在全局角度筛选出数据异常值,摒弃了之前只从单维度分析异常值的办法。
2、提升异常值检测效率,在海量数据层面上也具有可操作性,节省时间和技术人员的精力。
附图说明
图1为本发明的流程示意图。
具体实施方式
下面结合具体实施例来对本发明进行进一步说明,但并不将本发明局限于这些具体实施方式。本领域技术人员应该认识到,本发明涵盖了权利要求书范围内所可能包括的所有备选方案、改进方案和等效方案。
参见图1,本实施例提供了一种微波交通数据采集设备的异常值检测方法,包括:
(1)读取微波交通数据采集设备采集的数据,对同一维度的数据进行有效维度筛选,得到多维度有效数据;
本实施例所述多维度有效数据包括车流量、车辆过车速度、车辆占有率。原始微波数据包含很多维度,但是有些维度如纪录自增ID对异常值检测起不了显著作用,而有些维度虽然有字段存在但是并没有采集到任何数据,这些维度都可以视为无效维度,剔除这些之后剩下的即有意义的维度,即车流量,车辆过车速度,车辆占有率。具体的,现在可以拿到的一条微波数据包括以下维度:自增ID,设备编码,设备车道编码,车道流量,C1车型流量,C2车型流量,C3车型流量,C4车型流量,C5车型流量,车道平均速度,车头时距,车头间距,车道占有率,车道饱和度,采集时间,记录时间。对于异常值来说,类似自增ID这种纯粹标记性的维度是没有作用的。其次,虽然存在类似车头时距,车头间距这样的字段,但是这个字段中没有数据,即数据都为0,所以将其视为无效。在这些维度当中,车流量,平均速度,和占有率是比较有效的,首先因为这三个变量是交通理论中重要的概念,其次这三个维度的数据质量比较高,基本不存在0的情况。所以假如我们要判断一条微波数据是否为异常值,可以从这三个重要且有效的维度入手判断。
(2)对多维度有效数据进行标准化处理;
筛选出有效维度之后,因为各个维度的取值范围不一,在进行聚类的时候会对结果造成不良影响,所以需要对数据进行标准化处理。在对数据进行min-max标准化处理之后,各个维度对结果都具有同等的解释作用。
(3)对标准化后的多维有效数据进行Kmeans聚类,得到相应的聚类结果;
Kmeans聚类具体步骤包括:
1)设定聚类类别数目k,需要聚类的类别数目需要事先确定,最后的结果会把所有原始数据点分为k个类别,而每个类别中的数据都具有相似性;
2)随机分配k个聚类中心,根据给定的数目k自生成k个聚类中心;
3)对每一数据点都计算其与所有聚类中心的欧式距离,然后把数据点分配给最近的聚类中心;
4)重复2)和3)的内容,直到聚类中心收敛,即所有数据点不会随着聚类中心的变化改变其所依附的聚类中心。本发明最后得到的数据是按照以k个聚类中心所代表的k个类别组织的;每个类别内的成员都具有相似性而与其他类别不相似。
(4)对聚类结果中各个类别中的数据的表现通过业务员的经验和相关维度之间的相关性关系进行判定,确定哪些类别的数据为异常数据并打上数据异常的标签。具体的业务人员的经验指的是聚类出来的一个类别之中将近99%的数据表现为该类别中的数据包括多个维度是无效数据或者是仅有一个维度是有效数据,则该类别的数据是异常数据。比如:聚类出来的一个类别之中将近99%的数据表现为流量为0,占有率为0,速度为240,这些属于有意义的有效数据但是这显然属于异常值。相关维度之间的相关性关系指的是多维度有效数据中的各维度之间存在一定的线性关系,若某个类别中的数据在所述线性关系的置信区间内,则该类别的数据可以视作异常数据。比如车流量=车辆过车速度*密度;车辆占有率=车辆平均长度*密度,根据交通理论中这两个关于车流量,车辆过车速度,车辆占有率的公式,可以看到通过密度这个维度可以把三者直接放在一个公式下:车流量=车辆过车速度*车辆占有率/车辆平均长度,即三者存在一定的线性关系。当得到聚类结果的时候,可以通过对各个类别中数据的表现进行判定,若这个类别中的数据在这个线性关系的置信区间内,则这个类别的数据可以视作异常数据。
当异常数据的样本足够多时,步骤(3)中的聚类是可以根据已有的异常数据的样本值进行训练,对获得异常数据的特征聚类并建立异常数据模型,根据异常数据模型对在新数据中找出类似的异常数据。
本发明针对交通采集设备微波所采集的数据特点,运用聚类算法对其多维度特征进行分类,达到异常值检测的效果。同时,在判断聚类类别是否为异常值时,需要运用交通理论的知识。将聚类算法运用在交通数据上,同时结合了交通理论中流量速度和占有率的关系判断异常值。先通过聚类找出疑似异常值,再根据变量之间的相关性关系判断异常值,这一套流程在交通领域目前是具有先进性的。
Claims (9)
1.一种微波交通数据采集设备的异常值检测方法,其特征在于:包括:
(1)读取微波交通数据采集设备采集的数据,对同一维度的数据进行有效维度筛选,得到多维度有效数据;
(2)对多维度有效数据进行标准化处理;
(3)对标准化后的多维有效数据进行聚类,得到相应的聚类结果;
(4)对聚类结果中各个类别中的数据的表现进行判定,确定哪些类别的数据为异常数据。
2.根据权利要求1所述的一种微波交通数据采集设备的异常值检测方法,其特征在于:所述多维度有效数据包括车流量、车辆过车速度、车辆占有率。
3.根据权利要求1所述的一种微波交通数据采集设备的异常值检测方法,其特征在于:步骤(2)中的标准化处理采用min-max标准化处理。
4.根据权利要求1~3之一所述的一种微波交通数据采集设备的异常值检测方法,其特征在于:步骤(3)聚类采用Kmeans聚类,具体步骤包括:
1)设定聚类类别数目k,需要聚类的类别数目需要事先确定,最后的结果会把所有原始数据点分为k个类别,而每个类别中的数据都具有相似性;
2)随机分配k个聚类中心,根据给定的数目k自生成k个聚类中心;
3)对每一数据点都计算其与所有聚类中心的欧式距离,然后把数据点分配给最近的聚类中心;
4)重复2)和3)的内容,直到聚类中心收敛,即所有数据点不会随着聚类中心的变化改变其所依附的聚类中心。
5.根据权利要求1~3之一所述的一种微波交通数据采集设备的异常值检测方法,其特征在于:步骤(3)中的聚类是根据已有的异常数据的样本值进行训练,对获得异常数据的特征聚类并建立异常数据模型,根据异常数据模型对在新数据中找出类似的异常数据。
6.根据权利要求4所述的一种微波交通数据采集设备的异常值检测方法,其特征在于:步骤(4)中的判定方法包括业务员的经验和相关维度之间的相关性关系。
7.根据权利要求6所述的一种微波交通数据采集设备的异常值检测方法,其特征在于:业务人员的经验指的是聚类出来的一个类别之中将近99%的数据表现为该类别中的数据包括多个维度是无效数据或者是仅有一个维度是有效数据,则该类别的数据是异常数据。
8.根据权利要求6所述的一种微波交通数据采集设备的异常值检测方法,其特征在于:相关维度之间的相关性关系指的是多维度有效数据中的各维度之间存在一定的线性关系,若某个类别中的数据在所述线性关系的置信区间内,则该类别的数据可以视作异常数据。
9.根据权利要求1所述的一种微波交通数据采集设备的异常值检测方法,其特征在于:对步骤(4)中确定的异常数据打上数据异常的标签。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810585465.6A CN108847022B (zh) | 2018-06-08 | 2018-06-08 | 一种微波交通数据采集设备的异常值检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810585465.6A CN108847022B (zh) | 2018-06-08 | 2018-06-08 | 一种微波交通数据采集设备的异常值检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108847022A true CN108847022A (zh) | 2018-11-20 |
CN108847022B CN108847022B (zh) | 2021-04-16 |
Family
ID=64210324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810585465.6A Active CN108847022B (zh) | 2018-06-08 | 2018-06-08 | 一种微波交通数据采集设备的异常值检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108847022B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110007660A (zh) * | 2019-04-10 | 2019-07-12 | 河北工业大学 | 一种火电厂汽轮机组瞬态等效热应力在线软测量方法 |
CN110021164A (zh) * | 2019-03-02 | 2019-07-16 | 合肥学院 | 基于行驶时间数据的网约车路网占有率分析方法 |
CN110021163A (zh) * | 2019-03-02 | 2019-07-16 | 合肥学院 | 基于行驶里程数据的网约车路网占有率分析方法 |
CN111160401A (zh) * | 2019-12-09 | 2020-05-15 | 国网辽宁省电力有限公司电力科学研究院 | 一种基于均值漂移和XGBoost的异常用电判别方法 |
CN112434260A (zh) * | 2020-10-21 | 2021-03-02 | 北京千方科技股份有限公司 | 一种公路的交通状态检测方法、装置、存储介质及终端 |
CN114495498A (zh) * | 2022-01-20 | 2022-05-13 | 青岛海信网络科技股份有限公司 | 一种交通数据分布有效性判别方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727749A (zh) * | 2009-11-19 | 2010-06-09 | 北京世纪高通科技有限公司 | 异常车速数据的识别方法和装置 |
CN103839409A (zh) * | 2014-02-27 | 2014-06-04 | 南京大学 | 基于多断面视觉传感聚类分析的交通流状态判别方法 |
CN104809877A (zh) * | 2015-05-14 | 2015-07-29 | 重庆大学 | 基于特征参数加权gefcm算法的高速公路地点交通状态估计方法 |
CN106919953A (zh) * | 2017-02-23 | 2017-07-04 | 北京工业大学 | 一种基于轨道交通数据分析的异常出行群体识别方法 |
WO2017154947A1 (ja) * | 2016-03-10 | 2017-09-14 | 新日本理化株式会社 | 粉体状1,4-シクロヘキサンジカルボン酸 |
-
2018
- 2018-06-08 CN CN201810585465.6A patent/CN108847022B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727749A (zh) * | 2009-11-19 | 2010-06-09 | 北京世纪高通科技有限公司 | 异常车速数据的识别方法和装置 |
CN103839409A (zh) * | 2014-02-27 | 2014-06-04 | 南京大学 | 基于多断面视觉传感聚类分析的交通流状态判别方法 |
CN104809877A (zh) * | 2015-05-14 | 2015-07-29 | 重庆大学 | 基于特征参数加权gefcm算法的高速公路地点交通状态估计方法 |
WO2017154947A1 (ja) * | 2016-03-10 | 2017-09-14 | 新日本理化株式会社 | 粉体状1,4-シクロヘキサンジカルボン酸 |
CN106919953A (zh) * | 2017-02-23 | 2017-07-04 | 北京工业大学 | 一种基于轨道交通数据分析的异常出行群体识别方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110021164A (zh) * | 2019-03-02 | 2019-07-16 | 合肥学院 | 基于行驶时间数据的网约车路网占有率分析方法 |
CN110021163A (zh) * | 2019-03-02 | 2019-07-16 | 合肥学院 | 基于行驶里程数据的网约车路网占有率分析方法 |
CN110021163B (zh) * | 2019-03-02 | 2020-10-13 | 合肥学院 | 基于行驶里程数据的网约车路网占有率分析方法 |
CN110007660A (zh) * | 2019-04-10 | 2019-07-12 | 河北工业大学 | 一种火电厂汽轮机组瞬态等效热应力在线软测量方法 |
CN110007660B (zh) * | 2019-04-10 | 2020-06-16 | 河北工业大学 | 一种火电厂汽轮机组瞬态等效热应力在线软测量方法 |
CN111160401A (zh) * | 2019-12-09 | 2020-05-15 | 国网辽宁省电力有限公司电力科学研究院 | 一种基于均值漂移和XGBoost的异常用电判别方法 |
CN112434260A (zh) * | 2020-10-21 | 2021-03-02 | 北京千方科技股份有限公司 | 一种公路的交通状态检测方法、装置、存储介质及终端 |
CN114495498A (zh) * | 2022-01-20 | 2022-05-13 | 青岛海信网络科技股份有限公司 | 一种交通数据分布有效性判别方法及装置 |
CN114495498B (zh) * | 2022-01-20 | 2023-01-10 | 青岛海信网络科技股份有限公司 | 一种交通数据分布有效性判别方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108847022B (zh) | 2021-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108847022A (zh) | 一种微波交通数据采集设备的异常值检测方法 | |
WO2020238631A1 (zh) | 一种基于手机信令数据的人群类型识别方法 | |
CN108922174B (zh) | 一种快速路入口匝道周围交叉口群内路径动态分级方法 | |
CN109753499A (zh) | 一种运维监控数据治理方法 | |
CN105719045A (zh) | 留任风险确定器 | |
CN109191828B (zh) | 基于集成学习的交通参与者事故风险预测方法 | |
CN111640304B (zh) | 面向连续流交通设施的交通拥堵传播特征自动化量化提取方法 | |
CN113689594A (zh) | 一种新的汽车行驶工况构造方法 | |
CN104156403A (zh) | 一种基于聚类的大数据常态模式提取方法及系统 | |
CN106021386B (zh) | 面向海量分布式数据的非等值连接方法 | |
CN116013084B (zh) | 交通管控场景确定方法、装置、电子设备及存储介质 | |
CN106650113A (zh) | 一种基于模糊聚类的桥梁监测数据异常情况识别方法 | |
CN105049286A (zh) | 基于层次聚类的云平台测速数据判定方法 | |
Xu et al. | A novel algorithm for urban traffic congestion detection based on GPS data compression | |
CN116307682A (zh) | 一种考虑多元指标的机场航线网络失稳状态综合评价方法 | |
CN116894744A (zh) | 一种基于改进k均值聚类算法的电网用户数据分析方法 | |
CN107944487A (zh) | 一种基于混合协同过滤算法的作物育种品种推荐方法 | |
CN104361259B (zh) | 一种变压器越限运行特征分析方法 | |
CN107070748A (zh) | 一种通信大数据的处理系统及方法 | |
CN110610446A (zh) | 一种基于两步聚类思路的县域城镇分类方法 | |
CN109344729B (zh) | 一种识别道路中人员运动的方法 | |
CN106651630A (zh) | 关键用电行业识别方法和系统 | |
CN112241820A (zh) | 资金流动中关键节点的风险识别方法、装置及计算设备 | |
CN115618229A (zh) | 一种基于大数据的线损分摊算法 | |
CN114328513A (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 |