CN106553770A - 遥感卫星姿态运动补偿的成像试验方法 - Google Patents
遥感卫星姿态运动补偿的成像试验方法 Download PDFInfo
- Publication number
- CN106553770A CN106553770A CN201611083773.6A CN201611083773A CN106553770A CN 106553770 A CN106553770 A CN 106553770A CN 201611083773 A CN201611083773 A CN 201611083773A CN 106553770 A CN106553770 A CN 106553770A
- Authority
- CN
- China
- Prior art keywords
- satellite
- attitude
- image
- attitude motion
- remote sensing
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 29
- 238000010998 test method Methods 0.000 title claims abstract description 16
- 230000000694 effects Effects 0.000 claims abstract description 10
- 238000010835 comparative analysis Methods 0.000 claims abstract description 4
- 230000003287 optical effect Effects 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000000007 visual effect Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 abstract description 2
- 238000011161 development Methods 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 13
- 238000012360 testing method Methods 0.000 description 11
- 238000005259 measurement Methods 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种遥感卫星姿态运动补偿的成像试验方法，包括以下步骤：步骤一：获得卫星零姿态下的基准图像；步骤二：获得卫星存在姿态运动时的图像；步骤三：获得卫星存在姿态运动时对姿态运动进行补偿后的图像；步骤四：对比评估姿态运动补偿效果。本发明能够在卫星研制阶段，利用地面成像系统对卫星的卫星姿态运动补偿设计方案进行验证。
Description
技术领域
本发明涉及光学成像领域和精密测量领域，更具体的说，涉及一种遥感卫星姿态运动补偿的成像试验方法，能够模拟卫星在轨姿态偏差的状态，利用卫星遥感仪器成像功能在地面评估姿态运动补偿方案。
背景技术
随着遥感卫星应用技术水平的提高，通常要求高精度遥感卫星对地成像的图像被精确定位，并配准到地球标称网格。但卫星平台和遥感仪器受到空间力学环境、热环境变化等因素的影响，卫星会存在轨道漂移和姿态指向偏差，仪器会产生几何变形，影响遥感图像的定位与配准精度。卫星图像导航与配准技术可解决此类问题：卫星对遥感仪器的二维扫描机构进行实时角度增量补偿，将仪器观测点在地表的移动轨迹引导至预设路径，从而实现图像导航，使遥感图像配准到地球标称网格。
对遥感卫星在轨期间的姿态误差补偿测试通常采用在轨试验调试的方法，该方法风险大，代价高，参数调试、代码调试流程复杂。而在地面等比例对户外真实场景成像试验受场地大小限制，不可能完全模拟在轨工况。且室外场景实际尺寸不可测，受安全性防污染等因素考虑卫星也不便于长时间在外场测试。平行光管可以大幅缩短所需光程长度，使遥感仪器成像试验在室内进行等效模拟成为可能。而从常规思路考虑验证卫星导航与配准的姿态运动补偿需将整星实际姿态偏置一个固定小角度，且需对该角度做角秒级的精确测量，这将对整星的总装、精测带来极大挑战，且需要配备高精度的角度测量仪器。
袁修孝、余翔在测绘学报第3期第41卷(2012年6月)的文章《高分辨率卫星遥感影像姿态角系统误差检校》详细推导了卫星姿态误差对像质的影响，并给出了地面标校方法，但该方法仅能通过地面修正后处理，无法实现在轨实时修正，且地面验证手段也仅仅停留在仿真分析阶段。
专利公开号为CN102267575A、名称为“敏捷卫星成像沿任意方向推扫速度失配时的姿态补偿方法”的中国专利介绍了推扫式遥感仪器姿态补偿方法，该方法通过解算下传的地面图像配准精度来评估补偿效果，仍然不能做到地面的等效试验。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种遥感卫星姿态运动补偿的成像试验方法，其保持星体实际姿态不变，仅改变靶标位置和角度，等效在轨卫星的姿态变化导致的图像偏移，在实际操作中仅需对靶标的位移进行精测，通过光管等效焦距的放大，对位移的精测精度降低，且位移的精测难度比角度测量要小很多，这对试验操作流程带来极大便利。
根据本发明的一个方面，提供一种遥感卫星姿态运动补偿的成像试验方法，其特征在于，包括以下步骤：步骤一：获得卫星零姿态下的基准图像；步骤二：获得卫星存在姿态运动时的图像；步骤三：获得卫星存在姿态运动时对姿态运动进行补偿后的图像；步骤四：对比评估姿态运动补偿效果。
2、根据权利要求1所述的遥感卫星姿态运动补偿的成像试验方法，其特征在于，所述步骤一包括以下步骤：调整平行光管光轴对准卫星仪器视场中心，卫星在不开启姿态补偿的情况下对靶标成像，获得的图像作为基准图像。
3、根据权利要求1所述的遥感卫星姿态运动补偿的成像试验方法，其特征在于，所述步骤二包括以下步骤：靶标沿卫星本体系x轴或y轴方向平移一段距离，模拟卫星在轨姿态偏移造成的目标在视场中的移动效果；卫星在不开启姿态补偿的情况下对靶标成像，获得姿态运动时的图像。
4、根据权利要求1所述的遥感卫星姿态运动补偿的成像试验方法，其特征在于，所述步骤三包括以下步骤：在卫星内部姿轨控计算机软件中设置与靶标偏移量相对应的姿态偏置参数，卫星在开启姿态运动补偿的情况下对靶标成像，获得补偿后图像。
5、根据权利要求1所述的遥感卫星姿态运动补偿的成像试验方法，其特征在于，所述步骤四包括以下步骤：将基准图像的灰度矩阵分别与姿态运动时的图像、补偿后图像的灰度矩阵相减，从而观测图像的重合程度，两图像不重合的像元个数。
与现有技术相比，本发明具有如下的有益效果：本发明保持星体实际姿态不变，仅改变靶标位置和角度，等效在轨卫星的姿态变化导致的图像偏移，在实际操作中仅需对靶标的位移进行精测，通过光管等效焦距的放大，对位移的精测精度降低，且位移的精测难度比角度测量要小很多，这对试验操作流程带来极大便利。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为卫星姿态运动补偿成像试验原理示意图。
附图2为卫星地面成像试验示意图。
附图3为卫星姿态运动补偿成像试验流程示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
本发明根据卫星姿态运动后对地面成像造成的偏移，利用靶标等效平移模拟星体姿态运动对图像造成影响。
具体的，根据本发明提供的遥感卫星姿态运动补偿的成像试验方法，包括如下步骤：
步骤一：调整平行光管光轴对准卫星仪器视场中心，卫星在不开启姿态补偿的情况下对靶标成像，获得的图像作为基准图像。
步骤二：靶标沿卫星本体系x轴或y轴方向平移一段距离，模拟卫星在轨姿态偏移造成的目标在视场中的移动效果。卫星在不开启姿态补偿的情况下对靶标成像，获得姿态运动时的图像。
步骤三：在卫星内部姿轨控计算机软件中设置与靶标偏移量相对应的姿态偏置参数，卫星在开启姿态运动补偿的情况下对靶标成像，获得补偿后图像。
步骤四：将基准图像的灰度矩阵分别与姿态运动时的图像、补偿后图像的灰度矩阵相减，从而观测图像的重合程度，两图像不重合的像元个数即为图像偏差，从而对比评估姿态运动补偿效果。
所述步骤一包括以下步骤：
建立卫星、遥感仪器、平行光管、靶标的相对关系，其中平行光管光轴与卫星本体系z轴平行，并将遥感仪器视轴、平行光管同轴对准。卫星本体系y轴垂直指向地面，靶标固定，将靶标的该位置计为零位，遥感仪器不开启姿态运动补偿，获得的图像作为基准图像。
所述步骤二包括以下步骤：靶标沿卫星本体系x轴或y轴方向平移，并测量平移距离L，模拟卫星姿态偏差造成的成像效果。若靶标沿x轴或y轴正方向平移，则位移L为正，反之为负。遥感仪器在不开启姿态补偿的情况下靶标成像，获得偏移图像。
所述步骤三包括以下步骤：星体实际位置保持不动，在卫星内部姿轨控软件中设置姿态偏置参数。若靶标沿卫星本体系x轴方向偏移，则设置卫星俯仰角，角度设置满足如下关系式：
θ＝-arctan(L/f)
若靶标沿卫星本体系y轴方向偏移，则设置卫星滚动角，角度设置满足如下关系式：
其中L为靶标平移位移，f为平行光管等效焦距。开启姿态补偿后遥感仪器对偏移的靶标成像，获得的图像即为补偿图像。
所述步骤四包括以下步骤：根据获得的基准图像和模拟姿态偏移图像，分析靶标形状特征，提取图像中靶标边缘直线，并将两图像叠加，计算直线的偏移像素个数n，则补偿前的配准精度为n1τ，补偿后的配准精度为n2τ，其中τ为遥感仪器单个像元的视场大小，比较姿态补偿前后配准精度即可对补偿效果进行评估。
由于卫星在轨期间受到内外力矩的干扰，姿态控制系统不能将卫星平台姿态严格保持在参考坐标系，导致载荷视轴指向随平台姿态运动一起发生变化，姿态运动补偿利用仪器二维扫描镜增加补偿量的方法，将姿态运动后的载荷视轴实时指向预期方向，补偿后的效果相当于卫星一直保持在0姿态状态对地球成像(图1)。本发明能够在卫星研制阶段，利用地面成像系统对卫星的卫星姿态运动补偿设计方案进行验证。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。
Claims (5)
1.一种遥感卫星姿态运动补偿的成像试验方法，其特征在于，包括以下步骤：步骤一：获得卫星零姿态下的基准图像；步骤二：获得卫星存在姿态运动时的图像；步骤三：获得卫星存在姿态运动时对姿态运动进行补偿后的图像；步骤四：对比评估姿态运动补偿效果。
2.根据权利要求1所述的遥感卫星姿态运动补偿的成像试验方法，其特征在于，所述步骤一包括以下步骤：调整平行光管光轴对准卫星仪器视场中心，卫星在不开启姿态补偿的情况下对靶标成像，获得的图像作为基准图像。
3.根据权利要求1所述的遥感卫星姿态运动补偿的成像试验方法，其特征在于，所述步骤二包括以下步骤：靶标沿卫星本体系x轴或y轴方向平移一段距离，模拟卫星在轨姿态偏移造成的目标在视场中的移动效果；卫星在不开启姿态补偿的情况下对靶标成像，获得姿态运动时的图像。
4.根据权利要求1所述的遥感卫星姿态运动补偿的成像试验方法，其特征在于，所述步骤三包括以下步骤：在卫星内部姿轨控计算机软件中设置与靶标偏移量相对应的姿态偏置参数，卫星在开启姿态运动补偿的情况下对靶标成像，获得补偿后图像。
5.根据权利要求1所述的遥感卫星姿态运动补偿的成像试验方法，其特征在于，所述步骤四包括以下步骤：将基准图像的灰度矩阵分别与姿态运动时的图像、补偿后图像的灰度矩阵相减，从而观测图像的重合程度，两图像不重合的像元个数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611083773.6A CN106553770B (zh) | 2016-11-30 | 2016-11-30 | 遥感卫星姿态运动补偿的成像试验方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611083773.6A CN106553770B (zh) | 2016-11-30 | 2016-11-30 | 遥感卫星姿态运动补偿的成像试验方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106553770A true CN106553770A (zh) | 2017-04-05 |
| CN106553770B CN106553770B (zh) | 2019-01-08 |
Family
ID=58445512
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611083773.6A Active CN106553770B (zh) | 2016-11-30 | 2016-11-30 | 遥感卫星姿态运动补偿的成像试验方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106553770B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107292831A (zh) * | 2017-04-24 | 2017-10-24 | 上海卫星工程研究所 | 静止卫星成像导航与配准的姿态运动补偿方法 |
| CN107329191A (zh) * | 2017-05-17 | 2017-11-07 | 上海卫星工程研究所 | 静止气象卫星图像导航与配准的成像试验系统及试验方法 |
| CN109141368A (zh) * | 2018-07-19 | 2019-01-04 | 北京空间机电研究所 | 一种高定位精度星载成像系统及方法 |
| CN113063434A (zh) * | 2021-02-25 | 2021-07-02 | 上海卫星工程研究所 | 卫星指向恒星的精度评估方法及系统 |
| CN113063437A (zh) * | 2021-02-26 | 2021-07-02 | 上海卫星工程研究所 | 适用于全物理仿真的卫星视线指向模拟方法及系统 |
| CN113063435A (zh) * | 2021-02-25 | 2021-07-02 | 上海卫星工程研究所 | 卫星姿态稳定度和指向精度评估方法及系统 |
| RU2798008C2 (ru) * | 2019-02-11 | 2023-06-14 | Виасат, Инк. | Дельта-кодирование для дистанционного зондирования |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1094002A2 (de) * | 1999-10-18 | 2001-04-25 | Astrium GmbH | Regelungsanordnung und Regelungsverfahren für Satelliten |
| CN1789082A (zh) * | 2005-12-16 | 2006-06-21 | 西安电子科技大学 | 在轨地球同步自旋卫星红外弦宽差分定姿方法 |
| CN103162711A (zh) * | 2013-02-28 | 2013-06-19 | 北京航空航天大学 | 一种高动态星敏感器像增强器的误差补偿方法和装置 |
| US20150298827A1 (en) * | 2014-04-22 | 2015-10-22 | Massachusetts Institute Of Technology | Attitude determination using infrared earth horizon sensors |
| CN105572692A (zh) * | 2015-12-16 | 2016-05-11 | 上海卫星工程研究所 | 卫星图像导航与配准全物理试验装置及测试方法 |
| CN106027904A (zh) * | 2016-06-29 | 2016-10-12 | 中国人民解放军国防科学技术大学 | 视频卫星姿态运动图像补偿方法 |
-
2016
- 2016-11-30 CN CN201611083773.6A patent/CN106553770B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1094002A2 (de) * | 1999-10-18 | 2001-04-25 | Astrium GmbH | Regelungsanordnung und Regelungsverfahren für Satelliten |
| CN1789082A (zh) * | 2005-12-16 | 2006-06-21 | 西安电子科技大学 | 在轨地球同步自旋卫星红外弦宽差分定姿方法 |
| CN103162711A (zh) * | 2013-02-28 | 2013-06-19 | 北京航空航天大学 | 一种高动态星敏感器像增强器的误差补偿方法和装置 |
| US20150298827A1 (en) * | 2014-04-22 | 2015-10-22 | Massachusetts Institute Of Technology | Attitude determination using infrared earth horizon sensors |
| US9809328B2 (en) * | 2014-04-22 | 2017-11-07 | Massachusetts Institute Of Technology | Attitude determination using infrared earth horizon sensors |
| CN105572692A (zh) * | 2015-12-16 | 2016-05-11 | 上海卫星工程研究所 | 卫星图像导航与配准全物理试验装置及测试方法 |
| CN106027904A (zh) * | 2016-06-29 | 2016-10-12 | 中国人民解放军国防科学技术大学 | 视频卫星姿态运动图像补偿方法 |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107292831A (zh) * | 2017-04-24 | 2017-10-24 | 上海卫星工程研究所 | 静止卫星成像导航与配准的姿态运动补偿方法 |
| CN107329191A (zh) * | 2017-05-17 | 2017-11-07 | 上海卫星工程研究所 | 静止气象卫星图像导航与配准的成像试验系统及试验方法 |
| CN107329191B (zh) * | 2017-05-17 | 2020-04-21 | 上海卫星工程研究所 | 静止气象卫星图像导航与配准的成像试验系统及试验方法 |
| CN109141368A (zh) * | 2018-07-19 | 2019-01-04 | 北京空间机电研究所 | 一种高定位精度星载成像系统及方法 |
| CN109141368B (zh) * | 2018-07-19 | 2020-08-14 | 北京空间机电研究所 | 一种高定位精度星载成像系统及方法 |
| RU2798008C2 (ru) * | 2019-02-11 | 2023-06-14 | Виасат, Инк. | Дельта-кодирование для дистанционного зондирования |
| CN113063434A (zh) * | 2021-02-25 | 2021-07-02 | 上海卫星工程研究所 | 卫星指向恒星的精度评估方法及系统 |
| CN113063435A (zh) * | 2021-02-25 | 2021-07-02 | 上海卫星工程研究所 | 卫星姿态稳定度和指向精度评估方法及系统 |
| CN113063435B (zh) * | 2021-02-25 | 2022-12-13 | 上海卫星工程研究所 | 卫星姿态稳定度和指向精度评估方法及系统 |
| CN113063434B (zh) * | 2021-02-25 | 2023-03-24 | 上海卫星工程研究所 | 卫星指向恒星的精度评估方法及系统 |
| CN113063437A (zh) * | 2021-02-26 | 2021-07-02 | 上海卫星工程研究所 | 适用于全物理仿真的卫星视线指向模拟方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106553770B (zh) | 2019-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106553770B (zh) | 遥感卫星姿态运动补偿的成像试验方法 | |
| CN105572692B (zh) | 卫星图像导航与配准全物理试验装置及测试方法 | |
| Stöcker et al. | Quality assessment of combined IMU/GNSS data for direct georeferencing in the context of UAV-based mapping | |
| CN112629431B (zh) | 土木结构变形监测方法及相关设备 | |
| CN106767714B (zh) | 提高卫星图像定位精度的等效失配模型多级标校方法 | |
| US9194954B2 (en) | Method for geo-referencing an imaged area | |
| CN104006787B (zh) | 空间飞行器姿态运动模拟平台高精度姿态确定方法 | |
| CN107014399B (zh) | 一种星载光学相机-激光测距仪组合系统联合检校方法 | |
| CN107144293A (zh) | 一种视频卫星面阵相机的几何定标方法 | |
| CN111076880A (zh) | 一种考虑相机姿态变化的大跨桥梁多点挠度测量方法 | |
| CN105444778B (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
| CN105758400B (zh) | 静止卫星成像导航与配准恒星敏感东西参数提取方法 | |
| CN110779688B (zh) | 一种大视场面阵器件视场拼接精度测试方法 | |
| CN105812791B (zh) | 一种光学跟踪测量数据系统误差补偿方法 | |
| CN111473802A (zh) | 一种基于线阵推扫的光学传感器内方位元素定标方法 | |
| CN107633108A (zh) | 静止卫星成像导航与配准的轨道运动补偿方法 | |
| Liu et al. | Accurate mapping method for UAV photogrammetry without ground control points in the map projection frame | |
| CN110766757A (zh) | 一种具有二维指向镜的面阵相机几何成像模型标定方法 | |
| CN109855822A (zh) | 一种基于无人机的高铁桥梁竖向动扰度测量方法 | |
| Yuan et al. | A precise calibration method for line scan cameras | |
| CN107329191B (zh) | 静止气象卫星图像导航与配准的成像试验系统及试验方法 | |
| Albayrak et al. | The use of the automated digital zenith camera system in Istanbul for the determination of astrogeodetic vertical deflection | |
| Schwind et al. | An in-depth simulation of EnMAP acquisition geometry | |
| CN108956101A (zh) | 一种相机视轴变化的测量装置和测量方法 | |
| CN112097794B (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 |