CN115993667A - 光学系统闭环条件下的实时大气湍流强度分布测量方法 - Google Patents
光学系统闭环条件下的实时大气湍流强度分布测量方法 Download PDFInfo
- Publication number
- CN115993667A CN115993667A CN202211702651.6A CN202211702651A CN115993667A CN 115993667 A CN115993667 A CN 115993667A CN 202211702651 A CN202211702651 A CN 202211702651A CN 115993667 A CN115993667 A CN 115993667A
- Authority
- CN
- China
- Prior art keywords
- real
- time
- atmospheric turbulence
- data
- intensity distribution
- 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
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000003287 optical effect Effects 0.000 title claims abstract description 41
- 230000003044 adaptive effect Effects 0.000 claims abstract description 19
- 238000012937 correction Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000013479 data entry Methods 0.000 claims description 3
- 238000013500 data storage Methods 0.000 claims description 3
- 238000005311 autocorrelation function Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 11
- 239000010410 layer Substances 0.000 description 9
- 238000001514 detection method Methods 0.000 description 5
- 238000000691 measurement method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000005469 granulation Methods 0.000 description 1
- 230000003179 granulation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Telescopes (AREA)
Abstract
本发明公开了一种光学系统闭环条件下的实时大气湍流强度分布测量方法。该方法包括:在大视场自适应光学系统闭环条件下实时获取系统的各部分数据和参数;实时数据传输至数据预处理程序中进行开环重建;利用重建的实时开环数据测量大气湍流的强度分布信息;输出实时的大气湍流强度分布信息。利用本发明,可以在大视场自适应光学系统闭环条件下实时获取大气湍流强度分布信息。在闭环条件下实现大气湍流强度分布的测量可以保证在大视场自适应光学系统不同条件下全天候地监测大气湍流的状态,从而更好地调整大视场自适应光学系统参数。
Description
技术领域
本发明属于大气光学技术领域,具体涉及一种光学系统闭环条件下的实时大气湍流强度分布测量方法。
背景技术
自适应光学(Adaptive Optics,AO)目前已成为大口径地基望远镜实现高分辨率成像不可或缺的技术。基于单颗导引星的传统自适应光学技术,其校正视场非常有限,严重限制了其在天文观测中的应用。从大视场成像的需求出发,科学家们提出多层共轭自适应光学技术(Multi-ConjugateAdaptive Optics,MCAO),通过对大气湍流分层探测,并利用多块变形镜(DM)共轭校正相应高度层的大气湍流,从而实现校正视场的扩大。然而,保证MCAO系统达到优异校正效果的前提之一是:输入的大气湍流分层信息是准确的。因此,必须对望远镜站址上空的大气湍流分层信息进行实时监测与准确计算。同样,对于MCAO系统的特例-地表层自适应光学系统(Ground LayerAdaptive Optics,GLAO),大气湍流的分布信息同样有助于优化导星布局及其他相应的参数。
大气光学湍流随垂直高度的分布信息常用大气折射率结构常数表征。同时该参数有着非常广泛的应用领域,例如:用于评价望远镜站址好坏;优化系统参数,包括伺服闭环带宽、波前重构算法以及MCAO系统DM的共轭高度等。此外,大气视宁度也是评价大气湍流特性的重要参数,二者关系为:
因而,通常将大气层离散为有限数量的均匀薄层,并测量每层的视宁度参数,从而得到大气折射率结构常数。目前主流的基于波前斜率探测的大气湍流强度分层测量方法主要有:SLODAR(SLOpe Detection And Ranging,SLODAR,Butterley T,Wilson R W,SarazinM.Determination of the profile of atmospheric optical turbulence strengthfrom SLODAR data[J].Monthly Notices of the Royal Astronomical Society,2006,369(2):835-845.)、S-DIMM+(S-DIMM+,Scharmer G B,VanWerkhoven T I M.S-DIMM+height characterization of day-time seeing using solar granulation[J].Astronomy&Astrophysics,2010,513:A25.)、增加分层数目的S-DIMM+方法(ZhiyongWang,Lanqiang Zhang,Lin Kong,Hua Bao,Youming Guo,Xuejun Rao,Libo Zhong,LeiZhu,Changhui Rao,A modified S-DIMM+:applying additional height grids forcharacterizing daytime seeing profiles[J],Monthly Notices of the RoyalAstronomical Society,Volume 478,Issue 2,August 2018,Pages 1459–1467,https://doi.org/10.1093/mnras/sty1097)。
上述几种方法都采用大视场自适应光学系统开环状态下采集到的大视场夏克-哈特曼波前传感器数据进行相应的大气湍流强度分布测量。需要直接采集到开环状态下的大气湍流波前数据,但是在大视场自适应光学系统闭环条件下,波前传感器只能采集到实时的波前校正残差,并不是真实的大气湍流波前数据。这给大气湍流强度分布测量造成了很大的困难,如果不能在闭环条件下测量出大气湍流的强度分布,就不能做到全时段实时监测大气湍流,不能满足大视场自适应光学系统协同调整的要求。因此,为满足大视场自适应光学系统的性能需求,目前亟需寻找一种在大视场自适应光学系统闭环条件下的实时大气湍流分层强度测量方法。
针对上述问题,本发明提出一种光学系统闭环条件下的实时大气湍流强度分布测量方法,使之可以在大视场自适应光学系统闭环的条件下,测量出大气湍流的实时强度分布,为闭环条件下的大视场自适应光学系统提供实时准确的大气湍流强度分布信息输入。
发明内容
针对大视场自适应光学系统闭环条件下大气湍流强度分布测量的需求问题,提出一种光学系统闭环条件下的实时大气湍流强度分布测量方法。
本发明采用的技术方案为:
一种光学系统闭环条件下的实时大气湍流强度分布测量方法,所述方法通过实时获取大视场自适应光学系统的闭环数据并实时重建出开环数据,最终通过大气湍流强度分布测量算法得到实时的大气湍流强度分布信息,该方法包括如下步骤:
步骤(1):在大视场自适应光学系统闭环条件下,采集系统闭环数据,包括实时波前传感器探测到的波前数据、变形镜的校正电压、以及自适应系统中的各种参数,并将实时系统闭环数据传输至实时大气湍流测量程序中;
步骤(2):利用开环数据重建方法将系统闭环数据重建为系统开环数据;
步骤(3):利用实时大气湍流强度分布测量方法测量实时的大气湍流强度分布;
步骤(4):输出大气湍流分层强度信息并实时显示。
进一步地,步骤(1)中,在大视场自适应光学系统闭环情况下,获取大视场自适应光学系统中的实时波前传感器采集的波前斜率、变形镜的实时波前校正电压、以及系统的各种参数作为所述实时系统闭环数据,然后将所述实时系统闭环数据逐帧传输给后端的实时大气湍流测量程序中以用于大气湍流强度分布的实时计算。
进一步地,步骤(2)中,所述开环数据重建方法包括基于模式法校正的畸变波前重建以及基于直接斜率法校正的畸变波前重建。
进一步地,步骤(3)中,所述大气湍流强度分布测量方法通过基于大视场夏克-哈特曼波前传感器的大气湍流强度分布测量方法实现,主要包括SLODAR方法、S-DIMM+方法、增加分层网格的S-DIMM+方法以及适应真实湍流的大气湍流强度分布测量方法。
进一步地,步骤(3)中,如图2所示,采用“栈”式数据存储方式保留一定时间的波前斜率数据:采用先进先出的窗口滑动式数据进出规则,当新的数据进入时,将在“栈”中保存最久的数据踢出。
本发明原理在于:本发明一种大视场自适应光学系统闭环条件下的实时大气湍流强度分布测量方法,所述步骤(1)中,需要实时获取系统闭环数据。大视场自适应光学系统中,波前处理机作为系统的关键,承担着根据实时波前传感器探测的波前畸变控制变形镜的重要角色。利用波前处理机就可以获得所需的所有实时系统闭环数据,然后波前处理机将实时系统闭环数据一帧一帧的传输给大气湍流测量程序中。在步骤(2)中,需要将每一帧实时闭环数据重建为开环数据。波前处理机通过各种波前控制方法包括模式法和直接斜率法将实时得到的波前畸变信息转化为变形镜上的校正电压。利用这个逆过程可以将一帧系统闭环数据重建为系统开环数据。步骤(3)中,需要缓存一定帧数的开环波前数据,并对获取的开环波前数据进行计算。大气湍流强度信息是基于大气特性的统计结果,因此需要在积累一定时间后的计算结果才更具有可信度。在本发明中具体体现在:在积累一定帧数的开环波前数据后,再计算大气湍流强度分布信息。同时,为了更好的实时性效果,本发明中采用“滑窗”式的开环波前数据选择,即在一个时间点上,当该时间点的大气湍流墙布分布计算完成后,只将波前数据中时间相对较旧的部分数据踢出,并接收最新的波前斜率数据。
本发明预计取得的有益效果为:
(1)利用本发明大视场自适应光学系统闭环条件下,实时获取大气湍流强度分布信息,能够更好地监控大气湍流强度信息;
(2)利用本发明可以在闭环条件下为大视场自适应光学系统提供调整系统参数的重要依据,提升系统校正性能;
(3)利用本发明不需要增加其他硬件条件。本发明可以集成到大视场自适应光学系统的实时处理机中,且大部分大口径地基天文望远镜均配备大视场自适应光学系统。
附图说明
图1为本发明一种光学系统闭环条件下的实时大气湍流强度分布测量方法主要的框架;
图2为本发明使用的“窗口滑动”结构示意图;
图3为本发明大视场自适应光学系统中大视场夏克-哈特曼波前传感器的子孔径排布;
图4为本发明大视场自适应光学系统中变形镜的驱动单元排布;
图5为本发明实现的大气湍流强度分布测量结果显示示意图。
具体实施方式
下面结合例图以及具体实施进一步说明本发明。
如图1所示为本发明所述一种光学系统闭环条件下的实时大气湍流强度分布测量方法的主要框架,所述方法通过实时获取大视场自适应光学系统的闭环数据并实时重建出开环数据,最终通过大气湍流强度分布测量算法得到实时的大气湍流强度分布信息,该方法包括如下步骤:
步骤(1):在大视场自适应光学系统闭环条件下,采集系统闭环数据,包括实时波前传感器探测到的波前数据、变形镜的校正电压、以及自适应系统中的各种参数,并将实时系统闭环数据传输至实时大气湍流测量程序中;
主要的数据参数包括波前斜率数据S(导星数×子孔径数×2个浮点数)、变形镜的控制电压V(驱动器单元个数个浮点数)、波前传感器子孔径的排布方式(子孔径位置信息、子孔径大小等等)、导引星区域的参数(导星位置和导星大小)。这些数据从实时的波前处理机中传输至实时的大气湍流测量程序中。
步骤(2):利用开环数据重建方法将系统闭环数据重建为系统开环数据;
利用自适应光学系统中变形镜的驱动电压和大视场波前传感器斜率之间的交互矩阵Z(该矩阵中每一列对应一个驱动器在大视场波前传感器上的响应),可以重建出变形镜闭环时被校正掉的波前信息。
SDM=ZV (2)
式中SDM表示由变形镜校正掉的波前信息,Z为交互矩阵,V为变形镜的控制电压。将变形镜校正掉的波前信息和大视场波前传感器探测到的校正残余信息相加,即可得到真实的大气湍流信息。
Spseudo=SDM+S (3)
式中Spseudo为最终重建的伪开环数据,S为由大视场波前传感器探测的残余波前信息。
步骤(3):利用实时大气湍流强度分布测量方法测量实时的大气湍流强度分布;
大气湍流的强度分布信息属于统计信息,因此需要利用容器保存一定帧数的大气湍流信息,然后才能对大气湍流进行统计测量。如图2所示,本发明利用一个先进先出(FIFO)的容器来保存一定帧数的大气湍流数据。采用“栈”式数据存储方式保留一定时间的波前斜率数据:采用先进先出的窗口滑动式数据进出规则,当新的数据进入时,将在“栈”中保存最久的数据踢出。
在获得足够帧数的伪开环数据之后,便可以使用大气湍流强度分布的测量算法对数据进行分析。此处的方法主要包括SLODAR方法和S-DIMM+方法。
步骤(4):输出大气湍流分层强度信息并实时显示。
具体实施例:采用中国科学院云南天文台1米新真空太阳望远镜研制的MCAO系统。整套MCAO系统可以视为一个地表层大视场自适应光学系统(GLAO)然后加上一个高层湍流的探测校正系统。由于大气湍流测量需要真实的大气斜率信息,所以只能从GLAO系统中获取大气湍流信息。
该GLAO系统使用9×9哈特曼波前探测器,有效子孔径个数为48个,如图3所示,由于其六边形的子孔径形状,实际中一行中只有8个子孔径。GLAO哈特曼具有约四十角秒的大视场,可以在这个视场中选择合适的导星,目前采用3×3导星排布,每个导星具有约十角秒的视场。该系统通过控制一个127单元的变形镜来实现对畸变波前的实时校正。如图4所示,127个驱动器单元在镜面上规则排布。
波前实时处理机获取波前探测器的探测斜率数据,然后根据相应的方法(模式法或者直接斜率法)计算得到变形镜的电压。在系统闭环条件下,波前探测器只能探测到校正后的参与波前信息,这并不能用于大气湍流的强度分布测量。在闭环条件下,可以获取相应的波前探测器数据和变形镜电压数据,并结合波前控制方法,重建出相应的开环波前信息。每一帧重建的开环数据都会被传输至大气湍流强度分布测量程序中。本实施例设置每秒重建的开环数据的帧数为200帧。
在本实施例的大气湍流测量程序中,设置了一个容纳2000帧数据的栈空间。当栈被填满时,会进行一次计算,计算得到相应的大气湍流强度分布结果。完成一次计算之后,栈会弹出200帧数据,然后继续接收200帧数据用于下一次计算。每一次计算之后,程序会显示出相应的测量结果,如图5所示,图5中展示了利用本发明的方法对云南天文台1米新真空太阳望远镜站址所在大气湍流进行测量的一组实时结果。测量结果为大气湍流的归一化强度分布(Turbulence Intensity Distribution),横坐标为测量的高度(Height:单位为米)层。标题处显示了测量的大气相干长度(Total Seeing表示大气湍流的总强度),单位为米。
本实施例每次计算的实时计算结果可以用于实时调整大视场自适应光学系统的相应参数,便于达到更好的校正效果,实现衍射极限分辨率的大视场天文观测,也可以将相应的结果记录或者传输至显示系统进行展示。
以上所述,仅为本发明的具体实施方式,未详细阐述的部分属于本领域公知技术。但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解到的替换或增减,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (5)
1.一种光学系统闭环条件下的实时大气湍流强度分布测量方法,其特征在于,所述方法通过实时获取大视场自适应光学系统闭环条件下的系统数据,通过重建开环获取真实的大气数据,然后利用各种大气湍流强度分层测量方法,得到实时的大气湍流强度分布信息,该方法包括如下步骤:
步骤(1):在大视场自适应光学系统闭环条件下,采集系统闭环数据,包括实时波前传感器探测到的波前数据、变形镜的校正电压、以及自适应系统中的各种参数,并将实时系统闭环数据传输至实时大气湍流测量程序中;
步骤(2):利用开环数据重建方法将系统闭环数据重建为系统开环数据;
步骤(3):利用实时大气湍流强度分布测量方法测量实时的大气湍流强度分布;
步骤(4):输出大气湍流分层强度信息并实时显示。
2.根据权利要求1所述的一种光学系统闭环条件下的实时大气湍流强度分布测量方法,其特征在于,步骤(1)中,在大视场自适应光学系统闭环情况下,获取大视场自适应光学系统中的实时波前传感器采集的波前斜率、变形镜的实时波前校正电压、以及系统的各种参数作为所述实时系统闭环数据,然后将所述实时系统闭环数据逐帧传输给后端的实时大气湍流测量程序中以用于大气湍流强度分布的实时计算。
3.根据权利要求1所述的一种光学系统闭环条件下的实时大气湍流强度分布测量方法,其特征在于,步骤(2)中,所述开环数据重建方法包括基于模式法校正的畸变波前重建以及基于直接斜率法校正的畸变波前重建。
4.根据权利要求1所述的一种光学系统闭环条件下的实时大气湍流强度分布测量方法,其特征在于,步骤(3)中,所述大气湍流强度分布测量方法通过基于大视场夏克-哈特曼波前传感器的大气湍流强度分布测量方法实现,包括SLODAR方法、S-DIMM+方法、增加分层网格的S-DIMM+方法以及基于数据统计自相关函数的大气湍流强度分布测量方法。
5.根据权利要求4所述的一种光学系统闭环条件下的实时大气湍流强度分布测量方法,其特征在于,步骤(3)中,采用“栈”式数据存储方式保留一定时间的波前斜率数据:采用先进先出的窗口滑动式数据进出规则,当新的数据进入时,将在“栈”中保存最久的数据踢出。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211685624 | 2022-12-27 | ||
CN2022116856242 | 2022-12-27 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115993667A true CN115993667A (zh) | 2023-04-21 |
Family
ID=85989937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211702651.6A Pending CN115993667A (zh) | 2022-12-27 | 2022-12-29 | 光学系统闭环条件下的实时大气湍流强度分布测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115993667A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115016035A (zh) * | 2022-05-31 | 2022-09-06 | 中国科学院光电技术研究所 | 一种基于波前探测的实时大气湍流分层强度测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353514A (zh) * | 2011-08-25 | 2012-02-15 | 中国科学院光电技术研究所 | 一种用于实时大气湍流参数计算的信号处理系统 |
CN104391466A (zh) * | 2014-11-24 | 2015-03-04 | 中国科学院自动化研究所 | 一种快速自适应光学控制器的设计方法 |
CN115016035A (zh) * | 2022-05-31 | 2022-09-06 | 中国科学院光电技术研究所 | 一种基于波前探测的实时大气湍流分层强度测量方法 |
-
2022
- 2022-12-29 CN CN202211702651.6A patent/CN115993667A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353514A (zh) * | 2011-08-25 | 2012-02-15 | 中国科学院光电技术研究所 | 一种用于实时大气湍流参数计算的信号处理系统 |
CN104391466A (zh) * | 2014-11-24 | 2015-03-04 | 中国科学院自动化研究所 | 一种快速自适应光学控制器的设计方法 |
CN115016035A (zh) * | 2022-05-31 | 2022-09-06 | 中国科学院光电技术研究所 | 一种基于波前探测的实时大气湍流分层强度测量方法 |
Non-Patent Citations (2)
Title |
---|
吴碧琳;饶长辉;: "自适应光学系统中高精跟踪回路性能在线分析", 光电工程, no. 06, 28 June 2006 (2006-06-28), pages 1 - 3 * |
张昂;饶长辉;张雨东;姜文汉;: "基于自适应光学系统闭环数据的大气相干长度的计算", 光电工程, no. 12, 28 December 2006 (2006-12-28), pages 19 - 22 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115016035A (zh) * | 2022-05-31 | 2022-09-06 | 中国科学院光电技术研究所 | 一种基于波前探测的实时大气湍流分层强度测量方法 |
CN115016035B (zh) * | 2022-05-31 | 2023-12-22 | 中国科学院光电技术研究所 | 一种基于波前探测的实时大气湍流分层强度测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6570991B2 (ja) | 大口径望遠鏡における、非同一平面上の(anisoplanatic)画像形成のためのレンズレット、ビームウォーク(beamwalk)、及びチルトの多様化 | |
CN115993667A (zh) | 光学系统闭环条件下的实时大气湍流强度分布测量方法 | |
CN102889935A (zh) | 基于相位差法的自适应光学系统近场波前传感器标定装置及标定方法 | |
CN102879901B (zh) | 基于太阳地表层共轭自适应光学系统的大视场高分辨力成像装置 | |
Miklenic et al. | Coronal dimmings and the early phase of a CME observed with STEREO and Hinode/EIS | |
CN107942406B (zh) | 一种基于总大气相干长度约束的分层大气湍流强度测量方法 | |
CN115016035B (zh) | 一种基于波前探测的实时大气湍流分层强度测量方法 | |
CN104486562A (zh) | 基于固定积分时间的嵌入式红外图像超帧处理方法 | |
CN102967379A (zh) | 一种用于太阳自适应光学系统的波前传感器 | |
Rimmele et al. | Solar adaptive optics: a progress report | |
US6563572B1 (en) | Correlating shack-hartmann wavefront sensor | |
Ichimoto et al. | Measurement of the coronal electron temperature at the total solar eclipse on 1994 November 3 | |
Padovani et al. | The ESO's Extremely Large Telescope Working Groups | |
Goode et al. | Measuring seeing from solar scintillometry and the spectral ratio technique | |
Ziad et al. | First statistics of the turbulence outer scale at Dome C | |
Hart et al. | Multiframe blind deconvolution for imaging in daylight and strong turbulence conditions | |
CN108362697A (zh) | 一种增加分层数目的大气视宁度分层测量方法 | |
Schreiber et al. | The laser guide stars wavefront sensor prototype for the E-ELT: test results | |
Agapito et al. | NGSs acquisition in MORFEO | |
Wilson et al. | Turbulence profiler and seeing monitor for laser guide star adaptive optics | |
Wilson et al. | Development of a portable SLODAR turbulence profiler | |
Ziad et al. | The first optical characterization of the Oukaïmeden site with the Generalized Seeing Monitor (GSM) | |
Peter | Point-spread function reconstruction for the ground layer adaptive optics system ARGOS | |
Lehtonen et al. | Real-time turbulence profiling using particle filtering | |
Hart et al. | Quantitative measurements of daytime near infrared sky brightness at the AEOS 3.6 m telescope |
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 |