CN107644119A - 一种基于三维扫描点云的半孔率自动计算方法 - Google Patents

一种基于三维扫描点云的半孔率自动计算方法 Download PDF

Info

Publication number
CN107644119A
CN107644119A CN201710694502.2A CN201710694502A CN107644119A CN 107644119 A CN107644119 A CN 107644119A CN 201710694502 A CN201710694502 A CN 201710694502A CN 107644119 A CN107644119 A CN 107644119A
Authority
CN
China
Prior art keywords
blasthole
point
hole
point cloud
axis
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
Application number
CN201710694502.2A
Other languages
English (en)
Other versions
CN107644119B (zh
Inventor
李海波
徐昊
周家文
李孝文
杨兴国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan University
Original Assignee
Sichuan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan University filed Critical Sichuan University
Priority to CN201710694502.2A priority Critical patent/CN107644119B/zh
Publication of CN107644119A publication Critical patent/CN107644119A/zh
Application granted granted Critical
Publication of CN107644119B publication Critical patent/CN107644119B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明属于爆破技术领域,具体为一种基于三维扫描点云的半孔率自动计算方法,扫描获取爆破后岩壁面的空间信息,得到点云数据,拼合成三维模型;将炮孔的信息导入到模型中,生成炮孔轴线并修正到实际的轴线位置,识别并提取炮孔孔痕所在的点云;计算孔痕上的点到炮孔轴线的距离D,判断D与阈值D0的大小关系,得到所有半孔上的点;生成垂直于炮孔轴线且等间距的平面,面与半孔相交,根据所交的点,拟合其组成的圆弧,计算圆弧长度C,判断C与阈值C0的大小关系,得到所有完整半孔段;统计完整半孔段的个数,乘以平面间距的设定值,得到全部半孔长度L1,计算L1与炮孔总长度L的比值,得到该次爆破的半孔率。本发明工作量小、误差小、效率高。

Description

一种基于三维扫描点云的半孔率自动计算方法
技术领域
本发明属于爆破技术领域,具体为一种基于三维扫描点云的半孔率自动计算方法。
背景技术
在台阶爆破中,预裂爆破和光面爆破对于关键部位的开挖,起到非常重要的作用,其中,半孔率是评价预裂爆破和光面爆破效果的一项重要指标。
半孔率也称炮孔痕迹率,为开挖壁面上的炮孔痕迹总长与炮孔长的百分比。
在计算半孔率方面,传统的方法是手工皮卷尺测量,手工计算,这在工程应用中存在误差大、作业环境危险以及效率低下等缺点,测量的数据可靠度较低,不足以全面评价预裂爆破和光面爆破的效果。利用人工皮卷尺测量每一段完整半孔长度,然后手工将所有半孔长度相加,得到开挖壁面上的炮孔痕迹总长,除以炮孔长度,得到开挖壁面的半孔率。
该方法:1、工作量大:需要人工攀爬测量岩壁上的半孔长度,而在实际工作中,炮孔数量很多,导致需要测量的半孔工作量很大。
2、误差大:因为是人工用皮卷尺测量,测量工具落后,且存在人为操作上的的不准确,导致测量计算误差大。
3、作业环境危险:在实际测量工作中,需要人工攀爬至近乎垂直的岩壁上,且岩壁高度很高,给人工测量工作带来不便的同时也威胁着测量人员的安全。
申请号为201410467656.4,申请公布号为CN 104217199A的中国发明专利公开了一种基于近景摄影测量技术的半孔率统计方法:采用近景摄影测量技术对爆破后岩壁面进行非接触测量,获取壁面空间三维信息,建立壁面三维数字高程模型。通过影像识别技术,划分炮孔的孔痕,并根据孔痕在壁面上的不平整度来划分每一个炮孔半孔的长度。统计所有炮孔的半孔信息,计算得出评价此次爆破效果的半孔率。
该方法:1、测量精度差:由于该方法并非直接获取对象坐标,而是通过集合计算利用二维数码影像反算对象空间坐标,因此测量精度差。
2、工作环境要求高:因为是通过摄影机拍摄岩壁面的相片,所以对现场的光线条件要求较高,光线较差时无法进行工作。
3、对图像的处理繁琐:摄影测量技术在后期需要对图像经行人工拼接和矫正,此工程耗时较多。
4、自动化程度低:近景摄影测量技术对炮孔痕迹的识别和炮孔痕迹长度的统计都是手动完成,在实际半孔率统计中,炮孔数很多,手动操作费时费力。
发明内容
提出一种基于三维扫描点云的半孔率自动统计计算方法,解决传统测量计算过程中对于半孔率统计计算方法中存在的工作量大、误差大,效率低的缺陷。
具体技术方案如下:
一种基于三维扫描点云的半孔率自动计算方法,包括以下过程:
步骤1:用三维激光扫描仪分多站扫描获取爆破完成后整个岩壁面的空间信息,得到爆破后岩壁面的点云数据;
步骤2:将扫描得到的三维点云导入到点云处理软件中,通过点云拼接将各站扫描点云拼合成一幅完整的岩壁表面三维点云模型,并通过坐标转换将岩壁面三维点云模型转换成大地坐标;
步骤3:将炮孔的设计定位坐标、轴线方位及长度信息导入到三维点云模型中,并自动生成炮孔轴线,作为后续炮孔自动识别的参考;
步骤4:利用点云算法软件,通过最小二乘法将炮孔设计轴线方位修正到实际爆破完成后的轴线位置,并且以修正后的炮孔轴线为参考自动识别圈定并提取岩壁三维点云模型上的炮孔孔痕所在的点云区域;
步骤5:计算炮孔孔痕上的点到炮孔轴线的空间距离D,判断炮孔孔痕上点到炮孔轴线的距离D与设定的阈值D0之间的大小关系:当距离D大于设定的阈值D0时,排除该点;当距离D不大于设定的阈值D0时,将该点计入为半孔上的点。按上述方法依次判断所有炮孔孔痕上点是否为半孔上的点。
步骤6:生成一系列垂直于炮孔轴线且间距为设定值的平面,面与半孔相交,根据所交的点,拟合这些点所组成的圆弧,计算圆弧长度C,当相邻两个面所计算的C值均不小于设定的阈值C0时,认为这两个面之间的半孔段是完整的;否则,认为该半孔段是不完整的;
步骤7:统计所有完整半孔段的个数,乘以平面间距的设定值,得到所有的半孔长度L1,计算L1与炮孔总长度L的比值,得到该次爆破岩壁面的半孔率。
本发明提供的一种基于三维扫描点云的半孔率自动计算方法,带来的有益效果有:
1、工作量小:不用传统人工逐段手工测量计算,也没有近景摄影测量技术中对于图像的矫正、拼接以及逐个半孔长度计算累加等过程,本发明中的插件在输入点云数据后自动完成识别计算的全过程。
2、误差小:三维激光扫描的精确度可以到达100m的测距仅1mm的误差,人工测量时,每一个半孔长度的测量误差都是厘米级的,在累加的过程中误差会越积越大。
3、效率高:软件自动完成识别计算的全过程,只需将点云数据导入点云处理软件中,设定一系列的参数值就可以得到半孔率的结果。
附图说明
图1是实施例炮孔轴线的自动生成和矫正示意图;
图2是实施例以矫正后炮孔轴线为参考圈定的孔痕所在点云区域示意图;
图3是实施例炮孔点云的自动识别示意图。
具体实施方式
结合实施例说明本发明的具体技术方案。
如图1到图3所示,首先运用在某水电站地下厂房高边墙爆破开挖过程中的半孔率自动识别计算上。在打爆破孔之前,爆破设计人员会给出爆破孔的定位坐标、方位、直径、长度等信息。本案例中直径d=76mm,长度l=8m。
(1)用三维激光扫描仪分多次扫描某水电站地下厂房高边墙爆破开挖后的岩壁面,获取岩壁面的点云数据。
(2)将获取的点云数据导入到点云处理程序中,通过点云拼接将各站扫描点云拼合成一幅完整的岩壁表面三维点云模型,并通过坐标转换将岩壁面三维点云模型转换成大地坐标。
(3)将炮孔的设计定位点Pi的坐标、轴线方位及长度等信息导入到三维点云模型中,生成炮孔轴线Li。通过最小二乘法将炮孔设计轴线方位修正到实际爆破完成后的轴线位置,根据修正后的炮孔轴线自动识别出炮孔孔痕所在的点云。如图1所示,该次爆破岩壁面上,输入的炮孔轴线的起点为P1~P6,形成炮孔轴线6条,分别为L1~L6。炮孔总长度L为48m。
(4)以炮孔轴线为中轴线,半径r=100mm,作圆柱面并与点云相交,圆柱面之内的点云即是炮孔孔痕所在点云区域。
(5)计算炮孔孔痕所在区域上点到炮孔轴线Li的空间距离D,设计阈值D0=40mm,当距离D大于设定的阈值D0时,排除该点;当距离D不大于设定的阈值D0时,则将该点计入为半孔上的点。
(6)根据炮孔轴线的方程,生成一系列垂直于炮孔轴线且间距为0.05m的平面,面与半孔相交,根据所交的点,拟合这些点所组成的圆弧,计算圆弧长C,设定圆弧长的阈值C0=55mm,当相邻两个面所计算的C值均不小于设定的阈值C0时,认为这两个面之间的半孔段是完整的;否则,认为该半孔段是不完整的。图3中黑色线条部分为识别出的所有完整半孔段。
(7)统计所有完整半孔段的数量为624个,计算所有完整半孔段长度Lw=624*0.05m=31.2m,半孔率为Lw/L=31.2/48=65%。

Claims (1)

1.一种基于三维扫描点云的半孔率自动计算方法,其特征在于,包括以下过程:
步骤1:用三维激光扫描仪分多站扫描获取爆破完成后整个岩壁面的空间信息,得到爆破后岩壁面的点云数据;
步骤2:将扫描得到的三维点云导入到点云处理软件中,通过点云拼接将各站扫描点云拼合成一幅完整的岩壁表面三维点云模型,并通过坐标转换将岩壁面三维点云模型转换成大地坐标;
步骤3:将炮孔的设计定位坐标、轴线方位及长度信息导入到三维点云模型中,并自动生成炮孔轴线,作为后续炮孔自动识别的参考;
步骤4:利用点云算法软件,通过最小二乘法将炮孔设计轴线方位修正到实际爆破完成后的轴线位置,并且以修正后的炮孔轴线为参考自动识别圈定并提取岩壁三维点云模型上的炮孔孔痕所在的点云区域;
步骤5:计算炮孔孔痕上的点到炮孔轴线的空间距离D,判断炮孔孔痕上点到炮孔轴线的距离D与设定的阈值D0之间的大小关系:当距离D大于设定的阈值D0时,排除该点;当距离D不大于设定的阈值D0时,将该点计入为半孔上的点。按上述方法依次判断所有炮孔孔痕上点是否为半孔上的点。
步骤6:生成一系列垂直于炮孔轴线且间距为设定值的平面,面与半孔相交,根据所交的点,拟合这些点所组成的圆弧,计算圆弧长度C,当相邻两个面所计算的C值均不小于设定的阈值C0时,认为这两个面之间的半孔段是完整的;否则,认为该半孔段是不完整的;
步骤7:统计所有完整半孔段的个数,乘以平面间距的设定值,得到所有的半孔长度L1,计算L1与炮孔总长度L的比值,得到该次爆破岩壁面的半孔率。
CN201710694502.2A 2017-08-15 2017-08-15 一种基于三维扫描点云的半孔率自动计算方法 Active CN107644119B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710694502.2A CN107644119B (zh) 2017-08-15 2017-08-15 一种基于三维扫描点云的半孔率自动计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710694502.2A CN107644119B (zh) 2017-08-15 2017-08-15 一种基于三维扫描点云的半孔率自动计算方法

Publications (2)

Publication Number Publication Date
CN107644119A true CN107644119A (zh) 2018-01-30
CN107644119B CN107644119B (zh) 2020-08-04

Family

ID=61110740

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710694502.2A Active CN107644119B (zh) 2017-08-15 2017-08-15 一种基于三维扫描点云的半孔率自动计算方法

Country Status (1)

Country Link
CN (1) CN107644119B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111125901A (zh) * 2019-12-19 2020-05-08 中国水利水电第七工程局有限公司 一种基于激光扫描的三维爆破可视化设计方法
CN111798476A (zh) * 2020-06-08 2020-10-20 国网江西省电力有限公司电力科学研究院 一种高压隔离开关导电臂轴线提取方法
CN114186645A (zh) * 2022-02-15 2022-03-15 中国长江三峡集团有限公司 一种爆破半孔率检测方法、装置、电子设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008121105A (ja) * 2006-11-10 2008-05-29 China Steel Corp 高炉炉体のセメント吹付け厚さと跳ね返り量の測定方法
CN102798412A (zh) * 2012-07-31 2012-11-28 同济大学 一种基于三维激光扫描评定隧道钻爆施工质量的方法
CN104217199A (zh) * 2014-09-15 2014-12-17 武汉大学 一种基于近景摄影测量技术的半孔率统计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008121105A (ja) * 2006-11-10 2008-05-29 China Steel Corp 高炉炉体のセメント吹付け厚さと跳ね返り量の測定方法
CN102798412A (zh) * 2012-07-31 2012-11-28 同济大学 一种基于三维激光扫描评定隧道钻爆施工质量的方法
CN104217199A (zh) * 2014-09-15 2014-12-17 武汉大学 一种基于近景摄影测量技术的半孔率统计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
施富强等: "三维数字化爆破质量评价技术", 《工程爆破》 *
李海波等: "基于三维激光扫描的隧洞开挖衬砌质量检测技术及其工程应用", 《岩石力学与工程学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111125901A (zh) * 2019-12-19 2020-05-08 中国水利水电第七工程局有限公司 一种基于激光扫描的三维爆破可视化设计方法
CN111125901B (zh) * 2019-12-19 2023-06-16 中国水利水电第七工程局有限公司 一种基于激光扫描的三维爆破可视化设计方法
CN111798476A (zh) * 2020-06-08 2020-10-20 国网江西省电力有限公司电力科学研究院 一种高压隔离开关导电臂轴线提取方法
CN111798476B (zh) * 2020-06-08 2023-10-20 国网江西省电力有限公司电力科学研究院 一种高压隔离开关导电臂轴线提取方法
CN114186645A (zh) * 2022-02-15 2022-03-15 中国长江三峡集团有限公司 一种爆破半孔率检测方法、装置、电子设备及存储介质
CN114186645B (zh) * 2022-02-15 2022-04-26 中国长江三峡集团有限公司 一种爆破半孔率检测方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
CN107644119B (zh) 2020-08-04

Similar Documents

Publication Publication Date Title
CN110285792B (zh) 一种无人机倾斜摄影精细格网土方计量方法
WO2021128675A1 (zh) 一种基于bim技术的地下连续墙施工方法
CN102798412B (zh) 一种基于三维激光扫描评定隧道钻爆施工质量的方法
CN104792274B (zh) 一种圆形隧道收敛变形的测量方法
CN112282847B (zh) 一种煤矿井下巷道形变监测方法
CN107066758A (zh) 基于无人机摄影技术及bim技术下的室外施工方法
CN112560154A (zh) 一种基于Unity平台的公路隧道工程BIM模型动态创建方法
CN109508508B (zh) 一种露天矿山治理勘查设计方法
CN107644119A (zh) 一种基于三维扫描点云的半孔率自动计算方法
CN110276732B (zh) 一种顾及地形特征线要素的山区点云空洞修复方法
CN109635340B (zh) 一种基于倾斜摄影和bim的矿山加工系统设计方法
CN107862115A (zh) 三维辐射场模型构建方法、装置、计算机设备及存储介质
CN107833281A (zh) 一种快速建立地质临空面gocad三维地质模型的方法
CN111199066A (zh) 一种基于bim+gis的施工场地虚拟构建复原方法
CN111578916A (zh) 基于bim技术在装饰工程中智能放线施工方法
CN115560690B (zh) 一种基于三维激光扫描技术的结构物整体变形分析方法
CN112288848A (zh) 无人机航拍三维建模计算工程量的方法
CN109341603A (zh) 基于激光点云技术监测爆破开挖面的平整度的方法
CN109358353B (zh) 一种集电线路基础快速分坑放样方法
CN106052553A (zh) 基于三维激光扫描的土建临时工程的设计优化方法及装置
CN104197867A (zh) 一种基于近景摄影测量技术的台阶面平整度统计方法
CN102495879A (zh) 基于地面lidar的工业架空管线量测方法
CN104217199B (zh) 一种基于近景摄影测量技术的半孔率统计方法
CN110990923A (zh) 水利水电工程双曲拱坝的快速工程计算方法、系统及介质
CN116448080B (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