基于幀間差分累積的鐵路限界異物檢測提取算法

郭 碧,丁春平

(1.中鐵第四勘察設計院集團有限公司,武漢 430063； 2.浙江眾合科技股份有限公司,杭州 310051)

伴隨著軌道交通全自動運行系統技術的發展，對于無人值守列車前方障礙物探測的需求進一步加強，對列車行車前方限界內的障礙物檢測已成為亟待解決的課題[1-5]。當前，針對軌道交通限界環境下的障礙物檢測方法按感知方法可分成非接觸式和接觸式[6-7]。已投入運營的全自動運行線路北京燕房線即采用車底安裝的檢測橫梁通過與異物碰撞來實現對侵限異物的檢測與清除；非接觸式方法則包括激光、紅外、視覺識別等，其中視覺傳感器因其安裝簡便、信息豐富、結果直觀被廣泛采用[8-9]。

按照視頻設備和被監視場景的關系，利用車載視覺傳感器進行的異物檢測為動態背景下的前景目標檢測提取過程。目前，基于視覺傳感器技術的動態背景環境下前景對象的提取方法包括：一是根據空間坐標和成像坐標在運動中的對應關系，通過運動補償消減背景變化，再作背景減除提取前景目標。文獻[10]即以一維灰度投影結合高斯濾波來實現圖像快速去抖動，以統計分布方法完成背景更新，通過快速背景差分取得前景目標。文獻[11]基于高階奇異值分解和全變分思想，通過求解平穩背景、前景對象、動態背景之間的最小化約束方程，來實現對于運動目標的檢測。文獻[12]通過機器視覺理論和仿射原理構造空間模型，來實現對于障礙物入侵的判斷和測量。文獻[13]提出一種基于魯棒M估計全局像素點權值，通過Mean Shift聚類算法，實現不同運動點的分割動態背景下的目標檢測算法。二是將提取對象進行特征描述后，以重匹配方法實現目標檢測與跟蹤。如文獻[14]通過對檢測對象進行高斯金字塔分解建模后，來匹配動態背景中對象的位置。文獻[15]采用Random ferns和樸素貝葉斯方法，對目標對象進行特征描述匹配，通過背景補償引導檢測動態場景下的目標。文獻[16]以軌道線邊緣完整性及其特征來實現對于異物區域的檢測，通過提取異物特征以支持向量機對異物進行分類和辨識。三是以改進光流法等實現動態背景下目標的檢測方法[17]，考慮計算量及實時性，在此場景適用難度較大。鐵路限界內的障礙物檢測前景目標多樣，背景變化快，且夾雜由軌旁道砟、車輛振動、數據處理過程中的大量噪聲，加大了異物檢測提取的復雜度。

鑒于此，首先根據鐵路軌道交通線路環境變化緩和接近，針對動態背景環境，通過隔幀幀差方法來獲得多次差分結果，依托鐵路限界內路軌特征，采用數學形態學實現背景紋理重構來減除干擾及噪聲影響。隨后通過數值統計以側向差分灰度的累積投影值，來動態確定不同環境下的異物前景范圍，并通過最大類間方差法提取分割出前景目標。多幀隔幀差分相對于光流法具備實時性高的特點，相對于背景差分法能較好應對動態的背景條件，同時可提高對前景目標提取的實時性和魯棒性。

1 異物提取算法過程

本文所采用的軌道交通限界異物檢測提取方法的算法框圖如圖1所示。

圖1 異物檢測提取算法框圖

算法在初始化處理階段即提取軌道線性特征并以此為基礎確定前方軌道限界區域，隨后，依據列車前方限界區域范圍提取對應圖像數據作隔幀幀差處理，通過數學形態學來對鋼軌背景實現重構，來消除復雜背景中的道床枕木紋理特征。針對重建和濾除的差分結果，針對縱向和橫向作灰度累積投影變化策略來獲取差分突變位置作為前景目標區域。針對定位的圖像區域像素，以最大類間方差法來分割提取出前景目標，以方便作為模式識別的特征描述輸入參數。

2 初始化處理

2.1 圖像預處理

車載視覺傳感器采集的視頻序列幀中含有較多背景環境數據，同時因背景光照、行車振動、信號放大處理等摻雜干擾。本文以灰度直方圖均衡和中值濾波法來實現圖像數據平滑減輕噪聲干擾[18]。

2.2 軌道限界范圍標定

鐵路鋼軌因其材質和磨損，在背景圖像中呈現連續的線性特征，采用圖像邊緣提取算子可得到軌道線，考慮攝像機安裝位置與軌道相對固定，鐵路軌道曲率較小、軌面平坦，且軌道在圖像空間中呈現近景區域固定，遠景區域動態變化的特征，對圖像空間作近遠景劃分，通過分段曲線模型實現對于軌道線的檢測和跟蹤[19]。軌道限界范圍根據軌道線模型按圖像空間像素比例近似計算即可得到。根據鐵路限界要求，取最大邊距2 440 mm作為建筑限界外側距線路中心線距離，由鐵路標準軌距1 435 mm，可等比例得到對應的鐵路限界范圍，同時伴隨遠景區域軌道線跟蹤來動態更新限界范圍。根據空間坐標系中兩平行線在投影空間匯聚于一點，線路區域位于此消隱點(Pvanish(u，v))下方的原理[20-21]，本文以消隱點下二分之一位置vb處為分界點，標定出近景區和遠景區范圍。其中根據圖像數據編碼方式，文中圖像坐標系左上頂點為原點，以像素為單位，如圖2示。

圖2 軌道模型示意

3 隔幀幀間差分與背景補償

3.1 隔幀幀差法

將車載視頻傳感器采集的單幀圖像用一個二維函數I(u，v)表示，那么連續的序列幀可表示為I1，I2，I3，…，Ik…In。將序列中2組隔幀數為m的相鄰幀數據表示為：Ik-1(u，v)，Ik(u，v)和In-1(u，v)，In(u，v)，那么依據本文的思路，隔幀差分過程如下所述。

步驟1：數據參數及限界區初始化；取視頻序列幀數據Ik-1(u，v)，Ik(u，v)及In-1(u，v)和當前幀In(u，v)；

步驟2：差分Ik(u，v)與Ik-1(u，v)得Dk，k-1(u，v)；記錄并存儲為Dk，k-1(u，v)；將Ik(u，v)數據賦給變量Ik-1(u，v)；

步驟3：取In(u，v)與In-1(u，v)差分，記錄并存儲為Dn，n-1(u，v)，賦當前幀數據給變量In-1(u，v)；

步驟4：將兩次差分結果Dn，n-1(u，v)，Dk，k-1(u，v)以及當前Ik-1(u，v)，In-1(u，v)送入下一模塊，背景補償并作側向投影統計求差，以突變特征初定目標對象范圍；回到步驟2，重復此過程。

隔幀差分算法在視頻采集中同步進行，可應對異物出現在背景中無相對運動的場景，同時對丟幀具有適應性。本文取司機室內視覺傳感器采集的實驗場景為例，限界內前景目標為無源應答器。實驗結果如圖3所示。

圖3 兩組序列幀差分結果

對比幀組信息一致，結果中灰度數據信息少；當前幀差結果包含有因傳感器振動而出現的軌道輪廓；且差分結果中均存在干擾定位的枕木紋理等道床信息，需進一步進行補償處理。

3.2 基于數學形態學的背景補償

數學形態學開運算以結構元素構造的方法實現對于灰度圖像確定紋理的重構建。構造形態學結構元素b，以i(u，v)表示圖像I(u，v)的灰度值。形態學開運算過程為通過結構元素完成的先腐蝕后膨脹處理。其中灰度膨脹運算可表示為

(I⊕b)(u,v)=max{i(u-u′,v-v′)+

b(u′,v′)|(u′,v′)∈Db}

(1)

其過程表示結構元素b在I(u，v)定義域內灰度加運算遍歷過程中的最大值。為不更改原有像素點灰度值，本文取結構元素b的定義域Db內的所有坐標值為0，式(1)描述的算法過程為在Db對應領域范圍內確定最大值的過程。同樣地,灰度腐蝕運算可表示為

(IΘb)(u,v)=min{i(u+u′,v+v′)-

b(u′,v′)|(u′,v′)∈Db}

(2)

腐蝕運算則表示通過b與在I(u，v)定義域內取灰度最小值的過程。由上，灰度圖像開運算可表示為

I°b=(IΘb)⊕b

(3)

開運算通過結構元素腐蝕過程消除微細噪聲信息，再以膨脹重構出對象紋理結構。重構的精度受構造的結構元素的相似性影響。已知的紋理結構特征包括：橫向水平線性枕木邊緣結構、顯著軌道邊緣輪廓，以及其他相關設備結構。通過開運算可重構顯著紋理，來補償列車前方限界內枕木環境及抖動帶來的干擾。文中針對性構造結構元素b1，b2，b3，序列幀分辨率為521×480 pixel，其中：

b1：構造枕木邊緣，取30 pixel水平直線型結構元素；

b2：構造軌道結構輪廓，近景區軌道邊緣傾角由軌道線追蹤算法中的Hough變換提取得到，其中左右側斜率分別為kl，kr，設65 pixel線性結構元素。遠景區部分以縱向投影統計削減。

b3：本場景中車輛駕駛室外雨刷作特殊結構描述為60長25 pixel線性結構。

通過b1，b2，b3分別對當前幀In差分所得作開運算后，累加得到背景紋理結構的重構圖如圖4(a)所示。通過減除重構紋理實現補償的結果如圖4(b)所示，行車前方限界內道床、設備等因抖動產生的干擾得到削弱。

圖4 差分背景補償結果

4 基于側向灰度累積值的目標定位提取與分割

4.1 目標定位提取

車載視覺傳感器采集的視頻序列幀中包含的路軌道床信息會因光照和行車抖動而富含大量噪聲干擾，且不同場景下的信息具有較大的差異性，為應對不斷變化的環境條件，本文考慮采用統計思想，針對相鄰幀背景變化較小的特點，以限界區背景數據值的單方向投影來反映灰度信息在坐標中的分布，鄰幀的背景信息分布呈現較強的近似性。當有非路軌環境目標出現在限界范圍內時，單向投影統計值會呈現出部分范圍持續性的數據變化。文中通過單側統計差分結果的灰度和值對比得到異物疑似區域。

算法思路主要包括以下過程。

步驟1：參數初始化，求Dk，k-1(u，v)，Dn，n-1(u，v)遠景橫向投影累積值，記錄兩組結果相減的偏差向量ec[v]；

步驟2：以本幀橫向偏差向量ec[v]與前v-1項偏差均值Thr之差作為疑似區判定邊界。當ec>Thr，且持續有L個像素時，判別為疑似區域并記錄為[Gs，Gx]，若ec≤Thr，則返回繼續循環；

圖5 基于側向灰度累積值的目標定位算法示意

步驟3：在v軸[Gs，Gx]區域，作縱向差分灰度投影累積，并同步驟3方法完成u軸疑似區域[Kz，Ky]確定。實驗結果如圖5所示。因圖中限界內圖像數據差分后信號較為微弱，為方便說明此處以負片效果作為顯示，說明側向投影累積與前景目標的對應位置。

由目標出現而產生的與原枕木道床環境不同的灰度分布差異，可完成對于軌行區限界范圍內的前景目標定位，同時為方便目標分割及描述特征向量值的提取，在疑似區域外多擴3像素以加強背景參照信息。

4.2 目標分割與二值化

為區分目標和背景分界，利于后期對于異物對象特征向量值的描述和提取，針對疑似范圍內的圖像作分割運算以確定前景對象?？紤]基于最大類間方差的動態閾值分割方法相對于聚類分析、區域生長、主動輪廓等方法在速度、魯棒性、算法復雜度上的優勢，用最大類間方差法實現目標分割。最大類間方差通過統計不同灰度級別的分布的概率，以使方差最大來確定最佳閾值，采用的數學描述可表示為

Ts=Argmax[w1(t)(J1-J0)2+w2(t)(J2-J0)2]

(4)

式中，Ts表示灰度閾值；w1(t)和w2(t)分別表示前景目標和背景信息在圖像像素點中的比重；J0(t)為平均灰度值，J1(t)和J2(t)分別表示前景目標和背景的灰度均值；取限界內疑似范圍外方灰度值作背景劃分參量，本文舉例的實驗提取結果如圖6所示。

圖6 前景目標提取及分割結果

5 前景目標定位提取算法流程

本文基于幀間差分累積的鐵路限界異物檢測提取算法的主要過程包括初始化、隔幀幀差、背景補償、側向灰度累積定位與分割幾個主要部分，算法流程如圖7所示。算法自讀取圖像序列幀時，即開始循環做差分檢測，以避免漏檢車輛緩動場景下的異物目標，完成定位提取目標后，將目標及相應幀送入后續模式識別等模塊，作為危險辨識及預警輸入。

圖7 鐵路限界異物定位提取算法流程

6 實驗分析與驗證

本文通過室外場景實驗、TrainSimulator環境仿真、現場行車采集，共構造了47個不同的場景，用于算法場景測試，并將測試序列幀統一為512×480像素尺寸，每場景共包括45幀圖像。算法處理過程通過Matlab2013實現編寫和測試。

本次以現場行車采集的視頻序列幀場景為例，采用傳統幀間差分法作為對比說明。差分處理結果如圖8所示。

圖8 與幀間差分法結果對比

幀間差分法得到的差分結果如圖8(b)所示，圖中可見大量的背景干擾信息，較難區分出限界內前景目標，且差分結果中出現了目標的重影現象。本文算法結果如圖8(c)所示，通過對背景進行形態學重構補償，有效地消減了干擾信息。因幀間差分法無目標定位提取算法過程，本文后續通過對測試場景進行實驗，來驗證算法性能。

本文算法處理場景數據得到的疑似物目標定位提取結果(部分)如圖9所示。其中，圖9(a)表示光線充足下列車行車前軌行區域有未撤離人員環境場景；圖9(b)表示較低光線條件下限界內有山體落石出現場景；圖9(c)表示正常光線條件下，在道床上有懸空或跌落對象的場景。本文算法對異物目標的檢測提取結果如圖9(d)、圖9(e)和圖9(f)所示，可見前景異物目標于限界區域成像部分均實現提取和分割。

圖9 實驗場景定位及提取部分結果

針對構造實驗場景檢測統計實驗情況如表1所示。于行車前方限界內有異物的場景中，異物目標檢測率為96.87%，而于行車前方限界內無異物的場景出現了較高的誤檢率，誤檢率為26.67%。經實驗，算法具有較好的檢測效率，同時定位提取過程的平均耗時為0.137s。

表1 本文算法檢測數據統計

分析實驗過程及場景發現，造成誤檢的主要對象為路軌道床中間的安裝設備以及光照下的陰影和光斑等對象，此類不影響安全行車的非障礙物目標可在后續建立特征庫通過學習算法在目標識別過程中去除。

7 結論

對于車載視覺傳感器下的軌道交通線路動態背景環境，通過改進的隔幀差分法能較為有效地實現對于前景目標的提取。算法通過隔幀差分以獲得目標對照信息，并以鐵路限界背景紋理先驗知識，通過重構形態學背景紋理在一定程度上有效補償了限界范圍內的相似背景信息，且具有較好的實時性；算法通過側向差分灰度的累積投影值來動態確定不同環境下的異物對象范圍，可檢測隨機性對象，同時可較好應對車載視覺傳感器因抖動和數據處理中累積的干擾及噪聲。通過構建不同環境下的實驗，算法可完成對于動態路軌環境下的異物目標定位提取，檢測的實時性和準確率較好。后續可通過提取本算法檢測獲得的前景目標描述特征，通過模式識別算法實現對于異物目標的分類識別，輔助行車決策控制。