CN106935079B - 基于软件无线电平台的ads-b信号采集与可靠性验证方法 - Google Patents
基于软件无线电平台的ads-b信号采集与可靠性验证方法 Download PDFInfo
- Publication number
- CN106935079B CN106935079B CN201710107967.3A CN201710107967A CN106935079B CN 106935079 B CN106935079 B CN 106935079B CN 201710107967 A CN201710107967 A CN 201710107967A CN 106935079 B CN106935079 B CN 106935079B
- Authority
- CN
- China
- Prior art keywords
- ads
- signal
- message
- message information
- reliability
- 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
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000012795 verification Methods 0.000 title claims abstract description 9
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 238000001514 detection method Methods 0.000 claims abstract description 7
- 238000010219 correlation analysis Methods 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 6
- 230000009191 jumping Effects 0.000 claims description 4
- 102100026882 Alpha-synuclein Human genes 0.000 claims description 3
- 108090000185 alpha-Synuclein Proteins 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 4
- 238000012544 monitoring process Methods 0.000 abstract description 13
- 238000012545 processing Methods 0.000 abstract description 6
- 238000004891 communication Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0073—Surveillance aids
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Traffic Control Systems (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种基于软件无线电平台的ADS‑B信号采集与可靠性验证方法,属于信号检测与处理及安全监测领域,本发明的目的是解决现有民用航空器ADS‑B监视的准确性和可靠性不足的问题。本发明方法包括以下步骤:一、驱动软件无线电平台采集民航飞机的ADS‑B信号;二、对解调后的ADS‑B信号,设计一种基于的小波滤波器,滤除噪声,提取信号特征;三、设计一种基于相关性分析的报头检测算法,对ADS‑B数字信号解码,得到民航飞机的飞行参数信息;四、基于报文信息,对解析的民航飞机ADS‑B信息可靠性加以验证。本发明用于民航飞机ADS‑B监视及运行安全保障。
Description
技术领域
本发明涉及一种基于软件无线电平台的ADS-B信号采集与可靠性验证方法,属于民用航空器ADS-B监视及运行安全保障领域。
背景技术
在现代航空管理领域有许多亟待解决的问题,首要的问题是广泛使用的雷达系统无法工作于非陆地区域、雷达盲区以及机场跑道、滑行道上,使得像沙漠、海洋的监视管理工作很难进行;同样在大型机场的场面跑道、滑行道的雷达监视管理工作也很难进行。其次是航空通信方面,传统ATC系统中90%使用的都是语音通信,这种通信模式会导致飞机与基站之间的通信不足,存在安全隐患。在导航方面,如今的固定航线、航线分配模式难以实现最短距离、最佳航线行驶,这样不仅浪费时间、燃油,还造成严重污染。
为了应对日益增长的航空需求和空中交通管制相对落后的矛盾,国际民航组织提出一种新航行系统下的广播式自动相关监视技术即ADS-B。ADS-B是一种基于全球卫星定位系统和空-空、地-空数据链通信的航空器运行监视技术。ADS-B技术除可以有效的解决上述问题外,还将是实施自由飞行的奠基石,欧洲虽然并不提倡前景遥远、偏于理想、概念模糊的“自由飞行”,但也提出了“自由航路”概念,而ADS-B也将是促进实施自由航路的可行手段之一。
ADS-B和传统雷达监视系统相比可以获得更高精度的目标位置和高度信息,突破其无法覆盖荒漠远洋等地区、数据精度受限的问题。通过为飞行员提供更加快捷的飞行数据传输服务,使用了ADS-B可以优化航路设置,减小间隔标准,使空域得到更充分的利用。此外,ADS-B还可以为航空器提供更加全面天气、地形、周边交通情况的飞行信息,使机组更加了解航空器周边的情况,提高情境意识,从而实现空、天、地一体化信息共享和协同监视。另外,从经济角度来看,相较于传统的航管雷达需要建立地面基站,定期维护保养耗资巨大不同,ADS-B的建设保养费用都远远低于传统雷达,具有推广意义。因此,ADS-B被认为是空中交通管制的未来。
当前市场流行ADS-B设备是一个“黑箱”,用户只能得到解析后的民航飞机ADS-B航迹数据(高度、经纬度、速度等),且无法判断解析的数据是否可靠。因此空管运行单位只能将其作为辅助监视手段。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种基于软件无线电平台采集民航飞机ADS-B信号并进行可靠性验证的方法,用于解决现有技术中只能得到解析后的民航飞机ADS-B航迹数据且无法判断解析的数据是否可靠的技术问题。
技术方案:为实现上述目的,本发明采用的技术方案为:步骤一、驱动软件无线电平台采集民航飞机的ADS-B信号;
软件无线电(USRP)代替传统硬件无线电平台是现阶段信号处理的主流。软件无线电采用固定不变的硬件平台,通过软件升级来实现灵活多变的通信体制和通信功能的无线电系统。其特点是通用化、标准化、模块化,以及对信号波形的广泛适应性。软件无线电的核心是其驻留在FPGA内部的可升级、可重构的功能软件,以适应不同的技术标准、接口协议和信号波形。
在软件无线电平台内设置当前ADS-B基站的地理信息坐标、天线接收频率接收范围和拟采集信号的功率阈值以及信号采样率等参数,完成ADS-B信号采集环境配置。驱动USRP平台采集接收范围内ADS-B信号,存储待滤波和解码处理。
步骤二,基于小波滤波器的ADS-B信号带通滤波器设计
步骤2.1、按照(1)式对USRP采集经I/Q解调后的ADS-B信号x(n),对其进行K阶DB小波分解和重构:
其中,xk(n)为采用k阶小波系数重构的信号;n为x(n)的采样点;
步骤2.2、采用(2)式的相关算子计算xk(n)的频率:
其中,Ts为x(n)的采样周期,∠为复函数的相角,N为x(n)的采样点总数,上标*表示对偶;
步骤2.3、为抑制信道外干扰,设置频率阈值fth,按照(3)式设计带通滤波器进行滤波:
同时,为防止信道内脉冲干扰,按照(4)式计算xk(n)与标准1090MHz脉冲信号xs(n)之间的相关性ek
计算得到的ek需满足ek>eth,eth为预设的相关性阈值;
因此,得到综合滤波后的信号如(5)式所示
步骤三、设计一种基于相关性分析的报头检测算法解码ADS-B数字信号
对滤波处理后的ADS-B信号y(n)取模以得到信号包络。系统自动生成一个时长8μs固定格式的ADS-B报头与信号包络进行匹配处理。匹配得到的输出波形的起始处检测到一个相关峰值,而其临近区间不会有较大的峰值输出,即找到一个ADS-B信号起始位置。记录该起始位置后长度为112比特的ADS-B数字信号(由数据0或数据1组成)形成一条报文。
在现有技术中,按照ADS-B协议标准将数字信号解析为航班飞行参数,步骤为:
步骤c1、对每条报文预处理,内容包括:(1)判断该条报文是否为非空;(2)判断该条报文的存储格式是否为88位;(3)判断该条报文是否为ADS-B消息。对同时满足上述3条要求的报文进行下一步处理。
步骤c2、对预处理后的报文按相同的注册号即ICAO地址进行分组,并将相同注册号的报文存入同一分组中,即将同一飞机的报文存入同一分组中。
步骤c3、对同一分组中的报文进行如下操作:对当前一条报文与其之前的一条报文二者的第54位进行比较,判断这两条报文是否匹配,若匹配,则采用全局CPR算法求解出与当前报文相对应的飞行参数;若不匹配,采用局部CPR算法求解出与当前报文相对应的飞行参数。
步骤c4、依次对不同注册号的飞机对应的报文进行以上的解码操作,求解出不同飞机的飞机报文解析结果。
步骤四、基于报文的解析结果,计算如下4种精确性指标,对报文信息可靠性加以验证。
A、利用位置导航精度类别比率μNACp表示报文信息中位置信息的精度,位置导航精度类别比率μNACp按照(6)式计算:
NACp为位置导航精度类别指标Navigation Accuracy Category-Position的简称,报文解码后该指标是一个0到无穷大的数值,所以这里设定NACp>5表示导航精度高;
NoAT(H)为H时间内,所选定的一架飞机所监视到的所有轨迹点数目;H的选择范围在2-24小时之间;
NoAT(h)为第h秒时所选定的一架飞机所监视到的轨迹点;
B、利用速度导航精度类别比率μNACv表示报文信息中速度信息的精度,速度导航精度类别比率μNACv按照(7)式计算:
NACv为速度导航精度类别Navigation Accuracy Category-Velocity的简称,报文解码后该指标是一个0到无穷大的数值,所以这里设定NACv>5表示导航精度高;
C、按照(8)式计算报文信息中的航迹漏点率μmp
NoMP(H)为H时间内,所选定的一架飞机的所有漏点数目。对于选定任意一架航空器,连续5秒内未收到航迹点,则统计为报文信息中的1个航迹漏点。
D、按照(9)式计算ADS-B信息中的航迹跳点率μjp:
NoJP(H)为H时间内,选定某一架飞机所有跳点数目。考虑到各种因素的噪声干扰以及在数据解析过程存在的解析误差,导致实际航迹数据与真实航迹数据存在一定的偏差,即为航迹数据的位置离散度。计算航迹数据的位置离散度的时候,首先提取H时间段的报文信息,用卡尔曼滤波算法对报文信息中的位置信息进行滤波处理。计算滤波前后两个航迹点之间的误差,计算误差的均值和方差,当误差的均值和方差大于阈值,则认为此位置点为跳点。
依据USRP设备采集H时间内实时ADS-B航迹,对于某一架航班,利用上述4个精确性指标进行可靠性指标Ire判别,得到采用经验公式如下
Ire=w1(1-μNACp+1-μNACv)+w2μmp+w3μjp (10)
其中w1、w2、w3均为权重,满足w2>w1>w3且w1+w2+w3=1;
根据研究多组过去24小时内ADS-B解析的航迹信息(经度、纬度、速度)对管制的影响作用得到的经验可知,其中航迹漏点率最重要,直接影响航迹是否可用,其次是位置导航精度类别比率和速度导航精度类别比率,最后是航迹跳点率,因为航迹跳点率可能由高空风等因素影响。
当可靠性指标Ire处于如下阈值范围内,如表1所示,认为其在H时间段内的报文信息可靠性通过,否则,认为在H时间段内报文信息不可靠。
表1
有益效果:
本发明提供的基于软件无线电平台采集民航飞机ADS-B信号,并结合小波分析的带通滤波器设计和相关性分析的报头检测算法对ADS-B信号进行解码,实现对民航飞机飞行状态监视的同时,依据报文所包含的准确性信息特征,对ADS-B监视的可靠性能够进行验证。具体的:
1)本发明所提出的基于小波分析的滤波器能够抑制ADS-B接收的信道外和信道内的干扰。
2)本发明所提出的基于软件无线电的ADS-B采集系统能够获取ADS-B从采集-解码-解析-显示-验证多个模态下的内容,进而能够实现抑制干扰。
3)本发明所提出的验证方法能够验证民航飞机ADS-B航迹的可靠性,将本方法运用于空管运行单位,经过实践验证,当报文的可靠程度为可靠时,能够准确把握民航飞机位置,为空管做出正确决策提供保障,当报文的可靠程度为较可靠时,也能为空管做出正确决策给予一定的参考作用。
附图说明
图1为本发明的流程图;
图2为ADS-B信号采集图;
图3为ADS-B信号经过小波滤波器滤波后的频谱图;
图4为报头相关性检测后提取的ADS-B数字信号结果图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
本实施方式包括以下步骤:
步骤一、驱动USRP B210软件无线电平台,采集民航飞机的ADS-B信号,对软件无线电平台设置当前ADS-B基站的地理信息坐标;设置天线接收频率接收范围和拟采集信号的功率阈值以及信号采样率等参数,完成ADS-B信号采集环境配置,实现ADS-B信号的采集,如图2所示。
步骤二、基于的小波滤波器的ADS-B信号带通滤波器设计,采集经I/Q解调后的信号ADS-B信号,进行3阶DB小波分解和重构,如图3所示。
步骤三、设计一种基于相关性分析的报头检测算法解码ADS-B数字信号,对滤波处理后的ADS-B数字复基带信号取模以得到信号包络。然后,根据包络位置信息进行同步和码字的判决,即在这一串数字信号序列中找到相对应的ADS信号波形起始位置。结果如图4所示。
步骤四、基于解析报文的精确性信息,对ADS-B信息可靠性加以验证。依据4项精确性指标,对ADS-B信号解析出的数据进行精确性分析,根据大量数据和专家经验,考虑到航迹漏点率最为重要,接近一半,故占考虑因素的45%,取w2=0.45,导航和速度精度对ADS-B航迹能用的参考价值占三分之一左右,故考虑为30%,取w1=0.3,跳点率可能是风、误差的原因,因此占25%,取w3=0.25,从而得出数据的可靠性分析结果。结果如表2所示。
表2
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (2)
1.基于软件无线电平台的ADS-B信号采集与可靠性验证方法,其特征在于:包括顺序执行的以下步骤:
步骤一、利用软件无线电平台,设置天线工作频率和增益,采集并存储民航飞机的ADS-B信号,并进行I/Q解调;
步骤二、对解调后的ADS-B信号,采用基于小波分析的滤波器,依据信号频率特征,滤除高低频噪声;
步骤三、根据ADS-B数字信号编码协议,设计基于相关性分析的报头检测算法解码ADS-B信号所对应的报文信息,得到民航飞机的飞行参数信息;
步骤四、利用得到的飞行参数信息,计算报文信息中的位置导航精度类别比率μNACp、速度导航精度类别比率μNACv、航迹漏点率μmp、以及航迹跳点率μjp形成4个精确性指标,利用上述精确性指标组成如下经验公式计算可靠性指标Ire
Ire=w1(1-μNACp+1-μNACv)+w2μmp+w3μjp
其中:w1、w2、w3均为权重,满足w2>w1>w3且w1+w2+w3=1;
当可靠性指标Ire小于5%时认为报文信息可靠,当可靠性指标Ire为大于等于5%且小于等于10%时认为报文信息较可靠,当可靠性指标Ire大于10%时认为报文信息不可靠;将可靠性指标Ire和对应的飞行参数信息均发送给管制员供管制员参考;
所述位置导航精度类别比率μNACp的计算公式如下:
NACp为位置导航精度类别的简称,报文解码后该指标是一个0到无穷大的数值;NoAT(H)为H时间内,所选定的一架飞机所监视到的所有轨迹点数目;H的选择范围在2-24小时之间;NoAT(h)为第h秒时所选定的一架飞机所监视到的轨迹点;
所述速度导航精度类别比率μNACv的计算公式如下:
NACv为速度导航精度类别的简称,报文解码后该指标是一个0到无穷大的数值;
所述航迹漏点率μmp的计算公式如下:
NoMP(H)为H时间内,所选定的一架飞机的所有漏点数目;
所述航迹跳点率μjp的计算公式如下:
NoJP(H)为H时间内,选定某一架飞机所有跳点数目;
所述步骤二包括以下过程:
首先,将解调后的ADS-B信号x(n),按照(1)式进行K阶DB小波分解和重构:
其中,xk(n)为采用k阶小波系数重构的信号;n为x(n)的采样点;
然后,按照(2)式计算xk(n)的频率:
其中,Ts为x(n)的采样周期,∠为方括号中复函数的相角,N为x(n)的采样点总数,上标*表示对偶;
最后,按照(5)式的带通滤波器进行滤波,获得滤波信号:
2.根据权利要求1所述的基于软件无线电平台的ADS-B信号采集与可靠性验证方法,其特征在于:步骤三中的具体过程如下:
步骤3.1、利用报头的相关分析方法,确定ADS-B信号的位置并提取ADS-B数字信号形成报文;
步骤3.2、按照ADS-B协议标准解码ADS-B信号中的报文为民航飞机飞行参数信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710107967.3A CN106935079B (zh) | 2017-02-27 | 2017-02-27 | 基于软件无线电平台的ads-b信号采集与可靠性验证方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710107967.3A CN106935079B (zh) | 2017-02-27 | 2017-02-27 | 基于软件无线电平台的ads-b信号采集与可靠性验证方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106935079A CN106935079A (zh) | 2017-07-07 |
CN106935079B true CN106935079B (zh) | 2020-04-24 |
Family
ID=59424176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710107967.3A Active CN106935079B (zh) | 2017-02-27 | 2017-02-27 | 基于软件无线电平台的ads-b信号采集与可靠性验证方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106935079B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10650689B2 (en) | 2016-11-01 | 2020-05-12 | The Mitre Corporation | Waveform authentication system and method |
CN107369339A (zh) * | 2017-08-29 | 2017-11-21 | 中国民用航空总局第二研究所 | 基于cpci架构的ads‑b下行数据链分析验证系统 |
CN109327254A (zh) * | 2018-07-25 | 2019-02-12 | 南京航空航天大学 | 基于软件无线电的ads-b收发一体机及方法 |
CN114049796B (zh) * | 2021-11-09 | 2022-08-16 | 中国电子科技集团公司第二十八研究所 | 一种定向传输ads-b广播信号的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101110165A (zh) * | 2007-08-15 | 2008-01-23 | 民航数据通信有限责任公司 | 一种ads-b监视技术的数据评估系统 |
EP1884908A2 (en) * | 2006-07-25 | 2008-02-06 | Era Systems Corporation | Method and system to extend ADS performance metrics |
CN101575970A (zh) * | 2008-05-09 | 2009-11-11 | 高岩 | 一种随钻岩性和油藏特性识别方法 |
CN106027201A (zh) * | 2016-05-06 | 2016-10-12 | 电子科技大学 | 基于相关的星载ads-b报头检测方法 |
-
2017
- 2017-02-27 CN CN201710107967.3A patent/CN106935079B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1884908A2 (en) * | 2006-07-25 | 2008-02-06 | Era Systems Corporation | Method and system to extend ADS performance metrics |
CN101110165A (zh) * | 2007-08-15 | 2008-01-23 | 民航数据通信有限责任公司 | 一种ads-b监视技术的数据评估系统 |
CN101575970A (zh) * | 2008-05-09 | 2009-11-11 | 高岩 | 一种随钻岩性和油藏特性识别方法 |
CN106027201A (zh) * | 2016-05-06 | 2016-10-12 | 电子科技大学 | 基于相关的星载ads-b报头检测方法 |
Non-Patent Citations (4)
Title |
---|
ADS-B监视数据质量分析研究;钟建华;《西安航空学院学报》;20130531;第31卷(第3期);第72-75页 * |
Study of ADS-B Data Evaluation;ZHANG Jun 等;《Chinese Journal of Aeronautics 》;20110831;第461-466页 * |
一种ADS-B报头互相关检测方法;张涛 等;《电讯技术》;20160229;第56卷(第2期);第156-160页 * |
基于ADS-B数据监视性能评估技术研究;唐鹏;《中国优秀硕士学位论文全文数据库》;20160215;第15-23页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106935079A (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106935079B (zh) | 基于软件无线电平台的ads-b信号采集与可靠性验证方法 | |
US20190304313A1 (en) | Estimating aircraft operations at airports using transponder data | |
CN103199944B (zh) | 广播式自动相关监视信号检测方法及装置 | |
CN102854884B (zh) | 向飞行管理系统提供飞行器轨迹的气象信息的方法 | |
CN105501464B (zh) | 试飞测试航空机载无线电通讯与导航系统性能的方法 | |
CN105867272A (zh) | 一种通用航空飞行器及无人机的安全监控方法 | |
US6812890B2 (en) | Voice recognition landing fee billing system | |
EP2541529A2 (en) | Meteorological data selection along an aircraft trajectory | |
EP2541528A2 (en) | Meteorological modeling along an aircraft trajectory | |
CN108318868A (zh) | 基于ads-b数据源的雷达战术性能测试评估方法 | |
CN114295137B (zh) | 一种基于fpga的塔康信号识别方法及系统 | |
CN107944701B (zh) | 一种飞机着陆过程中冲出跑道风险的检测方法及装置 | |
CN103616665A (zh) | 一种利用通用照射源的无人机无源定位系统 | |
CN110060515B (zh) | 一种基于gnss的航空器监控系统及方法 | |
CN107329151B (zh) | 一种电力巡检无人机的gps欺骗检测方法 | |
CN108134649B (zh) | 一种星载ads-b信号突发检测方法及系统 | |
Verbraak et al. | Large-scale ADS-B data and signal quality analysis | |
CN104502923A (zh) | 一种机场gnss监测接收系统信号质量监测的方法 | |
CN103136033B (zh) | 基于ads-b测试信标设备的3d航迹模拟系统及方法 | |
CN102331579A (zh) | 机载导航定位系统以及应用该系统的卫星导航定位方法 | |
CN113031018A (zh) | 一种基于无人机的北斗卫星导航系统性能测试装置 | |
CN103546468A (zh) | 基于北斗卫星的无人机飞行诸元数据的传输方法及系统 | |
CN108737030B (zh) | 一种基于星载场景的ads-b信号低灵敏度接收方法 | |
Varga et al. | ADS-B based real-time air traffic monitoring system | |
CN104394107B (zh) | 一种基于二次拐点法的1090es信号解同频交叠方法 |
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 |