基于点云数据智能识别运营期隧道侵界的方法与流程

1.本发明涉及一种基于点云数据智能识别运营期隧道侵界的方法。适用于隧道工程测量技术领域。


背景技术：

2.在城市轨道交通飞跃式发展的同时，相应的安全事故也在频繁发生，隧道坍塌事故给国家带来重大的经济损失，尤其是已建成隧道的安全运营，时刻影响着人民群众的人身安全，隧道安全检测课题的研究是十分必要的。
3.隧道安全检测其中一个很重要的指标就是隧道限界，隧道限界指为了保证地铁上的车在地铁轨道上安全运行，从车的角度要限制车体的横断面，使其不向外超界，从设备和建筑的角度要限制设备和建筑的横断面，使其不向内侵界。其中，关于车的横断面轮廓的参数指标即为车辆限界，关于设备和建筑轮廓的参数指标即为设备限界和建筑限界。
4.隧道限界近年来发展迅速，人们通过大量的研究和时间，研究出很多设备，用来解决隧道检测各个方面的问题。他们所采用的方法各不相同，但最终大都殊途同归，较好的完成了隧道检测任务。由于隧道呈狭长的管状结构，人们往往利用机车牵引相关设备进行检测，然而极个别严重的端点处仍然需要特殊处理。
5.国外早期的限界研究采用带有可伸缩框架的检测车进行人工检测，在不同的路段对检测车框架进行调整，然后再进行相关的检测。我国地铁限界检测目前水平较低，最常用的是断面检测法，根据限界标准制作可伸缩框架的限界车，全线逐段检测，缺点是在不同的路段就要重新调整检测车各伸缩杆的长度，耗费时间；当出现限界问题时，需要人工进行多次测量和记录，效率低下；受人工情绪的影响较大，从而使检测结果的误差较大。


技术实现要素：

6.本发明要解决的技术问题是：针对上述存在的问题，提供一种基于点云数据智能识别运营期隧道侵界的方法。
7.本发明所采用的技术方案是：一种基于点云数据智能识别运营期隧道侵界的方法，其特征在于：
8.s1、获取隧道的点云数据；
9.s2、基于隧道各里程处对应的断面设计坐标，得到隧道断面；
10.s3、根据隧道点云数据到隧道断面的距离从隧道点云数据中提取相应里程处的断面点云数据；
11.s4、基于断面点云数据提取两侧轨道的轨道顶中心点；
12.s5、以两侧轨道的轨道顶中心点之间连线的中点为坐标原点，建立限界坐标系；
13.s6、获取车辆限界模型并导入所述限界坐标系；
14.s7、判断限界坐标系中相应断面点云数据是否有点位于车辆限界模型范围内，如有，则表明该隧道断面存在侵界。
15.所述根据隧道点云数据到隧道断面的距离从隧道点云数据中提取相应里程处的断面点云数据，包括：
16.s31、基于隧道点云数据提取隧道中轴线；
17.s32、基于隧道中轴线在相应里程处与断面设计坐标所处坐标系的夹角转换隧道点云数据坐标；
18.s33、根据坐标转换后的隧道点云数据计算其与隧道断面之间的距离；
19.s34、提取距离小于一定值的点云数据作为断面点云数据。
20.所述基于隧道点云数据提取隧道中轴线，包括：
21.以采集所述隧道点云数据的三维激光扫描仪为坐标原点o，建立xyz三轴坐标系，其中y轴为隧道走向方向；
22.将隧道点云数据分别投影至xyz三轴坐标系中的xoy面和yoz面上，提取隧道投影量测的边界点集；
23.基于边界点集拟合得到xoy面上的两侧边界线和yoz面上的两侧边界线；
24.对xoy面上一侧的边界线采样得到xy边界线点集ⅰ(x
i1
，y
i1
)，计算该边界线上点集ⅰ法方向与另一侧边界线的交点，得到xy边界线点集ⅱ(x
i2
，y
i2
)；
25.计算xy边界线点集ⅰ、ⅱ中每一对边界点的中点坐标(x
i
，y
i
)；
26.对yoz面上一侧的边界线采样得到yz边界线点集ⅰ(y
i1
，z
i1
)，计算该边界线上点集ⅰ法方向与另一侧边界线的交点，得到yz边界线点集ⅱ(y
i2
，z
i2
)；
27.计算yz边界线点集ⅰ、ⅱ中每一对边界点的中点坐标(y
i
，z
i
)；
28.基于xoy面上的中点坐标(x
i
，y
i
)和yoz面上的中点坐标(y
i
，z
i
)，拟合得到隧道中轴线。
29.所述基于隧道中轴线在相应里程处与断面设计坐标所处坐标系的夹角转换隧道点云数据坐标，包括：
30.提取相应里程处隧道中轴线与xoy面和yoz面的夹角；
31.基于隧道中轴线与xoy面和yoz面的夹角转换所述隧道点云数据的坐标，使隧道点云数据的中轴线平行于y轴。
32.所述基于断面点云数据提取两侧轨道的轨道顶中心点，包括：
33.从断面点云数据中提取两侧轨道的轨道部分点云；
34.提取出轨道部分点云中的最高点；
35.提取出最高点一定范围内的轨道部分点云；
36.计算提取出的轨道部分点云的平均值作为相应轨道的轨道顶中心点。
37.所述判断限界坐标系中相应断面点云数据是否有点位于车辆限界模型范围内，采用inpolygon函数进行判定。
38.还包括：
39.s8、计算断面点云数据至车辆限界模型的最短距离，并形成图表。
40.一种基于点云数据智能识别运营期隧道侵界的装置，其特征在于：
41.隧道点云获取模块，用于获取隧道的点云数据；
42.隧道断面获取模块，用于基于隧道各里程处对应的断面设计坐标，得到隧道断面；
43.断面点云提取模块，用于根据隧道点云数据到隧道断面的距离从隧道点云数据中
提取相应里程处的断面点云数据；
44.轨道中心提取模块，用于基于断面点云数据提取两侧轨道的轨道顶中心点；
45.限界坐标系构建模块，用于以两侧轨道的轨道顶中心点之间连线的中点为坐标原点，建立限界坐标系；
46.车辆限界获取模块，用于获取车辆限界模型并导入所述限界坐标系；
47.侵界判断模块，用于判断限界坐标系中相应断面点云数据是否有点位于车辆限界模型范围内，如有，则表明该隧道断面存在侵界。
48.一种存储介质，其上存储有能被处理器执行的计算机程序，其特征在于：所述计算机程序被执行时实现所述基于点云数据智能识别运营期隧道侵界的方法的步骤。
49.一种基于点云数据智能识别运营期隧道侵界的系统，其特征在于，包括：
50.三维激光扫描仪，用于采集隧道点云数据；
51.数据处理设备，具有存储器和处理器，存储器上存储有能被处理器执行的计算机程序，所述计算机程序被执行时实现所述基于点云数据智能识别运营期隧道侵界的方法的步骤。
52.本发明的有益效果是：本发明利用隧道点云数据结合车辆限界模型进行计算，判断是否存在侵界的情形，简单、快速而又准确。
53.本发明针对于隧道既有线，利用双向投影法得出隧道中轴线，基于隧道中轴线与某一坐标轴的角度，并进行旋转，用点到面的距离应用于断面点云的提取，基于断面点云进行隧道是否侵界的分析。
54.三维激光扫描技术是一种高精度、高效率、智能化的测量方式，通过本发明的后续处理可以快速准确的得到隧道任意里程的横断面以及判断出该断面是否侵界，与其它方法相比，该方法较为简单、快速、方便、充分利用海量点云原始数据、对于计算机配置无要求、高精度、高效率。
附图说明
55.图1为实施例的流程图。
56.图2为实施例中隧道中轴线提取过程图(a、某环附近点云的选取；b、双向投影到xoy平面和yoz平面；c、提取边界点后的yoz图形；d、提取边界点后的xoy图形；e、提取出的中轴线)。
57.图3为实施例中隧道旋转过程图(a、旋转前的投影图；b、旋转后的投影图)。
58.图4为实施例中隧道横断面示意图(a、提取断面图；b、投影后的断面图)。
59.图5为实施例中提取轨道部分点云示意图(a、提取左侧轨道点云；b、提取右侧轨道点云)。
60.图6为实施例中提取轨道顶中心点的示意图。
61.图7为实施例中坐标转换示意图。
62.图8为实施例中侵界判断示意图(a、限界模型导入图；b、判定侵界点)。
63.图9为实施例中的限界对比图(a、垂足点不在线段上点到线段最短距离示意图；b、垂足点在线段上点到线段最短距离示意图；c、断面点云到限界模型的距离)。
具体实施方式
64.本实施例为一种基于点云数据智能识别运营期隧道侵界的方法，包括以下步骤：
65.s1、获取隧道的点云数据。为获取较为精确的点云数据，需合理的选取测站间距和扫描分辨率，在测站间距和扫描分辨率确定后，进行现场扫描；在扫描的同时，用全站仪对扫描仪和标靶的位置进行测量，确定测站和标靶的绝对坐标。
66.三维激光扫描仪每次获取的点云数据都是以仪器中心为依据建立独立内部坐标系来获得点云坐标信息，对于多站点云数据只有通过各坐标系之间的精密转换，才能得到整个隧道统一点云坐标信息。隧道点云拼接一般基于标靶拼接，在相邻两站重叠区域布设多个标靶，并保证有三个以上非同一直线上的标靶位于两站之间的重叠区域。通过公共的标靶点参数转化完成整体点云的坐标系统一。在整体点云坐标系统一后，需要给标靶和测站点赋予用高精度全站仪测量的绝对坐标，使整体点云数据处于绝对坐标系下。
67.测站间距的主要影响因素是隧道的内径尺寸和激光扫描时的最大入射角。三维激光的入射角越大采集到的点云误差越大，而最大入射角θ
max
与测站间距s和隧道直径d有固定的几何关系，其建立的几何关系式为：
[0068][0069]
式(1)中：θ
max
为扫描范围内最大的入射角，s为测站间距，d为隧道的最大宽度。通过试验验证：当入射角大于60
°
时，误差急剧增加，由式(1)得，当θ
max
＝60
°
时，s＝1.7d。
[0070]
扫描仪设置的分辨率越高扫描得到的点云越密，测量精度就越高，扫描所需的时间也就更长，因此需要根据不同档次扫描分辨率的时间与点云数据精度的要求来合理的确定分辨率。在测站间距和扫描分辨率确定后，进行现场扫描，获取隧道点云数据。
[0071]
s2、基于隧道各里程处对应的断面设计坐标，得到隧道断面。
[0072]
根据隧道的设计文件，得到不同里程(环数)断面的设计坐标(x,y,z)，利用平面的法向量和设计坐标可得到该端面的平面方程，如公式(2)所示。
[0073][0074]
s3、根据隧道点云数据到隧道断面的距离从隧道点云数据中提取相应里程处的断面点云数据。
[0075]
s31、基于隧道点云数据提取隧道中轴线。隧道的中轴线是表示隧道走势的空间曲线，提取的隧道断面点云与隧道中轴线处处垂直，所以中轴线的准确性影响着隧道点云数据的后续处理，本实施例通过双向投影法提取特征点云，基于ransac算法并结合最小二乘法提取曲线隧道中轴线。
[0076]
本实施例中提取中轴线的具体步骤为：
[0077]
设以采集隧道点云数据的三维激光扫描仪为坐标原点o，建立xyz三轴坐标系，其中y轴为隧道走向方向，x轴与y轴垂直且两轴均位于横向扫描面内，z轴位于竖向扫描面内，三个轴互相垂直构成右手坐标系；
[0078]
将隧道点云数据分别投影至xoy面与yoz面上，提取隧道投影两侧的边界点集；
[0079]
通过边界点坐标集利用二次曲线方程(如公式(3))拟合边界线方程，利用ransac随机采样一致性算法结合最小二乘法拟合两条二次曲线。
[0080][0081]
对xoy面上一侧的边界线采样得到xy边界线点集ⅰ(x
i1
，y
i1
)，计算该边界线上点集ⅰ法方向与另一侧边界线的交点，得到xy边界线点集ⅱ(x
i2
，y
i2
)；
[0082]
计算xy边界线点集ⅰ、ⅱ中每一对边界点的中点坐标(x
i
，y
i
)；
[0083]
对yoz面上一侧的边界线采样得到yz边界线点集ⅰ(y
i1
，z
i1
)，计算该边界线上点集ⅰ法方向与另一侧边界线的交点，得到yz边界线点集ⅱ(y
i2
，z
i2
)；
[0084]
计算yz边界线点集ⅰ、ⅱ中每一对边界点的中点坐标(y
i
，z
i
)；
[0085]
基于xoy面上的中点坐标(x
i
，y
i
)和yoz面上的中点坐标(y
i
，z
i
)，拟合得到隧道中轴线。
[0086]
为提高中线拟合精度，将上述计算方法分别应用于两侧边界线从而计算出2组中点坐标进行基于ransac的参数估计的曲线拟合得到最终隧道中轴线。
[0087]
s32、基于隧道中轴线在相应里程处与断面设计坐标所处坐标系的夹角转换隧道点云数据坐标；
[0088]
提取相应里程处隧道中轴线与xoy面和yoz面的夹角α、β；
[0089]
基于隧道中轴线与xoy面和yoz面的夹角将原始隧道点云数据经过旋转、平移，使隧道点云数据的母线平行于坐标y轴。[x
i y
i z
i
]
t
和变化后的坐标[x
′
i y
′
i z
′
i
]
t
有式(4)中的转换关系：
[0090][0091]
s33、利用经过旋转、平移之后的隧道表面的点云数据到步骤s2所确定的相应里程处端面平面的距离进行截取隧道横断面，合理设置点到平面距离δ，根据公式(5)进行隧道的断面厚度的设置：
[0092]
|y
‑
y|<δ
ꢀꢀꢀ
(5)
[0093]
s34、提取隧道点云数据中满足公式(5)的点云数据作为断面点云数据。
[0094]
s4、基于断面点云数据提取两侧轨道的轨道顶中心点。
[0095]
提取里程方向对应某一环的断面点云数据进行分析，在运营期地铁隧道断面点云中可以确定轨道点云大致范围，采用坐标分割方法初步粗提取出轨道部分点云；
[0096]
在粗提取的一侧轨道部分点云中将提取出轨顶中心位置，方法为提取轨道部分点云中z坐标轴方向的最大值，设定一个阈值d，将z轴方向距离最高点阈值d范围内的点云提取出，计算此部分点云的平均值作为一侧轨道顶的中心点；同理提取另一侧轨道顶的中心点。
[0097]
s5、以两侧轨道的轨道顶中心点之间连线的中点为坐标原点，建立限界坐标系。
[0098]
通过提取出的两轨道顶中心点的坐标可以得知，两轨道顶部连线与原始坐标系夹角为α，将坐标系经过平移至两轨道顶部连线的中心点，通过公式(6)便可将坐标系转换到限界坐标系下。
[0099][0100]
s6、获取车辆限界模型并与相应的相应断面点云数据一同导入限界坐标系；
[0101]
s7、判断限界坐标系中相应断面点云数据是否有点位于车辆限界模型范围内，如有，则表明该隧道断面存在侵界。
[0102]
在车辆限界模型与断面点云形成的断面轮廓处于限界坐标系之下后，进行是否有侵界点的判断，即判断点是否在车辆限界模型的多边形内。
[0103]
本实施例采用inpolygon函数进行判定，inpolygon函数的作用是判断一个点是否存在于一个多边形内，它的用法是:
[0104]
[in on]＝inpolygon(x,y,xv,yv)
ꢀꢀꢀ
(8)
[0105]
其输入值x和y分别表示测试点的横坐标值和纵坐标值，而xv和yv则分别表示多边形定点的横坐标值和纵坐标值。若某一点在多边形内，则in返回值为1，否则in为0，若某一点在多边形的边界上则on为1，若不在边界上on为0。
[0106]
s8、计算断面点云数据至车辆限界模型的最短距离，并形成图表。计算点到车辆限界模型多边形距离，首先将多边形视为多个线段，计算点到每个线段的最短值。
[0107]
本实施例中计算方法为计算点云中一个点到每条线段的最短距离，利用点到直线的距离公式：
[0108][0109]
计算a点到直线bc的最短距离，然后判定d点是否在线段上，通过判断d点的横坐标与b点和d点横坐标关系实现。如果d点不在线段上，如图9(a)所示，计算bc线段两个端点到点a的距离，取其较小值为a点到线段bc的最小距离。如果d点在线段bc上，如图9(b)所示，则线段ad为点a到线段bc的最小距离。这样，依次计算每个点云到限界模型每个线段最短距离，取其最小值，此最小值为点到多边形最短距离。
[0110]
将断面点云按y轴向下为起始点，逆时针方向以隧道断面为中心点极坐标展开，画出如图9(c)所示的图形(横坐标为角度，纵坐标为距离)，即可较为方便的看出哪些方向点云侵界以及哪些点云容易侵界。
[0111]
本实施例还提供一种基于点云数据智能识别运营期隧道侵界的装置，包括隧道点云获取模块、隧道断面获取模块、断面点云提取模块、轨道中心提取模块、限界坐标系构建模块、车辆限界获取模块和侵界判断模块。
[0112]
本例中隧道点云获取模块用于获取隧道的点云数据；隧道断面获取模块用于基于隧道各里程处对应的断面设计坐标，得到隧道断面；断面点云提取模块用于根据隧道点云数据到隧道断面的距离从隧道点云数据中提取相应里程处的断面点云数据；轨道中心提取模块用于基于断面点云数据提取两侧轨道的轨道顶中心点；限界坐标系构建模块用于以两侧轨道的轨道顶中心点之间连线的中点为坐标原点，建立限界坐标系；车辆限界获取模块用于获取车辆限界模型并导入所述限界坐标系；侵界判断模块用于判断限界坐标系中相应断面点云数据是否有点位于车辆限界模型范围内，如有，则表明该隧道断面存在侵界。
[0113]
本实施例还提供一种存储介质，其上存储有能被处理器执行的计算机程序，该计算机程序被执行时实现本实施例中基于点云数据智能识别运营期隧道侵界的方法的步骤。
[0114]
本实施例还提供一种数据处理设备，具有存储器和处理器，存储器上存储有能被处理器执行的计算机程序，所述计算机程序被执行时实现本实施例中基于点云数据智能识别运营期隧道侵界的方法的步骤。
[0115]
本实施例还提供一种基于点云数据智能识别运营期隧道侵界的系统，包括三维激光扫描仪、全站仪和本例中的数据处理设备。