黃超
摘 要: 針對目前機場區(qū)域視覺監(jiān)視的需要,在對機場區(qū)域采集的視頻圖像進(jìn)行預(yù)處理的基礎(chǔ)上,利用背景差分法、幀間差分法對機場敏感區(qū)域靜態(tài)背景下的運動目標(biāo)進(jìn)行檢測;以昌北國際機場和哈爾濱太平國際機場區(qū)域為例,采用Kalman濾波跟蹤算法和Meanshift均值漂移算法檢測區(qū)域的飛機及引導(dǎo)車等移動目標(biāo)。通過Matlab仿真,分析并對比不同算法的適用情況并驗證其可行性。
關(guān)鍵詞: 機場區(qū)域; 圖像監(jiān)控; 目標(biāo)檢測; 濾波跟蹤
中圖分類號: TN911.23?34 文獻(xiàn)標(biāo)識碼: A 文章編號: 1004?373X(2017)14?0092?03
Abstract: For the demand of visual surveillance in airport area, the background difference method and inter?frame difference method are used to detect the moving target in static background of the airport sensitive area, which is on the basis of the preprocessing for the video image collected in the airport area. Taking the areas of Nanchang Changbei International Airport and Harbin Taiping International Airport as the examples, the Kalman filtering tracking algorithm and Meanshift algorithm are used to detect the moving targets in the area, such as the airplane and guided vehicle. With Matlab simulation, the application situation of different algorithms is analyzed and compared to verify their feasibilities.
Keywords: airport area; image monitoring; target detection; filtering tracking
0 引 言
隨著機場流量的日益增長,機場不同區(qū)域布局復(fù)雜性的逐漸提高,對機場區(qū)域監(jiān)視的需求也越來越大。機場區(qū)域的特點是,監(jiān)視環(huán)境是一個有限的、大尺度的空間。場面移動目標(biāo)繁多,有飛機、車輛以及機場工作人員等,場面為室外環(huán)境,易受天氣、光照影響。目前應(yīng)用于機場區(qū)域目標(biāo)監(jiān)視與跟蹤的手段主要有三種,場面監(jiān)視雷達(dá)(SMR)、ADS?B和多點定位(Multi?lateration)系統(tǒng)[1?2]。針對機場區(qū)域的上述特點,場面雷達(dá)獲取目標(biāo)距離、方位精度高,但作用距離短,造價昂貴;ADS?B為合作式監(jiān)視手段,場面移動目標(biāo)需安裝設(shè)備,實現(xiàn)難度大;多點定位系統(tǒng)需要地面基站與飛機應(yīng)答機形成小型收發(fā)網(wǎng)絡(luò),目前技術(shù)成熟度不夠。智能監(jiān)控系統(tǒng)可以利用計算機圖像識別技術(shù)實現(xiàn)在機場區(qū)域上人的視覺功能,同時降低成本,可實現(xiàn)度高[3?4]。本文采用一種對機場區(qū)域進(jìn)行監(jiān)控的新方法,利用計算機視覺技術(shù)對場面目標(biāo)進(jìn)行檢測和跟蹤。
1 采集機場圖像的噪聲抑制
在機場區(qū)域運用攝像機采集到的圖像受到的噪聲會直接影響到檢測目標(biāo)的準(zhǔn)確性,因為較大的噪聲導(dǎo)致目標(biāo)的丟失或者偽目標(biāo)的出現(xiàn)[5]。因此在對圖像進(jìn)行處理前需要先去噪。均值濾波也叫作線性濾波,該方法用像素灰度平均的方法平滑圖像,去除噪聲。
首先選取一個移動窗口,確定該窗口的中心,對窗口內(nèi)鄰域的像素點的灰度平均值進(jìn)行計算,計算出的值作為該點的新的灰度值。設(shè)有一幅圖像,是通過均值濾波方法得到的新圖像,除去中心點的鄰域的像素點共有S個,可得到:
在通過均值濾波對圖像進(jìn)行處理時,選取的窗口大小也會對處理產(chǎn)生影響。常用的有窗口和窗口,窗口的圖像比窗口的圖像更加清晰,但是去噪性能較差。均值濾波算法實現(xiàn)簡單,僅需要對窗口內(nèi)的像素灰度進(jìn)行平均計算。但由于各像素均采用了經(jīng)過計算的平均值,所以有些細(xì)小的特征就被忽略了。因此均值濾波算法可能會使圖像的輪廓變得圓滑,細(xì)小特征無法被識別。中值濾波能夠有效減弱圖像細(xì)節(jié)模糊,利用它對圖像掃描中產(chǎn)生的噪聲消除非常有效[6]。中值濾波仍選擇移動窗口并確定其中心,將該點鄰域的所有像素的灰度值按順序排列出來,針對窗口選取奇偶的不同,選取其中這些像素的中點,選取的值即為中心點的灰度值。設(shè)為的圖像,圖像為經(jīng)過中值濾波后的圖像, 除去中心點的鄰域的像素點共有S個。將圖像中各像素點按灰度值大小排成一列,如下:
2 背景差分法和幀間差分法融合的Kim算法
為了能夠?qū)煞N運動目標(biāo)檢測方法結(jié)合在一起,避免各自影響檢測目標(biāo)準(zhǔn)確性的因素,更好地應(yīng)用到機場區(qū)域運動目標(biāo)檢測中,使用與或運算將背景差分二值圖像和幀間差分二值圖像聯(lián)系起來,提出以下算法模型:
用表示幀差圖像,代表背景差分圖像。這兩幅差分圖像要在二值化之后進(jìn)行相與運算。二值化選取公式如下:
式中,和分別為幀間差分和背景差分圖像二值化的閾值。二值化圖像經(jīng)過相與運算后可以得到最終的二值運動模板,模板表示如下:
在Kim方法中,利用差分圖像進(jìn)行檢測時,既采用了連續(xù)的幀差圖像,也利用了當(dāng)前幀與事先存儲的背景幀的差分圖像。這樣既可以避免由于圖像重疊導(dǎo)致的目標(biāo)點檢測不完整,也防止了對背景更新的依賴性和對外界環(huán)境的敏感性。所以該方法可以比較準(zhǔn)確地檢測出運動目標(biāo),并且實時性強。
3 機場區(qū)域運動目標(biāo)跟蹤方法設(shè)計
3.1 Kalman濾波跟蹤方法設(shè)計
Kalman濾波器的本質(zhì)是一種遞歸的估計算法,該算法根據(jù)物體前一時刻的狀態(tài)和觀測到的狀態(tài)來反推物體當(dāng)前狀態(tài)[6?7]。Kalman濾波算法應(yīng)用于目標(biāo)跟蹤的主要過程是:首先獲取前幀的目標(biāo)跟蹤結(jié)果;利用該結(jié)果預(yù)測當(dāng)前幀的目標(biāo)位置;在該位置附近搜索目標(biāo)。
若目標(biāo)存在于搜索區(qū)域內(nèi),則繼續(xù)處理下一幀圖像;若目標(biāo)不存在于搜索區(qū)域內(nèi),則擴大搜索區(qū)域,在當(dāng)前幀重新搜索目標(biāo)。在本文中,設(shè)運動目標(biāo)的初始位置坐標(biāo)為, 同時設(shè)運動目標(biāo)在橫縱坐標(biāo)軸上的速度無限小并接近于0,則初始化的狀態(tài)向量為:
系統(tǒng)噪聲與觀測噪聲的協(xié)方差矩陣分別為四維和二維的單位矩陣。觀測狀態(tài)向量為,其中為物體在k時刻的橫縱位置坐標(biāo)。利用Kalman五個更新方程完成目標(biāo)位置的迭代搜索。
3.2 Meanshift跟蹤算法設(shè)計
Meanshift 算法也叫均值漂移算法,該算法應(yīng)用于目標(biāo)跟蹤的主要過程是:
(1) 獲取含有目標(biāo)的單幀圖像,在確定目標(biāo)區(qū)域后,建立目標(biāo)模板,確定運動目標(biāo)特征,且,初始位置為。
(2) 建立候選目標(biāo)模板,確定候選運動目標(biāo)特征,將目標(biāo)在前一幀的初始位置代入,確定一個候選目標(biāo)的初始特征。
(3) 計算Bhattacharyya系數(shù),所得的值對兩個模板的相似性進(jìn)行判斷。
(4) 為了保證不同像素點對確定目標(biāo)的貢獻(xiàn)不同,計算加權(quán)值。
(5) 利用Meanshift向量及迭代過程對目標(biāo)的新位置進(jìn)行計算并更新。
(6) 更新目標(biāo)位置后,再計算新目標(biāo)特征和目標(biāo)特征的相似性。
(7) 歐幾里得距離最大化也就是Bhattacharyya系數(shù)的最小化,因此對在和位置處兩個相似性系數(shù)進(jìn)行大小判斷。如果新位置處的相似性大于原位置處的相似性,那么就繼續(xù)循環(huán)尋找,直到找到某位置使得:
(8) 設(shè)定一個非負(fù)的容許參數(shù),對新位置與舊位置的坐標(biāo)距離進(jìn)行限定,因為幀間時間間隔較短,目標(biāo)移動的速度是有限的,如果超出該容許參數(shù)則認(rèn)為位置定義不準(zhǔn)確,需要重新進(jìn)行迭代計算。
(9) 如果新舊位置距離恰當(dāng),需要繼續(xù)跟蹤目標(biāo),則重復(fù)步驟(1)~(8),如果不需要進(jìn)行跟蹤,此時輸出結(jié)果,算法結(jié)束。
4 仿真實驗跟蹤對比
為了實現(xiàn)機場區(qū)域滑行飛機的Kalman濾波算法與Meanshift算法的仿真跟蹤對比,以下選取了在哈爾濱太平國際機場區(qū)域滑行的飛機進(jìn)行跟蹤。在視頻中截取第20幀、第25幀、第30幀和第40幀圖片,然后分別利用兩種跟蹤算法進(jìn)行仿真實驗。如圖1~圖4所示。
第20幀機場區(qū)域飛機行駛仿真結(jié)果:目標(biāo)真實位置為(796,235), Kalman濾波結(jié)果為(803,227), Meanshift算法目標(biāo)像素位置為(871,223)。
第25幀機場區(qū)域飛機行駛仿真結(jié)果:目標(biāo)真實位置為(779,226), Kalman濾波結(jié)果為(826,238), Meanshift算法目標(biāo)像素位置為(790,241)。
第30幀機場區(qū)域飛機行駛仿真結(jié)果:目標(biāo)真實位置為(751,214), Kalman濾波結(jié)果為(535,226),Meanshift算法目標(biāo)像素位置為(784,217)。
第40幀機場區(qū)域飛機行駛仿真結(jié)果:目標(biāo)真實位置為(706,205), Kalman濾波結(jié)果為(380,253),Meanshift算法目標(biāo)像素位置為(710,214)。
在對飛機跟蹤過程中,初始時Kalman濾波的跟蹤結(jié)果較好,但是隨著時間的遷移,跟蹤效果越來越差,最后在干擾下,甚至無法預(yù)測到飛機的正確位置。但是對于Meanshift算法,由于進(jìn)行配對的是候選目標(biāo)模板和真實目標(biāo)模板,所以對于部分遮擋的情況(如本例,指揮員進(jìn)行指揮),只要能夠搜索到部分候選目標(biāo)的模板并將其與真實目標(biāo)模板配對成功,這樣就可以跟蹤到目標(biāo)。通過上文的仿真實驗驗證可以看出,Kalman濾波跟蹤算法適用于天氣情況較差,無法對識別的模板做出正確判斷的情況下,由于該方法計算量小,所以實時性強,但隨著時間的推移,預(yù)測結(jié)果可能會出現(xiàn)較大偏差;而Meanshift 算法是一種迭代循環(huán)算法,它適用于部分遮擋的情況,只要能夠識別出目標(biāo)模板的一部分,那么就可以正確地追蹤到目標(biāo)。
5 結(jié) 論
在機場區(qū)域監(jiān)視領(lǐng)域中,少數(shù)的雷達(dá)及合作式設(shè)備已投入使用,但昂貴的造價及維修費用使得在中小機場無法大面積推廣[8?9]。采用計算機圖像識別技術(shù)可以有效降低成本,同時經(jīng)仿真驗證,在機場區(qū)域大尺度環(huán)境中,利用本文的檢測及跟蹤算法可以有效地識別并跟蹤機場區(qū)域移動目標(biāo),具有一定的現(xiàn)實意義和研究價值。
未來在對目標(biāo)的檢測和跟蹤可以從以下兩個方向進(jìn)一步完善:人工神經(jīng)網(wǎng)絡(luò)技術(shù),這種網(wǎng)絡(luò)系統(tǒng)具有很強的學(xué)習(xí)性和自適應(yīng)性,未來,多目標(biāo)跟蹤識別應(yīng)用非均勻網(wǎng)絡(luò)(接口技術(shù))將是一種發(fā)展趨勢;人機交互合作式技術(shù),計算機雖然具有運行準(zhǔn)確,處理數(shù)據(jù)速度快等優(yōu)點,但是它無法對各種隨機事件的發(fā)生做出正確的應(yīng)對,因此單單靠計算機進(jìn)行識別難度非常大[10]。所以在一定情況下,人的指導(dǎo)和參與是必不可少的。并且人的加入可以使得用戶與系統(tǒng)之間的聯(lián)系更加緊密,溝通更加快速、有效[11]。
參考文獻(xiàn)
[1] 顧春平.空中交通管制監(jiān)視新技術(shù)簡介[J].現(xiàn)代雷達(dá),2010,32(9):1?5.
[2] SAGE K, YOUNG S. Security applications of computer vision [J]. IEEE aerospace & electronic systems magazine, 1999, 22(8): 1090?1094.
[3] 張銘銘.視頻監(jiān)控新系統(tǒng)中圖像增強技術(shù)的研究[D].濟(jì)南:濟(jì)南大學(xué),2009.
[4] PIZURICA A, PHILIPS W W, LEMANHIEU I, et al. A versatile wavelet domain noise filtration technique for medical imaging [J]. IEEE transactions on medical imaging, 2003, 22(3): 323?331.
[5] 肖大雪.淺析數(shù)學(xué)形態(tài)學(xué)在圖像處理中的應(yīng)用[J].科技廣場,2013(5):10?19.
[6] 朱明旱,羅大庸.基于幀間差分背景模型的運動物體檢測與跟蹤[J].計算機測量與控制,2006,14(8):1004?1006.
[7] 林佳乙,于哲舟,張健,等.基于背景差分法和幀間差分法的視頻運動檢測[J].儀器儀表學(xué)報,2008,29(4):111?114.
[8] 崔錦泰,陳關(guān)榮.卡爾曼濾波及其實時應(yīng)用[M].北京:清華大學(xué)出版社,2013:10?14.
[9] 周尚波,李昆.核函數(shù)帶寬自適應(yīng)的均值漂移目標(biāo)跟蹤算法[J].計算機工程與技術(shù),2012,33(7):2776?2781.
[10] 胡波,陳懇,徐建瑜,等.基于Kalman預(yù)測和Meanshift算法的視頻目標(biāo)跟蹤[J].光電子·激光,2009,20(11):1517?1522.
[11] 馬翠霞,任磊,滕東興,等.云制造環(huán)境下的普適人機交互技術(shù)[J].計算機集成制造系統(tǒng),2011,17(3):504?510.