基于Hausdorff度量的高光譜異常目標檢測算法

趙春暉，尤偉，齊濱2，王玉磊

（1.哈爾濱工程大學(xué)信息與通信工程學(xué)院，黑龍江哈爾濱150001；2.哈爾濱工程大學(xué)水聲工程學(xué)院，黑龍江哈爾濱150001）

基于Hausdorff度量的高光譜異常目標檢測算法

趙春暉1，尤偉1，齊濱2，王玉磊1

（1.哈爾濱工程大學(xué)信息與通信工程學(xué)院，黑龍江哈爾濱150001；2.哈爾濱工程大學(xué)水聲工程學(xué)院，黑龍江哈爾濱150001）

高光譜遙感異常目標檢測中，目標和背景光譜信息難以精確地界定，導(dǎo)致目標檢測性能下降。針對經(jīng)典RX檢測算法存在虛警概率高、計算量大、過程繁瑣等缺點，將Hausdorff度量引入高光譜異常檢測，利用改進的Hausdorff距離（MHD）從光譜匹配程度的角度，進行了高光譜異常目標檢測，最大程度地將異常目標和背景分離。采用模擬數(shù)據(jù)和真實高光譜數(shù)據(jù)進行大量實驗，檢測性能大幅提升，算法的計算效率提高了60%，證明了本文算法比RX算法、因果RX算法和KRX算法檢測效果好，效率高等優(yōu)勢，算法的低復(fù)雜度特性為硬件實現(xiàn)提供了良好的算法支持。

高光譜遙感；異常檢測；RX算法；光譜匹配；Hausdorff度量

近年來，高光譜遙感在目標探測與識別等方面具有了廣泛的應(yīng)用。在目標檢測中，不需要先驗光譜信息的異常目標檢測技術(shù)具有較強的實用性，經(jīng)典的異常檢測算法是源于多光譜遙感的RX算法［1-2］，RX算法最早是由Reed和Xiaoli Yu于1990年提出的，是用于在高斯背景統(tǒng)計特性和空間白化的條件下目標光譜信息和背景光譜信息均未知的目標檢測［3］。隨后Chang等在RX算法基礎(chǔ)上對檢測算子進行了改進［4］，2013年，Acito等提出了一種RX局部異常檢測算子（L-RXD），提高了檢測性能，Riley等聯(lián)合了歐氏距離構(gòu)造了新檢測器［2］，有效地抑制了噪聲的干擾。由于RX算法具有較大的計算量，導(dǎo)致了算法計算效率不高［5］。隨后，Chang等在RX算法的基礎(chǔ)上進行了進一步改進［6-8］，提出了高光譜目標因果RX算子（CRRXD）和實時檢測RX算法，利用迭代思想大大加快了計算速度，提高了算法效率，但尚未從算法本身提高檢測精度。為降低虛警概率，Kwon等提出了一種基于核機器學(xué)習的非線性RX改進算法（KRX）［9］，該算法充分地利用了高光譜波段之間的非線性統(tǒng)計特性，然而，大量非線性核函數(shù)的運算使得算法效率降低。

在模式識別中，豪斯多夫（Hausdorff）度量被廣泛應(yīng)用于人臉檢測和圖像配準［10-12］，例如，何力提出一種基于 Hausdorff距離的運動目標跟蹤算法［13］，在運動物體的檢測、識別和跟蹤取得了較好的效果，H.Peyvandi提出了將基于Hausdorff相似度量的徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)應(yīng)用于聲吶信號的探測和分類［14］，在較低信噪比情況下獲得了很好的檢測效果，Manuel Grana等提出了基于Hausdorff距離的高光譜數(shù)據(jù)特征檢索算法［15］，提高了傳統(tǒng)算法的效率和準確率。

針對現(xiàn)有RX異常檢測算法虛警概率高、計算量較大的問題，本文提出了一種基于Hausdorff度量的高光譜異常目標檢測算法。異常目標檢測問題是一個當前像元與背景光譜信息匹配相似度問題，由于使用了Hausdorff距離作為距離測度無需考慮兩個點集中元素之間的對應(yīng)關(guān)系，可以有效地降低噪聲的干擾，最大程度地將異常目標和背景分離，提高目標邊緣檢測概率，檢測性能大幅提升，同時，計算復(fù)雜度大大降低，提高了計算效率。

1　異常目標檢測RX算子

1.1RX異常檢測算子

設(shè)有L個波段的高光譜像元的光譜向量可表示為一個L維列向量ri=［r1i，r2i，...，rLi］T，假設(shè)已知光譜信號為s=［s1，s2，...，sL］T，令a=［a1，a2，...，aL］T為信號豐度，則二元假設(shè)檢驗可表示如下

式中:H0代表目標不存在，H1代表目標存在；n是背景雜波信號，這里采用加性模型。通過廣義似然比檢驗并經(jīng)過一系列化簡可以得出RX檢測器，給出RX簡化形式如下

其中

式中:r為L維待檢測像元的光譜向量，μ為高光譜數(shù)據(jù)樣本均值向量，K為高光譜數(shù)據(jù)樣本協(xié)方差矩陣，η為檢測閾值，N為高光譜數(shù)據(jù)中觀測樣本的總數(shù)。

基于相關(guān)矩陣的因果RX異常檢測算子:

其中

式中:rn是當前待檢測像元 （nth）的光譜特性，R（n）是采樣數(shù)據(jù)的相關(guān)矩陣，該相關(guān)矩陣由當前像元之前的所有像元光譜求相關(guān)性而得，具有因果特性，因此稱之為因果RX算子（CRRXD）。

1.2KRX異常檢測算子

核RX（Kernel RX）是Kwon等人針對RX算子存在的問題提出的改進算法，采用非線性函數(shù)Φ將原始的高光譜遙感數(shù)據(jù)映射到高維特征空間中，得到特征空間中KRX算法可表示為

其中

化簡最終得到KRX算子的最終算子表達式:

其中:KB=K（ PB，PB）表示一個中心化的Gram矩陣，通過非中心化的Gram矩陣得到，是估計均值，KB=K（ PB，PB）為定義N×N中心化Gram核矩陣。

2　基于Hausdorff距離的異常檢測算子

HD（Hausdorff distance）為在度量空間中任意兩個集合之間定義的一種距離，可以用來描述兩組點集之間的相似程度。

設(shè)兩個點集:A=｛a1，a2，a3，...，am｝，B=｛b1，b2，b3，...，bn｝

定義Hausdorff距離為

式中:h（A，B）和h（B，A）為兩個點集A和B的單向距離，分別定義為

式（11）也被稱為雙向HD距離，是最簡單的形式，‖·‖表示某種距離范數(shù)，如歐氏距離。兩個單向HD距離，通常是不一樣的，表示了其中一個集合中最不匹配的點到最近的另一個集合中的點的距離。應(yīng)用HD距離需要證明其在集合族內(nèi)部是完備度量。

定理:用κ表示Rn中全體非空緊致子集構(gòu)成的集合族，即 HD:κ× κ →R，HD（ A，B）= max｛ h（A，B），h（ B，A）｝，HD是κ上一個完備度量，稱作Hausdorff度量。

證明:根據(jù)度量的定義，度量需要滿足三個條件，即非負性，對稱性和三角不等式。非負性和對稱性顯然，下面證明三角不等式。設(shè)A，B，C∈κ根據(jù)緊致性，顯然有HD（ A，B＜∞）設(shè)ε＞0，那么存在b∈B使得

因此，存在z∈C使得

這證明A包含在C的h（A，B）+h（B，A）+2ε鄰域的閉包內(nèi)。類似地，C包含在A的h（A，B）+ h（B，C）+2ε鄰域的閉包內(nèi)。所以

由于ε＞0的任意性，證得三角不等式性質(zhì)

HD距離作為一種度量，等于兩個單向HD距離的最大值，它可以解釋為保證兩個集合互相包含所需要的最小半徑，兩個集合越不相似，其互相包含的半徑越大，因此，它衡量了兩個點集的最大不相似程度，可以廣泛應(yīng)用于圖像配準、目標識別和檢測。

在HD距離原始定義中，它對噪聲干擾的集合輪廓非常敏感，要使得兩組點集之間的HD距離小，就要保證一個集合的每一個元素都需要位于另一個集合的附近。假設(shè)一個集合中存在一個或被噪聲干擾的點不滿足這個條件，則這兩個集合的HD距離將完全取決于噪聲點的位置，適應(yīng)性大大降低。為提高其魯棒性，采用改進的HD距離，表示為

其中:

NA=m是A點集中元素的數(shù)量，改進的mh（A，B）表示為B集合中元素最小值與A中所有元素非對應(yīng)歐氏距離的均值。改進的HD距離稱為MHD距離，MHD距離可以有效地緩解A中噪聲對距離測度的敏感程度。

RX檢測算子本質(zhì)是利用馬氏距離進行待檢測像元和背景相似度檢測，引入Hausdorff距離度量解決了RX檢測的灰度值圖像高虛警問題，對比于KRX算法，本方法從根本上降低了計算復(fù)雜度，使得計算效率大大提高。

Hausdorff距離用于衡量兩組點集之間的相似程度，在異常目標檢測任務(wù)中，提取像元光譜維數(shù)據(jù)。Kwon提出了一種基于雙層窗特征分解的異常檢測算法。該算法利用外窗和內(nèi)窗兩個窗口，最大化分離背景向量和目標向量的光譜特征。其具體局部檢測模型如圖1所示。外層窗截取檢測背景信息，而檢測點位于內(nèi)窗中央，內(nèi)窗的其他部分認為是目標不做操作。

圖1　高光譜異常目標局部檢測模型Fig.1 Hyperspectral anomaly target local detection model

由于空間上同一個位置附近的背景光譜曲線差異不大，對經(jīng)過雙窗模型得到的背景光譜數(shù)據(jù)取均值作為平均背景光譜向量B，同時，內(nèi)窗提取的像元作為待檢測像元光譜向量A，對兩個向量求改進的HD距離作為待檢測像元的檢測結(jié)果。用MHD距離來度量兩個向量的相似度，兩個向量的數(shù)據(jù)越相似，其MHD距離越小，兩個向量的數(shù)據(jù)差別越大，其MHD距離越大，檢測結(jié)果就越大，超過閾值顯示檢測到異常。另外，對背景光譜數(shù)據(jù)取均值作為平均背景光譜向量可以弱化MHD距離對B中噪聲的敏感程度。因此，此方法可以提升原始HD距離對A和B的抗噪聲性能。實際上，在高光譜異常目標檢測中，相鄰像元光譜數(shù)據(jù)具有相似分布特性，從原理上，應(yīng)用MHD距離可以充分抑制噪聲對檢測結(jié)果的干擾，對兩個向量進行匹配是可以對各個元素的微小差異進行積累，它會更好的檢測相似度，因此，可以獲得更好的檢測結(jié)果。

基于MHD距離的高光譜異常檢測的計算過程:

1）對于高光譜數(shù)據(jù)進行預(yù)處理，求均值做歸一化。

2）設(shè)待檢測像元xi，高光譜數(shù)據(jù)波段數(shù)L，利用雙窗模型在排除目標像元的前提下最大化提取待檢測像元附近背景W=｛xi（j）｝，i=1，2，...，N，j=1，2，...，M，N是待檢測像元總數(shù)，M是待檢測像元背景像元個數(shù)。

4）分別計算待檢測像元到匹配背景像元和匹配背景像元到待檢測像元的MHD距離:

5）比較mh（xi，xim）和mh（xim，xi）取其中的最大值作為MHD距離，將MHD距離做閾值判決輸出結(jié)果。

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

為驗證本文提出的MHD算法的有效性和高效性，分別用合成的模擬數(shù)據(jù)和真實的高光譜數(shù)據(jù)進行了實驗。計算機的硬件配置，處理器型號為Intel Core I5，主頻2.5 GHz，內(nèi)存為DDR3 8 GB，仿真實驗平臺為MATLAB2014版本。

3.1高光譜數(shù)據(jù)描述

3.1.1模擬數(shù)據(jù)介紹

基于本文MHD算法針對異常目標檢測，取用真實數(shù)據(jù)AVIRIS高光譜遙感中多種異常目標進行目標插入，該AVIRIS高光譜數(shù)據(jù)有126個波段，數(shù)據(jù)共有400×400像元。圖2中圈出了四種地物的地理分布，該四種地物分別是油桶 （gasoline can，G）、房屋（house，H）、樹木（tree，T）和飛機（plane，P）。合成90×90空間大小的含有126個波段的高光譜數(shù)據(jù)進行仿真實驗。

該模擬數(shù)據(jù)共有4×4=16個異常點。該16個3×3大小的像元生成如下:用圖2的四種物質(zhì)（G，H，T，P）的光譜特性模擬四種異常目標，每行異常目標由同種物質(zhì)構(gòu)成。其中第一列由純像元構(gòu)成，第二列包含了75%異常像元光譜特性和25%背景光譜特性，第三列包含了50%異常像元光譜特性和50%背景光譜特性，第四列包含了25%異常像元光譜特性和75%背景光譜特性。該圖可用于不同背景干擾下的多種異常目標探測。圖3即為模擬數(shù)據(jù)與目標分布。

圖2　四類異常目標在真實數(shù)據(jù)中的分布Fig.2 The distribution of four types of anomaly targets in real data

圖3　模擬數(shù)據(jù)第100波段和目標分布Fig.3 Synthetic data 100 band image and the target distribution

3.1.2真實數(shù)據(jù)

為了驗證本文提出的算法在實際應(yīng)用中的有效性，采用一幅真實的高光譜數(shù)據(jù)進行仿真實驗。該數(shù)據(jù)是美國圣地亞哥海軍基地AVIRIS高光譜數(shù)據(jù)截取的100×100空間大小的數(shù)據(jù)。該數(shù)據(jù)有224個波段，去除水的吸收帶和信噪比較低的波段后，余下的126個波段用于實驗，圖4即為原始高光譜的第100波段灰度顯示和真實地物分布。

圖4　真實數(shù)據(jù)第100波段灰度圖和真實地物分布Fig.4 Real data 100 band image and the target distribution

3.2模擬數(shù)據(jù)仿真實驗及分析

利用模擬數(shù)據(jù)的目的在于測試6種算法對于不同背景干擾下的多種異常目標的抗噪聲性能。實驗中，局部算子選擇目標檢測窗口設(shè)為3×3像元，背景信息提取窗口大小設(shè)為11×11像元。首先對原始高光譜數(shù)據(jù)歸一化，然后根據(jù)合成數(shù)據(jù)大小和目標分布，對圖2所示的合成數(shù)據(jù)分別進行全局RX算法（G-RXD）、局部RX算法（L-RXD）、全局因果RX算法（G-CRRXD）、局部因果RX算法（L-RXD）、KRX算法和本文提出的MHD算法進行檢測。KRX算法中的核函數(shù)采用的是高斯徑向基（RBF）核函數(shù)k（ x·xi）=exp（‖x-xi‖2/δ2），將原始數(shù)據(jù)映射到高維特征空間，通過多次仿真實驗比較最終確定徑向基函數(shù)的寬度δ2=100。

圖5　檢測結(jié)果二值圖像對比Fig.5 Results of binary image contrast

對比最終檢測結(jié)果的二值圖像，如圖5所示。圖5（a）～（d）描述了采用RX檢測算子在設(shè)定的閾值下得到的異常檢測結(jié)果，基于馬氏距離的全局RX算法檢測效果很不理想，全局檢測算子存在漏檢目標的情況嚴重，圖5（c）、（d）是局部算子，但由于選擇的外窗的背景信息并不是完全服從正態(tài)分布，因此在16個異常目標中每一個都有不同程度的漏警，且在異常目標周圍受背景干擾產(chǎn)生很多虛警點。圖5（e）、（f）可以看出，KRX算法較RX算法檢測效果好很多，沒有產(chǎn)生漏警，只是在個別目標周圍產(chǎn)生了虛警點，相比較而言，MHD算法效果更好，幾乎沒有產(chǎn)生漏警和虛警點，對多種異常目標的檢測具有通用性。

為了量化分析本文算法的優(yōu)越性，以不同算子檢測的結(jié)果中的目標所占像元數(shù)量、虛警所占像元數(shù)量和計算時間作為性能指標，對上述算法進行比較分析，如表1所示。

如表1所示，由于MHD算法的計算復(fù)雜度較RX算法小很多，因此計算效率大大提高，KRX算法取得了較好的檢測效果但計算時間也有所增加，而MHD算法進一步減少了計算，極大地減少了計算時間?；具_到了實時處理的要求，因此該算法更具實用性。

表1　不同算法性能量化分析Table 1 The quantitative analysis of different algorithms

3.3真實數(shù)據(jù)實驗結(jié)果與分析

實驗中，對高光譜數(shù)據(jù)去均值做歸一化后，采用與模擬數(shù)據(jù)實驗相同的局部窗口和核函數(shù)。為了更直觀地表現(xiàn)HD距離用于檢測的優(yōu)越性，實驗對比兩種算法繪制出了檢測灰度圖像的3D峰度對比圖如圖6所示。其中XY軸二維平面上的點代表像元在灰度圖中位置，圖中藍色到紅色對應(yīng)的為灰度依次由小到大，三維Z軸縱坐標表示對應(yīng)像元灰度值的強弱（即峰度），由圖中的峰度對比得出結(jié)論，和KRX算法相比，采用本文的MHD算法的檢測結(jié)果能夠產(chǎn)生較低的漏警概率，更容易的選擇最佳閾值進行目標和背景的分離。

圖6較直觀地表現(xiàn)了傳統(tǒng)RX、KRX和MHD算法對圣地亞哥海軍基地數(shù)據(jù)飛機場的檢測結(jié)果，其中RX算法檢測效果不理想，由于背景光譜信息的混擾在檢測到的目標附近產(chǎn)生較多虛警目標，KRX算法由于加入了非線性核函數(shù)，效果有所改善，但在少量目標附近仍不能區(qū)分目標和背景，采用的MHD算法是從光譜匹配程度的角度出發(fā)，檢測目標和背景光譜信息差異，并將這種差異作為不相似度量，將不相似度超過閾值判斷為目標，否則為判斷為背景。

ROC曲線又稱為接收機工作特性曲線，是一種表述檢測概率Pd與虛警概率Pf之間變化關(guān)系的曲線，是檢測性能定量分析中常用的手段。

圖6　六種算法不同檢測效果對比Fig.6 Kurtosis of MHD compared with other method

圖7給出了RX、KRX和MHD三種算法的ROC曲線，一般要求虛警概率不能大于10%，從圖中可以看出，在相同的虛警概率下，MHD的檢測概率要遠遠高于RX和KRX算法，從而可知MHD算法性能遠遠優(yōu)于其他兩種算法。

圖7　ROC曲線對比Fig.7 ROC curves contrast

與模擬數(shù)據(jù)實驗相同，為了量化分析本文算法的優(yōu)越性，與表1對比，增加了目標個數(shù)虛警個數(shù)作為性能指標對上述算法進行比較分析，如表2所示。

表2　不同算法性能量化分析Table 2 The quantitative analysis of different algorithms （Real data）

由表2可以看出，與其他檢測算法相比，MHD算法在提升檢測性能的同時極大地減少了計算時間，這是因為，在RX和KRX中存在著矩陣求逆運算和復(fù)雜的矩陣操作導(dǎo)致耗時較長，而在MHD中僅僅需要集合運算，這大大降低了計算復(fù)雜度，節(jié)省了計算時間，可以大量應(yīng)用到具有較低計算能力的硬件實現(xiàn)中，例如FPGA、ARM、DSP等。上述實驗充分證明了MHD算法的有效性、通用性和實用性。

4　結(jié)論

針對現(xiàn)有高光譜異常檢測的RX算法易受背景干擾檢測性能不佳和計算復(fù)雜度高這一問題，本文將Hausdorff度量引入高光譜異常檢測，利用改進的Hausdorff距離（MHD）從光譜匹配程度的角度，進行了高光譜異常目標檢測，最大程度地將異常目標和背景分離，同時降低了算法的計算復(fù)雜度。實驗結(jié)果表明:

1）本文提出的MHD算法可以更好地檢測異常目標和背景的光譜相似度，對不同復(fù)雜背景下的多異常目標具有很好的魯棒性；

2）算法時間顯著降低，接近實時高光譜異常目標檢測，較低的計算復(fù)雜度為硬件實現(xiàn)提供了良好的算法支持。

對于極弱小異常目標的檢測還存在一些不足，下一步將針對異常窗口的選擇提取背景光譜信息進行深入研究。

［1］REED I S，YU X.Adaptive multiple-band CFAR detection of an optical pattern with unknown spectral distribution［J］. IEEE transactions on acoustics，speech，and signal processing，1990，38（10）:1760-1770.

［2］RILEY R A，NEWSOM R K，ANDREWS A K.Anomaly detection in noisy hyperspectral imagery［C］//Proceedings of the SPIE 5546，Imaging Spectrometry X.United States: SPIE，2004:159-170.

［3］王立國，趙春暉.高光譜圖像處理技術(shù)［M］.北京:國防工業(yè)出版社，2013:110-169. WANG Liguo，ZHAO Chunhui.Processing techniques of hyperspectral imagery［M］.Beijing:National Defend Industry Press，2013:110-169.

［4］CHANG C I，CHIANG S S.Anomaly detection and classification for hyperspectral imagery［J］.IEEE transactions on geoscience and remote sensing，2002，40（6）:1314-1325.

［5］MATTEOLI S，VERACINI T，DIANI M，et al.A locally adaptive background density estimator:an evolution for RX-based anomaly detectors［J］.IEEE geoscience and remote sensing letters，2014，11（1）:323-327.

［6］ZHAO Chunhui，WANG Yulei，QI Bin，et al.Global and local real-time anomaly detectors for hyperspectral remote sensing imagery［J］.Remote sensing，2015，7（4）:3966-3985.

［7］CHEN S Y，WANG Yulei，WU Chaocheng，et al.Real-time causal processing of anomaly detection for hyperspectral imagery［J］.IEEE transactions on aerospace and electronic systems，2014，50（2）:1511-1534.

［8］WANG Yulei，CHEN S Y，LIU Chunhong，et al.Background suppression issues in anomaly detection for hyperspectral imagery［C］//Proceedings of the SPIE 9124，Satellite Data Compression，Communications，and Processing X.United States:SPIE，2014.

［9］GOLDBERG H，KWON H，NASRABADI N M.Kernel eigenspace separation transform for subspace anomaly detection in hyperspectral imagery［J］.IEEE geoscience and remote sensing letters，2007，4（4）:581-585.

［10］王靖，朱夢宇，趙保軍，等.基于小波和改進型Hausdorff距離的遙感圖像配準方法［J］.電子學(xué)報，2006，34 （12）:2167-2169. WANG Jing，ZHU Mengyu，ZHAO Baojun，et al.A remote sensing image registration method based on wavelet decomposition and the improved Hausdorff distance［J］.Acta electronica sinica，2006，34（12）:2167-2169.

［11］楊清夙，游志勝，張先玉.基于豪斯多夫距離的快速多人臉檢測算法［J］.電子科技大學(xué)學(xué)報，2004，33（4）: 407-409. YANG Qingsu，YOU Zhisheng，ZHANG Xianyu.Fast multi-face detection algorithm based on Hausdorff distance ［J］.Journal of university of electronic science and technology of China，2004，33（4）:407-409.

［12］YI J H，BHANU B，LI Ming.Target indexing in SAR images using scattering centers and the Hausdorff distance［J］. Pattern recognition letters，1996，17（11）:1191-1198.

［13］何力.一種基于Hausdorff距離的運動目標跟蹤算法［D］.合肥:中國科學(xué)技術(shù)大學(xué)，2010. HE Li.An object tracking algorithm based on hausdorff distance［D］.Hefei:University of Science and Technology of China，2010.

［14］PEYVANDI H.Evolving radial basis function neural network with Hausdorff similarity measure for SONAR signals detection/classification［C］//Proceedings of the Oceans 2009-Europe.Bremen:IEEE，2009.

［15］GRANA M，VEGANZONES M A.An endmember-based distance for content based hyperspectral image retrieval［J］. Pattern recognition，2012，45（9）:3472-3489.

本文引用格式：

趙春暉，尤偉，齊濱，等.基于Hausdorff度量的高光譜異常目標檢測算法［J］.哈爾濱工程大學(xué)學(xué)報，2016，37（7）:979-985.

ZHAO Chunhui，YOU Wei，QI Bin，et al.Hyperspectral anomaly target detection algorithm based on Hausdorff measure［J］.Journal of Harbin Engineering University，2016，37（7）:979-985.

Hyperspectral anomaly target detection algorithm based on Hausdorff measure

ZHAO Chunhui1，YOU Wei1，QI Bin2，WANG Yulei1
（1.College of Information and Communication Engineering，Harbin Engineering University，Harbin 150001，China；2.College of Underwater Acoustic Engineering，Harbin Engineering University，Harbin 150001，China）

In anomaly target detection in hyperspectral imagery，it can be difficult to accurately distinguish between the spectral information of the targets and background，which leads to a decline in target detection performance.The results of the classic RX detection algorithm have a high false alarm probability，and the process is characterized by a large amount of calculation and complexity.To address these issues，we introduce the Hausdorff metric to hyperspectral anomaly target detection，prove the usefulness of its application，and make a number of improvements to suppress noise interference.In terms of the spectral matching，we separate the targets and background to a greater degree based on the improved Hausdorff distance.Experiments were performed using both synthetic and real hyperspectral data.Moreover，the results show improved detection performance and an increase in computational efficiency of nearly 60%.These experimental results prove that this algorithm has lower computational complexity and better performance than the traditional RX algorithm，casual RX algorithm，and KRX algorithm and can better support the implementation of hardware.

hyperspectral remote sensing；anomaly target detection；RX algorithm；spectral matching；Hausdorff metric

10.11990/jheu.201506087

TN911.73

A

1006-7043（2016）07-979-08

2015-06-26.網(wǎng)絡(luò)出版日期：2016-05-27.

國家自然科學(xué)基金項目（61571145，61405041）；黑龍江省自然科學(xué)基金項目 （ZD201216）；哈爾濱市優(yōu)秀學(xué)科帶頭人基金項目 （RC2013XK009003）；中國博士后基金項目（2014M551221）；中央高校基礎(chǔ)研究基金項目（HEUCF1608）.

趙春暉（1965-），男，教授，博士生導(dǎo)師，博士.

趙春暉，E-mail:zhaochunhui@hrbeu.edu.cn.

網(wǎng)絡(luò)出版地址：http：//www.cnki.net/kcms/detail/23.1390.u.20160527.1446.028.html