周再文 王 建 朱 恰 劉星雨 馬紫雯 高賢君

(長(zhǎng)江大學(xué) 地球科學(xué)學(xué)院, 湖北 武漢 430100)

0 引言

建筑物作為地理信息研究的比較重要的地理要素之一,是現(xiàn)代城市的主要組成,也是用來(lái)衡量城市發(fā)展的比較重要的指標(biāo)。建筑物信息的提取對(duì)城市規(guī)劃、地理信息系統(tǒng)的實(shí)際應(yīng)用、數(shù)字城市應(yīng)用等方面具有重要的意義和作用[1]。隨著遙感技術(shù)的快速發(fā)展,各種遙感影像的空間分辨率也越來(lái)越高。由于高分辨率遙感影像具有地物紋理信息豐富、成像光譜波段多、重訪時(shí)間短等特點(diǎn),其已經(jīng)成為地物信息獲取的重要的數(shù)據(jù)來(lái)源[2]。但是高分遙感影像在提高地物光譜特征,突出地物結(jié)構(gòu)、紋理信息的同時(shí),會(huì)帶來(lái)嚴(yán)重的異物同譜的問(wèn)題。實(shí)際工作中建筑物周圍的各種其他無(wú)關(guān)地物的干擾也會(huì)導(dǎo)致建筑物信息最終的提取結(jié)果與原始影像中的建筑物的形狀和邊界輪廓不一致。因此基于已有的建筑物信息提取方法對(duì)建筑物的邊界輪廓進(jìn)一步優(yōu)化對(duì)于提高建筑物精度具有十分重要的研究意義。

目前比較主流的建筑物信息提取方法,按照遙感影像提取建筑物的原理主要分成三類:

(1)基于區(qū)域分割的提取方法:先利用監(jiān)督、非監(jiān)督或面向?qū)ο蟮姆诸惙椒▽?duì)影像進(jìn)行初步的分類;再對(duì)建筑物所在的類別根據(jù)建筑物自身的幾何形狀、空間信息等特征提取具體的建筑物目標(biāo)[3]。黃慧萍等[4]提出了以多尺度影像分割與面向?qū)ο笥跋穹治鰹橹饕夹g(shù),利用樣本多邊形對(duì)象的成員函數(shù)建立訓(xùn)練區(qū),自動(dòng)提取大慶市城市綠地覆蓋信息,達(dá)到清查城市綠地的目的。

(2)結(jié)合輔助知識(shí)的提取方法:利用遙感圖像中建筑物的陰影、紋理、空間位置、幾何結(jié)構(gòu)等特征來(lái)提取建筑物,提高建筑物信息提取精度。高賢君等[5]設(shè)計(jì)了基于偏移陰影分析的建筑物樣本自動(dòng)提取方法,結(jié)合支持向量機(jī)(SVM)分類模型對(duì)影像進(jìn)行初提取并進(jìn)行地物分類,在進(jìn)行相應(yīng)的后處理操作之后,利用相交邊界陰影比率篩選出無(wú)陰影的非建筑物等措施獲取最終結(jié)果。

(3)基于邊緣和角點(diǎn)監(jiān)測(cè)與匹配的提取方法:使用邊緣監(jiān)測(cè)的算法得到建筑物的邊緣,再根據(jù)建筑物的空間關(guān)系對(duì)影像中提取的邊緣線段進(jìn)行分組,搜索其平行線,通過(guò)經(jīng)驗(yàn)知識(shí)組成符合建筑物空間結(jié)構(gòu)和輪廓的矩形,實(shí)現(xiàn)對(duì)建筑物的提取。Taejung Kim[6]等提出了根據(jù)現(xiàn)實(shí)中建筑物的屋頂邊緣多呈相互平行的線段特點(diǎn),利用投票策略選取直線最多的方向?yàn)榉课葸吘壸呦?達(dá)到提取建筑物信息的目的。陶超等[7]在Harris算子提出角點(diǎn)的基礎(chǔ)上,加上約束準(zhǔn)則檢測(cè)居民地的角點(diǎn),根據(jù)影像中角點(diǎn)的分布情況,構(gòu)建函數(shù)度量影像中每一個(gè)像素屬于建筑物的概率,最后利用二值分割的方法提取影像中建筑物的信息。

以上三種方法在進(jìn)行建筑物信息提取的過(guò)程中或多或少都會(huì)受到地物周圍的陰影、植被等噪聲干擾,導(dǎo)致提取結(jié)果出現(xiàn)邊界鋸齒化、不規(guī)則等問(wèn)題。本文在利用深度學(xué)習(xí)提取的高分辨遙感影像的建筑物信息基礎(chǔ)上,提出一種解決建筑物圖像邊界問(wèn)題的思路。

1 數(shù)據(jù)來(lái)源

基于深度學(xué)習(xí)提取建筑物信息的方法是以谷歌地球上開源的高分遙感影像數(shù)據(jù)集作為基礎(chǔ)的實(shí)驗(yàn)數(shù)據(jù),因谷歌地球上的開源的高分遙感影像國(guó)外地區(qū)的建筑物結(jié)構(gòu)比較簡(jiǎn)單,同時(shí)遙感影像數(shù)據(jù)比較清晰,故實(shí)驗(yàn)數(shù)據(jù)集多選取國(guó)外的高分遙感影像進(jìn)行實(shí)驗(yàn)?；趖ensorflow深度學(xué)習(xí)的方法構(gòu)建建筑物識(shí)別的模型,通過(guò)增加訓(xùn)練樣本的數(shù)據(jù)量,不斷地修正訓(xùn)練模型,提高測(cè)試樣本的提取精度。最后將提取精度最高的樣本進(jìn)行建筑物后處理驗(yàn)證,得到建筑物信息的初始提取結(jié)果。上述基于深度學(xué)習(xí)的建筑物信息提取的方法參考了朱光亞等[8]基于深度學(xué)習(xí)的遙感影像建筑物信息提取研究的論文和高揚(yáng)等[9]基于卷積神經(jīng)網(wǎng)絡(luò)的高分辨率遙感影像建筑物信息提取的論文的部分思想。

UNet網(wǎng)絡(luò)作為比較早的卷積神經(jīng)網(wǎng)絡(luò)之一,主要用于圖像的處理,由用于捕捉語(yǔ)義特征的編碼端,精準(zhǔn)定位的解碼端以及結(jié)合底層信息和高層信息的跳躍連接這三部分組成。本實(shí)驗(yàn)中,輸入影像的影像塊大小為512×512,從而構(gòu)建特征金字塔(實(shí)驗(yàn)的數(shù)據(jù)集是新西蘭克賴斯特徹奇地區(qū)的高分辨率遙感影像)。本次實(shí)驗(yàn)的環(huán)境為:i7-9700k處理器,32 GB內(nèi)存,NVIDIA RTX2080Ti 8G顯卡,GPU加速庫(kù)采用CUDA10.0。DeepingLearn框架以Keras作為前端,TensorFlow為后端。圖1、圖2的(a)和(b)是部分原始遙感影像和利用UNet神經(jīng)網(wǎng)絡(luò)提取的灰度圖二值化的結(jié)果。

由圖1和圖2的(b)和(c)的對(duì)比,可以很明顯地看出,經(jīng)過(guò)后處理驗(yàn)證的優(yōu)化提取處理,錯(cuò)分建筑物的問(wèn)題得到了部分解決,提高了建筑物特征信息的提取精度。但是如圖1(c)和圖2(c)所示,部分建筑物的邊界輪廓還是存在不規(guī)整問(wèn)題,另有部分建筑物內(nèi)部存在細(xì)小孔洞等問(wèn)題,為此,需要繼續(xù)對(duì)其進(jìn)行建筑物邊界輪廓的優(yōu)化。

圖1 遙感影像YX1提取優(yōu)化對(duì)比

圖2 遙感影像YX2提取優(yōu)化對(duì)比

2 建筑物輪廓優(yōu)化原理

對(duì)在基于深度學(xué)習(xí)提取建筑物特征和改進(jìn)的馬爾科夫隨機(jī)場(chǎng)處理錯(cuò)分建筑物的問(wèn)題的基礎(chǔ)上所得到的提取結(jié)果仍然存在建筑物邊界不規(guī)整、邊界鋸齒化嚴(yán)重的問(wèn)題,本文提出的解決思路是:首先是利用目前最常用的線要素壓縮算法——垂距法[10]對(duì)提取到的建筑物輪廓?jiǎng)h除邊界輪廓的冗余點(diǎn)并對(duì)建筑物輪廓進(jìn)行多邊形擬合處理。然后通過(guò)對(duì)多邊形擬合之后的輪廓利用最小外接矩形的擬合算法[11],得到最小外接矩形的四個(gè)頂點(diǎn)坐標(biāo),最后再對(duì)擬合之后的最小外接矩形的長(zhǎng)寬比、輪廓面積、擬合出來(lái)的多邊形點(diǎn)數(shù)做篩選,選擇最合適的外接矩形輪廓作為最后邊界優(yōu)化的最終結(jié)果,使通過(guò)優(yōu)化后的邊界輪廓盡量與遙感影像原圖上的建筑物保持一致,提高研究后期進(jìn)行建筑物變化檢測(cè)的最終精度。

2.1 基于垂距法的多邊形擬合

經(jīng)過(guò)前期深度學(xué)習(xí)提取的建筑物邊界一般是由許多連續(xù)的像素點(diǎn)組成,會(huì)出現(xiàn)大量相鄰的點(diǎn)位于同一條直線上產(chǎn)生數(shù)據(jù)冗余,也是導(dǎo)致提取的建筑物邊界輪廓出現(xiàn)鋸齒化的一個(gè)原因。本文利用地理信息原理中常見(jiàn)的矢量數(shù)據(jù)壓縮的方法——垂距法來(lái)解決這個(gè)問(wèn)題,即將提取出來(lái)的建筑物輪廓進(jìn)行冗余點(diǎn)去除,將建筑物的邊界輪廓擬合成相對(duì)規(guī)則的多邊形，如圖3所示。

圖3 垂距法處理過(guò)程示意圖

垂距法的算法思想[10]是:每一次按照一定順序取曲線上連續(xù)的三個(gè)點(diǎn),計(jì)算中間的點(diǎn)到其余兩點(diǎn)連線的垂直距離d,并將d與給定的限差(閾值)D比較。求得的垂直距離d小于給定的閾值D時(shí),則去掉中間點(diǎn);反之,則中間點(diǎn)保留。然后按照順序依次取下三個(gè)點(diǎn)繼續(xù)處理,直到遍歷整條曲線。多邊形擬合原理的示意圖如圖4所示,具體過(guò)程如下:

圖4 垂距法處理結(jié)果對(duì)比圖[10]

(1)取任意建筑物輪廓邊界上的一段邊界P1～P4,分別連接P1P3和P3P4;

(2)按順序取P1P2P3前三個(gè)點(diǎn),計(jì)算中間點(diǎn)P2到線段P1P3的距離d;

(3)比較中間點(diǎn)到其余兩點(diǎn)的距離d與給定的閾值D的大小,若d≤D,則去掉P2點(diǎn),P1～P3線段處理結(jié)束。

不過(guò)，盡管短周期存在周期及基數(shù)壓力，空調(diào)依然是最為穩(wěn)固的家電細(xì)分板塊，格力美的的雙寡頭地位穩(wěn)固，維持市場(chǎng)大體穩(wěn)定，也讓行業(yè)均價(jià)處于緩慢提升通道。不過(guò)從競(jìng)爭(zhēng)格局來(lái)看，二線品牌還能保持一定市場(chǎng)空間，但其它雜牌生存空間則日益狹小。

(4)若d≥D,保留中間點(diǎn)P2,然后沿著直線繼續(xù)依次對(duì)剩下的曲線進(jìn)行處理,直到整條線處理結(jié)束。

(5)將經(jīng)過(guò)上面四個(gè)步驟保留下來(lái)的點(diǎn)依次連接起來(lái),得到的曲線就是通過(guò)垂距法對(duì)建筑物輪廓擬合的結(jié)果。

2.2 基于建筑物最小外接矩形的擬合

現(xiàn)實(shí)生活中建筑物是比較有規(guī)律的人造工程產(chǎn)物,建筑風(fēng)格和外形千差萬(wàn)別，但是,就建筑物的輪廓而言,大多數(shù)建筑物的邊界輪廓為矩形或者矩形的組合[11-12]。我們?cè)谕ㄟ^(guò)處理基于垂距法對(duì)建筑物輪廓進(jìn)行多邊形擬合后,可以繼續(xù)基于建筑物的最小外接矩形來(lái)實(shí)現(xiàn)對(duì)建筑物邊界輪廓的進(jìn)一步擬合逼近。

具體處理過(guò)程如下:

(1)首先利用輪廓監(jiān)測(cè)函數(shù)將多邊形擬合的邊界輪廓監(jiān)測(cè)出來(lái)并進(jìn)行保存。

(2)利用第一步得到的邊界輪廓點(diǎn)集,采用Sklansky算法,得到邊界輪廓點(diǎn)集的凸包。

(3)利用第二步得到的建筑物邊界的凸包,根據(jù)旋轉(zhuǎn)卡尺算法(Rotating Calipers),生成最小面積外接矩形。

(4)通過(guò)對(duì)滿足長(zhǎng)寬比、輪廓面積、擬合出的多邊形定點(diǎn)數(shù)三個(gè)條件的最小外接矩形進(jìn)行篩選,對(duì)滿足條件的輪廓畫出輪廓擬合的最小外接矩形,輸出處理后二值化影像。

3 實(shí)驗(yàn)結(jié)果分析

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

通過(guò)本文提供的建筑物輪廓優(yōu)化的方法,對(duì)遙感影像YX1和YX2進(jìn)行進(jìn)一步優(yōu)化,結(jié)果如圖5和圖6所示。

圖5 遙感影像YX1逐級(jí)優(yōu)化結(jié)果對(duì)比圖

圖6 遙感影像YX2逐級(jí)優(yōu)化結(jié)果對(duì)比圖

通過(guò)圖7和圖8的局部?jī)?yōu)化放大圖可以很明顯看出，本文所介紹的建筑物輪廓優(yōu)化方法對(duì)于部分建筑物邊界輪廓有一定的優(yōu)化效果,能夠較大限度地逼近實(shí)際建筑物邊界輪廓。但是對(duì)于較復(fù)雜的建筑物(如邊界輪廓為回字形建筑和階梯形建筑)邊界輪廓處理上還存在問(wèn)題。本文所介紹的方法只能處理由比較簡(jiǎn)單矩形組合的建筑物輪廓,有一定的應(yīng)用局限性。

圖7 遙感影像YX1局部?jī)?yōu)化放大圖

圖8 遙感影像YX2局部?jī)?yōu)化放大圖

3.2 實(shí)驗(yàn)精度評(píng)價(jià)

傳統(tǒng)的影像分類方法經(jīng)常采用總體精度(overall accuracy, OA)、平均精度(Average Accuracy,AA)和Kappa系數(shù)作為評(píng)價(jià)指標(biāo)[13]。本文選用其中的OA指標(biāo)和混淆矩陣中的用戶精度的概念對(duì)本試驗(yàn)的結(jié)果進(jìn)行精度評(píng)定。精度評(píng)定結(jié)果如表1所示。

表1 實(shí)驗(yàn)結(jié)果精度評(píng)價(jià)表

表1顯示,通過(guò)深度學(xué)習(xí)提取的建筑物特征信息經(jīng)過(guò)建筑物信息后處理驗(yàn)證,消除部分錯(cuò)分建筑物問(wèn)題和輪廓優(yōu)化,使最終建筑物特征信息提取的精度有了明顯提升,實(shí)現(xiàn)了建筑物信息的精提取。但是當(dāng)建筑物顏色和地面顏色很相近(如遙感影像YX2的情況),部分建筑物的信息經(jīng)過(guò)卷積神經(jīng)網(wǎng)絡(luò)提取和后處理之后,仍然存在相應(yīng)的誤差,導(dǎo)致建筑物用戶精度稍低。

4 結(jié)束語(yǔ)

本文在深度學(xué)習(xí)提取建筑物信息基礎(chǔ)上,提出一種優(yōu)化建筑物圖像邊界輪廓問(wèn)題的思路。對(duì)于深度學(xué)習(xí)提取建筑物圖像結(jié)果中出現(xiàn)的鋸齒化、建筑物邊界不規(guī)整等問(wèn)題,提出了一種利用垂距法進(jìn)行建筑物邊界輪廓的多邊形擬合,再繼續(xù)利用初提取的結(jié)果進(jìn)行建筑物輪廓的最小外接矩形選取,最后得到較符合實(shí)際的建筑物邊界輪廓。通過(guò)實(shí)驗(yàn)驗(yàn)證,本文所介紹的方法在對(duì)深度學(xué)習(xí)提取高分遙感影像的建筑物結(jié)果的后處理過(guò)程中,對(duì)建筑物的邊界輪廓進(jìn)行的擬合逼近,使得提取后的結(jié)果與實(shí)際影像中的建筑物的輪廓比較接近,為建筑物變化監(jiān)測(cè)提供可靠的基礎(chǔ)數(shù)據(jù)。這種方法為通過(guò)高分遙感影像提取建筑物信息的后處理的邊界優(yōu)化提供了一種可行的新思路。