CN110133695B - 一种双天线gnss位置延迟时间动态估计系统及方法 - Google Patents
一种双天线gnss位置延迟时间动态估计系统及方法 Download PDFInfo
- Publication number
- CN110133695B CN110133695B CN201910314091.9A CN201910314091A CN110133695B CN 110133695 B CN110133695 B CN 110133695B CN 201910314091 A CN201910314091 A CN 201910314091A CN 110133695 B CN110133695 B CN 110133695B
- Authority
- CN
- China
- Prior art keywords
- inertial navigation
- antenna
- gnss
- error
- module
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 230000003044 adaptive effect Effects 0.000 claims abstract description 12
- 238000005259 measurement Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 description 12
- 230000001133 acceleration Effects 0.000 description 4
- 230000004927 fusion Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/52—Determining velocity
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种双天线GNSS位置延迟时间动态估计系统及方法,本发明系统包括双天线RTK GNSS接收机模块,用于输出GNSS主天线处位置、GNSSS双天线、GNSSS双天线经度信息和纬度信息;惯导模块,用于获取车辆角速度并根据车辆电机转速和陀螺仪角速度对惯导处位置和航向进行推算;初始化模块,用于根据双天线RTK GNSS接收机模块输出的GNSSS双天线航向、GNSSS双天线经度信息、纬度信息以及惯导模块获取的角速度信息对惯导模块进行初始化;自适应卡尔曼滤波惯导误差估计模块,用于对惯导模块的惯导处位置误差、航向误差以及位置延迟时间误差进行估计。与现有技术相比,本发明实现了对GNSS位置延迟时间的准确估计,且对实现GNSS与惯导时间同步具有重要的工程意义。
Description
技术领域
本发明涉及一种车辆GNSS位置延迟时间估计方法,尤其是涉及一种双天线GNSS位置延迟时间动态估计系统及方法。
背景技术
自动驾驶汽车近年来飞速发展,自动驾驶技术已成为整车企业甚至是IT企业竞争的焦点。准确的车辆位置和航向对于自动的驾驶实现至关重要。融合多源信息对位置、航向和位置延迟时间进行估计成为了亟需解决的问题。
目前国内外在本领域的研究点大多集中在位置和航向误差估计,现有的位置误差和航向估计的主要方法有:1、卫星导航定位系统(GNSS)通过RTK载波相位差分实现精确位置和航向计算,但其鲁棒性较差,GNSS信号极易受干扰;2、惯性导航系统(INS)实现对位置和航向计算,但其存在累积误差问题;3、基于轮速和INS组合的位置和航向估计,由于没有全局位置修正,存在累积误差问题;4、基于激光雷达和视觉等传感器的组合位置和航向估计,受环境特征点分布影响严重,并且各个传感器性能要求较高;5、基于多源信息融合的位置、航向和GNSS位置延迟时间估计,该方法可通过估计GNSS位置延迟时间实现对位置和航向的准确计算,但目前尚未有一种有效的能够同时对位置和航向误差进行估计的GNSS位置延迟时间估计方法。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种双天线GNSS位置延迟时间动态估计系统及方法。
本发明的目的可以通过以下技术方案来实现:
一种双天线GNSS位置延迟时间动态估计系统,该系统包括初始化模块、惯性导航模块、双天线RTK GNSS接收机模块和自适应卡尔曼滤波误差估计模块。
双天线RTK GNSS接收机模块用于输出GNSS主天线处位置、GNSSS双天线、GNSSS双天线经度信息和纬度信息;
惯导模块用于获取车辆角速度并根据车辆电机转速和陀螺仪角速度对惯导处位置和航向进行推算;
初始化模块用于根据双天线RTK GNSS接收机模块输出的GNSSS双天线航向、GNSSS双天线经度信息、纬度信息以及惯导模块获取的角速度信息对惯导模块进行初始化;
自适应卡尔曼滤波惯导误差估计模块用于对惯导模块的惯导处位置误差、航向误差以及位置延迟时间误差进行估计。
一种双天线GNSS位置延迟时间动态估计方法,该方法包括下列步骤:
步骤一、在车辆静止条件下,在一段时间内采用双天线RTK GNSS接收机模块输出的GNSS主天线位置和GNSSS双天线航向获取车辆位置和航向的初始值。
步骤二、惯导模块根据获取的角速度对车辆航向角进行估计。
步骤三、根据车辆电机转速和陀螺仪角速度,利用自适应卡尔曼滤波模块对惯导处位置误差和航向误差进行推算,并获取GNSS位置延迟时间误差。
惯导模块根据车辆电机转速和陀螺仪角速度对惯导处位置和航向进行推算的具体步骤包括:
1)估计得到车辆航向角后,使用该航向角将电机转速换算得到的合速度分解到导航坐标系下,计算经纬度;
2)转向时,令车辆的转向中心处于后轴之上,获取后轴中心处的东、北方向速度;
3)将后轴中心处的东、北方向速度转换到GNSS天线处,使后轴中心与GNSS天线的位置作位置融合。
其中,惯导处航向推算的表达式为:
惯导处位置推算的表达式为:
其中,LWDec和λWDec分别为惯导推算得到的纬度、经度,LDec和λDec分别为真实的纬度、经度,δLDec和δλDec分别为惯导的纬度和经度测量误差;
GNSS位置的测量方程为:
取惯导处位置推算结果与GNSS位置的测量方程之差作为观测量,获取纬度、经度、航向、角速度零偏以及位置延迟时间误差的观测方程为:
式中,L为纬度,下标Dec表示纬度小数部分,RM为地球卯酉圈半径,RN为地球子午圈半径,h表示车辆距水平面高度,vE、vN分别为后轴中心处的东、北方向速度;δT即为GNSS位置延迟时间的估计结果。
步骤四、将步骤三获取的惯导处位置误差和航向误差反馈至惯导模块,惯导模块获取惯导航向,重新执行步骤二。
与现有技术相比,本发明具有以下优点:
(1)本发明方法在低速下引入车辆电机转速通过惯导航向分解得到东北方向速度,进而通过积分得到位置,可以削弱惯性元件的误差对于速度估计的影响,提高位置估计精度;
(2)本发明方法充分融合车辆信息、惯导信息以及GNSS信息,并通过自适应卡尔曼滤波方法使得算法适应与GNSS异常工况,使用轮速辅助对GNSS的位置延迟时间进行估计,实现了对GNSS位置延迟时间的准确估计,且对实现GNSS与惯导时间同步具有重要的工程意义。
附图说明
图1为本发明方法的流程示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。显然,所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范围。
如图1所示,本发明涉及一种双天线GNSS位置延迟时间动态估计系统,包括初始化模块、惯性导航模块、双天线RTK GNSS接收机模块和自适应卡尔曼滤波误差估计模块。
所述双天线RTK GNSS接收机模块用于输出GNSS主天线处位置和双天线航向;所述初始化模块根据双天线RTK GNSS接收机模块输出的航向、经度以及纬度信息以及惯导输出角速度信息对惯导进行初始化;所述惯导模块根据车辆电机转速和陀螺仪角速度对惯导处位置和航向进行推算;所述自适应卡尔曼滤波惯导误差估计模块对位置延迟时间误差、惯导模块的位置、航向以及进行估计。
一种双天线GNSS位置延迟时间动态估计方法,该方法基于上述系统,包括下列步骤:
步骤一、在车辆静止条件下,在一段时间内(约20~30秒),采用双天线RTK GNSS接收机模块输出的GNSS天线位置和航向计算得到位置和航向的初始值,计算方法如下:
步骤二、惯导模块根据惯导测得的角速度对航向进行估计:
航向角由四部分构成,即:
第一部分是上一时刻航向值第二部分是横摆角增量由惯导传感器直接获取;第三部分为航向角估计误差 为反馈系数,由误差估计模块中航向角误差估计结果提供;第四部分是角速度零偏由静止时初始化模块输出的角速度零偏获取,Δt为惯导采样时间间隔。
惯导模块根据电机转速以及估计所得航向对位置进行推算,推算方法如下:
当估计得到航向角后,使用该航向可将电机转速换算得到的合速度分解到导航坐标系下来计算经纬度。转向时,车辆的转向中心处于后轴之上,因此对于后轴中心而言,在忽略侧偏的情况下,后轴中心点的航向即为速度方向。令后轴轮速换算的车速为VW,则后轴中心处的东、北方向速度分别为:
而GNSS天线测得的位置和速度是GNSS天线安放处的位置和速度,为了和GNSS天线的位置作位置融合,需要将后轴中心处的速度转换到GNSS天线处。
假设在车辆前左上坐标系下,天线相对于后轴中心处的三个方向距离分别为lx,ly,lz,则由于该距离产生的速度误差为:
式中,vU为天向速度。
则经过该补偿后,后轴中心处的东、北方向速度公式变为:
基于该速度,位置经纬度的积分公式为:
该积分公式由三部分构成:第一部分为上一时刻纬度Lk-1和经度λk-1;第二部分为东向和北向速度产生的纬度增量vN·dt/(RM+h)和经度增量vE·dt/(RN+h)/cosLk-1;第三部分为惯导纬度估计误差和惯导经度估计误差kL为纬度误差反馈系数,kλ为经度误差反馈系数。
步骤三、利用自适应卡尔曼滤波模块对位置误差、航向误差和GNSS位置延迟时间误差进行估计。估计方法如下:
为了对GNSS接收机输出位置的延迟时间进行估计,使用车辆轮速进行辅助,且认为轮胎半径已知,轮速换算得到的车速准确,这里在小纵向动态和小侧向动态工况下,使用估计所得航向角替代轨迹角,用来将轮速分解至导航坐标系下,并认为估计所得的航向角存在一个误差角,则即惯导推算航向的表达式为:
假设真实的东、北方向速度分别为VE和VN:
则得到东向和北向速度误差方程分别为δVEW和δVNW:
其中,δT为GNSS的位置延迟时间误差,为角速度的零偏误差,δLDec和δλDec分别为惯导纬度误差、惯导经度误差,和wδT分别为纬度误差的过程噪声、经度误差的过程噪声、航向角误差的过程噪声、角速度零偏误差的过程噪声和位置延迟时间误差的过程噪声。
由于在位置误差动态方程中,已将经纬度的整数部分和小数部分拆开并做过放大1e7处理,测量方程中也要做相应处理。
由轮速推算所得位置方程为:
其中,LWDec和λWDec分别为惯导推算得到的纬度、经度,LDec和λDec分别为真实的纬度、经度,δLDec和δλDec分别为惯导的纬度和经度测量误差。
而GNSS位置测量方程为:
取二者之差作为观测量,得到纬度、经度、航向、角速度零偏以及位置延迟时间误差的观测方程为:
式中,L为纬度,下标Dec表示纬度小数部分,RM为地球卯酉圈半径,RN为地球子午圈半径,h表示车辆距水平面高度。
基于以上的状态方程和观测方程便可以通过自适应卡尔曼滤波算法对这些误差进行估计。
步骤四、在误差估计模块中所使用的是自适应卡尔曼滤波算法,其中标准卡尔曼滤波算法如下:
P0|0=Var(x0)
Pk|k=(I-GkCk|k-1)Pk|k-1
其中P0|0为初始时刻状态误差的协方差矩阵,x0为初始状态,Var(x0)为初始状态的方差,E(x0)为初始状态的期望,为k时刻的状态预测值,Ak-1为k-1时刻的系统矩阵,为k-1时刻的估计值,Pk|k-1为k时刻协方差矩阵预测值,Pk-1|k-1为k-1时刻的协方差矩阵,Γk-1为过程噪声的输入矩阵,Qk-1为过程噪声的协方差矩阵,Ck|k-1为测量矩阵,Gk为卡尔曼滤波增益矩阵,Rk为测量噪声的协方差矩阵,Pk|k为k时刻协方差矩阵,I为单位阵,为k时刻状态的最优估计值,zk为k时刻的测量。
基于标准卡尔曼滤波算法,通过对其残差进行处理从而实现测量噪声的自适应提升算法的动态性能,自适应方法如下:
从而实现了误差估计模块卡尔曼滤波算法的自适应。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的工作人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (3)
1.一种应用双天线GNSS位置延迟时间动态估计系统的估计方法,其特征在于,该系统包括:
双天线RTK GNSS接收机模块,用于输出GNSS主天线处位置、GNSSS双天线、GNSSS双天线经度信息和纬度信息;
惯导模块,用于获取车辆角速度并根据车辆电机转速和陀螺仪角速度对惯导处位置和航向进行推算;
初始化模块,用于根据双天线RTK GNSS接收机模块输出的GNSSS双天线航向、GNSSS双天线经度信息、纬度信息以及惯导模块获取的角速度信息对惯导模块进行初始化;
自适应卡尔曼滤波惯导误差估计模块,用于对惯导模块的惯导处位置误差、航向误差以及位置延迟时间误差进行估计;
所述估计方法包括下列步骤:
1)在车辆静止条件下,在一段时间内采用双天线RTK GNSS接收机模块输出的GNSS主天线位置和GNSSS双天线航向获取车辆位置和航向的初始值;
2)惯导模块根据获取的角速度对车辆航向角进行估计;
3)根据车辆电机转速和陀螺仪角速度,利用自适应卡尔曼滤波模块对惯导处位置误差和航向误差进行推算,并获取GNSS位置延迟时间误差;
4)将步骤3)获取的惯导处位置误差和航向误差反馈至惯导模块,惯导模块获取惯导航向,重新执行步骤2);
步骤3)中,惯导模块根据车辆电机转速和陀螺仪角速度对惯导处位置和航向进行推算的具体内容为:
1)估计得到车辆航向角后,使用该航向角将电机转速换算得到的合速度分解到导航坐标系下,计算经纬度;
2)转向时,令车辆的转向中心处于后轴之上,获取后轴中心处的东、北方向速度;
3)将后轴中心处的东、北方向速度转换到GNSS天线处,使后轴中心与GNSS天线的位置作位置融合。
3.根据权利要求1所述的一种双天线GNSS位置延迟时间动态估计方法,其特征在于,步骤3)中,惯导处位置推算的表达式为:
其中,LWDec和λWDec分别为惯导推算得到的纬度、经度,LDec和λDec分别为真实的纬度、经度,δLDec和δλDec分别为惯导的纬度和经度测量误差;
GNSS位置的测量方程为:
取惯导处位置推算结果与GNSS位置的测量方程之差作为观测量,获取纬度、经度、航向、角速度零偏以及位置延迟时间误差的观测方程为:
式中,L为纬度,下标Dec表示纬度小数部分,RM为地球卯酉圈半径,RN为地球子午圈半径,h表示车辆距水平面高度,vE、vN分别为后轴中心处的东、北方向速度;δT即为GNSS的位置延迟时间误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910314091.9A CN110133695B (zh) | 2019-04-18 | 2019-04-18 | 一种双天线gnss位置延迟时间动态估计系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910314091.9A CN110133695B (zh) | 2019-04-18 | 2019-04-18 | 一种双天线gnss位置延迟时间动态估计系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110133695A CN110133695A (zh) | 2019-08-16 |
CN110133695B true CN110133695B (zh) | 2023-04-28 |
Family
ID=67570417
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910314091.9A Active CN110133695B (zh) | 2019-04-18 | 2019-04-18 | 一种双天线gnss位置延迟时间动态估计系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110133695B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111025908B (zh) * | 2019-12-23 | 2023-03-31 | 上海理工大学 | 基于自适应机动加速度扩展卡尔曼滤波器的航姿参考系统 |
CN111197994B (zh) * | 2019-12-31 | 2021-12-07 | 深圳一清创新科技有限公司 | 位置数据修正方法、装置、计算机设备和存储介质 |
CN114394130B (zh) * | 2021-12-27 | 2022-11-11 | 中国矿业大学 | 一种煤矿辅助运输车辆定位方法及定位系统 |
CN114353835B (zh) * | 2022-01-21 | 2022-10-11 | 中国铁道科学研究院集团有限公司铁道建筑研究所 | 惯性轨道测量仪动态校准系统、方法及其应用 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8600660B2 (en) * | 2006-09-29 | 2013-12-03 | Honeywell International Inc. | Multipath modeling for deep integration |
CN101750066B (zh) * | 2009-12-31 | 2011-09-28 | 中国人民解放军国防科学技术大学 | 基于卫星定位的sins动基座传递对准方法 |
CN103245963A (zh) * | 2013-05-09 | 2013-08-14 | 清华大学 | 双天线gnss/ins深组合导航方法及装置 |
CN103743395B (zh) * | 2014-01-17 | 2016-09-14 | 哈尔滨工程大学 | 一种惯性重力匹配组合导航系统中时间延迟的补偿方法 |
RU2617565C1 (ru) * | 2015-12-02 | 2017-04-25 | Акционерное общество "Раменское приборостроительное конструкторское бюро" | Способ оценивания ошибок инерциальной информации и её коррекции по измерениям спутниковой навигационной системы |
CN106767787A (zh) * | 2016-12-29 | 2017-05-31 | 北京时代民芯科技有限公司 | 一种紧耦合gnss/ins组合导航装置 |
CN107037469A (zh) * | 2017-04-11 | 2017-08-11 | 北京七维航测科技股份有限公司 | 基于安装参数自校准的双天线组合惯导装置 |
CN108180925B (zh) * | 2017-12-15 | 2020-09-01 | 中国船舶重工集团公司第七0七研究所 | 一种里程计辅助车载动态对准方法 |
CN109459044B (zh) * | 2018-12-17 | 2022-09-06 | 北京计算机技术及应用研究所 | 一种gnss双天线辅助的车载mems惯导组合导航方法 |
-
2019
- 2019-04-18 CN CN201910314091.9A patent/CN110133695B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110133695A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110133695B (zh) | 一种双天线gnss位置延迟时间动态估计系统及方法 | |
CN110133694B (zh) | 基于双天线gnss航向和轮速辅助的车辆定位方法及系统 | |
CN110779521A (zh) | 一种多源融合的高精度定位方法与装置 | |
CN106840179B (zh) | 一种基于多传感器信息融合的智能车定位方法 | |
US8195392B2 (en) | Position detecting apparatus and method used in navigation system | |
US8510044B2 (en) | Position sensing device and method | |
US11079494B2 (en) | Positioning device | |
JP5419665B2 (ja) | 位置標定装置、位置標定方法、位置標定プログラム、速度ベクトル算出装置、速度ベクトル算出方法および速度ベクトル算出プログラム | |
JP5586994B2 (ja) | 位置標定装置、位置標定装置の位置標定方法および位置標定プログラム | |
CN109343095B (zh) | 一种车载导航车辆组合定位装置及其组合定位方法 | |
CN109917440B (zh) | 一种组合导航方法、系统及车辆 | |
JP6409346B2 (ja) | 移動距離推定装置 | |
US20100076681A1 (en) | Dead reckoning system | |
CA2989529C (en) | Positioning apparatus | |
WO2015104757A1 (ja) | 速度推定装置 | |
CN102508280B (zh) | 陀螺辅助双天线测量单元确定整周模糊度和航向的方法 | |
CN109470276B (zh) | 基于零速修正的里程计标定方法与装置 | |
CN111751857A (zh) | 一种车辆位姿的估算方法、装置、存储介质及系统 | |
CN105606093B (zh) | 基于重力实时补偿的惯性导航方法及装置 | |
JP3218876B2 (ja) | 車両用現在位置検出装置 | |
CN116337053A (zh) | 车辆导航方法、装置、电子设备及存储介质 | |
JP6531768B2 (ja) | センサ誤差補正装置及び方法 | |
US12019171B2 (en) | Method for determining the position of a vehicle as a function of the vehicle velocity | |
JP3628046B2 (ja) | 車両用現在位置検出装置 | |
CN113048987A (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 |