王建華　秦其明　高中靈　葉昕　孟晉杰

摘要：由于道路與建筑物等其他不透水層存在光譜相似性，導(dǎo)致僅利用光譜信息進(jìn)行道路提取的效果不佳。本文針對(duì)高等級(jí)城市道路目標(biāo)，提出了一種加入空間紋理信息的遙感圖像道路提取方法。首先，對(duì)圖像進(jìn)行空間自相關(guān)Moran指數(shù)計(jì)算，提取圖像空間紋理信息，并將其加入到原始光譜波段中；其次，通過(guò)建立知識(shí)模提取假設(shè)道路段，并對(duì)提取結(jié)果進(jìn)行假設(shè)驗(yàn)證；最后，采用數(shù)學(xué)形態(tài)學(xué)的方法對(duì)驗(yàn)證后的結(jié)果進(jìn)行后處理。以空間分辨率為0.1m的航空影像為數(shù)據(jù)源，對(duì)本方法進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)結(jié)果表明，加入空間紋理信息的遙感圖像道路提取精度總體達(dá)到88%，比不加入空間紋理的提取精度要提高約5%。

關(guān)鍵詞：空間紋理信息；Moran指數(shù)；圖像處理；假設(shè)驗(yàn)證；道路提取

中圖分類(lèi)號(hào)：TP751 文獻(xiàn)標(biāo)識(shí)碼：A

道路作為現(xiàn)代化建設(shè)的基礎(chǔ)設(shè)施，在人們?nèi)粘Ｉ詈蛙娛律隙季哂胁豢珊鲆暤淖饔茫詮倪b感圖像上提取道路具有重要意義，它可以為城市規(guī)劃提供參考，為地圖更新提供數(shù)據(jù)，為國(guó)土資源管理提供依據(jù)，為搶險(xiǎn)救災(zāi)提供指導(dǎo)。針對(duì)遙感圖像道路提取問(wèn)題，國(guó)內(nèi)外很多專(zhuān)家和學(xué)者都開(kāi)展了相應(yīng)的工作，Trinder等人在高分辨率圖像上提取道路段，并利用道路語(yǔ)義模型對(duì)漏提的道路段進(jìn)行推理、連接，從而得到道路網(wǎng)提取結(jié)果。Singh等人利用自適應(yīng)全局閾值和數(shù)學(xué)形態(tài)學(xué)操作從高分辨率遙感圖像中提取道路網(wǎng)。Shi等人運(yùn)用自適應(yīng)鄰域法和空間分析法在遙感圖像中準(zhǔn)確地提取道路中心線。Yuan等人利用局部興奮全局抑制振蕩網(wǎng)絡(luò)對(duì)遙感圖像進(jìn)行分割，從而進(jìn)行道路提取。Sent-hilnath等人通過(guò)歸一化割算法從全色遙感圖像中很好地提取出道路。Miao等人提出了一種半自動(dòng)化的道路中心線提取方法，該方法通過(guò)人工設(shè)置道路種子點(diǎn)，結(jié)合核密度估計(jì)和短程線法提取道路中心線。羅慶洲等人結(jié)合了道路的光譜特征和幾何特征從影像中很好地提取了道路信息。唐偉等人基于Snake模型的優(yōu)勢(shì)及缺陷，提出一種高分辨率多光譜圖像道路提取方法。以上方法中都用到了道路的光譜信息，但由于道路與建筑物等地物光譜相似性，從而影響道路提取精度。隨著遙感技術(shù)的發(fā)展，高分辨率遙感圖像包含了豐富的空間細(xì)節(jié)信息，使得地物的內(nèi)部結(jié)構(gòu)清晰可辨。這為加入空間紋理信息進(jìn)行道路提取提供了數(shù)據(jù)基礎(chǔ)。

本文針對(duì)上述問(wèn)題，提出了一種加入空間紋理信息的遙感圖像道路提取方法。將空間紋理信息加入到道路提取中，通過(guò)建立道路提取的假設(shè)和驗(yàn)證模型，基于知識(shí)的方法從遙感圖像中提取道路。本方法適用于遙感圖像中高等級(jí)城市道路提取，具有較高的精度。

1 研究方法

本文在研究遙感圖像中道路特征知識(shí)的基礎(chǔ)上，將圖像的空間紋理信息加入到遙感圖像道路提取中，通過(guò)知識(shí)模型提取道路。主要包括紋理信息提取、道路提?。僭O(shè)道路提取和假設(shè)驗(yàn)證）和提取結(jié)果后處理3個(gè)主要步驟。具體流程如圖1所示。

1.1 空間紋理信息提取

紋理是指遙感圖像中地物內(nèi)部色調(diào)有規(guī)則變化形成的一種影像結(jié)構(gòu)。高分辨率遙感圖像中包含著豐富的空間細(xì)節(jié)信息，能夠使地物內(nèi)部結(jié)構(gòu)更加清晰。為了解決道路和建筑物等其他不透水層混淆問(wèn)題，提高道路提取的準(zhǔn)確性，本文引入了空間紋理信息。采用空間統(tǒng)計(jì)（Moran指數(shù)）的方法提取圖像的紋理信息。首先利用主成分變化對(duì)多波段圖像進(jìn)行分析，然后利用Moran指數(shù)對(duì)第一主分量進(jìn)行統(tǒng)計(jì)，提取圖像的空間紋理。其中，Moran指數(shù)的定義如下：式中，x_i和x_j分別指空間單元i和j的屬性值，n是指空間單元的個(gè)數(shù)，ω_ij為空間權(quán)矩陣，表示空間單元i和j之間的影響程度。

1.2 基于知識(shí)的道路提取

道路在遙感圖像中一般呈現(xiàn)：1）寬度變化??；2）與周?chē)尘按嬖谝欢ǖ姆床?，表現(xiàn)出明顯的邊緣特征；3）道路內(nèi)部總體灰度比較均勻；4）道路對(duì)象具有較高的矩形度和長(zhǎng)寬比。本文在研究遙感圖像中道路的特征表達(dá)與抽取基礎(chǔ)上，建立道路提取的知識(shí)模型（假設(shè)模型和驗(yàn)證模型）。

道路特征知識(shí)定量表達(dá)的參數(shù)有亮度、標(biāo)準(zhǔn)差、矩形度、長(zhǎng)寬比和面積等。通過(guò)研究道路在遙感圖像上的表現(xiàn)特征，本文選取亮度、標(biāo)準(zhǔn)差建立道路提取假設(shè)模型，選取矩形度、長(zhǎng)寬比和面積建立驗(yàn)證模型。其中道路假設(shè)模型定義如下：式中，H_road是指假設(shè)道路，B_road和S_road分別為某對(duì)象的亮度和標(biāo)準(zhǔn)差，b₁和b₂是亮度閾值，s₁和s₂是標(biāo)準(zhǔn)差閾值。針對(duì)某對(duì)象，當(dāng)其亮度和標(biāo)準(zhǔn)差滿(mǎn)足B_road∈[b₁，b₂]或者S_road∈s₁，s₂]時(shí)，將其定義為假設(shè)道路段。由于異物同譜現(xiàn)象的存在，導(dǎo)致假設(shè)道路段中存在部分如車(chē)輛、房屋等虛假道路目標(biāo)。本研究通過(guò)建立假設(shè)驗(yàn)證模型，去除這些虛假道路目標(biāo)。道路對(duì)象的長(zhǎng)寬比和矩形度一般較大，所以本研究用道路的長(zhǎng)寬比和矩形度作為驗(yàn)證模型參數(shù)，并將面積屬性加入驗(yàn)證模型中去除零碎的細(xì)小對(duì)象。驗(yàn)證模型表示如下：式中，V_road是指驗(yàn)證后的道路；R，W和A分別指某對(duì)象的矩形度、長(zhǎng)寬比和面積；r₁，ω₁和a₁分別指矩形度閾值、長(zhǎng)寬比閾值和面積閾值。公式（3）的含義為，在假設(shè)道路段中，當(dāng)?shù)缆穼?duì)象的矩形度和長(zhǎng)寬比滿(mǎn)足R∈（0，r₁）或者W∈（0，ω₁）時(shí)，將其從假設(shè)道路段中去除。然后對(duì)假設(shè)道路段進(jìn)行區(qū)域合并，將面積小于a₁的對(duì)象從假設(shè)道路對(duì)象中去除，從而得到驗(yàn)證后的道路段。

1.3 提取結(jié)果后處理

由于車(chē)輛、行車(chē)線、建筑物等地物的影響很大，盡管知識(shí)模型在一定程度上可以減小影響，但是道路提取結(jié)果中仍存在很多噪聲和孔洞，本文采用數(shù)學(xué)形態(tài)學(xué)閉運(yùn)算對(duì)道路提取的結(jié)果進(jìn)行后處理來(lái)提高道路提取的準(zhǔn)確性，為了保證提取道路的完整性以及道路邊緣的平滑性，此處選取八鄰域作為數(shù)學(xué)形態(tài)學(xué)運(yùn)算的結(jié)構(gòu)元素進(jìn)行后處理。

2 實(shí)驗(yàn)與結(jié)果

為了驗(yàn)證本方法的有效性和優(yōu)越性，本文通過(guò)實(shí)測(cè)數(shù)據(jù)進(jìn)行實(shí)驗(yàn)，并對(duì)結(jié)果進(jìn)行精度評(píng)價(jià)。

2.1 數(shù)據(jù)源

本研究選取廣東省陽(yáng)江地區(qū)航空影像數(shù)據(jù)為實(shí)驗(yàn)數(shù)據(jù)，成像時(shí)間是2013年4月，圖像空間分辨率為0.1m，圖像大小為2808×2719像素。道路寬度為10m左右，道路類(lèi)型是高等級(jí)城市道路，實(shí)驗(yàn)數(shù)據(jù)如圖2（a）所示。圖像中包括道路、房屋、植被、車(chē)輛、行車(chē)線等目標(biāo)。圖像具有紅、綠、藍(lán)三個(gè)波段信息。

2.2 道路目標(biāo)提取

對(duì)原始圖像進(jìn)行主成分分析，第一主分量如圖2（b）所示，它包含了原始圖像的主要信息。利用Moran指數(shù)對(duì)圖2（b）進(jìn)行空間統(tǒng)計(jì)，提取圖像的空間紋理信息如圖2（c）所示，道路在圖中呈現(xiàn)低值聚集，比周?chē)匚锏牧炼让黠@要低，道路的邊緣突出，有助于更加完整地將道路從周?chē)匚镏蟹指铋_(kāi)。因此，將圖2（c）加入到原始圖像中，并采用多尺度分割的方法對(duì)加入空間紋理信息后的圖像進(jìn)行分割，設(shè)定分割尺度為200，分割結(jié)果如圖2（d）所示，從圖中可以看出道路的完整性保留得很好。利用道路提取假設(shè)模型進(jìn)行假設(shè)道路的提取，提取結(jié)果如圖2（e）所示（其中，b₁=74；b₂=102；s₁=12；s₂=15）。利用驗(yàn)證模型進(jìn)行假設(shè)驗(yàn)證，去除虛假道路目標(biāo)（r₁=0.6；ω₁=2；a₁=5000pixels），圖2（f）為驗(yàn)證結(jié)果。利用數(shù)學(xué)形態(tài)學(xué)閉運(yùn)算對(duì)提取結(jié)果進(jìn)行后處理，將結(jié)果與原始圖像進(jìn)行疊加，并與真實(shí)結(jié)果進(jìn)行對(duì)比。圖2（g）為疊加效果圖，其中紅色部分為正確提取的道路，藍(lán)色部分為漏提取的道路，綠色部分為錯(cuò)提取為道路。從圖中可以看出有一小段道路漏提，分析發(fā)現(xiàn)該小段道路左側(cè)邊緣受到了破壞，同時(shí)道路上塵土較多，導(dǎo)致與其他道路段的紋理和光譜特征產(chǎn)生了差異，故而被錯(cuò)分到其他地物中。為了證明本方法的優(yōu)越性，本文設(shè)計(jì)了對(duì)比實(shí)驗(yàn)，圖2（h）為不加入空間紋理信息提取的結(jié)果。通過(guò)對(duì)比分析，可以直觀地看出加入紋理信息提取的結(jié)果更精確。

2.3 精度評(píng)價(jià)

為了對(duì)道路提取結(jié)果進(jìn)行精度評(píng)價(jià)，本文選擇3個(gè)最廣泛被人認(rèn)可的評(píng)價(jià)指標(biāo)，具體如下：式中，E₁為生產(chǎn)者精度，E₂為用戶(hù)精度，E₃為提取質(zhì)量。k_TP為正確提取為道路的區(qū)域，k_FN為未被提取的道路區(qū)域，k_FP為錯(cuò)誤提取為道路的區(qū)域。利用選取的3個(gè)評(píng)價(jià)指標(biāo)，分別對(duì)上述實(shí)驗(yàn)中是否加入空間紋理信息提取的道路區(qū)域（面積）進(jìn)行評(píng)價(jià)，以此來(lái)突出本方法的優(yōu)越性。具體的精度評(píng)價(jià)結(jié)果如表1所示。

由表1可知，本方法能夠很好地提取道路，E₁，E₂和E₃分別達(dá)到了94.29%，88.62和84.11%，平均精度達(dá)到了88%，比不加入空間紋理信息提取精度提高了大約5%。然而，由于混合像元的影響，導(dǎo)致道路邊界不清楚，從而錯(cuò)誤地將道路邊界外側(cè)鄰近的部分區(qū)域提取為道路，這使得道路提取的用戶(hù)精度相對(duì)較低。

3 結(jié)論

本文針對(duì)高等級(jí)城市道路目標(biāo)，提出了一種加入空間紋理信息的遙感圖像道路提取的方法，有效地解決了道路與其他不透水層的光譜相似性問(wèn)題。該方法在深入研究道路特征知識(shí)表達(dá)的基礎(chǔ)上，建立了道路提取的假設(shè)模型和驗(yàn)證模型，減少了道路中車(chē)輛、行車(chē)線等地物的不利影響。最后利用廣東省陽(yáng)江地區(qū)影像數(shù)據(jù)進(jìn)行了實(shí)驗(yàn)，實(shí)驗(yàn)結(jié)果表明，本方法能夠克服道路周?chē)h(huán)境的影響，有效地從遙感圖像中提取城市主干道路。