【技术实现步骤摘要】
一种基于前视全景钻孔摄像的虚拟岩芯成像设备及方法
[0001]本专利技术涉及工程地质勘察领域,属于地质钻孔分析技术,具体涉及一种基于前视钻孔摄像的虚拟岩芯建模技术。
技术介绍
[0002]在施工过程中,常常需要了解工程所处地区岩体的力学特性及其稳定性信息。而岩体中结构面的类别属性、分布疏密和发育情况等特性是影响岩体结构稳定性、力学特性和工程特性的关键因素。许多重大工程事故的发生往往是由于在施工初期,不清晰岩体的内部构造,对存在的裂隙、破裂面和软弱夹层等薄弱结构面勘察不准导致的。因此,如何准确的掌握和描述岩体内部结构和力学特性是工程中不断精进的重要课题。
[0003]前视全景钻孔摄像系统是一项常用的钻孔勘察手段。其通过安装在探头顶端的CCD摄像头采集孔壁的影像数据,传输到监视器上供工程人员实时监测,并输出到记录仪上存档以备后续的调用和分析。工程人员可通过监测端/记录端观测有效光照区域内的孔壁视频录像,并以视频回放或屏幕截图的方式对岩体进行分析。该系统适用范围广,对钻孔质量要求较低,且设备轻巧,价格低廉,但缺点是图像数据量大,无法提供钻孔的全景展开图和“虚拟岩芯”模型,导致对岩体的定量分析能力较弱。
技术实现思路
[0004]本专利技术所要解决的技术问题是:优化现有设备的定量分析能力,实现钻孔岩壁的全景展开和“虚拟岩芯”构建,通过耦合图像处理、视觉重建和点云建模等领域的方法,实现从孔壁的影像数据中快速拼接融合出岩芯的全景展开图,并且实现岩芯三维模型的“升维”反演。
[0005]为了解决以上 ...
【技术保护点】
【技术特征摘要】
1.一种基于前视全景钻孔摄像的虚拟岩芯成像设备,其特征在于,包括:孔口支架(1)、深度编码器(2)、电动绞盘(3)、视频传输线(4)、数据传输线(5)、控制主机(6)、测量探头(7)、计算机系统(8);其中,所述孔口支架(1),放置于在钻孔的孔口上方;所述电动绞盘(3)包含步进电机、线缆圆盘,所述步进电机用于驱动线缆圆盘转动,提供平稳的下放或提拉速度;其中线缆圆盘通过一固定结构与孔口支架(1)相连接;所述深度编码器(2),安装在孔口支架(1)上,用于记录实时的测量探头(7)的深度信息,以及视频传输线(4)的收放长度;所述视频传输线(4),用于传输测量探头录制的影像数据,其一端与测量探头(7)连接,另一端通过深度编码器(2)后缠绕在线缆圆盘上;所述数据传输线(5),用于传输带有深度值的影像数据,其一端连接深度编码器(2),另一端通过电动绞盘(3)依次连接控制主机(6)以及计算机系统(8);所述测量探头(7),用于录制钻孔岩壁影像数据,其通过视频传输线(4)被电动绞盘(3)控制于钻孔中匀速下放至孔底后再提拉至孔口;所述控制主机(6),根据实时接收的影像数据进行抽帧形成图像序列,并将图像序列与探头深度信息进行时间匹配,生成带有探头深度信息的钻孔图像序列;计算机系统(8),用于根据带有探头深度信息的钻孔图像序列对钻孔岩壁进行三维重建与钻孔岩壁图像拼接,构建出虚拟岩芯。2.一种基于权利要求1所述虚拟岩芯成像设备的成像方法,其特征在于,具体步骤如下:步骤S1、将测量探头与孔口轴心对齐,操作控制主机启动设备,归零深度编码器,并检查控制主机是否有钻孔图像的实时反馈;步骤S2、调试成功后,启动步进电机,使测量探头匀速下降到孔底;稳定一段时间,启动步进电机上拉测量探头并实时传输钻孔岩壁影像数据,深度编码器将测量探头在钻孔中的实时深度传输给控制主机;控制主机对影像数据进行抽帧形成图像序列,并将图像序列与探头深度进行时间上的匹配,生成带有探头深度信息的钻孔图像序列;步骤S3、针对钻孔图像序列,采用霍夫圆检测法识别钻孔中心和岩环的最小半径R
min
;步骤S4、通过掩码序列对钻孔图像序列进行划分,生成有效图像序列;步骤S5、分别提取有效图像序列的DSP
‑
SIFT特征和SuperPoint特征,并进行顺序匹配;步骤S6、分别将所有匹配好的DSP
‑
SIFT特征和SuperPoint特征投影到图像上,并通过匹配对生成图像序列的特征关联网络V=(I,T,L),其中I为有效图像,T为特征点轨迹,L为轨迹的边,同时采用测地线一致性对特征关联网络进行误匹配剔除;步骤S7、通过初始搭建、深度计算、捆绑优化模块构建岩芯模型的归一化模型;步骤S8、使用统计去噪算法来去除模型噪点,并基于有效图像序列的深度值实现归一化模型的尺度恢复;步骤S9、采用FPFH特征进行粗匹配,然后采用ICP算法实现岩芯三维模型的精确耦合;步骤S10、通过基于随机一致性采样算法对岩芯三维模型进行圆柱拟合,获得钻孔轴线以及圆柱半径r,计算岩芯三维模型上3D点到圆柱轴线的垂直距离并与圆柱半径r相减,获得3D点轴向深度信息;
步骤S11、对有效图像序列进行抽稀和筛选,生成环形图像序列,通过矩形展开图的笛卡尔坐标系与环形图像序列的极坐标系之间的坐标变换关系,构建两者的像素映射关系;并结合双线性插值算法进行图像展开;步骤S12、通过模板匹配方法进行矩形展开图的拼接匹配:首先设置上一帧图像为参考帧,当前帧图像为目标帧,并在参考帧图像中选择模板图像与当前帧进行模板匹配;利用模板匹配算法获得两帧图像的匹配参数,以进行矩形展开图的拼接匹配,然后,将拼接生成的新图像再与下一帧进行模板匹配与拼接,最后拼接成为一幅完整的全景孔壁展开图;步骤S13、采用相互赋值的方式将岩芯三维模型和全景孔壁展开图相互关联:根据岩芯三维模型的轴向深度信息,将岩芯三维模型的每个3D点映射到全景孔壁展开图上,形成具有高程信息的点云模型,并根据全景孔壁展开图中对应位置的RGB信息对3D点赋予颜色,生成实景孔壁展开图;最后,将3D点的RGB信息赋予岩芯三维模型,生成岩芯实景三维模型。3.根据权利要求2所述的成像方法,其特征在于,步骤S3中所采用的霍夫圆检测法具体采用canny边缘检测算法处理钻孔图像,包括:(1)高斯滤波:采用双维度同标准差的高斯核,并将其与图像进行离散卷积:其中:G(x,y)为高斯核;σ2为预选方差;(2)像素梯度运算:采用Sobel算子进行像素梯度转换,并生成图像的梯度强度矩阵;Sobel算子是两个3
×
3的矩阵,分别为S
x
和S
y
,其中S
x
用于计算图像x方向的像素梯度矩阵g
x
,S
y
用于计算图像y方向的像素梯度矩阵g
y
,具体为:,具体为:,具体为:其中I
gray
代表钻孔图像序列的灰度矩阵,g
xy
(i,j)代表梯度强度矩阵,对x的梯度矩阵g
x
(i,j)和y的梯度矩阵g
y
(i,j)取二范数,i,j分别代表两个不同的像素点;(3)非极大值像素梯度抑制:采用梯度幅值的非极大值抑制方法NMS将当前像素梯度强度g
xy
(i,j)与正负梯度方向的相邻强度比较,抑制非边缘点,消除边缘检测的杂散响应。(4)阈值滞后处理:采用双阈值进一步对边缘点进行分类,低于低阈值的像素点被重新抑制;高于高阈值的像素点视为强边缘点;介于两者之间的定义为弱边缘,等待进一步筛分;(5)孤立弱边缘抑制:对弱边缘采用邻域索引的方式进行筛除,若邻域内含有强边缘点则保留;根据边缘点的梯度方向,根据投影不变量实现边缘线段的椭圆检测;通过椭圆拟合,求解最佳圆心,并根据点到圆心的距离,计算出最小岩环半径。4.根据权利要求2所述的成像方法,其特征在于,步骤S5中,DSP
‑
SIFT特征计算方法如下:
H
DSP
(θ|I)[p]=∫h
SIFT
(...
【专利技术属性】
技术研发人员:何磊,楚宇扬,徐浩若,卫海,肖怀广,王睿华,郑彦龙,
申请(专利权)人:东南大学,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。