CN108871373B - 一种基于俯仰滚转台和非线性优化的星敏感器标定方法 - Google Patents
一种基于俯仰滚转台和非线性优化的星敏感器标定方法 Download PDFInfo
- Publication number
- CN108871373B CN108871373B CN201710329219.XA CN201710329219A CN108871373B CN 108871373 B CN108871373 B CN 108871373B CN 201710329219 A CN201710329219 A CN 201710329219A CN 108871373 B CN108871373 B CN 108871373B
- Authority
- CN
- China
- Prior art keywords
- star
- star sensor
- coordinate system
- parameter
- spot
- 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
- 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 By Optical Means (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种基于俯仰滚转台和非线性优化的星敏感器标定方法,首先搭建基于俯仰滚转台和非线性优化的星敏感器标定系统,包括:参数初始值估计模块、参数记录模块、参数求取模块、参数优化模块和参数确定模块;然后依次通过估计星敏感器的焦距和主点的初始值;记录星点光斑在星敏感器测量坐标系下的质心坐标和平行光在转动后的星敏感器测量坐标系下的方向矢量;求取星点光斑对在不同坐标系下的角距余弦值;构建目标函数优化星敏感器的光学参数;通过星点光斑的重投影误差判断算法是否满足迭代终止条件确定最终的优化参数。本方法的优点在于能够对星敏感器的光学参数进行标定,同时可以估计俯仰滚转台和星敏感器之间的旋转矩阵。
Description
技术领域
本发明涉及一种星敏感器标定方法,特别是一种基于俯仰滚转台和非线性优化的星敏感器标定方法。
背景技术
星敏感器作为目前航天器中最重要的姿态测量敏感器,其精度直接影响航天器的姿态定位精度,因此对其光学参数进行标定显得尤为重要。星敏感器的光学参数主要包括其主点和焦距,通常,在星敏感器的研制过程中,会带来主点和焦距的误差,这些误差会影响星敏感器质心定位的精度,降低星图识别算法的准确率,进而影响星敏感器的姿态定位精度。
传统的通过旋转俯仰滚转台来标定星敏感器的方法,通常忽略星敏器与转台之间的旋转矩阵的影响,估计出来的星敏感器光学参数误差比较大。
发明内容
本发明目的在于提供一种基于俯仰滚转台和非线性优化的星敏感器标定方法,解决传统星敏感器光学参数误差比较大的问题。
一种基于俯仰滚转台和非线性优化的星敏感器标定方法的具体步骤为:
第一步 搭建基于俯仰滚转台和非线性优化的星敏感器标定系统
基于俯仰滚转台和非线性优化的星敏感器标定系统,包括:参数初始值估计模块、参数记录模块、参数求取模块、参数优化模块和参数确定模块。
参数初始值估计模块的功能为:估计星敏感器的焦距和主点的初始值;
参数记录模块的功能为:记录星点光斑在星敏感器测量坐标系下的质心坐标和平行光在转动后的星敏感器测量坐标系下的方向矢量;
参数求取模块的功能为:求取星点光斑对在不同坐标系下的角距余弦值;
参数优化模块的功能为:构建目标函数优化星敏感器的光学参数;
参数确定模块的功能为:通过星点光斑的重投影误差判断算法是否满足迭代终止条件确定最终的优化参数;
第二步 参数初始值估计模块估计星敏感器的焦距和主点的初始值
参数初始值估计模块根据星敏感器的图像分辨率和视场角估计星敏感器的焦距f,见公式(1)所示,主点(x0,y0)的初始值为星敏感器成像平面的中心位置,式中Fx为星载相机在x方向的视场宽度,设定转台坐标系为XtYtZt,星敏感器测量坐标系为XsYsZs,两个坐标系之间存在转动矩阵Mst,Mst初始值为3×3阶单位矩阵I;
第三步 参数记录模块记录星点光斑在星敏感器测量坐标系下的质心坐标和平行光在转动后的星敏感器测量坐标系下的方向矢量
调节俯仰滚转台采集平行光管在星敏感器的星点光斑,参数记录模块记录星点光斑在星敏感器测量坐标系下的质心坐标(xs,ys),则在星敏感器测量坐标系下平行光矢量的测量值为b=(xs ys -f)T,令平行光矢量在XtYtZt坐标系下为(0,0,-1)T,当转台先绕Yt轴转动了角度θ,又绕Xt轴转动了角度γ时,则平行光在转动后的转台坐标系下矢量为:
结合所述的矩阵Mst,则平行光在转动后的星敏感器测量坐标系下方向矢量为:
第四步 参数求取模块求取星点光斑对在不同坐标系下的角距余弦值
参数求取模块利用星敏感器得到平行光管所呈现的星点光斑在星敏感器成像坐标系的坐标值,利用公式(4)求取星点光斑对在星敏感器测量坐标系下的角距余弦值,两两不同的星点光斑组成星点光斑对,利用星点光斑在转动后的星敏感器测量坐标系的方向矢量,结合公式(5)计算星点光斑对方向矢量之间的角距余弦值;
fij=ri Trj (4)
其中
式中,i和j分别表示第i个和第j个星点光斑,bi和ri分别表示第i个星点光斑在星敏感器测量坐标系和转台坐标系下的方向矢量,(x0,y0)为星敏感器的主点,f为星敏感器焦距,fij表示第i个和第j个星点光斑在转动后的星敏感器测量坐标系下的星间角距余弦值,表示第i个星点光斑与第j个星点光斑在星敏感器测量坐标系下的星间角距余弦值。
第五步 参数优化模块通过构建目标函数优化星敏感器的光学参数
式中,表示第i个星点光斑与第j个星点光斑在星敏感器测量坐标系下的星间角距余弦值,fij表示第i个和第j个星点光斑在转动后的星敏感器测量坐标系下的星间角距余弦值,利用Levenberg-Marquardt非线性优化算法实现目标函数的优化,获得光学参数的估计值。
第六步 参数确定模块通过星点光斑的重投影误差判断算法是否满足迭代终止条件,确定最终的优化参数
参数确定模块根据估计得到的主点和焦距通过转动转台使星点光斑移动到主点位置,记录下转台转动的角度,解算出星敏感器的测量坐标系和转台坐标系的转动矩阵估计值将解算的重新代入公式(3)得到平行光在转动后的星敏感器测量坐标系下方向矢量为根据小孔成像原理,得到平行光在星敏感器测量坐标下的光斑重投影质心坐标为
式中,i表示第i个星点光斑,N表示星点光斑的总个数,(xsi,ysi)和分别表示第i个星点光斑在星敏感器测量坐标下的质心坐标和重投影质心坐标;当重投影误差εp小于容许误差ξp时,则算法结束,得到最终的星敏感器主点和焦距估计值;当重投影误差εp大于等于容许误差ξp时,则利用估计值 以及更新星敏感器的主点、焦距以及星敏感器的测量坐标系和转台坐标系之间的转动矩阵,重复执行步骤三至五,直到重投影误差小于容许误差为止。
至此,实现一种基于俯仰滚转台和非线性优化的星敏感器的标定。
本发明对传统俯仰滚转台标定星敏感器的方法的改进,主要体现在以下两个方面:一是利用星点光斑的角距在不同坐标系下的不变性,构建目标函数,利用文献1:More,J.J.:The Levenberg-Marquardt algorithm:Implementation andtheory.Number.Analysis 630,106-115提出的Levenberg-Marquardt非线性优化方法完成光学参数的优化;二是结合俯仰滚转台和Levenberg-Marquardt非线性优化方法能够完成对俯仰滚转台和星敏感器之间的旋转矩阵的估计。本发明的优点在于能够对星敏感器的光学参数进行标定,同时能够估计俯仰滚转台和星敏感器之间的旋转矩阵。
附图说明
图1一种基于俯仰滚转台和非线性优化的星敏感器标定方法的流程图;
图2一种基于俯仰滚转台和非线性优化的星敏感器标定方法的试验系统的示意图。
1.平行光管 2.星敏感器 3.转台
具体实施方式
一种基于俯仰滚转台和非线性优化的星敏感器标定方法的具体步骤为:
第一步 搭建基于俯仰滚转台和非线性优化的星敏感器标定系统
基于俯仰滚转台和非线性优化的星敏感器标定系统,包括:参数初始值估计模块、参数记录模块、参数求取模块、参数优化模块和参数确定模块。
参数初始值估计模块的功能为:估计星敏感器的焦距和主点的初始值;
参数记录模块的功能为:记录星点光斑在星敏感器测量坐标系下的质心坐标和平行光在转动后的星敏感器测量坐标系下的方向矢量;
参数求取模块的功能为:求取星点光斑对在不同坐标系下的角距余弦值;
参数优化模块的功能为:构建目标函数优化星敏感器的光学参数;
参数确定模块的功能为:通过星点光斑的重投影误差判断算法是否满足迭代终止条件确定最终的优化参数;
第二步 参数初始值估计模块估计星敏感器的焦距和主点的初始值
参数初始值估计模块根据星敏感器的图像分辨率和视场角估计星敏感器的焦距f,见公式(1)所示,主点(x0,y0)的初始值为星敏感器成像平面的中心位置,式中Fx为星载相机在x方向的视场宽度,设定转台坐标系为XtYtZt,星敏感器测量坐标系为XsYsZs,两个坐标系之间存在转动矩阵Mst,Mst初始值为3×3阶单位矩阵I;
第三步 参数记录模块记录星点光斑在星敏感器测量坐标系下的质心坐标和平行光在转动后的星敏感器测量坐标系下的方向矢量
调节俯仰滚转台采集平行光管在星敏感器的星点光斑,参数记录模块记录星点光斑在星敏感器测量坐标系下的质心坐标(xs,ys),则在星敏感器测量坐标系下平行光矢量的测量值为b=(xs ys -f)T,令平行光矢量在XtYtZt坐标系下为(0,0,-1)T,当转台先绕Yt轴转动了角度θ,又绕Xt轴转动了角度γ时,则平行光在转动后的转台坐标系下矢量为:
结合所述的矩阵Mst,则平行光在转动后的星敏感器测量坐标系下方向矢量为:
第四步 参数求取模块求取星点光斑对在不同坐标系下的角距余弦值
参数求取模块利用星敏感器得到平行光管所呈现的星点光斑在星敏感器成像坐标系的坐标值,利用公式(4)求取星点光斑对在星敏感器测量坐标系下的角距余弦值,两两不同的星点光斑组成星点光斑对,利用星点光斑在转动后的星敏感器测量坐标系的方向矢量,结合公式(5)计算星点光斑对方向矢量之间的角距余弦值;
fij=ri Trj (4)
其中
式中,i和j分别表示第i个和第j个星点光斑,bi和ri分别表示第i个星点光斑在星敏感器测量坐标系和转台坐标系下的方向矢量,(x0,y0)为星敏感器的主点,f为星敏感器焦距,fij表示第i个和第j个星点光斑在转动后的星敏感器测量坐标系下的星间角距余弦值,表示第i个星点光斑与第j个星点光斑在星敏感器测量坐标系下的星间角距余弦值。
第五步 参数优化模块通过构建目标函数优化星敏感器的光学参数
式中,表示第i个星点光斑与第j个星点光斑在星敏感器测量坐标系下的星间角距余弦值,fij表示第i个和第j个星点光斑在转动后的星敏感器测量坐标系下的星间角距余弦值,利用Levenberg-Marquardt非线性优化算法实现目标函数的优化,获得光学参数的估计值。
第六步 参数确定模块通过星点光斑的重投影误差判断算法是否满足迭代终止条件,确定最终的优化参数
参数确定模块根据估计得到的主点和焦距通过转动转台使星点光斑移动到主点位置,记录下转台转动的角度,解算出星敏感器的测量坐标系和转台坐标系的转动矩阵估计值将解算的重新代入公式(3)得到平行光在转动后的星敏感器测量坐标系下方向矢量为根据小孔成像原理,得到平行光在星敏感器测量坐标下的光斑重投影质心坐标为
式中,i表示第i个星点光斑,N表示星点光斑的总个数,(xsi,ysi)和分别表示第i个星点光斑在星敏感器测量坐标下的质心坐标和重投影质心坐标;当重投影误差εp小于容许误差ξp时,则算法结束,得到最终的星敏感器主点和焦距估计值;当重投影误差εp大于等于容许误差ξp时,则利用估计值 以及更新星敏感器的主点、焦距以及星敏感器的测量坐标系和转台坐标系之间的转动矩阵,重复执行步骤三至五,直到重投影误差小于容许误差为止。
至此,实现一种基于俯仰滚转台和非线性优化的星敏感器的标定。
Claims (2)
1.一种基于俯仰滚转台和非线性优化的星敏感器标定方法,其特征在于具体步骤为:
第一步 搭建基于俯仰滚转台和非线性优化的星敏感器标定系统
基于俯仰滚转台和非线性优化的星敏感器标定系统,包括:参数初始值估计模块、参数记录模块、参数求取模块、参数优化模块和参数确定模块;
参数初始值估计模块的功能为:估计星敏感器的焦距和主点的初始值;
参数记录模块的功能为:记录星点光斑在星敏感器测量坐标系下的质心坐标和平行光在转动后的星敏感器测量坐标系下的方向矢量;
参数求取模块的功能为:求取星点光斑对在不同坐标系下的角距余弦值;
参数优化模块的功能为:构建目标函数优化星敏感器的光学参数;
参数确定模块的功能为:通过星点光斑的重投影误差与容许误差比较,判断算法是否满足迭代终止条件确定最终的优化参数;
第二步 参数初始值估计模块估计星敏感器的焦距和主点的初始值
参数初始值估计模块根据星敏感器的图像分辨率和视场角估计星敏感器的焦距f,见公式(1)所示,主点(x0,y0)的初始值为星敏感器成像平面的中心位置,式中Fx为星载相机在x方向的视场宽度,设定转台坐标系为XtYtZt,星敏感器测量坐标系为XsYsZs,两个坐标系之间存在转动矩阵Mst;
第三步 参数记录模块记录星点光斑在星敏感器测量坐标系下的质心坐标和平行光在转动后的星敏感器测量坐标系下的方向矢量
调节俯仰滚转台采集平行光管在星敏感器的星点光斑,参数记录模块记录星点光斑在星敏感器测量坐标系下的质心坐标(xs,ys),则在星敏感器测量坐标系下平行光矢量的测量值为b=(xs,ys,-f)T,令平行光矢量在XtYtZt坐标系下为(0,0,-1)T,当转台先绕Yt轴转动了角度θ,又绕Xt轴转动了角度γ时,则平行光在转动后的转台坐标系下矢量为:
结合所述的矩阵Mst,则平行光在转动后的星敏感器测量坐标系下方向矢量为:
第四步 参数求取模块求取星点光斑对在不同坐标系下的角距余弦值
参数求取模块利用星敏感器得到平行光管所呈现的星点光斑在星敏感器成像坐标系的坐标值,利用公式(4)求取星点光斑对在星敏感器测量坐标系下的角距余弦值,两两不同的星点光斑组成星点光斑对,利用星点光斑在转动后的星敏感器测量坐标系的方向矢量,结合公式(5)计算星点光斑对方向矢量之间的角距余弦值;
fij=ri Trj (4)
其中
式中,i和j分别表示第i个和第j个星点光斑,bi和ri分别表示第i个星点光斑在星敏感器测量坐标系和转台坐标系下的方向矢量,(x0,y0)为星敏感器的主点,f为星敏感器焦距,fij表示第i个和第j个星点光斑在转动后的星敏感器测量坐标系下的星间角距余弦值,表示第i个星点光斑与第j个星点光斑在星敏感器测量坐标系下的星间角距余弦值;
第五步 参数优化模块通过构建目标函数优化星敏感器的光学参数
式中,表示第i个星点光斑与第j个星点光斑在星敏感器测量坐标系下的星间角距余弦值,fij表示第i个和第j个星点光斑在转动后的星敏感器测量坐标系下的星间角距余弦值,利用Levenberg-Marquardt非线性优化算法实现目标函数的优化,获得光学参数的估计值;
第六步 参数确定模块通过星点光斑的重投影误差与容许误差比较,判断算法是否满足迭代终止条件,确定最终的优化参数
参数确定模块根据估计得到的主点和焦距通过转动转台使星点光斑移动到主点位置,记录下转台转动的角度,解算出星敏感器的测量坐标系和转台坐标系的转动矩阵估计值将解算的重新代入公式(3)得到平行光在转动后的星敏感器测量坐标系下方向矢量为根据小孔成像原理,得到平行光在星敏感器测量坐标下的光斑重投影质心坐标为
式中,i表示第i个星点光斑,N表示星点光斑的总个数,(xsi,ysi)和分别表示第i个星点光斑在星敏感器测量坐标下的质心坐标和重投影质心坐标;当重投影误差εp小于容许误差ξp时,则算法结束,得到最终的星敏感器主点和焦距估计值;当重投影误差εp大于等于容许误差ξp时,则利用估计值以及更新星敏感器的主点、焦距以及星敏感器的测量坐标系和转台坐标系之间的转动矩阵,重复执行步骤三至五,直到重投影误差小于容许误差为止;
至此,实现一种基于俯仰滚转台和非线性优化的星敏感器的标定。
2.如权利要求1所述的基于俯仰滚转台和非线性优化的星敏感器标定方法,其特征在于所述Mst初始值为3×3阶单位矩阵I。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710329219.XA CN108871373B (zh) | 2017-05-11 | 2017-05-11 | 一种基于俯仰滚转台和非线性优化的星敏感器标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710329219.XA CN108871373B (zh) | 2017-05-11 | 2017-05-11 | 一种基于俯仰滚转台和非线性优化的星敏感器标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108871373A CN108871373A (zh) | 2018-11-23 |
CN108871373B true CN108871373B (zh) | 2021-06-15 |
Family
ID=64319422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710329219.XA Active CN108871373B (zh) | 2017-05-11 | 2017-05-11 | 一种基于俯仰滚转台和非线性优化的星敏感器标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108871373B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110702098B (zh) * | 2019-10-14 | 2022-12-09 | 中国科学院新疆理化技术研究所 | 一种基于星对角距测量精度的星敏感器辐射损伤实验室评估方法 |
CN111402300B (zh) * | 2020-04-21 | 2022-09-20 | 中国科学院光电技术研究所 | 一种基于双谱域主成分分析的高动态星敏感器运动参数估计方法 |
CN112683307B (zh) * | 2020-12-14 | 2022-08-19 | 中国科学院光电技术研究所 | 一种星敏感器星点质心位置精度测量方法 |
CN113218577A (zh) * | 2021-04-15 | 2021-08-06 | 中国科学院光电技术研究所 | 一种星敏感器星点质心位置精度的外场测量方法 |
CN113177987B (zh) * | 2021-04-30 | 2022-03-25 | 北京航空航天大学 | 一种视觉跟踪测量系统外场全局标定方法及系统 |
CN114370866B (zh) * | 2021-11-23 | 2024-01-30 | 上海航天控制技术研究所 | 一种星敏感器主点和主距测量系统及方法 |
CN116399233B (zh) * | 2023-06-07 | 2023-10-27 | 北京融为科技有限公司 | 一种星载激光载荷伺服转台安装误差标定装置和方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1949002A (zh) * | 2005-10-12 | 2007-04-18 | 北京航空航天大学 | 一种星敏感器内外方元素校准方法 |
CN101013033A (zh) * | 2006-03-21 | 2007-08-08 | 北京航空航天大学 | 一种基于无偏差带的星敏感器地面校准方法 |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 |
-
2017
- 2017-05-11 CN CN201710329219.XA patent/CN108871373B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1949002A (zh) * | 2005-10-12 | 2007-04-18 | 北京航空航天大学 | 一种星敏感器内外方元素校准方法 |
CN101013033A (zh) * | 2006-03-21 | 2007-08-08 | 北京航空航天大学 | 一种基于无偏差带的星敏感器地面校准方法 |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 |
Non-Patent Citations (2)
Title |
---|
"接近段目标双目视觉测量中的精度分析与像机标定方法研究";周致远;《中国优秀硕士学位论文全文数据库信息科技辑》;20150115;正文第46、56页 * |
"星敏感器模型参数分析与校准方法研究";郝雪涛;《光电工程》;20050331;第32卷(第3期);正文第5-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108871373A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108871373B (zh) | 一种基于俯仰滚转台和非线性优化的星敏感器标定方法 | |
CN110390697B (zh) | 一种基于lm算法的毫米波雷达与摄像头联合标定方法 | |
CN108413988B (zh) | 机器人末端经纬仪坐标系快速标定方法 | |
CN111220120B (zh) | 一种动平台双目测距自校准方法及装置 | |
CN109459058B (zh) | 一种基于三轴转台的多视场星敏感器的标定方法 | |
CN108648242B (zh) | 基于激光测距仪辅助无公共视场的两相机标定方法及装置 | |
CN112683307B (zh) | 一种星敏感器星点质心位置精度测量方法 | |
CN110398208A (zh) | 基于摄影全站仪系统的大数据变形监测方法 | |
CN109087355A (zh) | 基于迭代更新的单目相机位姿测量装置与方法 | |
CN111457911B (zh) | 一种基于偏振二维残差信息的仿生偏振罗盘标定方法 | |
CN106525054A (zh) | 一种采用星上推扫遥感图像信息的单星自主测定轨方法 | |
CN110793542A (zh) | 一种基于广义探元指向角的面阵光学遥感卫星在轨几何定标方法 | |
CN110068313B (zh) | 一种基于投影变换的数字天顶仪定向方法 | |
CN108225305B (zh) | 一种星敏感器多参数优化方法及系统 | |
CN117197241B (zh) | 一种基于多目视觉的机器人末端绝对位姿高精度跟踪方法 | |
CN103685936A (zh) | 宽视场摄像机图像的校准和扭曲恢复 | |
CN111998823B (zh) | 基于双目异光源测距装置的目标测距方法 | |
CN106023146B (zh) | 用于摄影测量中的场相关单边自标定光束平差方法 | |
CN112116665A (zh) | 一种结构光传感器标定方法 | |
CN113781581B (zh) | 基于靶标松姿态约束的景深畸变模型标定方法 | |
CN115797463A (zh) | 一种fod雷达与相机的神经网络训练方法及全局标定方法 | |
CN112485773B (zh) | 激光雷达与倾角传感器的外参信息标定方法 | |
CN110715795B (zh) | 一种光电跟踪系统中快速反射镜的标定和测量方法 | |
CN109506645B (zh) | 一种星敏感器安装矩阵地面精测方法 | |
CN112927299B (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 |