CN118224999A - 一种基于结构光的钢轨磨耗识别方法及设备 - Google Patents
一种基于结构光的钢轨磨耗识别方法及设备 Download PDFInfo
- Publication number
- CN118224999A CN118224999A CN202410650265.XA CN202410650265A CN118224999A CN 118224999 A CN118224999 A CN 118224999A CN 202410650265 A CN202410650265 A CN 202410650265A CN 118224999 A CN118224999 A CN 118224999A
- Authority
- CN
- China
- Prior art keywords
- image
- camera
- exposure time
- feature map
- rail
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005299 abrasion Methods 0.000 title claims abstract description 23
- 229910000831 Steel Inorganic materials 0.000 claims abstract description 41
- 239000010959 steel Substances 0.000 claims abstract description 41
- 238000001514 detection method Methods 0.000 claims abstract description 35
- 238000010586 diagram Methods 0.000 claims abstract description 24
- 238000005259 measurement Methods 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 description 11
- 230000007547 defect Effects 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 101100512352 Enterococcus faecalis (strain TX4000 / JH2-2) mapP gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61K—AUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
- B61K9/00—Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
- B61K9/08—Measuring installations for surveying permanent way
-
- 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
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/54—Extraction of image or video features relating to texture
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
- G06V10/758—Involving statistics of pixels or of feature values, e.g. histogram matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/761—Proximity, similarity or dissimilarity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/98—Detection or correction of errors, e.g. by rescanning the pattern or by human intervention; Evaluation of the quality of the acquired patterns
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Artificial Intelligence (AREA)
- Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Mechanical Engineering (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于结构光的钢轨磨耗识别方法及设备,属于钢轨磨耗检测技术领域。本发明通过检测小车获取钢轨的第一图像、第二图像、第三图像、第四图像。根据振动传感器的振动信号确认伸缩缝的位置,根据伸缩缝的位置从第一图像和第三图像中提取第一特征图。根据紧固件的纹理特征,从第二图像和第四图像中提取第二特征图。根据第一特征图调整第一摄像机和第三摄像机的第一曝光时间和第三曝光时间。根据第二特征图调整第二摄像机和第四摄像机的第二曝光时间和第四曝光时间。通过调整曝光时间,显著提高图像的清晰度。根据第一图像、第二图像、第三图像以及第四图像,生成钢轨的截面轮廓,根据截面轮廓获得更精确的磨耗值。
Description
技术领域
本发明涉及钢轨磨耗检测技术领域,尤其涉及一种基于结构光的钢轨磨耗识别方法及设备。
背景技术
由于列车长时间高负荷运转,对钢轨带来日益递增的磨耗,对列车正常运行带来极大的安全风险,准确而高效的检测钢轨磨耗极其重要。目前,由于结构光技术的快速发展,特别是在帧率上的提高,其非接触性、高精度、高效率动态检测的优点被广泛应用于轨道行业的检测。公开号为CN108955576A的中国专利申请公开了一种钢轨廓形动态检测中的多线结构光自标定方法及系统,所述方法从钢轨的场景图像中根据角提取算法获得多线结构光条上的特征点,根据特征点及钢轨自身的几何形状,并结合结构光视觉模型搭建钢轨三维廓形,并根据钢轨三维廓形获得投影廓形,对投影廓形和标准廓形在轨鄂、轨头内侧和轨底三个非磨损区域上的尺寸和形状偏差,利用曲线相似性度量方法进行量化,根据量化偏差迭代优化多线结构光平面参数,实现多线结构光的自标定。该方法解决了依赖标定模板不能自主迅速对结构光在线自标定的问题,提高测量精度,但是没有考虑到检测环境中环境光的强度和检测时的震动对检测结果的影响。由于轨道处于自然场景下,结构光获取的数据经常干扰特别多,导致数据的稳定性不好,所以现有技术有进一步改进的必要。
发明内容
为了解决上述现有技术存在的缺陷,本发明提出了一种基于结构光的钢轨磨耗识别方法及设备。本发明的基于结构光的钢轨磨耗识别方法通过调节第一曝光时间和第二曝光时间、第三曝光时间和第四曝光时间提高第一图像、第二图像、第三图像以及第四图像的清晰度,根据检测小车的振动情况确认伸缩缝的第一特征图,结合紧固件的纹理特征识别紧固件的第二特征图,可以显著提高图像质量。
本发明的技术方案是这样实现的:
一种基于结构光的钢轨磨耗识别方法,包括以下步骤:
步骤1:在检测小车上安装由第一激光器、第一摄像机和第二摄像机组成的第一测量单元和由第二激光器、第三摄像机和第四摄像机组成的第二测量单元以及振动传感器;
步骤2:钢轨经紧固件安装在基座上,相邻的钢轨具有伸缩缝,读取伸缩缝的第一基准特征图,读取紧固件的第二基准特征图,所述检测小车沿钢轨移动;
步骤3:第一激光器向钢轨的一个侧面发射线结构光,第一摄像机和第二摄像机接收该线结构光的反射光,第二激光器向钢轨的另一个侧面发射线结构光,第三摄像机和第四摄像机接收该线结构光的反射光;
步骤4:第一摄像机经第一曝光时间生成第一图像,第二摄像机经第二曝光时间生成第二图像,第三摄像机经第三曝光时间生成第三图像,第四摄像机经第四曝光时间生成第四图像;
步骤5:根据振动传感器获取的振动信号预测伸缩缝的位置,识别第一图像和第三图像中伸缩缝的第一特征图,若第一特征图与第一基准特征图的结构相似度低于伸缩缝阈值,基于该第一特征图调整第一曝光时间和第三曝光时间;
步骤6:基于所述第二基准特征图识别第二图像和第四图像中的紧固件的第二特征图,基于该第二特征图调整第二曝光时间和第四曝光时间;
步骤7:重复步骤3至步骤4,再次生成第一图像、第二图像、第三图像以及第四图像,分别根据第一图像、第二图像、第三图像以及第四图像生成第一点云、第二点云、第三点云以及第四点云;
步骤8:合并第一点云、第二点云、第三点云以及第四点云生成钢轨的截面轮廓,识别该截面轮廓的磨耗值,返回至步骤3。
在本发明中,在步骤5中,提取振动信号的频谱特征和幅值特征,根据该频谱特征和幅值特征确定检测小车的车轮经过伸缩缝的时刻t1,根据时刻t1计算第一摄像机经过伸缩缝的时刻t2,提取第一摄像机在时刻t2生成的第一图像。
在本发明中,在步骤5中,结构相似度QSSIM=Q1×Q2×Q3,亮度参数,对比度参数,相关性参数,其中,/>为第一特征图的像素均值,/>为第一基准特征图的像素均值,σ1为第一特征图的像素标准差,σ2为第一基准特征图的像素标准差,σ12为第一特征图与第一基准特征图的像素协方差,C1、C2、C3均为常数。
在本发明中,在步骤5中,计算第一图像中第一特征图的曝光像素参数W1和曝光质量参数W2,并根据曝光像素参数和曝光质量参数计算第一图像参数A1,A1=(W1+W2)/2,若A1<A0,调整后的第一曝光时间E=E0+(1-QSSIM)ΔE,若A1>A0,调整后的第一曝光时间E=E0-(1-QSSIM)ΔE,E0为第一曝光时间,ΔE为曝光调节步长,A0为第一基准特征图的图像参数。
在本发明中,在步骤5中,曝光像素参数,曝光质量参数,p3为第一特征图的像素极差。
在本发明中,在步骤6中,在第二图像中提取多个预选窗口,计算每一预选窗口与第二基准特征图的直方图相似度,根据直方图相似度识别第二特征图,根据第二特征图计算第二图像参数并确定第二曝光时间。
在本发明中,在步骤6中,预选窗口和第二基准特征图的直方图相似度,p4i为预选窗口的直方图中像素值i的概率,p5i为第二基准特征图的直方图中像素值i的概率,/>为预选窗口的直方图中像素值的平均概率,/>为第二基准特征图的直方图中像素值的平均概率。
在本发明中,在步骤7中,标定第一摄像机的内部参数以及第一摄像机相对于第一激光器的外部参数,根据该内部参数和外部参数变换第一图像的每一像素点的位置坐标并生成第一点云,根据刚体变换估计将第一点云匹配至基准坐标系。
在本发明中,在步骤8中,所述磨耗值由垂直磨耗、内圆周磨耗以及外圆周磨耗组成。
一种应用于所述基于结构光的钢轨磨耗识别方法的钢轨磨耗识别设备,包括:
沿钢轨移动的检测小车,检测小车设有振动传感器;
数据库,数据库被配置为存储伸缩缝的第一基准特征图以及紧固件的第二基准特征图;
第一测量单元,第一测量单元被配置为生成第一图像和第二图像,第一测量单元由第一激光器、第一摄像机和第二摄像机组成;
第二测量单元,第二测量单元被配置为生成第三图像和第四图像,第二测量单元由第二激光器、第三摄像机和第四摄像机组成;
控制单元,控制单元被配置为分别根据第一图像、第二图像、第三图像以及第四图像调整第一曝光时间、第二曝光时间、第三曝光时间和第四曝光时间;
数据处理单元,数据处理单元被配置为生成钢轨的截面轮廓并识别该截面轮廓的磨耗值。
实施本发明的这种基于结构光的钢轨磨耗识别方法及设备,具有以下有益效果:针对钢轨不同角度的光线差异,本发明从现场实测图像中提取钢轨的第一特征图和第二特征图,预测现场实测的第一特征图和第二特征图的曝光像素参数和曝光质量参数,进而调整各个角度的摄像机的曝光时间,提高图像的质量,使得磨耗值的测量更准确。进一步的,针对钢轨顶部的第一图像和第三图像,结合检测小车的振动情况确认伸缩缝的位置,可以快速识别伸缩缝的第一特征图。针对钢轨侧面的第二图像和第四图像,伸缩缝被钢轨连接件阻挡,结合紧固件的纹理特征识别紧固件的第二特征图,可以显著提高图像的识别与曝光调节速度。
附图说明
图1为列车钢轨的示意图;
图2为本发明基于结构光的钢轨磨耗识别方法的流程图;
图3为本发明检测小车的示意图;
图4为本发明第一测量单元与第二测量单元的安装示意图;
图5为本发明提取第一特征图的示意图;
图6为本发明根据第二图像调整第二曝光时间的示意图;
图7为本发明截面轮廓的示意图;
图8为本发明根据第一图像调整第一曝光时间的示意图;
图9为本发明钢轨磨耗识别设备的框图。
附图中附图标记:101左侧钢轨、102右侧钢轨、103伸缩缝、104紧固件、105推行装置、106行进压紧装置、107车轮行进装置、108测量装置安装面板、109第一激光器、110第二激光器、111第一摄像机、112第二摄像机、113第三摄像机、114第四摄像机、115轨道接缝引起的周期性幅值特征、116轨道缺陷引起的非周期性幅值特征、117检测小车沿钢轨滚动引起的周期性幅值特征、118第一图像、119截面轮廓、120标准钢轨轮廓。
具体实施方式
为更清楚地理解本申请的目的、技术方案和优点,下面结合附图和实施例,对本申请进行了描述和说明。
由于轨道处于自然场景下,线结构光获取的数据受到的干扰特别多,导致数据的稳定性不好。本发明根据检测小车的的振动信号识别第一图像和第三图像中伸缩缝的第一特征图,基于该第一特征图调整第一曝光时间和第三曝光时间。根据第二基准特征图识别第二图像和第四图像中的紧固件的第二特征图,基于该第二特征图调整第二曝光时间和第四曝光时间。通过调节第一曝光时间和第二曝光时间、第三曝光时间和第四曝光时间提高第一图像、第二图像、第三图像以及第四图像的清晰度。根据第一图像、第二图像、第三图像以及第四图像生成第一点云、第二点云、第三点云以及第四点云,合并第一点云、第二点云、第三点云以及第四点云生成钢轨的截面轮廓,识别该截面轮廓的磨耗值。
实施例一
如图1至图7所示,本发明基于结构光的钢轨磨耗识别方法包括以下步骤。
步骤1:在检测小车上安装由第一激光器、第一摄像机和第二摄像机组成的第一测量单元和由第二激光器、第三摄像机和第四摄像机组成的第二测量单元以及振动传感器。如图3所示,检测小车包括推行装置105、行进压紧装置106、车轮行进装置107、测量装置安装面板108、光电旋转编码器。推行装置105用于提供检测小车推动力,车轮行进装置107用于保证检测小车沿着钢轨滚动,行进压紧装置106用于保证检测小车沿着钢轨滚动时不发生大的扭转或者滑移等现象,光电旋转编码器用于精确记录检测小车的行进里程。
步骤2:钢轨经紧固件安装在基座上,相邻的钢轨具有伸缩缝,读取伸缩缝的第一基准特征图,读取紧固件的第二基准特征图,所述检测小车沿钢轨移动。如图1所示,钢轨每间隔相同的距离都会出现伸缩缝103,在两个相邻的紧固件104之间。预设检测小车沿钢轨移动的速度,根据两个相邻的紧固件的距离和检测小车的速度设置检测小车的检测频率。
步骤3:第一激光器向钢轨的一个侧面发射线结构光,第一摄像机和第二摄像机接收该线结构光的反射光,第二激光器向钢轨的另一个侧面发射线结构光,第三摄像机和第四摄像机接收该线结构光的反射光。如图4所示,第一激光器109、第一摄像机111和第二摄像机112安装在检测小车的右侧。具体安装角度通过现场调试确定,以第一摄像机和第二摄像机能获取预定位置的图像为准。第二激光器110、第三摄像机113和第四摄像机114以对称于第一激光器109、第一摄像机111和第二摄像机112的位置安装在检测小车的左侧。
步骤4:第一摄像机经第一曝光时间生成第一图像,第二摄像机经第二曝光时间生成第二图像,第三摄像机经第三曝光时间生成第三图像,第四摄像机经第四曝光时间生成第四图像。反射光在第一摄像机内生成条纹图,将连续的条纹图组成第一图像。在较佳的曝光时间下,条纹图结构清晰,反射光和背景色在第一图像、第二图像、第三图像和第四图像中呈现鲜明的对比,图像质量高。
步骤5:根据振动传感器获取的振动信号预测伸缩缝的位置,识别第一图像和第三图像中伸缩缝的第一特征图,若第一特征图与第一基准特征图的结构相似度低于伸缩缝阈值,基于该第一特征图调整第一曝光时间和第三曝光时间。提取振动信号的周期特征和幅值特征,根据该周期特征和幅值特征确定检测小车的车轮经过伸缩缝的时刻t1,根据时刻t1计算第一摄像机经过伸缩缝的时刻t2,第一摄像机在时刻t2生成的第一图像即为第一特征图。
如图5所示,波形曲线为等时间间隔所采集的振动信号,横坐标为周期,纵坐标为幅值,曲线中出现了轨道接缝引起的周期性幅值特征115,轨道缺陷引起的非周期性幅值特征116,检测小车沿钢轨滚动引起的周期性幅值特征117,其中,轨道接缝引起的周期性幅值特征115大于检测小车沿钢轨滚动引起的周期性幅值特征117,从图中根据幅值特征提取相应的周期特征,根据周期特征获得检测小车的车轮经过伸缩缝的时刻t1,根据检测小车沿钢轨移动的速度v和第一摄像机到检测小车的车轮的距离d确定时间差Δt=d/v,时刻t2=t1+Δt,根据时刻t2在第一图像118中找到第一特征图。
结构相似度QSSIM=Q1×Q2×Q3,亮度参数,对比度参数,相关性参数,其中,/>为第一特征图的像素均值,/>为第一基准特征图的像素均值,σ1为第一特征图的像素标准差,σ2为第一基准特征图的像素标准差,σ12为第一特征图与第一基准特征图的像素协方差,C1、C2、C3均为常数,可随机选取接近0的正数,防止亮度参数、对比度参数、相关性参数的分母为零。基于第一特征图调整第一曝光时间和第三曝光时间的具体步骤在实施例二中进行详细说明。
步骤6:基于所述第二基准特征图识别第二图像和第四图像中的紧固件的第二特征图,基于该第二特征图调整第二曝光时间和第四曝光时间。如图6所示,选取扣件区域作为一个窗口,使用滑动窗口方法,在第二图像中提取多个预选窗口,计算每一预选窗口与第二基准特征图的直方图相似度,选择直方图相似度最大的预选窗口作为第二特征图,根据第二特征图计算第二图像参数并确定第二曝光时间。根据第二特征图确定第二曝光时间的具体步骤在实施例二中进行说明。预选窗口和第二基准特征图的直方图相似度,p4i为预选窗口的直方图中像素值i的概率,p5i为第二基准特征图的直方图中像素值i的概率,/>为预选窗口的直方图中像素值的平均概率,/>为第二基准特征图的直方图中像素值的平均概率。
步骤7:重复步骤3至步骤4,再次生成第一图像、第二图像、第三图像以及第四图像,分别根据第一图像、第二图像、第三图像以及第四图像生成第一点云、第二点云、第三点云以及第四点云。标定第一摄像机的内部参数以及第一摄像机相对于第一激光器的外部参数,根据该内部参数和外部参数变换第一图像的每一像素点的位置坐标并生成第一点云,根据刚体变换估计将第一点云匹配至基准坐标系。按照从第一图像到第一点云的操作,根据第二图像、第三图像以及第四图像得出第二点云、第三点云以及第四点云。具体来说,例如根据内部参数和外部参数确定第一图像的二维坐标系与钢轨三维坐标系的变换矩阵,将第一图像的每一像素点映射至三维坐标系,再填充丢失像素,得到第一点云。基于钢轨的形状确定刚体变换估计的参数,最后将四组点云拼在一起。
步骤8:合并第一点云、第二点云、第三点云以及第四点云生成钢轨的截面轮廓,识别该截面轮廓的磨耗值,返回至步骤3。合并第一点云、第二点云、第三点云以及第四点云生成钢轨的三维模型,提取三维模型的截面生成钢轨的截面轮廓。如图7所示,所述磨耗值由垂直磨耗、内圆周磨耗以及外圆周磨耗组成,将钢轨的截面轮廓119和标准钢轨轮廓120进行对比测量得出垂直磨耗L1、内圆周磨耗L2以及外圆周磨耗L3。磨耗值L=β1L1+β2L2+β3L3,β1为垂直磨耗的权重,β2为内圆周磨耗的权重,β3为外圆周磨耗的权重,一般来说,内圆周磨耗的权重大于垂直磨耗的权重,垂直磨耗的权重大于外圆周磨耗的权重。
实施例二
如图8所示,本发明进一步公开了调整第一曝光时间、第二曝光时间、第三曝光时间、第四曝光时间的方法。
根据振动信号从第一图像中选取第一特征图,计算第一特征图和第一基准特征图的结构相似度QSSIM,根据结构相似度QSSIM调整第一曝光时间。计算第一图像中第一特征图的曝光像素参数W1和曝光质量参数W2,并根据曝光像素参数和曝光质量参数计算第一图像参数A1,A1=(W1+W2)/2,若A1<A0,调整后的第一曝光时间E=E0+(1-QSSIM)ΔE,若A1>A0,调整后的第一曝光时间E=E0-(1-QSSIM)ΔE,E0为第一曝光时间,ΔE为曝光时间调节步长,可以设置为1毫秒,A0为第一基准特征图的图像参数,通常为0.9到1.2之间,可以为1.1。曝光像素参数,曝光质量参数/>,p3为第一特征图的像素极差,像素极差是指最大像素值与最小像素值之间的差值。第三曝光时间的调整方法参照第一曝光时间。
在第二图像中提取多个预选窗口,计算每一预选窗口与第二基准特征图的直方图相似度,选择直方图相似度最大的预选窗口作为第二特征图,根据QHIS调整第二曝光时间。计算第二图像中第二特征图的曝光像素参数W3和曝光质量参数W4,参照第一特征图的曝光像素参数和曝光质量参数的计算方法,并根据曝光像素参数W3和曝光质量参数W4计算第二图像参数A2,A2=(W3+W4)/2,若A2<A3,调整后的第二曝光时间E2=E1+(1-QHIS)ΔE,若A2>A3,调整后的第二曝光时间E2=E1-(1-QHIS)ΔE,A3为第二基准特征图的图像参数,E1为第二曝光时间。第四曝光时间的调整方法参照第二曝光时间。
实施例三
如图9所示,本发明应用于所述基于结构光的钢轨磨耗识别方法的钢轨磨耗识别设备,包括:检测小车、数据库、第一测量单元、第二测量单元、控制单元、数据处理单元。检测小车沿钢轨移动,用于检测钢轨的磨耗值,检测小车设有振动传感器、第一测量单元、第二测量单元、控制单元、数据库。数据库被配置为存储伸缩缝的第一基准特征图以及紧固件的第二基准特征图。第一测量单元由第一激光器、第一摄像机和第二摄像机组成,第一摄像机生成第一图像,第二摄像机生成第二图像。第二测量单元由第二激光器、第三摄像机和第四摄像机组成,第三摄像机生成第三图像,第四摄像机生成第四图像。控制单元被配置为分别根据第一图像、第二图像、第三图像以及第四图像调整第一曝光时间、第二曝光时间、第三曝光时间和第四曝光时间。数据处理单元被配置为根据第一图像、第二图像、第三图像以及第四图像生成第一点云、第二点云、第三点云、第四点云以及钢轨的截面轮廓,识别该截面轮廓的垂直磨耗、内圆周磨耗以及外圆周磨耗,根据垂直磨耗、内圆周磨耗以及外圆周磨耗计算磨耗值。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改,等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于结构光的钢轨磨耗识别方法,其特征在于,包括以下步骤:
步骤1:在检测小车上安装由第一激光器、第一摄像机和第二摄像机组成的第一测量单元和由第二激光器、第三摄像机和第四摄像机组成的第二测量单元以及振动传感器;
步骤2:钢轨经紧固件安装在基座上,相邻的钢轨具有伸缩缝,读取伸缩缝的第一基准特征图,读取紧固件的第二基准特征图,所述检测小车沿钢轨移动;
步骤3:第一激光器向钢轨的一个侧面发射线结构光,第一摄像机和第二摄像机接收该线结构光的反射光,第二激光器向钢轨的另一个侧面发射线结构光,第三摄像机和第四摄像机接收该线结构光的反射光;
步骤4:第一摄像机经第一曝光时间生成第一图像,第二摄像机经第二曝光时间生成第二图像,第三摄像机经第三曝光时间生成第三图像,第四摄像机经第四曝光时间生成第四图像;
步骤5:根据振动传感器获取的振动信号预测伸缩缝的位置,识别第一图像和第三图像中伸缩缝的第一特征图,若第一特征图与第一基准特征图的结构相似度低于伸缩缝阈值,基于该第一特征图调整第一曝光时间和第三曝光时间;
步骤6:基于所述第二基准特征图识别第二图像和第四图像中的紧固件的第二特征图,基于该第二特征图调整第二曝光时间和第四曝光时间;
步骤7:重复步骤3至步骤4,再次生成第一图像、第二图像、第三图像以及第四图像,分别根据第一图像、第二图像、第三图像以及第四图像生成第一点云、第二点云、第三点云以及第四点云;
步骤8:合并第一点云、第二点云、第三点云以及第四点云生成钢轨的截面轮廓,识别该截面轮廓的磨耗值,返回至步骤3。
2.根据权利要求1所述的基于结构光的钢轨磨耗识别方法,其特征在于,在步骤5中,提取振动信号的频谱特征和幅值特征,根据该频谱特征和幅值特征确定检测小车的车轮经过伸缩缝的时刻t1,根据时刻t1计算第一摄像机经过伸缩缝的时刻t2,提取第一摄像机在时刻t2生成的第一图像。
3.根据权利要求1所述的基于结构光的钢轨磨耗识别方法,其特征在于,在步骤5中,结构相似度QSSIM=Q1×Q2×Q3,亮度参数,对比度参数/>,相关性参数,其中,/>为第一特征图的像素均值,/>为第一基准特征图的像素均值,σ1为第一特征图的像素标准差,σ2为第一基准特征图的像素标准差,σ12为第一特征图与第一基准特征图的像素协方差,C1、C2、C3均为常数。
4.根据权利要求3所述的基于结构光的钢轨磨耗识别方法,其特征在于,在步骤5中,计算第一图像中第一特征图的曝光像素参数W1和曝光质量参数W2,并根据曝光像素参数和曝光质量参数计算第一图像参数A1,A1=(W1+W2)/2,若A1<A0,调整后的第一曝光时间E=E0+(1-QSSIM)ΔE,若A1>A0,调整后的第一曝光时间E=E0-(1-QSSIM)ΔE,E0为第一曝光时间,ΔE为曝光调节步长,A0为第一基准特征图的图像参数。
5.根据权利要求4所述的基于结构光的钢轨磨耗识别方法,其特征在于,在步骤5中,曝光像素参数,曝光质量参数/>,p3为第一特征图的像素极差。
6.根据权利要求1所述的基于结构光的钢轨磨耗识别方法,其特征在于,在步骤6中,在第二图像中提取多个预选窗口,计算每一预选窗口与第二基准特征图的直方图相似度,根据直方图相似度识别第二特征图,根据第二特征图计算第二图像参数并确定第二曝光时间。
7.根据权利要求6所述的基于结构光的钢轨磨耗识别方法,其特征在于,在步骤6中,预选窗口和第二基准特征图的直方图相似度,p4i为预选窗口的直方图中像素值i的概率,p5i为第二基准特征图的直方图中像素值i的概率,/>为预选窗口的直方图中像素值的平均概率,/>为第二基准特征图的直方图中像素值的平均概率。
8.根据权利要求1所述的基于结构光的钢轨磨耗识别方法,其特征在于,在步骤7中,标定第一摄像机的内部参数以及第一摄像机相对于第一激光器的外部参数,根据该内部参数和外部参数变换第一图像的每一像素点的位置坐标并生成第一点云,根据刚体变换估计将第一点云匹配至基准坐标系。
9.根据权利要求1所述的基于结构光的钢轨磨耗识别方法,其特征在于,在步骤8中,所述磨耗值由垂直磨耗、内圆周磨耗以及外圆周磨耗组成。
10.一种应用于如权利要求1所述基于结构光的钢轨磨耗识别方法的钢轨磨耗识别设备,其特征在于,包括:
沿钢轨移动的检测小车,检测小车设有振动传感器;
数据库,数据库被配置为存储伸缩缝的第一基准特征图以及紧固件的第二基准特征图;
第一测量单元,第一测量单元被配置为生成第一图像和第二图像,第一测量单元由第一激光器、第一摄像机和第二摄像机组成;
第二测量单元,第二测量单元被配置为生成第三图像和第四图像,第二测量单元由第二激光器、第三摄像机和第四摄像机组成;
控制单元,控制单元被配置为分别根据第一图像、第二图像、第三图像以及第四图像调整第一曝光时间、第二曝光时间、第三曝光时间和第四曝光时间;
数据处理单元,数据处理单元被配置为生成钢轨的截面轮廓并识别该截面轮廓的磨耗值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410650265.XA CN118224999B (zh) | 2024-05-24 | 2024-05-24 | 一种基于结构光的钢轨磨耗识别方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410650265.XA CN118224999B (zh) | 2024-05-24 | 2024-05-24 | 一种基于结构光的钢轨磨耗识别方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN118224999A true CN118224999A (zh) | 2024-06-21 |
CN118224999B CN118224999B (zh) | 2024-07-23 |
Family
ID=91505276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410650265.XA Active CN118224999B (zh) | 2024-05-24 | 2024-05-24 | 一种基于结构光的钢轨磨耗识别方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118224999B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118379640A (zh) * | 2024-06-26 | 2024-07-23 | 山东光安智能科技有限公司 | 一种粮仓智能巡检方法、系统、设备、介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000136988A (ja) * | 1998-10-30 | 2000-05-16 | East Japan Railway Co | レール波状摩耗検出手法 |
EP1496343A2 (de) * | 2003-07-08 | 2005-01-12 | Holger Dipl.-Ing. Wente | Verfahren und Messeinrichtung zur Aufzeichnung der Bewegungsspur von Beobachtungspunkten an schwingenden Gegenständen |
CN101576375A (zh) * | 2009-05-21 | 2009-11-11 | 北京航空航天大学 | 一种钢轨磨耗激光视觉图像快速处理方法 |
CN106595596A (zh) * | 2016-12-14 | 2017-04-26 | 中南大学 | 一种铁路轨道轨向检测的摄影测量方法 |
KR101742981B1 (ko) * | 2016-06-03 | 2017-06-14 | 대구도시철도공사 | 철도레일 마모 측정장치 |
CN111122598A (zh) * | 2019-12-16 | 2020-05-08 | 北京冶自欧博科技发展有限公司 | 一种钢轨表面缺陷三维检测系统及方法 |
CN111623868A (zh) * | 2020-05-12 | 2020-09-04 | 西南交通大学 | 用于钢轨波磨识别的卷积神经网络构建方法 |
US20200334895A1 (en) * | 2017-12-22 | 2020-10-22 | Zivid Labs As | Method and system for generating a three-dimensional image of an object |
CN114575205A (zh) * | 2022-04-28 | 2022-06-03 | 中铁第四勘察设计院集团有限公司 | 一种基于图像数据处理的水射流钢轨廓形智能化打磨系统 |
-
2024
- 2024-05-24 CN CN202410650265.XA patent/CN118224999B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000136988A (ja) * | 1998-10-30 | 2000-05-16 | East Japan Railway Co | レール波状摩耗検出手法 |
EP1496343A2 (de) * | 2003-07-08 | 2005-01-12 | Holger Dipl.-Ing. Wente | Verfahren und Messeinrichtung zur Aufzeichnung der Bewegungsspur von Beobachtungspunkten an schwingenden Gegenständen |
CN101576375A (zh) * | 2009-05-21 | 2009-11-11 | 北京航空航天大学 | 一种钢轨磨耗激光视觉图像快速处理方法 |
KR101742981B1 (ko) * | 2016-06-03 | 2017-06-14 | 대구도시철도공사 | 철도레일 마모 측정장치 |
CN106595596A (zh) * | 2016-12-14 | 2017-04-26 | 中南大学 | 一种铁路轨道轨向检测的摄影测量方法 |
US20200334895A1 (en) * | 2017-12-22 | 2020-10-22 | Zivid Labs As | Method and system for generating a three-dimensional image of an object |
CN111122598A (zh) * | 2019-12-16 | 2020-05-08 | 北京冶自欧博科技发展有限公司 | 一种钢轨表面缺陷三维检测系统及方法 |
CN111623868A (zh) * | 2020-05-12 | 2020-09-04 | 西南交通大学 | 用于钢轨波磨识别的卷积神经网络构建方法 |
CN114575205A (zh) * | 2022-04-28 | 2022-06-03 | 中铁第四勘察设计院集团有限公司 | 一种基于图像数据处理的水射流钢轨廓形智能化打磨系统 |
Non-Patent Citations (1)
Title |
---|
赵静: "机器视觉下的钢轨廓形提取与磨耗测量", 《测绘通报》, no. 12, 31 December 2023 (2023-12-31), pages 159 - 163 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118379640A (zh) * | 2024-06-26 | 2024-07-23 | 山东光安智能科技有限公司 | 一种粮仓智能巡检方法、系统、设备、介质 |
Also Published As
Publication number | Publication date |
---|---|
CN118224999B (zh) | 2024-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN118224999B (zh) | 一种基于结构光的钢轨磨耗识别方法及设备 | |
US6909514B2 (en) | Wheel profile inspection apparatus and method | |
CN108986082A (zh) | 一种基于epnp的钢轨廓形检测方法及系统 | |
CN107014294B (zh) | 一种基于红外图像的接触网几何参数检测方法及系统 | |
EP3196863A1 (en) | System and method for aircraft docking guidance and aircraft type identification | |
WO2016015547A1 (zh) | 一种基于机器视觉的飞机入坞引导和机型识别的方法及系统 | |
CN108725511B (zh) | 钢轨波磨弦测点的实时位置修正方法 | |
CN115482195B (zh) | 一种基于三维点云的列车部件变形检测方法 | |
CN111539436B (zh) | 一种基于一字模板匹配的轨道扣件定位方法 | |
Xing et al. | Online detection system for wheel-set size of rail vehicle based on 2D laser displacement sensors | |
Liu et al. | Recognition and calibration of rail profile under affine-distortion-based point set mapping | |
CN109507198B (zh) | 基于快速傅里叶变换和线性高斯的口罩检测系统及方法 | |
CN113324478A (zh) | 一种线结构光的中心提取方法及锻件三维测量方法 | |
US10408610B1 (en) | Method and system for displacement measurement of surfaces on a moving vehicle | |
CN110954005B (zh) | 一种基于图像处理的中低速磁悬浮列车悬浮间隙检测方法 | |
CN110458785B (zh) | 一种基于图像传感的磁悬浮球悬浮间隙检测方法 | |
CN112408153B (zh) | 电梯制动过程钢丝绳移动量监测方法 | |
CN112950562A (zh) | 一种基于线结构光的扣件检测算法 | |
Molleda et al. | A profile measurement system for rail manufacturing using multiple laser range finders | |
CN111551122A (zh) | 一种基于激光雷达的火车车皮数量和长度测量系统及方法 | |
CN114426039B (zh) | 一种基于城轨列车车轮踏面的点云数据处理方法 | |
RU2650857C1 (ru) | Система определения геометрических параметров трехмерных объектов | |
CN115187536A (zh) | 一种基于机器视觉的钢轨波磨光截图像测量方法及设备 | |
Xu et al. | On-line visual inspection system for backside weld of tailored blanks laser welding | |
Wang et al. | An accurate method for the extraction of line structured light stripe. |
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 |