CN102425990B - 工作空间测量定位系统网络状态监测方法 - Google Patents
工作空间测量定位系统网络状态监测方法 Download PDFInfo
- Publication number
- CN102425990B CN102425990B CN201110270929.2A CN201110270929A CN102425990B CN 102425990 B CN102425990 B CN 102425990B CN 201110270929 A CN201110270929 A CN 201110270929A CN 102425990 B CN102425990 B CN 102425990B
- Authority
- CN
- China
- Prior art keywords
- theta
- cell site
- prime
- receiver
- monitoring
- 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
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及工业现场大尺寸三维坐标测量方法。为提供一种工作空间测量定位系统网络状态监测方法,本发明采取的技术方案是,工作空间测量定位系统网络状态监测方法:包括下列步骤:1、布置好发射站;2、在工作空间内安装两个以上用于监控发射站状态的接收器;3、WMPS测量网络定向完成后,可通过状态监控接收器采集发射站空间角度初始数据θ10或θ20并储存在上位机内,工作时每隔一定监控周期接收器多次连续采集发射站的扫描角度θ1、θ2并与上位机存储的初始数据对比、判定发射站工作状态;4、具体发射站工作状态采用测量结果统计值考核。本发明主要应用于三座标测量场合。
Description
技术领域
本发明涉及工业现场大尺寸三维坐标测量方法,特别涉及一种采用光信号接收器的工作空间测量定位系统网络状态监测方法,具体讲,涉及工作空间测量定位系统网络状态监测方法。
背景技术
本发明所涉及的方法是基于论文《扫描平面激光空间定位系统测量网络的构建》中所述的工作空间测量定位系统(WMPS:Workspace Measurement Positioning System),该系统为针对大型制造业测量需求及全局测量控制网的特点发展一种新型多站网络式室内测量定位系统,可实现工业现场大尺度空间坐标的高精度自动测量。如图1所示,此类系统借鉴了全球定位系统的测量思想,使用多个激光发射装置(发射站)组成测量网络,采用基于光电扫描的空间角度交会自动测量方法对单个接收器进行定位,发射站在工作时不负责解算接收器坐标,而是通过向外发射带有角度信息的光信号,为测量空间内的光电接收器提供定位服务。激光发射站结构如图2所示,发射站工作时转台绕固定轴匀速旋转并向外发射两束随转台一同旋转的平面激光,同时每当转台转至一个预定位置时基座上的激光器发出全向光脉冲作为单周旋转起点的同步标记。接收器接收到同步标记光信号及扫描平面光信号并通过内部计时器记录此时时间值,解算发射站转过角度。处理电路在获得发射站角度信息后通过无线网络将角度值发送给解算工作站。在已知发射站方位信息后,只要有两个以上基站的角度测量值,解算工作站便可以使用角度交汇方法计算出此时接收器的精确坐标。这种发射站网络单向广播的坐标测量方式使WMPS系统可同时对多个接收器进行跟踪定位并为现场的自动化组装设备(如机械手、自动引导小车等)提供导航信息,非常适于工业应用。
由测量原理可知,欲使WMPS系统测量结果稳定精确不但需要保证网络内发射站转速稳定在某一固定值,而且还要保证其位置姿(坐标及方向角)固定不变。虽然在现场使用时发射站一般安装在稳定的地基或墙壁侧立面上,但是大型装配设备长时间运动产生的震动和厂房框架受温度影响产生的变形仍会不可避免的使某些发射站位姿发生轻微改变,导致测量精度降低。同时,由于转台轴系本身寿命限制网络中个别长时间使用的发射站旋转稳定性降低同样会影响系统测量精度。鉴于WMPS系统的单向广播式测量方式,欲保证测量网络的稳定可靠需要解决以下两方面关键技术问题:
1、及时获得网络中各个发射站的工作状态并发现有问题的设备,对测量网络的工作装态实施有效监控。
2、对发生位姿改变但仍能使用的发射站通过辅助设施实现位置及姿态参数的在线补偿,实现可自愈的测量网络,提高系统的可靠性。
发明内容
为克服现有技术的不足,提供一种工作空间测量定位系统网络状态监测方法,本发明采取的技术方案是,工作空间测量定位系统网络状态监测方法:包括下列步骤:
1、布置好发射站,实现发射站全局定向,完成WMPS测量网络的构建;
2、在工作空间内安装两个以上用于监控发射站状态的接收器,并使每个接收器可接收到该测量区域内的所有发射站信息,用于安装接收器的支架由殷钢或陶瓷等温度系数较小的材料制成并固连在工位地基上;
3、WMPS测量网络定向完成后,可通过状态监控接收器采集发射站空间角度初始数据θ10或θ20并储存在上位机内,工作时每隔一定监控周期接收器多次连续采集发射站的扫描角度θ1、θ2并与上位机存储的初始数据对比、判定发射站工作状态;
4、具体发射站工作状态采用以下测量结果统计值考核:
(1)单个监控周期内的θ1或θ2的多次测量结果标准差:
该结果反映发射站轴系运行状态,当标准差大于某一阈值时,可被判断为发射站转速失稳,轴系结构存在问题;
(2)多个监控周期内的θ1或θ2的多次测量结果均值与初始数据θ10或θ20之差:
该结果反映发射站与接收器间空间位置关系,当δm变化超过某一预订阈值即可判断为发射站与监控接收器间的相对位姿发生变化,工作现场中当多个接收器测量同一发射站得到的统计量δ1、δ2同时发生变化时,可判定为发射站发生旋转或位移,需对发射站重新定向;另一方面如果现场只有一个监控接收器测量结果超阈值,其他结果保持在阈值内时,可判定为该接收器坐标变化,需要重新固定接收器,并更新空间角度初始数据测量值θ10或θ20。
所述采集发射站的扫描角度θ1、θ2是,
1、激光发射站的数学模型被抽象为绕公共转轴上一点的两个旋转的非平行半平面及一个以固定频率发射脉冲光的点光源,而接收器简化为以靶面光学中心为中心的质点,组装完毕后定义发射站坐标系可定义如下,以激光平面1与转轴的交点为激光发射站原点,旋转轴为Z轴,X轴为初始时刻,即转台转至固定位置,发射器发射脉冲光时,激光器1光轴所在位置,Y轴遵循右手定则,出厂前,对发射站旋转角速度ω进行设定,并标定发射站激光平面在初始位置处的平面ABCD参数,则发射站坐标系下光平面方程可表示为:
a′1、a′2、b′1、b′2、c′1、c′2、d′2、d′2即为ABCD参数,当两平面激光分别扫过接收器时其在发射站坐标系下的两平面方程系数变为转台旋转角度θ1、θ2的函数:
在经过标定建立测量坐标系后,每个发射站在测量坐标系下位姿已知,即确定了每台激光发射站原点在全局坐标系下的坐标PTXG及激光发射站坐标系到测量坐标系的旋转矩阵RTXG,发射站的两个扫描光分别到达接收器时,发射站转台转过角度分别为θ1、θ2时易知发射站旋两个转光平面方程在测量坐标系下的方程参数(aTXG1(θ1),bTXG1(θ1),cTXG1(θ1),dTXG1(θ1))T及(aTXG2(θ2),bTXG2(θ2),cTXG2(θ2),dTXG2(θ2))T可表示为:
2、当发射站光平面扫过接收器时,由于接收器在平面内,方程(1)在测量坐标系下仍然成立,如此时接收器在测量坐标系下坐标为(xRXG yRXG zRXG),可得:
观察(5)式可知,当发射站制造完成后其光平面参数(a′m b′m c′m d′m)在工作现场一般不会产生变化,此时,如果接收器坐标(xRXG yRXG zRXG)被视为定值,其测量到的转台旋转角度θ1、θ2直接反映了发射站坐标PTXG旋转矩阵RTXG的变化情况,T表示转置。
本发明具有如下技术效果:
本发明充分利用WMPS系统接收器的计时测角原理,将其作为状态监控用传感器,通过在测量空间内布置监控传感器实时监控网络内发射站轴系状况及位姿变化情况而无需解算传感器坐标,方法简单实用同时增加了测量系统的可靠性。此外,本发明还具有以下优势:
1、监控过程中无需解算接收器坐标,计算负担较小,方法简单有效;
2、监控过程与测量过程同时进行,互不干扰;
3、可及时获得网络中各个发射站的工作状态并发现有问题的设备,极大增加了测量系统的可靠性。
附图说明
图1工作空间测量定位系统组成示意图。图中,1为发射站,2为接收器。
图2激光发射站结构示意图。图中,3为旋转平台,4为线激光器,5为脉冲光源,6为光栅码盘,7为电机驱动,8为固定基座,9为激光平面1,10为激光平面2。
图3激光发射数学模型。
图4通过接收器监控发射站状态。图中,12为测量用接收器,13为监测用接收器,14为接收器支架,15为地基。
具体实施方式
本发明的技术方案如下:
1、激光发射站的数学模型如图3所示,可被抽象为绕公共转轴上一点的两个旋转的非平行半平面及一个以固定频率发射脉冲光的点光源,而接收器可简化为以靶面光学中心为中心的质点。组装完毕后定义发射站坐标系可定义如下,以激光平面1与转轴的交点为激光发射站原点,旋转轴为Z轴,X轴为初始时刻(即转台转至固定位置,发射器发射脉冲光时)激光器1光轴所在位置,Y轴遵循右手定则。出厂前,对发射站旋转角速度ω进行设定,并标定发射站激光平面在初始位置处的平面ABCD参数,则发射站坐标系下光平面方程可表示为:
当两平面激光分别扫过接收器时其在发射站坐标系下的两平面方程系数变为转台旋转角度θ1、θ2的函数:
在经过标定建立测量坐标系后,每个发射站在测量坐标系下位姿已知,即确定了每台激光发射站原点在全局坐标系下的坐标PTXG及激光发射站坐标系到测量坐标系的旋转矩阵RTXG。发射站的两个扫描光分别到达接收器时,发射站转台转过角度分别为θ1、θ2时易知发射站旋两个转光平面方程在测量坐标系下的方程参数(aTXG1(θ1),bTXG1(θ1),cTXG1(θ1),dTXG1(θ1))T及(aTXG2(θ2),bTXG2(θ2),cTXG2(θ2),dTXG2(θ2))T可表示为:
2、当发射站光平面扫过接收器时,由于接收器在平面内,方程(1)在测量坐标系下仍然成立,如此时接收器在测量坐标系下坐标为(xRXG yRXG zRXG),可得:
观察(5)式可知,当发射站制造完成后其光平面参数(a′m b′m c′m d′m)在工作现场一般不会产生变化。此时,如果接收器坐标(xRXG yRXG zRXG)被视为定值,其测量到的转台旋转角度θ1、θ2直接反映了发射站坐标PTXG旋转矩阵RTXG的变化情况,因此无需解算接收器坐标,可通过在现场每个一段时间采集发射站的旋转角度值直接监测发射站位姿变化。T的含义是转置,例如
3、在工作空间内安装两个以上用于监控发射站状态的接收器,并使每个接收器可接收到该测量区域内的所有发射站信息,如图4所示。用于安装接收器的支架由殷钢或陶瓷等温度系数较小、刚性较好不易变形的材料制成并固连在工位地基上,因此接收器在工作时可被认为静止。
4、在WMPS系统定向完成后,可通过状态监控接收器采集发射站空间角度初始数据θ10或θ20并储存在上位机内,系统工作时每隔一定监控周期Tmonitor接收器多次连续采集发射站的扫描角度θ1、θ2并与上位机存储的初始数据对比、判定发射站工作状态。具体发射站工作状态采用以下测量结果统计值考核:
(1)单个监控周期内的θ1或θ2的多次测量结果标准差:
该结果反映发射站轴系运行状态,当标准差大于某一阈值时,可被判断为发射站转速失稳,轴系结构存在问题。此时,需更换发射站。
(2)多个监控周期内的θ1或θ2的多次测量结果均值与初始数据θ10或θ20之差:
该结果反映发射站与接收器间空间位置关系,当δm变化超过某一预订阈值即可判断为发射站与监控接收器间的相对位姿发生变化。工作现场中当多个接收器测量同一发射站得到的统计量δ1、δ2同时发生变化时,可判定为发射站发生旋转或位移,需对发射站重新定向;另一方面如果现场只有一个监控接收器测量结果超阈值,其他结果保持在阈值内时,可判定为该接收器坐标变化,需要重新固定接收器,并更新空间角度初始数据测量值θ10、θ20。
本发明充分利用WMPS系统接收器的计时测角原理,将其作为状态监控用传感器,通过在测量空间内布置监控用冗余接收器实现了测量网络内发射站轴系状况及位姿变化情况的实时监控。具有以下优势:
1、监控过程中无需解算接收器坐标,计算负担较小,方法简单有效。
2、监控过程与测量过程同时进行,互不干扰。
3、可及时获得网络中各个发射站的工作状态并发现有问题的设备,极大增加了测量系统的可靠性。
下面结合附图和具体实施方式进一步详细说明本发明。
1、布置好发射站,采用《扫描平面激光空间定位系统测量网络的构建》中所述发射站定向方法实现发射站全局定向,完成WMPS测量网络的构建。
2、在工作空间内安装两个以上用于监控发射站状态的接收器,并使每个接收器可接收到该测量区域内的所有发射站信息,如图4所示。用于安装接收器的支架由殷钢或陶瓷等温度系数较小的材料制成并固连在工位地基上,因此接收器在工作时可被认为静止。
3、WMPS系统定向完成后,可通过状态监控接收器采集发射站空间角度初始数据θ10或θ20并储存在上位机内,系统工作时每隔一定监控周期Tmonitor接收器多次连续采集发射站的扫描角度θ1、θ2并与上位机存储的初始数据对比、判定发射站工作状态。
4、具体发射站工作状态采用以下测量结果统计值考核:
(1)单个监控周期内的θ1或θ2的多次测量结果标准差:
该结果反映发射站轴系运行状态,当标准差大于某一阈值时,可被判断为发射站转速失稳,轴系结构存在问题。此时,需更换发射站。
(2)多个监控周期内的θ1或θ2的多次测量结果均值与初始数据θ10或θ20之差:
该结果反映发射站与接收器间空间位置关系,当δm变化超过某一预订阈值即可判断为发射站与监控接收器间的相对位姿发生变化。工作现场中当多个接收器测量同一发射站得到的统计量δ1、δ2同时发生变化时,可判定为发射站发生旋转或位移,需对发射站重新定向;另一方面如果现场只有一个监控接收器测量结果超阈值,其他结果保持在阈值内时,可判定为该接收器坐标变化,需要重新固定接收器,并更新空间角度初始数据测量值θ10、θ20。
Claims (2)
1.一种工作空间测量定位系统网络状态监测方法,其特征是,包括如下步骤:
1)、布置好发射站,实现发射站全局定向,完成WMPS测量网络的构建;
2)、在工作空间内安装两个以上用于监控发射站状态的接收器,并使每个接收器可接收到该测量区域内的所有发射站信息,用于安装接收器的支架由温度系数较小的材料殷钢或陶瓷制成并固连在工位地基上;
3)、WMPS测量网络定向完成后,可通过状态监控接收器采集发射站空间角度初始数据测量值θ10或θ20并储存在上位机内,工作时每隔一定监控周期接收器多次连续采集发射站的扫描角度θ1、θ2并与上位机存储的初始数据对比、判定发射站工作状态;
4)、具体发射站工作状态采用以下测量结果统计值考核:
(1)单个监控周期内的θ1或θ2的多次测量结果标准差:
该结果反映发射站轴系运行状态,当标准差大于某一阈值时,可被判断为发射站转速失稳,轴系结构存在问题;
(2)多个监控周期内的θ1或θ2的多次测量结果均值与初始数据测量值θ10或θ20之差:
该结果反映发射站与接收器间空间位置关系,当δm变化超过某一预定阈值即可判断为发射站与监控接收器间的相对位姿发生变化,工作现场中当多个接收器测量同一发射站得到的统计量δ1、δ2同时发生变化时,可判定为发射站发生旋转或位移,需对发射站重新定向;另一方面如果现场只有一个监控接收器测量结果超阈值,其他结果保持在阈值内时,可判定为该接收器坐标变化,需要重新固定接收器,并更新空间角度初始数据测量值θ10、θ20。
2.如权利要求1所述方法,其特征是,所述采集发射站的扫描角度θ1、θ2是:
1)、激光发射站的数学模型被抽象为绕公共转轴上一点的两个旋转的非平行半平面及一个以固定频率发射脉冲光的点光源,而接收器简化为以靶面光学中心为中心的质点,组装完毕后定义发射站坐标系可定义如下,以激光平面(9)与转轴的交点为激光发射站原点,旋转轴为Z轴,X轴为初始时刻,即转台转至固定位置,发射器发射脉冲光时,激光器(4)光轴所在位置,Y轴遵循右手定则,出厂前,对发射站旋转角速度ω进行设定,并标定发射站激光平面在初始位置处的平面ABCD参数,则发射站坐标系下光平面方程可表示为:
a'1、a'2、b'1、b'2、c'1、c'2、d'1、d'2即为ABCD参数,当两平面激光分别扫过接收器时其在发射站坐标系下的两平面方程系数变为转台旋转角度θ1、θ2的函数:
在经过标定建立测量坐标系后,每个发射站在测量坐标系下位姿已知,即确定了每台激光发射站原点在全局坐标系下的坐标PTXG及激光发射站坐标系到测量坐标系的旋转矩阵RTXG,发射站的两个扫描光分别到达接收器时,发射站转台转过角度分别为θ1、θ2时易知发射站两个旋转光平面的平面方程在测量坐标系下的方程参数(aTXG1(θ1),bTXG1(θ1),cTXG1(θ1),dTXG1(θ1))T及(aTXG2(θ2),bTXG2(θ2),cTXG2(θ2),dTXG2(θ2))T可表示为:
2)、当发射站光平面扫过接收器时,由于接收器在平面内,方程(1)在测量坐标系下仍然成立,如此时接收器在测量坐标系下坐标为(xRXG yRXG zRXG),可得:
观察(5)式可知,当发射站制造完成后其光平面参数(a′m b′m c′m d′m)在工作现场一般不会产生变化,此时,如果接收器坐标(xRXG yRXG zRXG)被视为定值,其测量到的转台旋转角度θ1、θ2直接反映了发射站坐标PTXG旋转矩阵RTXG的变化情况,T表示转置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110270929.2A CN102425990B (zh) | 2011-09-14 | 2011-09-14 | 工作空间测量定位系统网络状态监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110270929.2A CN102425990B (zh) | 2011-09-14 | 2011-09-14 | 工作空间测量定位系统网络状态监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102425990A CN102425990A (zh) | 2012-04-25 |
CN102425990B true CN102425990B (zh) | 2014-03-19 |
Family
ID=45959996
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110270929.2A Active CN102425990B (zh) | 2011-09-14 | 2011-09-14 | 工作空间测量定位系统网络状态监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102425990B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103512499A (zh) * | 2013-10-22 | 2014-01-15 | 天津大学 | 一种基于光电扫描的单站式三维坐标测量方法 |
CN103616662A (zh) * | 2013-11-07 | 2014-03-05 | 天津大学 | 一种用于角度交会测量的可顶装倒置发射基站 |
CN105222718B (zh) * | 2015-09-21 | 2017-05-17 | 天津大学 | 室内空间测量定位网络动态坐标测量多站数据同步方法 |
CN105607034A (zh) | 2015-12-23 | 2016-05-25 | 北京凌宇智控科技有限公司 | 一种三维空间检测系统、定位方法及系统 |
CN107144273B (zh) * | 2017-06-06 | 2018-05-22 | 天津大学 | 基于倾角传感的室内测量定位系统基站姿态自动补偿方法 |
CN107388974B (zh) * | 2017-09-01 | 2019-10-15 | 浙江华东工程安全技术有限公司 | 光电式双向位移量测方法 |
CN114413754B (zh) * | 2022-01-05 | 2023-11-07 | 天津大学 | 工作空间测量定位系统动态误差建模与自主补偿方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007124902A2 (de) * | 2006-04-28 | 2007-11-08 | Prueftechnik Dieter Busch Ag | Vorrichtung und verfahren zur beurteilung der relativen raumlage zweier gegenstände |
CN100575865C (zh) * | 2008-07-18 | 2009-12-30 | 西安交通大学 | 基于双旋转激光平面发射机网络的空间定位方法 |
CN101373130A (zh) * | 2008-10-10 | 2009-02-25 | 西安交通大学 | 空间测量定位系统的相关时间提取及特征时间计算方法 |
US8803055B2 (en) * | 2009-01-09 | 2014-08-12 | Automated Precision Inc. | Volumetric error compensation system with laser tracker and active target |
CN101504275A (zh) * | 2009-03-11 | 2009-08-12 | 华中科技大学 | 一种基于空间无线定位的手持式线激光三维测量系统 |
CN101726257B (zh) * | 2009-12-22 | 2011-03-09 | 西安交通大学 | 多目大范围激光扫描测量方法 |
TR201208632T2 (tr) * | 2010-01-27 | 2012-10-22 | Engi̇n Tuncer Temel | Anten dizilerinin aynı anda pozisyon, karşılıklı bağlaşım ve kazanç/faz kalibrasyonu için bir sistem ve yöntem. |
-
2011
- 2011-09-14 CN CN201110270929.2A patent/CN102425990B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102425990A (zh) | 2012-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102425990B (zh) | 工作空间测量定位系统网络状态监测方法 | |
CN102435177B (zh) | 室内测量定位系统单发射站位姿参数在线修正方法 | |
Yi et al. | Recent research and applications of GPS‐based monitoring technology for high‐rise structures | |
CN101672913B (zh) | 激光三点动态定位方法及系统 | |
CN101962925B (zh) | 基于轨道精密控制网的轨道三维坐标高效测量方法 | |
US8077098B2 (en) | Antenna test system | |
CN102980551B (zh) | 一种基于激光超声融合扫描的无线定位平整度检测系统 | |
CN103064098B (zh) | 卫星导航数字多波束发射阵列天线指向偏差修正方法 | |
CN105425208A (zh) | 一种用于无人机精确导航的定位系统及定位方法 | |
CN110001710B (zh) | 一种磁悬浮轨道三维场景快速重构系统、方法及应用 | |
CN104101786A (zh) | 全空域有源多波束球面相控阵天线方向图测量系统 | |
CN104880204B (zh) | 利用gps及自动跟踪与测量系统对高精度激光测距仪的校准方法 | |
CN106842229A (zh) | 一种激光云台智能空间定位系统及空间定位方法 | |
CN102997846A (zh) | 基于工作空间测量定位系统的飞机全机水平测量方法 | |
CN102176003B (zh) | 一种机载激光雷达航测参数的优化设计方法 | |
CN105157687A (zh) | 一种基于wMPS的动态物体的位置姿态测量方法 | |
CN101832773B (zh) | 三维地貌观测装置 | |
CN204495330U (zh) | 重载铁路路基沉降自动监测及预警系统 | |
CN103512499A (zh) | 一种基于光电扫描的单站式三维坐标测量方法 | |
CN103278132A (zh) | 自行走式地下掘进机器人姿态测量装置及其测量方法 | |
CN102636774A (zh) | 一种基于光电扫描与超声测距的全站式空间测量定位方法 | |
CN104777452B (zh) | 一种移动设备的定位系统及定位方法 | |
CN112305570B (zh) | 一种导航定位设备的动态性能测试系统及测试方法 | |
CN105157729A (zh) | 一种用于室内空间测量定位系统的调平与校准方法 | |
CN202583889U (zh) | Dcs监控激光标校太阳能聚光反射镜定位装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |