CN116299562A - 一种高度计测距电离层误差校正滤波处理方法 - Google Patents
一种高度计测距电离层误差校正滤波处理方法 Download PDFInfo
- Publication number
- CN116299562A CN116299562A CN202310601000.6A CN202310601000A CN116299562A CN 116299562 A CN116299562 A CN 116299562A CN 202310601000 A CN202310601000 A CN 202310601000A CN 116299562 A CN116299562 A CN 116299562A
- Authority
- CN
- China
- Prior art keywords
- error correction
- signal
- threshold
- frequency
- ionosphere
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 58
- 238000001914 filtration Methods 0.000 title claims abstract description 34
- 239000005433 ionosphere Substances 0.000 title claims abstract description 30
- 238000003672 processing method Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000012804 iterative process Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000007873 sieving Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 6
- 230000007547 defect Effects 0.000 description 5
- 238000011282 treatment Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000013450 outlier detection Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
- G01S19/072—Ionosphere corrections
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种高度计测距电离层误差校正滤波处理方法,属于卫星图像处理技术领域。本发明将小波去噪与类EEMD算法相结合,以分解得到的IMF1作为基函数,使用小波去噪的方法分离高频噪声与可能的地球物理信号,且因为IMF1是由数据本身所分解得到的,其有直观、后验和自适应的优点。本发明采用了一种改进的阈值函数,将其应用于IMF1的小波去噪过程,并将上述组合算法应用于双频电离层误差校正的滤波处理,以消除仪器精度差异造成的噪声,使校正结果更加接近于实际值。
Description
技术领域
本发明属于卫星图像数据处理技术领域,具体涉及一种高度计测距电离层误差校正滤波处理方法。
背景技术
卫星高度计测高通过发射已知功率的短脉冲,记录脉冲信号发射后经过粗糙海平面反射回到高度计的传播时间来计算高度计到海面的距离。在脉冲信号穿越电离层的过程中,由于传播路径上带电粒子的影响,信号传播路径发生弯曲进而造成电磁波传播速度减慢,使得基于信号传播时间和光速计算得到的传播路径并非信号出发点到信号接收点的真实几何距离,由此出现的偏差即为电离层路径延迟误差,之后简称为电离层误差。
卫星高度计得到电离层误差的方式分为两种,以高度计Sentinel-3A为例,其一是基于美国喷气推进实验室(Jet Propulsion Laboratory,简称JPL)提供的电子含量数据,经过电离层误差校正模型计算得到沿轨方向的电离层误差校正;第二种是基于高度计本身两个频段(Ku波段和C波段)的观测数据,包括基于SWH(Significant Wave Height)计算的SSB(Sea State Bias)数据以及测距数据(Range),经双频校正模型得到电离层误差校正。
然而,由于高度计C波段的测距精度低于Ku波段,所以基于两种频段测量数据得到的双频电离层误差校正会包含较大的噪声,为降低由于两频段仪器噪声引入的误差,需要对双频电离层误差校正进行低通滤波处理。目前对高度计使用的滤波处理方法一般为在一定大小的移动窗口内取平均值,如Jason系列卫星在用户手册中提供的滤波方法:根据当地时间分段滤波,在当地时间6-24时之间,使用100-150km尺度滤波;在当地时间0-6时之间,采用150-200km尺度滤波。
目前的电离层误差校正方法存在的主要不足是:滑动窗口取均值的方法变相降低了观测的分辨率,并且不能将潜在的异常值剔除,导致校正的准确性下降,同时在海冰覆盖等观测不连续处会出现空白。
发明内容
本发明的目的是一种卫星高度计双频电离层误差校正滤波处理方法,该方法基于类EEMD与改进小波去噪的联合算法,以弥补现有技术的不足。
集合经验模态分解算法(Ensemble Empirical Mode Decomposition,简称EEMD)是经验模态分解算法(Empirical Mode Decomposition,简称EMD)的改进方法,EMD被开发用于分析非平稳和非线性的信号,可以将离散的信号按比例分解为有限个调频函数,其中各个模态函数中包含着高频噪声,但是EMD方法也有相应的不足,如末端效应、模态混叠问题等,其中模态混叠会导致错误的IMF分量,从而使IMF丧失具体的物理意义,EEMD针对EMD的模态混叠问题进行了改进,通过对信号多次引入均匀分布的白噪声,对多次重构信号进行均值处理降低了模态混叠问题产生的影响。
小波去噪的关键在于小波基和阈值函数的选择,在EEMD与小波去噪的结合算法中,通常将EEMD分解后的模态函数作为小波去噪的小波基,故阈值函数的选择对滤波结果具有极大的影响,目前常见的阈值函数主要有硬阈值函数、软阈值函数、半软阈值函数、Garrote阈值函数等,常用的阈值函数有各自的缺陷,硬阈值函数的结果不连续导致去噪信号有突变、软阈值函数可能导致有用信号丢失、半软阈值函数计算阈值较为困难、Garrote改进程度有限,去噪效果仍不够理想。
EEMD算法利用白噪声均值为0的特性,通过在信号中多次引入均匀分布的白噪声,将信号本身的噪声通过多次人为添加的噪声掩盖过去,从而得到更加精准的上下包络线,再依据分解模型来进行信号分解,将复杂信号分解为有限个本征模态函数(IntrinsicMode Function,简称IMF),所分解出来的各IMF分量包含了原始信号的不同特征尺度的局部特征信息。EEMD算法均在原始信号中加入不同的白噪声,本发明基于EEMD算法,将白噪声添加在迭代得到的重构信号,称为类EEMD算法。
本发明将小波去噪与类EEMD算法相结合,以分解得到的IMF1作为基函数,使用小波去噪的方法分离高频噪声与可能的地球物理信号,且因为IMF1是由数据本身所分解得到的,这种方法有直观、后验和自适应的优点。
传统的小波去噪阈值函数有各自的缺陷,为得到去噪效果更优的阈值函数需要考虑到函数的连续性、在阈值附近能否平滑过渡、阈值函数高阶可导等条件,本发明基于上述条件,采用了一种改进的阈值函数,将其应用于IMF1的小波去噪过程,并将上述组合算法应用于双频电离层误差校正的滤波处理,以消除仪器精度差异造成的噪声,使校正结果更加接近于实际值。
为达到上述目的,并基于上述原理,本发明采取的具体技术方案为:
一种高度计测距电离层误差校正滤波处理方法,包括如下步骤:
S1:获取电离层双频校正信号数据,并对其进行预处理,得到原始信号;
S2:对原始信号使用集合经验模态分解算法将其分解为多个本征模态函数IMF和残差信号;
S3:将分解得到的IMF1作为基函数,基于小波去噪的原理首先进行阈值的选取,使用改进的阈值函数计算得到;
S4:基于改进的阈值函数进行小波去噪,低于阈值的信号波动被认为与噪声关联,故将其去除;
S5:通过将去噪后的IMF1与剩余的高阶IMF、残差信号相加,执行信号的重构;
S6:在重构信号中加入均匀分布的白噪声后将其作为原始信号输入,多次迭代过程S2到S5后将得到的所有重构信号取均值,最终得到双频电离层误差校正的滤波结果。
进一步的,所述S1中预处理为剔除缺失值、离群值检索等。
进一步的,所述S2中,确定初始双频电离层误差校正的极大值和极小值,通过插值分别对极值点进行三次样条拟合,得到基于极大、极小值的上下包络线/>、,将初始校正与上下包络线的均值相减,得到第一阶的第一剩余/>:
进一步的,所述S3中,改进的阈值函数表达式如下所示:
改进的阈值计算公式为:
当然该信号重构方法可以不止上述一种,能够达到重构效果即可。
本发明的优点和技术效果:
本发明基于EEMD算法原理,在此基础上对白噪声的添加进行了不同的处理,以达到更优的去噪效果;另外,在小波去噪的过程中采用一种改进的阈值函数进行降噪处理,弥补了传统阈值函数的不足,使去噪结果更加接近于真实信号。
本发明将类EEMD算法与改进的小波去噪相结合,应用于双频电离层误差校正的滤波处理过程,为高度计测距过程中双频电离层误差校正提供了一种新的滤波方法。
附图说明
图1为本发明的整体流程图。
图2为重构信号用类EEMD算法分解结果。
图3为双频电离层误差校正初始值与基于类EEMD算法滤波处理后结果对比图。
具体实施方式
以下通过具体实施例并结合附图对本发明进一步解释和说明。
实施例1:
一种基于类EEMD与改进的小波去噪结合的高度计双频电离层误差校正的滤波处理方法,如图1所示,包括如下步骤:
S1:自高度计2级GDR产品信号数据中提取电离层双频校正值,对其进行预处理,如剔除缺失值、离群值检索等,得到原始信号。
S2:对原始信号使用集合经验模态分解算法将其分解为多个本征模态函数(IMF)和残差信号(Residual Signal)。
S3:将分解得到的IMF1作为基函数,基于小波去噪的原理首先进行阈值的选取,使用改进的阈值函数公式计算得到。
S4:基于改进的阈值函数进行小波去噪,低于阈值的信号波动被认为与噪声关联,故将其去除。
S5:通过将去噪后的IMF1与剩余的高阶IMF、残差信号相加,执行信号的重构。
S6:在重构信号中加入均匀分布的白噪声后将其作为原始信号输入,多次迭代过程S2到S5后将得到的所有重构信号取均值,最终得到双频电离层误差校正的滤波结果。
具体如下:
在一实施例中,所述S2中:对预处理后的数据增加均匀分布的白噪声,之后再进行类EEMD分解,分解的结果如图2所示。
在一实施例中,基于上述分解的IMF1使用改进的阈值函数计算阈值。
改进的阈值函数表达式如下所示:
改进的阈值计算公式为:
计算公式为:
在一实施例中,所述S5中,对信号进行重构:
在一应用实例中,所述S1中:使用高度计Jason-3F版本2019年8月1日中位于海面的双频电离层误差校正数据(使用200个星下点),剔除其中缺失值和质量标识为1(bad)的数据,使用大小为5的滑动窗口进行离群值检测,对于窗口中的每个数据点,测试其与相邻值的差异,若该差异大于数据标准差的4.5倍,则用窗口内数据的均值代替。
以该实施例所获得的具体处理信号数据为基础,将上述方法中的处理步骤重复20次后,得到最终双频电离层误差校正的滤波结果,图3所示结果为200个连续星下点原始双频电离层误差校正数据与类EEMD与改进小波去噪算法滤波结果之间的对比图,其中实线为通过类EEMD与改进小波去噪算法滤波后的双频电离层误差校正,虚线为高度计Jason-3滤波后的双频电离层误差校正,带点的虚线为高度计Jason-3原始的双频电离层误差校正。
结果说明:由图中可以看出,Jason-3双频电离层误差校正数据由于仪器误差导致校正结果中存在正值,经过类EEMD与改进小波去噪算法滤波处理后的结果均为负值,结果更加合理。且经过滤波处理后的结果较原始数据波动范围更小,但在总体趋势上仍与原始数据保持一致,说明滤波处理既去除了原始数据中的噪声,也保留了可用信息。
最后,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (6)
1.一种高度计测距电离层误差校正滤波处理方法,其特征在于,包括如下步骤:
S1:获取电离层双频校正信号数据,并对其进行预处理,得到原始信号;
S2:对原始信号使用集合经验模态分解算法将其分解为多个本征模态函数IMF和残差信号;
S3:将分解得到的IMF1作为基函数,基于小波去噪的原理首先进行阈值的选取,使用改进的阈值函数计算得到;
S4:基于改进的阈值函数进行小波去噪,低于阈值的信号波动被认为与噪声关联,故将其去除;
S5:通过将去噪后的IMF1与剩余的高阶IMF、残差信号相加,执行信号的重构;
S6:在重构信号中加入均匀分布的白噪声后将其作为原始信号输入,多次迭代过程S2到S5后将得到的所有重构信号取均值,最终得到双频电离层误差校正的滤波结果。
2.如权利要求1所述的高度计测距电离层误差校正滤波处理方法,其特征在于,所述S1中预处理包括剔除缺失值、离群值检索。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310601000.6A CN116299562B (zh) | 2023-05-26 | 2023-05-26 | 一种高度计测距电离层误差校正滤波处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310601000.6A CN116299562B (zh) | 2023-05-26 | 2023-05-26 | 一种高度计测距电离层误差校正滤波处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116299562A true CN116299562A (zh) | 2023-06-23 |
CN116299562B CN116299562B (zh) | 2023-08-04 |
Family
ID=86794597
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310601000.6A Active CN116299562B (zh) | 2023-05-26 | 2023-05-26 | 一种高度计测距电离层误差校正滤波处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116299562B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100179974A1 (en) * | 2009-01-10 | 2010-07-15 | Industrial Technology Research Institute | Signal Processing Method for Hierarchical Empirical Mode Decomposition and Apparatus Therefor |
CN107121705A (zh) * | 2017-04-28 | 2017-09-01 | 中南大学 | 一种基于自动反相校正和峰度值比较的探地雷达回波信号去噪算法 |
CN109145729A (zh) * | 2018-07-13 | 2019-01-04 | 杭州电子科技大学 | 基于改进小波阈值与eemd的肌电信号去噪方法 |
CN110542406A (zh) * | 2019-09-03 | 2019-12-06 | 兰州交通大学 | 基于emd-mpf改进的陀螺仪信号去噪方法 |
CN110688964A (zh) * | 2019-09-30 | 2020-01-14 | 哈尔滨工程大学 | 一种基于稀疏分解的小波阈值与emd联合降噪方法 |
CN112180408A (zh) * | 2020-09-29 | 2021-01-05 | 中山大学 | 一种基于智能终端的多径误差提取方法和相关装置 |
CN113156468A (zh) * | 2021-04-30 | 2021-07-23 | 中国电建集团成都勘测设计研究院有限公司 | 一种大坝gnss多路径改正模型优化方法 |
CN114079503A (zh) * | 2021-10-26 | 2022-02-22 | 广东工业大学 | 一种ceemdan-小波阈值去噪方法、装置、设备及光时域反射仪 |
CN114970602A (zh) * | 2022-04-13 | 2022-08-30 | 浙大宁波理工学院 | 基于改进的经验模态分解与小波阀值函数的信号去噪方法与系统 |
CN115700544A (zh) * | 2022-10-25 | 2023-02-07 | 山东科技大学 | 一种联合经验模态分解及小波软阈值的色谱信号去噪方法 |
-
2023
- 2023-05-26 CN CN202310601000.6A patent/CN116299562B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100179974A1 (en) * | 2009-01-10 | 2010-07-15 | Industrial Technology Research Institute | Signal Processing Method for Hierarchical Empirical Mode Decomposition and Apparatus Therefor |
CN107121705A (zh) * | 2017-04-28 | 2017-09-01 | 中南大学 | 一种基于自动反相校正和峰度值比较的探地雷达回波信号去噪算法 |
CN109145729A (zh) * | 2018-07-13 | 2019-01-04 | 杭州电子科技大学 | 基于改进小波阈值与eemd的肌电信号去噪方法 |
CN110542406A (zh) * | 2019-09-03 | 2019-12-06 | 兰州交通大学 | 基于emd-mpf改进的陀螺仪信号去噪方法 |
CN110688964A (zh) * | 2019-09-30 | 2020-01-14 | 哈尔滨工程大学 | 一种基于稀疏分解的小波阈值与emd联合降噪方法 |
CN112180408A (zh) * | 2020-09-29 | 2021-01-05 | 中山大学 | 一种基于智能终端的多径误差提取方法和相关装置 |
CN113156468A (zh) * | 2021-04-30 | 2021-07-23 | 中国电建集团成都勘测设计研究院有限公司 | 一种大坝gnss多路径改正模型优化方法 |
CN114079503A (zh) * | 2021-10-26 | 2022-02-22 | 广东工业大学 | 一种ceemdan-小波阈值去噪方法、装置、设备及光时域反射仪 |
CN114970602A (zh) * | 2022-04-13 | 2022-08-30 | 浙大宁波理工学院 | 基于改进的经验模态分解与小波阀值函数的信号去噪方法与系统 |
CN115700544A (zh) * | 2022-10-25 | 2023-02-07 | 山东科技大学 | 一种联合经验模态分解及小波软阈值的色谱信号去噪方法 |
Non-Patent Citations (2)
Title |
---|
蔡文航等: "CEEMD联合小波阈值去噪法及其在GPS多路径效应中的应用", 江西科学, vol. 36, no. 1, pages 73 - 78 * |
赵迎;乐友喜;黄健良;王姣;刘陈希;刘兵卿;: "CEEMD与小波变换联合去噪方法研究", 地球物理学进展, no. 06, pages 2870 - 2877 * |
Also Published As
Publication number | Publication date |
---|---|
CN116299562B (zh) | 2023-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111239697B (zh) | 低秩矩阵分解的多维域联合sar宽带干扰抑制方法 | |
CN109269497B (zh) | 基于auv切法向速度模型的多尺度无迹卡尔曼滤波估计方法 | |
CN109871733B (zh) | 一种自适应海杂波信号去噪方法 | |
CN108009122B (zh) | 一种改进的hht方法 | |
CN111982489B (zh) | 选择集成的改进局部特征分解的微弱故障特征提取方法 | |
CN111934711A (zh) | 一种时频混叠跳频信号的参数估计方法 | |
Wang et al. | Frequency down extrapolation with TV norm minimization | |
Lu et al. | A piecewise parametric method based on polynomial phase model to compensate ionospheric phase contamination | |
CN116299562B (zh) | 一种高度计测距电离层误差校正滤波处理方法 | |
Jicheng et al. | Seismic data random noise reduction using a method based on improved complementary ensemble EMD and adaptive interval threshold | |
CN103514600A (zh) | 一种基于稀疏表示的红外目标快速鲁棒跟踪方法 | |
Crawford et al. | Multipath artifact corrections in ultrasonic transmission tomography | |
CN116203634A (zh) | 一种基于低秩约束的鬼波去除方法 | |
CN115630295A (zh) | 双深度神经网络约束的无监督学习分离多震源数据的方法 | |
Wang et al. | Radar emitter signal intra-pulse feature extraction based on empirical mode decomposition | |
CN112149614B (zh) | 一种叠前井震联合智能去噪的方法 | |
CN110703089B (zh) | 一种用于低频振荡Prony分析的小波阈值去噪方法 | |
CN112379415A (zh) | 基于减采样的重构低频数据多尺度全波形反演方法及装置 | |
CN107315713B (zh) | 一种基于非局部相似性的一维信号去噪增强方法 | |
Venkatappareddy et al. | A Novel Thresholding Methodology Using WSI EMD and Adaptive Homomorphic Filter | |
Poczekajło et al. | Synthesis algorithm of low sensitivity non‐separable N‐D finite impulse response filters configured of fast orthogonal 1‐D rotation structures implemented in field programmable gate arrays IC | |
CN117706592A (zh) | 一种基于小波变换的北斗卫星周跳修复方法 | |
Chen et al. | GeoGCN: Geometric dual-domain graph convolution network for point cloud denoising | |
CN115184986B (zh) | 不依赖震源的全局包络互相关全波形反演方法 | |
Lin et al. | Research on optical fiber IMU signal processing based on wavelet algorithm |
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 |