一种基于摄像的火车轨道松动非接触实时检测方法技术

针对当前火车轨道松动主要靠人工巡检效率低的问题，本发明专利技术公开了一种基于摄像的火车轨道松动非接触实时检测方法。该方法用高速摄像机获取轨道在火车通过或者人为激励下振动视频，利用火车轨道和枕木间交错排布的特征，提出了基于图像像素灰度梯度和特征聚类算法的虚拟特征点检测方法，采用光流算法对虚拟特征点的光流进行计算，进而对轨道的时域振动进行测量。通过FFT分解得到轨道固有振动频率，进而通过固有振动频率的变化实时判断轨道是否发生了松动。所提方法测量装置简单、精度高、成本较低、易于操作。易于操作。易于操作。

【技术实现步骤摘要】
一种基于摄像的火车轨道松动非接触实时检测方法


[0001]本专利技术属于结构运动测量
，特别涉及一种基于摄像的火车轨道松动非接触实时检测方法。

技术介绍

[0002]火车轨道系统是交通运输系统中的重要部分，主要由轨道、枕木、扣件和路基组成。轨道由扣件固定在枕木上，在列车行驶过程中，其产生周期性冲击动载荷容易造成扣件发生振动，长此以往，会引发轨道松动。随着轨道松动的加剧，轨道结构的动力响应振幅明显变大，严重时会造成列车脱轨事故。因此，火车轨道松脱检测对保障列车行驶安全具有重大意义。
[0003]目前，常用的火车轨道松动检测方法主要为人工巡检，通常由经验丰富的巡检人员用肉眼判断扣件状态。该方法操作简单、效率低、成本高、漏检率高、安全隐患大。近年来，随着图像处理技术的发展，计算机视觉技术也逐渐被应用于火车轨道的结构健康监测。目前，使用的计算机视觉的检测方法需要对结构进行标记，通过跟踪标记点的方式对轨道振动进行测量。然而，由于火车轨道数量多，且长期处于室外环境，标记点容易脱落，给测量带来极大的挑战。

技术实现思路

[0004]为了克服上述现有技术的缺点，解决基于人工巡检的火车轨道松动检测方法效率、精度低的问题，本专利技术的目的在于提供一种基于摄像的火车轨道松动非接触实时检测方法，利用火车轨道和枕木间交错排布的特征，提出了基于图像像素灰度梯度和特征聚类算法的虚拟特征点检测方法，采用光流算法对虚拟特征点的光流进行计算，进而对轨道的时域振动进行测量，通过FFT分解得到轨道固有振动频率，进而通过固有振动频率的变化实时判断轨道是否发生了松动。
[0005]为了实现上述目的，本专利技术采用的技术方案是：
[0006]一种基于摄像的火车轨道松动非接触实时检测方法，包括以下步骤：
[0007]步骤1)，利用高速摄像机对火车轨道的振动进行视频记录；
[0008]步骤2)，对记录的视频，逐帧对图像进行灰度化处理；
[0009]步骤3)，计算第一帧图像像素的灰度梯度，采用K
‑
means聚类算法对灰度梯度进行聚类，确定轨道和扣件区域，将区域内像素灰度极值作为该区域的虚拟特征点；
[0010]步骤4)，通过图像多尺度分解技术构建图像多尺度金字塔，基于短时亮度恒定理论建立不同尺度上图像的光流方程，采用最小二乘算法计算不同尺度图像上虚拟特征点的光流；
[0011]步骤5)，利用金字塔图像间的尺度关系，对不同尺度图像上虚拟特征点的光流进行融合，得到虚拟特征点光流的计算结果，基于图像标定技术，获取轨道和扣件时域振动信号；
[0012]步骤6)，对轨道和扣件时域振动信号进行频域分析，得到轨道的固有频率，通过固有频率变化特征对轨道松动进行检测。
[0013]进一步的，所述步骤2具体为，为提高虚拟特征点的筛选效率，对相机获取的彩色图像逐帧进行灰度化处理，获取灰度图像：
[0014]I9x,y)＝0.299*R9x,y)+0.579*G9x,y)+0.114*B(x,y)
[0015]式中：I(x,y)为像素(x,y)的灰度值；R(x,y)、G(x,y)、B(x,y)分别为像素(x,y)的三个通道的像素值。
[0016]进一步的，所述步骤3具体为，基于图像像素灰度梯度和特征聚类算法对虚拟特征点进行筛选。
[0017]在本专利技术中，设第一帧图像像素所计算的水平方向上的梯度集合为采用K
‑
means聚类算法对水平方向上的梯度进行聚类分析：
[0018][0019]式中：E
x
为水平方向上分类平方和的计算结果，其值最小时，为最佳聚类；v为聚类的种类，在本专利技术中v为2，v＝1代表轨道区域，v＝2代表其他区域；表示水平方向梯度均值；l为图像像素的数目。
[0020]设第一帧图像像素所计算的竖直方向上的梯度集合为采用K
‑
means聚类算法对竖直方向上的梯度进行聚类分析：
[0021][0022]式中：E
y
为竖直方向上分类平方和的计算结果，其值最小时，为最佳聚类；h为聚类的种类，在本专利技术中h为2，h＝1代表扣件和枕木区域，h＝2代表其他区域；表示竖直方向梯度均值；l为图像像素的数目。
[0023]本专利技术中，设检测出的轨道和扣件区域记为Ω＝(Ω1,Ω2,...,Ω
q
)，q为确定的轨道和扣件的区域数目。利用不同区域内像素灰度极值点作为该区域的虚拟特征点。对于区域Ω1，该区域内像素灰度极值表示为：
[0024][0025]式中：为Ω1区域内灰度值最大的点的像素坐标，在本专利技术中，将该点作为Ω1区域虚拟特征点，虚拟特征点的灰度值，I(x,y)为像素(x,y)位置处的灰度值，n2为Ω1区域内的像素数目。
[0026]对不同位置处的轨道和扣件区域内的灰度极值点进行检测，获得不同位置处的虚拟特征点。
[0027]进一步的，所述步骤4具体为，对于虚拟特征点(k
x
,k
y
)，基于短时亮度恒定理论及空间一致性假设得到不同尺度下图像的光流方程。选取虚拟特征点邻域窗口大小m
×
m，根据邻域内像素运动一致性原则，虚拟特征点邻域内m
×
m个像素满足：
[0028][0029]式中：u,v为虚拟特征点(k
x
,k
y
)在水平和竖直方向上的光流，(k
x
‑
m,k
y
‑
m),
…
,(k
x
+m,k
y
+m)分别为虚拟特征点(k
x
,k
y
)的邻域内像素坐标，I
y
(k
x
+m,k
y
+m)和I
y
(k
x
+m,k
y
+m)分别表示虚拟特征点邻域内像素的灰度在x和y方向上的梯度，I
t
(k
x
+m,k
y
+m)表示虚拟特征点邻域内像素的灰度关于时间t的导数。
[0030]考虑上述方程为超静定方程，在本专利技术中，采用最小二乘法对方程进行求解，可得任意尺度下虚拟特征点的光流信息。
[0031]进一步的，所述步骤5具体为，根据图像间的尺度关系，将不同尺度下虚拟特征点的光流信息进行融合，得到虚拟特征点光流。利用虚拟特征点的光流，计算像素坐标系下的轨道和扣件的振动。
[0032]设不同帧虚拟特征点的光流在水平和竖直方向上分别为{u
k
,v
k
|k＝1,2,3,...,K}，单位为像素/帧，K为视频的总帧数，通过虚拟特征点的光流可以得到结构的运动：
[0033][0034]式中：m
x
、m
y
分别为虚拟特征点在水平和竖直方向上的运动，f为拍摄帧率。
...

【技术保护点】

【技术特征摘要】
1.一种基于摄像的火车轨道松动非接触实时检测方法，其特征在于，包括以下步骤：步骤1)，利用高速摄像机对火车轨道的振动进行视频记录；步骤2)，对记录的视频，逐帧对图像进行灰度化处理；步骤3)，计算第一帧图像像素的灰度梯度，采用K
‑
means聚类算法对灰度梯度进行聚类，确定轨道和扣件区域，将区域内像素灰度极值作为该区域的虚拟特征点；步骤4)，通过图像多尺度分解技术构建图像多尺度金字塔，基于短时亮度恒定理论建立不同尺度上图像的光流方程，采用最小二乘算法计算不同尺度图像上虚拟特征点的光流；步骤5)，利用金字塔图像间的尺度关系，对不同尺度图像上虚拟特征点的光流进行融合，得到虚拟特征点光流的计算结果，结合图像标定技术，获取轨道和扣件时域振动信号；步骤6)，对轨道和扣件时域振动信号进行频域分析，得到轨道的固有频率，通过固有频率变化特征对轨道松动进行检测。2.根据权利要求1所述基于摄像的火车轨道松动非接触实时检测方法，其特征在于，所述步骤3)中，设第一帧图像像素所计算的水平方向上的灰度梯度集合为采用K
‑
means聚类算法对水平方向上的灰度梯度聚类分析：式中：E
x
为水平方向上分类平方和的计算结果，其值最小时，为最佳聚类；v为聚类的种类，v＝1代表轨道区域，v＝2代表其他区域；表示水平方向梯度均值；l为图像像素的数目；设第一帧图像像素所计算的竖直方向上的灰度梯度集合为采用K
‑
means聚类算法对竖直方向上的灰度梯度进行聚类分析：式中：E
y
为竖直方向上分类平方和的计算结果，其值最小时，为最佳聚类；h为聚类的种类，h＝1代表扣件和枕木区域，h＝2代表其他区域；表示竖直方向梯度均值；l为图像像素的数目；将检测出的轨道和扣件区域记为Ω＝(Ω1,Ω2,...,Ω
q
)，q为确定的轨道和扣件的区域数目，利用各个区域内像素灰度极值点作为该区域的虚拟特征点，对于区域Ω1，该区域内像素灰度极值表示为：式中：为Ω1区域内灰度值最大的点的像素坐标，将该点作为Ω1区域虚拟特征点；虚拟特征点的灰度值；I(x,y)为像素(x,y)位置处的灰度值；n2为Ω1区域内的像素数目；对不同位置处的轨道和扣件区域内的灰度极值点进行检测，获得不同位置处的虚拟特
征点。3.根据权利要求1所述基于摄像的火车轨道松动非接触实时检测方法，其特征在于，所述步骤4)中，对于虚拟特征点(k
x
,k
y
)，基于短时亮度恒定理论及空间一致性假设得到不同尺度图像的光流方程，选取虚拟特征点邻域窗口大小为m
×
m，根据邻域内像素运动一致性原则，虚拟特征点邻域内m...

【专利技术属性】
技术研发人员：徐自力，辛存，王存俊，李康迪，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：