張 巍 唐耿彪 林立宇

（1.中科芯集成電路有限公司 無錫 214062）（2.武漢大學(xué)測繪遙感信息工程國家重點實驗室 武漢 430079）

1 引言

紅外探測技術(shù)在紅外預(yù)警、紅外制導(dǎo)等軍事領(lǐng)域以及醫(yī)學(xué)病變細(xì)胞診斷、工業(yè)探傷等民用領(lǐng)域中的應(yīng)用日趨廣泛［1～3］。然而，在紅外目標(biāo)檢測任務(wù)中，遠(yuǎn)距離成像通常使得紅外小目標(biāo)沒有明顯的紋理特性和特定的形狀［4～7］。因此，在復(fù)雜背景下檢測紅外小目標(biāo)是一項極具挑戰(zhàn)性的任務(wù)［8］。

在早期的研究中，Max-Mean［9］和 Top-Hat［10］等濾波算法被直接應(yīng)用于紅外小目標(biāo)檢測，雖然算法實現(xiàn)簡單但實際效果不佳。近年來，基于HVS的方法被廣泛研究并取得令人矚目的進(jìn)展［11～14］。ILCM［11］利用固定尺寸的局部窗口結(jié)構(gòu)測量局部對比度，在檢測率上有所提升，但難以檢測未知大小的目標(biāo)。針對這個問題，多尺度的局部對比度測量方法（LCM）［12］被提出，它能有效提高對未知大小目標(biāo)的檢測率，但是在抑制雜波方面效果仍然不佳?；诰矸e神經(jīng)網(wǎng)絡(luò)的一些方法［15～16］因其樣本依賴過大而應(yīng)用場合有限。隨后，MPCM［13］被提出用于抑制邊緣背景，然而它難以進(jìn)一步抑制背景并降低虛警率。雖然MPCM算法采用多尺度的局部窗口結(jié)構(gòu)能適應(yīng)不同大小的目標(biāo)，但由于窗口的被設(shè)置為規(guī)則大?。ㄈ?×3或5×5），所以難以測量尺寸大小不規(guī)則的目標(biāo)，導(dǎo)致了較高的虛警率。

為解決上述問題，本文提出了一種基于無監(jiān)督聚類的局部窗口結(jié)構(gòu)用于測量局部對比度。大量實驗證明，相比于已有的檢測算法，該方法具有更好的檢測性能和背景抑制效果。

2 基于無監(jiān)督聚類的紅外小目標(biāo)檢測方法

本文提出的算法具體步驟如圖1所示。首先利用滑動窗口快速提取候選像素，然后用基于無監(jiān)督聚類的窗口結(jié)構(gòu)分離目標(biāo)與背景，接著計算局部對比度用于背景抑制及目標(biāo)增強(qiáng)，最后在增強(qiáng)圖中通過自適應(yīng)閾值分割獲得檢測結(jié)果。

圖1 基于無監(jiān)督聚類的紅外小目標(biāo)檢測算法流程圖

2.1 候選像素的快速提取

如圖2所示，采用一個3×3圖像塊組成的滑動窗口快速獲取候選像素。這里，T和Bi（i=1，2，…，8）分別對應(yīng)檢測目標(biāo)和背景區(qū)域。為了檢測任意大小的目標(biāo)（范圍從 2× 1到 9 × 9）［11］，本文將圖像塊大小設(shè)置為9×9像素。

然后，使用滑動窗口對原始紅外圖像I進(jìn)行處理，預(yù)處理后的圖像E可表示如下：

其中（x，y）表示當(dāng)前處理的像素坐標(biāo)，mi指的是背景塊Bi的灰度均值，在預(yù)處理后的圖像E中，亮度高的像素更可能是真實目標(biāo)，故而將圖像E用以下的閾值進(jìn)行分割：

其中μ和δ分別指預(yù)處理后的圖像E的灰度平均值和標(biāo)準(zhǔn)差，k是一個經(jīng)驗參數(shù)，文中設(shè)置為6。經(jīng)過閾值分割后的像素組成候選像素集合。

2.2 基于無監(jiān)督聚類窗口構(gòu)造

構(gòu)造一個基于無監(jiān)督聚類的窗口結(jié)構(gòu)是為了可以在抑制各種背景雜波的同時增強(qiáng)被檢測目標(biāo)。圖3展示了一個以候選像素（x，y）為中心的構(gòu)造窗口。此處T指的是目標(biāo)區(qū)域，B1和B2分別為亮背景和暗背景，它們的形狀取決于無監(jiān)督聚類的結(jié)果。（x1，y1）（x2，y2）（x3，y3）（x4，y4）分別對應(yīng)區(qū)域中最左、最上、最右、最下點的坐標(biāo)。本文將窗口設(shè)置為15×15像素。然后使用模糊C均值算法（Fuzzy C-Means，F(xiàn)CM）對這些像素進(jìn)行聚類。傳統(tǒng)的FCM的核心思想是將每個類中心值與每個像素灰度值之間的差值的和函數(shù)最小化。由于傳統(tǒng)方法沒有考慮空間信息，使得分割效果不理想。因此本文將目標(biāo)函數(shù)改寫為

圖3 無監(jiān)督窗口結(jié)構(gòu)

其中c表示類的數(shù)量，本文設(shè)置為3，Ij表示第j個像素的灰度值，Ω表示圖像的像素集，μjk表示隸屬權(quán)值，且滿足，q表示模糊器值，通常設(shè)置為2，vk則表示每個類的中心值。ηjk是添加的空間權(quán)重項，（x0，y0）是潛在的小目標(biāo)中心坐標(biāo)，（xj，yj）是局部區(qū)域像素j的坐標(biāo)。本文中，將局部區(qū)域分為三類，第一類小目標(biāo)，第二類和第三類分別為亮背景和暗背景。這一改動使得聚類有距離參數(shù)，能更加精確地分離小目標(biāo)和復(fù)雜背景。

2.3 局部對比度計算與自適應(yīng)閾值分割

獲得聚類結(jié)果以后，計算每個候選像素的局部對比度進(jìn)行目標(biāo)增強(qiáng)和背景抑制。因為紅外小目標(biāo)通常比其背景更亮，本文定義FCM_LCM來描述目標(biāo)與背景之間的亮度差異。

這里，MT、MB1和MB2分別指的是目標(biāo)、亮背景和暗背景的亮度平均值。經(jīng)過計算以后，通常會發(fā)現(xiàn)目標(biāo)的FCM_LCM值是顯著的，用閾值分割即可得到檢測結(jié)果。因此使用以下的自適應(yīng)閾值來提取目標(biāo)：

此處Fmax和Fmean分別指最終增強(qiáng)圖像F中的像素的灰度最大值和平均值，λ是一個經(jīng)驗值，在本文中設(shè)置為0.5。

3 實驗驗證分析

為驗證本文方法的有效性，選用不同場景下的真實紅外序列來進(jìn)行實驗，四個數(shù)據(jù)集共計超過1200幀圖像。實驗中的β值設(shè)置為50。對比算法中的參數(shù)均設(shè)置為他們文中的推薦值。實驗的計算機(jī)配置為內(nèi)存8GB的Intel i5 3.3GHz雙核CPU，軟件平臺為Matlab。

3.1 評估指標(biāo)

為了直觀比較不同算法的背景抑制效果，將各個算法的檢測結(jié)果圖像進(jìn)行展示。為了進(jìn)一步定量比較不同算法的性能，本文選擇檢測概率Pd和虛警率Fd作為主要評價指標(biāo)［15］，其定義如下：

3.2 實驗結(jié)果與分析

本文選用改進(jìn)的局部對比度測量ILCM［11］、LCM［12］、多尺度基于塊的對比度測量MPCM［13］等高性能的檢測算法作為對比方法。圖4展示了不同算法在復(fù)雜背景中檢測目標(biāo)的結(jié)果圖。從圖中可以直觀看到ILCM殘留大量背景雜波，且漏檢目標(biāo)。LCM能有效檢測目標(biāo)，但仍有大量背景未能移除。MPCM表現(xiàn)出更好的背景抑制效果，但依然存在少量背景雜波難以移除。相比之下，本文的算法不僅能精確檢測到目標(biāo)，且移除了幾乎所有的背景干擾，具有最佳的背景抑制效果。

圖4 不同算法的檢測結(jié)果圖對比

在表1中，本文定量比較了不同方法的檢測性能?？梢钥闯?，本文的算法在三個數(shù)據(jù)集中都有最高的檢測率和更低的虛警率，這表明提出的方法有最佳的檢測性能。

表1 不同算法的檢測性能對比

4 結(jié)語

基于無監(jiān)督聚類的窗口結(jié)構(gòu)實現(xiàn)了目標(biāo)與背景的精確分離，充分利用了目標(biāo)的分布特征，提出了FCM-LCM描述子，有效增強(qiáng)了目標(biāo)，抑制了背景。大量圖像的實驗結(jié)果表明，與其他基線方法相比，本文的方法具有更好的檢測性能和對復(fù)雜背景的抑制效果。