CN109884621B - 雷达高度表回波相参积累方法 - Google Patents
雷达高度表回波相参积累方法 Download PDFInfo
- Publication number
- CN109884621B CN109884621B CN201910153944.5A CN201910153944A CN109884621B CN 109884621 B CN109884621 B CN 109884621B CN 201910153944 A CN201910153944 A CN 201910153944A CN 109884621 B CN109884621 B CN 109884621B
- Authority
- CN
- China
- Prior art keywords
- echo
- angle
- wide
- altimeter
- signal
- 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
Abstract
本发明提供了一种雷达高度表回波相参积累方法,本发明构建了宽角度下的回波信号模型;利用多项式模型,采用非线性映射的方法,对接收到的宽角度回波信号进行相参累积,从而改善了回波信号的波形;利用相参处理后的信号,实现了宽角度SAR高度表的测高。
Description
技术领域
本发明涉及一种雷达高度表回波相参积累方法。
背景技术
雷达高度表可以通过接收和处理散射回波来获取观测的地面或海平面高度。传统的雷达高度表只能处理天线当前波束内的回波,这就导致测高的精度和方位向的空间分辨率较低,分别在10cm和1km的量级。合成孔径雷达高度表借鉴了传统合成孔径雷达(SAR)的虚拟合成孔径技术,将沿航迹向对地面目标多个波束观测的回波信号进行积累,能将测高精度和空间分辨率分别提高到2cm~5cm和200m~500m。然而,目前合成孔径雷达高度表在方位向的处理并不是完全聚焦的,通常发射若干个脉冲簇,脉冲簇内部进行相干合成而脉冲簇之间进行非相干线性叠加(类似于传统SAR中的子孔径处理)。针对这种情况,2017年A.Egido和W.Smith提出了完全聚焦合成孔径雷达高度表(fully focused SARaltimetry)。这一新的进展实现了脉冲簇间的精确相位补偿,将方位向空间分辨率提高到理论值即天线长度的一半,同时测高精度也得到进一步提升。从提高空间分辨率的角度去审视雷达高度表的发展脉络,与SAR技术的发展一脉相承,而近年来处于研究前沿的宽角度SAR不仅能获得超高的方位向分辨率(方位向360度合成的圆迹SAR理论分辨率为1/4雷达波长),而且能获得丰富的目标特性。对比类推,宽角度合成孔径雷达高度表这一现实应用需求也有望获得更高的空间分辨率和更高的测高精度,是雷达高度表的重要发展方向。
尽管如此,在宽角度SAR高度表的研究领域仍然有许多困难需要解决。首先,观测场景及目标的后向散射特性会随着方位角的变化,如果不考虑这种变化,将会导致回波的直接累积出现误差。因此,宽角度回波需要在某种补偿与变换之后才能累积。其次,地面场景的复杂性会进一步增加信号模型和处理算法的难度。而目前已有的宽角度SAR高度表回波的处理方法,要么没有充分考虑散射特性的影响,要么只能利用统计学方法研究海平面等场景的散射特性。
发明内容
本发明的目的在于提供一种雷达高度表回波相参积累方法。
为解决上述问题,本发明提供一种雷达高度表回波相参积累方法,包括:
步骤A,建立宽角度SAR雷达高度表回波模型;
步骤B,基于多项式模型的宽角度回波非线性映射与相参累积;
步骤C,基于参数估计的宽角度SAR雷达高度表进行高度估计。
进一步的,在上述方法中,所述步骤A包括:
步骤11,按照预设周期发送线性调频信号;
步骤12,获取预设角度范围内的目标散射回的宽角度SAR雷达高度表的回波信号;
步骤13,对接收到的回波信号进行解调,去掉其中的载频项;
步骤14,对所述回波信号的脉冲压缩,得到处理后的快时间/慢时间域下的回波信号Eout(t,tn),t,tn分别为快时间、慢时间。
进一步的,在上述方法中,所述步骤B,包括:
步骤21,将所述回波信号Eout(t,tn)离散化,并得到总功率Ps(t);
步骤22,对回波功率进行延迟补偿:
步骤23,求解优化问题:
其中,Θ为多项式系数构成的MN×1维向量,am,n为多项式系数,M和N分别是角度域和时间域的采样数;
步骤24,利用求解出来的多项式系数,将原回波信号与多项式函数相乘,从而得到非线性映射后的信号:
步骤25,将映射后的信号进行相参累积:
进一步的,在上述方法中,所述步骤C包括求解优化问题:
其中
gk(H)=pk(H)-sk(H),
s=(s1,s2,......sK)是离散化的回波功率向量,p=(p1,p2,......pK)是含有噪声的实测信号,依然采用高斯牛顿法求解该问题,以得到H的估计值。
与现有技术相比,针对SAR高度表现有回波信号处理方法在宽角度场景下的缺陷,本发明结采用多项式建模的方法,将多角度回波信号映射到同一角度再进行累积。随后利用累积后的信号估计高度表的高度等参数,从而实现高度表的测高功能。本发明构建了宽角度下的回波信号模型;利用多项式模型,采用非线性映射的方法,对接收到的宽角度回波信号进行相参累积,从而改善了回波信号的波形;利用相参处理后的信号,实现了宽角度SAR高度表的测高。
附图说明
图1为本发明一实施例的示意图;
图2a、2b为本发明一实施例的雷达平台运动示意图与慢时间-方位角示意图;
图3为本发明一实施例的圆柱体散射函数的示意图;
图4为不同多普勒段圆柱体散射回波的示意图;
图5a、5b为非线性映射前后的圆柱体散射回波的示意图;
图6为长方体散射函数的示意图;
图7为长方体散射回波的示意图;
图8a、8b为非线性映射前后的长方体散射回波的示意图;
图9a、9b为二面角散射函数及散射回波的示意图;
图10a、10b为非线性映射前后的二面角散射回波的示意图;
图11为相参积累后的圆柱体散射回波的示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
解决宽角度SAR高度表的回波信号处理问题,依然需要从SAR高度表回波信号模型的构建出发。Brown将雷达高度表回波建模为三项卷积的形式,即平坦表面脉冲响应函数(FSIR)、地面高度概率密度函数(PDF)和雷达系统目标响应函数(PTR)的卷积,此后有大量学者对该模型进行了改进以适用不同应用场景,以及将其推广到合成孔径雷达高度表中。
为了考虑散射特性随角度变化对回波信号处理的影响,基于物理光学和几何绕射理论,研究宽角度回波与角度的函数模型,并用散射函数替代三项卷积模型中恒定的散射系数,即可获得宽角度下的SAR高度表回波信号模型。
本发明提供一种基于电磁散射参数化模型的宽角度的雷达高度表回波相参积累方法,包括如下步骤:
步骤A,建立宽角度SAR雷达高度表回波模型;
步骤B,基于多项式模型的宽角度回波非线性映射与相参累积;
步骤C,基于参数估计的宽角度SAR雷达高度表进行高度估计。
在此,针对SAR高度表现有回波信号处理方法在宽角度场景下的缺陷,本发明结采用多项式建模的方法,将多角度回波信号映射到同一角度再进行累积。随后利用累积后的信号估计高度表的高度等参数,从而实现高度表的测高功能。本发明构建了宽角度下的回波信号模型;利用多项式模型,采用非线性映射的方法,对接收到的宽角度回波信号进行相参累积,从而改善了回波信号的波形;利用相参处理后的信号,实现了宽角度SAR高度表的测高。
本发明的雷达高度表回波相参积累方法一实施例中,步骤A具体为:
步骤11,按照预设周期发送线性调频信号;
步骤12,获取预设角度范围内的目标散射回的宽角度SAR雷达高度表的回波信号;
步骤13,对接收到的回波信号进行解调,去掉其中的载频项;
步骤14,对所述回波信号的脉冲压缩,得到处理后的快时间/慢时间域下的回波信号Eout(t,tn),t,tn分别为快时间,慢时间。
本发明的雷达高度表回波相参积累方法一实施例中,步骤B具体为:
步骤21,将所述回波信号Eout(t,tn)离散化,并得到总功率Ps(t);
步骤22,因为宽角度SAR高度表需要在飞行轨迹的不同观测点对同一目标进行多次观测,这会造成距离迁徙现象,因此需要对回波功率进行延迟补偿:
步骤23,求解优化问题:
其中Θ为多项式系数构成的MN×1维向量,am,n为多项式系数,M和N分别是角度域和时间域的采样数。
步骤24,利用求解出来的多项式系数,将原回波信号与多项式函数相乘,从而得到非线性映射后的信号:
步骤25,将映射后的信号进行相参累积:
本发明的雷达高度表回波相参积累方法一实施例中,步骤C具体为,求解优化问题:
其中
gk(H)=pk(H)-sk(H)
s=(s1,s2,......sK)是离散化的回波功率向量,p=(p1,p2,......pK)是含有噪声的实测信号。依然采用高斯牛顿法求解该问题,便可以得到H的估计值。
具体的,
一、宽角度SAR高度表回波模型的建立
如图1所示,步骤11,按照一定周期发送线性调频信号:
调制后的发射信号为:
步骤12,获取一定角度范围内的雷达高度表回波:
步骤13,对接收到的信号进行解调,去掉其中的载频项,即:
其中λ是发射信号的波长,Lp是双程传播损失系数,v是雷达平台运动速度,G(t)是天线增益,γ(tn)是目标散射函数,X是合成孔径长度,R是雷达高度表和目标之间的距离:
雷达天线方向图可近似为一个sinc函数:
方位向波束跨度βbw=0.886λ/L,L为方位向天线长度。由于雷达能量的双程传播,接收信号的强度由p(θ)的平方给出,即
步骤14,对回波信号脉冲压缩,其中匹配滤波器为:
之后,脉冲压缩处理后的信号表达式为:
传统的SAR高度表只考虑角度范围较小的情况,因此γ是常数,而宽角度SAR高度表中的散射系数随角度变化,是一个函数,并且是已知的。为了简化模型,如图2a、2b所示,假设雷达平台沿着x轴方向运动,当信号频率确定时,目标散射函数由方位角决定,与俯仰角无关。而又仅由慢时间tn决定。它们的关系为:
二、基于多项式模型的宽角度回波非线性映射与相参累积
步骤21,将回波信号离散化,总功率表达式为:
步骤12,对回波功率进行延迟补偿:
步骤13,将不同方位角对应的回波信号与多项式函数相乘并将其变换到零多普勒段上,求解优化问题:
其中Θ为多项式系数构成的MN×1维向量,am,n为多项式系数,M和N分别是角度域和时间域的采样数。高斯牛顿法求解上述优化问题的步骤如下所示:
初始化:
初始点设为Θ(0),最大允许误差为ε,设k=0
迭代过程:
while(true):
If G(Θ(k+1))-G(Θ(k))<ε:
break
else:
Θ(k+1)=Θ(k)-(JTJ)-1JTg(Θ)
k=k+1.
end
步骤14,利用求解出来的多项式系数,将原回波信号与多项式函数相乘,从而得到非线性映射后的信号:
步骤15,将映射后的信号进行相参累积:
三、基于参数估计的宽角度SAR高度表高度估计
求解优化问题:
其中
gk(H)=pk(H)-sk(H)
s=(s1,s2,......sK)是离散化的回波功率向量,p=(p1,p2,……pK)是含有噪声的实测信号。依然采用高斯牛顿法求解该问题,便可以得到H的估计值。
四、仿真与结果
1、参数设置
以星载雷达高度表为例,以圆柱体、长方体、二面角结构为目标进行仿真,仿真的内容包括两部分,第一部分利用多项式函数对回波进行非线性映射,完成对多项式系数的求解;第二部分利用模拟出的接收到的回波信号对目标参数和雷达高度进行求解。仿真参数如表1所示:
表1仿真参数
仿真中默认的方位角范围为(-π/2,π/2)。
2、散射回波非线性映射
1)圆柱体
一个高度为200m,半径为100m的圆柱体散射函数随方位角变化的规律如图3所示,其散射回波如图4所示。可以看出,与随角度变化的散射函数相乘后,相当于在多普勒域对波形做了加权处理。
在仿真中,选取的多项式系数M,N分别为3,利用高斯牛顿法迭代得到多项式系数后,可得到非线性映射后的信号波形。图5a、5b将映射前后的波形在多普勒域进行了比较,其中横坐标为多普勒段,对应着方位角,可以看出,非线性映射后的信号功率的分布更加集中于零多普勒点,这也就达到了进行非线性映射的最初目的。
2)长方体
一个长、宽、高分别为200m、150m、100m的长方体的散射函数和回波分别如图6、7所示。非线性映射前后的散射回波如图8a、8b所示。最终得到的结果与圆柱体仿真类似,这里就不再赘述。
3)二面角
对一个公共边为200m,垂直边分别为200m、150m的二面角进行仿真,从图9a、9b可以看出,二面角结构的散射函数形式较为复杂,因此有必要增加多项式的项数以提高模型的拟合能力。经过反复验证,在M和N的值都取6的时候拟合效果最好。非线性映射前后的散射回波如图10a、10b所示。
3、高度估计
以圆柱体为例,利用高斯牛顿法,分别用传统的SAR高度表和宽角度SAR高度表对高度参数进行估计。参数估计的结果和误差如表2所示:
表2高度估计的结果与误差
测高结果表明,宽角度SAR高度表的估计性能优于传统的延迟/多普勒SAR高度表。这是由于宽角度SAR高度表在高度估计的过程中相当于采集了更大角度范围内的回波样本,且鉴于模型中散射函数随角度的变化而变化,求解过程考虑到这一情况,不再把散射系数看作常数。
4回波相参积累
在完成了宽角度SAR回波的非线性映射以后,可以将映射到同一方位角的回波信号叠加,从而得到相参积累后的波形,这里以圆柱体为例,如图11所示。
与非相参积累相比,相参积累后回波解决了波形失真的问题,同时信号功率的分布更加集中于待估测高度附近,这也提高了信号功率的利用率。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
显然,本领域的技术人员可以对发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包括这些改动和变型在内。
Claims (3)
1.一种雷达高度表回波相参积累方法,其特征在于,包括如下步骤:
步骤A,建立宽角度SAR雷达高度表回波模型;
步骤B,基于多项式模型的宽角度回波非线性映射与相参累积;
步骤C,基于参数估计的宽角度SAR雷达高度表进行高度估计;
所述步骤B,包括:
步骤21,将回波信号Eout(t,tn)离散化,并得到总功率Ps(t);
步骤22,对回波功率进行延迟补偿:
步骤23,求解优化问题:
其中,Θ为多项式系数构成的MN×1维向量,am,n为多项式系数,M和N分别是角度域和时间域的采样数;
步骤24,利用求解出来的多项式系数,将原回波信号与多项式函数相乘,从而得到非线性映射后的信号:
步骤25,将映射后的信号进行相参累积:
2.如权利要求1所述的雷达高度表回波相参积累方法,其特征在于,所述步骤A包括:
步骤11,按照预设周期发送线性调频信号;
步骤12,获取预设角度范围内的目标散射回的宽角度SAR雷达高度表的回波信号;
步骤13,对接收到的回波信号进行解调,去掉其中的载频项;
步骤14,对所述回波信号的脉冲压缩,得到处理后的快时间/慢时间域下的回波信号Eout(t,tn),t,tn分别为快时间、慢时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910153944.5A CN109884621B (zh) | 2019-02-28 | 2019-02-28 | 雷达高度表回波相参积累方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910153944.5A CN109884621B (zh) | 2019-02-28 | 2019-02-28 | 雷达高度表回波相参积累方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109884621A CN109884621A (zh) | 2019-06-14 |
CN109884621B true CN109884621B (zh) | 2023-01-06 |
Family
ID=66930102
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910153944.5A Active CN109884621B (zh) | 2019-02-28 | 2019-02-28 | 雷达高度表回波相参积累方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109884621B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112558497B (zh) * | 2020-12-01 | 2022-06-10 | 中国人民解放军火箭军工程大学 | 一种雷达高度表抗干扰数字仿真方法及系统 |
CN113009473B (zh) * | 2021-02-03 | 2023-08-01 | 中山大学 | 一种多波束凝视雷达低仰角目标测高方法、装置及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101581779A (zh) * | 2008-05-14 | 2009-11-18 | 中国科学院电子学研究所 | 一种层析合成孔径雷达三维成像原始回波信号生成方法 |
CN103344949A (zh) * | 2013-06-18 | 2013-10-09 | 中国人民解放军海军航空工程学院 | 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法 |
CN106646447A (zh) * | 2017-01-18 | 2017-05-10 | 武汉雷博合创电子技术有限公司 | 基于线性调频连续波的雷达目标长时间积累检测方法 |
CN106896358A (zh) * | 2017-04-27 | 2017-06-27 | 电子科技大学 | 一种基于位置旋转变换的高速目标相参积累检测方法 |
CN107329138A (zh) * | 2017-08-07 | 2017-11-07 | 电子科技大学 | 一种pd雷达距离走动校正与相参积累检测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7728757B2 (en) * | 2008-06-09 | 2010-06-01 | The Boeing Company | Autofocus for minimum entropy through multi-dimensional optimization |
US10641885B2 (en) * | 2017-07-21 | 2020-05-05 | Honeywell International Inc. | Systems and methods for measuring velocity and acceleration with a radar altimeter |
-
2019
- 2019-02-28 CN CN201910153944.5A patent/CN109884621B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101581779A (zh) * | 2008-05-14 | 2009-11-18 | 中国科学院电子学研究所 | 一种层析合成孔径雷达三维成像原始回波信号生成方法 |
CN103344949A (zh) * | 2013-06-18 | 2013-10-09 | 中国人民解放军海军航空工程学院 | 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法 |
CN106646447A (zh) * | 2017-01-18 | 2017-05-10 | 武汉雷博合创电子技术有限公司 | 基于线性调频连续波的雷达目标长时间积累检测方法 |
CN106896358A (zh) * | 2017-04-27 | 2017-06-27 | 电子科技大学 | 一种基于位置旋转变换的高速目标相参积累检测方法 |
CN107329138A (zh) * | 2017-08-07 | 2017-11-07 | 电子科技大学 | 一种pd雷达距离走动校正与相参积累检测方法 |
Non-Patent Citations (10)
Title |
---|
A Space-Time Coding Scheme With Time and Frequency Comb-Like Chirp Waveforms for MIMO-SAR;Shangwen Liu et al.;《IEEE Journal of Selected Topics in Signal Processing》;20170331;第11卷(第2期);全文 * |
Coherent ranging with Envisat radar altimeter A new perspective in analyzing altimeter data using Doppler processing;R.Abileah et al.;《Remote Sensing of Environment》;20130907;全文 * |
Height Precision of SAR Altimeter and Conventional Radar Altimeter Based on Flight Experimental Data;Lingwei Shi et al.;《IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing》;20160630;第9卷(第6期);全文 * |
一种径向匀加速目标包络徙动补偿新方法;蒋千等;《雷达科学与技术》;20130215;第11卷(第01期);全文 * |
分布式全相参雷达相参性能分析;宋靖等;《电子与信息学报》;20150131;第37卷(第1期);全文 * |
多航迹圆迹SAR三维联合稀疏成像方法;朱仕恒等;《信息技术》;20161231(第12期);全文 * |
机载SAR图像快速经纬度计算及精度分析;彭代强等;《现代雷达》;20100315;第32卷(第03期);全文 * |
空间探测雷达多目标跟踪技术研究;张硕;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180615(第6期);全文 * |
米波雷达低仰角测高的地形补偿方法;马献德等;《现代雷达》;20160930;第38卷(第9期);全文 * |
脉冲压缩多普勒雷达高度表信号处理算法;刘建新;《探测与控制学报》;20090615;第31卷(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109884621A (zh) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106093870B (zh) | 高超声速飞行器下降段的sar-gmti杂波抑制方法 | |
CN105259552B (zh) | 一种基于非线性调频信号的合成孔径雷达成像方法和装置 | |
CN107976660B (zh) | 弹载多通道雷达超低空目标分析与多径回波建模方法 | |
CN105445701B (zh) | Ddma‑mimo雷达目标的单脉冲角度估计方法 | |
CN110426707B (zh) | 基于轨道角动量的涡旋sar成像方法及成像系统 | |
CN106405552B (zh) | 基于wvd—pga算法的sar雷达目标聚焦方法 | |
CN108107430A (zh) | 一种基于分数阶傅立叶变换的舰船目标isar成像方法 | |
CN109270528B (zh) | 基于全解析距离模型的一站固定式双站sar成像方法 | |
CN104035095A (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
CN107271977B (zh) | 基于移动激励源fdtd算法的高精度sar回波仿真方法 | |
CN104950305A (zh) | 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法 | |
CN104280566A (zh) | 基于空时幅相估计的低空风切变风速估计方法 | |
CN109884621B (zh) | 雷达高度表回波相参积累方法 | |
CN109085556B (zh) | 一种基于一二阶峰比值的高频地波雷达浪场形成方法 | |
CN105738887A (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
CN106054187A (zh) | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 | |
CN113376597A (zh) | 基于数字高程地图和gpu的复杂地形电磁散射快速仿真方法 | |
CN107607945A (zh) | 一种基于空间嵌入映射的扫描雷达前视成像方法 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN110244300B (zh) | 基于球体模型和fenlcs算法的弹载sar平飞段高分辨率成像方法 | |
CN106291547B (zh) | 基于天线方向图辅助的多普勒模糊分量自适应抑制方法 | |
CN105866743A (zh) | 基于混响室的雷达杂波频谱特性模拟系统及方法 | |
CN104020465B (zh) | 基于八单元小孔径圆阵天线的外辐射源雷达测角方法 | |
CN110297237B (zh) | 考虑天线方向图的探地雷达绕射叠加成像方法及系统 | |
CN110441771B (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 |