李洪川，楚 恒,2,3，霍英海

(1. 重慶郵電大學(xué)光通信與網(wǎng)絡(luò)重點實驗室，重慶 400065； 2. 西南大學(xué)地理科學(xué)學(xué)院，重慶 400715； 3. 重慶市勘測院，重慶 400020)

城市道路作為高分辨率遙感影像中重要的地物信息，其精確提取在地理國情普查與監(jiān)測、城市規(guī)劃、交通導(dǎo)航等領(lǐng)域應(yīng)用中具有重要意義。由于高分辨率遙感影像的空間分辨率較中低分辨率影像顯著提高，地物細節(jié)信息更加豐富，紋理更加清晰，“同物異譜”“同譜異物”現(xiàn)象廣泛存在，道路提取更加困難[1]。如何從擁有復(fù)雜地物的高分辨率遙感影像中精確提取道路信息成為目前研究的熱點與難點。

目前常用的道路提取方法有模板匹配法[2]、Snakes模型方法[3]、邊緣檢測法[4]、面向?qū)ο蠓╗5]等。其中面向?qū)ο蟮牡缆诽崛》椒ㄔ诳朔敖符}”噪聲的同時能夠利用影像豐富的空間信息，提高道路提取精度，成為目前道路提取首選的方法。SVM在解決小樣本、非線性及高維數(shù)據(jù)的模式識別問題具有特有的優(yōu)勢，被越來越多的國內(nèi)外學(xué)者應(yīng)用到面向?qū)ο蟮牡缆诽崛≈?。ZHU等提出了一種改進SVM的道路提取方法。該方法通過聚類減少SVM分類數(shù)目，然后再訓(xùn)練SVM實現(xiàn)道路提取，取得了不錯的效果[6]。MIAO等提出了一種面向?qū)ο蟮腟VM道路提取方法。該方法在對象級影像上利用SVM并結(jié)合張量投票實現(xiàn)道路提取[7]。林鵬等提出一種面向?qū)ο蟮某擎?zhèn)道路提取方法。該方法利用SVM實現(xiàn)了農(nóng)村道路的精確提取[5]。以上方法均是以單核SVM對影像進行提取，由于不同核函數(shù)的SVM分類器道路提取精度不同，同時，影像中不同的特征對不同核函數(shù)的適應(yīng)性也不相同，造成道路提取的精度還有待提高。

基于此，本文提出一種高分辨率影像城市道路的多特征多核SVM分層提取方法。首先，基于FCM算法將原始輸入高分辨率遙感影像粗分類為包含道路、建筑物等地物信息的建成區(qū)和非建成區(qū)，去除非建成區(qū)；運用分水嶺分割算法對建成區(qū)影像進行分割，同時提取分割后對象的光譜特征和空間特征；然后選取RBF核函數(shù)和POLY核函數(shù)加權(quán)組合成MKSVM并對建成區(qū)進行二次分類，獲取道路初始信息；最后經(jīng)過數(shù)學(xué)形態(tài)學(xué)處理，實現(xiàn)道路的最終提取。通過重慶市某地區(qū)的QuickBird影像試驗證明：該方法能夠較好地提取出城市道路信息，且道路的分類精度明顯高于單一核函數(shù)的SVM和其他對比方法。

1 道路提取過程

試驗所選影像地物信息復(fù)雜，既有水體、植被，又有與道路信息易于混淆的建筑物等地物，能夠很好地驗證本文所提方法的有效性。道路提取具體實現(xiàn)流程如圖1所示。

圖1 道路提取流程

1.1 FCM聚類算法

模糊C均值聚類(fuzzy C-means,FCM)，是一種基于目標(biāo)函數(shù)的模糊聚類算法。該算法采用隸屬度確定每個數(shù)據(jù)點屬于某個聚類的程度，是對傳統(tǒng)硬聚類算法的一種改進[8]。

假設(shè)數(shù)據(jù)集X={x1,x2,…,xn}，則它的模糊c劃分可用模糊矩陣U=uij表示，矩陣U的元素uij表示第j(j=1,2,…,n)個數(shù)據(jù)點屬于第i(i=1,2,…,c)類的隸屬度，其中uij滿足以下條件

(1)

(2)

1.2 分水嶺分割算法分割建成區(qū)影像

經(jīng)過FCM聚類處理后的影像被粗分為具有光譜相似性的包含道路、建筑物等地物的建成區(qū)和與道路光譜差異較大的非建成區(qū)，去除非建成區(qū)地物后，剩余基于像素FCM聚類的建成區(qū)，因此需要對建成區(qū)進行分割，以便對后續(xù)的基于對象的MKSVM道路提取進行處理。

分水嶺分割方法，是一種基于拓撲理論的數(shù)學(xué)形態(tài)學(xué)的分割方法，其基本思想是把圖像看作測地學(xué)上的拓撲地貌，圖像中每一點像素的灰度值表示該點的海拔高度，每一個局部極小值及其影響區(qū)域稱為集水盆，而集水盆的邊界則形成分水嶺[9]。由于受噪聲等因素影響分水嶺分割算法容易造成過分割現(xiàn)象。本文采用文獻[9]中的標(biāo)記分水嶺方法，能夠一定程度上抑制影像過分割，減少分割對象數(shù)目，減少后文工作的工作量，提高算法效率。

1.3 建成區(qū)分割對象的特征提取

分水嶺分割后的建成區(qū)影像對象既含有光譜特征，又含有大量的空間結(jié)構(gòu)特征。特征提取是通過MKSVM進行道路信息提取的關(guān)鍵步驟，特征提取的好壞直接影響后續(xù)道路提取的精度。

1.3.1 光譜特征

(1) 光譜均值：計算區(qū)域?qū)ο笾邢袼氐拿總€波段的光譜值，對每個波段光譜值求平均值作為區(qū)域?qū)ο蟮奶卣鳌?/p>

SAVE=SSUM/N

(3)

式中，SAVE為單波段對象光譜均值；SSUM為單波段區(qū)域?qū)ο髢?nèi)像素的光譜值之和；N為像素點個數(shù)。

(2) 均方差：高分辨率遙感影像譜內(nèi)差異較大，譜間差異較小，僅提取光譜均值作為光譜特征難以滿足分類需要，進一步提取均方差，與均值一起作為光譜特征。

1.3.2 紋理特征

區(qū)域形狀指數(shù)(regional shape index,RSI)是在PSI基礎(chǔ)上加以改進的區(qū)域空間特征，能夠很好地反映區(qū)域?qū)ο蟮募y理特征[10]。對于任一特定的區(qū)域?qū)ο?，從對象中心點沿不同方向發(fā)散一系列方向線，通過設(shè)定閾值T限制方向線條數(shù)。假設(shè)對象中心點為(x0,y0),方向角為θ，對于第tth條方向線上任意一坐標(biāo)點(mt,nt)，其方向線描述如下

nt=tan(t×θ)×(mt-x0)+y0

(4)

方向線長度計算公式如下

(5)

式中，(xt,yt)為方向線與邊界的交點。算出所有的方向線長度，并進行求和運算就可以得到RSI，計算公式如下

(6)

1.3.3 形狀特征

(1) 長寬比：城市道路通常表現(xiàn)為條帶狀結(jié)構(gòu)，長寬比較大。使用區(qū)域?qū)ο蟮淖钚⊥饨泳匦蔚拈L度與寬度比M來近似區(qū)域?qū)ο蟮拈L寬比[11]。

M=LR/W

(7)

式中，LR為外接矩形的長；W為外接矩形的寬。

(2) 矩形度：目標(biāo)區(qū)域?qū)ο竺娣e與包圍該對象的最小矩形面積的比值，適合檢測曲線道路。

R=S0/SR

(8)

式中，S0為目標(biāo)區(qū)域?qū)ο蟮拿娣e；SR為包圍該對象的最小矩形的面積；R的大小可以區(qū)分目標(biāo)對象與最小矩形的接近程度，其范圍是(0,1]。

(3) 線性特征指數(shù)(linear feature index,LFI)：道路是暗長區(qū)域，遍歷已得到道路的候選區(qū)域，利用一個最小外接矩形包圍連通區(qū)域，采用文獻[12]中的方法提取每個連通區(qū)域的中心線。連通區(qū)域轉(zhuǎn)換成矩形，該矩形滿足

LW=np

(9)

約束條件為：L=連通域中心線的長度;np=連通域的面積。其LFI可通過下式計算

(10)

式中，L和W分別是新矩形的長和寬；np是連通域的面積。道路有大的LFI值，LFI較小的連通區(qū)域被去除。

2 MKSVM道路提取及數(shù)學(xué)形態(tài)學(xué)處理

2.1 MKSVM基本原理

大量的試驗研究證明SVM對于小樣本、高維特征的分類具有較好的性能，SVM的成功運用促進了核方法的研究與發(fā)展[13-14]。MKSVM在SVM的基礎(chǔ)上以多核函數(shù)代替單核函數(shù)，將不同樣本的特征分量分別輸入不同的核函數(shù)進行映射，使數(shù)據(jù)在新特征空間中得到更好的表達，從而提高學(xué)習(xí)性能和分類精度[15]。其中核函數(shù)公式為

(11)

式中，dm為多核權(quán)系數(shù)；km為基核。

因此，與SVM問題類似，多核SVM類問題可描述為如下形式

(12)

約束條件為

(13)

將式(3)整理為以dm為變量的優(yōu)化問題

(14)

約束條件為

(15)

(16)

全局核函數(shù)泛化能力強，學(xué)習(xí)能力差；局部核函數(shù)學(xué)習(xí)能力強但泛化能力差。本文結(jié)合高分辨率遙感影像小樣本特征復(fù)雜的特點，以全局核函數(shù)POLY核函數(shù)與局部核函數(shù)RBF核函數(shù)線性加權(quán)組合，構(gòu)建MKSVM，采用文獻[16]中的蜂群算法求解MKSVM的權(quán)系數(shù)和相關(guān)核參數(shù)。人工選取訓(xùn)練樣本，將上文提取的特征進行歸一化處理后輸入MKSVM進行建成區(qū)影像的二次分類，提取出道路信息。

2.2 數(shù)學(xué)形態(tài)學(xué)處理

經(jīng)過上述MKSVM提取后的道路信息，仍存在由于車道線、車輛、陰影和其他噪聲干擾造成的孔洞現(xiàn)象。數(shù)學(xué)形態(tài)學(xué)利用具有一定形態(tài)的結(jié)構(gòu)元素去度量和提取圖像中的對應(yīng)形狀[17]，能夠很好地填充孔洞。在本文中用到的數(shù)學(xué)形態(tài)學(xué)算子是膨脹運算、腐蝕運算、擊中或不擊中變換。

設(shè)A和B為兩個非空集合，B為具有一定形態(tài)的結(jié)構(gòu)元素，那么：

腐蝕的定義為

AΘB=xBx?A

(17)

膨脹的定義為

(18)

擊中或不擊中變換

A?B=(AΘB1)∩(AcΘB2)

(19)

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

試驗選取重慶市某地區(qū)QuickBird影像為研究對象，實現(xiàn)城市道路信息的快速準(zhǔn)確提取，影像拍攝時間2016年8月，影像大小434×434。QuickBird影像含有R、G、B 3個波段，空間分辨率0.61 m。各試驗步驟結(jié)果如圖2所示。

圖2 道路提取步驟

從目視效果來看，本文所提方法可以較好地提取出道路信息。圖2(b)是經(jīng)FCM方法聚類后提取的建成區(qū)影像，建成區(qū)內(nèi)既含有道路，也含有與道路光譜相似的建筑物、水泥廣場等；圖2(c)經(jīng)過標(biāo)記分水嶺分割后，一定程度上抑制了過分割現(xiàn)象；圖2(d)利用MKSVM二次分類獲得了粗略的道路信息，從圖中可以看出，道路已經(jīng)基本上完整地被提取出來，但仍含有一些易混的噪聲；圖2(e)為經(jīng)過數(shù)學(xué)形態(tài)學(xué)處理后的道路提取結(jié)果，可以看出，道路提取效果比較好。

為了驗證本文所提方法的有效性，采用常用的道路提取算法邊緣檢測算法、多尺度提取及單一核函數(shù)SVM作為對比試驗，并采用客觀與主觀結(jié)合的方式對提取效果進行評價。具體的各方法道路提取效果圖如圖3所示。

圖3 不同方法道路提取結(jié)果

從目視效果來看，邊緣檢測的道路提取方法效果最差，存在一定程度上的道路誤分現(xiàn)象，一些建筑物被誤分成道路；多尺度提取及單核SVM的道路提取方法基本上可以準(zhǔn)確提取出道路，但仍有一定的凹洞、不連貫現(xiàn)象存在；本文所提方法經(jīng)過MKSVM提取和數(shù)學(xué)形態(tài)學(xué)處理后，能夠完整地提取出道路信息，提取效果相比其他3種方法有了明顯提高。

為了對道路的提取效果進行客觀評價，即將道路提取結(jié)果與真實的地物地面數(shù)據(jù)進行比較，本文通過建立誤差矩陣并進行相關(guān)的運算得到所分地物的分類精度；選用總體精度與Kappa系數(shù)對道路提取效果進行客觀評價[18]。

3.1 誤差矩陣

誤差矩陣是用來表示精度評價的一種標(biāo)準(zhǔn)格式，它是一個n×n的矩陣，其中n表示類別的數(shù)量。

(20)

式中，pij表示分類數(shù)據(jù)類型中第i類數(shù)據(jù)被劃分到j(luò)類中的總數(shù)。在誤差矩陣中，對角線上的值越大，說明分類精度越高，反之，若非對角線上的值越大，則分類精度越低。

3.2 總體精度(OA)

(21)

3.3 Kappa系數(shù)

Kappa系數(shù)是另外一種評價分類精度的方法，它采用離散的多元技術(shù)，考慮了誤差矩陣中的所有因素，是一個更為客觀的評價指標(biāo)，具體計算公式如下

(22)

為定量進行道路提取效果的評價，驗證本文所提方法的有效性，本文在原始圖像中隨機選取500個測試樣本點，并對道路和非道路樣本點進行個數(shù)統(tǒng)計。圖4為測試樣本點分布。道路提取結(jié)果評價見表1、表2。

圖4 測試樣本點分布

從表2中可以看出道路提取的總體精度達到了91.00%，Kappa系數(shù)達到了0.802 2，能夠有效地提取出道路信息。相比于邊緣檢測方法，多尺度方法及單核SVM方法都有一定程度的提高。

表1 測試樣本點道路提取結(jié)果

表2 不同方法道路提取綜合評價

為了更好地驗證本文所提方法的有效性，本文利用另外一幅WorldView-2影像進行道路提取試驗，各方法道路信息提取結(jié)果如圖5所示

圖5 不同方法提取效果

從表3中可以看出，本文所提方法與邊緣檢測方法和多尺度方法相比，總體精度和Kappa系數(shù)明顯提高，分別提高7.1%和4.4%，Kappa系數(shù)分別提高0.068 9和0.034 8。相比于單核SVM提取方法，提高幅度不大，總體分類精度和Kappa系數(shù)分別提高2.9%和0.024 9，但仍有一定效果的改善。由上述分析可知，本文方法道路提取效果更好。

表3 不同方法的道路提取評價結(jié)果

4 結(jié) 語

針對在高分辨率遙感影像中難以精確提取道路信息的特點及SVM分類性能的問題，本文提出了一種城市道路的多特征多核SVM分層提取的方法。綜合利用影像對象的多種特征，基于面向?qū)ο蠹夹g(shù)并結(jié)合多核學(xué)習(xí)和SVM分類器，實現(xiàn)了城市道路的有效提取?？傮w來看，本文所提方法能夠較精確地提取出城市道路信息，在地理國情普查與監(jiān)測、數(shù)字城市、智慧城市建設(shè)中具有一定的應(yīng)用價值。但本文由于利用多核SVM進行分類，勢必導(dǎo)致計算復(fù)雜度的增加，下一步將就如何降低計算復(fù)雜度進行研究。

[1] 劉純，洪亮，陳杰,等.融合像素—多尺度區(qū)域特征的高分辨率遙感影像分類算法[J].遙感學(xué)報，2015，19(2)：228-239.

[2] 孫可，張鈞萍.局部自適應(yīng)方向模板匹配的高分辨率遙感影像道路提取[C]∥2015年光學(xué)精密工程論壇論文集.長春：《光學(xué)精密工程》編輯部，2015：509-515.

[3] ZHAO L，WANG X.Road Extraction in High Resolution Remote Sensing Images Based on Mathematic Morphology and Snake Model[C]∥International Congress on Image and Signal Processing.[S.l.]：IEEE，2010:1436-1440.

[4] ZHAO W Z，LUO L Q，GUO Z，et al.Road Extraction in Remote Sensing Images Based on Spectral and Edge Analysis[J].Spectroscopy and Spectral Analysis，2015，35(10)：2814-2819.

[5] 林鵬,阮仁宗,王玉強.一種基于面向?qū)ο蟮某擎?zhèn)道路自動提取方法研究[J].地理與地理信息科學(xué)，2016，32(1)：49-54.

[6] ZHU E，SONG W，DAI J，et al.Road Extraction of High-resolution Remote Sensing Images Based on Improved SVM[J].Science of Surveying and Mapping，2016，41(12)：224-228.

[7] MIAO Z，SHI W，GAMBA P，et al.An Object-based Method for Road Network Extraction in VHR Satellite Images[J].IEEE Journal of Selected Topics in Applied Earth Observations & Remote Sensing，2016，8(10):4853-4862.

[8] JIA S，LI X，GUO F，et al.A Fast and Robust FCM Algorithm with Local Information for Image Segmentation[J].Icic Express Letters，2015，9(3):869-874.

[9] 賈永紅，馮在梅，沈庭樂.基于標(biāo)記分水嶺分割的高分辨率影像水體信息提取[J].測繪通報，2014(9)：31-33.

[10] CHU S，HONG L，LIU C，et al.A New Regional Shape Index for Classification of High Resolution Remote Sensing Images[C]∥Third International Workshop on Earth Observation and Remote Sensing Applications.Changsha：IEEE,2014.

[11] 陳杰，邵權(quán)斌，梅小明，等.結(jié)合尺度空間的面向?qū)ο蟾叻直媛视跋癯鞘械缆诽崛J].測繪工程，2016，25(12)：5-11.

[12] SHI W，MIAO Z，WANG Q，et al.Spectral-spatial Classification and Shape Features for Urban Road Centerline Extraction[J].IEEE Geoscience and Remote Sensing Letters，2014，11(4):788-792.

[13] ZHENG W，TANG H，HUANG T S.Emotion Recognition from Non-frontal Facial Images[M]∥Emotion Recogn-ition：A Pattern Analysis Approach.Hoboken：John Wiley and Sons，Inc.，2015:183-213.

[14] CUSANO C，NAPOLETANO P，SCHETTINI R.Remote Sensing Image Classification Exploiting Multiple Kernel Learning[J].IEEE Geoscience and Remote Sensing Letters，2015，12(11):2331-2335.

[15] 晁拴社，楚恒，王興.高光譜圖像數(shù)據(jù)的多尺度多核SVM分類[J].計算機與現(xiàn)代化，2016(2):11-14.

[16] 吳一全，殷駿，戴一冕，等.基于蜂群優(yōu)化多核支持向量機的淡水魚種類識別[J].農(nóng)業(yè)工程學(xué)報，2014，30(16):312-319.

[17] GONZALEZ R C，WOODS R E，EDDINS S L.數(shù)字圖像處理[M].阮秋琦，譯.北京：電子工業(yè)出版社，2014.

[18] 趙時英.遙感應(yīng)用分析原理與方法[M].北京：科學(xué)出版社，2003.