CN112344877A - 无人机测量大型岩体结构面三维形貌参数的装置及方法 - Google Patents
无人机测量大型岩体结构面三维形貌参数的装置及方法 Download PDFInfo
- Publication number
- CN112344877A CN112344877A CN202011252160.7A CN202011252160A CN112344877A CN 112344877 A CN112344877 A CN 112344877A CN 202011252160 A CN202011252160 A CN 202011252160A CN 112344877 A CN112344877 A CN 112344877A
- Authority
- CN
- China
- Prior art keywords
- point
- camera
- measurement
- rock mass
- aerial vehicle
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 46
- 239000013598 vector Substances 0.000 claims abstract description 49
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000001788 irregular Effects 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 90
- 239000011159 matrix material Substances 0.000 claims description 28
- 230000005540 biological transmission Effects 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000010008 shearing Methods 0.000 claims description 12
- 230000003287 optical effect Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 238000012876 topography Methods 0.000 claims description 9
- 238000012423 maintenance Methods 0.000 claims description 7
- 239000000725 suspension Substances 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 230000009471 action Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 6
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 230000002441 reversible effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 7
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012625 in-situ measurement Methods 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
Images
Classifications
-
- 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/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
-
- 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/30—Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开一种无人机测量大型岩体结构面三维形貌参数的装置及方法,通过数字光栅投影设备向岩体结构面出射一系列连续的宽度不同的光栅条纹,接触岩体表面后受非规则凸体的影响,条纹发生变形,被左右两台工业摄像机捕捉并记录,并将空间位置坐标、拍摄俯仰角集成至图像信息中存储在图片中,以光栅投影作为捕捉特征点对左右摄像机的集成图像信息进行计算和拼接,利用最小二乘法解算被测表面上点的实际空间三维坐标,通过计算相关任意不共线点所构成的空间曲面的法向向量,进行简化三角剖分,再利用计算公式,获得岩石结构面表面的形貌参数。本发明适应性强、简洁便携、精度高,可有效测量工程背景下的岩体结构面三维形貌参数。
Description
技术领域
本发明涉及岩石结构面测量技术领域,尤其涉及一种无人机测量大型岩体结构面三维形貌参数的装置及方法。
背景技术
岩体结构面的粗糙度表征是岩石力学基本物理行为特征的重要衡量指标,直接影响岩体结构面的峰值剪切强度和初始剪胀角。工程中常使用JRC-JCS模型表征岩体结构面粗糙度,但利用结构面剖切面线投影的办法衡量粗糙度是以大幅度牺牲精度为前提的。G.Grasselli通过光学扫描技术,对结构面的几何形貌特征进行了细观的三维定义,提出了最大接触面积比A0、最大表面视倾角分布参数C作为表面形貌参数,并建立了结构面视倾角分布函,形成了对应峰值剪切强度的Grasselli模型,通过YANG J、TATONE B S A.等人的进一步完善形成具有较高精度和完整物理释义的Grasselli系模型。基于三维形貌参数及其对应峰值剪切强度模型是有效提升测量精度的手段,该方法已在实验室及小型野外试验中取得成熟应用,但其参数测量需要精密仪器,应用半径无法满足工程大规模测量的要求。因此寻找一种简易测量方法对于工程应用推广具有重大意义。
发明内容
针对上述现有技术的不足,本发明提供一种无人机测量大型岩体结构面三维形貌参数的装置及方法。
为解决上述技术问题,本发明所采取的技术方案是:一种无人机测量大型岩体结构面三维形貌参数的装置,包括:
测量模块、动态维稳模块、动力模块和计算机分析模块;
所述测量模块包括:两个工业摄像头和一个数字光栅投影设备;
所述动态维稳模块包括:下置单云台和水平连杆;
所述动力模块由无人机的动力系统构成;
所述计算机分析模块内置数据采集卡和基于反算JRC的数据校准模块;
所述动态维稳模块将下置单云台与无人机的底部中心位置连接,下置单云台的另一端与水平连杆的中心位置连接,并在水平连杆的中心位置悬挂数字光栅投影设备,在水平连杆的两端分别通过摄像头旋动机械转盘悬挂一个工业摄像头。
所述数字光栅投影设备的角度不可调节,两个工业摄像头的机身角度可按档位机械调节。
所述下置单云台的水平旋转角度为0°-180°,水平旋转速度为3°-10°/s,垂直旋转速度为4°/s。
所述工业摄像头具有无线传输能力,将捕捉到岩体结构表面反射的光栅条纹信息通过无线通信发送到计算机分析模块中。
另一方面,本发明还提供一种采用上述无人机测量大型岩体结构面三维形貌参数的装置进行大型岩体结构面三维形貌参数测量的方法,包括如下步骤:
步骤1:标定空间坐标,由无人机内置GPS系统定位空间基点;
步骤2:标定内部参数,将标定块型号和误差限值输入系统,将左右摄像头拍摄的标定块特征点图像传输至系统,系统通过摄影测量学原理自动确定测量装置内部参数,并建立点云空间坐标和实际空间坐标的转换矩阵M;
步骤3:确定测量范围,依据裸露岩体整体倾向的一致性划分测量区域;
步骤4:在保证一次测量为完整基面积的前提下,以测量移动次数最少为原则,依据测量范围和重叠等级确定基测量面积;
步骤5:对无人机辅助测量结构面三维形貌参数的装置进行测量前的校准和摄像头的角度调整,过程如下:
步骤5.1:确定摄像头角度,依据基测量面积确定摄像头角度,利用机械转盘将摄像头旋至基面积指定档位,锁死;
步骤5.2:共点校准,试飞无人机,开启测量和数据传输,系统检索点云数据,将左右摄像头点云数据进行共点匹配,若无匹配项则返回步骤2,重新标定内部参数;
步骤5.3:单点校准,试飞无人机,开启测量和数据传输,系统检索点云数据,关闭测量和数据传输,将云台旋转180°,无人机定点悬停旋转180°,开启测量和数据传输,系统检索点云数据,系统计算误差,镜头传感器自动校准;
步骤5.4:GPS定位校准,飞行器自系统自动校准;
步骤5.5:悬停校准,飞行器自系统自动校准。
步骤6:路径设定,利用辅助无人机测绘岩体空间图像,输入基面积和重叠等级,系统根据基测量面积对岩石表面进行分块,并生成路径方案和悬停点坐标;
所述对岩石表面进行分块时相邻的两块的重叠区域不小于30%。
步骤7:设置无人机按设定路径自动巡航,或人工指引飞行,并开启测量;
步骤8:数字光栅投影将不同条纹间隔密度的光栅重叠,向岩体结构面表面出射一系列连续的宽度不同的光栅条纹,实现灰度二进制编码;接触岩体表面后因受到表面非规则凸体的影响,条纹发生变形,被左右两台工业摄像机捕捉并记录下来,并将空间位置坐标、拍摄俯仰角集成至图像信息中存储在图片中,通过无线通信发送到计算机分析模块;
步骤9:将左右摄像机的集成图像信息进行共点匹配,利用左右摄像机像点坐标和转换矩阵计算被测表面上点的实际三维坐标,通过计算相关任意不共线点所构成的空间曲面的法向向量,利用最小二乘法解算实现三角剖分,获取唯一解数据,再利用表面形貌参数的计算公式,获得岩石结构面表面的形貌参数;
所述利用左右摄像机像点坐标和转换矩阵计算被测表面上点的实际三维坐标的方法如下:
(1)拾取结构面凸体点:
投影装置将光栅条纹投射在岩体结构面表面上,经凸体变形,分别由左右摄像头捕捉,两侧摄像头分别以独立内置坐标完成测量;单个摄像头的测量原理如下:
将岩体结构面表面某一凸体与相邻凸体边界围合平面的绝对平整延伸面作为参考平面Ω,在被测结构面上选取任一点P,称为凸体拾取点,P在参考平面Ω上的投影为P’,线段PP’的长度即点P距离Ω平面的高度为h,数字光栅投影设备结构中心与点P连线延长线交参考面Ω的交点为点A,摄像机光学中心与点P连线延长线交参考平面Ω的交点为点B;
数字光栅投影设备出射光线至参考平面Ω上的点A,摄像机捕捉到A,经过结构面凸体变形作用,该光线照射到凸体表面的点P,此时摄像机捕捉到B,即经过凸体变形作用点A移动到新的位置点B,AB的距离a耦合了凸体拾取点的高程信息h(X,Y),以此表示受到了表面形状的调制;
其中,a与系统的几何参数有关,即工业摄像机光心到参考平面Ω的距离d、数字光栅投影设备结构中心与工业摄像机光学中心的距离L有关,即:
h(X,Y)=f(d,L,a) (2)
(2)通过共点匹配,确定被测点在左右摄像机中的像点坐标:
对于任意一点P,其在两台工业摄像机内置坐标系统中单独形成坐标P1(X1,Y1,Z1)和P2(X2,Y2,Z2),空间内任意点的两坐标向量符合同一个转换矩阵函数,形式如下:
(3)利用左右摄像机中的像点坐标,计算被测点实际空间坐标:
经由方程
其中,u,v为像点在摄像机坐标体系下的齐次坐标,即
u1=X1/Z1,v1=Y1/Z1,u2=X2/Z2,v2=Y2/Z2,消去Z1、Z2,利用最小二乘法解得:
其中,(X0,Y0,Z0)为所求被测点实际空间坐标,(X1,Y1)(X2,Y2)为左右摄像头像点,M为转换矩阵,左、右摄像头各存在一个转换矩阵,分别表示为:
由系统结合公式(2)自动标定。
所述利用表面形貌参数的计算公式,获得岩石结构面表面的形貌参数的过程如下:
(1)利用三角剖分将扫描的结构面形貌做精简网格划分,利用定点坐标确定凸体表面外法向向量 为在参考平面Ω的投影向量,为剪切方向向量,系统默认沿参考平面Ω间隔60°生成各向剪切方向向量,按下式计算{θ*}:
tanθ*=-tanθcosα
其中,为剪切方向向量,为凸体表面外法向向量,为参考平面法向向量,为凸体表面外法向向量在参考平面的投影,θ*为视倾角,{θ*}为结构面各个凸体视倾角的集合,θ为凸体表面外法向向量和参考平面Ω法向向量夹角,α为剪切方向向量和凸体表面外法向向量在参考平面的投影向量的夹角;
(2)将{θ*}制成曲线,曲线与坐标轴的交点G、F分别为A0、θmax *,其中,为面向剪切方向大于某一视倾角的凸体面积之和与粗糙表面实际面积的比值,A0为最大可能接触面积比,定义为面向剪切方向的微元面积和粗糙表面实际面积的比值,θmax *为视倾角集合{θ*}中的最大值;
步骤10:待路径上所有的悬停点都测量完毕后,采用测量拼接的方法将多次测量结果转换到同一个坐标系下,删除重合区域,完成所有测量区域的拼合;
所述测量拼接的方法如下:
在区域内设置至少3个不共线的特征点,对两个区域分别进行三维坐标测量,得出特征点在2次测量中的坐标值;不同次测量中,特征点之间间距是相同的,利用特征点间距不变原理,进行坐标的空间变换,将2次测量结果转换到同一个坐标系下,删除重合区域,即完成了两个区域的拼接过程。
本实施例中,为方便特征点的提取,特征点被设计成外方内圆、外黑内白的样式,装置在测量时通过识别黑白相交的边界捕捉圆心坐标。设置特征点时应使特征点位于2个视角的重叠区域。
步骤11:系统反算JRC,系统将所测各基面积结构面三维形貌参数按式(1)反算出JRC值,现场随机选取岩体结构面可测量的区域测算JRC值,与反算结果比对,误差容许值以工程需要确定,其中反算公式的系统误差为±2%;
步骤12:复测,无法满足工程精度要求的,需重新测量。
采用上述技术方案所产生的有益效果在于:本发明提供的装置及方法可对实际工程的大幅天然岩体结构面进行原位测量,获得三维形貌参数一手数据,突破了实验室测量的限制,可大规模用于工程测量;测量精度略低于机械触针式表面形貌仪的精度,但这个对于岩石表面形貌的测量足够反映出岩石表面的细节;一次测量的最大有效面积可以通过调整摄像头的角度进行定量调节;质量轻、体积小、系统模块化,其中测量模块可与满足性能要求的动力模块任意组合;除必要手动校准外,可实现全周期系统自动校准。
附图说明
图1为本发明实施例中无人机测量大型岩体结构面三维形貌参数的装置的结构图;
图2为本发明实施例中结构面三维形貌参数测量方法的流程图;
图3为本发明实施例中结构面凸体点拾取原理示意图;
图4为本发明实施例中共点匹配原理示意图;
图5为本发明实施例中结构面表面空间信息编码及测量原理示意图;
图6为本发明实施例中结构面形貌特征三角网格化及其参数几何关系示意图;
其中,1为无人机动力系统,2为距离深感探测器,3为红外感知系统,4为下置单云台,5为连杆,6为摄像头悬动机械转盘,7为工业摄像头,8为数字光栅投影设备,9为工业摄像头,10为副摄像头,11为镜头智能校准系统。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
如图1所示,本实施例所述的一种无人机测量大型岩体结构面三维形貌参数的装置,包括:
测量模块、动态维稳模块、动力模块和计算机分析模块;
所述测量模块包括:两个工业摄像头7、9和一个数字光栅投影设备8;
本实施例中优先选用禅思H20T-四传感器摄像头;
所述动态维稳模块包括:下置单云台4和水平连杆5;
所述动力模块由无人机的动力系统1构成;
本实施例动力模块优先选择经纬M300 RTK行业级无人机,应具有高精度GPS定位、10公里以上的图传距离、2.5kg以上的负载能力。悬停精度应满足以下要求:垂直向,±0.1m(视觉定位正常工作时)、±0.5m(GPS正常工作时)、±0.1m(RTK定位正常工作时);水平向,±0.3m(视觉定位正常工作时)、±1.5m(GPS正常工作时)、±0.1m(RTK定位正常工作时)。
所述计算机分析模块内置数据采集卡和基于反算JRC的数据校准模块;
所述动态维稳模块将下置单云台4与无人机的底部中心位置连接,下置单云台4的另一端与水平连杆5的中心位置连接,并在水平连杆5的中心位置悬挂数字光栅投影设备8,在水平连杆5的两端分别通过摄像头旋动机械转盘6悬挂工业摄像头7、9。
所述数字光栅投影设备8的角度不可调节,两个工业摄像头7、9的机身角度可按档位机械调节。
所述下置单云台4的水平旋转角度为0°-180°,水平旋转速度为3°-10°/s,垂直旋转速度为4°/s。
所述工业摄像头7、9具有无线传输能力,将捕捉到岩体结构表面反射的光栅条纹信息通过无线通信发送到计算机分析模块中。
另一方面,本实施例还提供一种采用上述无人机测量大型岩体结构面三维形貌参数的装置进行结构面三维形貌参数测量的方法,其流程如图2所示,包括如下步骤:
步骤1:标定空间坐标,由无人机内置GPS系统定位空间基点;
步骤2:标定内部参数,将标定块型号和误差限值输入系统,将左右摄像头拍摄的标定块特征点图像传输至系统,系统通过摄影测量学原理自动确定测量装置内部参数,并建立点云空间坐标和实际空间坐标的转换矩阵M;
步骤3:确定测量范围,依据裸露岩体整体倾向的一致性划分测量区域;
本实施例中,依据裸露岩体整体倾向的一致性划分实验测量区域为1000mm×1000mm。
步骤4:在保证一次测量为完整基面积的前提下,以测量移动次数最少为原则,依据测量范围和重叠等级确定基测量面积;
优先选择测量面积为200mm×200mm、300mm×300mm、500mm×500mm中的一种,在本实施例中选择测量面积为300mm×300mm。
步骤5:对无人机测量大型岩体结构面三维形貌参数的装置进行测量前的校准和摄像头的角度调整,过程如下:
步骤5.1:确定摄像头角度,依据基测量面积确定摄像头角度,利用机械转盘将摄像头旋至基面积指定档位,锁死;
步骤5.2:共点校准,试飞无人机,开启测量和数据传输,系统检索点云数据,将左右摄像头点云数据进行共点匹配,若无匹配项则返回步骤2,重新标定内部参数;
步骤5.3:单点校准,试飞无人机,开启测量和数据传输,系统检索点云数据,关闭测量和数据传输,将云台旋转180°,无人机定点悬停旋转180°,开启测量和数据传输,系统检索点云数据,系统计算误差,镜头传感器自动校准;
步骤5.4:GPS定位校准,飞行器自系统自动校准;
步骤5.5:悬停校准,飞行器自系统自动校准。
步骤6:路径设定,利用辅助无人机测绘岩体空间图像,输入基面积和重叠等级,系统根据基测量面积对岩石表面进行分块,并生成路径方案和悬停点坐标;
所述对岩石表面进行分块时相邻的两块的重叠区域不小于30%。
步骤7:设置无人机按设定路径自动巡航,或人工指引飞行,并开启测量;
步骤8:数字光栅投影将不同条纹间隔密度的光栅重叠,向岩体结构面表面出射一系列连续的宽度不同的光栅条纹,实现灰度二进制编码,其原理如图5所示;接触岩体表面后因受到表面非规则凸体的影响,条纹发生变形,被左右两台工业摄像机捕捉并记录下来,并将空间位置坐标、拍摄俯仰角集成至图像信息中存储在图片中,通过无线通信发送到计算机分析模块;
步骤9:将左右摄像机的集成图像信息进行共点匹配,利用左右摄像机像点坐标和转换矩阵计算被测表面上点的实际三维坐标,通过计算相关任意不共线点所构成的空间曲面的法向向量,利用最小二乘法解算实现三角剖分,获取唯一解数据,再利用表面形貌参数的计算公式,获得岩石结构面表面的形貌参数;
所述利用左右摄像机像点坐标和转换矩阵计算被测表面上点的实际三维坐标的方法如下:
(1)拾取结构面凸体点,其示意图如图3所示:
投影装置将光栅条纹投射在岩体结构面表面上,经凸体变形,分别由左右摄像头捕捉,两侧摄像头分别以独立内置坐标完成测量;单个摄像头的测量原理如下:
将岩体结构面表面某一凸体与相邻凸体边界围合平面的绝对平整延伸面作为参考平面Ω,在被测结构面上选取任一点P,称为凸体拾取点,P在参考平面Ω上的投影为P’,线段PP’的长度即点P距离Ω平面的高度为h,数字光栅投影设备结构中心与点P连线延长线交参考面Ω的交点为点A,摄像机光学中心与点P连线延长线交参考平面Ω的交点为点B;
数字光栅投影设备出射光线至参考平面Ω上的点A,摄像机捕捉到A,经过结构面凸体变形作用,该光线照射到凸体表面的点P,此时摄像机捕捉到B,即经过凸体变形作用点A移动到新的位置点B,AB的距离a耦合了凸体拾取点的高程信息h(X,Y),以此表示受到了表面形状的调制;
其中,a与系统的几何参数有关,即工业摄像机光心到参考平面Ω的距离d、数字光栅投影设备结构中心与工业摄像机光学中心的距离L有关,即:
h(X,Y)=f(d,L,a) (2)
(2)通过共点匹配,确定被测点在左右摄像机中的像点坐标,其示意图如图4所示:
对于任意一点P,其在两台工业摄像机内置坐标系统中单独形成坐标P1(X1,Y1,Z1)和P2(X2,Y2,Z2),空间内任意点的两坐标向量符合同一个转换矩阵函数,形式如下:
(3)利用左右摄像机中的像点坐标,计算被测点实际空间坐标:
经由方程
其中,u,v为像点在摄像机坐标体系下的齐次坐标,即
u1=X1/Z1,v1=Y1/Z1,u2=X2/Z2,v2=Y2/Z2,消去Z1、Z2,利用最小二乘法解得:
其中,(X0,Y0,Z0)为所求被测点实际空间坐标,(X1,Y1)(X2,Y2)为左右摄像头像点,M为转换矩阵,左、右摄像头各存在一个转换矩阵,分别表示为:
由系统结合公式(2)自动标定。
所述利用表面形貌参数的计算公式,获得岩石结构面表面的形貌参数的过程如下:
(1)利用三角剖分将扫描的结构面形貌做精简网格划分,利用定点坐标确定凸体表面外法向向量为在参考平面Ω的投影向量,为剪切方向向量,其示意图如图6所示,系统默认沿参考平面Ω间隔60°生成各向剪切方向向量,按下式计算{θ*}:
tanθ*=-tanθcosα
其中,为剪切方向向量,为凸体表面外法向向量,为参考平面法向向量,为凸体表面外法向向量在参考平面的投影,θ*为视倾角,{θ*}为结构面各个凸体视倾角的集合,θ为凸体表面外法向向量和参考平面Ω法向向量夹角,α为剪切方向向量和凸体表面外法向向量在参考平面的投影向量的夹角;
(2)将{θ*}制成曲线,如图7所示,曲线与坐标轴的交点G、F分别为A0、θmax *,其中,为面向剪切方向大于某一视倾角的凸体面积之和与粗糙表面实际面积的比值,A0为最大可能接触面积比,定义为面向剪切方向的微元面积和粗糙表面实际面积的比值,θmax *为视倾角集合{θ*}中的最大值;
步骤10:待路径上所有的悬停点都测量完毕后,采用测量拼接的方法将多次测量结果转换到同一个坐标系下,删除重合区域,完成所有测量区域的拼合;
所述测量拼接的方法如下:
在区域内设置至少3个不共线的特征点,对两个区域分别进行三维坐标测量,得出特征点在2次测量中的坐标值;不同次测量中,特征点之间间距是相同的,利用特征点间距不变原理,进行坐标的空间变换,将2次测量结果转换到同一个坐标系下,删除重合区域,即完成了两个区域的拼接过程。
步骤11:系统反算JRC,系统将所测各基面积结构面三维形貌参数按式(1)反算出JRC值,现场随机选取岩体结构面可测量的区域测算JRC值,与反算结果比对,误差容许值以工程需要确定,其中反算公式的系统误差为±2%;
步骤12:复测,无法满足工程精度要求的,需重新测量。
Claims (10)
1.一种无人机测量大型岩体结构面三维形貌参数的装置,其特征在于,包括测量模块、动态维稳模块、动力模块和计算机分析模块;
所述测量模块包括:两个工业摄像头和一个数字光栅投影设备;
所述动态维稳模块包括:下置单云台和水平连杆;
所述动力模块由无人机的动力系统构成;
所述计算机分析模块内置数据采集卡和基于反算JRC的数据校准模块;
所述动态维稳模块将下置单云台与无人机的底部中心位置连接,下置单云台的另一端与水平连杆的中心位置连接,并在水平连杆的中心位置悬挂数字光栅投影设备,在水平连杆的两端分别通过摄像头旋动机械转盘悬挂一个工业摄像头。
2.根据权利要求1所述的无人机测量大型岩体结构面三维形貌参数的装置,其特征在于,所述数字光栅投影设备的角度不可调节,两个工业摄像头的机身角度可按档位机械调节。
3.根据权利要求1所述的无人机测量大型岩体结构面三维形貌参数的装置,其特征在于,所述下置单云台的水平旋转角度为0°-180°,水平旋转速度为3°-10°/s,垂直旋转速度为4°/s。
4.根据权利要求1所述的无人机测量大型岩体结构面三维形貌参数的装置,其特征在于,所述工业摄像头具有无线传输能力,将捕捉到岩体结构表面反射的光栅条纹信息通过无线通信发送到计算机分析模块中。
5.采用权利要求1至4中任意一项所述的无人机测量大型岩体结构面三维形貌参数的装置进行大型岩体结构面三维形貌参数测量的方法,其特征在于,包括如下步骤:
步骤1:标定空间坐标,由无人机内置GPS系统定位空间基点;
步骤2:标定内部参数,将标定块型号和误差限值输入系统,将左右摄像头拍摄的标定块特征点图像传输至系统,系统通过摄影测量学原理自动确定测量装置内部参数,并建立点云空间坐标和实际空间坐标的转换矩阵M;
步骤3:确定测量范围,依据裸露岩体整体倾向的一致性划分测量区域;
步骤4:在保证一次测量为完整基面积的前提下,以测量移动次数最少为原则,依据测量范围和重叠等级确定基测量面积;
步骤5:对无人机测量大型岩体结构面三维形貌参数的装置进行测量前的校准和摄像头的角度调整;
步骤6:路径设定,利用辅助无人机测绘岩体空间图像,输入基面积和重叠等级,系统根据基测量面积对岩石表面进行分块,并生成路径方案和悬停点坐标;
步骤7:设置无人机按设定路径自动巡航,或人工指引飞行,并开启测量;
步骤8:数字光栅投影将不同条纹间隔密度的光栅重叠,向岩体结构面表面出射一系列连续的宽度不同的光栅条纹,实现灰度二进制编码;接触岩体表面后因受到表面非规则凸体的影响,条纹发生变形,被左右两台工业摄像机捕捉并记录下来,并将空间位置坐标、拍摄俯仰角集成至图像信息中存储在图片中,通过无线通信发送到计算机分析模块;
步骤9:将左右摄像机的集成图像信息进行共点匹配,利用左右摄像机像点坐标和转换矩阵计算被测表面上点的实际三维坐标,通过计算相关任意不共线点所构成的空间曲面的法向向量,利用最小二乘法解算实现三角剖分,获取唯一解数据,再利用表面形貌参数的计算公式,获得岩石结构面表面的形貌参数;
步骤10:待路径上所有的悬停点都测量完毕后,采用测量拼接的方法将多次测量结果转换到同一个坐标系下,删除重合区域,完成所有测量区域的拼合;
步骤11:系统反算JRC,系统将所测各基面积结构面三维形貌参数按式(1)反算出JRC值,现场随机选取岩体结构面可测量的区域测算JRC值,与反算结果比对,误差容许值以工程需要确定,其中反算公式的系统误差为±2%;
步骤12:复测,无法满足工程精度要求的,需重新测量。
6.根据权利要求5所述的大型岩体结构面三维形貌参数测量的方法,其特征在于,所述利用左右摄像机像点坐标和转换矩阵计算被测表面上点的实际三维坐标的方法如下:
(1)拾取结构面凸体点:
投影装置将光栅条纹投射在岩体结构面表面上,经凸体变形,分别由左右摄像头捕捉,两侧摄像头分别以独立内置坐标完成测量;单个摄像头的测量原理如下:
将岩体结构面表面某一凸体与相邻凸体边界围合平面的绝对平整延伸面作为参考平面Ω,在被测结构面上选取任一点P,称为凸体拾取点,P在参考平面Ω上的投影为P’,线段PP’的长度即点P距离Ω平面的高度为h,数字光栅投影设备结构中心与点P连线延长线交参考面Ω的交点为点A,摄像机光学中心与点P连线延长线交参考平面Ω的交点为点B;
数字光栅投影设备出射光线至参考平面Ω上的点A,摄像机捕捉到A,经过结构面凸体变形作用,该光线照射到凸体表面的点P,此时摄像机捕捉到B,即经过凸体变形作用点A移动到新的位置点B,AB的距离a耦合了凸体拾取点的高程信息h(X,Y),以此表示受到了表面形状的调制;
其中,a与系统的几何参数有关,即工业摄像机光心到参考平面Ω的距离d、数字光栅投影设备结构中心与工业摄像机光学中心的距离L有关,即:
h(X,Y)=f(d,L,a) (2)
(2)通过共点匹配,确定被测点在左右摄像机中的像点坐标:
对于任意一点P,其在两台工业摄像机内置坐标系统中单独形成坐标P1(X1,Y1,Z1)和P2(X2,Y2,Z2),空间内任意点的两坐标向量符合同一个转换矩阵函数,形式如下:
(3)利用左右摄像机中的像点坐标,计算被测点实际空间坐标:
经由方程
其中,u,v为像点在摄像机坐标体系下的齐次坐标,即
u1=X1/Z1,v1=Y1/Z1,u2=X2/Z2,v2=Y2/Z2,消去Z1、Z2,利用最小二乘法解得:
其中,(X0,Y0,Z0)为所求被测点实际空间坐标,(X1,Y1)(X2,Y2)为左右摄像头像点,M为转换矩阵,左、右摄像头各存在一个转换矩阵,分别表示为:
由系统结合公式(2)自动标定。
7.根据权利要求6所述的大型岩体结构面三维形貌参数测量的方法,其特征在于,所述步骤5的过程如下:
步骤5.1:确定摄像头角度,依据基测量面积确定摄像头角度,利用机械转盘将摄像头旋至基面积指定档位,锁死;
步骤5.2:共点校准,试飞无人机,开启测量和数据传输,系统检索点云数据,将左右摄像头点云数据进行共点匹配,若无匹配项则返回步骤2,重新标定内部参数;
步骤5.3:单点校准,试飞无人机,开启测量和数据传输,系统检索点云数据,关闭测量和数据传输,将云台旋转180°,无人机定点悬停旋转180°,开启测量和数据传输,系统检索点云数据,系统计算误差,镜头传感器自动校准;
步骤5.4:GPS定位校准,飞行器自系统自动校准;
步骤5.5:悬停校准,飞行器自系统自动校准。
8.根据权利要求6所述的大型岩体结构面三维形貌参数测量的方法,其特征在于,所述步骤6中对岩石表面进行分块时相邻的两块的重叠区域不小于30%。
9.根据权利要求6所述的大型岩体结构面三维形貌参数测量的方法,其特征在于,所述利用表面形貌参数的计算公式,获得岩石结构面表面的形貌参数的过程如下:
(1)利用三角剖分将扫描的结构面形貌做精简网格划分,利用定点坐标确定凸体表面外法向向量 为在参考平面Ω的投影向量,为剪切方向向量,系统默认沿参考平面Ω间隔60°生成各向剪切方向向量,按下式计算{θ*}:
tanθ*=-tanθcosα
其中,为剪切方向向量,为凸体表面外法向向量,为参考平面法向向量,为凸体表面外法向向量在参考平面的投影,θ*为视倾角,{θ*}为结构面各个凸体视倾角的集合,θ为凸体表面外法向向量和参考平面Ω法向向量夹角,α为剪切方向向量和凸体表面外法向向量在参考平面的投影向量的夹角;
(2)将{θ*}制成曲线,曲线与坐标轴的交点G、F分别为A0、θmax *,其中,为面向剪切方向大于某一视倾角的凸体面积之和与粗糙表面实际面积的比值,A0为最大可能接触面积比,定义为面向剪切方向的微元面积和粗糙表面实际面积的比值,θmax *为视倾角集合{θ*}中的最大值;
10.根据权利要求6所述的大型岩体结构面三维形貌参数测量的方法,其特征在于,所述测量拼接的方法如下:
在区域内设置至少3个不共线的特征点,对两个区域分别进行三维坐标测量,得出特征点在2次测量中的坐标值;不同次测量中,特征点之间间距是相同的,利用特征点间距不变原理,进行坐标的空间变换,将2次测量结果转换到同一个坐标系下,删除重合区域,即完成了两个区域的拼接过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011252160.7A CN112344877B (zh) | 2020-11-11 | 2020-11-11 | 无人机测量大型岩体结构面三维形貌参数的装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011252160.7A CN112344877B (zh) | 2020-11-11 | 2020-11-11 | 无人机测量大型岩体结构面三维形貌参数的装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112344877A true CN112344877A (zh) | 2021-02-09 |
CN112344877B CN112344877B (zh) | 2022-02-01 |
Family
ID=74363313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011252160.7A Active CN112344877B (zh) | 2020-11-11 | 2020-11-11 | 无人机测量大型岩体结构面三维形貌参数的装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112344877B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113533042A (zh) * | 2021-07-07 | 2021-10-22 | 北京科技大学 | 一种表征岩石应力与破裂的综合性指标计算方法及应用 |
CN113946982A (zh) * | 2021-11-19 | 2022-01-18 | 中铁二院工程集团有限责任公司 | 一种获取危岩体地形剖面的方法 |
CN114895367A (zh) * | 2022-04-27 | 2022-08-12 | 清华大学 | 岩体产状信息测量方法 |
CN118314192A (zh) * | 2024-06-07 | 2024-07-09 | 中国地质大学(武汉) | 三维激光扫描合理采样间距确定方法、设备及存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101033950A (zh) * | 2007-04-04 | 2007-09-12 | 同济大学 | 岩石表面三维形貌的测量方法 |
CN101038156A (zh) * | 2007-04-05 | 2007-09-19 | 同济大学 | 便携式岩石表面三维形貌的测量系统 |
CN101038158A (zh) * | 2007-04-05 | 2007-09-19 | 同济大学 | 大尺寸岩石表面三维形貌的测量方法 |
CN104006803A (zh) * | 2014-06-20 | 2014-08-27 | 中国人民解放军国防科学技术大学 | 自旋稳定空间飞行器转动运动参数的摄像测量方法 |
CN104501779A (zh) * | 2015-01-09 | 2015-04-08 | 中国人民解放军63961部队 | 基于多站测量的无人机高精度目标定位方法 |
CN206988775U (zh) * | 2017-07-05 | 2018-02-09 | 四川建筑职业技术学院 | 一种能够实现全景摄像的多机位稳定云台 |
CN108474658A (zh) * | 2017-06-16 | 2018-08-31 | 深圳市大疆创新科技有限公司 | 地面形态检测方法及系统、无人机降落方法和无人机 |
CN109900713A (zh) * | 2019-04-17 | 2019-06-18 | 中国人民解放军国防科技大学 | 摄像引导的无人机风电叶片缺陷动态检测系统及其方法 |
CN209258413U (zh) * | 2018-10-15 | 2019-08-16 | 广东电网有限责任公司 | 一种无人机单云台系统 |
CN110440763A (zh) * | 2019-08-15 | 2019-11-12 | 华南机械制造有限公司 | 一种航测装置 |
-
2020
- 2020-11-11 CN CN202011252160.7A patent/CN112344877B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101033950A (zh) * | 2007-04-04 | 2007-09-12 | 同济大学 | 岩石表面三维形貌的测量方法 |
CN101038156A (zh) * | 2007-04-05 | 2007-09-19 | 同济大学 | 便携式岩石表面三维形貌的测量系统 |
CN101038158A (zh) * | 2007-04-05 | 2007-09-19 | 同济大学 | 大尺寸岩石表面三维形貌的测量方法 |
CN104006803A (zh) * | 2014-06-20 | 2014-08-27 | 中国人民解放军国防科学技术大学 | 自旋稳定空间飞行器转动运动参数的摄像测量方法 |
CN104501779A (zh) * | 2015-01-09 | 2015-04-08 | 中国人民解放军63961部队 | 基于多站测量的无人机高精度目标定位方法 |
CN108474658A (zh) * | 2017-06-16 | 2018-08-31 | 深圳市大疆创新科技有限公司 | 地面形态检测方法及系统、无人机降落方法和无人机 |
CN206988775U (zh) * | 2017-07-05 | 2018-02-09 | 四川建筑职业技术学院 | 一种能够实现全景摄像的多机位稳定云台 |
CN209258413U (zh) * | 2018-10-15 | 2019-08-16 | 广东电网有限责任公司 | 一种无人机单云台系统 |
CN109900713A (zh) * | 2019-04-17 | 2019-06-18 | 中国人民解放军国防科技大学 | 摄像引导的无人机风电叶片缺陷动态检测系统及其方法 |
CN110440763A (zh) * | 2019-08-15 | 2019-11-12 | 华南机械制造有限公司 | 一种航测装置 |
Non-Patent Citations (2)
Title |
---|
朱承金 等: "无人机摄影测量边坡结构面及块体稳定分析", 《东北大学学报 (自然科学版)》 * |
陈曦 等: "基于 Grasselli 形貌参数的岩石节理初始剪胀角新模型", 《岩石力学与工程学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113533042A (zh) * | 2021-07-07 | 2021-10-22 | 北京科技大学 | 一种表征岩石应力与破裂的综合性指标计算方法及应用 |
CN113533042B (zh) * | 2021-07-07 | 2022-04-05 | 北京科技大学 | 一种表征岩石应力与破裂的综合性指标计算方法及应用 |
CN113946982A (zh) * | 2021-11-19 | 2022-01-18 | 中铁二院工程集团有限责任公司 | 一种获取危岩体地形剖面的方法 |
CN113946982B (zh) * | 2021-11-19 | 2023-05-23 | 中铁二院工程集团有限责任公司 | 一种获取危岩体地形剖面的方法 |
CN114895367A (zh) * | 2022-04-27 | 2022-08-12 | 清华大学 | 岩体产状信息测量方法 |
CN114895367B (zh) * | 2022-04-27 | 2023-11-21 | 清华大学 | 岩体产状信息测量方法 |
CN118314192A (zh) * | 2024-06-07 | 2024-07-09 | 中国地质大学(武汉) | 三维激光扫描合理采样间距确定方法、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112344877B (zh) | 2022-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112344877B (zh) | 无人机测量大型岩体结构面三维形貌参数的装置及方法 | |
CN111473739B (zh) | 一种基于视频监控的隧道塌方区围岩变形实时监测方法 | |
CN113607135B (zh) | 一种用于路桥施工领域的无人机倾斜摄影测量方法 | |
CN109764858B (zh) | 一种基于单目相机的摄影测量方法及系统 | |
CN110780307B (zh) | 基于电瓶车车载式激光点云移动测量系统获取道路横断面的方法 | |
CN103925872A (zh) | 一种用于获取目标物空间分布的激光扫描测量装置和方法 | |
CN110186383B (zh) | 基于目标点仰角变化的单目相机挠度测量方法 | |
CN113739720B (zh) | 一种融合声学与光学方法的水下高精度测量和缺陷检测方法 | |
CN111307046B (zh) | 基于半球图像的树高测量方法 | |
CN113034674B (zh) | 一种利用多设备配合的施工安全检查的方法和装置 | |
CN108830317B (zh) | 基于数字摄影测量的露天矿山边坡岩体节理产状快速精细取值方法 | |
CN111733673A (zh) | 基于移动平台的路面不平度非接触测量方法及其装置 | |
CN114018167A (zh) | 一种基于单目三维视觉的桥梁挠度测量方法 | |
CN109631841B (zh) | 一种基于激光投影的快速道路横断面测量方法及装置 | |
CN116839539A (zh) | 一种用于隧道工程建设的自动监控量测系统 | |
CN116859361A (zh) | 用于点云数据的精确匹配方法、激光雷达及系统 | |
CN115601517A (zh) | 岩体结构面信息采集方法、装置、电子设备及存储介质 | |
CN115393537A (zh) | 输电通道三维可视化建模的精准性评估系统及方法 | |
CN205300519U (zh) | 铁塔变形测量装置 | |
CN1405736A (zh) | 基于近景摄影测量的脑电电极空间定位方法 | |
CN114266835A (zh) | 一种非量测相机的变形监测控制方法与系统 | |
CN106091930B (zh) | 一种基于双相机测量系统及结构光传感器的实时在线测量方法 | |
CN114660586B (zh) | 一种雷达弹目偏差测量报告自动生成系统及方法 | |
Fa et al. | Improving Vibration Monitoring of Structures Using Theodolites with Built-In Image Sensors | |
CN113610925B (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 |