CN116047532A - 一种测距方法和测距系统 - Google Patents

一种测距方法和测距系统 Download PDF

Info

Publication number
CN116047532A
CN116047532A CN202111260507.7A CN202111260507A CN116047532A CN 116047532 A CN116047532 A CN 116047532A CN 202111260507 A CN202111260507 A CN 202111260507A CN 116047532 A CN116047532 A CN 116047532A
Authority
CN
China
Prior art keywords
ranging
spad
sequence
time
pulse signals
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
Application number
CN202111260507.7A
Other languages
English (en)
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.)
Ningbo Abax Sensing Electronic Technology Co Ltd
Original Assignee
Ningbo Abax Sensing Electronic Technology Co Ltd
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 Ningbo Abax Sensing Electronic Technology Co Ltd filed Critical Ningbo Abax Sensing Electronic Technology Co Ltd
Priority to CN202111260507.7A priority Critical patent/CN116047532A/zh
Priority to PCT/CN2022/126423 priority patent/WO2023071908A1/zh
Publication of CN116047532A publication Critical patent/CN116047532A/zh
Pending legal-status Critical Current

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
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06Systems determining position data of a target
    • G01S17/08Systems determining position data of a target for measuring distance only
    • G01S17/10Systems determining position data of a target for measuring distance only using transmission of interrupted, pulse-modulated waves
    • 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/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483Details of pulse systems
    • G01S7/486Receivers
    • G01S7/4865Time delay measurement, e.g. time-of-flight measurement, time of arrival measurement or determining the exact position of a peak

Landscapes

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

Abstract

申请提供一种测距方法和测距系统,包括:发射N个脉冲信号;存储所述N个脉冲信号及其发射时间,形成序列A;接收回波信号,所述回波信号触发SPAD雪崩;存储所述SPAD被触发的时刻,形成序列B;对所述序列A与所述序列B进行移动时间窗口卷积(Time correlated moving convolution,TCMC)计算,并根据所述计算结果得到测距结果。本申请提出的测距方法和测距系统,提高了测距系统的抗干扰性能,并且因为SPAD相比于传统的APD器件,可以连接更高精度的TDC,因此,本申请提出的测距系统和测距方法相比于传统的测距系统和测距方法具有更高的精度。

Description

一种测距方法和测距系统
技术领域
本申请涉及测距领域,更具体的,涉及一种测距方法和测距系统。
背景技术
随着三维成像信息技术的迅速发展,特别是建筑测量、室内定位和导航、立体影像和辅助生活环境应用对飞行时间(Time-of-Flight,TOF)成像提出了迫切需求。飞行时间测距法(Time of flight,TOF),其原理是通过给目标物连续发送光脉冲,然后用传感器接收从物体返回的光,通过探测光脉冲的飞行(往返)时间来得到目标物距离。
直接飞行时间探测(Direct Time of flight,DTOF)作为TOF的一种,DTOF技术通过计算光脉冲的发射和接收时间,直接获得目标距离,具有原理简单,信噪比好、灵敏度高、精确度高等优点,受到了越来越广泛的关注。
然而,在同一场景下,多台TOF设备同时工作时,会产生干扰现象,导致测距精度不准确,现有技术中主要的解决方法有:频率调制技术,时钟同步技术等。相应的,在TOF设备的接收端,需要根据发射端的频率或者同步时钟进行相应的调制,以达到测距的目的。此外,现有技术中的使用单光子雪崩二极管(Single Photon Avalanche Diode,SPAD)的TOF设备,也会使用脉冲激光编码方式进行测距过程中的抗干扰。
然而,传统的使用脉冲激光编码方式的虽然具有抗干扰性能,但是接收后直接进行统计,这样给后续处理电路带来极大的复杂度,同时,还需要更大量的存储开销,因此,需要提出一种可以节省存储空间的具有抗干扰性能的测距方法和测距系统。
发明内容
本申请的目的在于,针对上述现有技术中的不足,提供一种测距方法和测距系统。
为实现上述目的,本申请实施例采用的技术方案如下:
第一方面,本申请提供一种测距方法,其特征在于,包括:
发射N个脉冲信号;
存储所述N个脉冲信号及其发射时间,形成序列A;
接收回波信号,所述回波信号触发SPAD雪崩;
存储所述SPAD被触发的时刻,形成序列B;
对所述序列A与所述序列B进行移动时间窗口卷积计算,并根据所述计算结果得到测距结果。
可选地,所述N个脉冲信号由随机数生成。
可选地,所述N个脉冲信号的相邻脉冲的时间间隔随机。
可选地,所述N个脉冲的相邻脉冲的时间间隔与所述SPAD的死区时间成正比。
可选地,所述N个脉冲的个数与测距范围成正比。
第二方面,本申请提供一种测距系统,其特征在于,包括:
发射模块,用于发射N个脉冲信号;
SPAD或SPAD阵列,用于探测光信号;
接收模块,用于接收回波信号,所述回波信号触发所述SPAD雪崩;
存储模块,分别存储所述N个脉冲信号及其发射时间,形成序列A,所述SPAD被触发的时刻,形成序列B;
计算模块,根据所述序列A与所述序列B进行移动时间窗口卷积计算,并根据所述计算结果得到测距结果。
可选地,所述N个脉冲信号由随机数生成。
可选地,所述N个脉冲信号的相邻脉冲的时间间隔随机。
可选地,所述N个脉冲的相邻脉冲的时间间隔与所述SPAD的死区时间成正比。
可选地,所述N个脉冲的个数与测距范围成正比。
本申请提供一种测距方法和测距系统,提高了测距方法和测距系统的抗干扰性能,并且因为SPAD相比于传统的APD器件,可以连接更高精度的TDC,因此,本申请提出的测距系统和测距方法相比于传统的测距方法和测距系统具有更高的精度,此外,SPAD的状态可以直接用二进制数字0或1直接表示,相比于APD更加节省存储空间,也就是说,使用SPAD器件数字化之后,该测距系统可以减小数据量,支持更多并行处理的通道数,支持更大规模的阵列。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本申请的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为一种传统的探测系统示意图;
图2为一种传统的探测系统的直方图;
图3为本申请实施例提供的一种测距方法的流程图;
图4为本申请实施例提供的一种测距系统的模块框图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
飞行时间测距法(Time of flight,TOF),其原理是通过给目标物连续发送光脉冲,然后用传感器接收从物体返回的光,通过探测光脉冲的飞行(往返)时间来得到目标物距离。
而直接飞行时间探测(Direct Time of flight,DTOF)作为TOF的一种,DTOF技术通过计算光脉冲的发射和接收时间,直接获得目标距离,具有原理简单,信噪比好、灵敏度高、精确度高等优点,受到了越来越广泛的关注。
一般地,在一些DTOF测距应用中,可以使用包括单光子检测器(例如单光子雪崩二极管,Single Photon Avalanche Diode,SPAD,或单光子雪崩二极管阵列)在内的光电探测器阵列来执行光电探测的任务。一个或多个光电探测器可以限定阵列的探测器像素。SPAD阵列可以在可能需要高灵敏度和定时分辨率的成像应用中用作固态光电探测器。SPAD基于半导体结(例如,p-n结),例如,当通过或响应于具有期望脉冲宽度的选通信号而被偏置到其击穿区域之外时,该半导体结可以检测入射光子。高的反向偏置电压会产生足够大小的电场,从而使引入器件耗尽层的单个电荷载流子可以通过碰撞电离引起自持雪崩。
一般地,光在极其微弱时会离散成一个个的光子,称为单光子。单光子信号由于强度微弱且粒子性显著,常规技术难以对其检测,被认为是光电探测技术的极限,DTOF技术克服了光电探测技术的难点,实现单光子检测。
DTOF,是直接根据脉冲发射和接收的时间差来测算距离。激光发射的瞬间,电子时钟被激活。光束操纵单元将脉冲引导到所需方向。脉冲从探测目标反射回来,一部分被光电探测器所接收。在响应中,与前端电子器件连接的光电探测器产生电信号,从而使时钟生效。通过测量飞行时间Δt,计算出与反射物体之间的距离d,计算式为d=cΔt/2,其中c指光在介质中的速度。
图1为传统的一种探测系统示意图,图2为传统的测距结果的直方统计图,即直方图。如图1,示意了探测系统获取目标的基本原理,处理单元120控制光源110发出发射光,光源可以为LED或者激光源,此处为了考虑人眼安全等一般选择光源为具有近红外波长的激光源,激光源可以选择VSCEL阵列型激光源,此处并不限定,光源110至少部分单元发射探测光,当然也包含了全部发出发射光和部分发出发射光的场景此处不做限定,为了实现能量集中和系统工作的高效性,通常使用部分输出单元输出发射光,接收模块130包含SPAD单光子雪崩二极管探测单元,这样可以实现对于微弱能量下的准确探测。
在实际的探测过程中,光源发射具有一定脉宽的脉冲激光例如几纳秒级别,通过数以万次的发射由统计结果获得最终的目标距离信息,脉冲激光经过探测目标140反射返回处于包含雪崩状态SPAD的阵列型接收模块,其中处于雪崩状态的探测单元可以接收返回的信号,返回的信号可以是发射脉冲光束的回波并形成光子信号,当接收模块的SPAD被施加高的偏置电压时,接收模块的SPAD单元处于雪崩状态,因此可以感测背景光或者返回信号光的光子,再利用处理模块的统计功能,对于触发事件进行大量的统计,通过统计结果构造出直方图,如图2所示,再对于直方图的触发概率高的时间段信息输出即可实现对于探测目标140的距离信息获取。
然而,对于直方图统计的测距方法,需要消耗大量的内存,因此本申请提出了一种测距方法和测距系统。
如图3所示,是本申请提供的一种测距方法和测距系统的流程图。下面详细说明本申请提供的测距方法的具体步骤。
一般地,测距系统中会包括发射模块和接收模块,发射模块用于发射探测信号,接收模块用于接收回波信号,并根据该信号进行处理得到距离信息。
S1,发射N个脉冲信号。
首先,发射模块发射N个脉冲信号,该脉冲信号的脉冲间隔不是固定值,可以是随机值,例如,N可以是1024,即一次探测发射1024个脉冲,脉冲平均间隔为10ns左右,间隔分为8ns和12ns。
S2,存储N个脉冲信号及其发射时间,形成序列A。
将上一步骤中发射出的N个脉冲信号存储在电路中的存储模块,形成序列A。
例如,将上一步骤中1024个脉冲存储,脉冲间隔为8ns的代表0,脉冲间隔是12ns的代表1,这样会产生一列长度为1024个数字的伪随机序列,并将其存储至存储模块中。
S3,接收回波信号,所述回波信号触发SPAD雪崩。
接收模块接受回波信号,回波信号可以引起接收模块中的SPAD器件雪崩,产生雪崩信号,该信号用于记录SPAD是否被触发。
S4,存储SPAD被触发的时刻形成序列B。
当接收模块接受到回波信号以后,SPAD被触发,存储模块则记录SPAD被触发的时刻,并记录和存储该时刻。
例如,将一个探测周期的时间设置为11600ns,覆盖了1333.33ns的飞行时间,对应200m范围的测距,以及1024个光脉冲所占的10240ns时间。
时间从0时刻开始时,时间精度设置为1ns,记录SPAD是否被触发,0表示SPAD被触发,1表示SPAD没有被触发,记录到一个探测周期11600ns时间结束。
S5,对序列A与序列B进行移动时间窗口卷积(Time correlated movingconvolution,TCMC),得到测距结果。
对于序列A和B进行移动窗口卷积,移动的窗口时间为0~1333ns,即代表0~200m量程。每次移动一个最小时间(Δt)得到一个卷积结果,即
f(Δt)=A(t+Δt)*B(t)。
示例性地,卷积的操作如下流程:
第一步,设置Δt=0ns时,序列A的全部(10240bit)和序列B的第一个至10240个bit一一相乘,再相加,得到C1即:C_1=sum[A(1:10240).*B(1:10240)];
第二步,设置Δt=1ns时,列A的全部(10240bit)和序列B的第二个至10240个bit一一相乘,再相加,得到C2即:C_2=sum[A(1:10240).*B(2:10240)];
循环往复,可以得到整个代表卷积结果的C序列,即:C=[C_1,C_2,C_3,……C_n];C_n=sum[A(1:10240).*B(n:10240+n)]
最后,在C序列中找最大值Cm,则回波时间即为m ns。最终将时间转换为距离,就可以直接得到探测距离。
因此,本申请提供的测距方法可以应用在多台测量设备同时工作的场景下,提高了设备的抗干扰性能,又因为SPAD相比于传统的APD器件,可以连接更高精度的TDC,因此,本申请提出的测距系统和测距方法相比于传统的测距设备具有更高的精度,此外,SPAD的状态可以直接用二进制数字0或1直接表示,相比于APD更加节省存储空间,也就是说,使用SPAD器件数字化之后,该测距系统可以减小数据量,支持更多并行处理的通道数,支持更大规模的阵列。
此外,本申请还提供一种测距系统,如图4所示,包括发射模块401,接收模块402,和存储模块403,计算模块405。
发射模块401用于发射N个脉冲信号,该脉冲信号的脉冲间隔不是固定值,可以是随机值,例如,N可以是1024,即一次探测发射1024个脉冲,脉冲平均间隔为10ns左右,间隔分为8ns和12ns。
接收模块402包括SPAD器件,SPAD用于探测光信号,接收模块402接受回波信号,回波信号可以引起接收模块中的SPAD器件雪崩,产生雪崩信号,该信号用于记录SPAD是否被触发。
当接收模块402接受到回波信号以后,SPAD被触发,存储模块403则记录SPAD被触发的时刻,并记录和存储该时刻,形成序列B。
存储模块403,还用于存储发射出的N个脉冲信号的序列A。
计算模块405,用于根据序列A与序列B进行移动时间窗口卷积计算,并根据计算结果得到测距结果。
需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (10)

1.一种测距方法,其特征在于,包括:
发射N个脉冲信号;
存储所述N个脉冲信号及其发射时间,形成序列A;
接收回波信号,所述回波信号触发SPAD雪崩;
存储所述SPAD被触发的时刻,形成序列B;
对所述序列A与所述序列B进行移动时间窗口卷积计算,并根据所述计算结果得到测距结果。
2.如权利要求1所述的测距方法,其特征在于,所述N个脉冲信号由随机数生成。
3.如权利要求1或2所述的测距方法,其特征在于,所述N个脉冲信号的相邻脉冲的时间间隔随机。
4.如权利要求3所述的测距方法,其特征在于,所述N个脉冲的相邻脉冲的时间间隔与所述SPAD的死区时间成正比,且大于1。
5.如权利要求4所述的测距方法,其特征在于,所述N个脉冲的个数与测距范围成正比。
6.一种测距系统,其特征在于,包括:
发射模块,用于发射N个脉冲信号;
SPAD,用于探测光信号;
接收模块,用于接收回波信号,所述回波信号触发所述SPAD雪崩;
存储模块,分别存储所述N个脉冲信号及其发射时间,形成序列A,所述SPAD被触发的时刻,形成序列B;
计算模块,根据所述序列A与所述序列B进行移动时间窗口卷积计算,并根据所述计算结果得到测距结果。
7.如权利要求6所述的测距系统,其特征在于,所述N个脉冲信号由随机数生成。
8.如权利要求6或7所述的测距系统,其特征在于,所述N个脉冲信号的相邻脉冲的时间间隔随机。
9.如权利要求8所述的测距系统,其特征在于,所述N个脉冲的相邻脉冲的时间间隔与所述SPAD的死区时间成正比。
10.如权利要求8所述的测距系统,其特征在于,所述N个脉冲的个数与测距范围成正比。
CN202111260507.7A 2021-10-28 2021-10-28 一种测距方法和测距系统 Pending CN116047532A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202111260507.7A CN116047532A (zh) 2021-10-28 2021-10-28 一种测距方法和测距系统
PCT/CN2022/126423 WO2023071908A1 (zh) 2021-10-28 2022-10-20 一种测距方法和测距系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111260507.7A CN116047532A (zh) 2021-10-28 2021-10-28 一种测距方法和测距系统

Publications (1)

Publication Number Publication Date
CN116047532A true CN116047532A (zh) 2023-05-02

Family

ID=86129959

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111260507.7A Pending CN116047532A (zh) 2021-10-28 2021-10-28 一种测距方法和测距系统

Country Status (2)

Country Link
CN (1) CN116047532A (zh)
WO (1) WO2023071908A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116359935A (zh) * 2023-05-22 2023-06-30 苏州洞悉科技有限公司 一种门控成像的测距系统和测距方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3497477A1 (en) * 2016-08-12 2019-06-19 Fastree3D SA Method and device for measuring a distance to a target in a multi-user environment by means of at least one detector
US20180081041A1 (en) * 2016-09-22 2018-03-22 Apple Inc. LiDAR with irregular pulse sequence
IT201700000532A1 (it) * 2017-01-03 2018-07-03 St Microelectronics Srl Procedimento per rilevare oggetti, sistema ed apparecchiatura corrispondenti
DE102018126522A1 (de) * 2018-10-24 2020-04-30 Blickfeld GmbH Laufzeitbasierte Entfernungsmessung unter Verwendung von modulierten Pulsfolgen von Laserpulsen
JP2021096177A (ja) * 2019-12-18 2021-06-24 ソニーセミコンダクタソリューションズ株式会社 受光装置及び受光装置の制御方法、並びに、測距装置
CN111366944B (zh) * 2020-04-01 2022-06-28 浙江光珀智能科技有限公司 一种测距装置和测距方法
CN113093212A (zh) * 2021-03-30 2021-07-09 宁波飞芯电子科技有限公司 一种spad传感器与使用其的探测系统及电子设备

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116359935A (zh) * 2023-05-22 2023-06-30 苏州洞悉科技有限公司 一种门控成像的测距系统和测距方法
CN116359935B (zh) * 2023-05-22 2023-08-11 苏州洞悉科技有限公司 一种门控成像的测距系统和测距方法

Also Published As

Publication number Publication date
WO2023071908A1 (zh) 2023-05-04

Similar Documents

Publication Publication Date Title
US11754686B2 (en) Digital pixel
CN110632578B (zh) 用于时间编码时间飞行距离测量的系统及方法
US11644551B2 (en) Lidar systems with improved time-to-digital conversion circuitry
US11119196B2 (en) First photon correlated time-of-flight sensor
US10613223B2 (en) Method of detecting objects, corresponding system and apparatus
CN112859098B (zh) 一种光子数分辨测量增强单光子激光雷达系统及测距方法
US20210231782A1 (en) Dram-based lidar pixel
CN110741281B (zh) 采用迟锁盖格模式检测的LiDAR系统及方法
CN110632577B (zh) 时间编码解调处理电路及方法
US20220221562A1 (en) Methods and systems for spad optimizaiton
CN113093212A (zh) 一种spad传感器与使用其的探测系统及电子设备
US20220187430A1 (en) Time of flight calculation with inter-bin delta estimation
US20220099814A1 (en) Power-efficient direct time of flight lidar
WO2023071908A1 (zh) 一种测距方法和测距系统
CN111656219B (zh) 用于使用光信号确定至少一个对象的距离的装置和方法
CN114096879A (zh) 事件驱动的共享存储器像素
US11061137B2 (en) Proximity detection device and method comprising a pulse transmission circuit to transmit into a scene plural optical pulses with different pulse durations during a detection period
US20230196501A1 (en) Systems and Methods for Memory-Efficient Pixel Histogramming
Huseynzada et al. Innovative photodetector for LiDAR systems
US20230243928A1 (en) Overlapping sub-ranges with power stepping
WO2023150920A1 (en) Methods and apparatus for single-shot time-of-flight ranging with background light rejection
CN118112539A (zh) 测距方法、测距系统及测距装置
CN115685227A (zh) 一种探测方法
WO2023114253A1 (en) Systems and methods for memory-efficient pixel histogramming
CN114545371A (zh) 一种提高spad激光雷达精度系统及其方法

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