CN116542927A - 一种基于相机视觉的轨道板产品平面度检测方法 - Google Patents

一种基于相机视觉的轨道板产品平面度检测方法 Download PDF

Info

Publication number
CN116542927A
CN116542927A CN202310501340.1A CN202310501340A CN116542927A CN 116542927 A CN116542927 A CN 116542927A CN 202310501340 A CN202310501340 A CN 202310501340A CN 116542927 A CN116542927 A CN 116542927A
Authority
CN
China
Prior art keywords
image
camera
gray
track
flatness
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.)
Pending
Application number
CN202310501340.1A
Other languages
English (en)
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.)
Shenyang Institute of Technology
Original Assignee
Shenyang Institute of Technology
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 Shenyang Institute of Technology filed Critical Shenyang Institute of Technology
Priority to CN202310501340.1A priority Critical patent/CN116542927A/zh
Publication of CN116542927A publication Critical patent/CN116542927A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/194Segmentation; Edge detection involving foreground-background segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/593Depth or shape recovery from multiple images from stereo images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • G06T7/85Stereo camera calibration
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明一种基于相机视觉的轨道板产品平面度检测方法,涉及CRTSⅢ型轨道板生产制造领域。其特征在于:该检测方法按包括以下步骤:双目视觉***构建以及相机标定;轨道板图像采集;轨道板图像预处理;轨道板图像边缘检测;进行图像匹配获取视差图;平面度信息的提取与分析。通过相机与轨道板的相对运动采集图像,并通过中值滤波对图像进行预处理,增强图像效果。采用优化后的Canny算法对图像进行边缘检测,捕捉到图像更多的细节。利用像素灰度差的平方和函数(SSD)进行图像位置匹配,最后通过三角测量原理来提取平面度信息,得到深度图像,获得实验结果。本发明方法采用图像处理技术进行检测,在节约人力资源的同时,检测效率和准确率也大幅提高。

Description

一种基于相机视觉的轨道板产品平面度检测方法
技术领域
本发明涉及CRTSⅢ型轨道板生产制造领域,具体说是一种基于相机视觉的轨道板产品平面度检测方法。
背景技术
随着我国对交通运输领域的大规模建设,尤其是加快推进高速铁路建设,保障轨道的安全运输成为一项必需的任务。CRTSⅢ型轨道板是我国具有自主知识产权的板式无砟轨道,是带挡肩的新型单元板式无砟轨道结构。目前,CRTSⅢ轨道板作为轨道运输的关键,需要承受来自上方列车的压力,轨道板表面容易产生各种缺陷,这会导致十分严重的事故发生。所以如何对轨道板进行快速精确的平面度检测对于保障轨道安全运输越来越重要。
轨道板的平面度检测方法包含超声导波检测、游标卡尺法、数字摄影测量法和激光扫描仪法等。超声导波检测由于受到信号的限制可能会产生检测盲点,所以检测的准确度不高。游标卡尺法测量原理简单,但要使用多种检测工具且依靠人工记录数据,致使检测效率低。数字摄影测量法通过拍摄预先放置在特定位置的靶标,经过一定的数据处理得到结果。但需提前在轨道板上布设大量靶标,使得测量效率较低。目前使用普遍的轨道板平面度检测方法是全站仪法。通过全站仪结合特殊的工装,得到轨道板上设置的棱镜坐标,进而转换得到轨道板平面度信息。全站仪法对一块轨道板的平面度检测时间约15-30分钟,精度一般大于0.5mm,但是该测量方法检测时间过长,对测量场地环境的要求比较高,不适合用于大批量的轨道板平面度检测,不能较好的满足实际测量的需求。
随着无损检测技术的不断发展,结合计算机视觉技术,将其应用在轨道板平面度检测上,在节约人力资源的同时,检测效率和准确率也大幅提高。针对上诉检测方式存在的问题,本发明研制一种基于相机视觉的CRTSⅢ型轨道板产品平面度检测方法来解决上述问题。
发明内容
发明目的:针对现有的检测方法难以同时保证检测精度和检测效率的问题,本发明提供一种基于相机视觉的CRTSⅢ型轨道板产品平面度检测方法,实现对CRTSⅢ型轨道板产品的平面度检测。
技术方案:
本发明的技术方案是这样实现的:
一种基于相机视觉的轨道板产品平面度检测方法,其特征在于包括以下步骤:
S1:双目视觉***构建以及相机标定;
S2:轨道板图像采集;
S3:轨道板图像预处理;
S4:轨道板图像边缘检测;
S5:进行图像匹配获取视差图;
S6:平面度信息的提取与分析;
其中,所述轨道板主要是CRTSⅢ型轨道板。
步骤S1所述双目视觉***构建,具体为:使用两台相机构成的双目立体视觉***,两台相机同时拍摄同一个实际场景的二维图像,使用优化后的Canny算法进行图像边缘检测,以及像素灰度差的平方和函数SSD完成图像的匹配,进行三维坐标和深度信息,进而对场景进行轨道平面度信息的提取;所述相机标定,具体为:相机内参数的标定和相机外参数的标定;相机内参为dx、dy、u0、v0。dx和dy表示:x方向和y方向的一个像素分别占多少长度单位;u0和v0表示图像的中心像素坐标和图像原点像素坐标之间相差的横向和纵向像素数;相机的外部参数为三个轴的旋转参数ω、δ、θ,三个轴的平移参数Tx、Ty、Tz
步骤S2所述轨道板图像采集,具体为:使用带有两个传感器的双目视觉***,进行相机与轨道板之间的相对运动,进而采集高质量图像。
步骤S3所述图像预处理,具体为:采用中值滤波,将领域内像素灰度值相差较大的像素灰度值替换成和领域像素灰度值想接近的值,在消除噪声的同时,保证不丢失图像的细节。
步骤S4所述轨道板图像边缘检测,具体为:采用基于改进Canny算法进行轨道表面图像边缘检测,主要是采用改进的迭代算法对高低阈值进行计算,能够使检测图像拥有更多的边缘细节,具体流程如下:
1)利用直方图得到图像灰度的最大值Tmax和最小值Tmin,计算二者的均值
2)根据均值把原始图像分成高灰度区域和低灰度区域,灰度值≥K的区域为高灰度区域H,灰度值<K的区域为低灰度区域L;
3)求出H区域的灰度最大值T1和最小值T2,令根据K1把H区域分为前景和背景,求出二者的灰度值G1和G2
其中,E(i,j)表示图像中点(i,j)的灰度值;F(i,j)为该点的权重系数;
4)H区域内,灰度值≥K1的区域为高灰度区域H1,灰度值<K1的区域为高灰度区域H2
5)求出H1区域内的灰度最大值T3和最小值T4,令根据K2把H1区域分为前景和背景,求出二者的灰度值G3和G4
6)求出新的阈值
7)若T为预先设定的一个参数,那么K3是经过迭代法求出的高阈值,否则将K3的值赋给K1,再重复进行3-7,知道获得满足要求的高阈值;
8)同理,可在N区域内找到符合要求的低阈值。
步骤S5所述进行图像匹配获取视差图,具体为:
第一步:从两个相机中采集的图像选一幅作为参考图像,另外一副就是待匹配图像;
第二步:在参考图像中,结合轨道在三维场景中的实际结构,选取匹配特征;
第三步:在待匹配图像中,找到与参考图像中选区的实际物理结构的特征相同的对应图像特征;
第四步:采用像素灰度差的平方和函数SSD确定两个特征的相对位置。
步骤S6所述平面度信息的提取与分析,具体为:通过获取的视觉差,结合前面的相机的内外参数,使用三角测量原理来提取平面度信息,得到深度图像,获得实验结果。
第四步:采用基于灰度的误差平方和函数SSD确定两个特征的相对位置,计算两幅图像相似度,以左相机采集的图像的源匹配点为中心,定义一个窗口W,统计其窗口的灰度值的和,然后在右相机采集的图像中逐步计算其左右窗口的灰度差值的平方和;Il(u,v)和Ir(u+d,v)分别代表第一步中左右相机采集到的图像,公式如下:
C(u,v,d)=∑(i,j)∈W(Il(u+i,v+j)-Ir(u+i+d,v+j))2
有益效果:
本发明具有以下优点及有益效果:
1、本发明的检测方法采用优化后的Canny算法进行图像边缘检测有很好的自适应性,能更好地保留图像的细节。
2、使用带有两个图像传感器的双目视觉***,直接利用图像传感器对图像信息进行实时转换,提高了最终的检测精度和检测效率。
3、本发明的检测方法原理可靠,设备构造简单,具有良好的应用场景,可实现对CRTSⅢ型轨道板平面度的实时检测。
附图说明
图1为本发明方法的整体流程图;
图2为本发明方法的双目视觉***;
图3为本发明方法的轨道板图像;
图4为本发明方法的坐标标定转换图。
具体实施方式
下面结合附图对本发明方法作进一步的详细说明。
一种基于相机视觉的CRTSⅢ型轨道板产品平面度检测方法,参照图1-4,主要包括以下步骤:
1:双目视觉***构建以及相机标定:
如图2所示,是本方法的双目视觉***布置原理图:对于实际轨道板某位置点P对应于相机采集到的两幅图像中的位置点P1、P2,视觉***利用三角相似关系得到位置点P的三维坐标,获得该位置点的深度信息。要把两个相机的相对位置固定成这样。在此之前先要完成相机的标定,通过在棋盘标定板中选取标定点,标定选取的关键点在世界坐标系下和图像坐标系下的坐标,就可以建立世界坐标同左右相机坐标间正确的坐标转换关系,如图4所示。其中,图4所涉及的四个坐标系及其含义为:Ow-Xw-Yw-Zw是世界坐标系;Oc-Xc Yc-Zc是相机坐标系,光心为原点;O-x-y是图像坐标系,光心为原点;u-v是像素坐标系,原点为图像左上角;P代表世界坐标系中的一点,即为生活中真实的一点。p是P在图像中的成像点;f代表相机焦距。
2:轨道板图像采集:图像采集的重点在于如何确保相机的运动路线能够将目标轨道板的信息全面拍摄采集出来。轨道板随运转***到达双目***检测范围之内,并保持一个稳定的状态。然后通过PLC控制伺服电机,带动相机以45mm/s的速度按照规定路线匀速移动,实现相机与轨道板之间的相对运动采集高质量图像,进而传输到计算机中,轨道板中的承轨槽部分如图3所示;
3:轨道板图像预处理:采用中值滤波,将领域内像素灰度值相差较大的像素灰度值替换成和领域像素灰度值想接近的值,在消除噪声的同时,还尽量保证不丢失图像的细节;使用模板6*6大小的滤波模板,将模板中心与图中的一个像素重合,收集模板中所有像素的灰度值,将其进行大小排序,找到中间的灰度值,将它赋值给模板中心位置的像素,得到新的像素值,把图像的每一个像素点都进行这样的处理,就可以得到中值滤波处理之后的新的图像;
4:轨道板图像边缘检测(轨道表面提取):将中值滤波之后的新图像,采用基于改进Canny算法进行轨道表面图像边缘检测,主要是采用改进的迭代算法对高低阈值进行计算,能够使检测图像拥有更多的边缘细节;
5:进行图像匹配获取视差图:
第一步:从两个相机中采集的图像进行边缘检测之后,选左相机采集的图像作为参考图像,右相机采集图像为待匹配图像;
第二步:在参考图像中,结合轨道在三维场景中的实际结构,选取匹配特征点,例如选取一些轨道板特殊明显的角点、线等;
第三步:在待匹配图像中,找到与参考图像中选区的实际物理结构的特征相同的对应图像特征;
第四步:采用基于灰度的误差平方和函数(SSD)确定两个特征的相对位置,计算两幅图像相似度,以左相机采集的图像的源匹配点为中心,定义一个窗口W,统计其窗口的灰度值的和,然后在右相机采集的图像中逐步计算其左右窗口的灰度差值的平方和。Il(u,v)和Ir(u+d,v)分别代表第一步中左右相机采集到的图像,公式如下:
6:平面度信息的提取与分析:通过获取的视觉差,结合前面的相机的内外参数,使用三角测量原理来提取平面度信息,得到深度图像,获得实验结果。如图2所示,由三角测量原理,可得点P的深度信息为:Zc=Bf/d,其中B为两相机光心之间的距离,f为相机焦距,d为搜索匹配的视差:这样就可以得到轨道板图像的平面度信息。
为了评估这种方法对轨道板产品平面度检测的效果,按照轨道板标准平面度误差范围与其他方法进行比较,可以看到不管是检测精度还是检测时间上都比其他方法好,也证明了本发明方法的有效性。
如图1所示,首先要构建双目视觉***,也就是标定2台相机的相对位置关系。
如图2所示,点P在左右相机中立体校正前的位置是点P1、P2,c1、c2为左右相机的光心,2点之间的连线是基线,e1、e2为左右相机的成像平面与基线的交点,也就是极点。由极线的约束条件可知,左图像中点P1在右图像中的对应点必定位于直线e2P2上,极线约束令对应点的搜索由二维搜索降为一维搜索。但极线往往为一条斜线,为了进一步降低搜寻对应点的复杂度,还需要进行极线校正。校正后,两台相机光轴平行,极点位于无穷远处,像点在左右相机中高度相同,对应点搜索只需在同一行上进行。
表1为按照以上方法实施后,所得的对比数据;
表1检测方法对比数据

Claims (8)

1.一种基于相机视觉的轨道板产品平面度检测方法,其特征在于包括以下步骤:
S1:双目视觉***构建以及相机标定;
S2:轨道板图像采集;
S3:轨道板图像预处理;
S4:轨道板图像边缘检测;
S5:进行图像匹配获取视差图;
S6:平面度信息的提取与分析;
其中,所述轨道板主要是CRTSⅢ型轨道板。
2.根据权利要求1所述的一种基于相机视觉的轨道板产品平面度检测方法,其特征在于:步骤S1所述双目视觉***构建,具体为:使用两台相机构成的双目立体视觉***,两台相机同时拍摄同一个实际场景的二维图像,使用优化后的Canny算法进行图像边缘检测,以及像素灰度差的平方和函数SSD完成图像的匹配,进行三维坐标和深度信息,进而对场景进行轨道平面度信息的提取;所述相机标定,具体为:相机内参数的标定和相机外参数的标定;相机内参为dx、dy、u0、v0;dx和dy表示:x方向和y方向的一个像素分别占多少长度单位;u0和v0表示图像的中心像素坐标和图像原点像素坐标之间相差的横向和纵向像素数;相机的外部参数为三个轴的旋转参数ω、δ、θ,三个轴的平移参数Tx、Ty、Tz
3.根据权利要求1所述的一种基于相机视觉的轨道板产品平面度检测方法,其特征在于:步骤S2所述轨道板图像采集,具体为:使用带有两个传感器的双目视觉***,进行相机与轨道板之间的相对运动,进而采集高质量图像。
4.根据权利要求1所述的一种基于相机视觉的轨道板产品平面度检测方法,其特征在于:步骤S3所述图像预处理,具体为:采用中值滤波,将领域内像素灰度值相差较大的像素灰度值替换成和领域像素灰度值想接近的值,在消除噪声的同时,保证不丢失图像的细节。
5.根据权利要求1所述的一种基于相机视觉的轨道板产品平面度检测方法,其特征在于:步骤S4所述轨道板图像边缘检测,具体为:采用基于改进Canny算法进行轨道表面图像边缘检测,主要是采用改进的迭代算法对高低阈值进行计算,能够使检测图像拥有更多的边缘细节,具体流程如下:
1)利用直方图得到图像灰度的最大值Tmax和最小值Tmin,计算二者的均值
2)根据均值把原始图像分成高灰度区域和低灰度区域,灰度值≥K的区域为高灰度区域H,灰度值<K的区域为低灰度区域L;
3)求出H区域的灰度最大值T1和最小值T2,令根据K1把H区域分为前景和背景,求出二者的灰度值G1和G2
其中,E(i,j)表示图像中点(i,j)的灰度值;F(i,j)为该点的权重系数;
4)H区域内,灰度值≥K1的区域为高灰度区域H1,灰度值<K1的区域为高灰度区域H2
5)求出H1区域内的灰度最大值T3和最小值T4,令根据K2把H1区域分为前景和背景,求出二者的灰度值G3和G4
6)求出新的阈值
7)若T为预先设定的一个参数,那么K3是经过迭代法求出的高阈值,否则将K3的值赋给K1,再重复进行3-7,知道获得满足要求的高阈值;
8)同理,可在N区域内找到符合要求的低阈值。
6.根据权利要求1所述的一种基于相机视觉的轨道板产品平面度检测方法,其特征在于:步骤S5所述进行图像匹配获取视差图,具体为:
第一步:从两个相机中采集的图像选一幅作为参考图像,另外一副就是待匹配图像;
第二步:在参考图像中,结合轨道在三维场景中的实际结构,选取匹配特征;
第三步:在待匹配图像中,找到与参考图像中选区的实际物理结构的特征相同的对应图像特征;
第四步:采用像素灰度差的平方和函数SSD确定两个特征的相对位置。
7.根据权利要求1所述的一种基于相机视觉的轨道板产品平面度检测方法,其特征在于:步骤S6所述平面度信息的提取与分析,具体为:通过获取的视觉差,结合前面的相机的内外参数,使用三角测量原理来提取平面度信息,得到深度图像,获得实验结果。
8.根据权利要求6所述的一种基于相机视觉的轨道板产品平面度检测方法,其特征在于:第四步:采用基于灰度的误差平方和函数SSD确定两个特征的相对位置,计算两幅图像相似度,以左相机采集的图像的源匹配点为中心,定义一个窗口W,统计其窗口的灰度值的和,然后在右相机采集的图像中逐步计算其左右窗口的灰度差值的平方和;Il(u,v)和Ir(u+d,v)分别代表第一步中左右相机采集到的图像,公式如下:
C(u,v,d)=∑(i,j)∈W(Il(u+i,v+j)-Ir(u+i+d,v+j))2
CN202310501340.1A 2023-05-06 2023-05-06 一种基于相机视觉的轨道板产品平面度检测方法 Pending CN116542927A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310501340.1A CN116542927A (zh) 2023-05-06 2023-05-06 一种基于相机视觉的轨道板产品平面度检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310501340.1A CN116542927A (zh) 2023-05-06 2023-05-06 一种基于相机视觉的轨道板产品平面度检测方法

Publications (1)

Publication Number Publication Date
CN116542927A true CN116542927A (zh) 2023-08-04

Family

ID=87451753

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310501340.1A Pending CN116542927A (zh) 2023-05-06 2023-05-06 一种基于相机视觉的轨道板产品平面度检测方法

Country Status (1)

Country Link
CN (1) CN116542927A (zh)

Similar Documents

Publication Publication Date Title
KR100776649B1 (ko) 깊이 정보 기반 스테레오/다시점 영상 정합 장치 및 방법
JP3539788B2 (ja) 画像間対応付け方法
CN108986070B (zh) 一种基于高速视频测量的岩石裂缝扩展实验监测方法
CN104537707B (zh) 像方型立体视觉在线移动实时测量***
CN110334701B (zh) 数字孪生环境下基于深度学习和多目视觉的数据采集方法
CN106996748A (zh) 一种基于双目视觉的轮径测量方法
CN110288659B (zh) 一种基于双目视觉的深度成像及信息获取方法
CN111091076B (zh) 基于立体视觉的隧道限界数据测量方法
CN116342718B (zh) 一种线激光3d相机的标定方法、装置、存储介质及设备
CN111780678A (zh) 一种轨道板预埋套管直径的测量方法
CN112232319A (zh) 一种基于单目视觉定位的扫描拼接方法
CN114037762A (zh) 基于图像与高精度地图配准的实时高精度定位方法
CN115330684A (zh) 基于双目视觉与线结构光的水下构筑物表观缺陷检测方法
CN111583342A (zh) 一种基于双目视觉的目标快速定位方法及装置
CN111429571A (zh) 一种基于时空图像信息联合相关的快速立体匹配方法
CN104614372B (zh) 一种太阳能硅片检测方法
CN113393413B (zh) 基于单目与双目视觉协同的水域面积测量方法和***
CN110726402B (zh) 非正交轴系激光全站仪的激光点视觉引导方法
CN112525106A (zh) 基于三相机协同激光的3d探测方法及装置
CN116958218A (zh) 一种基于标定板角点对齐的点云与图像配准方法及设备
CN108090930A (zh) 基于双目立体相机的障碍物视觉检测***及方法
CN116542927A (zh) 一种基于相机视觉的轨道板产品平面度检测方法
CN114608522B (zh) 一种基于视觉的障碍物识别与测距方法
CN114199205B (zh) 基于改进四叉树orb算法的双目测距方法
Liu et al. A Method for Synchronous Automated Extrinsic Calibration of LiDAR and Cameras Based on a Circular Calibration Board

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