CN112818514B - 无线层析成像方法及结合水平距离衰减的椭圆权重模型 - Google Patents
无线层析成像方法及结合水平距离衰减的椭圆权重模型 Download PDFInfo
- Publication number
- CN112818514B CN112818514B CN202110007809.7A CN202110007809A CN112818514B CN 112818514 B CN112818514 B CN 112818514B CN 202110007809 A CN202110007809 A CN 202110007809A CN 112818514 B CN112818514 B CN 112818514B
- Authority
- CN
- China
- Prior art keywords
- wireless tomography
- wireless
- linear equation
- ith
- link
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/04—Architectural design, interior design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种无线层析成像方法及结合水平距离衰减的椭圆权重模型,无线层析成像方法包括如下步骤:根据检测区域内传感器节点的数量和位置进行建模,得到无线层析成像的线性方程;将正则约束项加入到无线层析成像的线性方程中,得到目标函数;对目标函数求导并将其导数设为零,得到无线层析成像线性方程的Tikhonov正则化解;在无线层析成像方法的自适应椭圆权重模型的基础上引入自适应水平距离衰减因子,通过对像素在链路上的投影到其起始节点和终节点的距离做归一化处理,使得像素权重成为归一化后距离的函数,保证目标位置信息与测量信息的一致性。
Description
技术领域
本发明涉及无线层析成像技术领域,具体涉及一种无线层析成像方法及结合水平距离衰减的椭圆权重模型。
背景技术
无线层析成像技术能够基于静态无线传感器之间的接收信号强度(RSS,ReceivedSignal Strength)测量值对传感器探测区域中的无线传播环境变化进行成像。因此,可以检测室内和室外环境中人体活动引起的RSS变化实现目标定位和跟踪,且具有低功率和低成本的特点。
在RTI系统中,椭圆权重模型建模的准确性是影响系统定位性能的关键因素之一。经典的椭圆权重模型由Wilson和Patwari首次提出,但是该模型存在一定的局限性,如椭圆内的像素权重相同,不能分辨椭圆覆盖区域内实际测量环境中每个像素反映的信号传输衰减的变化。
2014年,Hamilton提出了逆面积椭圆模型,在该模型中,信号传播路径越短,视距(LOS,Line of Sight)路径附近的像素权重越大,反之亦然。为了体现模型特性,每个像素权重建模为包含收发节点和该像素的最小椭圆面积的倒数。由于该模型可以区分每个像素在椭圆内的贡献,因此能够获得较其它模型更高的定位精度。然而,模型的计算复杂度太高导致无法实时定位。
在适用于短距离小覆盖的恒定离心率椭圆权重模型中,不同长度链路对应不同的椭圆短轴参数,可以用较少的像素重建图像,较大程度上降低了噪声干扰。然而,如何确定椭圆内的权重却没有被考虑。
基于几何的椭圆模型采用视距和非视距两种加权方式,将一个椭圆分割成多个不同的区域,并结合正交匹配追踪(OMP,Orthogonal Matching Pursuit)算法消除图像重建中的额外亮点。与经典权重模型相比,单目标的定位精度提高23.8%。但是,在视距路径和非视距路径通信中,代表障碍的模型参数需要通过实验来确定。当监测环境变化后,模型参数需要重新校正。
自适应椭圆权重模型同时考虑椭圆覆盖范围的调整与像素权重的选择,在垂直于两节点之间视距路径方向上,距离视距路径越远,权重越小;而且在链路较短时增大椭圆短轴参数得到更丰富的像素点衰减信息。室内仿真表明,自适应椭圆模型可以减弱噪声伪迹,提高定位性能。除此以外,还有相关以菲涅耳衍射理论为建模基础提出的椭圆模型,如马鞍模型和衍射模型。但这些模型给出的不再是单纯的权值,而是具体的RSS变化量,已无法直接用到现有无线层析成像方法中。
发明内容
本发明的目的是提供一种无线层析成像方法及结合水平距离衰减的椭圆权重模型,以解决现有技术中无线层析成像方法的目标位置信息与测量信息的一致性差的问题。
为实现上述目的,本发明的采用如下技术方案:
一种无线层析成像方法,包括如下步骤:
步骤一:根据检测区域内传感器节点的数量和位置进行建模,得到无线层析成像的线性方程;
步骤二:将正则约束项加入到无线层析成像的线性方程中,得到目标函数;
步骤三:对目标函数求导并将其导数设为零,得到无线层析成像线性方程的Tikhonov正则化解;
在无线层析成像的线性方程中,设其阴影权重矩阵中第i行第j列的元素为Wij,该元素为第i条链路中第j个像素点的权重值,则
Wij=e-h×e-s
否则
Wij=0
l1和l2分别是第j个像素点在第i条链路上的投影到其起始节点和终节点的距离,di是第i条链路的长度,dmax为所有链路长度中的最大值。
进一步的,所述无线层析成像的线性方程为
y=Wx+n
其中y表示监测网络内所有链路上的变化量,W表示阴影权重矩阵,x表示信号衰落值,n表示噪声。
进一步的,所述目标函数为
其中Q是Tikhonov正则化矩阵,α为正则化参数,y表示监测网络内所有链路上的变化量,W表示阴影权重矩阵,x表示信号衰落值。
进一步的,无线层析成像线性方程的Tikhonov正则化解为
x=(WTW+αQTQ)-1WTy
其中Q是Tikhonov正则化矩阵,α为正则化参数,y表示监测网络内所有链路上的变化量,W表示阴影权重矩阵。
结合水平距离衰减的椭圆权重模型,设其阴影权重矩阵中第i行第j列的元素为Wij,该元素为第i条链路中第j个像素点的权重值,则
Wij=e-h×e-s
否则
Wij=0
l1和l2分别是第j个像素点在第i条链路上的投影到其起始节点和终节点的距离,di是第i条链路的长度,dmax为所有链路长度中的最大值。本发明的有益效果:
本发明所提供的技术方案,在无线层析成像方法的自适应椭圆权重模型的基础上引入自适应水平距离衰减因子,通过对像素在链路上的投影到其起始节点和终节点的距离做归一化处理,使得像素权重成为归一化后距离的函数,保证目标位置信息与测量信息的一致性。
附图说明
图1是本发明实施例1中无线层析成像方法的流程图;
图2是本发明实施例1中检测区域设置传感器节点的示意图;
图3是本发明实施例1中椭圆模型的示意图;
图4是本发明实施例1中在室外环境下采用四种无线层析成像方法得到的位置1的重建图像;
图5是本发明实施例1中在室外环境下采用四种无线层析成像方法得到的位置2的重建图像;
图6是本发明实施例1中在室外环境下采用四种无线层析成像方法得到的位置3的重建图像;
图7为本发明实施例1中室外三个位置下的目标定位MSE曲线图;
图8是本发明实施例1中在室内环境下采用四种无线层析成像方法得到的位置1的重建图像;
图9是本发明实施例1中在室内环境下采用四种无线层析成像方法得到的位置2的重建图像;
图10是本发明实施例1中在室内环境下采用四种无线层析成像方法得到的位置3的重建图像;
图11为本发明实施例1中室内三个位置下的目标定位MSE曲线图。
具体实施方式
实施例1:
本实施例提供一种无线层析成像方法,在其自适应椭圆权重模型的基础上引入自适应水平距离衰减因子,通过对像素在链路上的投影到其起始节点和终节点的距离做归一化处理,使得像素权重成为归一化后距离的函数,保证目标位置信息与测量信息的一致性。
本实施例所提供的无线层析成像方法,其流程如图1所示,包括如下步骤
步骤一:根据检测区域内传感器节点的数量和位置进行建模,得到无线层析成像的线性方程。
设检测区域共有L个传感器节点,如图2所示,在任意两个传感器节点之间建立LOS路径,每条LOS路径为一条链路,得到M条链路,则
M=L(L-1)
得到的无线层析成像的线性方程为
y=Wx+n
其中y表示监测网络内所有链路上的变化量,W表示阴影权重矩阵,x表示信号衰落值,n表示噪声。
步骤二:将正则约束项加入到无线层析成像的线性方程中,得到目标函数。
设目标函数为f(x),则
其中Q是Tikhonov正则化矩阵,α为正则化参数。
步骤三:对目标函数求导并将其导数设为零,得到无线层析成像线性方程的Tikhonov正则化解。
x=(WTW+αQTQ)-1WTy
阴影权重矩阵W的获取方法为:
在无线层析成像的线性方程中,设其阴影权重矩阵中第i行第j列的元素为Wij,该元素为第i条链路中第j个像素点的权重值,则
Wij=e-h×e-s
否则
Wij=0
l1和l2分别是第j个像素点在第i条链路上的投影到其起始节点和终节点的距离,di是第i条链路的长度,dmax为所有链路长度中的最大值。
仿真实验:
设第一种无线层析成像方法所采用的椭圆模型为:
当dij(1)+dij(2)<di+λ时
否则
Wij=0
设第二种无线层析成像方法所采用的椭圆模型为:
当dij(1)+dij(2)<di/ε时
否则
Wij=0
其中ε是椭圆的离心率。
设第三种无线层析成像方法所采用的椭圆模型为
Wij=e-h
否则
Wij=0
设第四种无线层析成像方法为本实施例所提供的无线层析成像方法。
为了评估各种无线层析成像方法的定位性能,分别在室内环境和室外环境下选取多个测试位置进行定位性能测试,表1给出了椭圆参数λ、像素数、节点数和网格划分像素尺寸等实验参数设置,节点数即为监测区域中传感器的数量。为了减小地面反射的影响,测试时将传感器节点置于距地面高度为1m的支架上;为了抑制来自墙体,天花板等外部监测区域的反射信号,采用了水平波束宽度为110°,垂直波束宽度为30°的平面定向天线。节点完全支持IEEE802.15.4协议,最大传输功率为4.5dBm,足以覆盖整个监控区域。为了快速测量所有链路的RSS,本实施例采用了Wilson提出的令牌环式通信协议,在该协议中为了避免冲突,为每个节点分配唯一的ID号作为标识,发送顺序取决于节点的编号。当所有节点完成信号传输后,得到所有链路的RSS测量信息。
表1
室外环境下定位性能测试:
室外环境下监测区域为4.5m×4.5m的正方形,相邻传感器节点间隔1.5m。在实验中选取了(3.0,3.0),(2.3,2.3)和(1.5,0.9)三个不同位置分析定位精度,分别定义为位置1、位置2和位置3。
则三个位置下的均方根误差平均值为
对于给定的三个目标位置,分别使用各种无线层析成像方法进行图像重建,得到的结果分别如图4、图5和图6所示。从图4、图5和图6可以明显看出,对于所测试的三个位置,可以根据最亮区域能量的聚集度来评价定位性能,从这一点看,第四种无线层析成像方法对应的重建图像最亮区域能量聚集度最高,且由噪声引起的伪影区域小。
为了更好地体现所提模型的定位性能,图7描述了室外三个位置下的目标定位均方误差曲线。表2给出了相应的均方根误差。结合均方误差和均方根误差可知,第四种无线层析成像方法相比第三种无线层析成像方法,目标定位精度提升了约6%,提供了更高的目标定位精度。
表2
室内环境下定位性能测试:
室内环境下监测区域为6m×6m的正方形,相邻传感器节点间隔1.5m。在实验中选取了(3.0,4.8),(1.0,1.2)和(4.8,4.8)三个不同目标位置分析定位精度,分别定义为位置1、位置2和位置3。
对于给定的三个目标位置,分别使用各种无线层析成像方法进行图像重建,得到的结果分别如图8、图9和图10所示。从图8可以明显看出,与图4所示的室外环境成像结果相比,因室内环境测试时的多径效应影响,室内不同位置重建图像中的伪影增强,少数图像伪影区域的平均灰度值较高,这在实际应用时可能会被误判为有多目标存在或增加单目标定位的偏差。此外第四种无线层析成像方法对应的目标重建图像中的伪影区域的高灰度值像素较少,无线层析成像方法性能相对稳定。图11描述了不同椭圆权重模型在室内三个目标位置下的MSE曲线,表3给出了定位均方误差和均方根误差。可见,与室外相比,室内环境下的定位性能有所下降,这是由于室内多径效应的影响造成的。与其他无线层析成像方法相比,第四种无线层析成像方法的定位性能分别改善了15.4%,10.2%和2.34%。
表3
实验结论:
本实施例所提供的无线层析成像方法,其采用了结合水平距离衰减的椭圆权重模型,该模型同时考虑了成像椭圆内像素点距离链路起始节点和终节点在垂直方向上和水平方向上的距离变化,并建立了像素权值与距离之间的函数关系,从而能够更细粒度的区分椭圆内每个像素对目标遮挡造成的路径损耗的贡献度,从而实现测量信息与目标位置之间的对应关系更加匹配实际的信号传输衰减特性。室内环境和室外环境的实测表明了本实施例的椭圆权重模型能够有效减少重建图像伪影,提高定位精度,而且在实际应用中不仅椭圆权重能够自适应选择,椭圆的大小也能够随着节点距离的变化而调整,椭圆范围可以控制进入椭圆内的噪声干扰的数量,并决定用于目标图像重建的椭圆内像素的数量。上述椭圆权重模型的综合设计将直接影响无线层析成像方法目标图像重建的质量和定位性能。
实施例2:
本实施例提供一种结合水平距离衰减的椭圆权重模型,与上述实施例1中无线层析成像方法中的椭圆权重模型相同,该模型已在上述实施例1中做了详细介绍,这里不多做说明。
以上公开的本发明的实施例只是用于帮助阐明本发明的技术方案,并没有详尽叙述所有的细节,也不限制该发明仅为所述的具体实施方式。显然,根据本说明书的内容,可作很多的修改和变化。本说明书选取并具体描述这些实施例,是为了更好地解释本发明的原理和实际应用,从而使所属技术领域人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。
本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不会使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (2)
1.一种无线层析成像方法,包括如下步骤:
步骤一:根据检测区域内传感器节点的数量和位置进行建模,得到无线层析成像的线性方程;
步骤二:将正则约束项加入到无线层析成像的线性方程中,得到目标函数;
所述无线层析成像的线性方程为
y=Wx+n
其中y表示监测网络内所有链路上的变化量,W表示阴影权重矩阵,x表示信号衰落值,n表示噪声;
所述目标函数为
其中Q是Tikhonov正则化矩阵,α为正则化参数,y表示监测网络内所有链路上的变化量,W表示阴影权重矩阵,x表示信号衰落值;
步骤三:对目标函数求导并将其导数设为零,得到无线层析成像线性方程的Tikhonov正则化解;
其特征在于,在无线层析成像的线性方程中,设其阴影权重矩阵中第i行第j列的元素为
Wij,该元素为第i条链路中第j个像素点的权重值,则
Wij=e-h×e-s
否则
Wij=0
l1和l2分别是第j个像素点在第i条链路上的投影到其起始节点和终节点的距离,di是第i条链路的长度,dmax为所有链路长度中的最大值。
2.根据权利要求1所述的无线层析成像方法,其特征在于,无线层析成像线性方程的Tikhonov正则化解为
x=(WTW+αQTQ)-1WTy
其中Q是Tikhonov正则化矩阵,α为正则化参数,y表示监测网络内所有链路上的变化量,W表示阴影权重矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110007809.7A CN112818514B (zh) | 2021-01-05 | 2021-01-05 | 无线层析成像方法及结合水平距离衰减的椭圆权重模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110007809.7A CN112818514B (zh) | 2021-01-05 | 2021-01-05 | 无线层析成像方法及结合水平距离衰减的椭圆权重模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112818514A CN112818514A (zh) | 2021-05-18 |
CN112818514B true CN112818514B (zh) | 2023-02-07 |
Family
ID=75857389
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110007809.7A Active CN112818514B (zh) | 2021-01-05 | 2021-01-05 | 无线层析成像方法及结合水平距离衰减的椭圆权重模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112818514B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104486833A (zh) * | 2014-11-24 | 2015-04-01 | 北京理工大学 | 基于运动预测删除干扰链路的室内无线层析成像增强方法 |
WO2017130018A1 (en) * | 2016-01-25 | 2017-08-03 | Adaptix Ltd | Medical imaging system with a fixed array of x-ray detectors and a fixed array of x-ray emitters for producing a digital 3-dimensional image |
CN108828510A (zh) * | 2018-05-29 | 2018-11-16 | 南京信息工程大学 | 一种基于渐变阴影权重模型的射频层析成像定位方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103024900B (zh) * | 2013-01-07 | 2015-03-11 | 北京理工大学 | 无线层析成像网络的分数阶傅立叶多径伪像修正定位方法 |
US10386449B2 (en) * | 2016-02-10 | 2019-08-20 | United States Of America As Represented By The Secretary Of The Air Force | System and method for radio tomographic image formation |
US9911208B2 (en) * | 2016-04-11 | 2018-03-06 | Toshiba Medical Systems Corporation | Apparatus and method of iterative image reconstruction using regularization-parameter control |
CN106231552B (zh) * | 2016-08-02 | 2019-09-27 | 南京信息工程大学 | 一种基于混合方式的非合作目标定位方法 |
CN106940895B (zh) * | 2017-01-12 | 2021-06-29 | 北京理工大学 | 一种应用于无线层析成像系统的降质函数的估计方法 |
CN108490387B (zh) * | 2018-03-08 | 2021-11-30 | 河南工业大学 | 一种无线层析成像用网格像素衰减值自适应确定方法 |
CN110766561A (zh) * | 2019-09-11 | 2020-02-07 | 河南工业大学 | 基于无线层析成像技术的异常粮情检测方法 |
-
2021
- 2021-01-05 CN CN202110007809.7A patent/CN112818514B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104486833A (zh) * | 2014-11-24 | 2015-04-01 | 北京理工大学 | 基于运动预测删除干扰链路的室内无线层析成像增强方法 |
WO2017130018A1 (en) * | 2016-01-25 | 2017-08-03 | Adaptix Ltd | Medical imaging system with a fixed array of x-ray detectors and a fixed array of x-ray emitters for producing a digital 3-dimensional image |
CN108828510A (zh) * | 2018-05-29 | 2018-11-16 | 南京信息工程大学 | 一种基于渐变阴影权重模型的射频层析成像定位方法 |
Non-Patent Citations (1)
Title |
---|
基于传感器网络的无线层析成像方法;刘珩等;《北京理工大学学报》;20131115(第11期);参见第1155-1160页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112818514A (zh) | 2021-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111954228B (zh) | 波束管理方法、装置、电子设备及计算机可读存储介质 | |
JP4065437B2 (ja) | 無線信号強度を使用する装置位置測定システムのキャリブレーション | |
US20150230100A1 (en) | System and method for wireless positioning in wireless network-enabled environments | |
WO2019007388A1 (en) | SYSTEM AND METHOD FOR DEEP LEARNING AND WIRELESS NETWORK OPTIMIZATION USING DEEP LEARNING | |
US20060055592A1 (en) | Antenna measurement systems | |
Will et al. | The membership degree min-max localization algorithm | |
CN111007546A (zh) | 基于北斗伪距差分与无线的室内外融合定位技术 | |
CN110361705A (zh) | 一种相控阵天线近场迭代校准方法 | |
CN112907690A (zh) | 无线层析成像系统及基于Tikhonov-lp范数正则化的图像重建方法 | |
CN112818514B (zh) | 无线层析成像方法及结合水平距离衰减的椭圆权重模型 | |
KR102211847B1 (ko) | 경로 손실 지수 예측 시스템 및 방법 | |
CN108169709B (zh) | 一种基于分布式天线的室内三维定位方法 | |
CN105072582B (zh) | 基于rss分布的距离自适应的无线传感网络被动式定位方法 | |
Sala et al. | Using neural networks and Active RFID for indoor location services | |
CN108490387B (zh) | 一种无线层析成像用网格像素衰减值自适应确定方法 | |
Pascacio et al. | A collaborative approach using neural networks for BLE-RSS lateration-based indoor positioning | |
Hou et al. | FCLoc: A novel indoor Wi-Fi fingerprints localization approach to enhance robustness and positioning accuracy | |
CN107422296B (zh) | 无线层析成像用网格像素衰减值确定方法 | |
CN108107421A (zh) | 一种室内测距方法和装置 | |
CN111313999A (zh) | 一种基于零稀疏数据驱动权重模型的射频层析成像方法 | |
Zhou et al. | An enhanced RSSI-based framework for localization of Bluetooth devices | |
CN110996248A (zh) | 一种基于rss的凸优化无线传感器网络定位方法 | |
CN114252871B (zh) | 一种基于机器学习的雷达测量精度补偿方法 | |
JP7468107B2 (ja) | 電波強度推定装置、位置推定システム、位置推定方法、電波強度推定方法およびプログラム | |
CN114900234A (zh) | 一种太赫兹频谱环境地图构建方法及设备 |
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 |