基于顯著條件隨機場的海面目標分割算法

宋文韜 胡勇 鞏彩蘭 李賽

摘 ?要： 海面目標的分割在軍事和民用領域有著重要的價值，但現階段對光學遙感下海面艦船等目標的研究還較少。針對強雜波干擾海面場景下目標分割精度低的問題，文中設計一種顯著條件隨機場的目標分割算法。該算法基于背景可分性度量準則，改進了譜殘差顯著性，并使用全連接的隨機場替代一般的鄰域模型，利用全圖的上下文信息對目標精確分割。實驗結果表明，該算法能夠有效抑制海面雜波，實現了對像素的精確分割。

關鍵詞： 圖像分割; 視覺顯著性; 條件隨機場; 海面目標; 遙感; 艦船目標檢測; 圖像處理; 目標分割

中圖分類號： TN911.73?34; TP751.1 ? ? ? ? ? ? 文獻標識碼： A ? ? ? ? ? ? ? ? ?文章編號： 1004?373X（2019）10?0128?04

Sea?surface target segmentation algorithm based on saliency condition random field

SONG Wentao1.2，3， HU Yong1，2， GONG Cailan1，2， LI Sai1，2，3

（1. Shanghai Institute of Technical Physics of Chinese Academy of Sciences， Shanghai 200083， China;

2. Key Laboratory of Infrared System Detection and Imaging Technology， Chinese Academy of Sciences， Shanghai 200083， China;

3. University of Chinese Academy of Sciences， Beijing 100049， China）

Abstract： The sea?surface target segmentation has an important value in military and civilian fields. However， there are few studies on sea?surface ship targets under optical remote sensing at present. In allusion to the low segmentation accuracy of targets in the strong clutter interfered sea scenario， a target segmentation algorithm based on the saliency condition random field is designed. In the algorithm， the spectral residual saliency is improved based on the background separability metrics criterion. The general neighborhood model is substituted by the fully?connected random field. The accurate segmentation of targets is conducted by using the context information of the whole image. The experimental results show that the algorithm can effectively suppress sea clutters and achieve accurate segmentation of pixels.

Keywords： image segmentation; visual saliency; condition random field; sea?surface target; remote sensing; ship target detection; image processing; target segmentation

光學遙感相比傳統的雷達探測[1] ，具有強隱蔽性、環境適應性好的特點，得到越來越廣泛的應用。

現在對光學遙感海面目標的檢測與分割運用較多的是基于顯著度的方法。顯著度是一種主動模仿人眼特征的設計方法，模擬人眼可快速集中于感興趣區域，而提出的目標檢測算法模型[2] 一般可分為3類：基于對比度信息[3] 、基于頻率域特征[4]、基于計算的特征融合的算法[5] 。但對不同場景適應性差，會提取一定的錯誤區域;經顯著度處理后都有一個高斯濾波的過程，難以做到對像元的精確提取。概率圖模型是圖像分割領域的一個研究熱點，條件隨機場模型（Condition Random Field，CRF）是經典的馬爾科夫隨機場模型（MRF）的發展，主要有兩個優點：首先，去除觀測序列的獨立假設分布約束[6]，增強馬爾科夫隨機場模型的適應能力;其次是能夠對長距離和短距離像素關系進行建模[7] ，特別是對于海面這種大面積背景，能夠充分利用圖像的上下文關系。本文針對遙感海面圖像的特點，對現有條件隨機場的檢測框架進行改進，提出一種改進視覺顯著性一元勢函數，同時利用全連接的隨機場替代鄰域檢測模型，實現對海面艦船目標的精確分割。

1 ?改進視覺顯著性模型

1.1 ?基于譜殘差顯著性的目標分割算法

人眼能夠快速和精確地定位自然場景中的目標，并且在不需要事先識別目標的情況下，轉移注意力并聚焦該目標，同時將目標從背景中分割出來。基于此，人眼會對圖像中的變化更加敏感，一個基本原則就是抑制常規與頻繁出現的特征，同時對非常規的特征保持敏感。依據此將圖像分為背景和異常兩部分[8] ：

[H（Image）=H（Innovation）+H（Backgroung）] （1）

給定一幅圖像，首先進行二維離散傅里葉變換，將其從空間域轉換到頻域，對幅值取對數，根據Hou x等人提出的理論，該對數曲線滿足局部線性條件[9]，基于此先驗得到譜殘差（Spectral Residual，SR）顯著性。其步驟計算為：對于給定圖像[I（x）]，[A（f）]表示其經過傅里葉變換的幅度譜，對其取對數譜，得到[L（f）]：

[L（f）=log A（f）] （2）

對對數譜[L（f）]進行均值濾波，再將濾波結果與對數譜相減，得到剩余譜[R（f）]：

[R（f）=L（f）-hn（f）*L（f）] （3）

利用剩余譜和相位譜進行反向傅里葉變換，進行高斯濾波器，得到顯著圖：

[S（x）=g（x）*F-1[e[R（f）+P（f）]]2] （4）

式中：F表示傅里葉變換;P（f）表示其相位譜。

對顯著度進行閾值分割，得到目標的初始分割結果，一般選擇閾值threshold為3倍均值。

[threshold=3E（S（x））] （5）

將大于閾值的賦值為1，其余為0。

1.2 ?基于相對熵的背景可分性度量準則

一般的視覺顯著性模型，輸入任意的圖像，都會相應輸出高頻或者異常區域。與自然圖像不同，遙感廣域海面只有少部分影像包含待分割目標，這些顯著性方法已不再適用。針對遙感目標的局限性，本文提出一種基于背景可分性度量準則，對視覺顯著性進行改進，并對像元進行初始分割。相對熵（Relative Entropy）又稱為KL散度，是兩個概率分布一種度量準則。 對一個假設背景的概率分布為P和目標的概率分布為Q來說，他們的KL散度定義為：

[DKLPQ=iP（i）logP（i）Q（i）] （6）

對提取的顯著區域和背景進行計算，在高斯分布假設下通過定義最大容許的貝葉斯錯誤率可以計算KL離散度的最低閾值[10] ，符合條件的影像認為存在目標，將顯著分割結果作為其初始分割結果。

2 ?全連接條件隨機場精確分割

條件隨機場借助圖像的上下文關系，可以實現對顯著性初始分割結果的精確分割。一般對圖像隨機變量定義兩個場：一個是觀測場Y，對應到圖像像元的觀測值（如DN）的集合;一個是標記場，對應于像元的狀態空間。其中L是標記場X的配置空間，對于目標分割問題，L對應二元空間，即L={0，1}?；跅l件隨機場的聯合概率可用吉布斯分布描述[11]，如下：[PXY=1ZexpE（xy） ? ? ? ? ? ? ? ? ? =1Zexp-c∈C-VCXcY] （7）

式中：E在統計力學中表示[xy]的吉布斯能量，是基團勢函數的疊加;C為基團的集合，每個基團c都將產生一個勢函數V;Z是歸一化函數。

一般而言只考慮單點和雙點基團，即一元勢函數和成對勢（二元）函數，為了方便表示，用X來代替[XY]，所以其吉布斯能量為：

[E=iΨuXi+i≠jΨi，jXi，Xj] （8）

一元勢函數設置為視覺顯著度分割模型，其置信度M，一般設為50%，可得：

[ΨuXi=-ln M， ? ? ? ? ? ? ? ? ?像元為目標-ln1-M2， ? ? ? ? ?像元為背景] （9）

二元勢函數用高斯核函數表示，i，j為圖中任意像元的序號。定義二元勢函數為：

[Ψi，jXi，Xj=μXi，Xjm=1Kwmkmfi，fj] （10）

式中：[μ（Xi，Xj）]是變分模型，一般選Ising/Potts模型，即只有當像素i和j的標簽不同時，才會為1，其他情況下分0;每一個[km]都代表一個高斯核函數;[f]是由每個像素的圖像觀測值和位置特征構造成的特征向量;[wm]為權重。所以高斯核函數為[7]：[kfi ，fj=w1exp-Pi-Pj22θ2α-Ii-Ij22θ2β+ ? ? ? ? ? ? ? ? ? w2exp-Pi-Pj22θ2γ] （11）

式中：I為每個像素值;P為像素在圖像上的坐標。

為了降低優化的復雜度，用最小化KL散度[D（QP）]尋找分布Q（X）代替P（X），分布Q可以寫成獨立的邊緣概率分布的乘積[12]：

[QX=iQiXi]

根據概率性質有：[XiQi（Xi）=1]，利用拉格朗日乘數法求[argminQiD（QP）]，根據Q計算結果將每個像素概率最大狀態作為該像素的類別。

3 ?實驗結果

本文算法的測試環境為Core（TM） ?i7?2600 CPU@2.90 GHz，編程語言C++，選擇High resolution ship collections 2016（HRSC2016） 數據集[13] 進行對比實驗，以驗證本文方法在圖像分割中的效果。本節給出本文所提模型與基于迭代計算的ISODATA聚類分割算法、基于頻域視覺顯著度的FT算法[14] 、基于自適應閾值的Otsu算法的實驗結果，主要從分割視覺效果與分割精度兩個方面進行比較，結果對比圖如圖1所示。

圖1 ?海面目標分割結果對比圖

將本文算法與三類算法做對比，從圖中可以看出，ISODATA算法受海面雜波影響最大，分割結果不符合實際需求;FT顯著性和Otsu分割也會受到復雜海面影響，分割精度降低;而本文方法在3幅圖像視覺上都能保持較好的結果，魯棒性強。

為了對本文方法進行定量評價，以像元檢測精度R作為評價指標。

[R=正確分類像元數+錯誤分類像元數像元總數×100%] （12）

其結果如表1所示。

從表中可以看出，傳統方法在海面雜波較多時，容易標記出不屬于目標的像素，這給實際工程處理帶來了不必要的麻煩，受圖像質量影響比較大，本文方法能夠避免這種現象，對三種不同圖像都能保持較高的分割精度，有效抑制海面雜波的影響。

4 ?結 ?論

本文提出一種結合視覺顯著度和條件隨機場的艦船目標分割算法，改進了基于頻域信息的顯著性算法，使其能夠更突出異常目標，同時結合條件隨機場模型，更有效利用廣域遙感圖像的空間上下文信息，實現了對像素的精確分割。這種算法在遙感圖像數據處理領域具有一定的潛在應用前景。

參考文獻

[1] APRILE A， GROSSI E， LOPS M， et al. Track?before?detect for sea clutter rejection： tests with real data [J]. IEEE transactions on aerospace and electronic systems， 2016， 52（3）： 1035?1045.

[2] 徐芳，劉晶紅，曾冬冬，等.基于視覺顯著性的無監督海面艦船檢測與識別[J].光學精密工程，2017，25（5）：1300?1311.

XU Fang， LIU Jinghong， ZENG Dongdong， et al. Detection and identification of unsupervised ships and warships on sea surface based on visual saliency [J]. Optics and precision engineering， 2017， 25（5）： 1300?1311.

[3] GU K， ZHAI G， YANG X， et al. Automatic contrast enhancement technology with saliency preservation [J]. IEEE transactions on circuits and systems for video technology， 2015， 25（9）： 1480?1494.

[4] LI J， DUAN L Y， CHEN X， et al. Finding the secret of image saliency in the frequency domain [J]. IEEE transactions on pattern analysis and machine intelligence， 2015， 37（12）： 2428?2440.

[5] BORJI A. Boosting bottom?up and top?down visual features for saliency estimation [C]// Proceedings of IEEE Computer Society Conference on Computer Vision and Pattern Recognition. [S.l.]： IEEE Computer Society， 2012： 438?445.

[6] LAFFERTY J， MCCALLUM A， PEREIRA F. Conditional random fields： probabilistic models for segmenting and labeling sequence data [C]// Proceedings of the 18th International Conference on Machine Learning. San Francisco： Morgan Kaufmann Publishers Inc.， 2001： 282?289.

[7] KR?HENB?HL P， KOLTUN V. Efficient inference in fully connected CRFs with Gaussian edge potentials [J]. Advances in neural information processing systems， 2011， 24： 109?117.

[8] 陳海永，徐森，劉坤，等.基于譜殘差視覺顯著性的帶鋼表面缺陷檢測[J].光學精密工程，2016，24（10）：2572?2580.

CHEN Haiyong， XU Sen， LIU Kun， et al. Surface defect detection of steel strip based on spectral residual visual saliency [J]. Optics and precision engineering， 2016， 24（10）： 2572?2580.

[9] HOU X， ZHANG L. Saliency detection： a spectral residual approach [C]// Proceedings of IEEE Conference on Computer Vision and Pattern Recognition. Minneapolis： IEEE， 2007： 1?8.

[10] 劉豐軼，胡勇，饒鵬，等.一種改進的MRF點目標檢測算法[J].紅外與毫米波學報，2018，37（2）：212?218.

LIU Fengyi， HU Yong， RAO Peng， et al. Modified point target detection algorithm based on Markov random field [J]. Journal of infrared and millimeter waves， 2018， 37（2）： 212?218.

[11] 劉彤，黃修添，馬建設，等.基于完全聯系的條件隨機場的圖像標注[J].計算機應用，2017，37（10）：2841?2846.

LIU Tong， HUANG Xiutian， MA Jianshe， et al. Image labeling based on fully?connected conditional random field [J]. Journal of computer applications， 2017， 37（10）： 2841?2846.

[12] KR?HENB?HL P， KOLTUN V. Parameter learning and convergent inference for dense random fields [C]// Proceedings of the 30th International Conference on Machine Learning. Atlanta： IEEE， 2013： 513?521.

[13] LIU Z， YUAN L， WENG L， et al. A high resolution optical satellite image dataset for ship recognition and some new baselines [C]// Proceedings of the 6th International Conference on Pattern Recognition Applications and Methods. Porto： IEEE， 2017： 324?331.

[14] ACHANTA R， HEMAMI S， ESTRADA F， et al. Frequency?tuned salient region detection [C]// Proceedings of IEEE Conference on Computer Vision and Pattern Recognition. Miami： IEEE， 2009： 1597?1604.