CN115100153B - 基于双目匹配的管内检测方法、装置、电子设备及介质 - Google Patents
基于双目匹配的管内检测方法、装置、电子设备及介质 Download PDFInfo
- Publication number
- CN115100153B CN115100153B CN202210752142.8A CN202210752142A CN115100153B CN 115100153 B CN115100153 B CN 115100153B CN 202210752142 A CN202210752142 A CN 202210752142A CN 115100153 B CN115100153 B CN 115100153B
- Authority
- CN
- China
- Prior art keywords
- image
- tube
- determining
- cost
- binocular
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 51
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 54
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 21
- 238000001914 filtration Methods 0.000 claims description 24
- 230000002776 aggregation Effects 0.000 claims description 19
- 238000004220 aggregation Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 17
- 230000009466 transformation Effects 0.000 claims description 8
- 230000002146 bilateral effect Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 238000003708 edge detection Methods 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 230000001276 controlling effect Effects 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 11
- 230000007547 defect Effects 0.000 description 5
- 230000007797 corrosion Effects 0.000 description 3
- 238000005260 corrosion Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001502 supplementing effect Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- HPTJABJPZMULFH-UHFFFAOYSA-N 12-[(Cyclohexylcarbamoyl)amino]dodecanoic acid Chemical compound OC(=O)CCCCCCCCCCCNC(=O)NC1CCCCC1 HPTJABJPZMULFH-UHFFFAOYSA-N 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/50—Depth or shape recovery
- G06T7/55—Depth or shape recovery from multiple images
- G06T7/593—Depth or shape recovery from multiple images from stereo images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
- G06T7/85—Stereo camera calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/97—Determining parameters from multiple pictures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20028—Bilateral filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20228—Disparity calculation for image-based rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Quality & Reliability (AREA)
- Image Processing (AREA)
Abstract
本申请公开了一种基于双目匹配的管内检测方法、装置、电子设备及介质,该方法包括:首先,根据两个相机的位置,确定双目相机的内部参数和外部参数;然后,利用双目相机拍摄管内的图像,得到双目图像对,并且将不同的拍摄相机拍摄的图像进行分组,确定左图像集和右图像集;接下来,利用图像代价算法,分别对左图像集和右图像集进行处理,再将处理结果综合,从而确定图像对匹配视差图;最后,利用深度计算公式,对图像对匹配视差图进行数据处理,确定管内图像深度图。通过双目相机采集管内图像,能够更好地定位管内图像,并且利用图像代价算法确定图像对匹配视差图,能够得到精度较高的代价值,从而实现提高管内图像的定位精度。
Description
技术领域
本发明涉及计算机视觉领域,尤其涉及一种基于双目匹配的管内检测方法、装置、电子设备及介质。
背景技术
随着城市的规模不断扩大,工业、排污、生活等各种管道铺设在城市管网中,与城市生活息息相关,管道检测维护的任务也越来越重要,成为一项不可忽视的重要工作。
目前,管道检测根据检测设备所处位置不同可分为管道外检测和管道内检测两大类。外检测是将检测设备放在管道外部检测管道的情况,内检测是指将检测器放在管道内进行检测。内检测通常采用双目立体视觉的方法确定管道缺陷位置,并且双目立体视觉中的半全局算法采用SGBM算法。然而,SGBM算法中的SGBM匹配方法采用Sobel X算子进行边缘检测,存在一定的误匹配率;SGBM算法中的代价计算主要为灰度绝对差值(SAD),匹配精度较差,且容易受到环境光照等因素的影响,难以获取到管道内部较为精准的深度信息。
因此,在管内检测的过程中,现有技术存在管内深度的检测精度较低的问题。
发明内容
有鉴于此,有必要提供一种基于双目匹配的管内检测方法、装置、电子设备及介质,用以解决现有技术中管内深度的检测精度较低的问题。
为了解决上述问题,本发明提供一种基于双目匹配的管内检测方法,包括:
获取双目相机的内部参数、外部参数以及双目相机拍摄的双目图像对,并基于内部参数和外部参数,确定左图像集和右图像集;
根据左图像集和右图像集,利用图像代价算法,确定图像对匹配视差图;
根据图像对匹配视差图,利用深度计算公式,确定管内图像深度图。
进一步地,获取双目相机的内部参数和外部参数,包括:
根据张正友棋盘格标定法计算,确定双目相机的内部参数和外部参数。
进一步地,基于内部参数和外部参数,确定左图像集和右图像集,包括:
获取双目相机拍摄的双目图像对,根据内部参数和外部参数,校正双目图像对,得到水平一致的图像对;
对水平一致的图像对进行灰度化处理,得到灰度图像对;
对灰度图像对进行分类,确定左图像集和右图像集。
进一步地,图像代价算法包括梯度代价算法和汉明距离代价算法,利用图像代价算法,确定图像对匹配视差图,包括:
根据左图像集和右图像集,基于梯度代价算法,确定管内图像的梯度代价;
基于汉明距离代价算法,确定管内图像的汉明距离代价;
基于梯度代价和汉明距离代价,根据路径聚合函数,确定图像对匹配视差图。
进一步地,基于梯度代价算法,确定管内图像的梯度代价,包括:
利用双边滤波分别对左图像集和右图像集进行处理,去除异常值,得到目标左图像集和目标右图像集;
根据Canny算法分别对目标左图像集和目标右图像集进行边缘检测,确定管内图像的梯度特征图;
对梯度特征图进行Census变换,确定管内图像的梯度代价。
进一步地,基于汉明距离代价算法,确定管内图像的汉明距离代价,包括:
分别获取左图像集和右图像集的代价计算区域;
对代价计算区域进行滤波,得到中值滤波,并将代价计算区域的中心值替换为中值滤波;
对左图像集和右图像集进行Census变换,确定管内图像的汉明距离代价。
进一步地,根据路径聚合函数,确定图像对匹配视差图,包括:
根据具有调整参数的动态规划法,调控梯度代价和汉明距离代价对应的惩罚参数,确定最优路径聚合函数;
根据最优路径聚合函数,确定图像对匹配视差图。
为了解决上述问题,本发明还提供一种基于双目匹配的管内检测装置,包括:
图像集获取模块,用于获取双目相机的内部参数、外部参数以及双目相机拍摄的双目图像对,并基于内部参数和外部参数,确定左图像集和右图像集;
图像对匹配视差图获取模块,用于根据左图像集和右图像集,利用图像代价算法,确定图像对匹配视差图;
管内深度图确定模块,用于根据图像对匹配视差图,利用深度计算公式,确定管内深度图。
为了解决上述问题,本发明还提供一种电子设备,包括处理器以及存储器,存储器上存储有计算机程序,计算机程序被处理器执行时,实现如前文所述的基于双目匹配的管内检测方法。
为了解决上述问题,本发明还提供一种计算机可读存储介质,存储介质存储有计算机程序指令,当计算机程序指令被计算机执行时,使计算机执行如前文所述的基于双目匹配的管内检测方法。
采用上述技术方案的有益效果是:本发明提供一种基于双目匹配的管内检测方法、装置、电子设备及介质,该方法包括:首先,通过双目相机采集管内图像,利用两个相机同时采集图像,能够实现立体匹配图像数据;然后,利用图像代价算法,确定图像对匹配视差图,不仅降低了外界光照条件的影响,并且能够有效抑制噪声,提高了管内图像的有效度;最后,利用深度计算公式,对管内图像数据进行处理,平衡了代价聚合中噪音的影响。本申请中,通过图像代价算法提取双目图像对中的数据信息,能够有效提高图像数据的精确度;由于是通过双目相机获取图像,能够从多个角度观测到管内的信息,并且生成立体图像,从而降低了低维度信息带来的精度影响,提高了管内图像深度的检测精度;最后,基于深度计算,有针对性地获取图像的深度图,提高了管内深度的检测精度。
附图说明
图1为本发明提供的基于双目匹配的管内检测方法一实施例的流程示意图;
图2为本发明提供的确定左图像集和右图像集一实施例的流程示意图;
图3为本发明提供的确定图像对匹配视差图一实施例的流程示意图;
图4为本发明提供的确定管内图像的梯度代价一实施例的流程示意图;
图5为本发明提供的确定管内图像的汉明距离代价一实施例的流程示意图;
图6为路径聚合示意图一实施例的角度图;
图7为管内图像深度图一实施例的结果图;
图8为本发明提供的基于双目匹配的管内检测装置的结构示意图;
图9为本发明提供的电子设备一实施例的结构框图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理,并非用于限定本发明的范围。
在陈述实施例之前,先对管内检测进行阐述:
管内检测是指利用管输介质驱动检测器在管道内运行,实时检测和记录管道的变形、腐蚀等损伤情况,并进行精准定位的工作。通过管道中检测器的检测结果,结合计算机等数据分析设备,对检测结果进行分析判断,从而了解管道的状况。通常可用于检测管道的变形、腐蚀和缺陷等。
进一步地,在通过管内检测的方法检测到管道的变形、腐蚀和缺陷等问题后,为了便于对管道进行后续的维修和保护,还需要对管内出现问题的部分进行定位。一般情况下,通过管内深度检测的方法,确定问题部位与摄像机的距离。
目前,一般采用双目立体视觉的方法确定管道缺陷位置,也就是说用双目立体视觉方法确定管道缺陷的位置。然而,双目立体视觉方法存在匹配率较低,匹配精度较差的问题,从而导致现有技术中存在管内深度的检测精度较低的问题。
为了解决现有技术中,管内深度的检测精度较低的问题,本发明提供了一种基于双目匹配的管内检测方法、装置、电子设备及介质,以下分别进行详细说明。
如图1所示,图1为本发明提供的基于双目匹配的管内检测方法一实施例的流程示意图,包括:
步骤S101:获取双目相机的内部参数、外部参数以及双目相机拍摄的双目图像对,并基于内部参数和外部参数,确定左图像集和右图像集。
步骤S102:根据左图像集和右图像集,利用图像代价算法,确定图像对匹配视差图。
步骤S103:根据图像对匹配视差图,利用深度计算公式,确定管内图像深度图。
本实施例中,首先,双目相机包括两个属性完全相同的相机,根据两个相机的位置,确定双目相机的内部参数和外部参数;然后,利用双目相机拍摄管内的图像,得到双目图像对,并且将不同的拍摄相机拍摄的图像进行分组,确定左图像集和右图像集;接下来,利用图像代价算法,分别对左图像集和右图像集进行处理,得到其对应的梯度代价和汉明距离代价,再将梯度代价和汉明距离代价综合,从而确定图像对匹配视差图;最后,利用深度计算公式,对图像对匹配视差图进行数据处理,确定管内图像深度图。
可以理解的是,上述实施例中通过双目相机采集管内图像,能够更好地定位管内图像,并且利用图像代价算法确定图像对匹配视差图,能够得到精度较高的代价值,从而实现提高管内图像的定位精度。
作为优选的实施例,在步骤S101中,为了确定双目相机的内部参数和外部参数,采用张正友棋盘格标定法进行计算。
在一具体实施例中,使用棋盘格为9*9、方格边长为25mm的黑白棋盘格作为标定板,在双目相机工作距离的范围内拍摄25组不同位姿的标定板双目图像,利用MATLABCamera Calibrator对拍摄双目图像实现标定。
在一具体实施例中,双目相机拍摄的图像分辨率为800*800。
在一具体实施例中,双目相机的内部参数包括:x方向上一个像素分别占多少长度单位,记作dx;y方向一个像素分别占多少长度单位,记作dy;图像的中心像素坐标和图像原点像素坐标之间相差的横向像素数,记作u0;图像的中心像素坐标和图像原点像素坐标之间相差的纵向像素数,记作v0;以及相机焦距,记作f。
其中,dx和dy是实现图像物理坐标系与像素坐标系转换的关键。
双目相机的外部参数包括:三个轴的旋转参数,记作ω、δ、θ;三个轴的平移参数,记作Tx、Ty、Tz。
其中,将每个轴的3*3旋转矩阵进行组合,也就是矩阵之间相乘,得到集合三个轴旋转信息,记作R,其大小还是3*3;再将R与三个轴的平移参数进行组合,得到一个3*4的矩阵,该矩阵是转换到标定纸坐标的关键。
在其他实施例中,还可以选择其他规格的标定板,还可以拍摄其他数量的标定板双目图像;另外,双目相机拍摄的图像分辨率还可以根据实际情况需要进行调整。
作为优选的实施例,在步骤S101中,为了根据双目相机拍摄的双目图像对,确定左图像集和右图像集,如图2所示,图2为本发明提供的确定左图像集和右图像集一实施例的流程示意图,包括:
步骤S111:获取双目相机拍摄的双目图像对,根据内部参数和外部参数,校正双目图像对,得到水平一致的图像对。
步骤S112:对水平一致的图像对进行灰度化处理,得到灰度图像对。
步骤S113:对灰度图像对进行分类,确定左图像集和右图像集。
本实施例中,为了获取双目图像对,利用标定好的双目相机对管内情况进行拍摄,由于双目相机包括两个相机,因此,能够成对获取管内图像,即双目图像对。
在步骤S111中,根据双目相机的内部参数和外部参数,对双目图像对进行校正,一方面,消除由于相机的影响导致的双目图像的畸变;另一方面,进行极线校正,使得每一对双目图像对中的图像内容都能够在水平位置保持一致,即得到水平一致的图像对。
在步骤S112中,将水平一致的图像对进行灰度化处理,使RGB空间中的图像数据转换到灰度空间,从而得到灰度图像对。
在步骤S113中,将灰度图像对按照拍摄相机的不同进行分类,确定左图像集和右图像集。
本实施例中,通过对由双目相机拍摄的双目图像对进行校正处理,消除了由于双目相机拍摄造成的外部影响;然后,通过灰度化处理双目图像对,突出了双目图像对的特征,并且根据灰度图像对能够更加有针对性地获取目标数据;最后,对双目图像对进行分组,有利于后续有序地对双目图像的数据进行处理。
作为优选的实施例,在步骤S102中,得到了左图像集和右图像集之后,还需要通过图像代价算法,才能够确定图像对匹配视差图,如图3所示,图3为本发明提供的确定图像对匹配视差图一实施例的流程示意图,包括:
步骤S121:根据左图像集和右图像集,基于梯度代价算法,确定管内图像的梯度代价。
步骤S122:基于汉明距离代价算法,确定管内图像的汉明距离代价。
步骤S123:基于梯度代价和汉明距离代价,根据路径聚合函数,确定图像对匹配视差图。
本实施例中,首先,通过梯度代价算法,将左图像集和右图像集的数据合并,确定管内图像的梯度代价;然后,通过汉明距离代价算法,确定管内图像的汉明距离代价;最后,通过路径聚合函数,将梯度代价和汉明距离代价融入,确定图像对匹配视差图。将左图像集和右图像集的数据进行整合,最终确定管内图像的对匹配视差图,以便于后续处理。
作为优选的实施例,如图4所示,图4为本发明提供的确定管内图像的梯度代价一实施例的流程示意图,包括:
步骤S1211:利用双边滤波分别对左图像集和右图像集进行处理,去除异常值,得到目标左图像集和目标右图像集。
步骤S1212:根据Canny算法分别对目标左图像集和目标右图像集进行边缘检测,确定管内图像的梯度特征图。
步骤S1213:对梯度特征图进行Census变换,确定管内图像的梯度代价。
本实施例中,首先,根据实际需要设置图像数据的阈值区间,双边滤波根据阈值区间,分别去除左图像集和右图像集中的异常值,得到目标左图像集和目标右图像集;然后,利用Canny算法分别对目标左图像集和目标右图像集进行边缘检测,使得目标左图像集和目标右图像集变得平滑,然后将平滑后的目标左图像集和目标右图像集组合到一起,提取出梯度特征,从而确定管内图像的梯度特征图;最后,根据Census变换,将管内图像的梯度特征图转化为梯度代价,从而确定管内图像的梯度代价。
本实施例中,一方面,利用双边滤波去除异常值,实现了从多方面提取需要的图像数据,降低了后续的计算难度;另一方面,利用Canny算法对图像进行平滑处理并确定梯度特征图,优化了图像数据的梯度特征,并且有效提取到了管内图像的梯度代价,实现了对管内图像的量化计算。
在一具体实施例中,双边滤波器像素输出值的计算公式为:
其中,P(i,j)为像素输出值;f(k,l)为像素输入值;ω(i,j,k,l)为加权系数。
作为优选的实施例,如图5所示,图5为本发明提供的确定管内图像的汉明距离代价一实施例的流程示意图,包括:
步骤S1221:分别获取左图像集和右图像集的代价计算区域。
步骤S1222:对代价计算区域进行滤波,得到中值滤波,并将代价计算区域的中心值替换为中值滤波。
步骤S1223:对左图像集和右图像集进行Census变换,确定管内图像的汉明距离代价。
本实施例中,首先,分别对左图像集和右图像集进行区域划分,确定其代价计算区域;然后,依次对代价计算区域进行滤波计算,得到每个代价计算区域的中值滤波,进一步地,将代价计算区域的中心值替换为中值滤波;最后,利用Census变换,确定管内图像的汉明距离代价。
本实施例中,通过将代价计算区域的中心值替换为中值滤波,有效抑制了噪声,提高了管内图像的处理精度。
在一具体实施例中,选择5×5窗口的正方形代价计算区域作为滤波窗口,对滤波窗口内像素进行中值滤波,得到中心像素值;然后,用中心像素值填充区域中心像素;最后,利用Census变换得到管内图像的变换码,从而确定管内图像的汉明距离代价。
本实施例中,通过滤波窗口的滤波,可以有效剔除中心像素点的噪声点,降低噪声干扰。
在一具体实施例中,在中心值附近遭遇噪声时,利用传统变换,将得到截然不同的变换码,而加入了滤波后,与未受到噪声干扰的结果仅相差1,实现了有效抑制噪声。
进一步地,在步骤S103中,获取到管内图像的梯度代价和汉明距离代价后,还需要根据路径聚合函数,才能确定图像对匹配视差图。
首先,根据具有调整参数的动态规划法,对梯度代价和汉明距离代价对应的惩罚参数进行调控,确定最优路径聚合函数;然后,根据最优路径聚合函数,确定图像对匹配视差图。
具体地,利用动态规划方法对每一个路径进行平滑处理时,都涉及到两个惩罚参数,这两个惩罚参数对于平滑项起着重要作用,直接影响到每个路径的聚合结果,并影响到最终的视差图。
本实施例中,为了避免由于两个惩罚参数的选择失误,导致最终结果的精度低,在两个惩罚参数的基础上添加调整参数,通过调整惩罚函数,提高惩罚参数的准确性,实现路径更加平滑,边界更加清晰。
在一具体实施例中,调整参数PN的计算公式为:
PN=(1cosθ)*t
其中,θ为p-r和p+r向量角,t为常数项,设定PN的阈值上限,通常取1。
如图6所示,图6为路径聚合示意图一实施例的角度图。
在一具体实施例中,确定调整参数后,路径聚合函数的计算公式为:
Lr(p,d)=C(p,d)+min(Lr(p-r,d),
Lr(p-r,d-1)+P1+PN,
Lr(p-r,d+1)+P1+PN,
其中,P1、P2为常数,并且规定P1的常数项要小于P2的常数项。
特别地,当视差变化较大时,会采用P2进行惩罚,当P2设置过大时,会导致边界模糊,PN参数会调整惩罚函数,达到减小P2的效果,从而平滑路径,清晰边界。
在一具体实施例中,基于WTA算法对匹配代价进行处理,处理过程中以标准欧式距离替换原算法欧氏距离,消除图像数据梯度不平衡带来的影响,最后输出图像对匹配视差图。
作为优选的实施例,在步骤S104中,根据图像对匹配视差图,确定管内图像深度图,其中,深度计算公式为:
d=Xl-Xr
其中,Z代表图像深度,d代表左右视差,f可以通过标定获取,T为双目相机光心距离。
在一具体实施例中,如图7所示,图7为管内图像深度图一实施例的结果图。
通过上述方法,根据梯度代价计算和汉明距离代价计算,确定双目图像对的图像对匹配视差图,然后根据图像对匹配视差图,确定管内图像深度图。由于梯度代价计算和汉明距离代价计算能够较好地平滑管内图像的噪声,从而实现提高管内图像的处理精度,从而得到精度较高的管内图像深度图。
为了检验本技术方案的效果,进行了对比实验。
对比试验中,在PC机上,对Middlebury数据集提供的标准图像对,进行本发明所述的双目匹配算法与SGBM匹配算法,测试环境为Ubantu20.04,C++配合CUDA与OpenCV编程实现,硬件为Intel corei9-10980XE 3.0GHz,内存空间为128G,处理结果如一下两表所示:
%
表1算法误匹配率
s
表2算法匹配时间
由上述对比试验的结果可知,与SGBM方法相比,本发明的误匹配率更低,匹配时间更短。也就是说,本技术方案不仅能够提高管内深度的检测精度,还能提高检测速度。
为了解决上述问题,本发明还提供了一种基于双目匹配的管内检测装置,如图8所示,图8为本发明提供的基于双目匹配的管内检测装置的结构示意图,基于双目匹配的管内检测装置800包括:
图像集获取模块801,用于获取双目相机的内部参数、外部参数以及双目相机拍摄的双目图像对,并基于内部参数和外部参数,确定左图像集和右图像集;
图像对匹配视差图获取模块802,用于根据左图像集和右图像集,利用图像代价算法,确定图像对匹配视差图;
管内深度图确定模块803,用于根据图像对匹配视差图,利用深度计算公式,确定管内深度图。
在一具体实施例中,由于管内存在光线不足的情况,因此,在对管内情况进行拍摄时,还需要根据实际情况进行补光操作,也就是说图像集获取模块801还包括补光模块。
本发明还相应提供了一种电子设备,如图9所示,图9为本发明提供的电子设备一实施例的结构框图。电子设备900可以是移动终端、桌上型计算机、笔记本、掌上电脑及服务器等计算设备。电子设备900包括处理器901以及存储器902,其中,存储器902上存储有基于双目匹配的管内检测程序903。
存储器902在一些实施例中可以是计算机设备的内部存储单元,例如计算机设备的硬盘或内存。存储器902在另一些实施例中也可以是计算机设备的外部存储设备,例如计算机设备上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(SecureDigital,SD)卡,闪存卡(Flash Card)等。进一步地,存储器902还可以既包括计算机设备的内部存储单元也包括外部存储设备。存储器902用于存储安装于计算机设备的应用软件及各类数据,例如安装计算机设备的程序代码等。存储器902还可以用于暂时地存储已经输出或者将要输出的数据。在一实施例中,基于双目匹配的管内检测程序903可被处理器901所执行,从而实现本发明各实施例的基于双目匹配的管内检测方法。
处理器901在一些实施例中可以是一中央处理器(Central Processing Unit,CPU),微处理器或其他数据处理芯片,用于运行存储器902中存储的程序代码或处理数据,例如执行基于双目匹配的管内检测程序等。
本实施例还提供了一种计算机可读存储介质,其上存储有基于双目匹配的管内检测程序,计算机该程序被处理器执行时,实现如上述任一技术方案所述的基于双目匹配的管内检测方法。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (9)
1.一种基于双目匹配的管内检测方法,其特征在于,包括:
获取双目相机的内部参数、外部参数以及所述双目相机拍摄的双目图像对,并基于所述内部参数和所述外部参数,确定左图像集和右图像集;
根据所述左图像集和所述右图像集,利用图像代价算法,确定图像对匹配视差图;
根据所述图像对匹配视差图,利用深度计算公式,确定管内图像深度图;
所述图像代价算法包括梯度代价算法和汉明距离代价算法,所述利用图像代价算法,确定图像对匹配视差图,包括:
根据所述左图像集和所述右图像集,基于所述梯度代价算法,确定所述管内图像的梯度代价;
基于所述汉明距离代价算法,确定所述管内图像的汉明距离代价;
基于所述梯度代价和所述汉明距离代价,根据路径聚合函数,确定图像对匹配视差图;
深度计算公式为:
其中,Z代表图像深度,d代表左右视差,f通过标定获取,T为双目相机光心距离。
2.根据权利要求1所述的基于双目匹配的管内检测方法,其特征在于,获取双目相机的内部参数和外部参数,包括:
根据张正友棋盘格标定法计算,确定所述双目相机的内部参数和外部参数。
3.根据权利要求1所述的基于双目匹配的管内检测方法,其特征在于,所述基于所述内部参数和所述外部参数,确定左图像集和右图像集,包括:
获取所述双目相机拍摄的双目图像对,根据所述内部参数和所述外部参数,校正所述双目图像对,得到水平一致的图像对;
对所述水平一致的图像对进行灰度化处理,得到灰度图像对;
对所述灰度图像对进行分类,确定左图像集和右图像集。
4.根据权利要求1所述的基于双目匹配的管内检测方法,其特征在于,所述基于所述梯度代价算法,确定管内图像的梯度代价,包括:
利用双边滤波分别对所述左图像集和所述右图像集进行处理,去除异常值,得到目标左图像集和目标右图像集;
根据Canny算法分别对所述目标左图像集和所述目标右图像集进行边缘检测,确定所述管内图像的梯度特征图;
对所述梯度特征图进行Census变换,确定所述管内图像的梯度代价。
5.根据权利要求4所述的基于双目匹配的管内检测方法,其特征在于,所述基于所述汉明距离代价算法,确定所述管内图像的汉明距离代价,包括:
分别获取所述左图像集和所述右图像集的代价计算区域;
对所述代价计算区域进行滤波,得到中值滤波,并将所述代价计算区域的中心值替换为所述中值滤波;
对所述左图像集和所述右图像集进行Census变换,确定所述管内图像的汉明距离代价。
6.根据权利要求5所述的基于双目匹配的管内检测方法,其特征在于,所述根据路径聚合函数,确定图像对匹配视差图,包括:
根据具有调整参数的动态规划法,调控所述梯度代价和所述汉明距离代价对应的惩罚参数,确定最优路径聚合函数;
根据所述最优路径聚合函数,确定图像对匹配视差图。
7.一种基于双目匹配的管内检测装置,其特征在于,包括:
图像集获取模块,用于获取双目相机的内部参数、外部参数以及所述双目相机拍摄的双目图像对,并基于所述内部参数和所述外部参数,确定左图像集和右图像集;
图像对匹配视差图获取模块,用于根据所述左图像集和所述右图像集,利用图像代价算法,确定图像对匹配视差图;
管内深度图确定模块,用于根据所述图像对匹配视差图,利用深度计算公式,确定管内深度图;
所述图像代价算法包括梯度代价算法和汉明距离代价算法,所述利用图像代价算法,确定图像对匹配视差图,包括:
根据所述左图像集和所述右图像集,基于所述梯度代价算法,确定所述管内图像的梯度代价;
基于所述汉明距离代价算法,确定所述管内图像的汉明距离代价;
基于所述梯度代价和所述汉明距离代价,根据路径聚合函数,确定图像对匹配视差图;
深度计算公式为:
其中,Z代表图像深度,d代表左右视差,f通过标定获取,T为双目相机光心距离。
8.一种电子设备,其特征在于,包括处理器以及存储器,所述存储器上存储有计算机程序,所述计算机程序被所述处理器执行时,实现如权利要求1-6任一项所述的基于双目匹配的管内检测方法。
9.一种计算机可读存储介质,其特征在于,所述存储介质存储有计算机程序指令,当所述计算机程序指令被计算机执行时,使计算机执行如权利要求1至6中任一所述的基于双目匹配的管内检测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210752142.8A CN115100153B (zh) | 2022-06-29 | 2022-06-29 | 基于双目匹配的管内检测方法、装置、电子设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210752142.8A CN115100153B (zh) | 2022-06-29 | 2022-06-29 | 基于双目匹配的管内检测方法、装置、电子设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115100153A CN115100153A (zh) | 2022-09-23 |
CN115100153B true CN115100153B (zh) | 2024-09-24 |
Family
ID=83295842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210752142.8A Active CN115100153B (zh) | 2022-06-29 | 2022-06-29 | 基于双目匹配的管内检测方法、装置、电子设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115100153B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117593349A (zh) * | 2023-11-24 | 2024-02-23 | 中山大学·深圳 | 图像变形监测方法、装置、设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114519713A (zh) * | 2022-02-24 | 2022-05-20 | 银牛微电子(无锡)有限责任公司 | 深度测量方法、双目立体视觉系统、装置及存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020037575A1 (zh) * | 2018-08-22 | 2020-02-27 | 深圳市大疆创新科技有限公司 | 图像深度估计方法、装置、可读存储介质及电子设备 |
CN110675449B (zh) * | 2019-09-02 | 2020-12-08 | 山东科技大学 | 一种基于双目相机的离岸流检测方法 |
CN110830781B (zh) * | 2019-10-30 | 2021-03-23 | 歌尔科技有限公司 | 一种基于双目视觉的投影图像自动校正方法及系统 |
CN112070821B (zh) * | 2020-07-31 | 2023-07-25 | 南方科技大学 | 一种低功耗立体匹配系统及获取深度信息的方法 |
-
2022
- 2022-06-29 CN CN202210752142.8A patent/CN115100153B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114519713A (zh) * | 2022-02-24 | 2022-05-20 | 银牛微电子(无锡)有限责任公司 | 深度测量方法、双目立体视觉系统、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN115100153A (zh) | 2022-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107633536B (zh) | 一种基于二维平面模板的相机标定方法及系统 | |
CN108230397B (zh) | 多目摄像机标定与校正方法和装置、设备、程序和介质 | |
WO2019233264A1 (zh) | 图像处理方法、计算机可读存储介质和电子设备 | |
CN108335331B (zh) | 一种钢卷双目视觉定位方法及设备 | |
CN114029946A (zh) | 一种基于3d光栅引导机器人定位抓取的方法、装置及设备 | |
CN108381549B (zh) | 一种双目视觉引导机器人快速抓取方法、装置及存储介质 | |
WO2019232793A1 (zh) | 双摄像头标定方法、电子设备、计算机可读存储介质 | |
CN109479082B (zh) | 图象处理方法及装置 | |
WO2012011227A1 (ja) | ステレオ測距装置及びステレオ測距方法 | |
CN113326787B (zh) | 指针仪表读数的自动识别方法、系统及设备 | |
CN109559353B (zh) | 摄像模组标定方法、装置、电子设备及计算机可读存储介质 | |
CN115100153B (zh) | 基于双目匹配的管内检测方法、装置、电子设备及介质 | |
CN102422321A (zh) | 摄像装置以及图像复原方法 | |
Kruger et al. | In-factory calibration of multiocular camera systems | |
CN109584312A (zh) | 摄像头标定方法、装置、电子设备和计算机可读存储介质 | |
CN105335959B (zh) | 成像装置快速对焦方法及其设备 | |
CN116012242A (zh) | 一种相机畸变矫正效果评价方法、装置、介质及设备 | |
CN114049304A (zh) | 一种3d光栅检测方法、装置、计算机设备及可读存储介质 | |
CN113298885A (zh) | 双目标定方法及装置、设备和存储介质 | |
CN114332092B (zh) | 缺陷图像检测方法、系统、智能终端和介质 | |
CN113592934B (zh) | 一种基于单目相机的目标深度与高度测量方法及装置 | |
CN106780425B (zh) | 一种vver蒸汽发生器传热管涡流检测系统定位方法 | |
CN115661258A (zh) | 标定方法及装置、畸变校正方法及装置、存储介质、终端 | |
JP6909090B2 (ja) | コンクリート表面上の損傷部の定量評価方法、定量評価装置、定量評価システム及びプログラム | |
EP2953096B1 (en) | Information processing device, information processing method, system and carrier means |
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 |