CN111796249A - 一种基于墙体前表面能量对消的墙体参数估计方法 - Google Patents

一种基于墙体前表面能量对消的墙体参数估计方法 Download PDF

Info

Publication number
CN111796249A
CN111796249A CN202010126142.8A CN202010126142A CN111796249A CN 111796249 A CN111796249 A CN 111796249A CN 202010126142 A CN202010126142 A CN 202010126142A CN 111796249 A CN111796249 A CN 111796249A
Authority
CN
China
Prior art keywords
wall
echo
front surface
processing
wall body
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
Application number
CN202010126142.8A
Other languages
English (en)
Other versions
CN111796249B (zh
Inventor
丁泽刚
龙腾
向寅
王岩
赵祎昆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202010126142.8A priority Critical patent/CN111796249B/zh
Publication of CN111796249A publication Critical patent/CN111796249A/zh
Application granted granted Critical
Publication of CN111796249B publication Critical patent/CN111796249B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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/418Theoretical aspects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/887Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
    • G01S13/888Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于墙体前表面能量对消的墙体参数估计方法。基于接收到的回波信号的构成,对该信号进行分别处理,将该信号的直达波、墙体前表面回波、墙后目标回波消除,获得墙体后表面回波信号和噪声混杂的信号。之后对该信号进行分析,获得能量值最高的时刻,即为计算墙体参数估计公式中,墙体后表面回波的时延精估计值。之后利用牛顿迭代法求解墙体参数,这种方法可以快速并且精准的估计出墙体参数,无需遍历搜索,且即使在墙后表面能量低的情况下,也不会失效。

Description

一种基于墙体前表面能量对消的墙体参数估计方法
技术领域
本发明涉及一种基于墙体前表面能量对消的墙体参数估计方法,属于穿墙 雷达技术领域。
背景技术
穿墙雷达成像是近些年快速发展起来的一种对墙后目标进行微波成像的技 术,他能够对墙体后的目标进行探测、成像、定位、跟踪与识别等,因此在城 市巷战、反恐安保、灾害及人质救援有着广泛的应用。
在电磁波信号穿透墙壁的过程中,电磁波会经历折射、反射和多径传播及 传播速度改变等过程,这意味着目标的回波信号将会存在显著的衰减、失真、 额外的时延以及相位的突变等。这些效应都会对后续的信号处理过程和对目标 区域进行准确的成像带来极大的挑战,需要对获得的信号进行细致的考虑和仔 细的处理。由于墙体的材料、建筑机构、厚度等特性往往各不相同,很难用一 个统一的模型来描述他们对电磁波传播的影响,因此,对于不同电磁属性的墙 体需要更为细致的研究。在一般意义上的透障成像中,墙体信息都被假定为已 知的。但在真实情况中,墙体的主要属性,例如介电常数和厚度,往往都是未 知的,而这会给成像处理带来极大的困扰,造成图像质量的下降。为此,需要 获得墙体的准确参数。
现有的墙体未知参数估计有以下几种方案:一是通过对墙体实测数据与建 模数据间的相关系数的计算,搜寻使相关系数达到最大时的模型,并将此模型 的参数作为墙体参数的估计。二是通过图像质量评价准则对成像结果进行评价, 通过不断改变假设的墙体参数,直到获得图像质量指标最高的成像结果,并将 此时假设的墙体参数作为真实参数的最优估计。上述两种方法都是通过搜索遍 历进行估计,计算速度慢。三是通过分析墙体参数误差与成像结果中目标位置 的关系,在一次观测中假设不同的墙体参数进行多次成像来获得目标位置的变 化轨迹,后通过不同天线轨迹的观测,将获得的多条目标轨迹的交叉点作为目 标的真实位置,并将此位置对应的墙体参数作为真实参数的估计值,但此方法需要多次成像,计算效率较低。四是通过改变天线之间的相对位置,测量每一 位置下墙体前后表面回波时延,并利用几何关系建立墙体参数的方程组,最后 通过最小二乘法求解墙体未知参数,此方法在墙后表面能量低的情况下失效。 五是利用支持向量机对未知墙体进行参数估计的方法,此方法受到目标位置、 大小和形状变化的影响较大并且无法预测多层墙体未知参数。
综上所述,目前没有一种快速并且精确的墙体参数估计方法。
发明内容
有鉴于此,本发明提供一种基于墙体前表面能量对消的墙体参数估计方法。 能够实现通过对回波信号的处理后,获得墙体参数。
一种基于墙体前表面能量对消的墙体参数估计方法,采用n对收发天线, 将一个接收天线接收到的回波信号进行脉冲压缩处理及去直达波处理,并以处 理后的信号回波能量最强的时刻作为中心点,将所述中心点左右对称的信号进 行对消处理,以消除回波信号中墙体前表面能量,凸显墙体后表面回波;对对 消处理后的信号做MUSIC处理,并取MUSIC处理后的信号回波能量最强的时 刻作为墙体后表面回波的时延τbn
获得两个不同接收天线的墙体后表面回波的时延,分别代入墙体参数估计 公式,利用牛顿迭代法求解即可得到墙体未知介电常数εr和墙体厚度d;
Figure BDA0002394431760000031
其中,xn为发射天线处发射信号到达墙体前表面的距离沿方位向的分量,yn为发射天线处发射信号从墙体前表面到达墙体后表面的距离沿方位向的分量, xn和yn为无需求解的未知量;r为天线与墙体前表面距离,Ln第n组收发天线间 距的一半,均由红外测距仪获得。
较佳地,适用于一发多收模式和多发多收模式进行数据采集。
较佳地,所述进行脉冲压缩处理及去直达波处理为:将接收天线接收到的 回波信号首先进行脉冲压缩处理,之后进行去直达波处理。
较佳地,所述进行脉冲压缩处理及去直达波处理为:首先利用空场景回波 信号与接收天线接收到的回波信号进行对消处理,之后对对消处理后的信号进 行脉冲压缩处理。
较佳地,通过设置距离门的方式,使墙体前表面的回波置零、墙后目标的 回波置零、中心点左右对称的回波数据互相抵消,实现信号的对消处理。
有益效果:
1、本发明基于接收到的回波信号的构成,对该信号进行分别处理,将该信 号的直达波、墙体前表面回波、墙后目标回波消除,获得墙体后表面回波信号 和噪声混杂的信号。之后对该信号进行分析,获得能量值最高的时刻,即为计 算墙体参数估计公式中,墙体后表面回波的时延精估计值。之后利用牛顿迭代 法求解墙体参数,这种方法可以快速并且精准的估计出墙体参数,无需遍历搜 索,且即使在墙后表面能量低的情况下,也不会失效。
附图说明
图1为墙体后表面回波示意图
图2为算法流程图
图3为实验场景图
图4(a)为一维距离向收发天线间距2m脉压结果图
图4(b)为一维距离向收发天线间距3m脉压结果图
图5(a)为空场景收发天线间距2m对消结果图
图5(b)为空场景收发天线间距2m对消结果图
图6(a)为墙体前表面能量收发天线间距2m对消结果图
图6(b)为墙体前表面能量收发天线间距3m对消结果图
图7(a)为MUSIC收发天线间距2m处理结果图
图7(b)为MUSIC收发天线间距3m处理结果图
具体实施方式
下面结合附图并举两个实施例,对本发明进行详细描述。
本发明提供一种基于墙体前表面能量对消的墙体参数估计方法,将原始回 波数据进行脉冲压缩处理,并利用空场景对消方法消除直达波,之后利用墙体 前表面能量对称性原则将信号进行抵消,凸显墙体后表面位置。然后利用MUSIC 方法对墙体后表面回波时延进行精确的估计,并利用该方法,获得两组墙体后 表面回波时延值,最后将两个数值代入墙体参数估计公式,利用牛顿迭代法求 解即可得到墙体未知介电常数和墙体厚度。墙体后表面回波信号示意图如图1 所示。
令第n个接收天线在第t时刻接收到的回波信号Sr(n,t)表示为:
Sr(n,t)=Srd(t-τdn)+Srf(t-τfn)+Srb(t-τbn)+Srt(t-τtn)+N(n,t) (1)
其中,Srd(t-τdn)为直达波,τdn为第n个接收天线处的直达波时延,Srf(t-τfn) 为墙体前表面回波,τfn为第n个接收天线处的墙体前表面回波时延,Srb(t-τbn)为 墙体后表面回波,τbn为第n个接收天线处的墙体后表面回波时延,Srt(t-τtn)为 墙后目标回波,τtn为第n个接收天线处的墙后目标回波时延,N(n,t)为第n个接 收天线处的噪声。由于墙体内多径和墙后目标多径能量较弱,此处不考虑其影 响。
具体实现方法为:
步骤一、为了实现回波信号的聚焦作用,对第n个接收天线在第t时刻接收 到的回波信号Sr(n,t)进行脉冲压缩处理:
Figure BDA0002394431760000051
其中,Spc(n,t)为脉压结果,u为卷积时间轴,Sref(n,t)为参考信号,一般取 发射信号的共轭。
步骤二、对步骤一中进行脉冲压缩处理后的信号进行去直达波处理,设空 场景回波信号可以表示为Srz(n,t),则步骤一中进行脉冲压缩处理后的信号与空 场景回波信号进行对消处理的方式如下:
Spcz(n,t)=Spc(n,t)-Srz(n,t) (3)
其中,Spcz(n,t)为第n个接收天线在第t时刻的空场景对消结果。
或者,本发明还可以对回波信号首先进行去直达波处理,之后对去直达波 处理后的信号进行脉冲压缩处理。
步骤三、对空场景对消后的信号进行墙体前表面能量对消
由于墙体前表面能量非常强,墙体后表面能量会被前表面能量的旁瓣淹没, 所以要进行前表面能量对消处理,从而凸显墙体后表面能量。为此,需要在步 骤二中获得的空场景对消结果中得到能量最强回波对应的时刻作为第n个接收 天线处的墙体前表面回波时延τfn,即:
Figure BDA0002394431760000061
通过以墙体前表面回波时延τfn为中心点,由于脉冲压缩后的墙体前表面能 量具有左右对称性,故采用左右对称使回波数据互相抵消,通过设置距离门使 墙体前表面的回波置零、使墙体后目标回波置零的方法,凸显出墙体后表面回 波,即:
Figure BDA0002394431760000062
其中,SFEC(n,t)为墙体前表面能量对消结果。针对不同场景可以选取合适的 tmax来抑制墙体后目标的影响,
Figure BDA0002394431760000063
d为墙体厚度,c为光速。
步骤四、将步骤三中获得的信号进行墙体后表面时延精估计,则需要对步 骤三中的结果做MUSIC处理
SMUSIC(n,t)=MUSIC[SFEC(n,t)] (6)
在上述结果中,取得峰值最大处对应的时刻即为墙体后表面回波时延精估 计值τbn,即:
Figure BDA0002394431760000064
步骤五、按照步骤一至步骤四的方法,获得两个不同接收天线的墙体后表 面回波的时延精估计值,之后代入墙体参数估计公式(8),利用牛顿迭代法求 解墙体参数。
由图1可以得到如下关系式
Figure BDA0002394431760000071
其中,εr为墙体相对介点常数,d为墙体厚度,xn为第n个发射天线处发射 信号到达墙体前表面的距离沿方位向的分量,yn为第n个发射天线处发射信号 从墙体前表面到达墙体后表面的距离沿方位向的分量。xn和yn为求解过程的中 间变量;r为天线与墙体前表面距离,Ln第n组收发天线间距的一半,由红外测 距仪即可获得。将步骤四得到的墙体后表面回波时延τbn代入上式,并利用牛顿 迭代法求解即可求得墙体厚度与相对介电常数。处理流程图如图2所示。
实施例1
利用所提技术对穿墙实验数据进行处理,以进一步证明所提技术的可行性 与有效性。穿墙实验参数如表1所示。
表1穿墙实验参数
参数/单位
信号波形 步进频连续波
发射功率/W 1
天线与墙体前表面距离/m 1.5
收发天线间距/m 2、3
墙体厚度/cm 29.5
墙体相对介电常数 5-7
实验场景图如图3所示,对穿墙实验数据进行脉冲压缩处理,得到的一维 距离向脉压结果如图4(a)和图4(b)所示。从图像可以看出,由于直达波的 影响,墙体前后表面回波均被淹没,对上述结果进行直达波去除处理得到的结 果如图5(a)和图5(b)所示。由图可知,墙体后表面回波被墙体前表面回波 和杂波淹没,对上述结果进行墙体前表面能量对消的结果如图6(a)和图6(b) 所示。由图可知,墙体后表面回波得以凸显,但位置估计存在误差。对上述结 果经过MUSIC处理得到的结果如图7(a)和图7(b)所示,由图可知墙体后 表面回波位置得以精确估计,将墙体后表面位置估计值代入式,即可得到估计 的墙体厚度及墙体相对介电常数如表2所示。
表2墙体参数估计算法估计结果
参数 墙体介电常数 墙体厚度(cm) 厚度估计误差(%)
参数值 6.2075 29.04 1.5
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保 护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等, 均应包含在本发明的保护范围之内。

Claims (5)

1.一种基于墙体前表面能量对消的墙体参数估计方法,采用n对收发天线,其特征在于,
将一个接收天线接收到的回波信号进行脉冲压缩处理及去直达波处理,并以处理后的信号回波能量最强的时刻作为中心点,将所述中心点左右对称的信号进行对消处理,以消除回波信号中墙体前表面能量,凸显墙体后表面回波;对对消处理后的信号做MUSIC处理,并取MUSIC处理后的信号回波能量最强的时刻作为墙体后表面回波的时延τbn
获得两个不同接收天线的墙体后表面回波的时延,分别代入墙体参数估计公式,利用牛顿迭代法求解即可得到墙体未知介电常数εr和墙体厚度d;
Figure FDA0002394431750000011
其中,xn为发射天线处发射信号到达墙体前表面的距离沿方位向的分量,yn为发射天线处发射信号从墙体前表面到达墙体后表面的距离沿方位向的分量,xn和yn为求解过程的中间变量;r为天线与墙体前表面距离,Ln第n组收发天线间距的一半,均由红外测距仪获得。
2.如权利要求1所述的墙体参数估计方法,其特征在于,适用于一发多收模式和多发多收模式进行数据采集。
3.如权利要求2所述的墙体参数估计方法,其特征在于,所述进行脉冲压缩处理及去直达波处理为:将接收天线接收到的回波信号首先进行脉冲压缩处理,之后进行去直达波处理。
4.如权利要求2所述的墙体参数估计方法,其特征在于,所述进行脉冲压缩处理及去直达波处理为:首先利用空场景回波信号与接收天线接收到的回波信号进行对消处理,之后对对消处理后的信号进行脉冲压缩处理。
5.如权利要求1或2或3或4所述的墙体参数估计方法,其特征在于,通过设置距离门的方式,使墙体前表面的回波置零、墙后目标的回波置零、中心点左右对称的回波数据互相抵消,实现信号的对消处理。
CN202010126142.8A 2020-02-28 2020-02-28 一种基于墙体前表面能量对消的墙体参数估计方法 Active CN111796249B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010126142.8A CN111796249B (zh) 2020-02-28 2020-02-28 一种基于墙体前表面能量对消的墙体参数估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010126142.8A CN111796249B (zh) 2020-02-28 2020-02-28 一种基于墙体前表面能量对消的墙体参数估计方法

Publications (2)

Publication Number Publication Date
CN111796249A true CN111796249A (zh) 2020-10-20
CN111796249B CN111796249B (zh) 2023-04-21

Family

ID=72806251

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010126142.8A Active CN111796249B (zh) 2020-02-28 2020-02-28 一种基于墙体前表面能量对消的墙体参数估计方法

Country Status (1)

Country Link
CN (1) CN111796249B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120235849A1 (en) * 2010-11-19 2012-09-20 Tatoian James Z Method and apparatus for through the wall radar imaging
US20150198713A1 (en) * 2014-01-13 2015-07-16 Mitsubishi Electric Research Laboratories, Inc. Method and System for Through-the-Wall Imaging using Compressive Sensing and MIMO Antenna Arrays
CN106772361A (zh) * 2016-11-30 2017-05-31 桂林电子科技大学 一种基于fpga的超宽带穿墙雷达成像算法的实现方法
CN107300694A (zh) * 2017-06-14 2017-10-27 南京信息工程大学 一种基于电磁波透射系数的未知墙体参数估计方法
CN109799499A (zh) * 2019-01-28 2019-05-24 沈阳航空航天大学 一种穿墙雷达墙体参数估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120235849A1 (en) * 2010-11-19 2012-09-20 Tatoian James Z Method and apparatus for through the wall radar imaging
US20150198713A1 (en) * 2014-01-13 2015-07-16 Mitsubishi Electric Research Laboratories, Inc. Method and System for Through-the-Wall Imaging using Compressive Sensing and MIMO Antenna Arrays
CN106772361A (zh) * 2016-11-30 2017-05-31 桂林电子科技大学 一种基于fpga的超宽带穿墙雷达成像算法的实现方法
CN107300694A (zh) * 2017-06-14 2017-10-27 南京信息工程大学 一种基于电磁波透射系数的未知墙体参数估计方法
CN109799499A (zh) * 2019-01-28 2019-05-24 沈阳航空航天大学 一种穿墙雷达墙体参数估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CHRISTOPHER THAJUDEEN 等: "Wall Parameters Estimation using a Hybrid Time-Delay-Only and Reflected Wave Ratio Technique", 《AP-S 2013》 *
FOAD FEREIDOONY 等: "Efficient Method for Estimation of the thicknesses and complex dielectric constant of wall", 《THE 8TH EUROPEAN CONFERENCE ON ANTENNAS AND PROPAGATION (EUCAP 2014)》 *
PAVEL PROTIVA 等: "An Improved Technique for Estimating Wall Parameters from Time-Delay Through-Wall Radar Measurements", 《22ND INTERNATIONAL CONFERENCE RADIOELEKTRONIKA 2012》 *
PAVEL PROTIVA 等: "Estimation of Wall Parameters From Time-Delay-Only Through-Wall Radar Measurements", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 *

Also Published As

Publication number Publication date
CN111796249B (zh) 2023-04-21

Similar Documents

Publication Publication Date Title
Jin et al. Image-domain estimation of wall parameters for autofocusing of through-the-wall SAR imagery
Ahmad et al. Synthetic aperture beamformer for imaging through a dielectric wall
CN107861123B (zh) 一种穿墙雷达在复杂环境下对多运动目标实时跟踪的方法
CN109799499B (zh) 一种穿墙雷达墙体参数估计方法
CN109444896B (zh) 水声阵列定位系统及其定位方法
JP2011505173A (ja) サーチボリュームの内容物を測定する方法および装置
Sakamoto et al. Fast imaging method for security systems using ultrawideband radar
KR101897763B1 (ko) 레이더를 이용한 거리 측정 방법 및 장치
CN111045006B (zh) 一种基于多成像字典融合的拐角隐蔽目标成像方法
CN109471097B (zh) 一种穿墙雷达信号优化处理方法及装置
CN109696672B (zh) 一种基于空间结构关联性的高分辨率穿墙雷达成像方法
Aftanas et al. Efficient method of TOA estimation for through wall imaging by UWB radar
CN113447915B (zh) 一种适用于复杂多径环境下的超宽带层析成像方法
CN113064160B (zh) 一种u型建筑布局下目标定位方法
CN107219522B (zh) 一种椭圆-双曲线联合的穿墙雷达目标定位方法
CN111796249B (zh) 一种基于墙体前表面能量对消的墙体参数估计方法
Lei A time-domain beamformer for UWB through-wall imaging
Singh et al. A novel technique for contrast target detection in through-the-wall radar images
CN113960558A (zh) 基于多输入多输出雷达的非视距目标定位方法及系统
CN114152943B (zh) 一种基于超宽带穿墙雷达的两阶段墙体参数估计方法
CN111443333A (zh) 一种基于频谱综合的多方位匹配信号生成方法
CN115267721B (zh) 一种基于双频sar的地面动目标径向速度估计方法
CN113341408B (zh) 一种基于穿墙雷达杂波抑制的成像方法及系统
Muqaibel et al. Dynamic joint reconstruction of walls and targets in through-the-wall radar imaging
Ahmad Multi-location wideband through-the-wall beamforming

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