CN104833361B - 复杂路况下的基于多权重值的地图匹配方法 - Google Patents

复杂路况下的基于多权重值的地图匹配方法 Download PDF

Info

Publication number
CN104833361B
CN104833361B CN201510222397.3A CN201510222397A CN104833361B CN 104833361 B CN104833361 B CN 104833361B CN 201510222397 A CN201510222397 A CN 201510222397A CN 104833361 B CN104833361 B CN 104833361B
Authority
CN
China
Prior art keywords
road section
vehicle
road
matched
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.)
Active
Application number
CN201510222397.3A
Other languages
English (en)
Other versions
CN104833361A (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.)
Nanjing Post and Telecommunication University
Original Assignee
Nanjing Post and Telecommunication 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 Nanjing Post and Telecommunication University filed Critical Nanjing Post and Telecommunication University
Priority to CN201510222397.3A priority Critical patent/CN104833361B/zh
Publication of CN104833361A publication Critical patent/CN104833361A/zh
Application granted granted Critical
Publication of CN104833361B publication Critical patent/CN104833361B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/26Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/28Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network with correlation of data from several navigational instruments
    • G01C21/30Map- or contour-matching

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Navigation (AREA)

Abstract

本发明公开了复杂路况下的基于多权重值的地图匹配方法,该方法解决了复杂路况下车辆卫星导航系统GNSS定位点地图匹配准确度不高的问题。该方法使用车辆前后GNSS定位点连线与路段夹角、定位点与路段投影距离以及航向角与路段夹角作为主要因素,利用GNSS浮动车辆的数据进行轨迹修正,在复杂路段处可以修正偏移的GNSS定位点并提供较高的地图匹配精确度。

Description

复杂路况下的基于多权重值的地图匹配方法
技术领域
本发明涉及车联网中车辆基于GNSS定位点的地图匹配技术,特别涉及一种复杂路况下的基于多权重值的地图匹配方法。
背景技术
当前,当车辆行驶在一般的道路情况下时,使用一般的地图匹配方法都可以对于车辆行驶轨迹进行准确修正,但是,当车辆行驶在复杂的道路情况下时,如U型弯道或者十字路口,由于车辆定位点的模糊以及普通地图匹配算法的局限性,很容易将车辆匹配到错误的道路上,因此就无法提供准确的车辆行驶轨迹。如图1所示,车辆连续采样时刻的定位点分别是P1→P2→P3→P4→P5→P6。图中有三条待匹配路段L1、L2和L3,如果使用最短距离投影算法进行匹配,P4很有可能会被匹配到路段L1上,这样的匹配结果显然是错误的。因此在面对复杂路况时,必须使用一种匹配准确度更好的地图匹配算法。
当车辆行驶在一般的道路情况下时,使用一般的地图匹配算法都可以对于车辆行驶轨迹进行准确修正。但是当车辆行驶在复杂的道路情况下时,如U型弯道或者十字路口,由于车辆定位点的模糊以及普通地图匹配算法的局限性,很容易将车辆匹配到错误的道路上,因此就无法提供准确的车辆行驶轨迹。而本发明能够很好地解决上面的问题。
发明内容
本发明目的在于提供了一种复杂路况下的基于多权重值的地图匹配方法,该方法解决了复杂路况下车辆卫星导航系统GNSS定位点地图匹配准确度不高的问题。该方法使用车辆前后GNSS定位点连线与路段夹角、定位点与路段投影距离以及航向角与路段夹角作为主要因素,利用GNSS浮动车辆的数据进行轨迹修正,在复杂路段处可以修正偏移的GNSS定位点并提供较高的地图匹配精确度。
本发明解决其技术问题所采取的技术方案是:一种复杂路况下的基于多权重值的地图匹配方法,该方法包括如下步骤:
步骤1:获取车辆前后GNSS定位点连线与路段夹角、定位点与路段投影距离以及航向角与路段夹角;车辆航向角(VCA,Vehicle Course Angle)可以通过GNSS定位数据得到,将VCA定义为车辆行驶方向与地图正北向沿顺时针所形成的夹角,范围是[0,360°)。路段方位角(SOA,Sections Of Azimuth)是指路段上相邻两个节点间沿着车辆行驶方向的路段与地图正北向沿顺时针所形成的的夹角,范围同样是[0,360°)。SOA可以通过电子地图数据库中的道路信息得到,也可以通过直行道路的起始点坐标求斜率得知。图2显示了VCA和SOA之间的关系。VCA与SOA之间的夹角记作也就是之间的差值;定位点垂直投影距离是指待匹配路段与车辆定位点的最近距离,因此首先必须将车辆定位点对于待匹配路段进行垂直投影,图3展示了垂直投影的过程。图中P点是当前车辆定位点,将P点对待匹配路段作垂直投影,得到垂直投影点P1,A点和B点分别是待匹配路段的起始端点。前后两个采样时刻的车辆定位点的连线与待匹配路段之间的夹角θ也是决定车辆行驶轨迹与路段相似度的重要因素。θ值越小,轨迹与路段相似度越高,反之则越低。图4清楚展示了前后定位点与路段之间的夹角。
步骤2:将获取的三个关键参数作为主要因素,利用GNSS浮动车辆的数据进行轨迹修正,在复杂路段处修正偏移的GNSS定位点并提供较高的地图匹配精确。基于多权重值的地图匹配算法设定了三个权重参量Wd和Wθ作为地图匹配的主要影响因素,分别通过权重因子λd和λθ对其进行约束,可以根据不同的路况选择不同的k值来调整三者的大小关系:当定位点与路段较近,即d较小时,可以将k值调小,此时W主要由Wd和Wθ决定;当定位点与路段较远,即d较大时,可以将k值调大,此时W主要由决定。
本发明使用当前采样时刻车辆的航向角与待匹配路段的夹角、当前采样时刻车辆定位点与待匹配路段的垂直投影距离以及当前时刻与前一时刻车辆定位点的连线与待匹配路段的夹角作为地图匹配的主要因素,并根据不同的路况选择不同的权重因子。
有益效果:
1、本发明能够根据不同的路况选择不同的权重因子,提高了地图匹配精确度。
2、当面对一般的道路情况时,使用一般的地图匹配算法都可以对于车辆行驶轨迹进行准确修正。
附图说明
图1为U型弯道车辆轨迹图。
图2为车辆航向角和路段方位角示意图。
图3为定位点垂直投影示意图。
图4为前后定位点连线与路段夹角示意图。
图5为基于多权重值的地图匹配算法在U型弯道的应用示意图。
图6为本发明的方法流程图。
具体实施方式
下面结合说明书附图对本发明创造作进一步的详细说明。
与本发明相关的缩略语和关键术语定义:
AGNSS Assisted Global Positioning System 位辅系助统全 球卫星定
GNSS Global Navigation Satellite System 全统 球卫星导航系
IOV Internet of Vehicle 车联网
SOA Sections Of Azimuth 路段方位角
V2V Vehicle to Vehicle 车与车之间的通信
VCA Vehicle Course Angle 车辆航向角
当面对一般的道路情况时,使用一般的地图匹配算法都可以对于车辆行驶轨迹进行准确修正,但是,当车辆行驶在复杂的道路情况下时,如U型弯道或者十字路口,由于车辆定位点的模糊以及普通地图匹配算法的局限性,很容易将车辆匹配到错误的道路上,因此就无法提供准确的车辆行驶轨迹。如图1所示,车辆连续采样时刻的定位点分别是P1→P2→P3→P4→P5→P6。图中有三条待匹配路段L1、L2和L3,如果使用最短距离投影算法进行匹配,P4很有可能会被匹配到路段L1上,这样的匹配结果显然是错误的。因此在面对复杂路况时,必须使用一种匹配准确度更好的地图匹配算法。
本发明提出的基于多权重值的地图匹配算法主要与三种因素相关,分别是当前采样时刻车辆的航向角与待匹配路段的夹角、当前采样时刻车辆定位点与待匹配路段的垂直投影距离以及当前时刻与前一时刻车辆定位点的连线与待匹配路段的夹角。接下来我们具体分析这三个组成要素,由于此算法需要解决复杂路况下GNSS定位点匹配准确度不高的问题,所以这里以使用GNSS进行定位作为前提。
(1)车辆航向角与路段夹角
由于这里是使用GNSS定位,所以车辆航向角(VCA,Vehicle Course Angle)可以通过GNSS定位数据得到,这里将VCA定义为车辆行驶方向与地图正北向沿顺时针所形成的夹角,范围是[0,360°)。同理,路段方位角(SOA,Sections Of Azimuth)是指路段上相邻两个节点间沿着车辆行驶方向的路段与地图正北向沿顺时针所形成的的夹角,范围同样是[0,360°)。SOA可以通过电子地图数据库中的道路信息得到,也可以通过直行道路的起始点坐标求斜率得知。图2显示了VCA和SOA之间的关系。
假设VCA与SOA之间的夹角记作也就是之间的差值。的变化范围为[0,180°],可以由公式(1)表示:
(2)定位点垂直投影距离
定位点垂直投影距离是指待匹配路段与车辆定位点的最近距离,因此首先必须将车辆定位点对于待匹配路段进行垂直投影,图3展示了垂直投影的过程。图中P点是当前车辆定位点,将P点对待匹配路段作垂直投影,得到垂直投影点P1,A点和B点分别是待匹配路段的起始端点。由公式(2)可知P点到待匹配路段的垂直投影距离d。
(3)前后定位点连线与路段夹角
前后两个采样时刻的车辆定位点的连线与待匹配路段之间的夹角θ也是决定车辆行驶轨迹与路段相似度的重要因素。θ值越小,轨迹与路段相似度越高,反之则越低。图4清楚展示了前后定位点与路段之间的夹角。图中,假设前后两个连续采样时刻的车辆定位点分别为P1(x1,y1)和P2(x2,y2),路段L的方位角是θroad,那么θ可以通过公式(3)计算得到:
本发明提出的基于多权重值的地图匹配算法对每条待匹配路段进行评判,通过每条路段的多权重值的大小选择出匹配度最高的路段[54][55][56]。多权重值W共由3个部分组成,由公式(4)表示:
其中,表示VCA和SOA之间的相似度,计算公式如下:
就是VCA和SOA的夹角,表示权重因子,会随着系数k而变化,d由公式(4.6)得到,0<k<1)。越小,VCA和SOA之间的相似度就越高,也就越大,当大于90°时,可以为负值。
Wd表示GNSS定位点与待匹配路段之间的接近程度,由公式(6)表示:
上述公式中的d表示车辆定位点到待匹配路段的投影距离,可以由公式(2)得到,λd是可选权重因子,在这里取作15。d越小,表示车辆定位点与待匹配路段越相近,因此Wd也就越大,反之亦然。
Wθ可以通过对多组相邻定位点连线与路段夹角求均值实现,这里选取路段周围满足Wd>1的n个定位点,分别设为(x1,y1),(x2,y2),…,(xn,yn),然后由公式(3)得到相邻两点的θ值,分别为θ1,θ2,…,θn。Wθ可以由公式(7)计算得到:
上述公式中,λθ为权重因子,λθ=(1-k)λd,θi是集合{θ12,...,θn-1}中的元素,当λθ不变,θi越小,Wθ就越大。
本发明提出的基于多权重值的地图匹配算法设定了三个权重参量Wd和Wθ作为地图匹配的主要影响因素,分别通过权重因子λd和λθ对其进行约束,可以根据不同的路况选择不同的k值来调整三者的大小关系:当定位点与路段较近,即d较小时,可以将k值调小,此时W主要由Wd和Wθ决定;当定位点与路段较远,即d较大时,可以将k值调大,此时W主要由决定。这里再回顾图1中的U型弯道车辆轨迹图,我们加入VCA、SOA以及定位点到路段的距离等因素,如图5所示。图中每个采样时刻的VCA分别用α1、α2、…、α6表示,L1、L2和L3的SOA分别用β1、β2和β3表示,P4到三个待匹配路段L1、L2和L3的距离分别为8米,11米和13米。因为目前车辆行驶在U型弯道处,所以VCA应该为主要影响因素,所以将k值设为0.6。我们通过基于多权重值的地图匹配算法对P4点进行权重值计算,得到如下表格:
表1:基于多权重值的地图匹配算法对于P4点的权重值计算
由表1可以看出,当车辆行驶在U行弯道处时,由于车辆定位点距离路段较远,所以在W中占比重较大,Wd占比重较小。经过权重值总和的计算,可以得知P4点与路段L3的匹配度最高,因此最后将P4成功匹配至L3

Claims (1)

1.一种复杂路况下的基于多权重值的地图匹配方法,其特征在于,所述方法使用当前采样时刻车辆的航向角与待匹配路段的夹角、当前采样时刻车辆定位点与待匹配路段的垂直投影距离以及当前时刻与前一时刻车辆定位点的连线与待匹配路段的夹角作为地图匹配的因素,并根据不同的路况选择不同的权重因子,包括如下步骤:
步骤1:获取车辆前后GNSS定位点连线与路段夹角、定位点与路段投影距离以及航向角与路段夹角;车辆航向角通过GNSS定位数据得到,将VCA定义为车辆行驶方向与地图正北向沿顺时针所形成的夹角,范围是[0,360°),路段方位角是指路段上相邻两个节点间沿着车辆行驶方向的路段与地图正北向沿顺时针所形成的夹角,范围同样是[0,360°),路段方位角SOA通过电子地图数据库中的道路信息得到或者通过直行道路的起始点坐标求斜率得知;VCA与路段方位角SOA之间的夹角记作也就是之间的差值;定位点垂直投影距离是指待匹配路段与车辆定位点的最近距离,首先将车辆定位点对于待匹配路段进行垂直投影;垂直投影的过程的P点是当前车辆定位点,将P点对待匹配路段作垂直投影,得到垂直投影点P1,A点和B点分别是待匹配路段的起始端点,前后两个采样时刻的车辆定位点的连线与待匹配路段之间的夹角θ也是决定车辆行驶轨迹与路段相似度的重要因素,θ值越小,轨迹与路段相似度越高,反之则越低;
步骤2:将获取的三个关键参数作为因素,利用GNSS浮动车辆的数据进行轨迹修正,在复杂路段处修正偏移的GNSS定位点并提供较高的地图匹配精确;基于多权重值的地图匹配算法设定了三个权重参量Wd和Wθ作为地图匹配的影响因素,分别通过权重因子λd和λθ对其进行约束,根据不同的路况选择不同的系数k值来调整三者的大小关系:当定位点与路段较近,即垂直投影距离d较小时,将k值调小,此时W由Wd和Wθ决定;当定位点与路段较远,即d较大时,将k值调大,此时W由决定;
所述的基于多权重值的地图匹配算法对每条待匹配路段进行评判,通过每条路段的多权重值的大小选择出匹配度最高的路段,多权重值W共由3个部分组成,由公式(1)表示:
其中,表示VCA和路段方位角SOA之间的相似度,计算公式如下:
就是VCA和路段方位角SOA的夹角,表示权重因子,会随着系数k而变化, 越小,VCA和路段方位角SOA之间的相似度就越高,也就越大,当大于90°时,为负值,
Wd表示GNSS定位点与待匹配路段之间的接近程度,由公式(3 )表示:
上述公式中的d表示车辆定位点到待匹配路段的投影距离,由得到,λd是可选权重因子,在这里取作15,d越小,表示车辆定位点与待匹配路段越相近,因此Wd也就越大,反之亦然;
GNSS定位点与待匹配路段之间的接近程度Wθ能通过对多组相邻定位点连线与路段夹角求均值实现,这里选取路段周围满足Wd>1的n个定位点,分别设为(x1,y1),(x2,y2),…,(xn,yn),然后由得到n-1个相邻两点的θ值,分别为θ1,θ2,…,θn,θroad表示路段L的方位角,Wθ由公式(4 )计算得到:
上述公式中,λθ为权重因子,λθ=(1-k)λd,θi是集合{θ12,...,θn-1}中的元素,当λθ不变,θi越小,Wθ就越大。
CN201510222397.3A 2015-05-04 2015-05-04 复杂路况下的基于多权重值的地图匹配方法 Active CN104833361B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510222397.3A CN104833361B (zh) 2015-05-04 2015-05-04 复杂路况下的基于多权重值的地图匹配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510222397.3A CN104833361B (zh) 2015-05-04 2015-05-04 复杂路况下的基于多权重值的地图匹配方法

Publications (2)

Publication Number Publication Date
CN104833361A CN104833361A (zh) 2015-08-12
CN104833361B true CN104833361B (zh) 2019-02-19

Family

ID=53811398

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510222397.3A Active CN104833361B (zh) 2015-05-04 2015-05-04 复杂路况下的基于多权重值的地图匹配方法

Country Status (1)

Country Link
CN (1) CN104833361B (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10145691B2 (en) 2016-05-18 2018-12-04 Here Global B.V. Ambiguity map match rating
CN106447090A (zh) * 2016-09-07 2017-02-22 平安科技(深圳)有限公司 旅游规划路线的异常提醒方法及服务器
CN107941229B (zh) * 2016-10-13 2021-08-03 斑马智行网络(香港)有限公司 车辆定位方法和车辆定位系统
CN108022432B (zh) * 2016-10-31 2020-09-04 杭州海康威视系统技术有限公司 车辆在运行线路上行驶方向的确定方法以及装置
CN108106619A (zh) * 2016-11-25 2018-06-01 厦门雅迅网络股份有限公司 主辅路识别方法及其系统
DE102017208867A1 (de) * 2017-05-24 2018-11-29 Bayerische Motoren Werke Aktiengesellschaft Darstellung einer Fahrtroute durch ein Navigationssystem
CN109579839B (zh) * 2017-09-29 2020-11-03 阿里巴巴(中国)有限公司 一种平行路识别方法、平行路相似度确定方法及装置
CN111380546A (zh) * 2018-12-28 2020-07-07 沈阳美行科技有限公司 基于平行道路的车辆定位方法、装置、电子设备和介质
CN111380540B (zh) * 2018-12-29 2024-04-05 阿里巴巴集团控股有限公司 地图匹配方法及装置、介质、终端
CN109813327B (zh) * 2019-02-01 2022-06-07 安徽中科美络信息技术有限公司 一种车辆行驶轨迹缺失补偿方法
CN110006442B (zh) * 2019-04-17 2021-06-01 北京百度网讯科技有限公司 导航方法、装置、设备和介质
CN111006680B (zh) * 2019-12-04 2020-12-08 无锡物联网创新中心有限公司 一种基于v2i技术的自动驾驶车辆路径规划系统及方法
CN111896024A (zh) * 2020-07-24 2020-11-06 北京汽车股份有限公司 导航显示的控制方法、装置及ar-hud显示系统
CN112991806B (zh) * 2021-02-18 2022-01-07 安徽中科美络信息技术有限公司 一种车辆轨迹监测方法及装置
CN113267196B (zh) * 2021-05-19 2022-04-08 中移智行网络科技有限公司 车辆轨迹修正方法、终端及计算机可读存储介质
CN114485686B (zh) * 2022-04-07 2022-07-01 北京盈通恒信电力科技有限公司 用于车辆导航定位的方法、装置、电子设备、存储介质
CN115985097B (zh) * 2022-12-29 2024-04-26 浪潮通信信息系统有限公司 一种高速用户运行轨迹有效性判断方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175254A (zh) * 2011-01-13 2011-09-07 北京超图软件股份有限公司 导航定位纠偏方法和装置,定位导航系统
CN102879003A (zh) * 2012-09-07 2013-01-16 重庆大学 基于gps终端的面向车辆位置跟踪的地图匹配方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183256B (zh) * 2011-02-28 2012-12-26 重庆大学 行进车队的地图匹配方法
CN102589557B (zh) * 2012-01-13 2014-12-03 吉林大学 一种基于驾驶员行为特性和分对数模型的交叉口处地图匹配方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175254A (zh) * 2011-01-13 2011-09-07 北京超图软件股份有限公司 导航定位纠偏方法和装置,定位导航系统
CN102879003A (zh) * 2012-09-07 2013-01-16 重庆大学 基于gps终端的面向车辆位置跟踪的地图匹配方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于WSN节点定位与道路方位信息的地图匹配改进算法;杨绍恭等;《云南民族大学学报》;20110131;第20卷(第1期);50-52,70
顾及多因子影响的自适应地图匹配算法研究;蔡靖宇等;《现代测绘:江苏省测绘地理信息学会2014年学术年会会议论文集》;20141219;39-41

Also Published As

Publication number Publication date
CN104833361A (zh) 2015-08-12

Similar Documents

Publication Publication Date Title
CN104833361B (zh) 复杂路况下的基于多权重值的地图匹配方法
CN101922939B (zh) 一种导航过程中的地图匹配方法和装置
US10074270B2 (en) Clustering observations of objects along roads for navigation-related operations
CN111380539B (zh) 车辆定位、导航方法和装置及相关系统
US20170097642A1 (en) Apparatus for correcting vehicle location
CN106104656B (zh) 地图信息生成系统、方法及程序
Tao et al. Lane marking aided vehicle localization
CN104197945B (zh) 一种基于低采样率浮动车数据的全局投票地图匹配方法
CN108827292A (zh) 基于gnss和地面基站的组合式导航精确测速定位方法及系统
CN104062671B (zh) 曲率约束的gnss浮动车地图匹配方法及装置
WO2012118207A1 (ja) 局所地図生成装置、局所地図生成システム、グローバル地図生成装置、グローバル地図生成システム、及びプログラム
CN107328410A (zh) 用于定位自动驾驶车辆的方法和汽车电脑
CN106469505B (zh) 一种浮动车轨迹纠偏方法和装置
CN110308470B (zh) 车辆定位方法及车辆定位系统
CN108981702A (zh) 一种多位置联合粒子滤波的车辆定位方法
EP3093620B1 (en) System and method for detecting roundabouts from probe data using vector fields
CN109443370A (zh) 一种检测轨迹偏离的方法
CN107860399A (zh) 一种基于地图匹配的车载捷联惯导行进间精对准方法
CN109813332B (zh) 添加虚拟引导线的方法和装置
CN104615897A (zh) 基于低频gps数据的路段行程时间估算方法
Blazquez et al. Simple map-matching algorithm applied to intelligent winter maintenance vehicle data
CN107328423A (zh) 基于地图数据的弯道识别方法及其系统
CN113450455A (zh) 用于生成停车场的道路链路的地图的方法、设备和计算机程序产品
US8401787B2 (en) Method of determining drive lane using steering wheel model
JP2013032953A (ja) 位置判定装置およびナビゲーション装置並びに位置判定方法,プログラム

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant