王海晏,江 濤,王 芳,鄧世杰
(空軍工程大學(xué),陜西 西安 710038)
高光譜成像技術(shù)是建立在成像光譜學(xué)基礎(chǔ)上的全新的綜合性遙感技術(shù),其集光譜分光技術(shù)、光電轉(zhuǎn)換技術(shù)、探測器技術(shù)、光學(xué)成像技術(shù)于一體[1]。其利用成像光譜儀納米級的光譜分辨率,在幾十或幾百個波段同時對地表地物成像,能夠?qū)崿F(xiàn)地物空間信息、輻射信息、光譜信息的同步獲取,因而在資源勘察,災(zāi)害檢測,海洋、陸地科學(xué)研究,軍事探測等方面有巨大的應(yīng)用價值和發(fā)展前景。
自光譜成像技術(shù)出現(xiàn)以來,目標(biāo)探測已經(jīng)發(fā)展了許多年。從光譜維的角度出發(fā)發(fā)展出了約束能量最小化算法(constrained energy minimizatio,CEM),自動色彩均衡算法(automatic color equalization,ACE),恒虛警檢測算法(reed-xiao detector,RXD),光譜角匹配算法(spectral angle mapper,SAM),光譜相關(guān)系數(shù)映射算法(spectral correlation coefficient mapping,SCM),光譜信息散度算法(spectral information divergence,SID)[2-4]等第一代探測算法,這些算法奠定了目標(biāo)探測的主基調(diào),利用信號處理原理,從光譜角度出發(fā)通過對背景信號的抑制從而突出目標(biāo)信號實現(xiàn)目標(biāo)探測。
在發(fā)展過程中又出現(xiàn)了一個新的分支——結(jié)合空間信息的探測算法。雖然高光譜空間分辨率受到了限制,但是空間信息還是有可挖掘之處。文獻(xiàn)[5]從4個方面論述了空間輔助信息對高光譜遙感數(shù)據(jù)分析的輔助作用。文獻(xiàn)[6]研究了鄰域信息對待測中心像元的影響,在圖像稀疏特性的基礎(chǔ)上將像元的鄰域信息應(yīng)用到分類中,提出了WJSRC算法。文獻(xiàn)[7]將光譜維和空間維結(jié)合成SAM_SURF算法,其中光譜特征識別所用的算法是SAM,空間特征識別所用的算法是SURF,最后顯示空譜聯(lián)合識別結(jié)果。
對空間信息的利用從本質(zhì)上看也是利用多角度信息減小不確定性,于是考慮將人的先驗知識用知識圖譜的方式在識別過程中加以利用。文獻(xiàn)[8]在雷達(dá)目標(biāo)識別應(yīng)用領(lǐng)域,為提升雷達(dá)目標(biāo)識別精度,構(gòu)建了一種基于目標(biāo)特性知識圖譜,主要描述面向雷達(dá)目標(biāo)識別領(lǐng)域的行業(yè)知識圖譜技術(shù)。文獻(xiàn)[9]提出了以知識為中心的地物目標(biāo)識別,這個框架能夠更好地利用知識,能夠在不同的算法中最大程度地實現(xiàn)知識共享。
雖然利用先驗知識進(jìn)行目標(biāo)識別已經(jīng)有了一定的發(fā)展,但是在高光譜目標(biāo)識別中,依然沒有具體模型,因此本文針對知識圖譜對高光譜數(shù)據(jù)進(jìn)行目標(biāo)識別存在的問題,提出一種基于融合光譜、空間和人先驗知識的識別模型。
由人的視覺處理系統(tǒng)可以知道,僅從視覺角度可以通過物體的顏色、形狀輪廓、運動和深度等信息再結(jié)合主觀經(jīng)驗進(jìn)行目標(biāo)識別。而基于人的主觀經(jīng)驗識別出的目標(biāo)特征是指人在實踐過程中,逐步抽象出來的更高級的物體的屬性。以飛機目標(biāo)為例,屬性主要包含外形輪廓、顏色、運動等。此外在高光譜圖像中,數(shù)據(jù)是一個三維立方體,包含了二維平面信息和第三維光譜信息,因此屬性中還應(yīng)當(dāng)包含光譜維信息。
要理解關(guān)聯(lián)關(guān)系,首先看人對目標(biāo)的識別過程。圖1是普通高光譜圖像的目標(biāo),數(shù)據(jù)是ABU數(shù)據(jù)集中Urban中的一幅圖像,目標(biāo)是船只,由美國的機載可視紅外成像光譜儀(AVIRIS)所拍攝。
圖1 示例圖像Fig.1 Example image
僅根據(jù)圖1(a)很難對圖中的目標(biāo)做出判斷,然而對照圖1(b)則能夠知道圖像中可能是船,或者說在人的大腦認(rèn)知中認(rèn)為是船的概率比較高。這和大腦中存儲的知識有關(guān),進(jìn)一步說就是利用了知識之間的關(guān)聯(lián)性。人們的認(rèn)知就是船是在水里的,于是問題可以變成對水體的尋找,水和船變成了關(guān)聯(lián)較大的一組概念,在高光譜圖像空間維中往往蘊含豐富的信息,如果能充分利用其中信息,再結(jié)合人的主觀經(jīng)驗,使得最終決策更高效。
關(guān)聯(lián)關(guān)系包括目標(biāo)與目標(biāo)之間以及目標(biāo)與地區(qū)之間的主要關(guān)聯(lián)關(guān)系類型等。這個關(guān)聯(lián)過程主要通過圖連接來實現(xiàn),如果實體之間存在直接連接,則說明其關(guān)聯(lián)關(guān)系緊密,反之則不然。例如在不同飛機目標(biāo)之間、飛機和跑道、飛機和航站樓、飛機和停機坪存在關(guān)聯(lián)關(guān)系。
為了尋找邊緣,在數(shù)字圖像處理中,使用梯度的概念來完成。一幅圖像f在位置(x,y)處的梯度用?f來表示,?f表達(dá)式為
(1)
式(1)中,?f是一個向量,大小用M(x,y)表示,其值為
(2)
例如Sobel、Canny等檢測算子都是在此基礎(chǔ)上發(fā)展起來的。
Sobel算子是一種用于邊緣檢測的離散微分算子,它結(jié)合了高斯平滑和微分求導(dǎo)。該算子用于計算圖像明暗程度近似值,根據(jù)圖像邊緣旁邊明暗程度把該區(qū)域內(nèi)超過某個數(shù)的特定點記為邊緣。其在x和y方向上的梯度計算方法為
(3)
(4)
對外形輪廓進(jìn)行模板匹配時,有以下幾點要求或問題:1) 匹配算法要簡單,準(zhǔn)確率要較高。由于還有其他識別步驟,該過程使用復(fù)雜的算法,會降低整個識別過程的實時性。2) 匹配過程不能受目標(biāo)位置平移的影響。由于目標(biāo)在圖像中的位置未知,因此要求匹配算法不能對目標(biāo)位置變化敏感。3) 匹配過程需要減小分辨率不對應(yīng)造成的影響。實際運用中大部分圖像都可能是異源圖像,在匹配過程中要求算法能在不同分辨率之間圖像進(jìn)行匹配。4) 算法還需要解決目標(biāo)旋轉(zhuǎn)問題。目標(biāo)具有顯著的輪廓特征,那么隨著目標(biāo)的旋轉(zhuǎn),算法結(jié)果差異很大,也會造成最終漏判誤判。針對以上要求及問題,設(shè)計了一種模板匹配算法——特征點匹配算法(feature point matching algorithm)。
步驟1 求得目標(biāo)區(qū)域和模板圖像的中心,設(shè)(x,y)為非零像元坐標(biāo),則其中心(x0,y0)計算公式為
(5)
步驟2 選定特征點位置。以中心(x0,y0)作為坐標(biāo)原點,在坐標(biāo)軸上C2距離處分別取4個特征點。設(shè)特征點坐標(biāo)為(xti,yti),則
(xti,yti)∈{(x0-C2,y0),(x0,y0-C2),
(x0+C2,y0),(x0,y0+C2)},
(6)
式(6)中,C2=min{x0,y0,b},b為設(shè)定的距離常數(shù)。
步驟3 以特征點為基礎(chǔ),計算非零像元與特征點之間的歐式距離和。設(shè)特征點(xti,yti)對應(yīng)的特征值為ti,則
(7)
步驟4 4個特征值組成目標(biāo)圖像特征向量Ttar= (t1,t2,t3,t4),利用上述過程同樣可以獲得模板圖像的特征向量Tsam,設(shè)兩個向量之間的相似度為S,則
(8)
根據(jù)cos-1函數(shù)值的特點S距離越接近0,表示目標(biāo)和模板越相似。
光譜角匹配是一種基本的目標(biāo)識別算法,它通過計算光譜向量之間的夾角來衡量待識別像素光譜與已知目標(biāo)光譜之間的相似程度,從而進(jìn)行判定該像素是否包含目標(biāo)。其數(shù)學(xué)表達(dá)式為
(9)
式(9)中,X是像元向量,Y是目標(biāo)向量。Dsam值越接近0則兩個光譜向量越相似。
實驗數(shù)據(jù)是美國加州圣地亞哥North Island海軍機場高光譜圖像,由AVIRIS機載成像光譜儀拍攝。原始圖像大小為 100 × 100,空間分辨率為3.5 m。在人們認(rèn)知中對飛機最敏感的是其外形輪廓,因此利用光譜和外形輪廓屬性進(jìn)行識別。圖像模板如圖2所示。
首先利用光譜角匹配方法(spectral angle mapper, SAM)利用參考光譜(參考光譜用目標(biāo)像元平均值代替)進(jìn)行匹配,然后對匹配結(jié)果利用輪廓識別進(jìn)行背景剔除,以10×10大小窗口進(jìn)行識別,在圖中會出現(xiàn)5個待識別目標(biāo),如圖3所示。
圖2 圖像模板Fig.2 Image template
圖3 待識別目標(biāo)Fig.3 Target to be identified
以圖2中模板為基礎(chǔ)對圖3中目標(biāo)和背景進(jìn)行識別,設(shè)置C1=2,由式(5)—(7)可得各圖像的中心和特征向量,模板 0、目標(biāo)1和背景1的中心和特征向量如表1所示。
表1 圖像中心和特征向量Tab.1 Image center and feature vector
由式(8)得目標(biāo)1和模板0的相似度為0.200,背景1和模板0的相似度為0.397,同理可得各目標(biāo)、背景圖像和模板圖像的匹配結(jié)果如表2。
表2 C1=2模板匹配結(jié)果Tab.2 Template matching results when C1=2
從表2中可以看出目標(biāo)與模板相似度值明顯低于背景與模板的相似度,從而很容易將目標(biāo)和模板分開。同時在目標(biāo)與模板中的最小值總出現(xiàn)在模板7附近,可以判斷目標(biāo)與模板7最為接近。從圖3中目標(biāo)輪廓可以看出匹配結(jié)果與實際情況相符。為了更直觀體現(xiàn)以表2數(shù)據(jù)為基礎(chǔ),畫出目標(biāo)相似度雷達(dá)圖如圖4。
圖4 相似度雷達(dá)圖Fig.4 Similarity radar image
由圖4可以看出目標(biāo)1,2,3在模板7附近雷達(dá)圖明顯有內(nèi)凹的趨勢,說明目標(biāo)輪廓與模板7附近最為相似,可以判斷目標(biāo)朝向與模板7朝向相似。驗證了提出的特征點匹配方法對物體識別的有效性以及對方位判別的能力。最終識別檢測概率和虛警率如表3所示,其中檢測概率和虛警率分別為
Pd=Nt/St,
Pf=Nb/Sb,
(10)
式(10)中,Nt為正確像元數(shù);St為目標(biāo)像元數(shù);Nb為錯誤識別像元數(shù);Sb為背景像元總數(shù),如圖3(a) 所示,Sb=100×100。
表3 檢測概率和虛警率Tab.3 Detection probability and false alarm rate
由此可知輪廓識別后可以在相同檢測概率的條件下,有效降低其虛警率,保證目標(biāo)識別的準(zhǔn)確性。
利用圖1中的ABU數(shù)據(jù)集,目標(biāo)是要尋找船只,通過知識圖譜之間的關(guān)聯(lián)關(guān)系可以知道,船只和水體的關(guān)聯(lián)程度較高,因此轉(zhuǎn)而尋找水體。由于水體的光譜特征比較顯著因此可在圖中優(yōu)先識別水體,以水的光譜作為匹配光譜識別圖中水體區(qū)域,在水體區(qū)域內(nèi)再進(jìn)行二次識別,結(jié)果如圖5所示。
由圖5(b)、(c)可以看出通過水體識別的方式能夠有效減少水體之外的地物的干擾,提高識別的準(zhǔn)確率。圖5(d)是ROC曲線,其橫軸表示虛警概率,其縱軸表示檢測概率,根據(jù)不同門限值可以得到一組虛警和檢測概率的坐標(biāo),進(jìn)而畫出接收機特性曲線(ROC)[10-11]。從曲線的分布可以看出,基于關(guān)聯(lián)關(guān)系的方法其ROC曲線始終在光譜角匹配識別的曲線上方,說明在同樣虛警率的情況下,基于關(guān)聯(lián)關(guān)系的方法檢測率更高,或者說在同樣檢測概率的情況下,基于關(guān)聯(lián)關(guān)系方法的虛警率更低。進(jìn)而說明了基于知識圖譜方法對提高檢測效率的有效性。
圖5 識別效果Fig.5 Recognition effect
本文在充分考慮人的先驗知識、空間信息和光譜信息的基礎(chǔ)上,引入知識圖譜作為識別框架,充分利用形狀輪廓以及地物之間關(guān)聯(lián)關(guān)系等信息,結(jié)合光譜匹配的方法進(jìn)行地物識別,在測試數(shù)據(jù)集上表現(xiàn)良好,能有效區(qū)分地物輪廓并提高檢測效率,表明該方法具有一定的可實踐性。在后續(xù)的研究中將擴展目標(biāo)屬性,構(gòu)建更加全面的知識圖譜和光譜庫以適應(yīng)更廣泛、更復(fù)雜的地物。