CN116541639A - 一种基于点云数据的隧道超欠挖快速计算方法 - Google Patents
一种基于点云数据的隧道超欠挖快速计算方法 Download PDFInfo
- Publication number
- CN116541639A CN116541639A CN202310813010.6A CN202310813010A CN116541639A CN 116541639 A CN116541639 A CN 116541639A CN 202310813010 A CN202310813010 A CN 202310813010A CN 116541639 A CN116541639 A CN 116541639A
- Authority
- CN
- China
- Prior art keywords
- point cloud
- tunnel
- cloud data
- axis
- super
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 55
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000010606 normalization Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 9
- 238000012217 deletion Methods 0.000 claims description 3
- 230000037430 deletion Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 abstract description 13
- 238000005516 engineering process Methods 0.000 abstract description 4
- 239000011435 rock Substances 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000005641 tunneling Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C7/00—Tracing profiles
- G01C7/06—Tracing profiles of cavities, e.g. tunnels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Operations Research (AREA)
- Multimedia (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种基于点云数据的隧道超欠挖快速计算方法,包括:获取隧道点云数据与相应的轴线方程;依据轴线方程将点云归一化至原点;在测区范围内以隧道轴线的长度作为相对桩号,并依据桩号边界对点云分段;对分段点云单轴归一化;基于归一化点云,计算方位角与超欠挖值。本申请相较于传统隧道三维点云基于判断点距离断面距离的断面切片超欠挖计算方法,减少计算时间约85.3%(经采用一般笔记本电脑进行测试计算,在100万点、约400个断面的情况下计算超欠挖值仅耗时约6.41秒),解决了传统技术计算速率低下的问题,大幅提高了实时监测的效率,使得隧道安全监测可以更加实时化,便于第一时间反馈隧道围岩三维形变。
Description
技术领域
本发明涉及工程检测与监测技术领域,尤其涉及一种基于点云数据的隧道超欠挖快速计算方法。
背景技术
隧道的监控量测和断面检测是隧道施工安全控制的重要环节,随着三维激光扫描技术的发展,越来越多人将三维激光扫描技术引入隧道检测与监测工作当中,如隧道断面检测、超欠挖检测、监控量测等。现有应用三维激光扫描的隧道检测和监测方法中,大部分的隧道三维轮廓计算方法均需要通过点云切片方法来获取隧道点云断面。该方法先通过计算所有点云与每个待测隧道断面的距离,然后通过筛选距离值较小的点来作为待测断面的点,达到了点云切片的目的,从而进一步进行围岩轮廓计算。上述方法操作简便,且方法较为直观,所以大部分学者均以此作为获取点云断面的手段,如:[1]CN110706153B一种基于原始点云数据的隧道断面快速提取方法、[2]景冬,张光煜,刘如飞,王旻烨.一种隧道点云三维中轴线与横断面提取方法[J].激光与光电子学进展,2022,59(08):532-539.、[3]程效军,贾东峰,刘燕萍,程小龙.基于中轴线的隧道点云去噪算法[J].同济大学学报(自然科学版),2015,43(08):1239-1245.[4]林景峰,俞家勇,田茂义,徐飞,周茂伦.利用隧道激光点云提取中轴线及进行整体变形分析[J].遥感信息,2021,36(01):94-101。
该方法获取每个断面都需要计算大量点距离断面的距离,当需要进行隧道全断面三维变形、三维超欠挖计算时,又需要一次获取大量的断面,此时该方法总共将花费大量计算时间。点云数据共200万点,为46MB点云文件,断面数为450个,计算所用笔记本电脑(16年产品)CPU型号 I7-6700HQ,内存16G DDR3L,分别计算5次,平均耗时87.22秒。这明显严重拖慢了实时检测与实时监测效率,对于行业发展来说是不可忽视的问题,且对从业人员电脑配置要求也较高。
发明内容
本发明提供了一种基于点云数据的隧道超欠挖快速计算方法,包括以下步骤:
步骤一、获取隧道点云数据与相应的轴线方程;
步骤二、依据轴线方程将点云归一化至原点;
步骤三、在测区范围内以隧道轴线的长度作为相对桩号,并依据桩号边界对点云分段;
步骤四、对分段点云单轴归一化;
步骤五、基于归一化点云,计算方位角与超欠挖值。
可选的,本发明提供了一种基于点云数据的隧道超欠挖快速计算方法,其具体包括以下步骤:
步骤一、通过扫描装置对隧道进行扫描,获取隧道点云数据,并对隧道点云数据/>进行整平;坐标系/>轴正向为隧道桩号前进方向;
步骤二、对已知的轴线方程进行点云处理,得到转换后的点云数据/>和转换后的轴线方程/>;
步骤三、在测区范围内以隧道轴线的长度作为相对桩号,并设起点桩号为、终点桩号为/>、各个测量桩号/>;
通过积分法获得各个测量桩号所对应的/>坐标/>;
将转换后的点云数据在/>范围内分段,以获得份点云段数据/>、每段点云数据段/>对应的桩号边界所在桩号分别为/>和;
分别计算出在处/>的斜率/>以及在/>处/>的斜率/>,则有:
;
即得到:
;
其中:为测量编号且/>;/>为/>对应的桩号边界与隧道轴线交点的最小/>坐标,/>为/>对应的桩号边界与隧道轴线交点的最大/>坐标,且/>和/>均满足/>,/>为测量桩号的步长;/>为隧道宽度;/>为桩号边界最小/>值距离/>的距离,/>为桩号边界最大/>值距离/>的距离;/>是/>处/>的斜率/>与/>平面之间的夹角,/>是/>处/>的斜率/>与/>平面之间的夹角;
步骤四、对每个点云数据段进行单轴归一化处理,得到规一化点云切片/>;
步骤五、基于规一化点云切片,计算方位角与超欠挖值。
可选的,所述步骤一中的扫描装置设置为三维激光扫描仪或固态激光雷达。
可选的,所述步骤二中对轴线方程进行点云处理的具体方法如下:
针对基于地面三维激光扫描所获取的隧道点云数据,对隧道点云数据/>进行零点值归一化,以点云零点值为方位基准归一化至坐标原点;
针对基于固态激光雷达扫描所获取的隧道点云数据,对隧道点云数据/>进行最小值归一化,以最小值为方位基准归一化至坐标原点;
归一化后,得到转换后的点云数据和转换后的轴线方程/>,且转换后的点云数据/>和转换后的轴线方程/>均以隧道轴线过原点且原点处的切线斜率为0。
可选的,所述步骤二中针对隧道点云数据所获得的轴线方程/>进行点云处理的具体步骤如下:
①、计算轴线方程中/>处的切线与/>坐标轴的夹角/>;
②、使用四元数坐标转换方法将隧道点云数据和轴线方程/>同时在/>方向平移/>距离、在/>方向平移/>距离;
③、以坐标轴原点为旋转点旋转角度,得到转换后的点云数据/>和转换后的轴线方程/>。
可选的,所述步骤四中得到规一化点云切片的具体方法为:
计算轴线方程中/>处的切线与/>坐标轴的夹角/>,并将每个点云数据段使用四元数坐标转换方法分别在/>处以坐标轴原点为旋转原点旋转/>角度、以及在方向平移/>距离、在/>方向平移/>距离,以将/>和/>分别转换为/>和/>,作为近似桩号边界;
直接删除近似桩号边界和/>外的点,得到规一化点云切片/>。
可选的,计算方位角与超欠挖值的具体方法为:
以规一化点云切片的/>轴作为竖轴、/>轴作为横轴、/>和/>作为原点,建立断面坐标系/>;基于断面坐标系/>,将所有规一化点云切片/>当做同一断面,计算所有/>的方位角和超欠挖值,一次性得到形状类似柱状的隧道超欠挖测量值集合;
其中:为各个断面测量桩号/>的集合向量,/>分别为/>各个断面点云所对应的/>坐标向量和/>坐标向量,/>为各个断面点云所对应的超欠挖值向量,/>为各个断面点云所对应的断面方位角,/>共包含/>个断面点云/>。
与现有技术相比,本发明具有以下有益效果:
本申请包括依据已知的轴线方程将点云归一化至原点后;在测区范围内以隧道轴线的长度作为相对桩号,并依据桩号边界对点云分段;对分段点云单轴归一化;基于归一化点云,计算方位角与超欠挖值;通过采用上述步骤,使其相较于传统隧道三维点云基于判断点距离断面距离的断面切片超欠挖计算方法,减少计算时间约85.3%(经采用一般笔记本电脑进行测试计算,在100万点、约400个断面的情况下计算超欠挖值仅耗时约6.41秒),解决了传统技术计算速率低下的问题,大幅提高了实时监测的效率,使得隧道安全监测可以更加实时化,便于第一时间反馈隧道围岩三维形变;且采用本申请中所记载的隧道超欠挖快速计算方法不仅适用于TLS的广水平视角的非重复扫描点云数据,也适用于激光雷达的窄水平视角的重复扫描点云数据。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本实施例中一种基于点云数据的隧道超欠挖快速计算方法的整体流程示意图;
图2是桩号边界示意图;
图3是规则化点云切片示意图;
图4是柱状隧道超欠挖测量值集合示意图;
图5是柱状隧道超欠挖测量值集合图示。
具体实施方式
以下是本发明的具体实施例并结合附图,对本发明的技术方案作进一步的描述,但本发明并不限于这些实施例。
除非另有定义,本文所使用的所有的技术和科学术语与本发明的技术领域技术人员通常理解的含义相同。本文中在本发明说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。
本实施例:
参见图1所示,本发明所提供的一种基于点云数据的隧道超欠挖快速计算方法,具体步骤包括以下:
步骤一、通过地面三维激光扫描(Terrestrial Laser Scanning, TLS)或固态激光雷达(固态激光雷达包括纯固态激光雷达和半固态激光雷达)对隧道进行扫描,获取隧道点云数据,并对隧道点云数据/>进行整平;坐标系/>轴正向为隧道桩号前进方向。
步骤二、对已知的轴线方程进行点云处理,具体为:
针对基于地面三维激光扫描所获取的隧道点云数据,对隧道点云数据/>进行零点值归一化,以点云零点值为方位基准归一化至坐标原点,得到转换后的点云数据/>和转换后的轴线方程/>,且转换后的点云数据/>和转换后的轴线方程/>均以隧道轴线过原点且原点处的切线斜率为0;
针对基于固态激光雷达扫描所获取的隧道点云数据,对隧道点云数据/>进行最小值归一化,以最小值为方位基准归一化至坐标原点;归一化后,得到转换后的点云数据和转换后的轴线方程/>,且转换后的点云数据/>和转换后的轴线方程/>均以隧道轴线过原点且原点处的切线斜率为0。此处优选:隧道轴线方程/>的计算方式参考现有技术(具体如:程效军,贾东峰,刘燕萍,程小龙.基于中轴线的隧道点云去噪算法[J].同济大学学报(自然科学版),2015,43(08):1239-1245.)。
步骤三、在测区范围内以隧道轴线的长度作为相对桩号,并设起点桩号为、终点桩号为/>、各个测量桩号/>;具体的,固态激光雷达数据起点桩号/>;由于三维激光扫描仪起点桩号不确定,需根据现场实际需求确定,因此,三维激光扫描起点桩号一般设置为-40到-20之间;固态激光雷达在测区范围的相对桩号范围一般设置为0-50m,三维激光扫描在测区范围的相对桩号范围一般设置为-40m-40m;
通过积分法获得各个测量桩号所对应的/>坐标/>;
将转换后的点云数据在/>范围内分段,以获得份点云段数据/>、每段点云数据段/>对应的桩号边界分别为/>和/>;
分别计算出在处/>的斜率/>以及在/>处/>的斜率/>,则有:
;
即得到:
;
其中:为测量编号且/>;/>为/>对应的桩号边界与隧道轴线交点的最小/>坐标,/>为/>对应的桩号边界与隧道轴线交点的最大/>坐标,且/>和/>均满足/>,/>为测量桩号的步长(一般取值0.1m);/>为隧道宽度;为桩号边界最小/>值距离/>的距离,/>为桩号边界最大/>值距离/>的距离;/>是/>处/>的斜率/>与/>平面之间的夹角,/>是/>处/>的斜率与/>平面之间的夹角。
步骤四、计算中/>处的切线与/>坐标轴的夹角/>,并将每个点云数据段/>使用四元数坐标转换方法分别在/>处以坐标轴原点为旋转原点旋转/>角度、以及在/>方向平移/>距离、在/>方向平移/>距离,以将/>和/>分别转换为和/>,作为近似桩号边界;
直接删除近似桩号边界和/>外的点,得到规一化点云切片/>(具体的,规则化点云切片示意图见图3所示)。因《盾构法隧道施工及验收规范》中指出,地铁隧道平面曲线半径一般不小于 300 m,即表示一般情况隧道平面曲线半径均大于300m,假设隧道平面曲线半径为300m,则由几何关系可以推出/>和/>处的/>坐标与/>坐标最大值与处的/>坐标与/>坐标最大值差值均小于/>m,则可以认为其他平面曲线半径大于300m的隧道计算得到差值比该值更小,该误差远小于1mm(一般施工监测精度要求),对于工程测量可以忽略不计,可以认为/>范围内的点均属于/>断面。
步骤五、以规一化点云切片的/>轴作为竖轴、/>轴作为横轴、/>和/>作为原点,以建立断面坐标系/>;基于断面坐标系/>,将所有规一化点云切片/>视为同一断面,同时计算/>的方位角和超欠挖值,且计算前后均无需再进行坐标转换步骤,得到形状类似柱状的隧道超欠挖测量值集合/>;其中:/>为各个断面测量桩号的集合向量即/>,/>分别为各个点的/>、/>坐标向量,/>为各个点的超欠挖值向量,/>为各个点的断面方位角,/>共包含/>个断面点云/>即。
可选的,所述步骤二中针对基于地面三维激光扫描所获取的隧道点云数据所获得的轴线方程/>进行点云处理的具体步骤如下:
①、计算轴线方程中/>处的切线与/>坐标轴的夹角/>;
②、使用四元数坐标转换方法将隧道点云数据和轴线方程/>同时在/>方向平移/>距离、在/>方向平移/>距离;
③、以坐标轴原点为旋转点旋转角度,得到转换后的点云数据/>和转换后的轴线方程/>。
可选的,所述步骤二中针对基于固态激光雷达扫描所获取的隧道点云数据获得的轴线方程/>进行点云处理的具体步骤如下:
Ⅰ、计算轴线方程在点云/>坐标最小值/>处中的切线与/>坐标轴的夹角;
Ⅱ、使用四元数坐标转换方法将隧道点云数据和轴线方程/>同时在/>方向平移/>距离、在/>方向平移/>距离以及在/>方向平移距离后,以坐标轴原点为旋转原点旋转/>角度,得到转换后的点云数据和转换后的轴线方程/>。
实施例一:
以固态激光雷达对隧道进行扫描为例,具体说明如下:
步骤1)、通过固态激光雷达获取某段已通过倾角仪将扫描点云整平后的隧道点云数据,基于该隧道点云数据/>获取相应隧道点云处的轴线方程/>,以隧道轴线为坐标系/>轴,坐标系/>轴正向为隧道桩号前进方向。
步骤2)、设定点云测量范围为距离仪器15m至55m范围,对隧道点云数据进行最小值归一化,将所有隧道点云数据/>以最小值为方位基准归一化至坐标原点,得到转换后的点云数据/>和转换后的轴线方程/>,且转换后的点云数据/>和转换后的轴线方程均以隧道轴线过原点且原点处的切线斜率为0。具体步骤为:计算轴线方程/>在点云/>处中的切线与/>坐标轴的夹角/>,使用四元数坐标转换方法将隧道点云数据和轴线方程/>同时在/>方向平移/>距离、在/>方向平移/>距离以及在/>方向平移/>距离后,以坐标轴原点为旋转原点旋转/>角度,得到转换后的点云数据/>和转换后的轴线方程/>。
步骤3)、在测区范围内以隧道轴线的长度作为相对桩号(即起点桩号),设终点桩号/>(相对桩号最大值与点云测量范围最大值相似)、测量桩号步长,则测量桩号/>,/>为测量编号且设。
通过积分获得各个测量桩号对应的/>坐标/>,之后分别在范围内分段,获得401份点云段/>,每段点云数据段/>对应的桩号边界即为桩号/>(即/>)和/>,其中:/>和分别为/>对应的桩号边界与隧道轴线交点的/>坐标,且其满足;/>和/>分别为桩号边界最小/>值距离/>的距离和桩号边界最大/>值距离/>的距离;
再计算出和/>处/>的斜率/>和/>,则有:
;
即得到:
;
其中:为隧道宽度,该隧道/>。桩号边界示意图如图2所示。
步骤4)、计算中/>处的切线与/>坐标轴的夹角/>,每个/>分别使用四元数坐标转换方法在/>处以坐标轴原点为旋转原点旋转/>角度,并在/>方向和/>方向分别平移/>和/>距离,/>和/>转换后为/>和/>,作为近似桩号边界;此时直接删除近似桩号边界/>和/>外的点,得到规一化点云切片/>,规则化点云切片示意图见图3所示。由几何关系可以推出/>和/>处的/>坐标与/>坐标最大值与/>处的坐标与/>坐标最大值差值均小于/>m,该误差远小于1mm(一般施工监测精度要求),对于工程测量可以忽略不计,可以认为/>范围内的点均属于/>断面。
步骤5)、以规一化点云切片的/>轴作为竖轴,/>轴作为横轴,/>和/>作为原点,以建立断面坐标系/>,将所有规一化点云切片/>视为同一断面,使用现有方法(CN201810768199.0一种基于激光点云的隧道超欠挖数值计算方法中断面超欠挖值计算方法)同时计算/>的方位角和超欠挖值,且计算前后均无需再进行坐标转换步骤,得到形状类似柱状的隧道超欠挖测量值集合/>,其中:/>分别为各个点的、/>坐标向量,/>为各个点的超欠挖值向量,/>为各个点的断面方位角,/>共包含/>个断面点云/>,每个断面对应桩号/>。柱状隧道超欠挖测量值集合示意图和该次扫描点云计算得到的柱状隧道超欠挖测量值集合图示见图4和图5所示。
接下来分别通过对比现有基于计算点距离断面的距离进行切片的方法和本文方法的计算速度,点云数据共200万点,为46MB点云文件,断面数为401个,计算所用笔记本电脑(16年产品)CPU型号 I7-6700HQ,内存16G DDR3L,每个方法分别计算5次:具体参数参见下表所示,
由上表可知,本文基于坐标转换的点云切片方法相对于目前流行的基于计算点云与断面距离的点云切片方法,200万点数据计算用时减少约74.4秒,减少幅度约85.3%,大幅提高计算效率,为基于激光雷达的隧道实时监测提供可能。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于点云数据的隧道超欠挖快速计算方法,其特征在于,包括以下步骤:
步骤一、获取隧道点云数据与相应的轴线方程;
步骤二、依据轴线方程将点云归一化至原点;
步骤三、在测区范围内以隧道轴线的长度作为相对桩号,并依据桩号边界对点云分段;
步骤四、对分段点云单轴归一化;
步骤五、基于归一化点云,计算方位角与超欠挖值。
2.根据权利要求1所述的基于点云数据的隧道超欠挖快速计算方法,其特征在于,具体包括以下步骤:
步骤一、通过扫描装置对隧道进行扫描,获取隧道点云数据,并对隧道点云数据/>进行整平;坐标系/>轴正向为隧道桩号前进方向;
步骤二、对已知的轴线方程进行点云处理,得到转换后的点云数据/>和转换后的轴线方程/>;
步骤三、在测区范围内以隧道轴线的长度作为相对桩号,并设起点桩号为、终点桩号为/>、各个测量桩号/>;
通过积分法获得各个测量桩号所对应的/>坐标/>;
将转换后的点云数据在/>范围内分段,以获得/>份点云段数据/>、每段点云数据段/>对应的桩号边界所在桩号分别为/>和/>;
分别计算出在处/>的斜率/>以及在/>处/>的斜率/>,则有:
;
即得到:
;
其中:为测量编号且/>;/>为/>对应的桩号边界与隧道轴线交点的最小/>坐标,/>为/>对应的桩号边界与隧道轴线交点的最大/>坐标,且/>和/>均满足,/>为测量桩号的步长;/>为隧道宽度;/>为桩号边界最小值距离/>的距离,/>为桩号边界最大/>值距离/>的距离;/>是/>处/>的斜率与/>平面之间的夹角,/>是/>处/>的斜率/>与/>平面之间的夹角;
步骤四、对每个点云数据段进行单轴归一化处理,得到规一化点云切片/>;
步骤五、基于规一化点云切片,计算方位角与超欠挖值。
3.根据权利要求2所述的基于点云数据的隧道超欠挖快速计算方法,其特征在于,所述步骤一中的扫描装置设置为三维激光扫描仪或固态激光雷达。
4.根据权利要求2所述的基于点云数据的隧道超欠挖快速计算方法,其特征在于,所述步骤二中对轴线方程进行点云处理的具体方法如下:
针对基于地面三维激光扫描所获取的隧道点云数据,对隧道点云数据/>进行零点值归一化,以点云零点值为方位基准归一化至坐标原点;
针对基于固态激光雷达扫描所获取的隧道点云数据,对隧道点云数据/>进行最小值归一化,以最小值为方位基准归一化至坐标原点;
归一化后,得到转换后的点云数据和转换后的轴线方程/>,且转换后的点云数据和转换后的轴线方程/>均以隧道轴线过原点且原点处的切线斜率为0。
5.根据权利要求4所述的基于点云数据的隧道超欠挖快速计算方法,其特征在于,所述步骤二中针对隧道点云数据所获得的轴线方程/>进行点云处理的具体步骤如下:
①、计算轴线方程中/>处的切线与/>坐标轴的夹角/>;
②、使用四元数坐标转换方法将隧道点云数据和轴线方程/>同时在/>方向平移/>距离、在/>方向平移/>距离;
③、以坐标轴原点为旋转点旋转角度,得到转换后的点云数据/>和转换后的轴线方程/>。
6.根据权利要求2所述的基于点云数据的隧道超欠挖快速计算方法,其特征在于,所述步骤四中得到规一化点云切片的具体方法为:
计算轴线方程中/>处的切线与/>坐标轴的夹角/>,并将每个点云数据段/>使用四元数坐标转换方法分别在/>处以坐标轴原点为旋转原点旋转/>角度、以及在/>方向平移/>距离、在/>方向平移/>距离,以将/>和/>分别转换为/>和/>,作为近似桩号边界;
直接删除近似桩号边界和/>外的点,得到规一化点云切片/>。
7.根据权利要求2所述的基于点云数据的隧道超欠挖快速计算方法,其特征在于,计算方位角与超欠挖值的具体方法为:
以规一化点云切片的/>轴作为竖轴、/>轴作为横轴、/>和/>作为原点,建立断面坐标系/>;基于断面坐标系/>,将所有规一化点云切片/>当做同一断面,计算所有的方位角和超欠挖值,一次性得到形状类似柱状的隧道超欠挖测量值集合;
其中:为各个断面测量桩号/>的集合向量,/>分别为/>各个断面点云所对应的坐标向量和/>坐标向量,/>为各个断面点云所对应的超欠挖值向量,/>为各个断面点云所对应的断面方位角,/>共包含/>个断面点云/>。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310813010.6A CN116541639B (zh) | 2023-07-05 | 2023-07-05 | 一种基于点云数据的隧道超欠挖快速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310813010.6A CN116541639B (zh) | 2023-07-05 | 2023-07-05 | 一种基于点云数据的隧道超欠挖快速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116541639A true CN116541639A (zh) | 2023-08-04 |
CN116541639B CN116541639B (zh) | 2023-09-19 |
Family
ID=87454452
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310813010.6A Active CN116541639B (zh) | 2023-07-05 | 2023-07-05 | 一种基于点云数据的隧道超欠挖快速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116541639B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017157967A1 (en) * | 2016-03-14 | 2017-09-21 | Imra Europe Sas | Processing method of a 3d point cloud |
CN110030972A (zh) * | 2019-04-24 | 2019-07-19 | 中铁八局集团第三工程有限公司 | 基于ExcelVBA的隧道超欠挖检测方法 |
CN114993203A (zh) * | 2022-05-26 | 2022-09-02 | 西南交通大学 | 一种基于初支不等厚的隧道变形监测方法 |
CN115830587A (zh) * | 2022-12-02 | 2023-03-21 | 四川大学 | 一种基于高精度点云数据的结构面快速自动识别方法 |
CN116304764A (zh) * | 2023-05-22 | 2023-06-23 | 湖南联智科技股份有限公司 | 一种施工隧道点云中轴线自动拟合方法 |
-
2023
- 2023-07-05 CN CN202310813010.6A patent/CN116541639B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017157967A1 (en) * | 2016-03-14 | 2017-09-21 | Imra Europe Sas | Processing method of a 3d point cloud |
CN110030972A (zh) * | 2019-04-24 | 2019-07-19 | 中铁八局集团第三工程有限公司 | 基于ExcelVBA的隧道超欠挖检测方法 |
CN114993203A (zh) * | 2022-05-26 | 2022-09-02 | 西南交通大学 | 一种基于初支不等厚的隧道变形监测方法 |
CN115830587A (zh) * | 2022-12-02 | 2023-03-21 | 四川大学 | 一种基于高精度点云数据的结构面快速自动识别方法 |
CN116304764A (zh) * | 2023-05-22 | 2023-06-23 | 湖南联智科技股份有限公司 | 一种施工隧道点云中轴线自动拟合方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116541639B (zh) | 2023-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111927558B (zh) | 动水软弱围岩隧道全断面掘进的安全预警方法及装置 | |
Mah et al. | 3D laser imaging for surface roughness analysis | |
US20190339210A1 (en) | Multi-sensor pipe inspection system and method | |
WO2018137103A1 (zh) | 一种基于多源遥感数据的流域污染检测方法及*** | |
CN110243300B (zh) | 基于机器视觉技术的盾尾间隙测量方法及*** | |
JP2018531402A6 (ja) | 斜面安定性ライダー | |
CN103940356A (zh) | 一种基于三维激光扫描技术的建筑物整体变形监测方法 | |
CN106767524A (zh) | 一种水力机械叶片曲面检测方法和装置 | |
Watson et al. | Estimating underground mine ventilation friction factors from low density 3D data acquired by a moving LiDAR | |
CN102768022A (zh) | 采用数码照相技术的隧道围岩变形检测方法 | |
CN113091639A (zh) | 盾构隧道管片错台提取方法 | |
KR102001789B1 (ko) | 구조물 손상 검사장치 및 그 제어방법 | |
CN112233233B (zh) | 一种隧道相对点云数据的三维线形恢复方法 | |
Weixing et al. | Extraction of tunnel center line and cross-sections on fractional calculus, 3D invariant moments and best-fit ellipse | |
CN116304764A (zh) | 一种施工隧道点云中轴线自动拟合方法 | |
CN106813588A (zh) | 一种外浮顶储罐群变形监测方法 | |
CN110006432B (zh) | 一种基于几何先验信息下的室内机器人快速重定位的方法 | |
Tan et al. | Combining mobile terrestrial laser scanning geometric and radiometric data to eliminate accessories in circular metro tunnels | |
Dong et al. | The method for accurate acquisition of pavement macro-texture and corresponding finite element model based on three-dimensional point cloud data | |
Mehrishal et al. | A semi-automatic approach for joint orientation recognition using 3D trace network analysis | |
CN116541639A (zh) | 一种基于点云数据的隧道超欠挖快速计算方法 | |
CN115880597B (zh) | 一种基于遥感技术的岩溶区落水洞提取方法 | |
CN104316960B (zh) | 一种基于vsp的储层裂缝预测方法及*** | |
CN105716565A (zh) | 一种基于分布式角度信息的架空输电线路舞动的重建方法 | |
CN110111376A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |