CN112414676B - 一种多自由曲面反射镜测量数据的评价计算方法 - Google Patents
一种多自由曲面反射镜测量数据的评价计算方法 Download PDFInfo
- Publication number
- CN112414676B CN112414676B CN202011225596.7A CN202011225596A CN112414676B CN 112414676 B CN112414676 B CN 112414676B CN 202011225596 A CN202011225596 A CN 202011225596A CN 112414676 B CN112414676 B CN 112414676B
- Authority
- CN
- China
- Prior art keywords
- free
- reg
- model
- coordinate
- measured data
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明涉及一种多自由曲面反射镜测量数据的评价计算方法,包括下列步骤:(1)使用多表面测量数据中所有表面的测量数据与多自由曲面反射系统的设计模型进行坐标配准,完成测量坐标系与模型坐标系的统一;(2)使用经过坐标配准后的各表面测量数据分别与各自面型的设计模型方程进行坐标配准,完成各表面面形与位置误差的评价。
Description
技术领域
本发明涉及一种反射镜测量数据的评价计算方法。
背景技术
多自由曲面反射系统指由一个以上反射镜表面共同完成其设计功能效果的系统,系统功能由这些表面的质量以及各个表面之间的相对位置关系共同决定。各表面间可能相互关联存在共同的定位基准,也可能互相独立于空间之中。多反射镜光学系统的一体制造极大程度的降低了多个表面装配组装的难度,进一步推进了自由曲面的广泛应用。高精度的完整测量与评价是多反射光学系统制造的基础保障,如今成熟的测量技术能够完成对单个表面的高精度测量,在曲面的实际评价过程中,测量曲面往往只相当于整个设计曲面上的一部分,并且测量曲面在设计曲面上的具体对应区域是未知的,需要确定测量曲面在设计曲面上的对应区域,通常要求具有较好的初值位置及几何约束,才能完成测量曲面和与之对应的设计曲面之间的匹配实现面形误差评价。而多自由曲面反射系统测量除了面形评价外还需要各表面相对位置误差的评价方法。
目前多自由曲面反射系统的评价方法往往需要添加额外的空间位置基准,或者使用表面本身的棱边、顶点、几何中心等几何特征作为相对位置定位基准,但这些做法与多自由曲面反射系统的设计理念相悖,同时实际测量过程中大多测量方法表面的测量范围有限难以完全包含测量表面的几何特征,鲜有能够完成多自由曲面反射系统相对位置测量的设备,并且大部分设备只能完成单一类型的多自由曲面反射系统的测量评价,缺少通用性。当下仍然缺少一种多自由曲面反射系统测量数据的位置误差与面形误差评价方法,能够不依赖特定测量仪器或测量系统,较为准确的完成多自由曲面反射系统的测量评价。
发明内容
针对上述问题,本发明的目的是提供一种多表面面形测量数据的评价算法,无需在被测系统内添加额外的空间位置参考基准,无需计算表面的棱边、顶点、几何中心等几何特征,完成各表面面形误差与相对位置的评价。为实现上述目的,本发明采取以下的技术方案。
一种多自由曲面反射镜测量数据的评价计算方法,包括下列步骤:
(1)使用多表面测量数据中所有表面的测量数据与多自由曲面反射系统的设计模型进行坐标配准,完成测量坐标系与模型坐标系的统一:用于坐标配准过程的刚体变换齐次矩阵H通过绕过点(xL,yL,zL),单位矢量方向为[ex,ey,ez]T的空间直线L与沿空间笛卡尔坐标系的X轴,Y轴和Z轴的平移量a,b和c表示,表达方式为:
其中,θ是绕单位矢量[ex,ey,ez]T旋转的角度;
设测量中,第j个表面的第i个测量数据点的坐标为P(i,j),则其空间坐标如下表示:
P(i,j)=(X(i,j),Y(i,j),Z(i,j))
设P(i,j)经过刚体变换后对应的模型点为Q(i,j),此对应的模型点根据点距最近原则搜索得到;坐标配准过程优化求解刚体变换齐次矩阵H,使得测量数据到模型对应点的几何距离最小化,优化目标函数如下所示:
(2)使用经过坐标配准后的各表面测量数据分别与各自面型的设计模型方程进行坐标配准,完成各表面面形与位置误差的评价。方法如下:
设第j个表面经过第一步坐标配准后的测量数据点为T(i,j),则该面测量数据与其理想模型坐标配准的变换关系按照以下方式进计算:
其中,Rj和tj分别为第j个表面测量数据坐标配准到其表面设计模型上的刚体变换旋转矩阵和平移矢量;D[·]为距离算子,表示测量数据点T(i,j)经过刚体变换旋转矩阵和平移矢量的计算后与其理想表面的距离;此处Rj和tj的表达式为:
tj=[Δxj,Δyj,Δzj]T
坐标配准后的测量数据记为Treg(i,j),该数据坐标通过以下变换得到:
Treg(i,j)=RjT(i,j)+tj=(Xreg(i,j),Yreg(i,j),Zreg(i,j))T
其中Xreg,Yreg和Zreg是坐标配准后的测量数据的笛卡尔坐标;
此时,将各表面测量数据坐标配准后各数据点到自由曲面反射镜模型的距离作为面形误差,设第j个自由曲面反射镜的面形方程为z=Fj(x,y),则各测量数据点对应的面形误差E(i,j)为:
E(i,j)=Zreg(i,j)-Fj(Xreg(i,j),Yreg(i,j))
第j个自由曲面反射镜位置误差定义为其与模型面型间的位移偏差量Tj和绕模型坐标系中某一矢量方向Vj的偏差角度θj;其中Tj与平移矩阵tj一致,即Δxj,Δyj和Δzj分别为第j个自由曲面反射镜沿着沿空间笛卡尔坐标系的X轴,Y轴和Z轴方向的位移偏差量;Vj与θj通过对旋转矩阵Rj使用下述公式计算得到:
附图说明
图1为本发明所涉及的多自由曲面反射系统测量数据误差评价方法流程图。
图2为本发明实施例中离轴二反自由曲面系统的面形与位置误差分析过程。其中,(a)是两个反射面表面测量数据与两反设计模型数据初步配准结果;(b)是精细坐标配准结果;(c)是自由曲面反射面一面测量数据与其设计模型的匹配结果;(d)自由曲面反射面二面测量数据与其设计模型的匹配结果。
图3本发明实施例中离轴二反自由曲面系统两个表面测量数据的面形误差。其中(a)是第一个自由曲面反射面面形误差;(b)是第二个自由曲面反射面面形误差。
具体实施方式
结合附图和实施例对本发明的一种多表面面形测量数据的评价计算方法加以说明。评价计算方法流程如图1所示。首先,使用多表面测量数据中所有表面的测量数据与多自由曲面反射系统的设计模型进行坐标配准,完成测量坐标系与模型坐标系的统一,然后,将坐标系统一后的各表面测量数据与各表面的设计模型进行单独的坐标配准,最后,根据坐标配准结果获得被测表面面形与位置误差评价结果。
(1)将多表面测量数据中所有表面的测量数据与多自由曲面反射系统的模型进行匹配,完成测量坐标系与模型坐标系的统一。其中,用于坐标配准过程的刚体变换齐次矩阵H通过绕过点(xL,yL,zL),单位矢量方向为[ex,ey,ez]T的空间直线L与沿空间笛卡尔坐标系的X轴,Y轴和Z轴的平移量a,b和c表示,表达方式为:
其中:
具体的,对于刚体变换齐次矩阵H的计算包括以下步骤:
第一步,进行初步坐标配准计算。设第j个反射镜的第i个测量数据点的坐标为P(i,j),则其空间坐标表示为:
P(i,j)=(X(i,j),Y(i,j),Z(i,j))
多自由曲面反射镜的第j个反射镜的第i个模型数据点的坐标为Q(i,j),并且该处的法矢为nq(i,j),根据自由曲面反射镜模型的方程计算得到。经过配准变换后的P(i,j)与Q(i,j)根据点距最近原则确定为匹配点对。其优化函数的目标为:
其中Hrough为测量数据转化的刚体变换齐次矩阵。此时,对于测量数据点P(i,j),经过平移旋转后所获得的第j个反射镜的第i个数据点的坐标Pr(i,j)的表达式为
Pr(i,j)=HroughP(i,j)
第二步,基于初步坐标配准计算结果进行精细坐标配准计算。根据多自由曲面反射镜模型的基准方程,计算测量数据点沿法矢方向到自由曲面反射镜的空间距离,设D[·]为距离计算算子,表示测量数据点与其理想表面的距离。则通过以下目标函数优化计算两坐标系间的刚体变换齐次矩阵,进一步提升多表面测量数据与设计模型间的匹配精度。
其中Hfine为精细坐标配准过程中测量坐标系与模型坐标系统一的刚体变换齐次矩阵。经过测量坐标系与模型坐标系统一后,所获得的第j个反射镜的第i个数据点的坐标T(i,j)是与模型坐标系统一的测量数据,其表达式为:
T(i,j)=HfinePr(i,j)
(2)将每个表面的测量数据分别与其表面的设计模型进行匹配,设第j个表面经过测量坐标系与模型坐标系统一后的测量数据点为T(i,j),则该测量数据点与其理想模型匹配的变换关系按照以下方式进计算:
其中,Rj和tj分别为第j个表面测量数据坐标配准到其表面设计模型上的刚体变换旋转矩阵和平移矢量。此处Rj和tj的计算方式为:
tj=[Δxj,Δyj,Δzj]T
坐标配准后的测量数据Treg(i,m)坐标通过变换得到:
Treg(i,j)=RjT(i,j)+tj=(Xreg(i,j),Yreg(i,j),Zreg(i,j))T
此时,将各表面测量数据坐标配准后各数据点到自由曲面反射镜模型的距离作为面形误差,设第j个自由曲面反射镜的面形方程为z=Fj(x,y),则各测量数据点对应的面形误差E(i,j)为:
E(i,j)=Zreg(i,j)-Fj(Xreg(i,j),Yreg(i,j))
第j个自由曲面反射镜位置误差定义为其与模型面型间的位移偏差量Tj和绕模型坐标系中某一矢量方向Vj的偏差角度θj。其中Tj与平移矩阵tj一致,即Δxj,Δyj和Δzj分别为第j个自由曲面反射镜沿着沿空间笛卡尔坐标系的X轴,Y轴和Z轴方向的位移偏差量。Vj与θj通过对旋转矩阵Rj使用下述公式计算得到:
根据上述的方式可以获得表面的面形误差评价结果为E(i,j);表面的位置误差的评价结果为:第j个自由曲面反射镜沿着笛卡尔坐标系的各个轴的位置平移误差为Δxj,Δyj,Δzj,旋转误差为绕着矢量方向Vj的轴旋转,旋转角度为θj。
根据上述方法和流程,可实现多自由曲面反射系统各表面面形与位置误差的准确评估。
表1两自由曲面反射面的最大误差评价结果数据
本发明方法可适用于各类光学多自由曲面反射系统的综合误差评估,下面以一离轴二反自由曲面系统为例,其系统结构如图2所示,以进一步具体描述本发明的评价方法。其系统中的两个自由曲面反射面均为自由曲面,系统材质为铝,通过相关测量设备获得两个自由曲面反射面测量数据。使用本发明的评价流程,对离轴系统进行统误差分析过程如图2所示。首先将两反射面表面测量数据与两反设计模型数据进行初始配准其结果如图2(a)所示,对测量数据进行刚体变换后,进行精细坐标配准操作,其结果如图2(b)所示,完成测量坐标系与模型坐标系的统一。随后将两表面的测量数据与各表面的设计模型进行配准,其结果如图2(c)、(d)所示,分析两表面面形与位置误差。面形误差如图3所示,各表面最大位置误差的评价结果如表1所示。
Claims (1)
1.一种多自由曲面反射镜测量数据的评价计算方法,包括下列步骤:
(1)使用多表面测量数据中所有表面的测量数据与多自由曲面反射系统的设计模型进行坐标配准,完成测量坐标系与模型坐标系的统一:用于坐标配准过程的刚体变换齐次矩阵H通过绕过点(xL,yL,zL),单位矢量方向为[ex,ey,ez]T的空间直线L与沿空间笛卡尔坐标系的X轴,Y轴和Z轴的平移量a,b和c表示,表达方式为:
其中,θ是绕单位矢量[ex,ey,ez]T旋转的角度;
设测量中,第j个表面的第i个测量数据点的坐标为P(i,j),则其空间坐标如下表示:
P(i,j)=(X(i,j),Y(i,j),Z(i,j))
设P(i,j)经过刚体变换后对应的模型点为Q(i,j),此对应的模型点根据点距最近原则搜索得到;坐标配准过程优化求解刚体变换齐次矩阵H,使得测量数据到模型对应点的几何距离最小化,优化目标函数如下所示:
(2)使用经过坐标配准后的各表面测量数据分别与各自面型的设计模型方程进行坐标配准,求得各测量数据点对应的面形误差,以及自由曲面反射镜位置误差,方法如下:
设第j个表面经过第一步坐标配准后的测量数据点为T(i,j),则该面测量数据与其理想模型坐标配准的变换关系按照以下方式进计算:
其中,Rj和tj分别为第j个表面测量数据坐标配准到其表面设计模型上的刚体变换旋转矩阵和平移矢量;D[·]为距离算子,表示测量数据点T(i,j)经过刚体变换旋转矩阵和平移矢量的计算后与其理想表面的距离;此处Rj和tj的表达式为:
tj=[Δxj,Δyj,Δzj]T
坐标配准后的测量数据记为Treg(i,j),该数据坐标通过以下变换得到:
Treg(i,j)=RjT(i,j)+tj=(Xreg(i,j),Yreg(i,j),Zreg(i,j))T
其中Xreg,Yreg和Zreg是坐标配准后的测量数据的笛卡尔坐标;
此时,将各表面测量数据坐标配准后各数据点到自由曲面反射镜模型的距离作为面形误差,设第j个自由曲面反射镜的面形方程为z=Fj(x,y),则各测量数据点对应的面形误差E(i,j)为:
E(i,j)=Zreg(i,j)-Fj(Xreg(i,j),Yreg(i,j))
第j个自由曲面反射镜位置误差定义为其与模型面型间的位移偏差量Tj和绕模型坐标系中某一矢量方向Vj的偏差角度θj;其中Tj与平移矩阵tj一致,即Δxj,Δyj和Δzj分别为第j个自由曲面反射镜沿着沿空间笛卡尔坐标系的X轴,Y轴和Z轴方向的位移偏差量;Vj与θj通过对旋转矩阵Rj使用下述公式计算得到:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011225596.7A CN112414676B (zh) | 2020-11-05 | 2020-11-05 | 一种多自由曲面反射镜测量数据的评价计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011225596.7A CN112414676B (zh) | 2020-11-05 | 2020-11-05 | 一种多自由曲面反射镜测量数据的评价计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112414676A CN112414676A (zh) | 2021-02-26 |
CN112414676B true CN112414676B (zh) | 2022-02-11 |
Family
ID=74827040
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011225596.7A Active CN112414676B (zh) | 2020-11-05 | 2020-11-05 | 一种多自由曲面反射镜测量数据的评价计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112414676B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113835221B (zh) * | 2021-10-11 | 2022-07-08 | 东北大学 | 一种多反离轴光学系统初始结构一体化建模方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001216344A (ja) * | 2000-02-03 | 2001-08-10 | Canon Inc | 非球面形状モデルの形状誤差評価装置、非球面形状モデルの形状誤差評価方法および記憶媒体 |
JP2006343234A (ja) * | 2005-06-09 | 2006-12-21 | Yaskawa Information Systems Co Ltd | 非球面形状誤差の評価方法および形状評価装置 |
DE202015102459U1 (de) * | 2014-05-13 | 2015-06-01 | Carl Mahr Holding Gmbh | Messeinrichtung zur Durchführung eines Kalibrierverfahrens |
CN104808581A (zh) * | 2015-04-20 | 2015-07-29 | 天津大学 | 一种复杂面型曲面制造的补偿加工方法 |
CN105572837A (zh) * | 2015-12-23 | 2016-05-11 | 中国科学院长春光学精密机械与物理研究所 | 一种高精度空间相机反射镜镜面面形的处理方法 |
CN107091627A (zh) * | 2017-04-18 | 2017-08-25 | 复旦大学 | 一种多表面系统的综合测量与评估方法 |
CN107941471A (zh) * | 2017-11-17 | 2018-04-20 | 全普光电科技(上海)有限公司 | 一种自由曲面的检测方法 |
CN109387938A (zh) * | 2017-08-14 | 2019-02-26 | 清华大学 | 光学系统自由曲面面形公差分布的分析方法 |
CN110245367A (zh) * | 2018-03-09 | 2019-09-17 | 中国科学院长春光学精密机械与物理研究所 | 光机集成分析面形误差处理方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
HK1095987A2 (en) * | 2006-02-03 | 2007-05-18 | Univ Hong Kong Polytechnic | Methods and apparatus of aligning surfaces |
-
2020
- 2020-11-05 CN CN202011225596.7A patent/CN112414676B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001216344A (ja) * | 2000-02-03 | 2001-08-10 | Canon Inc | 非球面形状モデルの形状誤差評価装置、非球面形状モデルの形状誤差評価方法および記憶媒体 |
JP2006343234A (ja) * | 2005-06-09 | 2006-12-21 | Yaskawa Information Systems Co Ltd | 非球面形状誤差の評価方法および形状評価装置 |
DE202015102459U1 (de) * | 2014-05-13 | 2015-06-01 | Carl Mahr Holding Gmbh | Messeinrichtung zur Durchführung eines Kalibrierverfahrens |
CN104808581A (zh) * | 2015-04-20 | 2015-07-29 | 天津大学 | 一种复杂面型曲面制造的补偿加工方法 |
CN105572837A (zh) * | 2015-12-23 | 2016-05-11 | 中国科学院长春光学精密机械与物理研究所 | 一种高精度空间相机反射镜镜面面形的处理方法 |
CN107091627A (zh) * | 2017-04-18 | 2017-08-25 | 复旦大学 | 一种多表面系统的综合测量与评估方法 |
CN109387938A (zh) * | 2017-08-14 | 2019-02-26 | 清华大学 | 光学系统自由曲面面形公差分布的分析方法 |
CN107941471A (zh) * | 2017-11-17 | 2018-04-20 | 全普光电科技(上海)有限公司 | 一种自由曲面的检测方法 |
CN110245367A (zh) * | 2018-03-09 | 2019-09-17 | 中国科学院长春光学精密机械与物理研究所 | 光机集成分析面形误差处理方法 |
Non-Patent Citations (4)
Title |
---|
A profile error evaluation method for freeform surface measured by sweep scanning on CMM;Gaiyun He 等;《Precision Engineering》;20181215;第56卷;全文 * |
Radial scan form measurement for freeform surfaces with a large curvature using stylus profilometry;Yingmo Wang 等;《Measurement Science and Technology》;20190314;第30卷(第4期);全文 * |
基于光学性能评价自由曲面离轴三反系统加工误差分析;刘现磊 等;《光学学报》;20170510;全文 * |
自由曲面精密铣削加工与面形误差评价的研究;占策;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》;20180115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112414676A (zh) | 2021-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110440692B (zh) | 激光跟踪仪与结构光3d扫描仪组合式测量标定方法 | |
Li et al. | Geometric algorithms for workpiece localization | |
US8452564B2 (en) | Method of determining geometric errors in a machine tool or measuring machine | |
CN106441117B (zh) | 基于多站etalon激光跟踪系统的转台误差检测方法 | |
CN109883443B (zh) | 一种线结构光传感器空间姿态标定方法 | |
CN112833786B (zh) | 一种舱段位姿测量及对准系统、控制方法及应用 | |
CN110202318B (zh) | 一种基于双侧超声滚压加工的航空叶片定位与姿态调节方法 | |
Menq et al. | Curve and surface approximation from CMM measurement data | |
CN110161965B (zh) | 一种大型航天机匣斜孔的在机测量方法 | |
CN107091627B (zh) | 一种多表面系统的综合测量与评估方法 | |
CN112414676B (zh) | 一种多自由曲面反射镜测量数据的评价计算方法 | |
CN112229321B (zh) | 基于lasso算法求解三坐标测量机21项几何误差的方法 | |
CN101750045A (zh) | 烟机圆柱凸轮轮廓曲线的等分度测量方法 | |
CN109323665B (zh) | 一种线结构光驱动全息干涉的精密三维测量方法 | |
Gou et al. | On the symmetric location problem | |
Peng et al. | Measurement and identification of geometric errors of translational axis based on sensitivity analysis for ultra-precision machine tools | |
Zhang et al. | A new contouring error estimation for the high form accuracy of a multi-axis CNC machine tool | |
Gąska et al. | Simulation model for uncertainty estimation of measurements performed on five-axis measuring systems | |
CN111046584A (zh) | 星载仪器设备的精密调整方法 | |
CN114485392A (zh) | 基于激光跟踪仪的航天器大尺寸机械基准建立方法和系统 | |
CN113916128A (zh) | 一种基于光笔式视觉测量系统的提高精度的方法 | |
Deng et al. | Hand-eye calibration of line structured-light sensor by scanning and reconstruction of a free-placed standard cylindrical target | |
Huang et al. | Adaptive positioning technology of film cooling holes in hollow turbine blades | |
CN114034244B (zh) | 一种基于Makima算法获得CMM测量范围内任意点几何误差的方法 | |
CN115415853A (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 |