CN111983612B - 一种sar滑动聚束模式方位去斜方法 - Google Patents
一种sar滑动聚束模式方位去斜方法 Download PDFInfo
- Publication number
- CN111983612B CN111983612B CN202010874366.7A CN202010874366A CN111983612B CN 111983612 B CN111983612 B CN 111983612B CN 202010874366 A CN202010874366 A CN 202010874366A CN 111983612 B CN111983612 B CN 111983612B
- Authority
- CN
- China
- Prior art keywords
- earth
- azimuth
- satellite
- coordinate system
- position vector
- 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
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9052—Spotlight mode
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种SAR滑动聚束模式方位去斜方法。方法包括:根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角,确定每个方位时刻对应的多普勒频率;根据每个方位时刻对应的多普勒频率,确定成像时间内的多普勒历程;根据多普勒历程,确定成像时间内的相对距离历程;根据相对距离历程,确定方位去斜函数;将回波信号与方位去斜函数相乘,以对回波信号进行方位去斜。本发明提供的方位去斜方法不依赖虚拟转动点,可以解决大方位扫描角滑动聚束模式虚拟转动点模型失效的问题,同时可以精确去除波束方位扫描引起的多普勒混叠。
Description
技术领域
本发明涉及合成孔径雷达(SAR)数据处理技术领域,具体涉及一种SAR滑动聚束模式方位去斜方法。
背景技术
滑动聚束模式是先进体制高分宽幅SAR最为常用的一种成像模式,在滑动聚束模式下,天线波束中心的交点在地平面下方的某一点处,这一点被称为虚拟转动点。现有的频域成像算法的第一步通常是对回波信号进行方位去斜处理,采用的方位去斜方法通常先计算虚拟转动点到卫星的距离历程,再利用该距离历程实现方位去斜。可以看出,现有的方位去斜方法依赖虚拟转动点。
但是,随着SAR方位幅宽和方位分辨率的提高,天线波束方位扫描角范围迅速增大,匀速扫描控制下波束中心会明显偏离虚拟转动点,导致虚拟转动点模型失效,此时,无法继续利用虚拟转动点进行方位去斜,使得基于方位去斜的频域成像算法也难以继续适用。
因此,亟需一种不依赖虚拟转动点的方位去斜方法,以解决大方位扫描角滑动聚束模式下的虚拟转动点模型失效的问题。
发明内容
本发明的目的是提供一种SAR滑动聚束模式方位去斜方法,以解决大方位扫描角滑动聚束模式下的虚拟转动点模型失效的问题。
为达到上述目的,本发明实施例提供一种SAR滑动聚束模式方位去斜方法,方法包括:
根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角,确定每个方位时刻对应的多普勒频率;
根据所述每个方位时刻对应的多普勒频率,确定所述成像时间内的多普勒历程;
根据所述多普勒历程,确定所述成像时间内的相对距离历程;
根据所述相对距离历程,确定方位去斜函数;
将回波信号与所述方位去斜函数相乘,以对所述回波信号进行方位去斜。
由以上本发明实施例提供的技术方案可见,本发明提供的方位去斜方法不再依赖虚拟转动点,通过波束中心地表交点计算瞬时多普勒历程,再积分得到相对距离历程,就可以构建方位去斜函数,从而实现方位去斜处理。此外,基于本发明提供的方位去斜方法,SAR系统设计也不再需要严格按照虚拟转动点模型,方位波束可以采用更加简单的匀速扫描,使得系统设计和实现变得更加灵活和简单,并且本发明提供的方位去斜方法还能够精确地去除波束方位扫描引起的多普勒混叠。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是依赖虚拟转动点的滑动聚束模式示意图;
图2是不依赖虚拟转动点的滑动聚束模式示意图;
图3是本公开实施例提供的方位去斜方法的流程图;
图4是本公开实施例提供的方位去斜方法的算法流程框图。
具体实施方式
下面将结合附图和具体实施方式,对本发明的技术方案作详细说明,应理解这些实施方式仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落入本发明所附权利要求所限定的范围内。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施方式的目的,不是旨在于限制本发明。
参考图1所示,为依赖虚拟转动点的滑动聚束模式示意图。参考图2所示,为不依赖虚拟转动点的滑动聚束模式示意图。现有的方位去斜方法都是针对图1所示的依赖虚拟转动点的滑动聚束模式,为了使SAR数据处理不再依赖虚拟转动点,从而可以在图2所示的滑动聚束模式下进行SAR系统设计与数据处理,本公开实施例提供了一种SAR滑动聚束模式方位去斜方法,参考图3所示,该方法可以包括:
S1:根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角,确定每个方位时刻对应的多普勒频率。
其中,参考图4所示,步骤S1还可以包括子步骤S101至S104:
S101:根据所述成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角,计算在地球惯性坐标系下,卫星到波束中心地表交点的单位位置矢量、地心到卫星的位置矢量以及卫星的速度矢量。
具体的,地球惯性坐标系下的卫星到波束中心地表交点的单位位置矢量,通过以下公式计算得到:
其中,
u1=-u-π/2
i1=i-π/2
上述公式中的为地球惯性坐标系下的卫星到波束中心地表交点的单位位置矢量,为天线坐标系下的卫星到波束中心地表交点的单位位置矢量,Aov为地球惯性坐标系到轨道坐标系的转换矩阵的逆矩阵,Ave为轨道坐标系到本体坐标系的转换矩阵的逆矩阵,Aea为本体坐标系到天线坐标系的转换矩阵的逆矩阵,θlook为天线视角,θsq为天线斜视角,q为升交点赤径,θyaw为卫星的偏航角,θpitch为卫星的俯仰角,θroll为卫星的横滚角,u为卫星纬度角,i为轨道的倾角。
具体的,地球惯性坐标系下的地心到卫星的位置矢量以及卫星的速度矢量,通过以下公式计算得到:
其中,
上述公式中的为所述地球惯性坐标系下的所述地心到卫星的位置矢量,Aov为地球惯性坐标系到轨道坐标系的转换矩阵的逆矩阵,为轨道坐标系下的所述地心到卫星的位置矢量,为所述地球惯性坐标系下的卫星的速度矢量,为轨道坐标系下的卫星的速度矢量,a为轨道的半长轴,e为轨道的离心率,u为卫星纬度角,w为近地点辐角,μ为地球引力常数。
S102:根据所述单位位置矢量以及所述地心到卫星的位置矢量,计算在所述地球惯性坐标系下,卫星到波束中心地表交点的位置矢量。
具体的,地球惯性坐标系下,卫星到波束中心地表交点的位置矢量,可以通过以下公式计算得到:
其中,
上述公式中的为所述地球惯性坐标系下的卫星到波束中心地表交点的位置矢量,为卫星到波束中心地表交点的位置矢量的模值,为所述地球惯性坐标系下的卫星到波束中心地表交点的单位位置矢量,以及分别为在地球惯性坐标系的x轴、y轴、z轴上的投影,为地球惯性坐标系下的所述地心到卫星的位置矢量,Ea为地球的长半轴,Eb为地球的短半轴, 以及分别为在地球惯性坐标系的x轴、y轴、z轴上的投影。
S103:根据所述卫星到波束中心地表交点的位置矢量以及所述地心到卫星的位置矢量,计算在所述地球惯性坐标系下,地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量。
具体的,在地球惯性坐标系下,地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量,通过以下公式计算得到:
上述公式中的为所述地球惯性坐标系下的地心到波束中心地表交点的位置矢量,为所述地球惯性坐标系下的所述地心到卫星的位置矢量,为所述地球惯性坐标系下的卫星到波束中心地表交点的位置矢量,为所述地球惯性坐标系下的波束中心地表交点的速度矢量,ωe为地球的自转角速度。
S104:根据所述地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量,以及所述地心到卫星的位置矢量与卫星的速度矢量,计算所述方位时刻对应的多普勒频率。
具体的,方位时刻对应的多普勒频率通过以下公式计算得到:
上述公式中的fa为所述多普勒频率,为所述地球惯性坐标系下的所述地心到卫星的位置矢量,为所述地球惯性坐标系下的地心到波束中心地表交点的位置矢量,为所述地球惯性坐标系下的卫星的速度矢量,为所述地球惯性坐标系下的波束中心地表交点的速度矢量,λ为SAR信号波长。
S2:根据所述每个方位时刻对应的多普勒频率,确定所述成像时间内的多普勒历程。
具体的,通过多项式拟合,计算多普勒历程的各阶系数,以得到所述多普勒历程:
上述公式中的fa(ta)为所述多普勒历程,fdc、fdr1、fdr2、fdr3以及fdr4分别为各阶系数,ta为SAR方位时刻。
S3:根据所述多普勒历程,确定所述成像时间内的相对距离历程。
具体的,根据相对距离历程和多普勒历程的关系表达式:
fa(ta)=-2ΔR′a(ta)/λ
通过积分可以计算得到相对距离历程:
上述公式中的ΔRa(ta)为所述相对距离历程,fdc、fdr1、fdr2、fdr3以及fdr4分别为所述多普勒历程中的各阶系数,ta为SAR方位时刻,λ为SAR信号波长。
S4:根据所述相对距离历程,确定方位去斜函数。
具体的,方位去斜函数为距离时域的方位去斜函数,或距离频域的方位去斜函数;其中,距离时域的方位去斜函数通过以下公式计算得到:
距离频域的方位去斜函数通过以下公式计算得到:
上述公式中的Hdechirp(ta)为距离时域的方位去斜函数,ta为SAR方位时刻,λ为SAR信号波长,ΔRa(ta)为所述相对距离历程,Hdechirp(fτ,ta)为距离频域的方位去斜函数,fτ为距离频率,c为光速,fc为SAR信号载频。
S5:将回波信号与所述方位去斜函数相乘,以对所述回波信号进行方位去斜。
由以上本发明实施例提供的技术方案可见,本发明提供的方位去斜方法不再依赖虚拟转动点,通过波束中心地表交点计算瞬时多普勒历程,再积分得到相对距离历程,就可以构建方位去斜函数,从而实现方位去斜处理。此外,基于本发明提供的方位去斜方法,SAR系统设计也不再需要严格按照虚拟转动点模型,方位波束可以采用更加简单的匀速扫描,使得系统设计和实现变得更加灵活和简单。并且本发明提供的方位去斜方法还能够精确地去除波束方位扫描引起的多普勒混叠。
本说明书中的上述各个实施方式均采用递进的方式描述,各个实施方式之间相同相似部分相互参照即可,每个实施方式重点说明的都是与其他实施方式不同之处。
以上所述仅为本发明的几个实施方式,虽然本发明所揭露的实施方式如上,但所述内容只是为了便于理解本发明的技术方案而采用的实施方式,并非用于限定本发明。任何本发明所属技术领域的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施方式的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附权利要求书所界定的范围为准。
Claims (9)
1.一种SAR滑动聚束模式方位去斜方法,其特征在于,包括:
根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角,确定每个方位时刻对应的多普勒频率;
根据所述每个SAR方位时刻对应的多普勒频率,确定所述成像时间内的多普勒历程;
根据所述多普勒历程,确定所述成像时间内的相对距离历程;
根据所述相对距离历程,确定方位去斜函数;
将回波信号与所述方位去斜函数相乘,以对所述回波信号进行方位去斜;
其中,所述根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态以及天线扫描角,确定每个方位时刻对应的多普勒频率,包括以下步骤:
根据所述成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角,计算在地球惯性坐标系下,卫星到波束中心地表交点的单位位置矢量、地心到卫星的位置矢量以及卫星的速度矢量;
根据所述单位位置矢量以及所述地心到卫星的位置矢量,计算在所述地球惯性坐标系下,卫星到波束中心地表交点的位置矢量;
根据所述卫星到波束中心地表交点的位置矢量以及所述地心到卫星的位置矢量,计算在所述地球惯性坐标系下,地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量;
根据所述地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量,以及所述地心到卫星的位置矢量与卫星的速度矢量,计算所述方位时刻对应的多普勒频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010874366.7A CN111983612B (zh) | 2020-08-26 | 2020-08-26 | 一种sar滑动聚束模式方位去斜方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010874366.7A CN111983612B (zh) | 2020-08-26 | 2020-08-26 | 一种sar滑动聚束模式方位去斜方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111983612A CN111983612A (zh) | 2020-11-24 |
CN111983612B true CN111983612B (zh) | 2022-04-15 |
Family
ID=73440586
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010874366.7A Active CN111983612B (zh) | 2020-08-26 | 2020-08-26 | 一种sar滑动聚束模式方位去斜方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111983612B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115792906B (zh) * | 2023-02-02 | 2023-04-11 | 中国科学院空天信息创新研究院 | 一种星载大斜视滑动聚束sar成像处理方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2743727A2 (en) * | 2014-01-16 | 2014-06-18 | Institute of Electronics, Chinese Academy of Sciences | Method for implementing high-resolution wide-swath spaceborne SAR system |
CN104062659A (zh) * | 2014-07-11 | 2014-09-24 | 哈尔滨工业大学 | 基于压缩感知的滑动聚束sar三维成像方法 |
EP2873988A1 (en) * | 2013-11-14 | 2015-05-20 | Korea Aerospace Research Institute | Method for SAR processing |
CN105974409A (zh) * | 2016-01-13 | 2016-09-28 | 合肥工业大学 | 基于多频子带并发的星载滑动聚束mimo-sar成像方法 |
CN106226768A (zh) * | 2016-08-09 | 2016-12-14 | 北京空间飞行器总体设计部 | 超高分辨率敏捷sar卫星滑动聚束模式系统参数设计方法 |
CN106526553A (zh) * | 2016-10-31 | 2017-03-22 | 北京空间飞行器总体设计部 | 一种通用的sar卫星方位模糊度性能精确分析方法 |
CN106872979A (zh) * | 2017-03-02 | 2017-06-20 | 中国科学院电子学研究所 | 基于海面起伏目标sar图像再聚焦的海浪参数获取方法 |
CN107741585A (zh) * | 2017-09-29 | 2018-02-27 | 西安电子科技大学 | 适用于大斜视sar成像的距离空变md自聚焦方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103149564B (zh) * | 2013-01-29 | 2014-12-31 | 内蒙古工业大学 | 斜视偏置相位中心方位多波束星载sar 成像方法 |
US9182483B2 (en) * | 2013-03-15 | 2015-11-10 | Mitsubishi Electric Research Laboratories, Inc. | Method and system for random steerable SAR using compressive sensing |
-
2020
- 2020-08-26 CN CN202010874366.7A patent/CN111983612B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2873988A1 (en) * | 2013-11-14 | 2015-05-20 | Korea Aerospace Research Institute | Method for SAR processing |
EP2743727A2 (en) * | 2014-01-16 | 2014-06-18 | Institute of Electronics, Chinese Academy of Sciences | Method for implementing high-resolution wide-swath spaceborne SAR system |
CN104062659A (zh) * | 2014-07-11 | 2014-09-24 | 哈尔滨工业大学 | 基于压缩感知的滑动聚束sar三维成像方法 |
CN105974409A (zh) * | 2016-01-13 | 2016-09-28 | 合肥工业大学 | 基于多频子带并发的星载滑动聚束mimo-sar成像方法 |
CN106226768A (zh) * | 2016-08-09 | 2016-12-14 | 北京空间飞行器总体设计部 | 超高分辨率敏捷sar卫星滑动聚束模式系统参数设计方法 |
CN106526553A (zh) * | 2016-10-31 | 2017-03-22 | 北京空间飞行器总体设计部 | 一种通用的sar卫星方位模糊度性能精确分析方法 |
CN106872979A (zh) * | 2017-03-02 | 2017-06-20 | 中国科学院电子学研究所 | 基于海面起伏目标sar图像再聚焦的海浪参数获取方法 |
CN107741585A (zh) * | 2017-09-29 | 2018-02-27 | 西安电子科技大学 | 适用于大斜视sar成像的距离空变md自聚焦方法 |
Non-Patent Citations (6)
Title |
---|
Investigation on the Space-Time Waveform Encoding Single Phase Center Multiple Azimuth Beam sliding spotlight SAR for ultra-high resolution and wide swath imaging;Liu Yadong;《IET International Radar Conference 2015》;20160430;1-9 * |
On the Processing of Very High Resolution Spaceborne SAR Data: A Chirp-Modulated Back Projection Approach;Dadi Meng;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING,》;20180131;第56卷(第1期);191-201 * |
卫星姿态误差对多通道SAR成像质量的影响;吴亮等;《测绘通报》;20150125(第01期);全文 * |
双基前视髙机动平台SAR系统特性及成像算法研究;孟自强;《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》;20170215;20-21 * |
变波门大斜视滑动聚束SAR成像关键技术分析;聂鑫;《电子与信息学报》;20161231;第38卷(第12期);3122-3128 * |
星载干涉合成孔径雷达信号处理研究;申艳;《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》;20081015;34 * |
Also Published As
Publication number | Publication date |
---|---|
CN111983612A (zh) | 2020-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108508918B (zh) | 一种静轨遥感卫星数传天线高精度实时对地指向控制方法 | |
Davidson et al. | Signal properties of spaceborne squint-mode SAR | |
Dai et al. | Estimating the yaw-attitude of BDS IGSO and MEO satellites | |
CN107689064B (zh) | 顾及光行差改正的卫星光学严格几何成像模型构建方法 | |
CN111381256B (zh) | 主动遥感卫星天线相位中心偏移误差计算的方法和系统 | |
US20090012662A1 (en) | Method and apparatus for determining a satellite attitude using crosslink reference signals | |
CN101226058A (zh) | 一种卫星侧摆机动图像实时传输的实现方法 | |
CN112130590B (zh) | 一种基于瞬时惯性系下速度补偿的星载天线对地指向确定方法 | |
CN111983612B (zh) | 一种sar滑动聚束模式方位去斜方法 | |
CN113589285A (zh) | 一种飞行器sar实时成像方法 | |
CN103675760A (zh) | 一种星载地球同步轨道合成孔径雷达姿态导引方法 | |
CN113834482A (zh) | 定位方法、装置、电子设备及存储介质 | |
CN103344958B (zh) | 基于星历数据的星载sar高阶多普勒参数估算方法 | |
CN112713922A (zh) | 一种多波束通讯卫星的可见性快速预报算法 | |
JP2003322676A (ja) | 衛星追跡方法及び衛星追跡装置 | |
CN112255606A (zh) | 基于单反射面天线的Geo-SAR卫星正侧视成像姿态角计算方法 | |
CN111541477B (zh) | 低轨星座系统内部频率干扰抑制方法及装置 | |
JP4489654B2 (ja) | 衛星追尾用のアンテナ制御装置 | |
CN108489483B (zh) | 一种船载星光定向仪单星次优修正算法 | |
Ward | Spacecraft attitude estimation using GPS: Methodology and results | |
US6283415B1 (en) | Simplified yaw steering method for satellite antenna beam control | |
JP4295618B2 (ja) | Gpsの初期化のための衛星姿勢調整 | |
Li et al. | Adaptively robust filtering algorithm for maritime celestial navigation | |
KR102195422B1 (ko) | 안테나 제어 방법 및 장치 | |
CN114578353A (zh) | 一种sar波束时空覆盖特性分析方法及装置 |
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 |