李百壽，李靈芝，張 強(qiáng)

(1. 桂林理工大學(xué)測(cè)繪地理信息學(xué)院，廣西 桂林 541004； 2. 廣西空間信息與測(cè)繪重點(diǎn)實(shí)驗(yàn)室，廣西 桂林 541004)

隨著遙感技術(shù)的快速發(fā)展，光學(xué)遙感影像的分辨率不斷提高，采用多層次目標(biāo)識(shí)別[1]，柵格與矢量處理結(jié)合[2]，利用影像顏色、形狀、紋理、空間關(guān)系等圖像特征，進(jìn)行建筑物的自動(dòng)識(shí)別和標(biāo)繪一直是國(guó)內(nèi)外研究的熱點(diǎn)問(wèn)題[3]。其中如何充分利用建筑物的幾何特征[4]進(jìn)行建筑物角點(diǎn)檢測(cè)[5-6]、邊緣線檢測(cè)[7-8]、建筑物輪廓構(gòu)建[9-10]及主方向求取[11]，成為了研究中的重點(diǎn)。

張正通過(guò)分層次目標(biāo)過(guò)程模型描述建筑物輪廓，并通過(guò)角點(diǎn)檢測(cè)方法自動(dòng)標(biāo)繪出目標(biāo)矢量形狀，建筑物標(biāo)繪精度達(dá)到91.2%[1]，可用于多類城區(qū)建筑物的自動(dòng)標(biāo)繪。還可以通過(guò)機(jī)器學(xué)習(xí)，建筑物輪廓構(gòu)建采用兩個(gè)反向神經(jīng)網(wǎng)絡(luò)分別用來(lái)生成初始化建筑物輪廓和刪除無(wú)用輪廓信息，根據(jù)感知器內(nèi)柵格排列的情況訓(xùn)練神經(jīng)網(wǎng)絡(luò)，確定出建筑物輪廓描述的規(guī)則，實(shí)現(xiàn)建筑物聚集區(qū)域的標(biāo)繪[9]。

利用圖像邊緣檢測(cè)技術(shù)獲取建筑物的邊緣輪廓信息，然后通過(guò)對(duì)邊緣輪廓進(jìn)行篩選、優(yōu)化及填充等操作提取出了完整的建筑物邊界[2]。孫金彥結(jié)合SVM多類分割算法和數(shù)學(xué)形態(tài)學(xué)方法粗略得到建筑物的候選斑塊，通過(guò)高斯邊緣模板匹配的矢量化方法細(xì)化建筑物輪廓信息。矢量化結(jié)果平均準(zhǔn)確度為89%、完整度為98%、幾何形狀相似性為87%、整體質(zhì)量為85%[2]。呂鳳華基于模式識(shí)別和場(chǎng)景分析思想，通過(guò)HOG識(shí)別建筑物中心部分，再由紋理光譜融合特征剔除干擾地物，這種多層次特征結(jié)合的方法可細(xì)化建筑物的邊緣信息，其中建筑物正確檢出率96%以上，錯(cuò)檢率為12%[12]。

基于細(xì)化的建筑物標(biāo)繪方法精度較高，但較費(fèi)時(shí)，需要利用種子區(qū)域生長(zhǎng)、hough直線檢測(cè)技術(shù)，不能直接提取圖形的線寬等信息[13-14]。例如QI Min等提出了一種從單一的遙感影像中自動(dòng)提取建筑物的直線輪廓信息。首先提出一種獲取初始種子點(diǎn)區(qū)域生長(zhǎng)的方法，隨后以Canny邊緣圖輔助區(qū)域生長(zhǎng)策略選擇分離閾值去除非建筑物部分，最后屋頂被細(xì)化后并配用直線擬合的算法來(lái)獲得直線邊界，將直線輪廓通過(guò)邊界跟蹤提取[15]。程效軍提出一種綜合多種數(shù)據(jù)提取建筑物輪廓，通過(guò)建筑物同名角點(diǎn)將機(jī)載雷達(dá)點(diǎn)云與航空影像進(jìn)行配準(zhǔn)，根據(jù)融合影像對(duì)建筑物輪廓進(jìn)行精確提取[16]。多源遙感數(shù)據(jù)的綜合運(yùn)用雖然可以進(jìn)行精細(xì)化提取，但是需同時(shí)獲取高分影像數(shù)據(jù)和機(jī)載雷達(dá)點(diǎn)云數(shù)據(jù)，不僅增加了數(shù)據(jù)獲取難度且花費(fèi)成本較高。

本文基于分層次目標(biāo)提取的基本方法，提出一種最小外包矩形(MBR)約束下的高分光學(xué)影像直角建筑物標(biāo)繪方法。先通過(guò)CART分類樹(shù)結(jié)合面向?qū)ο蟮姆椒ù痔崛〗ㄖ锏妮喞畔ⅲ偻ㄟ^(guò)構(gòu)建角點(diǎn)方向決策器確定建筑物的角點(diǎn)及其順序，最后連接各個(gè)角點(diǎn)自動(dòng)標(biāo)繪出光學(xué)影像中的建筑物。

1 MBR約束下的建筑物提取與標(biāo)繪

1.1 建筑物提取與標(biāo)繪方法流程

MBR約束下的建筑物提取和標(biāo)繪技術(shù)流程如圖1所示，詳細(xì)步驟如下：

(1) 將高分辨率遙感影像通過(guò)面向?qū)ο蟮姆椒ㄟM(jìn)行多尺度影像分割，根據(jù)影像特征將同質(zhì)區(qū)域分割成多個(gè)子區(qū)域。

(2) 從分割后影像選取一定量有代表性的數(shù)據(jù)作為樣本對(duì)象特征數(shù)據(jù)，導(dǎo)出其屬性信息，通過(guò)SPSS分類模塊中的CART分類方法，確定各屬性的重要性。

(3) 選取重要性前三的屬性進(jìn)行訓(xùn)練，獲取CART分類樹(shù)分裂節(jié)點(diǎn)閾值，建立規(guī)則信息。

(4) 根據(jù)在SPSS中建立起的CART規(guī)則信息，在eCognition中重新建立軟件能夠讀取的分類規(guī)則，提取出建筑物所在區(qū)域。

(5) 通過(guò)形態(tài)學(xué)變換修復(fù)提取的建筑物區(qū)域，首先利用數(shù)學(xué)形態(tài)學(xué)處理步驟(4)提取的建筑物區(qū)域，去除提取的建筑物的“毛刺”“空洞”等問(wèn)題區(qū)域。并通過(guò)Candy算子進(jìn)行建筑物輪廓粗提取。

(6) 根據(jù)步驟(5)提取的建筑物輪廓，計(jì)算最小面積的MBR，并確定建筑物的主導(dǎo)方向。

(7) 計(jì)算建筑物各邊到MBR的一對(duì)直角邊的距離，利用聚類的方法把所得到的距離值分成N類，每一類用其平均值D表示距離直線的距離，N表示建筑物的邊數(shù)。

(8) 將N個(gè)距離值對(duì)應(yīng)的MBR邊平移D個(gè)單位以此來(lái)擬合建筑物邊緣線；計(jì)算各個(gè)直線的交點(diǎn)，以此來(lái)表示建筑物的角點(diǎn)或者位于邊界線上的點(diǎn)。

(9) 確定角點(diǎn)的行走方向，然后根據(jù)行走方向?qū)⒏鼽c(diǎn)按指定的順序連接起來(lái)形成的直角多邊形就是標(biāo)繪的最終結(jié)果。

圖1 建筑物提取與標(biāo)繪技術(shù)流程

1.2 CART決策樹(shù)結(jié)合多尺度分割提取建筑物區(qū)

1.2.1 CART決策樹(shù)

CART決策樹(shù)在遙感影像分類中得以有效地應(yīng)用。其基本思想是根據(jù)一定的分類規(guī)則，把遙感影像數(shù)據(jù)一分為二的向下劃分，最終得到屬性不同的子類。

決策樹(shù)中每一個(gè)子節(jié)點(diǎn)表示一個(gè)決策，不同的決策選擇影響著葉節(jié)點(diǎn)表示的決策結(jié)果。CART算法采用的是二元遞歸分割的方法將樣本集逐級(jí)劃分形成二叉樹(shù)。分割的評(píng)價(jià)策略為

(1)

式中，T為決策樹(shù)的某個(gè)節(jié)點(diǎn)；m為類別集因變量數(shù)目；Ci為類別；pi為類i在T中的相對(duì)概率。當(dāng)i在T中為傾斜時(shí)，gini(T)最小。當(dāng)gini(T)=0時(shí)，所有樣本屬于同一類，當(dāng)所有類在節(jié)點(diǎn)中表現(xiàn)為相同的概率時(shí)，gini(T)最大化。衡量出某個(gè)屬性的GINI指數(shù)后，可以得到ginisplit(T)數(shù)據(jù)

(2)

式(2)被稱為GINI信息增益，在CART中，選擇其中最小的信息增益作為結(jié)點(diǎn)劃分決策樹(shù)。訓(xùn)練數(shù)據(jù)集T中的屬性F將T分成T1和T2，給定劃分的GINI指標(biāo)

(3)

對(duì)于離散值屬性，在算法中遞歸的選擇該屬性產(chǎn)生最小GINI指標(biāo)的子集作為它的分裂子集。

1.2.2 融合CART決策樹(shù)和多尺度分割方法提取建筑物區(qū)

根據(jù)影像的特征，對(duì)影像進(jìn)行多尺度分割將不同質(zhì)區(qū)域分割成多個(gè)子區(qū)。從分割后的遙感影像子區(qū)選擇一定量的、有代表性的數(shù)據(jù)作為樣本。從eCognition中導(dǎo)出樣本數(shù)據(jù)屬性信息。利用SPSS分類模塊中的CART分類方法檢測(cè)各個(gè)屬性的重要性。根據(jù)選取重要性排在前三的屬性，獲取CART分類樹(shù)分裂節(jié)點(diǎn)的閾值。利用SPSS中CART功能，選擇相應(yīng)的屬性進(jìn)行訓(xùn)練，就可以直接獲取閾值。根據(jù)閾值做出CART分類樹(shù)的示意圖如圖2所示。

圖2 CART規(guī)則樹(shù)建立

根據(jù)在SPSS中建立起的CART規(guī)則信息，在eCognition中重新建立軟件能夠讀取的分類規(guī)則，進(jìn)行建筑物區(qū)域的粗提取。

1.3 建筑物的自動(dòng)化提取與形狀標(biāo)繪

1.3.1 建筑物MBR計(jì)算

在Matlab平臺(tái)，首先使用Candy算子對(duì)粗提取建筑物區(qū)域進(jìn)行輪廓提取。然后將建筑物輪廓信息轉(zhuǎn)化成輪廓點(diǎn)集，通過(guò)輪廓點(diǎn)集運(yùn)算，計(jì)算建筑物的MBR，詳細(xì)的計(jì)算步驟如下：

(1) 輸入輪廓點(diǎn)集Ei(xi,yi)，并計(jì)算輪廓的形心Center(p,q)，形心的x坐標(biāo)為

形心的y坐標(biāo)為

Center(p,q)、Ei(xi,yi)連線與x軸的夾角為a1。

(2) 找出Ei(xi,yi)中xi和yi的最小和最大值minx,miny,maxx,maxy，以(minx,miny)，(maxx,maxy)，(minx,maxy)，(maxx,miny)作為外接矩形的4個(gè)頂點(diǎn)確定一個(gè)矩形。計(jì)算矩形面積A1,將A1賦給min(Area)，頂點(diǎn)坐標(biāo)賦給

Rect=(minx,miny,maxx,maxy)

(4)

(4) 調(diào)用步驟(2)，計(jì)算旋轉(zhuǎn)β后的外接矩形的面積Ai，判斷(Ai)≤min(Area)是否成立;若成立,則min(Area)=Ai，同時(shí)將獲得的矩形4個(gè)頂點(diǎn)坐標(biāo)順時(shí)針旋轉(zhuǎn)β，更新

Rect=(minx′,miny′,maxx′,maxy′)

(5) 重復(fù)步驟(4)，直至程序結(jié)束，最后的矩形頂點(diǎn)坐標(biāo)集即構(gòu)成最小外包矩形。

1.3.2 輪廓分段擬合與優(yōu)化

提取建筑物的粗輪廓往往呈不規(guī)則的波浪形邊界，然而實(shí)際中建筑物邊界為規(guī)則直線。將不規(guī)則的輪廓轉(zhuǎn)成規(guī)則直線最基本常用的方法為直線擬合。設(shè)一組坐標(biāo)點(diǎn)Ni(xi,yi)∈E(i=1,2,…,M)之間的線性關(guān)系為：y=ax+b,通過(guò)最小二乘估計(jì)參數(shù)

(5)

建筑物邊界線由多個(gè)直線段組成，因此需要對(duì)建筑物粗輪廓進(jìn)行分段。分段的方法為將輪廓坐標(biāo)投影到該輪廓的主方向或副主方向上。

傳統(tǒng)直線擬合結(jié)果在確定建筑物角點(diǎn)時(shí)會(huì)出現(xiàn)很大偏差，擬合出的直線破壞了邊界線之間的正交關(guān)系和平行關(guān)系。為了解決此偏差，需要對(duì)最小二乘擬合的結(jié)果進(jìn)行優(yōu)化。本文在MBR約束條件下調(diào)整擬合直線斜率，使直線與主方向平行或正交，然后根據(jù)最小二乘法求得直線方程截距。設(shè)優(yōu)化后的直線方程為

y=a′x+b′

(6)

式中，α′=tanθ；b′=(∑yi-α′∑xi)/M。

1.3.3 交點(diǎn)方向決策器

擬合直線的交點(diǎn)并非全部位于建筑輪廓上，故先要設(shè)定一個(gè)半徑為R的圓，判斷圓域內(nèi)是否有輪廓片段數(shù)據(jù)。如果沒(méi)有，則不做決策，進(jìn)入下一個(gè)交點(diǎn)；如果有輪廓片段，則根據(jù)圓形區(qū)域內(nèi)輪廓片段的特征給該點(diǎn)的方向做出決策。最后將這些點(diǎn)按照順序連接，即可成功標(biāo)繪出建筑物。為此，本文設(shè)計(jì)了交點(diǎn)方向決策器。

交點(diǎn)方向決策器是一個(gè)以交點(diǎn)為圓心，半徑為R的圓形區(qū)域，具有決策該點(diǎn)方向的功能。為了讓決策器做出準(zhǔn)確判斷，首先定義出交點(diǎn)的方向關(guān)系。以建筑物主方向和副主方向?yàn)樽鴺?biāo)軸，交點(diǎn)坐標(biāo)Pt為坐標(biāo)原點(diǎn)，建立局部坐標(biāo)系如圖3所示，新建立的坐標(biāo)系把平面分成了4個(gè)區(qū)域，位于局部坐標(biāo)原點(diǎn)上方的區(qū)域定義為第一區(qū)，位于局部坐標(biāo)原點(diǎn)左邊的區(qū)域?yàn)榈诙^(qū)，位于局部坐標(biāo)原點(diǎn)下方的為第三區(qū)，位于局部坐標(biāo)原點(diǎn)右方的區(qū)域?yàn)榈谒膮^(qū)。

圖3 局部坐標(biāo)系

根據(jù)交點(diǎn)處輪廓的特征可以得出6種方向關(guān)系，如圖4所示，因此一個(gè)決策器可以根據(jù)輪廓片段特征做出6種方向關(guān)系的決策。圖4(a)表示第一區(qū)走向，(b)表示第四區(qū)走向，(c)表示第三區(qū)走向，(d)表示第二區(qū)走向，(e)表示主方向走向，(f)表示副主方向走向。

2 試驗(yàn)與結(jié)果分析

為了驗(yàn)證提出的基于MBR約束的光學(xué)影像直角建筑物標(biāo)繪方法的有效性，本文進(jìn)行了相關(guān)試驗(yàn)。試驗(yàn)1為從影像中截取了單一建筑物進(jìn)行提取和標(biāo)繪；試驗(yàn)2為對(duì)某一區(qū)域所有建筑物進(jìn)行標(biāo)繪；最后計(jì)算標(biāo)繪后得到的建筑物的周長(zhǎng)和面積，進(jìn)行精度評(píng)價(jià)。

圖4 決策器方向模式

2.1 試驗(yàn)1

本文所用研究數(shù)據(jù)為廣西某地高分遙感影像截取的單一建筑物，影像大小為219×219像素，試驗(yàn)環(huán)境為Matlab2013a 語(yǔ)言編程環(huán)境，使用計(jì)算機(jī)處理器為Intel(R) Xeon(R)CPU E5645，處理速度為2.40 GHz，內(nèi)存為4 GB。利用Matlab2013a編程平臺(tái)計(jì)算該建筑物的MBR，如圖5所示。

圖5 原始影像及其最小MBR計(jì)算

該MBR由直線L12、L23、L34、L41組成，根據(jù)點(diǎn)到直線的距離公式，計(jì)算輪廓上的坐標(biāo)點(diǎn)Ei(xi,yi)到直線L12、L23的距離，如圖6所示。

由圖6可知，曲線的走勢(shì)具有一定的穩(wěn)定性，即可以利用距離將輪廓分段。本文根據(jù)聚類波峰波谷點(diǎn)實(shí)現(xiàn)輪廓分段。當(dāng)某個(gè)點(diǎn)的二次差分值大于0，則該點(diǎn)表示波峰；當(dāng)某個(gè)點(diǎn)的二次微分值小于0，則該點(diǎn)表示波谷。波峰波谷聚集情況如圖7所示。

由圖7(b)可知，波峰波谷聚集在4處，而且每一處的波峰波谷分布在同一個(gè)距離層次上，說(shuō)明這4處的波峰波谷代表了4條平行于L12的直線；另外，共有波峰波谷聚集在4處，其中有2處分布在同一個(gè)距離層次上，說(shuō)明這2處的輪廓出自同一條直線，因此表明組成建筑物的邊界中有3條邊平行于L23。

根據(jù)波峰波谷聚集的情況，以及在距離層次上的分布，計(jì)算出L12和L23平移的距離情況。最終的平移距離以每一個(gè)段中距離的平均值來(lái)表示，見(jiàn)表1。

圖6 輪廓坐標(biāo)到MBR一對(duì)直角邊的距離

圖7 根據(jù)波峰波谷的輪廓分段

上一個(gè)過(guò)程中共獲得7條直線，分別代表著建筑物的7條邊。而直線的交點(diǎn)則表示建筑物的角點(diǎn)或建筑物邊上的點(diǎn)。7條直線彼此正交，產(chǎn)生了12個(gè)交點(diǎn)，交點(diǎn)的分布情況如圖8所示。

圖8 MBR擬合直線的交點(diǎn)

計(jì)算輪廓片段的曲率圓心平均坐標(biāo)，交點(diǎn)坐標(biāo)曲率圓的結(jié)果如圖9所示。

圖9 輪廓片段曲率圓

從圖9可以看到，非建筑物角點(diǎn)的交點(diǎn)，其區(qū)域內(nèi)輪廓片段的平均曲率圓的圓心位于該交點(diǎn)局部坐標(biāo)系的一對(duì)主方向軸上或在很小范圍內(nèi)波動(dòng)。若位于副主方向軸上(或附近)，則其行走方向?yàn)橹鞣较蜃呦?；若位于主方向軸線上(或附近)則其行走方向?yàn)楦敝鞣较颍蝗魣A心落于該點(diǎn)相應(yīng)的象限內(nèi)則該點(diǎn)的行走方向?yàn)樵撓笙薜男凶叻较?。根?jù)各交點(diǎn)對(duì)應(yīng)的平均曲率圓的圓心位置確定各交點(diǎn)的行走方向，如圖10所示。

圖10 交點(diǎn)行走方向確定

根據(jù)各角點(diǎn)的行走方向，把這些點(diǎn)依次連接起來(lái)形成一個(gè)封閉的正交多邊形，這就是基于MBR標(biāo)繪方法的建筑物標(biāo)繪結(jié)果，如圖11所示。

圖11 基于MBR標(biāo)繪結(jié)果

將標(biāo)繪結(jié)果與人工標(biāo)繪結(jié)果作對(duì)比，以人工獲取的形狀作為真實(shí)值，求出基于MBR標(biāo)繪的結(jié)果精度(見(jiàn)表2)。將標(biāo)繪結(jié)果的面積與周長(zhǎng)作為評(píng)價(jià)指標(biāo),面積的相對(duì)誤差為3.5%，周長(zhǎng)的相對(duì)誤差為4.6%。

表2 基于最小MBR標(biāo)繪的誤差

2.2 試驗(yàn)2

試驗(yàn)2中選取一幅包含多個(gè)直角建筑的高分遙感影像，采用本文的技術(shù)方法進(jìn)行建筑物標(biāo)繪，并將標(biāo)繪的結(jié)果與人工標(biāo)繪結(jié)果作對(duì)比,如圖12、圖13 所示。

圖12 直角建筑物自動(dòng)標(biāo)繪結(jié)果

圖13 本文方法標(biāo)繪結(jié)果與人工標(biāo)繪真實(shí)結(jié)果對(duì)比

從圖13可以看出，選取的9個(gè)建筑物整體的標(biāo)繪結(jié)果與實(shí)際值相差不多，其中面積自動(dòng)標(biāo)繪誤差比周長(zhǎng)的自動(dòng)標(biāo)繪誤差稍大。然而面積自動(dòng)標(biāo)繪誤差比周長(zhǎng)的自動(dòng)標(biāo)繪誤差稍大，引起這個(gè)現(xiàn)象的主要原因是因?yàn)樵几叻直媛视跋裰幸恍┓墙ㄖ锛?xì)節(jié)干擾信息很容易被誤提取，從而引起粗輪廓提取精度降低。而周長(zhǎng)的精度相對(duì)較高的原因是本文采取了MBR約束條件，保證了標(biāo)繪過(guò)程中直線邊界的幾何特征與實(shí)際建筑物特征相符。最后統(tǒng)計(jì)全部9個(gè)建筑物自動(dòng)標(biāo)繪的面積和周長(zhǎng)的平均誤差，得出面積的平均精度為93.3%，周長(zhǎng)的平均精度為96.1%。

3 結(jié) 語(yǔ)

本文基于建筑物的直角特征提出了一種高分辨率遙感影像中直角建筑物提取、標(biāo)繪的新方法。主要針對(duì)直角建筑物進(jìn)行標(biāo)繪，并將MBR標(biāo)繪方法用于實(shí)際建筑物提取與標(biāo)繪中。本方法的流程中，CART分類樹(shù)與多尺度分割結(jié)果包含誤差，通過(guò)對(duì)建筑物輪廓分段擬合優(yōu)化，可以有效去除粗提取建筑物結(jié)果中的邊緣“毛刺”，并由建筑物輪廓線特征提取出建筑物角點(diǎn)，利用角點(diǎn)決策器確定各個(gè)角點(diǎn)的行走方向，通過(guò)對(duì)建筑物角點(diǎn)的精細(xì)化提取，成功標(biāo)繪出建筑物信息，提高了高分遙感影像中建筑物提取的精度。

試驗(yàn)結(jié)果表明，本文所述標(biāo)繪方法僅采用高分光學(xué)數(shù)據(jù)而無(wú)需激光雷達(dá)數(shù)據(jù)，研究的方法復(fù)雜程度低，建筑物提取速度快，易于實(shí)現(xiàn)。建筑物提取過(guò)程中影像質(zhì)量、干擾地類等影響因素較多，目前光學(xué)影像中建筑物提取模型存在擾動(dòng)性，邊緣易產(chǎn)生錯(cuò)分和漏分現(xiàn)象，通過(guò)幾何約束條件可以有效減少提取誤差。因此，如何通過(guò)幾何約束更加快速有效地提取并標(biāo)繪出不同形狀的建筑物，提高算法的穩(wěn)健性為后續(xù)工作的研究方向。

參考文獻(xiàn)：

[1] 張正,王宏琦,寧忠磊,等.高分辨率遙感圖像中城區(qū)建筑物自動(dòng)測(cè)繪方法研究[J].國(guó)外電子測(cè)量技術(shù),2010,29(6):49-53.

[2] 孫金彥,黃祚繼,周紹光,等.高分辨率遙感影像中建筑物輪廓信息矢量化[J].遙感學(xué)報(bào),2017,21(3):396-405.

[3] BIOUCASDIAS J M，PLAZA A，CAMPSVALLS G，et al.Hyperspectral Remote Sensing Data Analysis and Future Challenges[J].IEEE Geoscience and Remote Sensing Magazine，2013，1(2): 6-36.

[4] 張祖勛,張劍清,胡翔云.基于物方空間幾何約束最小二乘匹配的建筑物半自動(dòng)提取方法[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2001,26(4):290-295.

[5] CUI S，YAN Q，REINARTZ P，et al.Complex Building Description and Extraction Based on Hough Transforma-tion and Cycle Detection[J].Remote Sensing Letters，2011，3(2): 151-159.

[6] TURKER M,KOC-SAN D.Building Extraction from High-resolution Optical Spaceborne Images Using the Integration of Support Vector Machine (SVM) Classific-ation，Hough Transformation and Perceptual Grouping[J].Applied Earth Observation and Geoinformation,2015,34(2):58-69.

[7] WANG J，YANG X，QIN X，et al.An Efficient Approach for Automatic Rectangular Building Extraction from Very High Resolution Optical Satellite Imagery[J].IEEE Geoscience and Remote Sensing Letters，2015，12(3): 487-491.

[8] FAZAN A J，DAL A P.Rectilinear Building Roof Contour Extraction Based on Snakes and Dynamic Programming[J].International Journal of Applied Earth Observation and Geoinformation，2013，25(8)：1-10.

[9] CHENG B，LIU Q，LI X，et al.Local Perception-based Intelligent Building Outline Aggregation Approach with Back Propagation Neural Network[J].Neural Processing Letters，2015，41(2): 273-292.

[10] CHAUDHURI D,KUSHWAHA N K,SAMAL A，et al.Automatic Building Detection from High-resolution Satellite Images Based on Morphology and Internal Gray Variance[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing，2016，9(5): 1767-1779.

[11] 程亮,龔健雅.LiDAR輔助下利用超高分辨率影像提取建筑物輪廓方法[J].測(cè)繪學(xué)報(bào)，2008，37(3)：391-393,399.

[12] 呂鳳華,舒寧,龔龑,等.利用多特征進(jìn)行航空影像建筑物提取[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2017,42(5):656-660.

[13] 吳秀蕓,李艷,周華.基于角點(diǎn)檢測(cè)的建筑物輪廓矢量化方法[J].遙感信息,2011,17(5):95-99.

[14] HUANG X，ZHANG L.Morphological Building/Shadow Index for Building Extraction from High-resolution Imagery over Urban Areas[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing，2012，5(1): 161-172.

[15] QI Min，JIA Jiankui.Automatic Extraction of the Building Rectilinear Rooftops Based on Region over Growing[J].Journal of Computers,2014，9(10)：2365-2370.

[16] 程效軍,程小龍,胡敏捷,等.融合航空影像和LIDAR點(diǎn)云的建筑物探測(cè)及輪廓提取[J].中國(guó)激光,2016,43(5):253-261.