琚麗君,田豐華,曾朝平

(1.河南測繪職業(yè)學(xué)院,鄭州 451464;2.鄭州華潤燃?xì)夤煞萦邢薰?鄭州 451464)

0 引言

隨著對地觀測技術(shù)的快速發(fā)展,影像的空間分辨率達(dá)到亞米級別。高分辨率遙感影像具有信息豐富、邊界輪廓清晰、光譜波段多等優(yōu)點,該優(yōu)點擴(kuò)展了遙感影像的應(yīng)用領(lǐng)域,但也增加了高分辨率遙感影像中有效信息的提取難度[1]。影像分割通過影像光譜、顏色以及紋理等信息,將整景影像分割成符合人眼視覺要求的多個區(qū)域,以增強(qiáng)區(qū)域內(nèi)特征的一致性,同時突出不同區(qū)域差異性[2-3]。

傳統(tǒng)的基于像元的分割方法未考慮影像紋理細(xì)節(jié)信息,僅依靠影像光譜信息導(dǎo)致分割結(jié)果存在很大的不確定性,加上影像存在噪聲,分割結(jié)果出現(xiàn)“椒鹽噪聲”現(xiàn)象,給后續(xù)處理帶來新的問題[4-5]?；谏疃葘W(xué)習(xí)的影像分割方法在樣本數(shù)據(jù)充足的情況下能獲取良好的分割結(jié)果,但不同季節(jié)、不同緯度區(qū)域、不同傳感器獲取影像分割效果不理想,而且基于深度學(xué)習(xí)的分割方法獲取地物分割邊界需要后期人工編輯[6-7]。在樣本數(shù)據(jù)不足情況下,為了獲取較為理想的分割結(jié)果,相關(guān)學(xué)者利用影像空間幾何、紋理特征以及光譜信息關(guān)系,將超像素分割方法[8-9]應(yīng)用于遙感影像分割中,該類方法以聚類為基本思想,將影像分割成接近均勻、互斥的超像素塊,提高后期影像處理效率[10]。董志鵬等[11]提出簡單線性迭代聚類與最小生成樹融合的高分辨率遙感影像分割方法,該方法通過區(qū)域動態(tài)約束策略對超像素分割結(jié)果進(jìn)行合并,獲取良好的分割結(jié)果;周佳超等[12]針對簡單線性迭代聚類(simple linear iterative cluster,SLIC)算法對遙感影像分割準(zhǔn)確率、效率不足等問題,采用Wasserstein距離優(yōu)化SLIC方法,在確保分割準(zhǔn)確性的前提下提升算法運算效率;張銳等[13]提出標(biāo)記分水嶺與改進(jìn)超像素融合的遙感影像分割方法,該方法在超像素分割時充分應(yīng)用影像顏色、空間、紋理等多種特征,對于過分割現(xiàn)象通過超像素合并的方法獲取最終結(jié)果。上述可知,由于遙感影像場景相比常規(guī)影像復(fù)雜得多,將SLIC超像素分割算法直接應(yīng)用于高分辨率遙感影像分割鄰域時,較多的微小特征與紋理會導(dǎo)致超像素塊分割獨立地物現(xiàn)象,分割性能較差。李嬌嬌等[14]提出深度網(wǎng)絡(luò)模型的遙感影像分割方法,但在遙感場景復(fù)雜、目標(biāo)地物的尺度不一致時,會導(dǎo)致算法處理效率低而且一景影像上不同地物的分割精度差別較大,因此需要較多樣本數(shù)據(jù)對模型進(jìn)行優(yōu)化。線性譜聚類(linear spectral cluster,LSC)[15]算法是一種較為新穎的超像素分割方法,其分割性能較高。

鑒于此,文中提出歸一化轉(zhuǎn)動慣量特征(normalized moment of inertia,NMI)優(yōu)化邊界敏感的LSC遙感影像分割算法。首先,為了增強(qiáng)影像邊緣特征,在LSC初始化聚類中心時,對包含地物邊界初始超像素進(jìn)行邊界約束,更新初始像素塊;然后,利用LSC算法對高分辨率遙感影像進(jìn)行分割,獲取候選超像素結(jié)果,并統(tǒng)計微小超像素塊;最后,計算微小超像素塊與相鄰超像素塊的歸一化轉(zhuǎn)動慣性量相似性度量值,并將微小超像素塊合并到相似性最小的超像素塊上,有效避免分割出微小超像素塊。

1 改進(jìn)的LSC超像素分割算法

1.1 邊緣敏感的線性譜聚類分割算法

線性譜聚類算法通過設(shè)計核函數(shù)將影像投影到高維特征空間,進(jìn)行簡單加權(quán)的K-means聚類數(shù)據(jù)處理,較大程度地降低算法復(fù)雜程度,提高算法效率。LSC算法將原始RGB空間影像數(shù)據(jù)轉(zhuǎn)換到高維特征空間,每個像素采用所在位置坐標(biāo)(x,y)和顏色值(l,α,β)共同構(gòu)成一個五維的特征相鄰表示。其中,l表示顏色亮度,α表示深綠色到亮粉紅色顏色通道值,β表示藍(lán)色到黃色顏色通道值。LSC生成超像素的具體實現(xiàn)過程如下。

1)將影像上任意像素點p,通過φ(p)核函數(shù)表示為多維特征向量,φ(p)核函數(shù)的表達(dá)如式(1)所示。

(1)

式中:(lp、αp、βp、xp、yp)、(lq、αq、βq、xq、yq)表示像素點p和像素點q的五維向量;V表示待分割影像數(shù)據(jù)集;p、q表示影像上任意兩個像素點;w(p)表示像素點p與像素點q在高維特征空間的距離;Cc、Cs表示調(diào)整顏色與空間權(quán)重的參數(shù)因子,Cc=Cs/r,當(dāng)r取值較小時,顏色空間相似的像素聚集,反之,空間相近的像素聚集,文中r=0.075。

2)在影像行列方向分別設(shè)置水平間隔vx、垂直間隔vy,生成K個初始化超像素種子點。

3)設(shè)置迭代次數(shù)T、初始化聚類中心ck與類內(nèi)加權(quán)均值mk,表達(dá)如式(2)所示。

(2)

式中:Mk為第k個超像素數(shù)據(jù)集,k=1,2,3,…,K。

4)統(tǒng)計ck類中心局部鄰域(τvx×τvy,τ>0.5)內(nèi)像素點p和mk的相似性值,并將像素點p更新到差異最小值類上,表達(dá)如式(3)所示。

(3)

5)更新式(2)獲取第k個超像素的加權(quán)均值mk,直至到達(dá)迭代次數(shù)或超像素加權(quán)均值mk收斂,否則迭代步驟1)～步驟3)。

在LSC初始化聚類中心時,存在包含地物邊界的跨邊界初始超像素塊。若直接進(jìn)行LSC處理,生成超像素塊的邊界難以有效地貼合地物邊界導(dǎo)致部分地物邊界信息丟失。跨邊界超像素塊如圖1所示。

圖1 跨邊界超像素塊

為了獲取地物邊界完整的超像素塊,提出邊緣敏感的LSC高分辨率遙感影像分割算法。該算法通過邊緣檢測算子提取地物邊界信息,結(jié)合初始超像素塊判斷是否為跨邊界初始超像素塊,若是跨邊界超像素塊,需用跨邊界初始超像素的區(qū)域質(zhì)心更新原始聚類中心。新的聚類中心確定規(guī)則如下。

1)若初始聚類中心局部S×S區(qū)域內(nèi)存在地物邊界,剔除初始聚類中心,并將地物邊界與超像素塊邊界閉合,將區(qū)域質(zhì)心更新為處理聚類中心,如圖2中Bi所示。

圖2 聚類中心更新示意圖

2)若初始聚類中心局部S×S區(qū)域內(nèi)不存在地物邊界,保留初始聚類中心不變,如圖2中Ci所示。

本文選擇LoG邊緣檢測算子對原始影像進(jìn)行邊緣檢測,該算子能夠有效地克服噪聲對提取結(jié)果的影響。然后,將邊界檢測結(jié)果與LSC分割產(chǎn)生的初始超像素塊融合,更新超像素塊聚類中心生成超像素分割結(jié)果。通過引入邊緣信息優(yōu)化LSC算法,可以充分利用高分辨率的地物邊緣細(xì)節(jié)信息,增強(qiáng)邊緣信息對LSC算法的邊界約束力,有效區(qū)分影像弱邊緣與噪聲信息,確保分割出的地物信息更加完整、吻合。

1.2 基于NMI特征的超像素合并

由上文分析可知,為避免分割結(jié)果中出現(xiàn)欠分割情況,超像素塊數(shù)K不能設(shè)置過小;但隨著超像素數(shù)K值的增加,分割結(jié)果也會存在過分割現(xiàn)象。針對邊緣敏感的LSC分割產(chǎn)生的微小超像素塊,本文通過計算影像微小超像素與鄰域超像素塊的繞質(zhì)心的歸一化轉(zhuǎn)動慣量(normalized moment of inertia,NMI)相似性,將其合并到鄰近超像素塊上。

對光譜紋理相似的區(qū)域,僅利用LSC算法計算影像局部區(qū)域的光譜與空間特征的相似性度量,難以有效進(jìn)行區(qū)分。若引入?yún)^(qū)域的整體特征,將提高超像素塊的分割精度。NMI特征是影像二階矩特征,通過影像空間信息表示影像的光譜復(fù)雜分布程度,具有平移、旋轉(zhuǎn)和縮放不變性的特點,已在影像匹配中有所應(yīng)用。本文將該特征作為計算超像素塊間的相似性值。對于一幅M×N大小的遙感影像,每個像素點可以視為平面上的一個質(zhì)點,每個波段的灰度值I(i,j)可以表示為質(zhì)點的質(zhì)量,影像所有像素灰度之和為影像質(zhì)量m,表達(dá)如式(4)所示。

(4)

影像質(zhì)心(i0,j0)如式(5)所示。

(5)

假如影像灰度分布比較均勻,那么影像的質(zhì)心應(yīng)該為圖形的平面形心。因此,影像繞質(zhì)心(i0,j0)的轉(zhuǎn)動慣量可以表示為J(i0,j0),如式(6)所示。

(6)

依據(jù)影像質(zhì)量、影像質(zhì)心以及轉(zhuǎn)動慣量定義歸一化轉(zhuǎn)動慣量,表達(dá)如式(7)所示。

(7)

由式(7)可知,歸一化轉(zhuǎn)動慣量與影像灰度值分布、質(zhì)心位置和灰度分布相對質(zhì)心的分布均有相關(guān)。因此,可利用歸一化轉(zhuǎn)動慣量區(qū)別紋理差異。

由于LSC算法分割影像生成超像素存在微小像素塊,因此需計算微小像素塊與其鄰域的像素塊的NMI值,并對微小超像素進(jìn)行合并處理。

LSC算法分割影像生成的超像素用S表示,表達(dá)如式(8)所示。

(8)

對于一幅M×N的遙感影像,通過LSC算法分割生成K個超像素,由文獻(xiàn)[16]可知,當(dāng)超像素塊滿足N′

1.3 算法流程

文中采用NMI特征的相似性度量方法合并微小超像素塊,提高影像分割精度,流程圖如圖3所示。首先,通過LOG算子提取影像邊界信息,并以邊界信息為引導(dǎo),結(jié)合初始超像素塊判斷是否為跨邊界初始超像素塊,若是跨邊界超像素塊,需用跨邊界初始超像素的區(qū)域質(zhì)心更新原始聚類中心,避免初始超像素塊內(nèi)存在地物邊界信息;然后,利用邊緣敏感的LSC算法分割高分辨率遙感影像生成候選超像素塊,并統(tǒng)計微小超像素塊;最后,為了避免超像素過分割結(jié)果存在,通過計算微小超像素塊與相鄰超像素塊的歸一化轉(zhuǎn)動慣性量相似性度量值,將微小超像素塊合并到相似性最小的超像素塊上,有助于后期影像處理效率。

圖3 算法流程圖

算法具體流程如下。

1)初始化超像素塊數(shù)K,計算初始聚類中心。

2)利用LOG算子提取影像邊界信息,生成邊界圖像。

3)以邊界影像為引導(dǎo),更新初始聚類中心。

4)獲取影像候選超像素分割結(jié)果,并標(biāo)記微小超像素塊。

5)計算微小超像素塊與局部鄰域超像素塊NMI特征相似性值。

6)將微小超像素與NMI特征差異最小的超像素塊合并。

7)影像分割完成。

2 實驗與分析

2.1 算法流程

為了驗證本文方法的有效性、可行性,文中對兩組不同場景的高分辨率遙感影像數(shù)據(jù)進(jìn)行分割處理,并設(shè)計兩種方法與本文方法對比:一是LSC算法,二是eCognition軟件中的多尺度分割方法。本文方法的實驗運行系統(tǒng)環(huán)境為Windows 10系統(tǒng),CPU為Core(TM) i5-124900f,內(nèi)存為32 GB,在實驗平臺MATLAB2016、eCognition Developer 9.0軟件上進(jìn)行仿真實驗。

選用定性與定量方法進(jìn)行分割結(jié)果精度評定。其中,定性方法通過視覺對比方法分析分割結(jié)果;定量方法通過查全率R、查準(zhǔn)率P以及R、P的平均值F1等指標(biāo)對分割結(jié)果進(jìn)行精度評價。定義如式(9)至式(11)所示。

(9)

(10)

(11)

式中:TP表示獨立地物分割屬于相同對象的數(shù)目;FN表示獨立地物被分割為其他地物的數(shù)目;FP表示獨立地物分割中出現(xiàn)其他對象的數(shù)目。其中,R、P值越大,表示分割精度越高。當(dāng)R與P難以描述分割結(jié)果時,通過F1進(jìn)行綜合對比分析。

2.2 實驗結(jié)果

1)實驗1結(jié)果與分析。選用高景一號01星全色影像與多光譜影像經(jīng)過最鄰近融合(nearest neighbor diffuse,NND)生成的多光譜數(shù)據(jù),影像包含近紅外、紅、綠以及藍(lán)4個波段,實驗輸入數(shù)據(jù)為紅、綠、藍(lán)3個波段,影像大小為1 200像素×1 200像素,影像分辨率為0.5 m。該區(qū)域包含水體、植被、裸地等典型地物信息,如圖4(a)所示。由圖4(b)可知,由于高分辨率遙感影像中的地物繁雜,采用LOG邊緣檢測算子在原始遙感影像中進(jìn)行邊緣檢測,會出現(xiàn)諸如模糊邊緣紋理、可能產(chǎn)生虛假邊界和對曲線邊界定位不準(zhǔn)確等問題,導(dǎo)致一定程度的分割性能降低,但邊緣敏感的線性譜聚類分割算法僅在初始化聚類中心時對存在邊界信息的初始像素塊進(jìn)行區(qū)域聚類中心更新,在后期NMI特征超像素合并時會對上一步結(jié)果進(jìn)行優(yōu)化。

圖4 影像數(shù)據(jù)

其中,圖5(a)為LSC算法分割高分辨率遙感影像結(jié)果。由圖5(a)可知,整個高分辨率遙感影像的分割結(jié)果存在嚴(yán)重的過分割現(xiàn)象,特別是植被和水體區(qū)域分割過于破碎,且存在微小的超像素塊區(qū)域。圖5(b)、圖5(c)分別為多尺度分割方法下不同尺度因子產(chǎn)生的分割結(jié)果,由于圖5(b)分割尺度因子為60,在水體內(nèi)部存在光譜變化較大的區(qū)域產(chǎn)生大量過分割結(jié)果。由圖5(c)可知,當(dāng)尺度因子增加到100時,過分割現(xiàn)象有所緩解,但在植被、水體區(qū)域的過分割現(xiàn)象未能消除,同時隨著尺度因子的增加,部分區(qū)域出現(xiàn)了欠分割現(xiàn)象,見圖5(c)橢圓形區(qū)域。由于文獻(xiàn)[12]采用Wasserstein距離優(yōu)化K-means算法中距離度量工具,難以有效區(qū)分耕地區(qū)域間的差異性,導(dǎo)致圖5(d)耕地區(qū)域存在過分割問題。由圖5(e)可知,本文方法有效解決水體、植被區(qū)域存在的過分割問題,通過邊緣特征優(yōu)化初始LSC聚類中心,避免初始像素塊穿過物體邊界,同時通過引入NMI特征相似度參數(shù),有效合并微小超像素塊獲取的分割結(jié)果基本與地物邊界一致,較好地保留了植被、水體、道路等區(qū)域邊界信息。

圖5 實驗數(shù)據(jù)及分割結(jié)果

由圖5(a)可知,LSC算法分割結(jié)果存在嚴(yán)重的過分割現(xiàn)象,導(dǎo)致分割精度很差。因此,本文不對LSC算法結(jié)果進(jìn)行精度評定。表1給出多尺度分割方法與本文方法的分割結(jié)果精度評定。由表1可知:本文分割結(jié)果的查準(zhǔn)率P為0.916,均高于多尺度分割方法(0.815、0.853),說明本文方法的分割結(jié)果優(yōu)于多尺度分割方法;本文方法分割結(jié)果的查全率R優(yōu)于0.9,說明本文分割結(jié)果與地物實際邊界基本保持一致;本文方法相比較不同尺度參數(shù)的多尺度分割方法平均值F1分別提高0.091、0.059,說明本文方法通過邊緣特征優(yōu)化LSC初始聚類中心,能有效避免后期分割超像素塊通過線性特征區(qū)域,通過引入NMI特征計算不同超像素塊間的相似性值,將微小超像素塊合并到相似性區(qū)域,避免過分割現(xiàn)象,進(jìn)一步提升影像分割精度。

表1 實驗1分割結(jié)果精度評定

圖6給出超像素數(shù)量K=[200,400,600,800,1 000,1 200,1 400]時生成的影像分割精度。由圖6可知,當(dāng)K取1 000時,算法分割精度的查全率最高,說明算法分割影像邊界與實際地物貼近,同一地物基本被分割到一個對象區(qū)域,有效避免過分割現(xiàn)象,分割結(jié)果更加符合人眼視覺效果;隨著K值增加,查全率與準(zhǔn)確率有所下降,因此邊緣貼合程度有所下降。

圖6 實驗1影像超像素參數(shù)初始值對分割精度影響

2)實驗2結(jié)果與分析。實驗2為吉林一號高分04A星數(shù)據(jù)全色影像與多光譜影像經(jīng)過NND融合生成的多光譜數(shù)據(jù),影像包含近紅外、紅、綠以及藍(lán)4個波段,實驗輸入數(shù)據(jù)為紅、綠、藍(lán)3個波段,影像大小為1 000像素×1 000像素,影像分辨率為0.3 m。該區(qū)域包含建筑物、道路、植被等典型地物信息,如圖7(a)所示,邊緣檢測結(jié)果如圖7(b)所示。

圖7 影像數(shù)據(jù)

圖8(a)為LSC方法過分割結(jié)果,由圖8(a)可知,該結(jié)果依然存在嚴(yán)重過分割現(xiàn)象。圖8(b)、圖8(c)為多尺度分割方法分割結(jié)果,圖8(b)為尺度因子70時分割結(jié)果,圖上的矩形是獨立建筑物被分割為多個對象,由于建筑物光譜特征差別較大,導(dǎo)致過分割現(xiàn)象。隨著尺度因子的增大,圖8(c)為尺度因子為100時的分割結(jié)果,由圖8(c)可知,該區(qū)域建筑物的過分割現(xiàn)象明顯好轉(zhuǎn),但隨著尺度因子的增加,圖上橢圓區(qū)域存在明顯的欠分割。由圖8(d)可知,由于建筑物頂部存在明顯彩鋼瓦拼接縫,導(dǎo)致建筑物頂部被分割為多個對象,通過采用改進(jìn)的距離相似度計算方法難以將該區(qū)合并,導(dǎo)致過分割問題。由圖8(e)可知,本文方法可以很好地區(qū)分建筑物、道路、植被等地物,較好地保留地物的邊緣信息,避免過分割與欠分割情況。

圖8 實驗數(shù)據(jù)及分割結(jié)果

表2為多尺度分割方法與本文方法的分割結(jié)果精度評定。由表2可知:本文分割結(jié)果的查準(zhǔn)率P為0.932,均高于多尺度分割方法(0.854、0.891),說明本文方法的分割結(jié)果優(yōu)于多尺度分割方法;本文方法分割結(jié)果的查全率R為0.966,說明本文分割結(jié)果與地物實際邊界基本保持一致;本文方法相比較不同尺度參數(shù)的多尺度分割方法平均值F1分別提高0.073、0.047,說明本文方法能有效避免后期分割超像素塊通過線性特征區(qū)域,分割結(jié)果更加合理有效。

表2 實驗2分割結(jié)果精度評定

圖9給出超像素數(shù)量K=[200,400,600,800,1 000,1 200,1 400]時生成的影像分割精度。由圖9可知,當(dāng)K取1 200時,算法分割影像查全率與查準(zhǔn)率最高,說明影像分割邊界與實際地物貼近,有效避免過分割現(xiàn)象,分割結(jié)果更加符合人眼視覺效果。

圖9 實驗2影像超像素參數(shù)初始值對分割精度影響

3 結(jié)束語

針對線性譜聚類方法處理復(fù)雜場景的高分辨率遙感影像時存在地物邊界丟失、過分割問題,本文提出NMI特征優(yōu)化邊緣敏感的LSC高分辨率遙感影像分割方法。利用邊緣信息引導(dǎo)LSC初始聚類中心,能有效避免初始超像素塊內(nèi)部存在線性地物,減少過分割問題。同時,結(jié)合NMI特征計算微小超像素塊與相鄰超像素塊相似性,合并微小超像素塊進(jìn)一步改善過分割問題。利用兩組高分辨率遙感影像數(shù)據(jù)進(jìn)行實驗,實驗結(jié)果表明,該方法可有效改善復(fù)雜區(qū)域的過分割現(xiàn)象,提升影像分割效果。在實驗過程中發(fā)現(xiàn),高分辨率遙感影像中的地物繁雜,在遙感影像中使用LOG邊緣檢測算子會存在模糊邊緣紋理、虛假邊界等問題,導(dǎo)致分割性能降低,這也是算法的下一步改進(jìn)方向。