CN112581521B - 一种磁浮轨道中心线提取方法 - Google Patents
一种磁浮轨道中心线提取方法 Download PDFInfo
- Publication number
- CN112581521B CN112581521B CN202011536006.2A CN202011536006A CN112581521B CN 112581521 B CN112581521 B CN 112581521B CN 202011536006 A CN202011536006 A CN 202011536006A CN 112581521 B CN112581521 B CN 112581521B
- Authority
- CN
- China
- Prior art keywords
- track
- point
- point cloud
- model
- magnetic levitation
- 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
- 239000000725 suspension Substances 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000003044 adaptive effect Effects 0.000 claims abstract description 19
- 238000013519 translation Methods 0.000 claims abstract description 15
- 238000005339 levitation Methods 0.000 claims description 50
- 239000013598 vector Substances 0.000 claims description 32
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000013461 design Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000011960 computer-aided design Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000011513 prestressed concrete Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 210000003781 tooth socket Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/68—Analysis of geometric attributes of symmetry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/12—Geometric CAD characterised by design entry means specially adapted for CAD, e.g. graphical user interfaces [GUI] specially adapted for CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
- G06V10/751—Comparing pixel values or logical combinations thereof, or feature values having positional relevance, e.g. template matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Architecture (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Human Computer Interaction (AREA)
- Artificial Intelligence (AREA)
- Multimedia (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Vehicles With Linear Motors And Vehicles That Are Magnetically Levitated (AREA)
Abstract
本发明涉及一种磁浮轨道中心线提取方法,该方法包括以下步骤:获取轨道断面点云数据和磁浮轨道断面模型点集,磁浮轨道断面模型点集包括轨道模型平面中心点;从轨道断面点云数据中提取点云轨道断面初始匹配中心点;依据轨道模型平面中心点至点云轨道断面初始匹配中心点所需的平移参数平移磁浮轨道断面模型点集,基于距离最小化原则,利用轨道断面点云数据重构磁浮轨道断面模型点集,得到适应性模型点集;利用最小二乘改进的二维ICP算法计算适应性模型点集与轨道断面点云数据的匹配参数,得到最终点云轨道中心点;基于最终点云轨道中心点,利用拟合算法得到磁浮轨道的中心线。与现有技术相比,可有效提高磁浮轨道中心线提取的效率和精度。
Description
技术领域
本发明涉及磁浮轨道测量领域,尤其是涉及一种磁浮轨道中心线提取方法。
背景技术
高速磁悬浮列车是一种利用电磁力实现车辆悬浮、导向,直线电机实现牵引驱动的新型地面交通工具,其轨道由预应力混凝土梁、连接件、功能件以及具有齿槽结构分段拼接的长定子组成,列车环抱滑行面、导向面、定子面“飞行”。高速列车的安全、平稳、高效运行根本在于铁路轨道是否满足设计的型和位,而轨道的空间位置、几何状态以及线路附属物需要定期的测量和维护才能准确地反应出来。磁浮轨道作为磁浮基础设施对象测绘定位中的关键要素,由于其轨道呈面状且路线长等特性,目前现有方法获得的数据比较离散且低效,无法准确、快速反映轨道各要素,难以满足现代化铁路运营维护的监测要求。
传统的轨道中心线测量主要通过以下几种方式:1)人工上线测量,钢尺丈量里程,用全站仪矢距法或偏角法对曲线测设,使用全站仪、GPS RTK的方法直接采取各要素坐标。然而这些方法由于需要人工上线安全性差,并且工作量大,需要大量的人力物力造成了资源的损失,同时存在着收集信息不全面、作业效率低等方面的缺点。2)使用高分辨率影像,基于栅格图像处理方法能够有效地提取铁路区域,减少铁路轨道线虚警目标的出现,微分几何算法实现准确的铁路轨道线提取,但目前提取率未达到百分百,测量精度相对传统方式较低。3)基于机载LiDAR的铁路轨道提取方法,利用机载LiDAR和影像数据,建立了一种面向对象的铁路区域提取算法,通过数学形态学处理消除干扰因素的影响,获取铁轨区域的LiDAR点云数据,然后根据铁路轨道点云高度特征进一步分离轨道点云,最后通过最小二乘法对轨道点云数据进行精确拟合,获取轨道信息。但目前该方法存在精度不能满足工程需求且在轨道曲线方向上的识别效果较差。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种磁浮轨道中心线提取方法。
本发明的目的可以通过以下技术方案来实现:
一种磁浮轨道中心线提取方法,该方法包括以下步骤:
S1:获取轨道断面点云数据和磁浮轨道断面模型点集,所述磁浮轨道断面模型点集包括轨道模型平面中心点;
S2:从轨道断面点云数据中提取点云轨道断面初始匹配中心点;
S3:依据轨道模型平面中心点至点云轨道断面初始匹配中心点所需的平移参数平移磁浮轨道断面模型点集,基于距离最小化原则,利用轨道断面点云数据重构磁浮轨道断面模型点集,得到适应性模型点集;
S4:利用最小二乘改进的二维ICP算法计算适应性模型点集与轨道断面点云数据的匹配参数,得到最终点云轨道中心点;
S5:基于最终点云轨道中心点,利用拟合算法得到磁浮轨道的中心线。
S1中,获取轨道断面点云数据和磁浮轨道断面模型点集包括以下步骤:
获取磁浮轨道扫描点云数据,基于磁浮轨道设计参数得到磁浮轨道断面模型点集;
基于移动扫描轨迹数据,对磁浮轨道扫描点云数据利用矢量投影进行切片,生成切片点云数据;
对切片点云数据进行投影、旋转和平移生成轨道断面点云数据。
所述轨道断面点云数据表示为:
基于磁浮轨道设计参数得到磁浮轨道断面模型点集的过程包括:
根据磁浮轨道设计参数绘制磁浮轨道断面CAD图;
在磁浮轨道断面CAD图中将磁浮轨道断面线段定距离等分,得到磁浮轨道断面模型点集。
S2中,提取点云轨道断面初始匹配中心点的过程包括:
从轨道断面点云数据中提取大于高度阈值的数据,得到轨道平面点集;
计算轨道平面点集的均值,得到点云轨道断面初始匹配中心点。
S3中得到适应性模型点集的过程包括:
在平移后的磁浮轨道断面模型点集中提取距离轨道断面点云数据中每一个点云数据点最近的模型点,重构磁浮轨道断面模型点集,再利用平移参数反向平移重构的磁浮轨道断面模型点集,得到适应性模型点集。
S4中,最小二乘改进的二维ICP算法的目标函数为:
其中,i为第i个数据点,x′i,y′i为参考数据点,b1和b2为平移参数,α为旋转矩阵参数,n为数据点个数。
S4中,计算匹配参数α,b1和b2的公式为:
其中,X′表示为:
Y′表示为:
X表示为:
Y表示为:
S4中,得到匹配参数α,b1和b2后,判断b1和b2与前一次迭代的匹配参数的差值是否大于等于阈值,若是,以b1和b2作为平移参数,执行S3,若否,计算最终点云轨道中心点最终点云轨道中心点的计算公式为:
与现有技术相比,本发明具有以下优点:
(1)提出应用模板点云匹配进行磁浮轨道的中心线提取,通过对轨道断面点云数据的匹配搜寻,可获取密集的点云轨道平面中心点,利用拟合的方式拟合磁浮轨道的中心线方程,克服了人工检测的工作强度大、效率低等缺点,对磁浮的安全运营维护提供数据基础。
(2)在实际扫描的磁浮轨道扫描点云数据中,因扫描仪的位置和扫描线是等角分布的原因,得到的轨道断面点云数据密度随位置分布而产生不均匀现象,若直接采用分布密度均匀的模型点集,其匹配效果必然受限,本发明中模型点集采用了基于距离最小化的适应性密度建立适应性模型点集,基于迭代和距离最小原则,不断对模型点集进行适应性密度重构,适应于轨道断面点云分布不均匀的场景,可有效提高模板点云匹配的精度。
附图说明
图1为本发明的流程图;
图2为本发明的利用矢量投影进行切片示意图;
图3为本发明的轨道断面点云数据示意图;
图4为本发明的磁浮轨道结构示意图;
图5为本发明的磁浮轨道断面模型点集示意图;
图6为本发明的适应性模型点集与轨道断面点云数据匹配结果图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例
本实施例提供一种基于模板点云匹配的磁浮轨道中心线提取方法,如图1所示,包括以下步骤:
S1:获取轨道断面点云数据和磁浮轨道断面模型点集,所述磁浮轨道断面模型点集包括轨道模型平面中心点;
S2:从轨道断面点云数据中提取点云轨道断面初始匹配中心点;
S3:依据轨道模型平面中心点至点云轨道断面初始匹配中心点所需的平移参数平移磁浮轨道断面模型点集,基于距离最小化原则,利用轨道断面点云数据重构磁浮轨道断面模型点集,得到适应性模型点集;
S4:利用最小二乘改进的二维ICP算法计算适应性模型点集与轨道断面点云数据的匹配参数,得到最终点云轨道中心点;
S5:基于最终点云轨道中心点,利用拟合算法得到磁浮轨道的中心线。
与现有技术相比,基于模板点云匹配的磁浮轨道中心线提取方法具有定位精确、效率高、安全性高等特点。
具体而言:
1)S1中,获取轨道断面点云数据和磁浮轨道断面模型点集包括以下步骤:
获取磁浮轨道扫描点云数据,基于磁浮轨道设计参数得到磁浮轨道断面模型点集;基于移动扫描轨迹数据,对磁浮轨道扫描点云数据利用矢量投影进行切片,生成切片点云数据;对切片点云数据进行投影、旋转和平移生成轨道断面点云数据。
磁浮轨道扫描点云数据包括移动扫描轨迹数据,将移动扫描轨迹数据插值生成间隔为1m的轨迹点集,利用相邻轨迹点,采用矢量投影的方法进行点云切片,具体原理如图2所示,A、B为相邻轨迹点,组成轨迹方向向量磁浮轨道扫描点云数据中任一点为C点,与A点组成向量则C点距离经过A点且垂直轨迹方向向量的平面的距离为L,其关系式为:
将切片点云数据经投影、旋转、平移生成二维点云(即轨道断面点云数据),根据点到点法式平面投影方法,将切片中的每一个点投影至经过轨迹点且法向量为轨迹方向向量的平面中,依据平面的法向量计算出将平面旋转至xoy平面的旋转矩阵R,以及旋转后在xoy平面中点坐标的平均值则生成的轨道断面点云坐标(xi,yi,0)为:
依据磁浮轨道的结构和相关设计参数,如图4所示,绘制其断面CAD图,如图5所示,其中,断面CAD图的中心点位于图中原点位置,在CAD中经绘图处理将磁浮轨道断面线段“定距等分”化为点集,在CAD工具中利用“数据提取”将点集坐标提出制成磁浮轨道模型点集文件。
2)提取点云轨道断面初始匹配中心点的过程包括:针对轨道断面点云数据,对比计算轨道断面点云数据中的最高点集,即轨道平面点集,以一定的高度阈值粗略提取轨道平面点集,基于磁浮轨道平面的对称性,以轨道平面点集的坐标平均值代替此断面的轨道平面中心,得到点云轨道断面初始匹配中心点。
3)利用步骤2)中计算出的点云轨道断面初始匹配中心点,对磁浮轨道断面模型点集进行平移,再利用距离最小化的原则,在磁浮轨道断面模型点集中计算离轨道断面点云数据中每一个点云数据点最近的模型点重构模型点集,寻求结束后,再利用平移参数反向平移模型点集,得到适应性模型点集。
4)为了精确计算适应性模型点集与轨道断面点云数据的匹配参数,采用基于最小二乘改进的二维ICP算法,具体步骤如下:4.1设定参考数据集(轨道断面点云数据)和目标数据集(适应性模型点集);4.2对目标数据集中的每个点,在参考数据集中寻找一个与之对应的最短距离的点;4.3建立匹配目标函数,对目标函数进行优化,求出目标函数最优解,得到新的目标数据集。
建立目标函数的具体过程如下:
其中,i为第i个数据点,x′i,y′i为参考数据点,b1和b2为平移参数,α为旋转矩阵参数,n为数据点个数,求解匹配参数α,b1,b2,使得f最小。
其中三个匹配参数解为:
则:
其中,R-1为xoy平面旋转至经过轨迹点且法向量为轨迹方向向量的平面的旋转矩阵。
5)基于提取的最终点云轨道中心点,利用拟合算法计算出磁浮轨道的中心线。
Claims (9)
1.一种磁浮轨道中心线提取方法,其特征在于,该方法包括以下步骤:
S1:获取轨道断面点云数据和磁浮轨道断面模型点集,所述磁浮轨道断面模型点集包括轨道模型平面中心点;
S2:从轨道断面点云数据中提取点云轨道断面初始匹配中心点;
S3:依据轨道模型平面中心点至点云轨道断面初始匹配中心点所需的平移参数平移磁浮轨道断面模型点集,基于距离最小化原则,利用轨道断面点云数据重构磁浮轨道断面模型点集,得到适应性模型点集;
S4:利用最小二乘改进的二维ICP算法计算适应性模型点集与轨道断面点云数据的匹配参数,得到最终点云轨道中心点;
S5:基于最终点云轨道中心点,利用拟合算法得到磁浮轨道的中心线;
S3中得到适应性模型点集的过程包括:
在平移后的磁浮轨道断面模型点集中提取距离轨道断面点云数据中每一个点云数据点最近的模型点,重构磁浮轨道断面模型点集,再利用平移参数反向平移重构的磁浮轨道断面模型点集,得到适应性模型点集。
2.根据权利要求1所述的一种磁浮轨道中心线提取方法,其特征在于,S1中,获取轨道断面点云数据和磁浮轨道断面模型点集包括以下步骤:
获取磁浮轨道扫描点云数据,基于磁浮轨道设计参数得到磁浮轨道断面模型点集;
基于移动扫描轨迹数据,对磁浮轨道扫描点云数据利用矢量投影进行切片,生成切片点云数据;
对切片点云数据进行投影、旋转和平移生成轨道断面点云数据。
5.根据权利要求2所述的一种磁浮轨道中心线提取方法,其特征在于,基于磁浮轨道设计参数得到磁浮轨道断面模型点集的过程包括:
根据磁浮轨道设计参数绘制磁浮轨道断面CAD图;
在磁浮轨道断面CAD图中将磁浮轨道断面线段定距离等分,得到磁浮轨道断面模型点集。
6.根据权利要求1所述的一种磁浮轨道中心线提取方法,其特征在于,S2中,提取点云轨道断面初始匹配中心点的过程包括:
从轨道断面点云数据中提取大于高度阈值的数据,得到轨道平面点集;
计算轨道平面点集的均值,得到点云轨道断面初始匹配中心点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011536006.2A CN112581521B (zh) | 2020-12-22 | 2020-12-22 | 一种磁浮轨道中心线提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011536006.2A CN112581521B (zh) | 2020-12-22 | 2020-12-22 | 一种磁浮轨道中心线提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112581521A CN112581521A (zh) | 2021-03-30 |
CN112581521B true CN112581521B (zh) | 2023-04-07 |
Family
ID=75138893
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011536006.2A Active CN112581521B (zh) | 2020-12-22 | 2020-12-22 | 一种磁浮轨道中心线提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112581521B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114396892B (zh) * | 2021-12-02 | 2023-08-25 | 重庆交通大学 | 轨道交通曲线轨道曲率测量方法 |
CN114187342A (zh) * | 2021-12-13 | 2022-03-15 | 中南大学 | 既有线路的钢轨点云提取方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156972A (zh) * | 2014-08-25 | 2014-11-19 | 西北工业大学 | 基于激光扫描测距仪与多相机融合的透视成像方法 |
CN109657698A (zh) * | 2018-11-20 | 2019-04-19 | 同济大学 | 一种基于点云的磁悬浮轨道障碍物检测方法 |
CN109670557A (zh) * | 2019-01-03 | 2019-04-23 | 中国矿业大学 | 一种基于杆状构筑物的公路点云自动配准方法 |
CN111489390A (zh) * | 2020-03-25 | 2020-08-04 | 山东科技大学 | 一种针对圆形隧道的中轴线提取方法 |
-
2020
- 2020-12-22 CN CN202011536006.2A patent/CN112581521B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156972A (zh) * | 2014-08-25 | 2014-11-19 | 西北工业大学 | 基于激光扫描测距仪与多相机融合的透视成像方法 |
CN109657698A (zh) * | 2018-11-20 | 2019-04-19 | 同济大学 | 一种基于点云的磁悬浮轨道障碍物检测方法 |
CN109670557A (zh) * | 2019-01-03 | 2019-04-23 | 中国矿业大学 | 一种基于杆状构筑物的公路点云自动配准方法 |
CN111489390A (zh) * | 2020-03-25 | 2020-08-04 | 山东科技大学 | 一种针对圆形隧道的中轴线提取方法 |
Non-Patent Citations (4)
Title |
---|
基于点云数据的隧道中轴线自动提取方法研究;李兵等;《重庆建筑》;20160125(第01期);全文 * |
基于盾构隧道环缝点云的中轴线提取方法;陈林等;《河海大学学报(自然科学版)》;20200125(第01期);全文 * |
对高速磁浮(轨道)平面缓和曲线线形的探讨;吴超平等;《交通科技》;20150615(第03期);全文 * |
徕卡MS50全站扫描仪在地铁隧道断面中的应用;李豪等;《测绘地理信息》;20180515(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112581521A (zh) | 2021-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113870123B (zh) | 基于车载移动激光点云的接触网导高与拉出值自动检测方法 | |
CN110986878B (zh) | 基于移动测量系统自动化提取铁轨断面的方法 | |
CN107301648B (zh) | 基于重叠区域边界角度的冗余点云去除方法 | |
CN109272524B (zh) | 一种基于阈值分割的小尺度点云噪声去噪方法 | |
CN112581521B (zh) | 一种磁浮轨道中心线提取方法 | |
CN106643551A (zh) | 一种叶片形状快速扫描装置及方法 | |
CN103106632B (zh) | 一种基于均值漂移的不同精度三维点云数据的融合方法 | |
CN107610223A (zh) | 基于LiDAR点云的电力塔三维重建方法 | |
CN111524103B (zh) | 一种基于三维激光点云的圆型隧道中轴线提取方法 | |
CN102798362A (zh) | 一种基于点云数据的铸件加工余量估计方法 | |
CN111179321B (zh) | 一种基于模板匹配的点云配准方法 | |
CN102750449B (zh) | 基于分步三维空间-特征域映射的点云直线特征提取方法 | |
CN111612902A (zh) | 一种基于雷达点云数据的煤矿巷道三维模型构建方法 | |
CN116204990B (zh) | 一种三维实测数据驱动的飞机大型骨架精准协调修配方法 | |
CN113689394A (zh) | 一种基于三维点云的隧道初支平整度检测方法 | |
CN110986877A (zh) | 基于高精度车载激光移动测量系统铁路工程限界检测方法 | |
CN112857315B (zh) | 基于三维激光扫描的钢管柱垂直度测量方法 | |
CN114049457A (zh) | 一种基于点云数据提取运营期隧道接触网高度和拉出值的方法 | |
CN115620021A (zh) | 一种基于三维激光点云的地下建筑物特征提取方法 | |
CN114897777A (zh) | 顾及空间关系的接触网支撑设施激光点云全自动提取方法 | |
CN116681912A (zh) | 铁路道岔的轨距检测方法及装置 | |
Li et al. | Research on Key Technologies for Container Ship Loading Test Based on 3D Laser Scanning | |
CN114742765A (zh) | 一种基于激光点云测量的隧道断面特征点精确提取方法 | |
CN118691836B (zh) | 基于种子引导的uav激光点云铁路轨道线提取方法及系统 | |
Zhao et al. | A Rapid Segmentation Method of Highway Surface Point Cloud Data Based on a Supervoxel and Improved Region Growing Algorithm |
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 |