高 飛，徐才進，梁 璟，倪辰辰，楊 玲

（中國電子科技集團公司第二十九研究所，四川 成都 610036）

0 引言

隨著現(xiàn)代戰(zhàn)場環(huán)境日益復(fù)雜以及戰(zhàn)場目標(biāo)活動愈發(fā)頻繁，傳統(tǒng)目標(biāo)情報中包括的基本屬性、信號參數(shù)等靜態(tài)信息已經(jīng)難以滿足聯(lián)合作戰(zhàn)對支援情報的需求。基于歷史偵察數(shù)據(jù)的目標(biāo)活動規(guī)律分析可以挖掘出重點目標(biāo)的典型航路和頻繁活動區(qū)域，不僅能為聯(lián)合作戰(zhàn)指揮控制、作戰(zhàn)籌劃等業(yè)務(wù)提供戰(zhàn)場目標(biāo)數(shù)據(jù)支撐，還能在戰(zhàn)場態(tài)勢感知中用于目標(biāo)活動異常檢測和預(yù)警等。受益于戰(zhàn)場偵察能力的提升和情報信息系統(tǒng)能力的增強，戰(zhàn)場目標(biāo)的歷史偵察數(shù)據(jù)在系統(tǒng)中持續(xù)積累，為目標(biāo)活動分析提供了良好的條件。作為戰(zhàn)場情報分析的主要組成部分，目標(biāo)活動分析依賴戰(zhàn)場偵察數(shù)據(jù)開展工作，傳統(tǒng)的目標(biāo)活動分析方法[1]以偵察數(shù)據(jù)統(tǒng)計、挖掘為主，輔以人工標(biāo)注等干預(yù)手段，雖然能得到較為完整準(zhǔn)備的目標(biāo)活動分析結(jié)論，但是存在挖掘算法設(shè)計困難、大數(shù)據(jù)計算耗時較長以及人工標(biāo)注效率較低等問題。蓬勃發(fā)展的圖像處理技術(shù)為目標(biāo)活動分析提供了新的技術(shù)思路，將偵察數(shù)據(jù)在地圖上進行圖像化，使用圖像處理方法進行目標(biāo)活動分析，可以有效減少算法設(shè)計難度、提升數(shù)據(jù)分析效率并增強分析結(jié)論可視化水平。

本文創(chuàng)新性地將圖像處理方法引入目標(biāo)活動規(guī)律分析。由于環(huán)境和裝備能力的影響，不同目標(biāo)的歷史偵察數(shù)據(jù)情況有所不同，根據(jù)數(shù)據(jù)量的數(shù)據(jù)分布情況可以分為歷史稀疏航跡點和歷史密集航跡點2 種情況：對于稀疏航跡點通過偵察數(shù)據(jù)圖像化、二值化以及圖像濾波實現(xiàn)了目標(biāo)活動航路提?。粚τ诿芗桔E點通過偵察數(shù)據(jù)熱度圖生成、圖像去噪、圖像增強實現(xiàn)了目標(biāo)活動區(qū)域提取。最后，使用仿真數(shù)據(jù)驗證了將圖像處理技術(shù)引入目標(biāo)活動分析的有效性。

1 稀疏航跡點目標(biāo)活動航路分析

針對稀疏型航跡點數(shù)據(jù)，主要研究基于歷史稀疏型航跡點的航路分析，具體處理分析流程如圖1 所示。

圖1 稀疏航跡點目標(biāo)活動航路分析流程

1.1 原始數(shù)據(jù)解析與航跡圖平滑化

通過戰(zhàn)場情報偵察等手段獲取的目標(biāo)位置信息主要是經(jīng)緯度信息，在進行基于數(shù)字圖像處理的活動航路分析前，需要將這種地理坐標(biāo)系統(tǒng)轉(zhuǎn)換為平面坐標(biāo)系統(tǒng)，形成平面坐標(biāo)航跡圖?？紤]到偵察時的環(huán)境影響、裝備狀況和目標(biāo)混疊等情況，原始偵察數(shù)據(jù)的經(jīng)緯度信息中可能存在一定誤差和噪聲，可以采用數(shù)字圖像處理中平滑模糊濾波方法[2]，使航跡圖像平滑模糊化并具有一定的容錯性。

圖像空間濾波器一般由一個鄰域和對鄰域所包圍圖像像素執(zhí)行的預(yù)定義操作組成。對每個像素濾波器產(chǎn)生一個新像素，其坐標(biāo)等于鄰域中心坐標(biāo)，其像素值是濾波的響應(yīng)結(jié)果。濾波器依次訪問輸入圖像中的每個像素后，就生成了濾波后的圖像。對圖像中的任意一點(x，y)，其像素值為f(x，y)，濾波器的響應(yīng)g(x，y)是濾波器系數(shù)與由該濾波器所包圍的圖像像素的乘積之和。濾波器鄰域一般為奇數(shù)尺寸，最小為3×3。對大小為m×n的鄰域，設(shè)奇數(shù)m=2a+1，奇數(shù)n=2b+1，其中a、b為正整數(shù)。使用大小為m×n的濾波器對大小為M×N的圖像進行線性空間濾波，可表示為：

例如，對于3×3 鄰域的線性空間濾波器：

1.2 航跡圖二值化

進行目標(biāo)活動航路分析時，圖像中只存在航路區(qū)域和非航路區(qū)域，是一個典型的二值化問題。圖像的二值化就是通過確定一個閾值，把灰度圖像的每個像素按其灰度值分為2 類，使得原圖像變?yōu)槎祱D像，即灰度圖像的二值化[3]。

設(shè)灰度圖像P尺寸為m×n，經(jīng)過上文濾波后處理的每一像素灰度值為x（i，j）。圖像像素灰度均值為：

設(shè)整數(shù)min、max，其中min 為最小灰度像素灰度值向下取整，max 為最大灰度像素灰度值向上取整。

設(shè)整數(shù)k∈［min，max］，對圖像進行遍歷，得到圖像中灰度值近似為k的像素點個數(shù)Ak。遍歷開始時Ak=0，對任一像素點，若k-0.5≤x（i，j） ＜k+0.5，則Ak=Ak+1。遍歷后得到序列Amin，Amin+1，…，Amax，對數(shù)列A進行排序，由大到小得到新的序列B1，B2，…，Bm，與B1相對應(yīng)的An的角標(biāo)n即為圖像中像素出現(xiàn)頻率最高的灰度，由此可得到灰度分布中一個較大的峰值處的灰度值。通常圖像的灰度均值應(yīng)在兩峰值之間，故可用下述方法確定另一峰值的位置，即按順序?qū)2，B3，…對應(yīng)An＇的角標(biāo)n＇代入式|n-n＇|，當(dāng)出現(xiàn)第一個滿足條件的n＇時，認為該n＇為另一波峰灰度值。所要取的閾值T應(yīng)在n與n＇之間 ：

取得閾值后，對圖像進行二值化運算得到邊緣清晰的二值化圖像：

1.3 航路區(qū)域提取

提取出二值化后的航跡圖像后，研究采用數(shù)字圖像處理中的數(shù)學(xué)形態(tài)學(xué)濾波方法，分析提取出航路區(qū)域。數(shù)學(xué)形態(tài)學(xué)可以很好地對圖像中存在的一些空洞或斷點進行處理，得到更加本質(zhì)的圖像形態(tài)，主要包括膨脹、腐蝕、細化、粗化、開運算和閉運算等。

膨脹是將目標(biāo)區(qū)域的背景點合并到該目標(biāo)物中，使目標(biāo)物邊界向外部擴張的處理。選取結(jié)構(gòu)元素B，將圖像X中的每一個像素點擴大為B+x，膨脹可以使圖像斷點連合，成為一個整體。

腐蝕是將目標(biāo)區(qū)域的邊界收縮，直觀上就是使得目標(biāo)區(qū)域變小，與膨脹是互為對偶的過程，但是并非經(jīng)過2 種運算后圖像可以還原。選取結(jié)構(gòu)元素B，將圖像中的X的B+x收縮為點x。腐蝕可以消除邊界處的干擾點。

另外，開運算和閉運算是為了解決膨脹和腐蝕運算改變圖像大小而提出來的。開運算是先腐蝕后膨脹的過程，可以消除圖像中細小噪聲，平滑物體邊界。閉運算是先膨脹后腐蝕的過程，可以消除圖像中存在的空洞，平滑物體邊界。

2 密集航跡點目標(biāo)活動區(qū)域分析

針對密集型航跡點數(shù)據(jù)，主要研究基于歷史密集型航跡點的區(qū)域分析，具體處理分析流程如圖2 所示。

圖2 密集航跡點目標(biāo)活動區(qū)域分析流程

2.1 原始數(shù)據(jù)解析與航跡點熱度圖生成

通過解析原始數(shù)據(jù)，將經(jīng)緯度信息轉(zhuǎn)換為平面坐標(biāo)系下的航跡圖，便于進行圖像處理。統(tǒng)計每個像素點區(qū)域目標(biāo)出現(xiàn)的熱度（像素點位置出現(xiàn)一次航跡點，便將像素點灰度值加1 或者其他固定值），分析處理后，目標(biāo)出現(xiàn)越頻繁的像素點區(qū)域灰度值越高，圖像對應(yīng)區(qū)域就越高亮顯示。

2.2 航跡點熱度圖去噪

通過圖像學(xué)中值濾波去除熱度圖中的“非頻繁噪聲點”。中值濾波的基本原理是把數(shù)字圖像或數(shù)字序列中一點的值用該點的一個鄰域中各點的中值代替，從而消除孤立的噪聲點[4]。

設(shè)二維矩［xij］表示一副數(shù)字化的圖像（其中i、j表示各像素點的位置，xij表示各像素點的灰度值），W[xij] 表示對圖像［xij］中的點xij做窗口操作，Med{xij}表示對窗口W[xij]內(nèi)的所有點取中值，若［yij］為圖像［xij］經(jīng)過中值濾波后的輸出圖像，濾波窗口的大小為m×n（其中m、n為奇整數(shù)），那么：

2.3 航跡點熱度圖增強

在數(shù)字圖像處理中，直方圖均衡化（HE）是一種最常用的空域直方圖修正技術(shù)，其以概率理論為基礎(chǔ)，運用灰度點運算來實現(xiàn)對圖像直方圖的變換，從而達到增強圖像對比度的目的。HE 是把給定圖像的直方圖分布改造成均勻直方圖分布，使輸出像素灰度的概率密度均勻分布?；驹硎钱?dāng)圖像中所有的灰度級出現(xiàn)的概率是一個均勻分布時，圖像所暴露出來的信息量最大[5]。

空域圖像增強的基本方法是對圖像做灰度變換T：f(x，y)→g(x，y)。直方圖修正試圖尋找灰度變換函數(shù)T：它能把f(x，y)映射為g(x，y)，使得g(x，y)的直方圖Pg(r)具有指定的形式qr，因此它不是事先指定變換函數(shù)，而是從Pf(r)及qr中計算得到T。任一圖像f(x，y)的直方圖Pf(r)在變換s=T(r)下具有以下形式：

式中，T(r)在區(qū)間［0，1］上是單調(diào)上升的單值函數(shù)，對于0 ≤r≤1，有0 ≤T(r)≤1。實際應(yīng)用當(dāng)中，處理的圖像都是數(shù)字圖像，因此需要離散化，HE 的離散化表達式為：

式中，NL、NS分別是圖像行數(shù)和列數(shù)，從而可以得到：

3 仿真驗證

基于上述目標(biāo)活動航路和活動區(qū)域分析方法，構(gòu)造對應(yīng)的稀疏、密集仿真數(shù)據(jù)，驗證其理論算法和分析流程的有效性。

3.1 稀疏航跡點目標(biāo)活動航路分析驗證

構(gòu)建稀疏的仿真目標(biāo)活動位置信息，形成原始的平面坐標(biāo)點跡圖，如圖3 所示。

圖3 原始點跡圖

基于時間序列繪制其航路，中間缺失的航跡點采用直線連接表示，如圖4 所示。

圖4 原始航跡圖

采用高斯濾波進行原始航跡圖平滑化處理。根據(jù)航跡點在圖像上的成像分辨率大小，選擇一定半徑的高斯核對航跡圖像進行平濾波滑，針對仿真數(shù)據(jù)選擇的高斯核大小為5×5，濾波效果如圖5 所示。

圖5 航跡圖平滑模糊濾波效果

航跡圖平滑模糊濾波后，進行圖像二值化處理，便于進行航路區(qū)域分析提取，試驗效果如圖6 所示。

圖6 航跡模糊圖像二值化

在二值化航跡圖像的基礎(chǔ)上，對二值化圖像重復(fù)進行3×3 或5×5 等區(qū)域大小的圖像學(xué)膨脹與腐蝕處理，重復(fù)進行2 至3 次，試驗效果如圖7 所示。

圖7 航路區(qū)域提取效果

為驗證航路分析提取效果，將其與原始航跡數(shù)據(jù)疊加進行驗證，試驗效果如圖8 所示。

圖8 航跡點航路驗證

3.2 密集航跡點目標(biāo)活動區(qū)域分集驗證

構(gòu)建密集的仿真目標(biāo)活動位置信息，并將地球經(jīng)緯度轉(zhuǎn)換為平面像素坐標(biāo)系下的絕對坐標(biāo)值，形成如圖9 所示的原始航跡點圖像。

圖9 原始航跡點圖

在原始航跡點數(shù)據(jù)基礎(chǔ)上，通過灰度圖亮度值累加，生成航跡點熱度圖效果，如圖10 所示。

圖10 航跡熱度圖

在航跡點熱度圖基礎(chǔ)上，利用圖像中值濾波原理，去除非頻繁活動航跡點，從而實現(xiàn)航跡點熱度圖去噪，效果如圖11 所示。

圖11 航跡熱度圖中值濾波效果

在通過中值濾波實現(xiàn)航跡熱度圖去噪后，采用直方圖均衡化方法，拉升圖像對比度，實現(xiàn)圖像增強，將航跡熱度圖顯性化，效果如圖12 所示。

圖12 航跡熱度圖增強效果

在航跡熱度圖增強效果基礎(chǔ)上，進行圖像二值化處理、圖像濾波，實現(xiàn)航跡熱度圖區(qū)域的擴張和活動區(qū)域顯性化，效果如圖13 所示。

圖13 航跡熱度圖頻繁活動區(qū)域效果

將提取的頻繁活動區(qū)域與原始數(shù)據(jù)進行疊加驗證，可以看出頻繁活動區(qū)域呈現(xiàn)出一些特定形狀區(qū)域、直線區(qū)域或點狀區(qū)域等，如圖14 所示。

圖14 頻繁活動區(qū)域效果驗證

將通過仿真數(shù)據(jù)提取到的活動航路、區(qū)域結(jié)果與原始航跡點進行疊加驗證，可以看出該方法得到的結(jié)論能有效地描述不同數(shù)據(jù)分布情況下的目標(biāo)典型航路和區(qū)域信息，結(jié)果邊界光滑且對比顯示清晰。

4 結(jié)束語

本文將數(shù)字圖像處理的相關(guān)方法引入目標(biāo)活動分析，對于稀疏航跡點通過航跡平滑化、二值化以及形態(tài)濾波實現(xiàn)了目標(biāo)活動航路提取，對于密集航跡點通過熱度圖生成、去噪、增強、濾波實現(xiàn)了目標(biāo)活動區(qū)域提取。并通過仿真數(shù)據(jù)驗證了方法的有效性。相對于傳統(tǒng)的基于統(tǒng)計和人工標(biāo)注的方法，本文提出的方法能有效消除數(shù)據(jù)噪聲和非頻繁活動點的影響，通過形態(tài)學(xué)濾波消除由于偵察數(shù)據(jù)缺失導(dǎo)致的邊界不平滑問題，并實現(xiàn)提取結(jié)果的增強顯示。后續(xù)可以基于本文的活動分析結(jié)果，同樣利用圖像處理的方法判斷目標(biāo)當(dāng)前是否在典型活動航路或區(qū)域內(nèi)，判斷其是否異?；顒?。■