郭雪松,王 勝
(中國電子科技集團公司第三研究所,北京 100015)
水下小目標(biāo)跟蹤技術(shù)是探測水下小目標(biāo)的關(guān)鍵技術(shù)之一[1]。本文在目標(biāo)檢測結(jié)果的基礎(chǔ)上(本文不討論目標(biāo)檢測結(jié)果的獲取過程),通過對聲吶圖像連續(xù)幀的處理,得到目標(biāo)運動軌跡,并根據(jù)目標(biāo)運動特征過濾掉干擾軌跡,最終獲得真實目標(biāo)的位置、方位等信息。
水下小目標(biāo)跟蹤技術(shù)可分為3 個步驟[2]:一是起始軌跡的建立(包含軌跡頭的確定);二是對聲吶圖像連續(xù)幀的處理,確認(rèn)軌跡并進行軌跡管理;三是軌跡撤銷/終止。
每個軌跡的第一個點稱為軌跡頭。軌跡生成前,以第一個點為中心,建立一個環(huán)形波門。此波門的范圍與目標(biāo)的最大速度、最小速度有關(guān)。
每一個軌跡頭與在波門內(nèi)的點關(guān)聯(lián)后,都要根據(jù)目標(biāo)的最大速度和最小速度推算出第三個點可能的范圍,若第三次關(guān)聯(lián)成功,則可以生成一條起始軌跡。一般假定目標(biāo)勻速直線運動得到起始軌跡。
采用如下方法進行軌跡的確認(rèn):
(1)若第四次關(guān)聯(lián)成功,則確認(rèn)是一條軌跡,放入軌跡庫進行管理;
(2)若第四次未關(guān)聯(lián)成功,則采用外推的方法,與假定的目標(biāo)點進行關(guān)聯(lián),第五次與真實檢測點關(guān)聯(lián)成功,則確認(rèn)是一條軌跡,放入軌跡庫進行管理;
(3)若第四次未關(guān)聯(lián)成功,采用外推的方法,與假定的目標(biāo)進行關(guān)聯(lián),第五次仍沒有真實檢測點進行關(guān)聯(lián),則清除起始軌跡點。
1.4.1 軌跡預(yù)測
軌跡預(yù)測的方法步驟如下。
(1)根據(jù)目標(biāo)的運動特征建立直角坐標(biāo)系下的運動模型[3]。目標(biāo)的模型反映在狀態(tài)矩陣上可表示為:
式中:A1表示加速運動模型,A2表示勻速運動模型,T表示運動時間。
(2)k+1 時刻預(yù)測模型概率μj(k+1|k),如式(3)所示:
式中:πij為狀態(tài)轉(zhuǎn)移矩陣。
交互權(quán)重計算如式(4)所示:
(3)確定k+1 時刻模型j對應(yīng)的濾波器估計算法的狀態(tài)輸入及其協(xié)方差:
式中:Xi(k)表示運動模型i對應(yīng)的濾波算法在時刻k的輸出,Pi(k)是其對應(yīng)的協(xié)方差。
(4)采用kalman 濾波對各個模型的運動狀態(tài)進行估計。
(5)利用濾波算法獲得各個模型的信息及其協(xié)方差S,計算各個模型的似然函數(shù)L。
(6)利用似然函數(shù)更新模型概率:
(7)根據(jù)更新后的模型概率μj(k+1),對各模型對應(yīng)的估計結(jié)果進行輸出交互,最終完成k+1 時刻的目標(biāo)狀態(tài)及其協(xié)方差的估計。
由于水下環(huán)境的限制及目標(biāo)在遠距離上的回波非常微弱,目標(biāo)檢測的概率會相應(yīng)地降低,此時目標(biāo)的軌跡就會不連續(xù)。多模型的方法對目標(biāo)位置的估計更準(zhǔn)確,更有利于目標(biāo)軌跡的管理。
1.4.2 軌跡擬合
目標(biāo)軌跡的測量點在空間上的離散程度很高,不易判斷目標(biāo)的真實運動方向。在對目標(biāo)進行估計和預(yù)測的同時,對其軌跡進行擬合,能夠反映目標(biāo)的真實運行方向。
(1)在直角坐標(biāo)系中,將坐標(biāo)系劃分為n×n的棋盤;
(2)對包含軌跡的方格求均值;
(3)如果方格中的點不連續(xù),則可以判定軌跡點在出方格后又進了方格,求兩個方格的均值;
(4)將所有的方格點按照時間順序排列。得到的新軌跡即為擬合后的軌跡。
1.4.3 檢測點關(guān)聯(lián)
采用最鄰近法進行關(guān)聯(lián)。具體為通過定義以下兩個度量來進行數(shù)據(jù)與目標(biāo)的關(guān)聯(lián)。
(1)基于歐氏距離的數(shù)據(jù)有效性度量。在波門中,只有小于門限值的測量點才有可能是目標(biāo),由于檢測門限設(shè)置得較寬松,可能會有多個點符合條件。如果沒有其他的度量方法,只能通過增加軌跡分支的方法,通過后續(xù)周期的測量與預(yù)測值的關(guān)聯(lián)來確定是否為真實軌跡,這會大大增加計算量。因此,可以采用目標(biāo)尺寸和反射強度的相似性度量來進一步提高關(guān)聯(lián)準(zhǔn)確度和關(guān)聯(lián)速度。
華堂村桃形李種植歷史悠久,但普遍種植、銷售還不到10年時間.村里在產(chǎn)業(yè)引導(dǎo)上也做了一些工作:一是有統(tǒng)一印制的禮盒包裝,在水果市場上打出桃形李品牌;二是建有一個桃形李市場,供農(nóng)戶和外地客商交易;三是組織過一些桃形李種植技術(shù)的培訓(xùn).另據(jù)村支書介紹,村里還與臺灣客商合作開發(fā)過桃形李果酒,但因為桃形李淀粉含量少沒有成功.
(2)基于目標(biāo)尺寸和反射強度信息的相似性度量。在有效性度量的基礎(chǔ)上,對符合條件的測量點進行相似性度量。相似性度量基于以下假設(shè):在空間域和時間序列中,真實目標(biāo)的尺寸(聲吶圖像中以像素個數(shù)表示)不會發(fā)生劇烈變化,即像素個數(shù)不會突然增加或減少很多,應(yīng)該有一個漸變過程。在反目標(biāo)探測上,基于尺寸的相似性度量需要實驗驗證,可能依據(jù)尺寸的變化不能起作用。在空間域和時間序列中,真實目標(biāo)的反射強度(聲吶圖像中以灰度級表示)不會發(fā)生突變。在相鄰的兩個周期中,目標(biāo)的亮度可能發(fā)生變化,變化的程度(或變化速度)應(yīng)該有一個上限和下限。設(shè)置一個門限,只要點在這個門限內(nèi),可以認(rèn)為其是目標(biāo)的可能性很大。
如果仍然有不止一個點符合條件,可以形成分支軌跡,然后依據(jù)軌跡撤銷規(guī)則對虛假軌跡進行撤銷。
針對以下幾種情況,需要分別進行處理。
(1)對只有軌跡頭的情況,只要其后的第一個周期內(nèi)沒有點跡出現(xiàn),就將其撤銷。
(2)對一條初始化的軌跡來說,如果在以后連續(xù)三個掃描周期內(nèi)沒有出現(xiàn)任何點跡(這里指沒有關(guān)聯(lián)上任何軌跡點),這條初始軌跡就可以在數(shù)據(jù)庫中被消去。
(3)對于一條已確認(rèn)的軌跡,可以設(shè)定連續(xù)4~6 個掃描周期內(nèi)沒有點跡落入相應(yīng)的關(guān)聯(lián)門內(nèi)作為軌跡撤銷準(zhǔn)則。此時,需要多次利用盲目外推的方法,擴大波門,對丟失目標(biāo)進行再捕獲。也可以采用小滑窗檢測器,設(shè)立一個軌跡結(jié)束準(zhǔn)則。只要滿足該準(zhǔn)則,就可以對該軌跡進行撤銷。令滑窗寬度為5,軌跡撤銷門限為4,這就意味著在連續(xù)5 個周期中,只要有4 個周期沒有點跡存在,就宣布該軌跡被撤銷。也可以采用連續(xù)3 個周期沒有點跡時宣布該軌跡被撤銷。對已經(jīng)確認(rèn)的軌跡進行撤銷時,可以考慮加大滑窗寬度,以便放寬撤銷條件。
(4)對于一條跟蹤了很長時間的、穩(wěn)定的軌跡所對應(yīng)的目標(biāo),在目標(biāo)離開探測器工作范圍時,也需要對其進行撤銷。
本文算法的整體流程如圖1 所示。
圖1 算法整體流程
利用Matlab 仿真所提出的跟蹤算法,結(jié)果如圖2 所示。
圖2 利用Matlab 仿真跟蹤算法
利用本文跟蹤算法的預(yù)測方法與傳統(tǒng)的kalman 預(yù)測方法[4-5]相比,兩者的跟蹤誤差如圖3 所示。由圖3 可以看出,本文跟蹤方法在不影響實時性的前提下,位置誤差相比傳統(tǒng)的kalman 預(yù)測方法有了很大的改善。
圖3 跟蹤誤差比較
針對多目標(biāo)跟蹤的仿真結(jié)果如圖4 所示。
圖4 對多目標(biāo)跟蹤算法的仿真結(jié)果
利用Matlab 對跟蹤算法進行仿真,輸出的軌跡如圖5 所示。
圖5 單目標(biāo)軌跡仿真
添加噪聲后的仿真軌跡如圖6 所示。噪聲為隨機在圖像中添加的稀疏噪聲點,輸出的軌跡如圖6 所示。
圖6 加噪聲的單目標(biāo)軌跡仿真
在單目標(biāo)跟蹤算法的基礎(chǔ)上,對多目標(biāo)跟蹤算法進行仿真,對多目標(biāo)跟蹤的結(jié)果如圖7 所示。
圖7 多目標(biāo)跟蹤仿真
水下小目標(biāo)的跟蹤分為軌跡起始、軌跡管理以及軌跡撤銷/終止3 大部分,其中的核心是數(shù)據(jù)關(guān)聯(lián)問題。由仿真試驗數(shù)據(jù)驗證,本文提出的跟蹤方法能夠?qū)崿F(xiàn)對水下小目標(biāo)的有效跟蹤。在此基礎(chǔ)上,可以繼續(xù)研究復(fù)雜背景下水下小目標(biāo)的穩(wěn)定跟蹤方法。