CN107036627A - 环形激光测角装置的地速投影分量误差自校准方法 - Google Patents
环形激光测角装置的地速投影分量误差自校准方法 Download PDFInfo
- Publication number
- CN107036627A CN107036627A CN201710204728.XA CN201710204728A CN107036627A CN 107036627 A CN107036627 A CN 107036627A CN 201710204728 A CN201710204728 A CN 201710204728A CN 107036627 A CN107036627 A CN 107036627A
- Authority
- CN
- China
- Prior art keywords
- angle
- measuring equipment
- loop laser
- error
- laser angle
- 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
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于环形激光测角技术,涉及一种环形激光测角装置的地速投影分量误差自校准方法。本发明环形激光测角装置的地速投影分量误差自校准方法采用最小二乘法,将来自于被测转台上所安装的圆感应同步器的角度信息与环形激光测角装置的角度信息进行融合,并通过离散化处理分离出误差项,从而消除因地速投影分量所引入的测角误差。本发明充分利用了环形激光测角装置短时间内精度高及转台角度误差不随时间发散的优点,可以成功地降低地速投影分量所引入的测角误差,得到精度更高、噪声更小的角度及角速度信息。
Description
技术领域
本发明属于环形激光测角技术,涉及一种环形激光测角装置的地速投影分量误差自校准方法。
背景技术
角度测量是几何量计量技术的重要组成部分,特别是高精度动态角度标定,例如在数控机床、高精度转台、卫星天线、远程望远镜等领域,都存在着既要求静态定位精度,也要求稳定的转动动态特性,也就是在相对较高速率的转动下其角位置精度的问题。此前,对转台特性的高精度测量通常会选择使用光学码盘或基于自准直仪的方法,但前者存在着精度有限,无法对角度转动的动态特性进行测量,而后者则存在测角范围小,无法整周测量等缺点。此外,目前的大多数角度测量方法无法实现很高数据更新率(例如,105Hz量级)的角增量输出,对深入研究转台控制系统的动态特性会有一定的影响。
在实际角度测量过程中,因为机械加工及机械安装精度的限制,环形激光测角装置的敏感轴与被测转台之间通常会存在数角分的夹角。而环形激光测角装置所感应的是相对于惯性坐标系中的角速度,因此,上述夹角会将地球自转角速度引入到环形激光测角装置的最终输出中,从而产生测角误差。
发明内容
本发明的目的:提供一种用于消除环形激光测角装置中的因地速投影分量所引入的测角误差的自校准方法。
本发明的技术方案:环形激光测角装置的地速投影分量误差自校准方法,其硬件部分包括环形激光测角装置、圆感应同步器、过渡板、同步采样电路及处理器,其中,环形激光测角装置和圆感应同步器分别装在过渡板两侧,且二者均与同步采样电路连接,而同步采样电路与处理器相连,其采用最小二乘法,将来自于被测转台上所安装的圆感应同步器的角度信息与环形激光测角装置的角度信息进行融合,并通过离散化处理分离出误差项,从而消除因地速投影分量所引入的测角误差。
所述的环形激光测角装置的地速投影分量误差自校准方法,其具体步骤如下:
步骤1:高速同步采集
由同步采样电路高速同步采集环形激光测角装置的角速度信号和圆感应同步器的角位置信号;
步骤2:误差分离
根据环形激光测角装置的相对敏感的为惯性空间的角速度,而圆感应同步器相对敏感的角位置为相对于转台自身的角位置的特征,得到环形激光测角装置的角速度信号输出
ωg=ωNsinθcosψsinΔψ+ωNsinθsinψcosΔψ+(ωr+ωU)cosθ+εg (1.6)
其中,ωg为环形激光测角装置输出,ψ为转台转角,ωN为地球自转角速度ωie在地理正北方向上的投影量,ωU为地球自转角速度ωie在地理垂线的反方向上的投影量,ωN=ωiecos L,ωU=ωiesin L,L为当地纬度;θ是被测转台的转台与环形激光测角装置之间的夹角;Δψ是环形激光测角装置的敏感轴在安装平面上的投影与地理正东方向之间的夹角;ωr是圆感应同步器的角速度,由圆感应同步器的角位置信息进行运算后得到;εg是环形激光测角装置自身的误差;
当外部同步时钟的周期为ΔT,对(1.1)式进行离散化,得到
ωg(k)=ωNsinθcosψ(k)sinΔψ+ωNsinθsinψ(k)cosΔψ+(ωr(k)+ωU)cosθ+εg
(1.7)
对(1.2)式进行矢量化,得到
利用最小二乘法求解出以下变量的最优估计值:(ωNsinθsinΔψ),(ωNsinθcosΔψ),(cosθ)和(ωUcosθ+εg),至此实现对环形激光测角装置的角速度信号输出的误差分离;
步骤3:误差消除
根据步骤2所计算得到的cosθ,计算出圆感应同步器的转台(12)与环形激光测角装置(13)之间的夹角θ,根据计算得到的θ,就能够根据下式对环形激光测角仪中因地速投影分量所引入的测角误差进行消除,得到修正后的、更准确的角度输出信息。
其中,ωr′(k)是修正后的第k个时钟周期内的角速度信息,对从m·ΔT到n·ΔT时间段内的ωr′(k)进行累加即可得到该时间段内的角度增量信息。
所述最小二乘法包括常规最小二乘法或加权最小二乘法或递推最小二乘法。
本发明的优点和有益效果是:本发明地速投影分量误差自校准方法,能够将来自于被测转台上所安装的圆感应同步器的角度信息与环形激光测角装置的角度信息进行融合,充分利用了环形激光测角装置短时间内精度高及转台角度误差不随时间发散的优点,采用最小二乘法,可以成功地降低地速投影分量给测角过程所引入的误差,得到精度更高、噪声更小的角度及角速度信息。该方法能够实现自标校,减少了人工干预,提高了产品可靠性和易用性。
附图说明
图1是环形激光测角装置及其安装示意图;
其中,1-环形激光测角装置,2-过渡板,3-圆感应同步器(与被测转台同心安装),4-同步采样电路,5-处理器(可以是数字处理器芯片,包括单片机、现场可编程门阵列,也可以是计算机),11-被测转台;
图2是环形激光测角装置敏感轴与被测转台之间夹角的示意图,
其中,11-被测转台,12-被测转台的转台,13-环形激光测角装置的敏感轴,θ是12与13之间的夹角。
图3是环形激光测角装置敏感轴在水平面内的投影(s)与北向(N)和东向(E)之间的夹角关系示意图,其中Δψ是s与东向的夹角。
具体实施方式
下面结合附图对本发明作进一步的说明:
请同时参阅图1、图2、图3,本发明环形激光测角装置的地速投影分量误差自校准方法的硬件部分包括环形激光测角装置、圆感应同步器、过渡板、同步采样电路及处理器,其中,环形激光测角装置和圆感应同步器分别装在过渡板两侧,且二者均与同步采样电路连接,而同步采样电路与处理器相连。本发明环形激光测角装置的地速投影分量误差自校准方法采用最小二乘法,将来自于被测转台上所安装的圆感应同步器的角度信息与环形激光测角装置的角度信息进行融合,并通过离散化处理分离出误差项,从而消除因地速投影分量所引入的测角误差。
所述的环形激光测角装置的地速投影分量误差自校准方法的具体实施过程如下:
步骤1:高速同步采集
由同步采样电路高速同步采集环形激光测角装置的角速度信号和圆感应同步器的角位置信号;
同步采样电路通常采用现场可编程门阵列器件实现,外同步信号通常采用TTL电平的数字脉冲信号,外同步信号的频率通常在1MHz左右;
步骤2:误差分离
根据环形激光测角装置的相对敏感的为惯性空间的角速度,而圆感应同步器相对敏感的角位置为相对于转台自身的角位置的特征,得到环形激光测角装置的角速度信号输出
ωg=ωNsinθcosψsinΔψ+ωNsinθsinψcosΔψ+(ωr+ωU)cosθ+εg (1.11)
其中,ωg为环形激光测角装置输出,ψ为转台转角,ωN为地球自转角速度ωie在地理正北方向上的投影量,ωU为地球自转角速度ωie在地理垂线的反方向上的投影量,ωN=ωiecos L,ωU=ωiesin L,L为当地纬度;θ是被测转台的转台与环形激光测角装置之间的夹角;Δψ是环形激光测角装置的敏感轴在安装平面上的投影与地理正东方向之间的夹角;ωr是圆感应同步器的角速度,由圆感应同步器的角位置信息进行运算后得到;εg是环形激光测角装置的量测误差;
当外部同步时钟的周期为ΔT,对(1.1)式进行离散化,得到
ωg(k)=ωNsinθcosψ(k)sinΔψ+ωNsinθsinψ(k)cosΔψ+(ωr(k)+ωU)cosθ+εg
(1.12)
对(1.2)式进行矢量化,得到
利用最小二乘法求解出以下变量的最优估计值:(ωNsinθsinΔψ),(ωNsinθcosΔψ),(cosθ)和(ωUcosθ+εg),从而实现对环形激光测角装置的角速度信号输出的误差分离;
步骤3:误差消除
根据cosθ精确地计算出被测转台的转台12与环形激光测角装置13之间的夹角θ,根据计算得到的θ,就能够根据式(1.2)对环形激光测角仪中因地速投影分量所引入的测角误差进行消除,得到修正后的、更准确的角速度输出信息。
其中,ωr′(k)是修正后的第k个时钟周期内的角速度信息,对从m·ΔT到n·ΔT时间段内的ωr′(k)进行累加即可得到该时间段内的角度增量信息。
另外,本发明环形激光测角装置的地速投影分量误差自校准方法中的最小二乘法可以为常规最小二乘法或加权最小二乘法或递推最小二乘法。
本发明环形激光测角装置的地速投影分量误差自校准方法能够有效地消除因环形激光测角装置的敏感轴与被测转台之间的夹角所引入的地速投影分量测角误差。相比于传统的地速投影分量误差消除算法,本发明所介绍的方法可以将地速分量投影误差、环形激光测角装置的自身误差进行分离,简单易用,且精度更高。
Claims (4)
1.环形激光测角装置的地速投影分量误差自校准方法,其硬件部分包括环形激光测角装置、圆感应同步器、过渡板、同步采样电路及处理器,其中,环形激光测角装置和圆感应同步器分别装在过渡板两侧,且二者均与同步采样电路连接,而同步采样电路与处理器相连,其特征在于,采用最小二乘法,将来自于被测转台上所安装的圆感应同步器的角度信息与环形激光测角装置的角度信息进行融合,并通过离散化处理分离出误差项,从而消除因地速投影分量所引入的测角误差。
2.根据权利要求1所述的环形激光测角装置的地速投影分量误差自校准方法,其特征在于,具体步骤如下:
步骤1:高速同步采集
由同步采样电路高速同步采集环形激光测角装置的角速度信号和圆感应同步器的角位置信号;
步骤2:误差分离
根据环形激光测角装置的相对敏感的为惯性空间的角速度,而圆感应同步器相对敏感的角位置为相对于转台自身的角位置的特征,得到环形激光测角装置的角速度信号输出
ωg=ωNsinθcosψsinΔψ+ωNsinθsinψcosΔψ+(ωr+ωU)cosθ+εg (1.1)
其中,ωg为环形激光测角装置输出,ψ为转台转角,ωN为地球自转角速度ωie在地理正北方向上的投影量,ωU为地球自转角速度ωie在地理垂线的反方向上的投影量,ωN=ωiecosL,ωU=ωiesinL,L为当地纬度;θ是被测转台的转台与环形激光测角装置之间的夹角;Δψ是环形激光测角装置的敏感轴在安装平面上的投影与地理正东方向之间的夹角;ωr是圆感应同步器测得的角速度,由圆感应同步器的角位置信息进行运算后得到;εg是环形激光测角装置自身的误差;
当外部同步时钟的周期为ΔT,对(1.1)式进行离散化,得到
ωg(k)=ωNsinθcosψ(k)sinΔψ+ωNsinθsinψ(k)cosΔψ+(ωr(k)+ωU)cosθ+εg
(1.2)
对(1.2)式进行矢量化,得到
利用最小二乘法求解出以下变量的最优估计值:(ωNsinθsinΔψ),(ωNsinθcosΔψ),(cosθ)和(ωUcosθ+εg),至此实现对环形激光测角装置的角速度信号输出的误差分离;
步骤3:误差消除
根据步骤2所计算得到的cosθ,计算出圆感应同步器的转台(12)与环形激光测角装置(13)之间的夹角θ,根据计算得到的θ,就能够根据下式对环形激光测角仪中因地速投影分量所引入的测角误差进行消除,得到修正后的、更准确的角度输出信息。
3.根据权利要求2所述的环形激光测角装置的地速投影分量误差自校准方法,其特征在于,步骤3中的角度输出信息中的角速度信号输出:
其中,ωr′(k)是修正后的第k个时钟周期内的角速度信息,对从m·ΔT到n·ΔT时间段内的ωr′(k)进行累加即可得到该时间段内的角度增量信息:
4.根据权利要求2所述的环形激光测角装置的地速投影分量误差自校准方法,其特征在于,所述最小二乘法包括常规最小二乘法或加权最小二乘法或递推最小二乘法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710204728.XA CN107036627B (zh) | 2017-03-30 | 2017-03-30 | 环形激光测角装置的地速投影分量误差自校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710204728.XA CN107036627B (zh) | 2017-03-30 | 2017-03-30 | 环形激光测角装置的地速投影分量误差自校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107036627A true CN107036627A (zh) | 2017-08-11 |
CN107036627B CN107036627B (zh) | 2020-06-30 |
Family
ID=59534186
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710204728.XA Active CN107036627B (zh) | 2017-03-30 | 2017-03-30 | 环形激光测角装置的地速投影分量误差自校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107036627B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108444433A (zh) * | 2018-03-07 | 2018-08-24 | 太原理工大学 | 基于面型基准的转台转角误差检测装置及方法 |
CN111336952A (zh) * | 2020-02-26 | 2020-06-26 | 北京空间机电研究所 | 一种感应同步器角度修正的系统、方法及装置 |
CN113624255A (zh) * | 2021-07-30 | 2021-11-09 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种环形激光测角仪敏感轴修正装置及方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050235507A1 (en) * | 2004-04-23 | 2005-10-27 | Badami Vivek G | Measurement of form of spherical and near-spherical optical surfaces |
CN101509785A (zh) * | 2009-03-16 | 2009-08-19 | 浙江大学 | 一种光纤陀螺仪输入轴失准角评测方法 |
CN102128645A (zh) * | 2010-12-13 | 2011-07-20 | 中国科学院光电技术研究所 | 用于光电轴角编码器动态测角误差检测仪 |
CN102175207A (zh) * | 2011-03-01 | 2011-09-07 | 东华大学 | 一种基于反向和多传感器法结合的高精度圆度检测方法 |
US20120222465A1 (en) * | 2009-11-26 | 2012-09-06 | Leica Geosystems Ag | Calibration method and angle measuring method for an angle measuring device, and angle measuring device |
CN102985793A (zh) * | 2010-05-17 | 2013-03-20 | 法罗技术股份有限公司 | 自补偿角度编码器 |
CN103017791A (zh) * | 2012-12-12 | 2013-04-03 | 中国科学院长春光学精密机械与物理研究所 | 高精度判别测角仪器的误差源的方法 |
CN104296644A (zh) * | 2014-09-26 | 2015-01-21 | 中国科学院光电研究院 | 一种测角误差校准装置及其安装调节方法 |
CN105026885A (zh) * | 2013-01-09 | 2015-11-04 | 卡尔蔡司工业测量技术有限公司 | 用于求取旋转装置旋转误差的检验体 |
CN105091835A (zh) * | 2014-05-07 | 2015-11-25 | 中国航空工业第六一八研究所 | 一种转台动态角位置误差测量精度提高装置及校正方法 |
CN105091844A (zh) * | 2014-05-07 | 2015-11-25 | 中国航空工业第六一八研究所 | 一种高精度动态测角装置及方法 |
CN105091789A (zh) * | 2014-05-07 | 2015-11-25 | 中国航空工业第六一八研究所 | 一种基于空间四频差动激光陀螺的高精度测角装置及其安装标定方法 |
CN106482668A (zh) * | 2016-11-21 | 2017-03-08 | 班戈设备系统(苏州)有限公司 | 一种角度高精确度测量系统及测角方法 |
-
2017
- 2017-03-30 CN CN201710204728.XA patent/CN107036627B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050235507A1 (en) * | 2004-04-23 | 2005-10-27 | Badami Vivek G | Measurement of form of spherical and near-spherical optical surfaces |
CN101509785A (zh) * | 2009-03-16 | 2009-08-19 | 浙江大学 | 一种光纤陀螺仪输入轴失准角评测方法 |
US20120222465A1 (en) * | 2009-11-26 | 2012-09-06 | Leica Geosystems Ag | Calibration method and angle measuring method for an angle measuring device, and angle measuring device |
CN102985793A (zh) * | 2010-05-17 | 2013-03-20 | 法罗技术股份有限公司 | 自补偿角度编码器 |
CN102128645A (zh) * | 2010-12-13 | 2011-07-20 | 中国科学院光电技术研究所 | 用于光电轴角编码器动态测角误差检测仪 |
CN102175207A (zh) * | 2011-03-01 | 2011-09-07 | 东华大学 | 一种基于反向和多传感器法结合的高精度圆度检测方法 |
CN103017791A (zh) * | 2012-12-12 | 2013-04-03 | 中国科学院长春光学精密机械与物理研究所 | 高精度判别测角仪器的误差源的方法 |
CN105026885A (zh) * | 2013-01-09 | 2015-11-04 | 卡尔蔡司工业测量技术有限公司 | 用于求取旋转装置旋转误差的检验体 |
CN105091835A (zh) * | 2014-05-07 | 2015-11-25 | 中国航空工业第六一八研究所 | 一种转台动态角位置误差测量精度提高装置及校正方法 |
CN105091844A (zh) * | 2014-05-07 | 2015-11-25 | 中国航空工业第六一八研究所 | 一种高精度动态测角装置及方法 |
CN105091789A (zh) * | 2014-05-07 | 2015-11-25 | 中国航空工业第六一八研究所 | 一种基于空间四频差动激光陀螺的高精度测角装置及其安装标定方法 |
CN104296644A (zh) * | 2014-09-26 | 2015-01-21 | 中国科学院光电研究院 | 一种测角误差校准装置及其安装调节方法 |
CN106482668A (zh) * | 2016-11-21 | 2017-03-08 | 班戈设备系统(苏州)有限公司 | 一种角度高精确度测量系统及测角方法 |
Non-Patent Citations (8)
Title |
---|
SWAMINATHAN, RAHUL等: ""Nonmetric calibration of wide-angle lenses and polycameras"", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》 * |
XUEJUN ZHANG等: ""Gain and phase error-free LINC transmitter"", 《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》 * |
丁国清: ""基于误差分离技术的超精密测量及校正方法研究"", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
任顺清等: ""圆感应同步器测角误差的分离技术"", 《中国电机工程学报》 * |
何征宇等: ""一种基于CMOS图像传感器的大视场角度测量装置的误差校正方法"", 《光学技术》 * |
傅鑫等: ""动态角度的精密测量技术"", 《第十八届中国科协年会-分6 军民融合高端论坛论文集》 * |
刘学峰: ""转台校准设备研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
吕沧海,等: "《国防工业出版社》", 31 March 2014 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108444433A (zh) * | 2018-03-07 | 2018-08-24 | 太原理工大学 | 基于面型基准的转台转角误差检测装置及方法 |
CN111336952A (zh) * | 2020-02-26 | 2020-06-26 | 北京空间机电研究所 | 一种感应同步器角度修正的系统、方法及装置 |
CN113624255A (zh) * | 2021-07-30 | 2021-11-09 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种环形激光测角仪敏感轴修正装置及方法 |
CN113624255B (zh) * | 2021-07-30 | 2023-09-05 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种环形激光测角仪敏感轴修正装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107036627B (zh) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11409001B2 (en) | Method for tilt measurement and compensation of surveying instrument based on GNSS receiver and IMU sensor | |
US10704902B2 (en) | Surveying pole | |
US20120245850A1 (en) | Azimuth initialization and calibration of wellbore surveying gyroscopic and inertial instruments by means of an external navigation system | |
CN102168990B (zh) | 惯性定向设备的高精度检测标定装置及其检测标定方法 | |
CN109459008B (zh) | 一种小型中高精度光纤陀螺寻北装置及方法 | |
CN106772493A (zh) | 基于北斗差分定位的无人机航向测算系统及其测算方法 | |
CN107036627A (zh) | 环形激光测角装置的地速投影分量误差自校准方法 | |
US11226203B2 (en) | Low cost INS | |
CA2738665A1 (en) | Method and apparatus for precision azimuth measurement | |
CN201955097U (zh) | 惯性定向设备的检测标定系统 | |
CN111504288B (zh) | 一种寻北装置 | |
Zaitsev et al. | Study of systems error compensation methods based on molecular-electronic transducers of motion parameters | |
CN109470274B (zh) | 一种车载光电经纬仪载车平台变形测量系统及方法 | |
Fang et al. | A new inclination error calibration method of motion table based on accelerometers | |
CN201983799U (zh) | 一种惯性定向设备的高精度检测标定装置 | |
CN207528248U (zh) | 一种星载二维指向机构测角精度检测装置 | |
CN109681189A (zh) | 一种井径扇区固井质量及轨迹一体化测量仪 | |
CN108775899A (zh) | 基于伪卫星和惯性信息的深部开采井上下坐标系连接方法 | |
CN102410845B (zh) | 误差校正、角速度检测及鼠标控制方法和装置,空间鼠标 | |
CN108803373B (zh) | 一种三轴转台的地速消除方法 | |
CN109470275A (zh) | 一种机动布站的光电经纬仪高精度自主定向方法 | |
CN107179066B (zh) | 旋转定标装置及其定标方法 | |
CN105044757A (zh) | 卫星信号遮蔽区域gnss差分与惯性测量组合测图方法 | |
CN112461268A (zh) | 天向陀螺的代位标定方法与双陀螺式导航设备的标定方法 | |
CN206756137U (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 |