作者簡(jiǎn)介：李培云（1995—），女，安徽滁州人，碩士生，主要研究方向?yàn)闄C(jī)器視覺(jué)，E-mail：1125924814@qq.com。

摘要：實(shí)驗(yàn)結(jié)果表明：該算法能在大多數(shù)情況下實(shí)時(shí)、準(zhǔn)確地檢測(cè)跟蹤接觸線，為城市軌道安全運(yùn)行提供了保障。

0   引言

隨著我國(guó)電氣化鐵路不斷地往高速、重載方向發(fā)展，使接觸網(wǎng)- 受電弓在高速運(yùn)行工作狀態(tài)下越來(lái)越受到關(guān)注^[1-2]。由于接觸網(wǎng)長(zhǎng)期保持和受電弓滑板相互接觸，在機(jī)車(chē)行駛過(guò)程中時(shí)刻磨損，從而給機(jī)車(chē)取電的性能帶來(lái)了很多不可靠的因素，直接影響機(jī)車(chē)的行駛狀態(tài)，嚴(yán)重時(shí)會(huì)導(dǎo)致發(fā)生安全事故。因此需要實(shí)時(shí)檢測(cè)接觸網(wǎng)的運(yùn)行狀態(tài)，以保障機(jī)車(chē)安全行駛。

近年來(lái)，國(guó)內(nèi)外對(duì)接觸網(wǎng)運(yùn)行狀態(tài)監(jiān)測(cè)受到了研究學(xué)者的廣泛關(guān)注。文獻(xiàn)^[3]基于識(shí)別出接觸線的導(dǎo)高、拉出值以及硬點(diǎn)部位，提出了一種基于圖像處理的鐵路接觸網(wǎng)檢測(cè)系統(tǒng)。文獻(xiàn)^[4]利用高精度傳感器采集線夾溫度等參數(shù)，提出了一種基于張力與溫度的接觸網(wǎng)安全狀態(tài)在線監(jiān)測(cè)系統(tǒng)。文獻(xiàn)^[5] 研究了一套牽引供電接觸網(wǎng)溫度在線監(jiān)測(cè)系統(tǒng)，但只是對(duì)接觸網(wǎng)某些關(guān)鍵部位進(jìn)行溫度檢測(cè)。文獻(xiàn)^[6]西安理工大學(xué)熊磊開(kāi)發(fā)了一套基于紅外熱像儀的接觸網(wǎng)載流溫度監(jiān)測(cè)管理系統(tǒng)，對(duì)實(shí)現(xiàn)自動(dòng)化監(jiān)測(cè)具有重要意義。

由于機(jī)車(chē)在行駛過(guò)程中遇到道路岔口需要頻繁換線，為了解決在行車(chē)過(guò)程中弓網(wǎng)脫離引發(fā)的安全事故，本文提出了一種接觸網(wǎng)目標(biāo)檢測(cè)與跟蹤算法，能夠?qū)崟r(shí)、快速地檢測(cè)接觸網(wǎng)，并能準(zhǔn)確地跟蹤接觸線。

1   算法概述

本文算法主要分為接觸線檢測(cè)和參數(shù)提取并跟蹤兩部分。檢測(cè)部分首先對(duì)輸入的幀做灰度化處理，然后分別運(yùn)用OTSU 自適應(yīng)閾值分割和邊緣檢測(cè)兩種方法提取接觸網(wǎng)的邊緣特征，最后利用Hough 變換連接直線邊緣。由于在紅外視頻圖像序列中，接觸線和承力索絕大多數(shù)都是成雙成對(duì)出現(xiàn)的，且在視頻畫(huà)面上都是呈直線型的，所以接觸線和承力索在視覺(jué)上無(wú)顯著特征，從而對(duì)目標(biāo)接觸線的跟蹤具有一定的難度。因此，跟蹤部分主要是利用Hough 變換提取目標(biāo)接觸線的參數(shù)和利用卡爾曼濾波方法跟蹤參數(shù)兩個(gè)部分。參數(shù)提取部分首先利用Hough 變換法記錄候選直線參數(shù)，根據(jù)采集到的紅外視頻，發(fā)現(xiàn)采集的視頻圖像中接觸網(wǎng)導(dǎo)線出現(xiàn)在畫(huà)面中近似直線，且一般是2~4 根，所以需要根據(jù)直線數(shù)目分情況討論，然后進(jìn)行判別并提取接觸網(wǎng)目標(biāo)直線參數(shù)；參數(shù)跟蹤部分通過(guò)卡爾曼濾波器預(yù)測(cè)紅外圖像中接觸線的參數(shù)，再求取實(shí)際接觸線的測(cè)量參數(shù)，并將參數(shù)與接觸線進(jìn)行一一配對(duì)，最終實(shí)現(xiàn)了紅外視頻中接觸線的檢測(cè)與跟蹤。一種接觸網(wǎng)目標(biāo)檢測(cè)與跟蹤算法流程如圖1 所示。

2   接觸線檢測(cè)

2.1 圖像預(yù)處理

首先，將紅外圖像（見(jiàn)圖2a）轉(zhuǎn)換為灰度圖像（見(jiàn)圖2b），可以增強(qiáng)圖像中的特征信息，并提高圖像處理的工作效率，同時(shí)也便于對(duì)圖像進(jìn)一步處理（如：圖像去噪、圖像分割等）。其次， 利用大津法自適應(yīng)求得圖像閾值，進(jìn)行分割后得到二值圖像（見(jiàn)圖2c）。大津法是圖像分割中閾值選取的最佳算法^[7]，運(yùn)算簡(jiǎn)單，不受圖像亮度和對(duì)比度的影響，因此在數(shù)字圖像處理上得到了廣泛的應(yīng)用。最后，為了獲取圖像的邊緣信息， 采用prewitt邊緣檢測(cè)算法^[8]提取接觸網(wǎng)的邊緣特征信息，得到二值化的邊緣圖像（見(jiàn)圖2d）。Prewitt 邊緣檢測(cè)算子對(duì)噪聲具有一定的抑制能力，對(duì)噪聲較多和灰度漸變的圖像處理產(chǎn)生了較好的邊緣效果。充分利用OTSU 自適應(yīng)閾值分割和Prewitt 算子在紅外圖像邊緣檢測(cè)中的各自?xún)?yōu)勢(shì)進(jìn)行特征提取，相比傳統(tǒng)的邊緣檢測(cè)算法，具有邊緣細(xì)節(jié)更清晰、背景噪聲更小等優(yōu)點(diǎn)，能夠更完整、更精確地提取邊緣信息。

（a）紅外圖像

（b）灰度圖像

（c）Otsu分割

（d）Prewitt算子

（e）log算子

（f）sobel算子

（g）kirsch算子

（h）roberts算子

圖2 紅外圖像預(yù)處理效果圖示

2.2 接觸線的邊緣連接

圖像采集的視頻是由紅外攝像機(jī)拍攝的，采集的圖像容易受自然界中的光照和機(jī)車(chē)運(yùn)行過(guò)程中的強(qiáng)電流、電壓磁場(chǎng)等干擾，導(dǎo)致獲取的圖像具有很多噪聲。由于接觸網(wǎng)導(dǎo)線和受電弓滑板相互摩擦產(chǎn)生磨損等原因，使得邊緣檢測(cè)后的直線邊緣出現(xiàn)中斷或者產(chǎn)生了很多偽邊緣。因此，常見(jiàn)的做法是在邊緣檢測(cè)后緊跟著做邊緣連接^[9]，形成連續(xù)的、有意義的邊緣。由于接觸線和承力索基本上都是成雙成對(duì)出現(xiàn)，且絕大部分情況都是呈直線型地呈現(xiàn)在視頻畫(huà)面中，因此本文選用Hough 變換[10-12]連接接觸線的邊緣，優(yōu)點(diǎn)是對(duì)噪聲不敏感，能較好地處理圖像中局部被遮擋和覆蓋的情況。Hough 變換法邊緣連接結(jié)果如圖3 所示。

（a）紅外圖像

（b）Otsu+Prewitt邊緣檢測(cè)

（c）邊緣連接結(jié)果

（d）邊緣連接的局部放大圖

圖3 邊緣連接實(shí)驗(yàn)結(jié)果圖

3   接觸線跟蹤

視頻圖像中接觸網(wǎng)安全運(yùn)行是機(jī)車(chē)行駛的可靠性和安全性的關(guān)鍵，在電氣化鐵路高速、重載發(fā)展中占據(jù)非常重要的位置。由于接觸線和承力索絕大多數(shù)都是成雙成對(duì)出現(xiàn)且判別性較弱，在接觸線跟蹤過(guò)程中出現(xiàn)道岔口需要換線或者橫支架部分遮擋等情況時(shí)，對(duì)目標(biāo)接觸線的跟蹤具有一定的難度。因此，提出在紅外視頻的各幀中利用Hough 變換檢測(cè)直線的優(yōu)勢(shì)來(lái)提取接觸線的參數(shù)，進(jìn)而使用Kalman 濾波實(shí)現(xiàn)對(duì)接觸線參數(shù)一幀幀地持續(xù)跟蹤，從而彌補(bǔ)接觸線和承力索在視覺(jué)上無(wú)顯著特征的不足。

3.1 接觸線的參數(shù)提取

Hough 變換是一種估計(jì)直線參數(shù)的方法，其抗噪性能力強(qiáng) , 對(duì)直線斷裂、局部遮擋等缺陷不敏感^[13]。接觸線的檢測(cè)問(wèn)題，經(jīng)過(guò)Hough 變換變成了對(duì)參數(shù)空間中半徑ρ 參數(shù)和角度θ 參數(shù)的求取過(guò)程。通過(guò)對(duì)候選直線參數(shù)的ρ 和θ 進(jìn)行分析，發(fā)現(xiàn)圖像中直線對(duì)應(yīng)的點(diǎn)皆在第一象限中，且目標(biāo)接觸線的ρ 和θ 在參數(shù)空間（極坐標(biāo)系）中的數(shù)值是最大值。根據(jù)采集到的紅外視頻中受電弓與接觸線的運(yùn)行狀態(tài)，發(fā)現(xiàn)采集的視頻圖像中接觸網(wǎng)導(dǎo)線出現(xiàn)在畫(huà)面中呈直線型，且直線數(shù)目一般是2-4根，所以需要根據(jù)直線數(shù)目分情況討論，然后進(jìn)行判別并提取目標(biāo)直線參數(shù)即接觸線參數(shù)。由此可知，可以將圖像空間中較為困難的全局檢測(cè)問(wèn)題轉(zhuǎn)化為參數(shù)空間中相對(duì)容易解決的局部峰值檢測(cè)問(wèn)題^[14]。

3.2 接觸線的參數(shù)跟蹤

本文利用卡爾曼濾波器^[15]預(yù)測(cè)接觸線的參數(shù)來(lái)估計(jì)出其目標(biāo)接觸線的軌跡位置并對(duì)得到的參數(shù)數(shù)據(jù)進(jìn)行修正。在采集的紅外視頻幀中，將目標(biāo)接觸線理解為一個(gè)包含參數(shù)空間中半徑和角度的向量，從當(dāng)前幀通過(guò)協(xié)方差矩陣推測(cè)出下一幀的參數(shù)預(yù)測(cè)值，即估計(jì)出下一幀中接觸線的軌跡位置。確定卡爾曼增益系數(shù)對(duì)參數(shù)的估計(jì)值和預(yù)測(cè)誤差進(jìn)行修正，接著進(jìn)入下一次迭代過(guò)程?？柭鼮V波預(yù)測(cè)接觸線軌跡算法流程如圖4 所示。

假設(shè)紅外視頻中目標(biāo)接觸線在s 幀的參數(shù)可以用一個(gè)向量Peak(s,:) 來(lái)表示，其中rho_max 表示其接觸線參數(shù)的極坐標(biāo)系下半徑最大值，theta_max 表示其接觸線參數(shù)的極坐標(biāo)系下角度最大值。

如果已知s-1 幀目標(biāo)接觸線的參數(shù)預(yù)測(cè)矩陣即將向量Peak 測(cè)量值賦值給x，s 幀接觸線的參數(shù)觀測(cè)矩陣可以用以下公式來(lái)表示

式中：x1(s,:) 是在第s 幀接觸線參數(shù)的預(yù)測(cè)矩陣，A 表示由s-1 幀最優(yōu)估計(jì)x(s-1) 推測(cè)出第s 幀接觸線軌跡的參數(shù)觀測(cè)矩陣，w 為預(yù)測(cè)模型的過(guò)程噪聲，B 為控制矩陣，表示控制量w 如何作用于當(dāng)前狀態(tài)。由于是對(duì)紅外視頻中目標(biāo)接觸線參數(shù)的預(yù)測(cè)，這里通過(guò)一個(gè)均值為0 的高斯噪聲來(lái)模擬預(yù)測(cè)時(shí)產(chǎn)生的誤差。視頻中每一幀的不確定性都是通過(guò)狀態(tài)噪聲協(xié)方差矩陣來(lái)表示，預(yù)測(cè)參數(shù)時(shí)必定是不準(zhǔn)確的，本文通過(guò)狀態(tài)噪聲協(xié)方差矩陣Q 來(lái)表示預(yù)測(cè)過(guò)程中產(chǎn)生的誤差，該誤差包含一些不確定因素，預(yù)測(cè)誤差協(xié)方差矩陣為

式中，p1(s,:) 為下一幀接觸線位置參數(shù)的預(yù)測(cè)誤差協(xié)方差矩陣，p(s-1,:) 為前一幀位置參數(shù)最優(yōu)估計(jì)的觀測(cè)誤差協(xié)方差矩陣，A 為上一幀最優(yōu)估計(jì)預(yù)測(cè)第s 幀接觸線軌跡的參數(shù)轉(zhuǎn)移矩陣的轉(zhuǎn)置矩陣。在紅外視頻中，接觸線在實(shí)際檢測(cè)中無(wú)論如何都避免不了誤差，通過(guò)Hough 變換檢測(cè)的接觸線嚴(yán)格來(lái)說(shuō)也只是對(duì)接觸線參數(shù)的預(yù)測(cè)估計(jì)。Kalman 由預(yù)測(cè)過(guò)程和測(cè)量修正兩部分組成，估計(jì)過(guò)程中包含接觸線當(dāng)前位置參數(shù)的預(yù)測(cè)和誤差協(xié)方差的預(yù)測(cè)。濾波器修正部分主要是對(duì)預(yù)測(cè)狀態(tài)的更新，包括Kalman 濾波增益的更新，利用K(s) 對(duì)位置參數(shù)值和協(xié)方差修正。

卡爾曼增益表達(dá)式為

式中，k(s) 為第s 幀的卡爾曼濾波常數(shù)，R 是參數(shù)觀測(cè)噪聲協(xié)方差矩陣，即參數(shù)的測(cè)量誤差，H 是真實(shí)參數(shù)到預(yù)測(cè)參數(shù)的預(yù)估觀測(cè)值的參數(shù)轉(zhuǎn)移矩陣。

推測(cè)出接觸線下一幀位置參數(shù)的最佳估計(jì)值的數(shù)學(xué)表達(dá)式為

式中，x1(s,:) 是第s 幀接觸線實(shí)際軌跡的參數(shù)推測(cè)值，y(s,:) 是第s 幀接觸線參數(shù)的實(shí)際測(cè)量值即參數(shù)觀測(cè)值。

最后，更新紅外視頻圖像中接觸線的實(shí)際位置參數(shù)與檢測(cè)位置參數(shù)之間的誤差協(xié)方差矩陣，即更新預(yù)測(cè)誤差協(xié)方差矩陣，數(shù)學(xué)表達(dá)為

由上面分析可知，卡爾曼濾波所用信息都是時(shí)域內(nèi)的量, 所以卡爾曼濾波器是在時(shí)域內(nèi)設(shè)計(jì)的, 且適用于多維情況。本文采用卡爾曼濾波器預(yù)測(cè)當(dāng)前視頻幀接觸線的參數(shù)來(lái)估計(jì)接觸線的軌跡位置，如果在紅外視頻某一幀中沒(méi)有檢測(cè)到接觸線的參數(shù)，同時(shí)又沒(méi)有預(yù)測(cè)到接觸線的參數(shù)時(shí)，則舍棄記錄該視頻幀接觸線軌跡位置。

4   實(shí)驗(yàn)結(jié)果

本文實(shí)驗(yàn)采用MATLAB R2018a 進(jìn)行仿真，運(yùn)行環(huán)境為：CPU：Intel Core i5-8250U@ 1.6 GHz 1.8 GHz，內(nèi)存：8G，Windows10 操作系統(tǒng)。使用MATLAB 仿真實(shí)現(xiàn)MP4 接觸網(wǎng)視頻文件到序列幀的轉(zhuǎn)換，如圖5 所示。以此序列圖片為例，成功實(shí)現(xiàn)了對(duì)接觸線的跟蹤。實(shí)驗(yàn)結(jié)果表明該算法準(zhǔn)確性好，跟蹤效果顯著，在電氣化鐵路高速發(fā)展及接觸線檢測(cè)和跟蹤研究中有較好的應(yīng)用價(jià)值。

5   結(jié)語(yǔ)

由于機(jī)車(chē)在行駛過(guò)程中需要頻繁換線從而保證機(jī)車(chē)高速安全運(yùn)行，本文提出了一種高速鐵路接觸線檢測(cè)與跟蹤方法，以解決在行車(chē)過(guò)程中弓網(wǎng)脫離而引發(fā)的安全事故，并克服了多個(gè)相似運(yùn)動(dòng)目標(biāo)同時(shí)運(yùn)動(dòng)復(fù)雜情況問(wèn)題。為了能夠?qū)崿F(xiàn)對(duì)紅外視頻中接觸線的準(zhǔn)確跟蹤，利用Hough 變換提取直線參數(shù)和卡爾曼濾波跟蹤參數(shù)，但尚不能做到自適應(yīng)檢測(cè)。Hough 變換檢測(cè)速度仍然有上升空間，另外卡爾曼濾波跟蹤參數(shù)部分，在接觸線運(yùn)動(dòng)特征變化不大的情況下跟蹤較為穩(wěn)定，但在目標(biāo)位置出現(xiàn)遮擋嚴(yán)重等背景復(fù)雜情況下跟蹤效果不佳，需要在后續(xù)工作中進(jìn)一步研究。

參考文獻(xiàn)：

[1] SINHA A,KIRUBARAJAN T,BAR-SHALOM Y.Application of the Kalman-levy Filter for Tracking Maneuvering Targets[J].IEEE Trans. on Aerospace and Electronic Systems,2007,43(3):1099-1107.

[2] 楊德江.高速鐵路接觸網(wǎng)檢測(cè)技術(shù)研究[J].建材與裝飾,2020(02):253-254.

[3] 陳高華,卓東風(fēng),夏錫瑞.基于Hough變換直線檢測(cè)中的峰值提取[J].太原科技大學(xué)學(xué)報(bào),2006(04):256-258.

[4] 何正友,程宏波.高速鐵路牽引供電系統(tǒng)健康管理及故障預(yù)警體系研究[J].電網(wǎng)技術(shù),2012,36(10):259-264.

[5] 閻玲玲,金立軍,張文豪,等.基于表面溫度法的接觸線磨耗紅外診斷[J].機(jī)電工程,2011,28(12):1426-1429.

[6] 熊磊.基于紅外熱像儀的接觸網(wǎng)載流溫度監(jiān)測(cè)管理系統(tǒng)的開(kāi)發(fā)[D].西安:西安理工大學(xué),2017.

[7] CAO J,CHEN L,WANG M,et al.Implementing a Parallel Image Edge Detection Algorithm Based on the Otsu-Canny Operator on the Hadoop Platform[J].Computational Intelligence ＆ Neuroscience,2018(3):1-12.

[8] 安建堯,李金新,孫雙平.基于Prewitt算子的紅外圖像邊緣檢測(cè)改進(jìn)算法[J].杭州電子科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2018,38(05):18-23+39.

[9] 吳銳,黃劍華,唐降龍,等.基于灰度直方圖和譜聚類(lèi)的文本圖像二值化方法[J].電子與信息學(xué)報(bào).2009,31(10):2460-2464.

[10] SATZODA R K,SATHYANARAYANA S,SRIKANTHAN T,et al.Hierarchical additive hough transform for lane detection[J].IEEE Embedded Systems Letters,2010,2(2):23-26.

[11] 王紹霖,付永生.Hough變換邊緣參數(shù)提取算法[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),1996(04):471-474.

[12] 劉春閣.基于Hough變換的直線提取與匹配[D].阜新:遼寧工程技術(shù)大學(xué),2009.

[13] 許龍.基于設(shè)備維修策略的接觸網(wǎng)管理信息系統(tǒng)的研究[D].成都:西南交通大學(xué),2013.

[14] 王紹霖,付永生.Hough變換邊緣參數(shù)提取算法[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),1996(04):471-474.

[15] 郭玲紅,李亞立.基于卡爾曼濾波的動(dòng)目標(biāo)運(yùn)動(dòng)參數(shù)跟蹤測(cè)量[J].自動(dòng)化技術(shù)與應(yīng)用,2012,31(08):20-23.

（本文來(lái)源于《電子產(chǎn)品世界》雜志2021年1月期）