基于直线段检测和LT描述符的矿井图像线特征匹配算法

朱代先, 秋强, 孔浩然, 胡其胜, 刘树林

朱代先,秋强,孔浩然,等. 基于直线段检测和LT描述符的矿井图像线特征匹配算法[J]. 工矿自动化,2024,50(2):72-82. DOI: 10.13272/j.issn.1671-251x.2023090045
引用本文: 朱代先,秋强,孔浩然,等. 基于直线段检测和LT描述符的矿井图像线特征匹配算法[J]. 工矿自动化,2024,50(2):72-82. DOI: 10.13272/j.issn.1671-251x.2023090045
ZHU Daixian, QIU Qiang, KONG Haoran, et al. A line feature matching algorithm for mine images based on line segment detection and LT descriptors[J]. Journal of Mine Automation,2024,50(2):72-82. DOI: 10.13272/j.issn.1671-251x.2023090045
Citation: ZHU Daixian, QIU Qiang, KONG Haoran, et al. A line feature matching algorithm for mine images based on line segment detection and LT descriptors[J]. Journal of Mine Automation,2024,50(2):72-82. DOI: 10.13272/j.issn.1671-251x.2023090045

基于直线段检测和LT描述符的矿井图像线特征匹配算法

基金项目: 陕西省重点研发计划项目(2021GY-338);西安市碑林区科技计划项目(GX2333)。
详细信息
    作者简介:

    朱代先(1970—),男,安徽安庆人,副教授,博士研究生,主要从事智能机器人、嵌入式系统方面的研究工作,E-mail:zhudaixian@xust.edu.cn

    通讯作者:

    秋强(1996—),男,陕西咸阳人,硕士研究生,研究方向为计算机视觉,E-mail:qiuqiang1168@163.com

  • 中图分类号: TD67

A line feature matching algorithm for mine images based on line segment detection and LT descriptors

  • 摘要: 图像匹配是同步定位与地图构建(SLAM)技术中极为重要的一环,用于根据图像之间的变换关系确定相机位姿。基于线特征的图像匹配方法具有较强的鲁棒性和抗噪能力,更加适用于井下图像匹配,基于深度学习的线描述符对线段遮挡等场景具有较高的鲁棒性,性能优于传统描述符,但卷积神经网络架构的描述符将可变长度线段抽象为固定维进行描述,不利于线段长度及视差变化较大图像的匹配。针对上述问题,提出一种基于直线段检测和线描述符的矿井图像线特征匹配算法。在频域利用单参数同态滤波降低图像的照射分量,并增强反射分量,提升亮度及对比度;在YUV空间利用对比度受限的自适应直方图均衡化(CLAHE)算法对亮度分量进行均衡,使亮度分布更加均匀;变换至RGB空间提取直线段检测(LSD)线,引入一种基于Transformer架构的LT描述符构建LSD线的特征向量,最后完成线特征匹配。实验结果表明:该算法结合了同态滤波和CLAHE算法的优点,增强后图像的亮度适中,对比度良好,灰度分布均匀,增强效果优于单参数同态滤波算法、EnlightenGAN算法;该算法提取的线特征数较原图平均提升了32.92%,在不同相似纹理占比、不同程度旋转与平移变化的井下图像匹配中鲁棒性好,平均正确匹配数为61.75对,平均精度为86.83%,优于线二进制描述符(LBD)算法、LBD_NNDR算法、LT算法,能够满足矿井图像稳健匹配的需求。
    Abstract: Image matching is an extremely important part of simultaneous localization and mapping (SLAM) technology. It is used to determine camera position and posture based on the transformation relationship between images. The image matching method based on line features has strong robustness and noise resistance, making it more suitable for underground image matching. The line descriptors based on deep learning have high robustness to scenes such as line segment occlusion, and their performance is better than traditional descriptors. However, the descriptors of convolutional neural network architecture abstract variable length line segments into fixed dimensions for description, which is not conducive to matching images with large changes in line segment length and parallax. In order to solve the above problems, a line feature matching algorithm for mine images based on line segment detection and line transformers (LT) is proposed. The algorithm uses single parameter homomorphic filtering in the frequency domain to reduce the lighting component of the image, enhance the reflection component, and improve brightness and contrast. The algorithm uses contrast limited adaptive histogram equalization (CLAHE) algorithm in YUV space to balance brightness components and make brightness distribution more even. The algorithm transforms to RGB space to extract line segment detection (LSD) lines. A LT descriptor based on Transformer architecture is introduced to construct the feature vector of LSD lines, and finally complete line feature matching. The experimental results show that the algorithm combines the advantages of homomorphic filtering and CLAHE algorithm. After image enhancement, the brightness of the image is moderate, the contrast is good, the grayscale distribution is even. The enhancement effect is better than the single parameter homomorphic filtering algorithm and EnlightenGAN algorithm. The number of line features extracted by this algorithm has increased by an average of 32.92% compared to the original image. It has good robustness in matching underground images with different proportions of similar textures, varying degrees of rotation and translation changes. The average correct matching number is 61.75 pairs, with an average precision of 86.83%. It is superior to the line binary descriptor (LBD) algorithm, LBD_NNDR algorithm, and LT algorithm. It can meet the requirements of robust matching of mine images.
  • 准确获取结构面信息是岩质边坡稳定性分析的重要前提。传统的结构面接触式测量方法有测线法和测窗法,这2种方法操作简单,但对测量者的专业技术水平要求较高,易受测量者主观影响,对于复杂结构面需要进行多次测量和综合分析,工作量大且效率低;尤其对于高陡边坡,测量者无法直接进行人工编录,作业危险性大,难以获取全面的结构面信息[1-2]。因此,高效准确地采集岩质高边坡结构面信息仍是一个关键难题。

    针对传统接触式测量技术存在的各种问题,研究者提出了非接触式测量技术,包括三维激光扫描技术和无人机摄影测量技术[3]。三维激光扫描技术对测量环境要求较高且相关设备价格昂贵。无人机摄影测量技术以无人机为云台,通过机身搭载高精度数码相机和传感器对被测物体拍摄高清影像,然后结合计算机图像处理软件进行后处理与分析[4],具有机动灵活、性价比高、维护简单和适用性强等优势,因而逐渐被应用于数字化地形测绘、建筑物测量等领域[5]

    近年来,研究者将无人机摄影测量技术应用于岩质高边坡结构面识别中。赵明宇等[6]通过无人机和全站仪对边坡进行联合测量,获取了边坡结构面迹线和产状信息。熊开治等[7]开展了基于无人机摄影测量技术的随机裂隙网络模型构建方法研究。李水清等[8]基于无人机和“三点法”提出了一种岩体结构面产状统计的半自动方法,克服了传统方法的局限性。虽然无人机摄影测量技术可以获取岩质高边坡表面的高清影像,但是缺少高效准确的影像后处理方法,导致结构面识别仍然存在准确性差和效率低等缺点。一些学者尝试采用不同数据处理算法进行结构面识别。沙鹏等[9]在采用无人机航测获取边坡图像的基础上,结合运动恢复结构法(Structure from Motion,SfM)获取边坡岩体点云模型,进而建立包含离散裂隙网络的实际边坡三维离散元模型。贾曙光等[10]使用无人机采集北京市某采石场地形数据,结合计算机视觉原理与SfM方法重建了三维模型测量结构面产状。张恺等[11]采用无人机和SfM法对浙江某露天采石场边坡进行三维重构,提取岩体结构面产状。Liu Yongqiang等[12]采用无人机拍摄高分辨率影像,开发了一种同时考虑产状、迹长和起伏度的结构面分组新方法。然而,当前算法在构建高密度三维点云模型、精准识别与拟合复杂结构面存在局限性,且没有充分考虑识别结构面产状信息统计特征的不确定性,亟需构建高效准确的融合高边坡摄影、后处理算法及统计分析的一体化结构面识别与产状信息采集方法。

    本文以江西省南昌市某露天矿高边坡为研究背景,提出一种岩质高边坡结构面识别及产状统计信息采集方法。采用无人机摄影测量技术和计算机视觉原理,使用无人机对边坡进行航测并获取基础影像数据;采用Context Capture系列图像处理软件、Cloud Compare、DSE(Discontinuity Set Extractor)和多种算法进行边坡结构面识别,进而提取结构面产状信息,进行产状统计特征分析。

    露天矿边坡现场如图1所示。通过现场踏勘得知,边坡长约200 m,坡顶坡底落差约100 m,边坡出露岩体节理裂隙发育。将研究区岩体划分为上层和下层2个部分,每层岩体高陡,难以直接采用传统的接触式测量方法进行结构面产状测量。选择地点开阔、人员稀少且不受磁场干扰影响的地区作为无人机飞行起降点,并通过云台设置覆盖整个研究区的飞行路线,进行无人机摄影测量工作。

    图  1  露天矿边坡现场
    Figure  1.  Open pit slope site

    考虑到研究区地形复杂、障碍物多等因素,选择体积小、控制性能好、操作简单的Phantom 4 Pro V2.0四旋翼无人机进行测量。无人机系统主要包括飞行器、云台、视觉系统、红外感知系统、相机、遥控器、智能飞行电池和存储器等组件。

    采用无人机获取高分辨率边坡影像后,需要选择合适的计算机图像处理软件基于图像二维坐标构建三维点云模型。目前,常用的无人机航测影像后处理软件主要包括Pix4D、大疆智图、Photo Scan、Context Capture、Open Drone Map等。其中, Context Capture具有快速、简单、全自动等优势,能够在无人工干预情况下,从简单连续的二维影像还原出真实的三维实景模型和超高密度点云模型。此外,该软件还具有良好的输入数据兼容性和广泛的输出数据格式,能够满足绝大多数的工作需要。因此,本文选择Context Capture软件进行无人机拍摄影像的后处理。

    由于无人机配备有惯性测量装置传感器及全球定位系统(Global Positioning System,GPS)与全球轨道导航卫星系统(Global Orbiting Navigation Satellite System,GLONASS),在飞行作业中拍摄的影像均保留了传感器信息、相机参数、照片姿态参数、照片位置参数等元数据,从而保证能够输出色彩真实、几何形态及细节精细的三维模型。Context Capture构建三维点云模型的流程如图2所示。生成的边坡上层和下层岩体高密度三维点云模型分别如图3图4所示,可看出模型与边坡实物图吻合度较高。

    图  2  三维点云模型构建流程
    Figure  2.  Construction process of 3D point cloud model
    图  3  上层岩体三维点云模型
    Figure  3.  3D point cloud model of rock mass at the upper part of the slope
    图  4  下层岩体三维点云模型
    Figure  4.  3D point cloud model of rock mass at the lower part of the slope

    在完成原始点云数据重采样后,采用A. J. Riquelme等[13]基于Matlab编写的开源软件DSE进行结构面识别,流程如图5所示。

    图  5  结构面识别流程
    Figure  5.  Identification process of structural plane

    1) 将带有空间点坐标(XYZ)和颜色信息的重采样点云数据输入到结构面提取程序中。

    2) 设置K近邻(K-Nearest Neighbor,KNN)算法中需要的参数−近邻点数量K和偏差参数${\eta _{\max }}$,并根据确定好的参数寻找共面点。

    3) 在搜寻到所有共面点子集后,分别进行最优平面拟合,计算拟合平面的法向量并作为对应目标点的法向量。

    4) 在法向量密度图的基础上,采用Riquelme算法计算极点密度,并确定主极点。

    5) 采用基于密度的聚类(Density-Based Spatial Clustering of Applications with Noise,DBSCAN)算法进行聚类分析和结构面识别,并通过最小二乘法计算结构面的平面方程。

    6) 输出三维点云及结构面数据。

    在边坡结构面识别过程中,寻找特征向量是进行后续结构面分组与提取的前提。通常选择点法向量作为特征向量,但是每个单独的点无法计算法向量,需要利用每个点和其相似点构成平面的法向量作为该点的法向量。KNN是一种非参数、有监督的学习分类器,通过邻近度指标对单个数据点进行分类或预测[14]。KNN算法主要包括近邻点距离度量法和确定近邻点数量法。① 近邻点距离度量法:为了确定最接近给定目标点的近邻点,需要计算该点与其他数据点之间的距离,以该点为球心,将所有满足给定距离的相邻点都归属于同一个点集。② 确定近邻点数量法:规定目标点需要寻找的近邻点数量为一定值,从目标点出发,将距离目标点最近的K个点归属于同一个点集。

    由于点云密度存在不均匀性,采用近邻点距离度量法进行搜索时,可能会因为边界处包含错误点而形成错误平面,导致后续计算流程出现错误[15]。而选择确定近邻点数量法进行搜索时,可通过适当缩小搜索半径降低结果出错率。因此,本文选择确定近邻点数量法构建相似点集。

    需要指出的是,设置合适的近邻点数量K是保证KNN算法性能的关键。如果K值过小(如K < 15),那么在计算极点密度过程中会保留噪声;若K值过大(如K > 30),则会降低结果的准确性。同样,对于偏差参数${\eta _{\max }}$,如果${\eta _{\max }} < 15{\text{%}} $,会导致符合的共面点被丢弃;而当${\eta _{\max }} > 25{\text{%}} $时,则会产生边缘点。根据A. C. Rencher和W. F. Christensen[16]的研究,如果1组主成分的方差贡献率不小于80%,则可认为其能够正确地表示数据。因此,本文中取K = 30,${\eta _{\max }}{\text{ = }}20{\text{%} }$。

    确定点法向量后,采用Riquelme算法计算极点的密度函数。该方法为点云中每个点指定主方向,但由于曲面的读取误差和极点分散等问题,使得局部极大值极点中只有少数点才是主极点。为此,通过设置参数${n_{\mathrm{p}}}$和$\gamma $来寻找局部最大值点并作为主极点,${n_{\mathrm{p}}}$为可建立结构面组的最大数量,$\gamma $为点法向量与指定主极点法向量之间的最大允许角度。本文设置${n_{\mathrm{p}}}{\text{ = }}10$,$\gamma $取20~30°。

    为了找到每个结构面组的数据子集并实现结构面识别,本文采用文献[17]提出的DBSCAN算法进行聚类分析。文献[18]证明了该算法在点云数据处理中的实用性。该算法需要设置2个输入参数:ε和 min-pts。ε 表示可视为相邻点的两点间最大距离,一般根据点云密度确定;min-pts为DBSCAN算法中的阈值参数,即在1个目标点的邻域内至少包含的点数。min-pts过大会导致碎片化,min-pts过小会导致子集连接过多。本文设置ε为4,min-pts为5。

    聚类分析原理:在数据集中任选1个初始点,如果在距离该点ε 距离内有超过min-pts个点,则认为这些点都属于同一个聚类。以所有新添加的点为起始点,搜索距离该点ε 距离内是否也有超过min-pts个点。如果有超过min-pts个点,则把这些点也增加到与初始点同一聚类中,依此循环,直至没有符合条件的可添加点为止,便可得到1个聚类的所有子集。以该聚类外的任一点为初始点重复上述过程,直到所有的聚类分组都被找到。

    完成聚类分析后,需要计算结构面的平面方程。设结构面组i的点集为${R_i}$,属于该组结构面聚类j的点集为${R_{ij}}$,寻找${R_{ij}}$的最佳拟合平面,其平面方程为

    $$ {A_{ij}}x + {B_{ij}}y + {C_{ij}}{\textit{z}} + {D_{ij}} = 0 $$ (1)

    式中:($x$,$y$,${\textit{z}}$)为点数据的坐标信息;${A_{ij}}$,${B_{ij}}$,${C_{ij}}$分别为平面单位法向量在xyz方向的分量;${D_{ij}}$为原点到平面的垂直距离。

    经过Context Capture软件处理获得的点云密度通常很大,直接用于结构面识别需要大量计算,时间成本大,效率低。为此,对原始点云数据进行前处理,在保证计算结果准确性的前提下,最大限度地提高计算效率。Cloud Compare软件是一款开源的三维点云编辑与处理软件,支持许多点云处理算法(包括点云配准、重采样、统计计算、自动分割等),并包含大量插件可供选择。因此,本文选择Cloud Compare软件进行点云数据重采样处理。重采样方法包括空间距离重采样、随机重采样、八叉树重采样3种。综合比较,本文选择空间距离重采样方法。由于研究区域范围较大,原始点云数据过于庞大,为提高计算效率,在Cloud Compare软件中进行重采样时,设置两点间的最小距离为5 cm。

    为了验证边坡结构面识别结果的准确性,将边坡结构面识别结果与野外实地勘测的结构面信息进行对比分析。

    图  6  上层岩体法向量密度
    Figure  6.  Normal vector density of rock mass at the upper part of the slope

    分别对图1中高边坡的上层岩体和下层岩体进行结构面识别,得到岩体法向量密度图和三维点云,如图6图9所示。边坡上层岩体共找到6个主极点J1—J6(倾向/倾角):J1(156.26°/66.07°),J2(177.52°/57.46°),J3(147.54°/34.49°),J4(197.45°/82.27°),J5(127.88°/78.31°),J6(221.05°/88.77°)。边坡下层岩体共找到5个主极点J7—J11:J7(155.00°/43.01°),J8(177.51°/57.45°),J9(135.02°/75.48°),J10(209.04°/72.65°),J11(354.48°/73.11°)。

    识别的三维点云数据与实际边坡结构面的对比如图10所示。可看出,用本文方法识别的结构面与边坡工程现场实际情况基本吻合。

    图  7  下层岩体法向量密度
    Figure  7.  Normal vector density of rock mass at the lower part of the slope
    图  8  上层岩体三维点云
    Figure  8.  3D point cloud of rock mass at the upper part of the slope
    图  9  下层岩体三维点云
    Figure  9.  3D point cloud map of rock mass at the lower part of the slope
    图  10  结构面识别结果与实际边坡结构面的对比
    Figure  10.  Comparison between structural plane recognition results and actual slope structural plane

    在聚类分析过程中对6组结构面的每个子聚类进行平面拟合,分别得到其平面方程参数,利用这些参数可计算得到每个结构面的倾角和走向,进而得到其倾向[19-20]。因为聚类结果中存在很多小的聚类,难以形成有效的结构面,所以设定一个阈值500,每个子聚类中的所有点数量不小于500时才在输出结果中显示。利用Python编写的计算程序得到边坡稳定性定量评价需要的结构面产状信息。

    对边坡上层岩体进行分析,得到拟合平面参数及结构面产状信息(表1),J1—J6所在结构面组中出露的结构面数量依次为70,66,103,49,55,20。通过统计分析,得到各组结构面倾角和倾向的概率分布直方图及函数拟合曲线,如图11所示,μ为均值,σ为标准差。

    表  1  上层岩体拟合平面参数及结构面产状信息
    Table  1.  Parameters of fitting plane and orientation information of structural plane of rock mass at the upper part of the slope
    主极点聚类号点数量拟合平面方程参数倾角/(°)倾向/(°)
    AijBijCijDij
    J1162 356−0.4070.813−0.416−49.92965.43153.40
    234 6810.268−0.8770.39942.26866.46163.01
    J2154 135−0.0100.860−0.510−33.92059.33179.36
    237 2590.2530.810−0.529−39.05558.04197.33
    J3162 0730.334−0.2880.898−58.57026.17130.82
    232 6840.323−0.4920.808−10.03636.07146.75
    J4188 3280.2690.947−0.177−44.65979.83195.89
    27 9730.2410.961−0.136−47.63782.17194.10
    J5137 880−0.9020.420−0.100−11.17984.26114.98
    227 7620.748−0.6200.23733.93176.31129.69
    J6120 0510.6430.766−0.020−31.94188.88220.03
    26 6600.6000.792−0.1160.26783.32217.13
    下载: 导出CSV 
    | 显示表格
    图  11  结构面组J1—J6倾角和倾向概率分布直方图与拟合曲线的比较
    Figure  11.  Comparison of probability distribution histograms and fitting curves for dip and dip angles of structural plane groups J1-J6

    边坡下层岩体拟合平面参数及结构面产状信息见表2。J7—J11所在结构面组中出露的结构面数量依次为102,92,41,42,12。同理,通过统计分析可得到各组结构面倾角和倾向的概率分布直方图、函数拟合曲线和结构面产状分布特征。结构面产状统计特征见表3。边坡结构面产状概率统计分布特征可以为高边坡概率稳定性的定量评价提供重要的数据来源。

    表  2  下层岩体拟合平面参数及结构面产状信息
    Table  2.  Parameters of fitting plane and orientation information of structural plane of rock mass at the lower part of the slope
    主极点聚类号点数拟合平面方程参数倾角/(°)倾向/(°)
    AijBijCijDij
    J7152 786−0.2970.687−0.664−13.62548.42156.60
    223 2170.317−0.6450.69622.36045.93153.80
    J8143 8800.001−0.8250.56530.00155.62179.93
    228 671−0.1970.846−0.495−26.15360.35166.89
    J9190 5170.742−0.6220.25115.80875.44129.97
    260 8420.723−0.6910.00156.81389.96133.68
    J10127 2160.4820.793−0.37215.41768.19211.28
    210 3960.2520.928−0.274−42.15274.07195.16
    J11
    17 6020.158−0.9870.02450.75388.61350.88
    21 6900.137−0.9890.04953.16287.16352.14
    下载: 导出CSV 
    | 显示表格
    表  3  边坡上层和下层岩体结构面产状统计特征
    Table  3.  Statistical characteristics of rock mass structural plane orientation at the upper and lower part of the slope
    上层岩体 下层岩体
    主极点 产状 分布特征 均值$\mu $ 标准差$\sigma $ 主极点 产状 分布特征 均值$\mu $ 标准差$\sigma $
    J1 倾角 对数正
    态分布
    68.32 2.92 J7 倾角 正态分布 43.39 4.42
    倾向 对数正
    态分布
    156.21 6.78 倾向 正态分布 154.31 8.50
    J2 倾角 正态分布 56.97 5.31 J8 倾角 对数正
    态分布
    59.33 5.36
    倾向 正态分布 184.40 8.40 倾向 对数正
    态分布
    170.74 5.73
    J3 倾角 正态分布 34.95 6.64 J9 倾角 正态分布 72.74 7.67
    倾向 正态分布 142.94 15.87 倾向 对数正
    态分布
    138.26 12.24
    J4 倾角 正态分布 81.39 2.93 J10 倾角 正态分布 72.09 16.48
    倾向 对数正
    态分布
    201.24 3.16 倾向
    J5 倾角 正态分布 77.58 4.47 J11 倾角
    倾向 正态分布 127.87 6.79 倾向
    J6 倾角 正态分布 83.22 4.97
    倾向 对数正
    态分布
    216.16 3.17
    下载: 导出CSV 
    | 显示表格

    对比分析研究区高边坡上层和下层岩体结构面的概率分布直方图及拟合结果可知,大部分结构面的函数拟合效果较好,仅有少数结构面的拟合效果较差。其可能原因主要如下:① 无人机获取的图像数量不够多,导致部分结构面细节丢失,进而增大了点云数据提取过程中的误差。② 研究区边坡下层岩体表面局部存在植被,干扰点云数据聚类分析,影响结构面产状的计算结果。③ 一些结构面组出露的结构面数量很少,无法对其进行准确拟合。

    相较于当前结构面识别研究[21-22],本文方法具有以下优势。

    1) 相较于三维激光扫描方法,本文方法具有效率高、操作简单和成本低等优势,尤其在一些人力无法到达和不利于架设观测站的情况(如水电站库岸边坡等),更能突出无人机机动灵活的优势。

    2) 相较于其他数据处理算法(K-Means算法等),DBSCAN算法不需要事先指定聚类个数,能够自动识别簇的数量和形状,同时能识别离群点,有效过滤噪声数据,使得结构面识别过程更高效、结果更准确,且无需过多人工干预。

    1) 提出了基于无人机摄影测量技术和Context Capture系列软件的边坡结构面点云数据获取方法,具有效率高、操作简单和成本低等优势,为解决高边坡结构面勘测难题提供了有效途径。

    2) 采用KNN算法和DBSCAN算法进行边坡结构面识别,实现了复杂边坡结构面的准确高效识别和产状信息快速提取,结构面识别结果与工程实际吻合。

    3) 岩质高边坡结构面识别及产状统计信息采集方法可获取高边坡结构面数量、产状信息及其统计特征,大部分结构面倾角和倾向概率分布与实测数据拟合较好,为高边坡裂隙网络模型构建及稳定性分析提供了重要数据来源。

    4) 下一步将继续开展岩质高边坡结构面网络模型构建和边坡稳定性与可靠度分析,并揭示结构面分布对边坡稳定性的影响机制,最终实现科学指导水利水电工程、公路和露天矿边坡的防护与安全运行,降低边坡支护成本和失稳破坏造成的潜在损失。

  • 图  1   基于LSD和LT描述符的矿井图像线特征匹配算法流程

    Figure  1.   Process of mine image line feature matching algorithm based on line segment detector(LSD) and line transformers(LT) descriptor

    图  2   高斯型同态滤波传递函数

    Figure  2.   Transfer function of gaussian homomorphic filtering

    图  3   单参数同态滤波传递函数

    Figure  3.   Transfer function of single parameter homomorphic filter

    图  4   CLAHE原理

    Figure  4.   Principle of contrast limited adaptive histogram equalization(CLAHE)

    图  5   基于LSD和LT描述符的矿井图像线特征匹配算法模型

    Figure  5.   Model of mine image line feature matching algorithm based on LSD and LT descriptor

    图  6   图像增强结果

    Figure  6.   Image enhancement results

    图  7   灰度直方图对比结果

    Figure  7.   Comparison results of gray histogram

    图  8   LSD提取对比

    Figure  8.   LSD extraction comparison

    图  9   实验图像

    Figure  9.   Experimental images

    图  10   图像1匹配结果

    Figure  10.   Image 1 matching results

    图  11   图像2匹配结果

    Figure  11.   Image 2 matching results

    图  12   图像3匹配结果

    Figure  12.   Image 3 matching results

    图  13   图像4匹配结果

    Figure  13.   Image 4 matching results

    图  14   平均精度统计

    Figure  14.   Average accuracy statistics

    图  15   平均正确匹配数统计

    Figure  15.   Statistics of the average number of correct matches

    表  1   LT描述符的训练参数

    Table  1   Training parameters of the LT descriptor

    参数
    学习率 0.001
    训练轮次 1 000
    图像大小 640×480
    线长度最小阈值 16
    最大Token数 21
    Token间距 8
    描述符维度 256
    注意力头数量 4
    编码器特征维数 [32 64 128 256]
    线段描述层数量 12
    前馈层内部维度 1 024
    签名网络层数 7
    Transformer编码器层数 12
    下载: 导出CSV

    表  2   图像增强结果统计

    Table  2   Statistics of image enhancement results

    算法 标准差 均值 信息熵 PSNR
    改进同态滤波算法 63.583 123.820 7.712 10.014
    EnlightenGAN算法 63.312 153.775 7.683 7.598
    本文算法 64.108 129.174 7.798 10.178
    下载: 导出CSV

    表  3   LSD线段提取数量

    Table  3   LSD line segment extraction quantity

    采集点帧数LSD线平均数/条增长率/%
    原图本文算法增强图像
    巷道48094.8125.632.49
    水房218120.5166.438.09
    工作面324138.8181.530.85
    避难硐室780184.4244.832.75
    平均值144.6192.232.92
    下载: 导出CSV

    表  4   图像属性

    Table  4   Image attributes

    图像序号 采集位置 旋转与平移程度 相似纹理占比
    1 巷道 较小 较小
    2 避难硐室 较小 较大
    3 巷道 较大 较小
    4 避难硐室 较大 较大
    下载: 导出CSV

    表  5   线特征匹配实验数据统计

    Table  5   Statistics of experimental data of line feature matching

    图像
    序号
    变化
    程度
    相似纹
    理占比
    算法 线特征数量/条 同名直
    线数/对
    正确匹
    配数/对
    匹配
    精度/%
    1 较小 较小 LBD 125 105 62 53 85.48
    LBD_NNDR 125 105 58 52 89.65
    LT 125 105 51 46 90.19
    本文算法 149 130 63 58 92.06
    2 较小 较大 LBD 180 182 50 28 56.00
    LBD_NNDR 180 182 33 25 75.76
    LT 180 182 62 52 83.87
    本文算法 327 335 123 100 85.47
    3 较大 较小 LBD 100 181 20 12 60.00
    LBD_NNDR 100 181 15 11 73.33
    LT 100 181 37 34 91.89
    本文算法 137 258 42 39 92.86
    4 较大 较大 LBD 287 252 60 9 15.00
    LBD_NNDR 287 252 18 8 44.44
    LT 287 252 55 38 69.09
    本文算法 350 385 65 50 76.92
    下载: 导出CSV
  • [1] 王国法,任世华,庞义辉,等. 煤炭工业“十三五”发展成效与“双碳”目标实施路径[J]. 煤炭科学技术,2021,49(9):1-8.

    WANG Guofa,REN Shihua,PANG Yihui,et al. Development achievements of China' s coal industry during the 13th Five-Year Plan period and implementation path of "dual carbon" target[J]. Coal Science and Technology,2021,49(9):1-8.

    [2] 程德强,钱建生,郭星歌,等. 煤矿安全生产视频AI识别关键技术研究综述[J]. 煤炭科学技术,2023,51(2):349-365.

    CHENG Deqiang,QIAN Jiansheng,GUO Xingge,et al. Review on key technologies of AI recognition for videos in coal mine[J]. Coal Science and Technology,2023,51(2):349-365.

    [3] 苗升,刘小雄,黄剑雄,等. 无人机视觉SLAM环境感知发展研究[J]. 计算机测量与控制,2021,29(8):1-6,41.

    MIAO Sheng,LIU Xiaoxiong,HUANG Jianxiong,et al. Research on development of UAV visual SLAM environment perception[J]. Computer Measurement & Control,2021,29(8):1-6,41.

    [4] 孔二伟,张亚邦,李佳悦,等. 面向煤矿井下低光照图像的增强方法[J]. 工矿自动化,2023,49(4):62-69,85.

    KONG Erwei,ZHANG Yabang,LI Jiayue,et al. An enhancement method for low light images in coal mines[J]. Journal of Mine Automation,2023,49(4):62-69,85.

    [5] 赵良玉,金瑞,朱叶青,等. 基于点线特征融合的双目惯性SLAM算法[J]. 航空学报,2022,43(3):363-377.

    ZHAO Liangyu,JIN Rui,ZHU Yeqing,et al. Stereo visual-inertial SLAM algorithm based on merge of point and line features[J]. Acta Aeronautica et Astronautica Sinica,2022,43(3):363-377.

    [6]

    ZHU Daixian,JI Kangkang,WU Dong,et al. A coupled visual and inertial measurement units method for locating and mapping in coal mine tunnel[J]. Sensors,2022,22(19). DOI: 10.3390/s22197437.

    [7] 谢晓佳. 基于点线综合特征的双目视觉SLAM方法[D]. 杭州:浙江大学,2017.

    XIE Xiaojia. Stereo visual SLAM using point and line features[D]. Hangzhou:Zhejiang University,2017.

    [8]

    WANG Wei,GAO Wei,CUI Hainan,et al. Reconstruction of lines and planes of urban buildings with angle regularization[J]. ISPRS Journal of Photogrammetry and Remote Sensing,2020,165:54-66. DOI: 10.1016/j.isprsjprs.2020.04.013

    [9] 宋佳璇,范大昭,董杨,等. 神经网络学习与灰度信息结合的跨视角影像线特征匹配算法[J]. 测绘学报,2023,52(6):990-999. DOI: 10.11947/j.AGCS.2023.20220468

    SONG Jiaxuan,FAN Dazhao,DONG Yang,et al. Line matching algorithm for cross-view images combining neural network learning with grayscale information[J]. Acta Geodaetica et Cartographica Sinica,2023,52(6):990-999. DOI: 10.11947/j.AGCS.2023.20220468

    [10]

    CHEN Min,YAN Shaohua,QIN Rongjun,et al. Hierarchical line segment matching for wide-baseline images via exploiting viewpoint robust local structure and geometric constraints[J]. ISPRS Journal of Photogrammetry and Remote Sensing,2021,181:48-66. DOI: 10.1016/j.isprsjprs.2021.09.002

    [11]

    CHEN Min,LI Wen,FANG Tong,et al. An adaptive feature region-based line segment matching method for viewpoint-changed images with discontinuous parallax and poor textures[J]. International Journal of Applied Earth Observation and Geoinformation,2023,117. DOI: 10.1016/j.jag.2023.103209.

    [12]

    LANGE M,SCHWEINFURTH F,SCHILLING A. DLD:a deep learning based line descriptor for line feature matching[C]. IEEE/RSJ International Conference on Intelligent Robots and Systems,Macau,2019:5910-5915.

    [13] 王竞雪,刘肃艳,王伟玺. 联合共线约束与匹配冗余的组直线匹配结果检核算法[J]. 测绘学报,2020,49(6):746-756.

    WANG Jingxue,LIU Suyan,WANG Weixi. A checking algorithm for pair-wise line matching based on collinearity constraint and matching redundancy[J]. Acta Geodaetica et Cartographica Sinica,2020,49(6):746-756.

    [14]

    LI Gang,ZENG Yawen,HUANG Huilan,et al. A multi-feature fusion slam system attaching semantic invariant to points and lines[J]. Sensors,2021,21(4). DOI: 10.3390/s21041196.

    [15]

    ZHENG Xianwei,YUAN Zhuang,DONG Zhen,et al. Smoothly varying projective transformation for line segment matching[J]. ISPRS Journal of Photogrammetry and Remote Sensing,2022,183:129-146. DOI: 10.1016/j.isprsjprs.2021.10.017

    [16]

    WANG Qiang,ZHANG Wei,LIU Xiaolong,et al. Line matching of wide baseline images in an affine projection space[J]. International Journal of Remote Sensing,2020,41(2):632-654. DOI: 10.1080/01431161.2019.1646937

    [17]

    SHEN Liang,ZHU Jiahua,XIN Qin,et al. Robust line segment mismatch removal using point-pair representation and Gaussian-uniform mixture formulation[J]. ISPRS Journal of Photogrammetry and Remote Sensing,2023,203:314-327. DOI: 10.1016/j.isprsjprs.2023.08.003

    [18] 刘肃艳,王竞雪,沈昭宇,等. 结合线对几何特征及单线描述符约束的直线匹配算法[J]. 武汉大学学报(信息科学版),2023,48(6):936-949.

    LIU Suyan,WANG Jingxue,SHEN Zhaoyu,et al. Line matching algorithm based on pair-wise geometric features and individual line descriptor constraints[J]. Geomatics and Information Science of Wuhan University,2023,48(6):936-949.

    [19] 张珊,张卡,赵立科,等. 结合网状描述符和单应约束的近景影像直线匹配[J]. 地球信息科学学报,2022,24(11):2186-2197.

    ZHANG Shan,ZHANG Ka,ZHAO Like,et al. Close-range image line matching based on mesh descriptor and homography constraint[J]. Journal of Geo-information Science,2022,24(11):2186-2197.

    [20]

    VAKHITOV A,LEMPITSKY A. Learnable line segment descriptor for visual SLAM[J]. IEEE Access,2019,7:39923-39934. DOI: 10.1109/ACCESS.2019.2901584

    [21]

    LANGE M,RAISCH C,SCHILLING A. WLD:a wavelet and learning based line descriptor for line feature matching[M]//KRÜGER J,NIESSNER M,STÜCKLER J. Vision,modeling,and visualization. Eindhoven:The Eurographics Association,2020:39-46.

    [22]

    PAUTRAT R,LIN J T,LARSSON V,et al. SOLD2:self-supervised occlusion-aware line description and detection[C]. IEEE/CVF Conference on Computer Vision and Pattern Recognition,Nashville,2021:11368-11378.

    [23]

    YOON S,KIM A. Line as a visual sentence:context-aware line descriptor for visual localization[J]. IEEE Robotics and Automation Letters,2021,6(4):8726-8733. DOI: 10.1109/LRA.2021.3111760

    [24]

    GUO Chunle,LI Chongyi,GUO Jichang,et al. Zero-reference deep curve estimation for low-light image enhancement[C]. IEEE/CVF Conference on Computer Vision and Pattern Recognition,Seattle,2020:1777-1786.

    [25] 郭永坤,朱彦陈,刘莉萍,等. 空频域图像增强方法研究综述[J]. 计算机工程与应用,2022,58(11):23-32.

    GUO Yongkun,ZHU Yanchen,LIU Liping,et al. Research review of space-frequency domain image enhancement methods[J]. Computer Engineering and Applications,2022,58(11):23-32.

    [26] 王智奇,李荣冰,刘建业,等. 基于同态滤波和直方图均衡化的图像增强算法[J]. 电子测量技术,2020,43(24):75-80.

    WANG Zhiqi,LI Rongbing,LIU Jianye,et al. Image enhancement algorithm based on homomorphic filtering and histogram equalization[J]. Electronic Measurement Technology,2020,43(24):75-80.

    [27]

    HANA F M,MAULIDA I D. Analysis of contrast limited adaptive histogram equalization (CLAHE) parameters on finger knuckle print identification[J]. Journal of Physics:Conference Series,2021,1764. DOI: 10.1088/1742-6596/1764/1/012049.

    [28]

    JIANG Yifan,GONG Xinyu,LIU Ding,et al. EnlightenGAN:deep light enhancement without paired supervision[J]. IEEE Transactions on Image Processing,2021,30:2340-2349. DOI: 10.1109/TIP.2021.3051462

    [29] 高宇彤. 基于透视变换与LBD描述子约束的特征线匹配算法[D]. 阜新:辽宁工程技术大学,2022.

    GAO Yutong. Feature line matching algorithm based on perspective transformation and LBD descriptor constraint[D]. Fuxin:Liaoning Technical University,2022.

  • 期刊类型引用(3)

    1. 雷富坤. 基于面保护原理的智能分布式煤矿供电防越级跳闸系统的研究. 电工技术. 2025(02): 109-113+119 . 百度学术
    2. 鹿伟强,于铄航,程卫健,穆艳祥. 实用的过零点电网频率计算方法. 煤矿机械. 2024(12): 4-6 . 百度学术
    3. 王世侠. 煤矿变电站自动化控制系统. 工矿自动化. 2023(S2): 150-151 . 本站查看

    其他类型引用(8)

图(15)  /  表(5)
计量
  • 文章访问数:  127
  • HTML全文浏览量:  36
  • PDF下载量:  19
  • 被引次数: 11
出版历程
  • 收稿日期:  2023-09-13
  • 修回日期:  2024-02-20
  • 网络出版日期:  2024-03-03
  • 刊出日期:  2024-02-24

目录

/

返回文章
返回