CN108535710A - 一种基于目标环境特征向量的干扰抑制与目标识别方法 - Google Patents
一种基于目标环境特征向量的干扰抑制与目标识别方法 Download PDFInfo
- Publication number
- CN108535710A CN108535710A CN201810182215.8A CN201810182215A CN108535710A CN 108535710 A CN108535710 A CN 108535710A CN 201810182215 A CN201810182215 A CN 201810182215A CN 108535710 A CN108535710 A CN 108535710A
- Authority
- CN
- China
- Prior art keywords
- target
- feature vector
- environment feature
- data
- target environment
- 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.)
- Pending
Links
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于目标环境特征向量的对空目标识别方法。其主要流程是:首先通过对雷达回波检测、点迹凝聚处理,找到目标回波的中心点,提取目标回波数据;通过提取目标环境特征向量对目标数据进行预处理;最后利用目标环境特征向量及预处理后的目标回波数据,提取目标调制谱特征,判断目标是喷气式飞机、螺旋桨飞机或是直升机。本发明所提供的方法主要针对工程应用实现,能够很好的剔除杂波干扰、雷达副瓣干扰的影响,识别正确率达到85%以上,并具有很好的稳健性和适应性。
Description
技术领域
本发明属于雷达目标信号分类识别技术。
背景技术
空中飞机目标主要分为喷气式飞机、螺旋桨飞机和直升机三类。这三类飞机上都存在旋转部件,如直升机的主旋翼和尾旋翼,螺旋桨飞机的螺旋桨叶片和喷气式飞机发动机的压缩叶片,它们在飞机飞行的同时,自身也存在周期性的旋转运动。目标的微动会对雷达回波产生附加的调制,使得目标回波多普勒是时变的,表现在频谱上式存在频谱旁瓣或多普勒频率展宽,折算目标结构部件与目标主体之间相互作用的结果。微动产生的微多普勒特征,即调制谱特征,是目标独特的特征,它反映了目标的几何结构和运动特征,可以用来确定目标的性质,为雷达目标识别提供了新的途径,具有重要的研究价值和意义。
基于微多普勒特征的技术研究是雷达目标识别领域的研究热点之一。如在加拿大Defence R&D Canada公司在研究报告《Micro-Doppler radar signatures forintelligent target recognition》中,对直升机叶片的频域特征进行分析,并且提出了相应的叶片频谱提取方法;在2009年2月的《指挥控制与仿真》中《多普勒频谱及JEM效应在空中目标识别中的应用分析》提出利用相位和幅度综合的方法试图对三类目标进行分类识别;在2006年《现代雷达》中《基于微多普勒特征的空中目标识别》提出利用多普勒频谱能量的不对称性和能量对比实现了固定翼飞机和直升机的分类识别。
与其他文献中提出的方法不同,本发明针对实际雷达的数据特点,通过对目标回波的数据质量进行评估,自适应调整数据预处理和特征提取策略和参数,有效提取目标的调制谱特征,能够有效降低因杂波干扰、副瓣干扰而产生的识别错误,具有更为稳健、抗干扰能力更强的分类识别能力。
发明内容
本发明的目的在于提供一种解决常规相参体制警戒雷达系统的空目标识别方法。通过本发明,能够在警戒雷达系统中实现喷气式飞机、螺旋桨飞机和直升机的分类识别能力,识别正确率达到85%以上。
实现本发明的技术解决方案为:
首先通过对雷达回波检测、点迹凝聚处理,找到目标回波的中心点,提取若干能够用于处理的目标回波数据;通过计算目标环境特征向量,对目标数据进行预处理;最后利用目标环境特征向量计算分割门限,提取目标调制谱特征,判断目标是喷气式飞机、螺旋桨飞机或是直升机。本发明所提供的方法主要针对工程应用实现,具有很好的抗干扰、抗杂波的特点,并且对空目标正确率达到85%以上。
本发明与现有技术相比,其显著优点为:
通过本发明的基于目标环境特征向量的数据预处理方法,能够有效抑制背景杂波、目标副瓣等因素产生的干扰,能准确、有效地提取出目标回波中存在的真实调制谱特征信息,该方法与现有调制谱特征识别处理方法相比,对强杂波干扰、强副瓣干扰、噪声干扰的抑制能力更好,预处理后的目标数据用于特征提取和分类识别的准确率更高;
采用目标环境特征向量进行调制谱特征提取和分类器设计,能够有效对空目标进行分类识别。该方法具有实时性好、检测概率高的特点,更为重要的是,该分类识别方法的适应性更强,能够适应各种强度杂波环境及噪声干扰环境,并且能够适应不同型号雷达,它的提出及工程实现在雷达目标检测与识别领域具有很高的推广应用价值。用本专利的方法,在实际雷达中在各种环境下进行现场测试,其中螺旋桨类246个批次、喷气式飞机900个批次,直升机480个批次,测试结果如下表所示:
表1现场测试结果表
下面结合附图对本发明作进一步详细描述。
附图说明
图1是本发明的工作流程图。
图2是背景杂波剔除处理流程。
图3是背景杂波未剔除原始结果1。
图4是背景杂波剔除效果1。
图5是背景杂波未剔除原始结果2。
图6背景杂波剔除效果2。
图7是副瓣抑制方法处理流程图。
图8是未进行副瓣抑制原始结果1。
图9是副瓣抑制处理效果1。
图10是未进行副瓣抑制原始结果2。
图11是副瓣抑制处理效果2。
具体实施方式
本发明实现方法具体实施步骤为,参见附图1:
(1)根据位置信息提取目标粗略范围根据已经建立的目标航迹信息,得到目标的方位αc和距离dc,利用方位αc和距离dc,从接收机的回波IQ数据截取目标区域的波门数据(二维数组,其大小为M×N,其中M为所取的方位脉冲个数,N为所取的距离单元个数,一般M可取150,N取40)。
(2)目标位置检测。利用MTD检测方法,获得目标所在的距离单元,得到目标回波数据fTarget。
(3)目标距离单元FFT处理。对目标回波数据fTarget进行FFT处理,得到目标回波频谱数据FTarget。
(4)计算目标环境特征向量。利用目标波门数据和FFT处理后的数据FTarget计算目标环境特征向量STB,具体计算方法如下:
a)将目标波门数据投影到脉冲格式方向:
N为所取的距离单元个数,目标区域的波门数据
b)对pPl进行谱变换,得到PPl:
M为所取的方位脉冲个数,v为目标估计的径向运动速度。
c)计算目标环境特征向量STB:
(5)背景杂波剔除。利用目标环境特征向量STB,从0频率开始判断当前频率是否属于杂波成分,逐步剔除目标回波频谱数据FTarget中的背景杂波成分,得到背景剔除后的目标频谱F′Target,具体流程如图2所示,分别是对目标频谱F′Target进行正频率方向和负频率方向的背景杂波剔除处理,具体流程说明如下:
●正频率方向背景杂波剔除处理
a)对目标回波频谱数据FTarget,从0频率开始,向正频率方向开始,判定当前频率点i的幅度FTarget(i)与噪声均值比值,得到频点幅度比kTarget(i);
b)如果kTarget(i)<STB(i),则判定当前频率点i符合噪声特性,则认为当前频率点i即为背景杂波的截止频率,则正频率方向杂波抑制结束;
c)如果kTarget(i)≥STB(i),则计算该频率点的导数F′Target(i);如果导数F′Target(i)<0,则认为当前频率点i属于背景杂波,将i加1,跳转到步骤a);
d)如果导数F′Target(i)≥0,则计算FTarget(i)与FTarget(0)比值,即频率点i幅度与频率点0幅度的比值是否大于KSobel(KSobel一般取0.1),判定是否是副瓣;
e)如果判定是副瓣,则将i加1(i=i+1),跳转回步骤a);否则认为当前频率点i即为背景杂波的截止频率,则正频率方向杂波抑制结束。
●负频率方向背景杂波剔除处理
a)对目标回波频谱数据FTarget,从0频率开始,向负频率方向开始,判定当前频率点i的幅度FTarget(i)与噪声均值比值,得到频点幅度比kTarget(i);
b)如果kTarget(i)<STB(i),则判定当前频率点i符合噪声特性,则认为当前频率点i即为背景杂波的截止频率,则负频率方向杂波抑制结束;
c)如果kTarget(i)≥STB(i),则计算该频率点的导数F′Target(i),如果导数F′Target(i)>0,则认为当前频率点i属于背景杂波,将i加1,跳转到步骤a);
d)如果导数F′Target(i)≤0,则计算FTarget(i)与FTarget(0)比值,即频率点i幅度与频率点0幅度的比值是否大于KSobel(KSobel一般取0.1),判定是否是副瓣;
e)如果判定是副瓣,则将i减1(i=i-1),跳转到步骤a);否则认为当前频率点i即为背景杂波的截止频率,则负频率方向杂波抑制结束。
经过背景杂波剔除前后的对比如图3、图4、图5、图6所示,可见该方法能够非常完整的剔除不同杂波谱宽的背景杂波,而不影响目标频谱;
(6)多普勒速度补偿。根据目标航迹的历史航迹点,用最小二乘法进行直线拟合,估计目标的运动航向和航速信息,进而得到目标的径向速度(目标沿雷达射线方向运动速度)vd,并估计目标自身的多普勒频率fd,其中:
fd=2·vd/λ,λ为雷达发射波长。
利用多普勒频率fd对F′Target进行多普勒速度补偿,将F′Target进行-fd的圆周平移,得到F″Target。
(7)副瓣抑制处理。利用目标环境特征向量STB对F″Target进行副瓣抑制,得到预处理结果具体处理流程如图7所示,具体流程说明如下:
a)估计当前频率点的副瓣幅度,初始i=-N/2(N是FFT后的数据长度),对副瓣抑制处理后数据F″Target,利用当前频率点i附近频率点的幅度估计当前频率点i的副瓣幅度FSobel:
b)比较估计的频率点i副瓣幅度FSobel(i)与实际的FTarget(i)大小,如果FSobel(i)≤FTarget(i),则跳转到步骤c;否则,根据对当前频率点i进行副瓣衰减处理,衰减方法为:
其中,F″Target为目标频谱经过背景杂波剔除、多普勒速度补偿后的结果。
c)i=i+1,跳转到步骤a,直到完成所有频率点副瓣抑制。
经过副瓣抑制前后的对比如图8、图9、图10、图11所示,可见目标频谱的副瓣都得到了很好的抑制,对后面分割处理不再会产生影响;
(8)门限分割处理。对副瓣抑制后的数据,需要通过设定一定的门限,判定具体哪些是噪声点和调制谱,采用数据质量评估结果动态调整分割门限,以能对各种条件下有效的提取目标调制谱。主要处理流程如下:
a)首先根据目标环境特征向量STB,选择门限计算策略,预估基础门限;
分别是用目标回波距离单元数据估计出的噪声均值和方差;分别是用背景噪声估计出的噪声均值和方差,n1、n2为常数。
b)根据目标环境特征向量STB的模,对门限再次进行调整。
k1、k2为常数。
c)根据目标的信噪比(SNR)大小,对门限进行调整;
h1、h2为常数。
(9)调制特征提取。经过门限分割后,剩余的频谱主要是目标机身产生的主频和由旋转部件运动产生的调制谱,针对飞机的调制谱特征分布,主要提取7种特征,主要包括:调制谱个数、最小谱间隔、最大谱间隔、平均谱间隔、邻域内积累谱数量、对称调制谱个数,去对称性最小谱间隔。
(10)分类识别。通过所提取的调制谱特征,设计对3类目标适用的分类器,主要分类准侧如下表所示:
表2不同类型目标特征分布特点
Claims (2)
1.一种基于目标环境特征向量的干扰抑制与目标识别方法,其特征在于:通过检测、凝聚方法,找到目标回波的中心点,提取能够用于处理的目标回波数据;对目标数据进行质量评估,计算目标环境特征向量,根据目标环境特征向量,对目标数据进行预处理,包括背景杂波剔除和副瓣抑制;依据目标环境特征向量,对预处理后的目标数据进行调制谱特征和分类处理,判断目标类型。
2.根据权利要求1所述的一种基于目标环境特征向量的干扰抑制与目标识别方法,其特征在于目标数据预处理方法:计算目标环境特征向量,利用该向量构建变化矩阵对目标数据进行预处理,剔除背景杂波、抑制信号副瓣。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810182215.8A CN108535710A (zh) | 2018-03-06 | 2018-03-06 | 一种基于目标环境特征向量的干扰抑制与目标识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810182215.8A CN108535710A (zh) | 2018-03-06 | 2018-03-06 | 一种基于目标环境特征向量的干扰抑制与目标识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108535710A true CN108535710A (zh) | 2018-09-14 |
Family
ID=63486753
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810182215.8A Pending CN108535710A (zh) | 2018-03-06 | 2018-03-06 | 一种基于目标环境特征向量的干扰抑制与目标识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108535710A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109633584A (zh) * | 2018-10-31 | 2019-04-16 | 上海无线电设备研究所 | 雷达lprf工作条件下的直升机目标识别方法 |
CN109975780A (zh) * | 2019-04-17 | 2019-07-05 | 西安电子工程研究所 | 基于脉冲多普勒雷达时域回波的直升机型号识别算法 |
CN113534093A (zh) * | 2021-08-13 | 2021-10-22 | 北京环境特性研究所 | 飞机目标的螺旋桨叶片数量反演方法及目标识别方法 |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB8823753D0 (en) * | 1987-10-28 | 1998-05-27 | Licentia Gmbh | Method of type classification of a target |
CN102043148A (zh) * | 2011-01-10 | 2011-05-04 | 西安电子科技大学 | 基于广义匹配滤波的地面运动目标分类方法 |
CN102176000A (zh) * | 2011-01-27 | 2011-09-07 | 电子科技大学 | 一种船用雷达海杂波的抑制方法 |
CN102540163A (zh) * | 2011-12-12 | 2012-07-04 | 中国船舶重工集团公司第七二四研究所 | 一种基于常规相参雷达的直升机分类识别方法 |
CN102540162A (zh) * | 2011-12-12 | 2012-07-04 | 中国船舶重工集团公司第七二四研究所 | 基于海杂波的低空电磁波传播特性估计方法 |
CN102565784A (zh) * | 2012-01-04 | 2012-07-11 | 西安电子科技大学 | 基于vsar系统的动目标重定位与速度解模糊方法 |
CN102621534A (zh) * | 2011-01-26 | 2012-08-01 | 中国科学院电子学研究所 | 基于噪声失效奇异值分解的地基雷达杂波抑制方法 |
CN102930276A (zh) * | 2012-11-05 | 2013-02-13 | 中国船舶重工集团公司第七二四研究所 | 基于线性回归和聚类分析理论的高速机动目标跟踪方法 |
CN104133211A (zh) * | 2014-07-07 | 2014-11-05 | 中国船舶重工集团公司第七二四研究所 | 一种多普勒频率变换雷达目标分类识别方法 |
CN104198993A (zh) * | 2014-07-29 | 2014-12-10 | 北京理工大学 | 适合参数估计的认知雷达最优波形设计方法 |
CN104515975A (zh) * | 2014-12-12 | 2015-04-15 | 中国电子科技集团公司电子科学研究院 | 一种面向杂波抑制的相干mimo雷达波形设计方法 |
CN104931949A (zh) * | 2015-06-15 | 2015-09-23 | 中国民航大学 | 雷达扫描模式下风轮机杂波背景下飞机目标检测方法 |
CN105242254A (zh) * | 2015-10-22 | 2016-01-13 | 中国船舶重工集团公司第七二四研究所 | 一种基于数据质量评估的对空目标识别方法 |
CN105388465A (zh) * | 2015-12-17 | 2016-03-09 | 西安电子科技大学 | 基于海浪谱模型的海杂波模拟方法 |
CN106772263A (zh) * | 2016-12-30 | 2017-05-31 | 无锡市雷华科技有限公司 | 对地监视雷达 |
WO2017188905A1 (en) * | 2016-04-29 | 2017-11-02 | Aselsan Elektronik Sanayi Ve Ticaret Anonim Sirketi | A method for motion classification using a pulsed radar system |
-
2018
- 2018-03-06 CN CN201810182215.8A patent/CN108535710A/zh active Pending
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB8823753D0 (en) * | 1987-10-28 | 1998-05-27 | Licentia Gmbh | Method of type classification of a target |
CN102043148A (zh) * | 2011-01-10 | 2011-05-04 | 西安电子科技大学 | 基于广义匹配滤波的地面运动目标分类方法 |
CN102621534A (zh) * | 2011-01-26 | 2012-08-01 | 中国科学院电子学研究所 | 基于噪声失效奇异值分解的地基雷达杂波抑制方法 |
CN102176000A (zh) * | 2011-01-27 | 2011-09-07 | 电子科技大学 | 一种船用雷达海杂波的抑制方法 |
CN102540163A (zh) * | 2011-12-12 | 2012-07-04 | 中国船舶重工集团公司第七二四研究所 | 一种基于常规相参雷达的直升机分类识别方法 |
CN102540162A (zh) * | 2011-12-12 | 2012-07-04 | 中国船舶重工集团公司第七二四研究所 | 基于海杂波的低空电磁波传播特性估计方法 |
CN102565784A (zh) * | 2012-01-04 | 2012-07-11 | 西安电子科技大学 | 基于vsar系统的动目标重定位与速度解模糊方法 |
CN102930276A (zh) * | 2012-11-05 | 2013-02-13 | 中国船舶重工集团公司第七二四研究所 | 基于线性回归和聚类分析理论的高速机动目标跟踪方法 |
CN104133211A (zh) * | 2014-07-07 | 2014-11-05 | 中国船舶重工集团公司第七二四研究所 | 一种多普勒频率变换雷达目标分类识别方法 |
CN104198993A (zh) * | 2014-07-29 | 2014-12-10 | 北京理工大学 | 适合参数估计的认知雷达最优波形设计方法 |
CN104515975A (zh) * | 2014-12-12 | 2015-04-15 | 中国电子科技集团公司电子科学研究院 | 一种面向杂波抑制的相干mimo雷达波形设计方法 |
CN104931949A (zh) * | 2015-06-15 | 2015-09-23 | 中国民航大学 | 雷达扫描模式下风轮机杂波背景下飞机目标检测方法 |
CN105242254A (zh) * | 2015-10-22 | 2016-01-13 | 中国船舶重工集团公司第七二四研究所 | 一种基于数据质量评估的对空目标识别方法 |
CN105388465A (zh) * | 2015-12-17 | 2016-03-09 | 西安电子科技大学 | 基于海浪谱模型的海杂波模拟方法 |
WO2017188905A1 (en) * | 2016-04-29 | 2017-11-02 | Aselsan Elektronik Sanayi Ve Ticaret Anonim Sirketi | A method for motion classification using a pulsed radar system |
CN106772263A (zh) * | 2016-12-30 | 2017-05-31 | 无锡市雷华科技有限公司 | 对地监视雷达 |
Non-Patent Citations (2)
Title |
---|
徐利刚,等: "一种基于调制谱的空中目标识别方法", 《指挥控制与仿真》 * |
管志强: "红外搜索系统中目标探测与识别技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109633584A (zh) * | 2018-10-31 | 2019-04-16 | 上海无线电设备研究所 | 雷达lprf工作条件下的直升机目标识别方法 |
CN109633584B (zh) * | 2018-10-31 | 2020-11-03 | 上海无线电设备研究所 | 雷达lprf工作条件下的直升机目标识别方法 |
CN109975780A (zh) * | 2019-04-17 | 2019-07-05 | 西安电子工程研究所 | 基于脉冲多普勒雷达时域回波的直升机型号识别算法 |
CN113534093A (zh) * | 2021-08-13 | 2021-10-22 | 北京环境特性研究所 | 飞机目标的螺旋桨叶片数量反演方法及目标识别方法 |
CN113534093B (zh) * | 2021-08-13 | 2023-06-27 | 北京环境特性研究所 | 飞机目标的螺旋桨叶片数量反演方法及目标识别方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108490410B (zh) | 一种两坐标雷达对海目标联合检测跟踪方法 | |
US7675458B2 (en) | Dual beam radar system | |
CN105242254A (zh) | 一种基于数据质量评估的对空目标识别方法 | |
CN109975780B (zh) | 基于脉冲多普勒雷达时域回波的直升机型号识别算法 | |
US20070024494A1 (en) | Classification system for radar and sonar applications | |
CN104267379A (zh) | 一种基于波形设计的主被动雷达协同抗干扰方法 | |
CN104133211B (zh) | 一种多普勒频率变换雷达目标分类识别方法 | |
CN110231615B (zh) | 基于海面目标空间约束的距离测量异常判别方法 | |
CN108535710A (zh) | 一种基于目标环境特征向量的干扰抑制与目标识别方法 | |
CN105116396B (zh) | 一种连续波雷达多普勒回波检测方法 | |
WO2016194036A1 (ja) | レーダ信号処理装置 | |
CN104215938B (zh) | 一种对地形反弹干扰的极化检测识别方法和装置 | |
EP3617740B1 (en) | Target detection in rainfall and snowfall conditions using mmwave radar | |
CN105785340A (zh) | 单脉冲雷达对主瓣内目标和诱饵干扰的测角与识别方法 | |
CN105842688A (zh) | 单脉冲雷达对空目标快速捕获方法 | |
CN107632291A (zh) | 一种基于极化雷达抗角反射器干扰的方法 | |
Molchanov et al. | Aerial target classification by micro-Doppler signatures and bicoherence-based features | |
CN104199002B (zh) | 抑制风力发电旋叶虚警目标的雷达信息处理方法 | |
CN106168661B (zh) | 随机捷变极化有源假目标的鉴别方法 | |
Xu et al. | Long-time coherent integration for radar target detection base on Radon-Fourier transform | |
CN111123269B (zh) | 用于无人机避障雷达的地面杂波抑制方法、模块及装置 | |
CN109683158A (zh) | 基于干涉仪信息融合的雷达地面静止防空单元检测方法 | |
CN111337920B (zh) | 一种防止云雾干扰的弹载雷达对地探测方法及装置 | |
McDonald et al. | Track-before-detect using swerling 0, 1, and 3 target models for small manoeuvring maritime targets | |
Tikkinen et al. | Joint utilization of incoherently and coherently integrated radar signal in helicopter categorization |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180914 |