王曉寧

摘? 要： 針對(duì)傳統(tǒng)城市空間尺度要素特征提取系統(tǒng)沒有分類處理遙感影像多尺度紋理特征，導(dǎo)致提取效果差、精度低等問題，設(shè)計(jì)一種基于機(jī)器視覺的城市空間尺度要素特征提取系統(tǒng)?；跈C(jī)器視覺技術(shù)將數(shù)字圖像處理與計(jì)算機(jī)理論相結(jié)合，采用灰度共生矩陣法提取遙感影像多尺度紋理特征，并利用決策樹分類法對(duì)其展開分類后，獲取城市遙感影像分類結(jié)果。通過(guò)直方圖閾值分割方法提取水體特征，利用決策樹分類法提取建筑物特征，通過(guò)紋理參數(shù)分割、數(shù)字形態(tài)學(xué)處理、特征矢量化等步驟提取植被特征，運(yùn)用數(shù)學(xué)形態(tài)學(xué)和邊緣檢測(cè)提取道路、陰影以及裸地特征。實(shí)驗(yàn)結(jié)果表明，該系統(tǒng)空間尺度要素特征提取，準(zhǔn)確率高達(dá)95.02%，精度高;并且速度快、效率高，說(shuō)明該系統(tǒng)實(shí)用性較好。

關(guān)鍵詞： 城市空間尺度; 要素特征提取; 機(jī)器視覺; 紋理特征提取; 決策樹分類; 系統(tǒng)設(shè)計(jì)

中圖分類號(hào)： TN911?34; TP75? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼： A? ? ? ? ? ? ? ? ? ? ? 文章編號(hào)： 1004?373X（2020）22?0168?05

Abstract： In allusion to the problem that traditional urban spatial scale element feature extraction system does not classify and process the multi?scale texture features of remote sensing images， which results in the poor extraction effect and low precision， an urban spatial scale elements feature extraction system based on machine vision is designed. The digital image processing is combined with computer theory based on machine vision technology， the multi?scale texture features of remote sensing images are extracted by means of the gray scale co?existing matrix method， and the remote sensing images are classified by means of the decision tree classification to obtain the classification results of urban remote sensing images. The water features are extracted by means of the histogram threshold segmentation， the building features are extracted by means of the decision tree classification， the vegetation features are extracted by the texture parameter segmentation， digital morphological processing and feature vectorization， and the features of the road， shadow and bare ground are extracted by the mathematical morphology and edge detection. The experimental results show that the extraction accuracy of spatial scale element feature of the system is 95.02%， the extraction precision is high， and the extraction efficiency and speed is excellent， which indicates that the system is practical.

Keywords： urban spatial scale; element feature extraction; machine vision; textural features extraction; decision tree classification; system design

0? 引? 言

伴隨著新政策的施行，城市空間布局發(fā)展體系逐漸完善，城市魅力呈現(xiàn)集中展示趨勢(shì)，城市化范圍擴(kuò)大是發(fā)展的必然趨勢(shì)，為降低城市景觀格局遭到破壞的概率，需要精確提取城市空間尺度要素特征[1]。常利用遙感技術(shù)對(duì)城市的各項(xiàng)數(shù)據(jù)進(jìn)行采集，將不斷擴(kuò)大的城市情況實(shí)時(shí)捕捉。城市空間一般具備集中或零散的公用設(shè)施與市政公共設(shè)施，城市空間尺度要素主要包括建筑物、植被、水體、道路、裸地、陰影等，高分辨率遙感影像能夠精準(zhǔn)捕捉這些尺度要素的細(xì)節(jié)[2]。

遙感影像的采集和處理需要利用到機(jī)器視覺技術(shù)，機(jī)器視覺技術(shù)將數(shù)字圖像處理與計(jì)算機(jī)理論相結(jié)合，對(duì)采集到的影像信息進(jìn)行處理，提取出有用的影像信息展開識(shí)別，即在對(duì)影像進(jìn)行處理時(shí)，使用相關(guān)數(shù)字技術(shù)或數(shù)字計(jì)算機(jī)，對(duì)影像進(jìn)行處理與運(yùn)算，獲取所需影像信息，例如使用相關(guān)算法對(duì)影像進(jìn)行紋理特征處理，獲得城市空間尺度要素特征[3]。紋理特征能夠精準(zhǔn)展現(xiàn)高分辨率遙感影像細(xì)節(jié)信息，為此結(jié)合遙感影像中的紋理信息可精準(zhǔn)提取城市空間尺度要素特征[4]。城市空間尺度要素分為多種尺度，具有差異性的空間分辨率，可以將不同的尺度要素進(jìn)行分類。有研究表明，多尺度紋理要素特征提取可以對(duì)尺度要素特征進(jìn)行更完善的描述，使分類精度得到提高，如果在尺度要素特征提取時(shí)直接使用原始圖像，噪聲干擾大、數(shù)據(jù)量龐雜的原始圖像會(huì)對(duì)提取結(jié)果造成影響[5]。

在城市空間環(huán)境中，尺度要素分布情況復(fù)雜，每種尺度要素都有適合自己的尺度，在進(jìn)行尺度要素信息提取時(shí)如果使用統(tǒng)一的尺度層次就無(wú)法得到精準(zhǔn)的影像信息[6]。以建筑物為例，不同的屋頂顏色及材質(zhì)會(huì)造成不同的光譜差異。在影像上會(huì)呈現(xiàn)出不同的建筑物尺度要素。

基于上述分析，本文利用機(jī)器視覺設(shè)計(jì)城市空間尺度要素特征提取系統(tǒng)，實(shí)現(xiàn)空間尺度重要要素的特征提取。

1? 城市空間尺度要素特征提取系統(tǒng)設(shè)計(jì)

使用高分辨率遙感影像，把城市空間尺度要素分別劃分為植被、水體、陰影、建筑物、道路、裸地六部分，系統(tǒng)對(duì)這六部分進(jìn)行尺度要素特征提取，提取流程見圖1。

使用灰度發(fā)生矩陣提取紋理特征，采用決策樹分類法展開監(jiān)督分類，根據(jù)已有研究，對(duì)高分辨率遙感影像中的大區(qū)域進(jìn)行劃分，將城市的空間尺度要素特征提取出來(lái)。

1.1? 多尺度紋理特征提取

在城市空間中，灰度的分布不斷交替變化形成紋理，在遙感圖像中空間灰度呈現(xiàn)相關(guān)特性，這是由于在空間中，一定距離間隔的兩個(gè)像素存在灰度關(guān)系。在對(duì)紋理特征進(jìn)行描述時(shí)常常使用灰度共生矩陣，遙感影像內(nèi)任意兩個(gè)點(diǎn)之間的灰度空間相關(guān)特性通過(guò)矩陣呈現(xiàn)，該矩陣是由影像內(nèi)的灰度值之間的二階聯(lián)合條件概率密度組成的[7?8]。使用灰度共生矩陣分析紋理時(shí)需要對(duì)移動(dòng)方向、移動(dòng)窗口大小以及移動(dòng)步長(zhǎng)這三個(gè)參數(shù)加以考慮。多尺度紋理特征提取時(shí)移動(dòng)方向選取30°，移動(dòng)步長(zhǎng)取值為2，移動(dòng)窗口大小選取4×4，6×6，…，26×26。在對(duì)紋理特征進(jìn)行統(tǒng)計(jì)時(shí)，以均勻性、非相似度、均值、熵、方差、相關(guān)度、均勻性、角二階矩等統(tǒng)計(jì)量來(lái)描述。

不同地物類型的紋理特征伴隨著紋理窗口大小的變化呈現(xiàn)的表達(dá)程度不同，地物類型紋理尺度對(duì)分類的影響很大。為了紋理尺度在優(yōu)化后的分類效果能夠達(dá)到良好的程度，對(duì)各類別地物之間的可分離性進(jìn)行表達(dá)時(shí)運(yùn)用平均J?M（Jeffries?Matusita）距離。J?M距離表達(dá)式如下：

式中：[Wx]表示類別[x]的矩陣樣本協(xié)方差;[Wy]表示類別[y]的矩陣樣本協(xié)方差;[Mx]和[My]分別表示[x]和[y]的樣本均值向量;[JM]表示平均J?M距離;[qx]和[qy]分別表示[x]和[y]的先驗(yàn)概率;[JMxy]表示這兩個(gè)類別之間的J?M距離;[Bxy]表示紋理特征因子;[N]表示閾值。

按照以上公式，選取樣本進(jìn)行訓(xùn)練，隨著紋理尺度的變化，平均J?M距離也產(chǎn)生變化。以下紋理窗口中尺度要素類型平均J?M距離最大：建筑物和裸地的紋理窗口是6×6，道路和植被的紋理窗口是8×8，水體的紋理窗口是16×16，陰影的紋理窗口是4×4，若要各種尺度要素的平均可分離性呈現(xiàn)較小的狀態(tài)，則紋理窗口值需達(dá)到21×21。

1.2? 決策樹分類法

決策樹分類法在遙感影像分類領(lǐng)域優(yōu)勢(shì)明顯，運(yùn)算效率高，清晰且靈活，應(yīng)用范圍廣泛，每個(gè)決策樹的組成都需要一整個(gè)內(nèi)部節(jié)點(diǎn)、單一根節(jié)點(diǎn)、大量分支節(jié)點(diǎn)和無(wú)數(shù)葉節(jié)點(diǎn)。在整個(gè)決策過(guò)程中，任何需要測(cè)試的屬性都來(lái)自決策樹的每個(gè)內(nèi)部節(jié)點(diǎn)，測(cè)試的結(jié)果可由各分支表示，圖像分類結(jié)果以各葉節(jié)點(diǎn)表示。

以CART算法作為決策樹構(gòu)建算法，該算法結(jié)構(gòu)清晰、運(yùn)算簡(jiǎn)單、處理數(shù)據(jù)時(shí)高效迅速。節(jié)點(diǎn)的分割規(guī)則是以各屬性最佳分割點(diǎn)作為各節(jié)點(diǎn)的選擇，在此節(jié)點(diǎn)分割出的兩個(gè)節(jié)點(diǎn)基礎(chǔ)上再次分割，按照這個(gè)規(guī)律不停分割，直到無(wú)法分割為止，通過(guò)該算法得到的決策樹稱為二叉樹，原因是各節(jié)點(diǎn)都存在兩個(gè)分支。

使用多尺度紋理特征提取和決策樹分類法得到的城市遙感影像分類結(jié)果，對(duì)城市空間尺度要素中的水體、建筑、植被等特征進(jìn)行提取。

1.3? 尺度要素特征提取

1.3.1? 水體特征提取

若要得到城市空間水體特征，需要在城市遙感影像分類結(jié)果基礎(chǔ)上，通過(guò)大區(qū)域分割方法獲得。使用直方圖閾值分割方法對(duì)水體進(jìn)行特征提取以解決水體灰度值過(guò)低不利于提取的問題。水體分割閾值的確定需要使用直方圖判定，該直方圖通過(guò)人機(jī)交互水體分割閾值獲得，閾值是與灰度直方圖內(nèi)第一個(gè)峰相對(duì)應(yīng)的波谷[9]。通過(guò)求出曲線極小值確定波谷，極小值的確定條件為：同時(shí)符合一階差分和二階差分分別等于0和大于0。綜上所述，特征提取方法步驟如下：

1） 以直方圖統(tǒng)計(jì)原始影像數(shù)據(jù)，將得到的數(shù)據(jù)以N級(jí)灰度劃分。

2） 對(duì)步驟1）中得到的原始直方圖展開中值濾波處理，去除直方圖中的毛刺[10]。

3） 為排除局部抖動(dòng)導(dǎo)致的誤差，對(duì)已經(jīng)去除毛刺的直方圖繼續(xù)進(jìn)行平滑處理。

4） 以一階差分[d′c=dc+1-dz]處理經(jīng)平滑處理后的直方圖，對(duì)一階直方圖從左到右檢測(cè)，其中，第一個(gè)從負(fù)轉(zhuǎn)變到正的躍階點(diǎn)即為波谷點(diǎn)，閾值[N]是在一階直方圖內(nèi)檢測(cè)到的與極小點(diǎn)相對(duì)應(yīng)的灰度值。

（5） 對(duì)原始圖像進(jìn)行分割處理，將小于閾值[N]的像元設(shè)置為水體，大于閾值[N]的像元判定為水體以外的尺度要素。

經(jīng)過(guò)以上5個(gè)步驟可以提取出城市空間尺度要素中的水體特征，經(jīng)矢量化處理得到多邊形水體區(qū)域。

1.3.2? 建筑物特征提取

按照建筑物屋頂顏色不同可以劃分為藍(lán)色、灰色、紅色3種。將每?jī)蓚€(gè)相鄰影像中的相對(duì)邊界和距離添加到紋理特征中，與光譜特征、形狀特征組成建筑物影像分類的原始特征子集。決策樹分類法的學(xué)習(xí)器有兩個(gè)重要的調(diào)整參數(shù)[11?13]，分別為修剪葉節(jié)點(diǎn)的置信系數(shù)[K]和節(jié)點(diǎn)需求實(shí)例中的最小數(shù)量參數(shù)[O]，這兩個(gè)調(diào)整參數(shù)的取值對(duì)分類結(jié)果起到重要作用，所以取值時(shí)要經(jīng)過(guò)參數(shù)尋優(yōu)處理。最小數(shù)量參數(shù)[O]和置信系數(shù)[K]取值分別為1～9和0.1～0.6，步長(zhǎng)分別為1.5和0.5。優(yōu)化決策樹得到最優(yōu)特征集，見表1。優(yōu)化后的最小數(shù)量參數(shù)[O]和置信系數(shù)[K]分別為4和0.2。

1.3.3? 植被特征提取

植被區(qū)域?qū)儆趬K狀基元，提取植被區(qū)域特征時(shí)一般需要使用影像中的紋理信息，特征提取方法如下：

1） 紋理參數(shù)分割。灰度共生矩陣通過(guò)高分辨率遙感影像生成，遙感影像和矩陣計(jì)算得出的5個(gè)紋理參數(shù)一起進(jìn)行分割計(jì)算[14?15]。

2） 數(shù)字形態(tài)學(xué)處理。為了消除影像中的斑塊和斑塊周圍的毛刺使得影像中的背景和斑塊形成鮮明對(duì)比，對(duì)分割后的二值圖像進(jìn)行數(shù)學(xué)形態(tài)算子處理。

3） 特征矢量化。以矢量化處理目標(biāo)圖像，單元形狀設(shè)置為多邊形，計(jì)算形狀指數(shù)，提取出植被區(qū)域。

1.3.4? 道路、陰影、裸地特征提取

運(yùn)用數(shù)學(xué)形態(tài)學(xué)和邊緣檢測(cè)對(duì)道路、陰影、裸地基元進(jìn)行提取，采用多尺度紋理特征提取高分辨率遙感影像的道路、陰影、裸地特征，利用決策樹分類法得出道路、陰影、裸地斑塊，使用啟發(fā)式連接規(guī)律連接道路的軸線，提取出道路、陰影、裸地區(qū)域特征。

2? 實(shí)驗(yàn)分析

為了驗(yàn)證基于機(jī)器視覺的城市空間尺度要素特征提取系統(tǒng)的整體有效性，需要對(duì)該系統(tǒng)進(jìn)行對(duì)比實(shí)驗(yàn)。實(shí)驗(yàn)前提取遼寧省沈陽(yáng)市城區(qū)內(nèi)某區(qū)域的城市空間尺度要素特征數(shù)據(jù)，用于對(duì)比驗(yàn)證本文系統(tǒng)的有效性，得到的實(shí)驗(yàn)結(jié)果如下。

2.1? 尺度要素特征提取效果對(duì)比

首先檢測(cè)不同系統(tǒng)方法的空間尺度要素特征提取效果，圖像越清晰，提取效果越好。實(shí)驗(yàn)結(jié)果如圖2所示。

從圖2a）中可以看出，在采用本文系統(tǒng)提取到的空間尺度要素中，可以清晰地提取出道路區(qū)域，同時(shí)其他小面積區(qū)域也能清楚地得到體現(xiàn)，這是因?yàn)樵撓到y(tǒng)采用多尺度紋理特征提取方法和決策樹分類方法，能夠更加清晰地獲取道路等不同事物的紋理特征，并保證道路特征不會(huì)與光譜相似的建筑物混淆，分類精度得到提高，總體效果較好。通過(guò)圖2b）可知，城市地理信息系統(tǒng)提取的空間尺度要素特征破碎，圖像十分模糊，道路特征展現(xiàn)出斷斷續(xù)續(xù)的形態(tài)，且與建筑物特征出現(xiàn)大面積混淆，陰影和水體不能明顯區(qū)分，大量建筑物被判定為道路和植被。

通過(guò)以上數(shù)據(jù)可以證明，本文系統(tǒng)提取的空間尺度要素特征精度較高，提取性能更好。

2.2? 尺度要素特征提取精度對(duì)比

為進(jìn)一步驗(yàn)證不同系統(tǒng)方法的特征提取有效性，對(duì)兩種系統(tǒng)處于不同信噪比下的特征提取精度進(jìn)行了檢測(cè)，結(jié)果如表2所示。

2.3? 尺度要素特征提取時(shí)間對(duì)比

由圖3可知，城市地理信息系統(tǒng)的空間尺度要素特征提取時(shí)間在3～6.8 s之間，波動(dòng)范圍較為明顯，而本文系統(tǒng)的空間尺度要素特征提取時(shí)間始終低于2 s，時(shí)間較為平穩(wěn)。特征提取時(shí)間越短，表明方法的提取速度越快，效率越高。通過(guò)以上數(shù)據(jù)可以證明，本文系統(tǒng)的空間尺度要素特征提取效率更高，實(shí)際使用價(jià)值更高。

3? 結(jié)? 論

將本文系統(tǒng)空間尺度要素特征提取與城市地理信息系統(tǒng)空間尺度要素特征提取進(jìn)行對(duì)比，城市地理信息系統(tǒng)所提取的尺度要素十分模糊，且斑點(diǎn)更多，混淆信息更多。在對(duì)兩個(gè)系統(tǒng)添加不同信噪比后對(duì)比特征提取精度，本文系統(tǒng)提取的各尺度要素平均精度均在95%以上，而城市地理信息系統(tǒng)提取的各尺度要素平均精度均在90%以下;本文系統(tǒng)的空間尺度要素特征提取時(shí)間始終低于2 s，而城市地理信息系統(tǒng)的空間尺度要素特征提取時(shí)間在3～6.8 s之間，證明本文系統(tǒng)實(shí)際應(yīng)用性能更優(yōu)秀。在今后的研究中，可從解決尺度要素混淆的角度出發(fā)，使用輔助或者優(yōu)化方法，提升提取精度。

參考文獻(xiàn)

[1] 吳志強(qiáng)，葉鍾楠.基于百度地圖熱力圖的城市空間結(jié)構(gòu)研究：以上海中心城區(qū)為例[J].城市規(guī)劃，2016，40（4）：33?40.

[2] 王釗，楊山，王玉娟，等.基于最小阻力模型的城市空間擴(kuò)展冷熱點(diǎn)格局分析：以蘇錫常地區(qū)為例[J].經(jīng)濟(jì)地理，2016，36（3）：57?64.

[3] 張少康，劉沛，魏冀明.基于風(fēng)環(huán)境分析的珠三角地區(qū)城鎮(zhèn)空間規(guī)劃引導(dǎo)[J].規(guī)劃師，2016，32（9）：118?122.

[4] 殷明，龐紀(jì)勇，魏遠(yuǎn)遠(yuǎn)，等.結(jié)合NSDTCT和稀疏表示的遙感圖像融合[J].光子學(xué)報(bào)，2016，45（1）：10?17.

[5] 周福送，霍宏，萬(wàn)衛(wèi)兵，等.基于強(qiáng)監(jiān)督部件模型的遙感圖像目標(biāo)檢測(cè)[J].計(jì)算機(jī)應(yīng)用，2016，36（6）：1714?1718.

[6] 李建，周屈，陳曉玲，等.近岸/內(nèi)陸典型水環(huán)境要素定量遙感空間尺度問題研究[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2018，43（6）：937?942.

[7] 徐磊，董捷，李璐，等.基于功能分區(qū)視角的長(zhǎng)江中游城市群國(guó)土空間特征及優(yōu)化[J].經(jīng)濟(jì)地理，2017，37（6）：76?83.

[8] 侯俊劍，吳艷靈，何文斌，等.基于三維聲場(chǎng)空間特征的故障診斷方法研究[J].振動(dòng)與沖擊，2018，37（13）：1?6.

[9] 李建，周屈，陳曉玲，等.近岸/內(nèi)陸典型水環(huán)境要素定量遙感空間尺度問題研究[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2018，43（6）：937?942.

[10] 段建民，石慧，戰(zhàn)宇辰.基于機(jī)器視覺篩選GPS衛(wèi)星信號(hào)的無(wú)人駕駛汽車組合導(dǎo)航方法[J].電子技術(shù)應(yīng)用，2016（1）：111?114.

[11] 呂志強(qiáng)，鄧睿，卿珊珊.大型山地城市建設(shè)用地空間擴(kuò)展及地形分異[J].水土保持研究，2017，24（1）：232?238.

[12] 唐燕妮.基于機(jī)器視覺的陶瓷大花瓶?jī)?nèi)表面裂紋自動(dòng)檢測(cè)系統(tǒng)研究[J].陶瓷學(xué)報(bào)，2016（4）：417?422.

[13] 和偉康，蘇向輝，馬瑛，等.烏魯木齊市產(chǎn)業(yè)結(jié)構(gòu)與城市空間形態(tài)耦合及響應(yīng)關(guān)系研究[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2017，38（9）：98?105.

[14] 張純，夏海山，宋彥.軌道交通帶動(dòng)下的城市形態(tài)演變：以北京為例[J].城市發(fā)展研究，2016，23（9）：107?112.

[15] 李在軍，張雅倩，胡美娟，等.基于最低成本?周期模型的昆明市空間形態(tài)演變研究[J].長(zhǎng)江流域資源與環(huán)境，2016，25（5）：708?714.