CN114494704A - 二值化图像抗噪声提取骨架的方法及系统 - Google Patents
二值化图像抗噪声提取骨架的方法及系统 Download PDFInfo
- Publication number
- CN114494704A CN114494704A CN202210140719.XA CN202210140719A CN114494704A CN 114494704 A CN114494704 A CN 114494704A CN 202210140719 A CN202210140719 A CN 202210140719A CN 114494704 A CN114494704 A CN 114494704A
- Authority
- CN
- China
- Prior art keywords
- pixel
- skeleton
- binary image
- noise
- sample
- 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 49
- 239000013598 vector Substances 0.000 claims description 40
- 238000006073 displacement reaction Methods 0.000 claims description 28
- 238000000605 extraction Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 238000000513 principal component analysis Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 238000007670 refining Methods 0.000 claims description 6
- 210000000988 bone and bone Anatomy 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000001133 acceleration Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Processing (AREA)
Abstract
本发明属于图像骨架化技术领域,具体公开了一种二值化图像抗噪声提取骨架的方法及系统,该方法利用二值图像的前景像素,计算每个前景像素在其搜索域中的向量之和,得到其初始位移矢量,进行主成分分析,对初始位移矢量进行优化,沿优化后的位移矢量移动每个像素点,并根据移动点的集合重新绘制图像,重复步骤,直到满足终止条件时停止,得到原始骨架,对原始骨架进行细化,删除不改变图连通性的非必要像素点,得到最终的骨架。采用本技术方案,收集更多的局部信息,使提取的骨架更加平滑,而且对噪声具有很强的鲁棒性。
Description
技术领域
本发明属于图像骨架化技术领域,涉及一种二值化图像抗噪声提取骨架的方法及系统。
背景技术
骨架化是将二值图像上的前景区域减少为骨架残余的过程,它在很大程度上保持了原始区域的连通性,从而提供了一种实用和紧凑的图像对象表示方式。提取的骨架可以作为一种有用的形状描述工具,也是一种著名的图像分析和形状识别的预处理方法。
现有技术中有许多在二维空间中获取骨骼的算法,可以将它们分类如下:
基于细化的骨架化方法:以直观的方式对形状进行剥离,得到一组单像素宽度的连通点,保持了形状的拓扑结构。即细化是一种旨在以并行或顺序的方式删除非终结简单点的操作,这些算法的主要优点是保留了形状拓扑。
基于距离图的骨架化方法:目标是识别距离图上的关键点,其中每个像素都用值标记与最接近的背景像素的距离。不同的距离映射近似或精确计算欧几里得距离。例如,倒角(通过局部掩模逼近欧几里得距离),平方欧几里得距距离,带符号的欧几里得距离和蜂窝(基于六边形网格)。
而现有技术基于固定4邻域和8邻域的算法,只能考虑相当局部的信息,特别是当噪声接近边界时,局部信息的丢失将导致生成的骨架质量严重退化,算法鲁棒性不足。
发明内容
本发明的目的在于提供一种二值化图像抗噪声提取骨架的方法及系统,以解决算法鲁棒性不足的问题。
为了达到上述目的,本发明的基础方案为:一种二值化图像抗噪声提取骨架的方法,包括如下步骤:
S1,获取二值图像的前景像素,计算每个前景像素在其搜索域中的向量之和,得到其初始位移矢量;
S2,进行主成分分析,对初始位移矢量进行优化;
S3,沿优化后的位移矢量移动每个像素点,并根据移动点的集合重新绘制图像;
S4,重复步骤S1-S3,直到满足终止条件时停止,得到原始骨架;
S5,对原始骨架进行细化,删除不改变图连通性的非必要像素点,得到最终的骨架。
本基础方案的工作原理和有益效果在于:本方案计算每个前景像素在其搜索域中的向量之和,得到其初始运动方向,然后对之前计算的向量进行优化。沿着优化后的位移矢量移动每个像素点,并根据移动点的集合重新绘制图像。直到满足终止条件,进行后处理,得到最终的骨架。这样基于最近邻的抗噪声骨架化算法,收集了更多的局部信息,使提取的骨架不仅更加平滑,而且对噪声具有很强的鲁棒性。
进一步,得到初始位移矢量的方法如下:
设二值图像S用大小为m×n的矩阵P表示;
给定一组前景像素,通过计算每个样本像素周围向量的和,得到初始位移向量ui为:
其中,i为待计算样本像素标识,j为i的邻域内样本像素标识,xi为待计算样本像素坐标,xj为邻域内样本像素坐标,ki为i的自适应的邻域内所含像素点个数;Ni,k是样本像素xi的k近邻的集合。
计算简单,便于操作。
进一步,对初始位移矢量进行优化的方法如下:
对每个前景像素在其Ni,k上执行主成分分析,得到一个对称的2×2半正定矩阵M:
为了明确区分哪些像素应该移动,哪些像素不应该移动,添加一个高斯权值来调整最终的位移矢量,得到优化后的位移矢量如下:
利用优化的位移矢量,确定需要移动的像素和不需要移动的像素,避免误移动,利于使用。
设置高斯权值,调整最终的位移矢量,利于使用。
进一步,为了能够提取不同尺度图像的骨架,需要自适应更新最近邻的个数,即k,计算如下:
式中,Abb为原始图像中目标的包围框面积,c是比例因子。
基于k近邻进行骨架提取,可在不同的尺度上提取准确、平滑的图像,便于使用。
进一步,对原始骨架进行细化的方法如下:
基于以上设定,过滤掉中心轴端点;
对于结果端点集,对于每个连通域,计算其质心,并将该质心视为实际端点;
对于前景像素,越靠近边界,|ui|越大,首先按照|ui|的值降序对像素进行排序,然后依次删除它们,确保总是先删除最外层的像素。
细化骨架,同时保证骨架质量,利于使用。
进一步,细化过程中减少计算量的方法如下:
当被计算的ui等于0,即样本像素的邻域是一个完全对称区域时,将其Δxi设为(0,0)。
减少计算量,加快计算速度,利于使用。
进一步,细化过程中减少计算量的方法如下:
基于高斯圆问题:
其中,r为圆的半径,N(r)为该半径下邻域中所含像素点个数;当半径达到一定的阈值时,在数组的下一次迭代中提前保存可能的邻居点的相对位置,在下一次迭代中,根据阵列中的相对位置直接计算出点的绝对位置。
在最终的程序中实施加速策略,完成计算加速,便于使用。
本发明还提供一种二值化图像抗噪声提取骨架的系统,包括图像获取模块和处理模块,所述图像获取模块用于获取二值图像,图像获取模块的输出端与处理模块连接,所述处理模块执行本发明所述方法,进行骨架提取。
该系统利用图像获取模块获取所需图像,并通过处理模块完成骨架提取,操作简单,利于使用。
附图说明
图1是本发明二值化图像抗噪声提取骨架的方法的流程示意图;
图2是本发明二值化图像抗噪声提取骨架的方法的样本及其邻域的坐标示意图;
图3是本发明二值化图像抗噪声提取骨架的方法的原始骨架的结构示意图;
图4是本发明二值化图像抗噪声提取骨架的方法的不均匀变薄的骨架的结构示意图;
图5是本发明二值化图像抗噪声提取骨架的方法的细化骨架的结构示意图;
图6是本发明二值化图像抗噪声提取骨架的方法的一种优选方案的边界像素的结构示意图;
图7是本发明二值化图像抗噪声提取骨架的方法的一种优选方案的边界像素的结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明的描述中,除非另有规定和限定,需要说明的是,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是机械连接或电连接,也可以是两个元件内部的连通,可以是直接相连,也可以通过中间媒介间接相连,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
本发明公开了一种二值化图像抗噪声提取骨架的方法,用于解决现有技术中基于固定4邻域和8邻域的算法鲁棒性不足的问题。如图1所示,本发明的二值化图像抗噪声提取骨架方法包括如下步骤:
S1,获取二值图像的前景像素,计算每个前景像素在其搜索域中的向量之和,得到其初始位移矢量;
S2,进行主成分分析,对初始位移矢量进行优化;
S3,沿优化后的位移矢量移动每个像素点,并根据移动点的集合重新绘制图像;
S4,重复步骤S1-S3,直到满足终止条件时停止,得到原始骨架;
S5,对原始骨架进行细化,删除不改变图连通性的非必要像素点,得到最终的骨架。
本发明的一种优选方案中,得到初始位移矢量的方法如下:
设二值图像S用大小为m×n的矩阵P表示;
给定一组前景像素,通过计算每个样本像素周围向量的和,得到初始位移向量ui为:
其中,i为待计算样本像素标识,j为i的邻域内样本像素标识,xi为待计算样本像素坐标,xj为邻域内样本像素坐标,ki为i的自适应的邻域内所含像素点个数;Ni,k是样本像素xi的k近邻的集合。位移向量ui决定了移动的主要方向,但当面对一个具有特殊局部信息的像素时,如图2中所示右边的样本像素,这个像素处于极端位置,它不应该被移动,但它的ui在栅格化后被计算为(0,-1)。因此,在移动前,需要对每个前景像素在其Ni,k上执行主成分分析(PCA)。
本发明的一种优选方案中,对初始位移矢量进行优化的方法如下:
对每个前景像素在其Ni,k上执行主成分分析,得到一个对称的2×2半正定矩阵M:
为了明确区分哪些像素应该移动,哪些像素不应该移动,添加一个高斯权值来调整最终的位移矢量,得到优化后的位移矢量如下:
其中,表示跳跃函数;细节因子γ表示提取骨架的细节程度,γ值越高,骨架的细节就越多;符号a是比例因子,默认值为10;表示样本像素xi的主方向;e为自然底数。例如,基于8个样本像(如图6所示)素的σi构造一个跳跃函数,σi本身包含了对对象方向不敏感的特征向量之间的比率信息。具体说明如下:
J中的偏移量0.7723为第一个矩形和第二个矩形边界像素的σ,即0.7875和0.7571的均值。通过这种改进,可以使边界正常移动,改进结果如图7所示。
所有线性测量的平均值σi表示为表示骨骼的线性程度;随着萃取的进行,将趋于1;当连续三次不变时,满足终止条件,停止提取。在每次迭代结束时,由于移动会有许多重复的点,这些点可以通过简单地按位置重新绘制来删除。
本发明的一种优选方案中,为了能够提取不同尺度图像的骨架,需要自适应更新最近邻的个数,即k,计算如下:
式中,Abb为原始图像中目标的包围框面积,c是比例因子,默认值为10。利用该更新方案,能够在不同的尺度上提取准确、平滑的图像。
经过上述提取,得到一个较厚但拓扑保存完好的原始骨架,原始骨架需要进一步减薄。其基本思想是,如果一个像素被删除,而不改变图的连通性,那么这个像素实际上是可以被删除的。但是端点应该进行预先保留,因为如果沿着端点移除,虽然可以保持图的连接性,但会导致整个骨架消失。因此,本发明的一种优选方案中,对原始骨架进行细化的方法如下:
根据这两个事实,过滤掉中心轴端点;
结果端点集可能是连接域的集合,对于每个连通域,计算其质心,并将该质心视为实际端点;然后对其进行稀释:由于原始骨架不是严格光滑的,如图3所示。如果简单地从一个方向或多个方向删除像素,最终会得到一个不均匀变薄的骨架,如图4所示。
根据ui长度的性质,即|ui|,对于前景像素,它越靠近边界,它的|ui|越大;
因此,首先按照|ui|的值降序对像素进行排序,然后依次删除它们,确保总是先删除最外层的像素,得到结果如图5所示。
细化过程中减少计算量的优选方法如下:
当被计算的ui等于0,即样本像素的邻域是一个完全对称区域时,将其Δxi设为(0,0)。
细化过程中减少计算量的另一种优选方法如下:
基于高斯圆问题:
其中,r为圆的半径,N(r)为该半径下邻域中所含像素点个数;当半径达到一定的阈值时,可能会出现新的邻居,例如,当半径在2到2.3之间时,不会出现新的邻居,只有当半径达到2.3时,才会导致邻居数量的可能变化。这揭示了一种加速策略,在这种策略中,不需要一点点地增加半径来搜索邻居,在数组的下一次迭代中提前保存可能的邻居点的相对位置,在下一次迭代中,根据阵列中的相对位置直接计算出点的绝对位置。随着半径r增大,所包含的潜在邻居点越来越多,且其坐标相对固定。如r在(0,1)内没有新增,在r=1时新增相对位置(0,1)(1,0)(0,-1)(-1,0)四个邻居点。在最终的程序中实施这个加速策略,可带来4倍的加速。
本发明还提供一种二值化图像抗噪声提取骨架的系统,包括图像获取模块和处理模块,图像获取模块用于获取二值图像,图像获取模块的输出端与处理模块电性连接,所述处理模块执行本发明所述方法,进行骨架提取。本方案基于最近邻的抗噪声骨架化算法,比4邻域或8邻域收集了更多的局部信息,提取的骨架不仅更加平滑,而且对噪声具有很强的鲁棒性。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。
Claims (9)
1.一种二值化图像抗噪声提取骨架的方法,其特征在于,包括如下步骤:
S1,获取二值图像的前景像素,计算每个前景像素在其搜索域中的向量之和,得到其初始位移矢量;
S2,进行主成分分析,对初始位移矢量进行优化;
S3,沿优化后的位移矢量移动每个像素点,并根据移动点的集合重新绘制图像;
S4,重复步骤S1-S3,直到满足终止条件时停止,得到原始骨架;
S5,对原始骨架进行细化,删除不改变图连通性的非必要像素点,得到最终的骨架。
7.如权利要求1所述的二值化图像抗噪声提取骨架的方法,其特征在于,细化过程中减少计算量的方法如下:
当被计算的ui等于0,即样本像素的邻域是一个完全对称区域时,将其Δxi设为(0,0)。
9.一种二值化图像抗噪声提取骨架的系统,其特征在于,包括图像获取模块和处理模块,所述图像获取模块用于获取二值图像,图像获取模块的输出端与处理模块连接,所述处理模块执行权利要求1-8之一所述方法,进行骨架提取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210140719.XA CN114494704A (zh) | 2022-02-16 | 2022-02-16 | 二值化图像抗噪声提取骨架的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210140719.XA CN114494704A (zh) | 2022-02-16 | 2022-02-16 | 二值化图像抗噪声提取骨架的方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114494704A true CN114494704A (zh) | 2022-05-13 |
Family
ID=81481313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210140719.XA Pending CN114494704A (zh) | 2022-02-16 | 2022-02-16 | 二值化图像抗噪声提取骨架的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114494704A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115311507A (zh) * | 2022-10-12 | 2022-11-08 | 南通宏梁建筑科技有限公司 | 基于数据处理的建筑板材分类方法 |
CN117218091A (zh) * | 2023-09-19 | 2023-12-12 | 徐州医科大学 | 面向骨折地图构建的骨折线提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156730A (zh) * | 2014-07-25 | 2014-11-19 | 山东大学 | 一种基于骨架的抗噪声汉字特征提取方法 |
US20160328855A1 (en) * | 2015-05-04 | 2016-11-10 | Siemens Aktiengesellschaft | Method and System for Whole Body Bone Removal and Vascular Visualization in Medical Image Data |
CN107330901A (zh) * | 2017-06-29 | 2017-11-07 | 西安理工大学 | 一种基于骨架的物体构件分解方法 |
CN110992304A (zh) * | 2019-10-30 | 2020-04-10 | 浙江力邦合信智能制动系统股份有限公司 | 二维图像深度测量方法及其在车辆安全监测中的应用 |
CN111968089A (zh) * | 2020-08-15 | 2020-11-20 | 晋江市博感电子科技有限公司 | 一种基于最大内切球机制的l1中值骨架提取方法 |
-
2022
- 2022-02-16 CN CN202210140719.XA patent/CN114494704A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156730A (zh) * | 2014-07-25 | 2014-11-19 | 山东大学 | 一种基于骨架的抗噪声汉字特征提取方法 |
US20160328855A1 (en) * | 2015-05-04 | 2016-11-10 | Siemens Aktiengesellschaft | Method and System for Whole Body Bone Removal and Vascular Visualization in Medical Image Data |
CN107330901A (zh) * | 2017-06-29 | 2017-11-07 | 西安理工大学 | 一种基于骨架的物体构件分解方法 |
CN110992304A (zh) * | 2019-10-30 | 2020-04-10 | 浙江力邦合信智能制动系统股份有限公司 | 二维图像深度测量方法及其在车辆安全监测中的应用 |
CN111968089A (zh) * | 2020-08-15 | 2020-11-20 | 晋江市博感电子科技有限公司 | 一种基于最大内切球机制的l1中值骨架提取方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115311507A (zh) * | 2022-10-12 | 2022-11-08 | 南通宏梁建筑科技有限公司 | 基于数据处理的建筑板材分类方法 |
CN117218091A (zh) * | 2023-09-19 | 2023-12-12 | 徐州医科大学 | 面向骨折地图构建的骨折线提取方法 |
CN117218091B (zh) * | 2023-09-19 | 2024-03-29 | 徐州医科大学 | 面向骨折地图构建的骨折线提取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10127675B2 (en) | Edge-based local adaptive thresholding system and methods for foreground detection | |
CN107038717B (zh) | 一种基于立体栅格自动分析3d点云配准误差的方法 | |
CN106981068B (zh) | 一种联合像素与超像素的交互式图像分割方法 | |
CN108537239B (zh) | 一种图像显著性目标检测的方法 | |
CN111340824B (zh) | 一种基于数据挖掘的图像特征分割方法 | |
CN101916446A (zh) | 基于边缘信息和均值移位的灰度目标跟踪算法 | |
CN110992381A (zh) | 一种基于改进Vibe+算法的运动目标背景分割方法 | |
CN114494704A (zh) | 二值化图像抗噪声提取骨架的方法及系统 | |
CN108537816B (zh) | 一种基于超像素和背景连接先验的显著物体分割方法 | |
CN106294705A (zh) | 一种批量遥感影像预处理方法 | |
CN111583279A (zh) | 一种基于pcba的超像素图像分割方法 | |
CN109741358B (zh) | 基于自适应超图学习的超像素分割方法 | |
CN114818851B (zh) | 一种基于特征差异性的农村道路矢量数据纠正方法 | |
CN113409332A (zh) | 一种基于三维点云的建筑物平面分割方法 | |
CN105405152B (zh) | 基于结构化支持向量机的自适应尺度目标跟踪方法 | |
CN112766340A (zh) | 基于自适应空间模式的深度胶囊网络图像分类方法及系统 | |
CN114842308B (zh) | 一种基于全特征融合建立目标预判优化模型的方法 | |
CN111091071B (zh) | 基于探地雷达双曲波拟合的地下目标检测方法及系统 | |
CN115471682A (zh) | 一种基于SIFT融合ResNet50的图像匹配方法 | |
CN106815851B (zh) | 一种基于视觉测量的栅格圆形油位计自动读数方法 | |
CN114862765A (zh) | 一种细胞分层图像处理方法 | |
CN109829511B (zh) | 基于纹理分类的下视红外图像中云层区域检测方法 | |
CN104881669A (zh) | 一种基于颜色对比度的局部区域检测子提取方法及系统 | |
CN113379802A (zh) | 一种多特征自适应融合的相关滤波目标跟踪方法 | |
CN111724409A (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 |