毛藝帆，　張多林，　王　路

(空軍工程大學 防空反導學院， 陜西 西安 710051)

?

基于三角形拓撲相似性的航跡關(guān)聯(lián)算法

毛藝帆，張多林，王路

(空軍工程大學 防空反導學院， 陜西 西安 710051)

針對多個傳感器量測航跡在系統(tǒng)誤差下關(guān)聯(lián)融合困難的問題，提出基于拓撲4交差模型和三角形拓撲相似度的航跡關(guān)聯(lián)方法，將傳感器探測目標分成若干個拓撲三角形結(jié)構(gòu)，使用4交差模型排除不同傳感器探測中相關(guān)可能性小的三角形；然后定義三角形相似度，對歸屬于不同探測源的拓撲三角形相似程度進行度量，從而得出目標之間的相關(guān)程度，獲取航跡關(guān)聯(lián)結(jié)果。仿真結(jié)果表明，在目標密集編隊時，該算法具有有效性。

航跡融合；航跡關(guān)聯(lián)；拓撲4交差模型；三角形相似度

航跡關(guān)聯(lián)是航跡融合中的關(guān)鍵問題，解決了不同航跡是否屬于同一目標的問題，是后續(xù)信息能夠進行可靠融合的基礎(chǔ)，是跟蹤成功的關(guān)鍵。在分布式信息處理系統(tǒng)中，信息處理子節(jié)點提供目標的狀態(tài)估計航跡；在信息處理中心，對不同節(jié)點的航跡信息關(guān)聯(lián)對準，對關(guān)聯(lián)成功的航跡對再進行航跡融合，得到最終的融合結(jié)果?，F(xiàn)階段，已有大量的關(guān)于航跡融合及航跡關(guān)聯(lián)的研究成果，主要分為基于統(tǒng)計的方法[1-2]、基于模糊理論的方法[3-6]和基于灰色理論的方法[7-8]。這些方法大多只考慮探測的隨機誤差，未考慮系統(tǒng)誤差，且在目標密集編隊時，無法在短時間內(nèi)獲取良好的關(guān)聯(lián)效果。因此，需要在現(xiàn)有研究的基礎(chǔ)上，發(fā)掘新的方法，解決上述問題。

本文首先推導存在系統(tǒng)誤差的目標粗相關(guān)波門，再將航跡目標進行分割，劃分成最簡單的拓撲目標——三角形目標進行研究，并通過基于空間拓撲模型進行篩選，通過篩選的目標進行基于三角形相似度比較的精相關(guān)，最后結(jié)合雙門限準則和多義性處理，得到最終的航跡關(guān)聯(lián)結(jié)果。

1　目標粗相關(guān)波門計算

目標粗相關(guān)的思想是根據(jù)關(guān)聯(lián)波門，對待關(guān)聯(lián)的目標進行篩選，以便在進一步的關(guān)聯(lián)操作中，摒棄關(guān)聯(lián)波門以外的目標。粗關(guān)聯(lián)可有效減小計算量，提高計算效率[9]。在粗關(guān)聯(lián)中，關(guān)聯(lián)波門的選擇是尤為重要的，因此，本文對存在系統(tǒng)誤差時的關(guān)聯(lián)波門計算方法進行闡述。

假設(shè)有傳感器S1、S2對規(guī)定區(qū)域進行目標探測和跟蹤，探測到的目標分別為m和n個；探測到的航跡信息傳遞到信息融合中心，以融合中心為坐標原點，建立笛卡爾坐標系，取其中兩維進行研究，設(shè)S1和S2在該坐標系下的坐標為(xs1,ys1)和(xs2,ys2)；將傳感器S1所在位置設(shè)為融合中心，則其與傳感器及任一目標位置關(guān)系如圖1所示。圖中，S1和S2分別表示2個傳感器，Target表示目標真實位置，TS1和TS2分別表示傳感器的測量值?？紤]傳感器的系統(tǒng)誤差，其測距以及測角系統(tǒng)誤差為(Δrs1,Δθs1)和(Δrs2,Δθs2)。

圖1　融合中心、傳感器及任一目標位置關(guān)系示意圖

(1)

(2)

(3)

(4)

(5)

(6)

2　基于三角形相似性的航跡相關(guān)

2.1目標三角形分割

由于三角形是形狀比較簡單的拓撲結(jié)構(gòu)，可以很好地反映局部的拓撲結(jié)構(gòu)，而且對大規(guī)模的運算而言，運算負擔不會很重，其在模式識別方面已經(jīng)得到廣泛應(yīng)用[10]，因此，本文選擇將三角形結(jié)構(gòu)作為拓撲結(jié)構(gòu)元，研究基于目標之間拓撲信息的航跡關(guān)聯(lián)方法。

a) 傳感器1中選中的三角形示意圖

b) 傳感器2中目標點跡三角形組合

a) 傳感器2中目標點跡三角形組合

b) 傳感器2中目標點跡三角形組合

2.2目標三角形篩選

(7)

式中，Aο表示集合A的內(nèi)部；Bο表示集合B的內(nèi)部；?A表示集合A的邊界；?B表示集合B的邊界。該模型可描述多種點/點、點/線、點/面、線/線、線/面、面/面關(guān)系。

相對于其他模型(如9交模型等)，4交模型只需計算4個元素，且可表征比較完備的空間目標拓撲關(guān)系，其應(yīng)用比較廣；然而要進行4個交運算，運算復(fù)雜性比較高。因此，鄧敏等[11]46-49提出了一種混合模型——4交差模型，如式(8)所示。較之4交模型，其只需進行2個交運算和2個差運算，顯然，差運算比交運算的計算復(fù)雜性低；而且，該模型同樣可以很好地表征空間目標的拓撲關(guān)系。

(8)

若2個拓撲目標相關(guān)，則其在空間中的距離不應(yīng)相距太遠，而彼此之間的距離衡量需要定義另外的函數(shù)。假設(shè)有一個映射函數(shù)σ(Δ),其定義如式(9)。進而，定義距離函數(shù)，如式(10)所示。

(9)

(10)

例如目標相等和目標相離之間拓撲距離計算過程為

已知

可以看出，目標相等和目標相交、包含、包含于、覆蓋、覆蓋于這幾種關(guān)系的距離較近，由于傳感器探測到的目標都不是百分百準確的目標信息，所以具有不同程度的模糊性。因此本文認定，當空間目標處于相等、相交、包含、包含于、覆蓋、覆蓋于這幾種拓撲關(guān)系時，目標是有相關(guān)的可能性的，只是可能性大小不同；反之，當處于相接和相離關(guān)系時，認定其是不相關(guān)的。根據(jù)此方法，排除相關(guān)可能性較低的拓撲三角形目標，使用Arcview軟件即可實現(xiàn)，并可與Matlab進行數(shù)據(jù)互通。

2.3三角形相似度度量

已經(jīng)經(jīng)過相關(guān)篩選的三角形目標，需要進行進一步的處理，才能獲取目標關(guān)聯(lián)結(jié)果。本文根據(jù)目標之間的拓撲相似性來判斷目標相關(guān)程度。具體來說，當拓撲目標之間的結(jié)構(gòu)相似度越大，則可知構(gòu)成三角形拓撲結(jié)構(gòu)的目標點之間的關(guān)聯(lián)程度就越高，繼而問題轉(zhuǎn)化為如何判斷三角形目標相似度。參考文獻[14]中對三角形相似度的闡述，定義三角形拓撲相似度為

(11)

2.4雙門限準則和多義性處理

某一時刻的目標關(guān)聯(lián)結(jié)果并不能作為整體的航跡關(guān)聯(lián)結(jié)果，需要經(jīng)過多次的關(guān)聯(lián)檢驗，才有可能形成固定的航跡關(guān)聯(lián)對，因此，需要引入雙門限準則和多義性處理幫助進行航跡關(guān)聯(lián)判定。

雙門限準則[9,15]：選擇整數(shù)I和R(一般根據(jù)經(jīng)驗選值)，?k=1,…,R時刻，若關(guān)聯(lián)成功，則計數(shù)器mij(k)=mij(k-1)+1，mij(0)置0；若經(jīng)過R次關(guān)聯(lián)檢測，計數(shù)器值大于門限值I,即mij(R)≥I,則判定目標i和目標j為同一目標，否則，需進入下一檢驗周期進行關(guān)聯(lián)檢測。已經(jīng)判定為關(guān)聯(lián)的目標對則不需要進行進一步的關(guān)聯(lián)檢測。

3　算法流程

5) 對目標進行基于I/R規(guī)則的航跡相關(guān)檢驗，并在必要的時候進行多義性處理，得到最終航跡關(guān)聯(lián)對。對已關(guān)聯(lián)航跡對，使用基本加權(quán)融合法，得到航跡融合結(jié)果。

4　仿真分析

為說明本文方法的有效性，設(shè)置如下的基本仿真環(huán)境：假設(shè)在以融合中心為原點的坐標系中有2個傳感器，坐標分別為(0,0)和(0,50)。在2個傳感器的重疊觀測區(qū)域內(nèi)，有30批目標。傳感器的重疊觀測區(qū)域一般是不規(guī)則圖形，為了方便對關(guān)聯(lián)方法進行討論，設(shè)置該區(qū)域為矩形，目標初始位置隨機分布在30km×30km的矩形區(qū)域內(nèi)，目標的初速和初始方位角分別位于20～200 m/s和0～2π rad范圍內(nèi)，目標進行勻速直線運動。隨機測量距離誤差為50 m，隨機測量角度誤差為0.5°。假設(shè)共同觀測概率為1，進行50次蒙特卡羅仿真，每次20步，使用本文的方法和傳統(tǒng)的最近鄰(Nearest Neighbor,NN)算法進行仿真。為對航跡關(guān)聯(lián)的效果進行度量，定義航跡正確關(guān)聯(lián)概率Pr，Pr=Cr/(Cr+Cw),其中Cr表示正確關(guān)聯(lián)航跡對數(shù)目，Cw表示錯誤關(guān)聯(lián)航跡對數(shù)目。I/R邏輯選擇3/5規(guī)則。設(shè)置不同的仿真條件，進行仿真實驗。

1) 系統(tǒng)誤差范圍固定。假設(shè)傳感器1的系統(tǒng)測量距離和角度偏差最大值分別為100 m和1°，傳感器2的系統(tǒng)測量距離和角度偏差最大值分別為100 m和1°。使用本文方法和基本NN關(guān)聯(lián)法以及文獻[5]1694-1697中所提的模糊關(guān)聯(lián)法進行仿真實驗。在起始時刻，目標真實位置及2個傳感器的測量值如圖4所示，使用本文的方法，其最終航跡關(guān)聯(lián)結(jié)果的局部如圖5所示。

圖4　初始時刻目標示意圖

圖5　航跡關(guān)聯(lián)結(jié)果局部圖

根據(jù)本文算法和NN算法，得到正確關(guān)聯(lián)概率如圖6。在經(jīng)過20步的關(guān)聯(lián)檢驗之后，3種方法的關(guān)聯(lián)正確率均可達到0.9以上，說明在檢驗步數(shù)一定的情況下，3種方法均有效。

2) 系統(tǒng)測角誤差不固定。假設(shè)傳感器1的系統(tǒng)測量距離和角度偏差最大值分別為100m和1°，傳感器2的系統(tǒng)測量距離最大值為100 m，其系統(tǒng)測量角度偏差最大值為1°～5°，分別為1°、1.5°、2°、2.5°、3°、3.5°、4°、4.5°和5°。使用本文方法和基本NN關(guān)聯(lián)法以及文獻[5]1694-1697中的模糊關(guān)聯(lián)法進行仿真實驗，得到正確關(guān)聯(lián)概率如圖7。由圖可得，隨著測角誤差增大，NN算法性能下降嚴重，文獻[5]1694-1697的方法也有所下降，且在誤差達到5°時關(guān)聯(lián)正確率已經(jīng)降至0.7，本文方法性能比較穩(wěn)定，關(guān)聯(lián)正確率仍舊在0.8以上。

3) 系統(tǒng)測距誤差不固定。假設(shè)傳感器1的系統(tǒng)測量距離和角度偏差最大值分別為100 m和1°，傳感器2的系統(tǒng)測量角度偏差最大值1°，其系統(tǒng)測量距離偏差最大值為100～500 m，分別為100 m、150 m、200 m、250 m、300 m、350 m、400 m、450 m和500 m。使用本文方法和基本NN關(guān)聯(lián)法以及文獻[5]1694-1697中的模糊關(guān)聯(lián)法進行仿真實驗，得到正確關(guān)聯(lián)概率如圖8。隨著測距誤差增大，NN算法和文獻[5]1694-1697的方法性能均下降嚴重，說明以上算法對距離誤差是比較敏感的，本文算法性能也有所下降，但正確關(guān)聯(lián)概率穩(wěn)定在0.8以上。

圖6　系統(tǒng)誤差固定時的關(guān)聯(lián)正確率

圖7　系統(tǒng)測角誤差變化時的關(guān)聯(lián)正確率

圖8　系統(tǒng)測距誤差變化時的關(guān)聯(lián)正確率

由以上不同條件的仿真結(jié)果可知，首先，在系統(tǒng)誤差相差比較大的情況下，能夠獲取較好的關(guān)聯(lián)效果；進而，在系統(tǒng)誤差產(chǎn)生變化時，算法性能穩(wěn)定，在系統(tǒng)測量誤差逐漸增大時，算法的正確關(guān)聯(lián)概率雖然也會逐漸減小，但比較平穩(wěn)，變化不大，在可接收的范圍內(nèi)。

5　結(jié) 束 語

針對現(xiàn)有航跡關(guān)聯(lián)算法對編隊飛行及密集目標關(guān)聯(lián)效果較差的問題，提出基于三角形拓撲相似性的航跡關(guān)聯(lián)方法，該方法有別于傳統(tǒng)的方法，是解決航跡關(guān)聯(lián)問題的一個有益的嘗試，仿真結(jié)果也表明本文方法的有效性。然而，本文算法雖然可以解決目標密集的問題，但運算量是一個不得不考慮的問題。在今后的研究中，應(yīng)著重提升算法的效率，以便在實時性要求比較高的場合能夠得到應(yīng)用。

References)

[1]盧大威, 高彥釗, 晏行偉, 等.基于角度觀測關(guān)聯(lián)對的NN-EKF紅外/毫米波融合跟蹤算法[J]. 彈箭與制導學報, 2011,31(5):4-8.

[2]耿利祥,陳錢,錢惟賢.改進的聯(lián)合概率數(shù)據(jù)關(guān)聯(lián)算法(JPDA)對紅外目標與誘餌的辨別[J].紅外與激光工程,2013,42(2):305-310.

[3]劉德浩,王國宏,陳中華.系統(tǒng)誤差下的模糊航跡關(guān)聯(lián)方法[J].現(xiàn)代防御技術(shù), 2012,40(3):128-131；143.

[4]郭徽東,章新華.基于模糊綜合函數(shù)的航跡關(guān)聯(lián)算法及其應(yīng)用[J].系統(tǒng)工程與電子技術(shù),2003,25(11):1401-1403.

[5]丁鋒,薛磊.基于模糊理論的證據(jù)航跡關(guān)聯(lián)算法分析[J].計算機應(yīng)用研究,2015,32(6):1694-1697.

[6]毛藝帆,王睿,張金成.基于空天信息的異類傳感器模糊航跡關(guān)聯(lián)算法[J].空軍工程大學學報(自然科學版),2012,13(5):35-39.

[7]閆子豪.基于灰關(guān)聯(lián)分析-模糊推理的目標關(guān)聯(lián)技術(shù)[J].華東理工大學學報(自然科學版),2007,33(5):712-715.

[8]李鵬,劉思峰.基于灰色關(guān)聯(lián)分析和D-S證據(jù)理論的區(qū)間直覺模糊決策方法[J].自動化學報,2011,37(8):993-998.

[9]宋強,熊偉,何友.基于復(fù)數(shù)域拓撲描述的航跡對準關(guān)聯(lián)算法[J]. 宇航學報, 2011, 32(3):560-566.

[10]趙鳳,劉漢強.不同形狀鄰域空間信息的模糊聚類圖像分割[J].計算機工程與應(yīng)用, 2015,51(10):12-15；35.

[11]鄧敏,李志林,吳靜.空間關(guān)系理論與方法[M].北京:科學出版社,2013.

[12]鄧敏.矢量GIS中拓撲關(guān)系的擴展模型、理論和方法[D].武漢:武漢大學,2003：85-90.

[13]鄧敏,劉文寶,馮學智.GIS面目標間拓撲關(guān)系描述的形式化模型[J].測繪學報, 2005,34(1):85-91.

[14]楊哲,韓崇昭,李晨,等.基于目標之間拓撲信息的數(shù)據(jù)關(guān)聯(lián)方法[J].系統(tǒng)仿真學報,2008,20(9):2357-2360.

[15]宋強,熊偉,馬強.基于目標不變信息量的模糊航跡對準關(guān)聯(lián)算法[J]. 系統(tǒng)工程與電子技術(shù),2011, 33(1):190-195.

(編輯：李江濤)

A Track Association Approach Based on Topological Triangle Similarity Degree

MAO Yifan，ZHANG Duolin，WANG Lu

(Air and Missile Defense College, Air Force Engineering University, Xi’an Shaanxi 710051, China)

To solve out the challenges over the track association and fusion with system error when measuring track with multiple sensors, the paper proposes a track association method based on topological 4-intersection-and-difference model and triangle similarity degree. By dividing the detected target of multiple sensors into several topological triangle similarity structures, using 4-intersection-and-difference model, the paper excludes those triangles with small possibility of association , then the paper defines the triangle similarity degree and measures the topological triangle similarity degree so as to draw the relevance among targets and the result of track association. The simulation result shows, the algorithm is effective in case of intensive formation targets.

track fusion; track association; topological 4-intersection-and-difference model; triangle similarity degree

2016-01-25

毛藝帆(1988—)，女，博士研究生，主要研究方向為作戰(zhàn)建模與仿真。myf1210@126.com

TP212

2095-3828(2016)04-0069-06

A

10.3783/j.issn.2095-3828.2016.04.015

張多林，男，教授，博士生導師。