朱添翼，范 強(qiáng)，杜漫飛

(遼寧工程技術(shù)大學(xué)測(cè)繪與地理科學(xué)學(xué)院，遼寧 阜新 123000)

?

一種遙感圖像建筑物直線特征提取算法

朱添翼，范 強(qiáng)，杜漫飛

(遼寧工程技術(shù)大學(xué)測(cè)繪與地理科學(xué)學(xué)院，遼寧 阜新 123000)

隨著遙感技術(shù)的不斷發(fā)展，高效、快速的遙感影像處理成為人們不斷研究的熱點(diǎn)。作為內(nèi)業(yè)處理的基礎(chǔ)，人們對(duì)影像直線特征提取的重視不斷加強(qiáng)。本次試驗(yàn)采用遼寧科技大學(xué)1∶2000正射影像圖，針對(duì)教學(xué)區(qū)內(nèi)建筑物、綠化區(qū)和道路混合的特點(diǎn)，提出了一種自動(dòng)識(shí)別并高效提取建筑物直線特征的算法。該算法將直線擬合和Hough變換進(jìn)行結(jié)合，加入角度先驗(yàn)方法。試驗(yàn)結(jié)果表明：該算法的直線提取性能優(yōu)于傳統(tǒng)Hough變換，能很好地適應(yīng)教學(xué)區(qū)域，對(duì)建筑物的識(shí)別能力、直線提取效率及精度有很大的提高。

Hough變換；直線擬合；教學(xué)區(qū)域；建筑物；直線提取

遙感技術(shù)的增強(qiáng)使得航測(cè)和遙感衛(wèi)星提供的影像快速增多，空間、光譜與時(shí)間分辨率不斷提高，為測(cè)繪學(xué)科提供了豐富的基礎(chǔ)信息數(shù)據(jù)源。建筑物的特征識(shí)別和提取是遙感影像處理的基礎(chǔ)之一，在目標(biāo)跟蹤[1]、圖像處理[2]和目標(biāo)辨識(shí)[3]等領(lǐng)域都具有重要意義。如何高效地利用遙感影像進(jìn)行建筑物特征提取，已成為國內(nèi)外遙感領(lǐng)域研究的熱點(diǎn)。

直線作為遙感影像的基礎(chǔ)特征之一，對(duì)其進(jìn)行識(shí)別可反映影像的大量高層信息。對(duì)建筑物進(jìn)行識(shí)別、提取和三維重建時(shí)，直線是其最重要、最穩(wěn)定的特征。直線提取效果直接影響后續(xù)工作的準(zhǔn)確性。對(duì)于直線提取，人們提出了許多算法：如Hough變換及其改進(jìn)方法、相位編組法、鏈碼法、最小特征根分析法和主元分析法(PCA)[4-8]等。其中，Hough變換及其改進(jìn)方法是提出較早、最常用且應(yīng)用最為廣泛的算法，與其他4種方法相比，其抗噪能力強(qiáng)，直線識(shí)別能力強(qiáng)；但在復(fù)雜地形下，如教學(xué)區(qū)域，由于道路和綠化區(qū)的影響，應(yīng)用Hough變換將會(huì)把大量密集且不屬于同一直線上的短直線誤識(shí)別成同一條直線，從而產(chǎn)生錯(cuò)誤。

本文針對(duì)Hough變換缺點(diǎn)，研究并提出一種新算法。該方法通過對(duì)Hough變換進(jìn)行改進(jìn)，加入直線擬合[9]和角度先驗(yàn)方法[10]，解決了Hough變換生成大量短直線和弧線的問題，并使其對(duì)教學(xué)區(qū)域內(nèi)建筑物的識(shí)別能力和提取效率得到了很大的提高。

1 直線提取基本原理

1.1 提取流程

試驗(yàn)先對(duì)影像進(jìn)行預(yù)處理(如灰度化等)，建立降噪濾波模型降噪；接著對(duì)處理結(jié)果進(jìn)行初步直線檢測(cè)和直線擬合，并提取建筑物所在最小矩形范圍；最后對(duì)提取區(qū)域使用改進(jìn)后的Hough變換進(jìn)行直線提取并輸出結(jié)果。提取流程如圖1所示。

圖1 提取流程

1.2 去噪模型

由于影像獲取流程中外界條件和相機(jī)本身傳感器限制的影響，每張遙感影像圖都會(huì)存在噪聲。噪聲的存在將對(duì)影像特征提取產(chǎn)生很大的影響，去噪是遙感影像處理不可或缺的步驟。本次試驗(yàn)采用高斯濾波模型[11-12]，其濾波函數(shù)如下

(1)

式中，σ為方差，也可稱其為模型的平滑尺度，它決定了高斯函數(shù)的寬度。σ過大會(huì)造成影像細(xì)節(jié)模糊和圖像輪廓變形等問題，經(jīng)多次試驗(yàn)發(fā)現(xiàn)，式(1)中取σ=0.1并建立5×5的高斯去噪聲模型對(duì)試驗(yàn)影像的噪聲去除效果最佳。

1.3 邊緣檢測(cè)與直線提取

邊緣檢測(cè)：直線特征提取前應(yīng)對(duì)影像進(jìn)行邊緣檢測(cè)，邊緣檢測(cè)的效果直接影響后續(xù)直線提取的結(jié)果。由于Canny算子擁有較好的檢測(cè)精度和高信噪比[13]，挑選其作為邊緣檢測(cè)算法加入試驗(yàn)。

直線提?。杭碒ough變換，其原理是利用圖像空間和Hough參數(shù)空間的點(diǎn)線對(duì)偶性，把圖像空間中的點(diǎn)線關(guān)系轉(zhuǎn)換到參數(shù)空間中，將邊緣提取結(jié)果中有一定關(guān)系的像元進(jìn)行累加，設(shè)置直線判斷閾值，超過閾值則判斷為一條直線并輸出[14-16]。

Hough變換公式如下

ρ=xcosθ+ysinθ0≤θ<π

(2)

式中，ρ為原點(diǎn)到直線交點(diǎn)的距離；θ為垂直于圖像中直線的角度。

Hough變換原理如圖2所示。

圖2 Hough變換原理

2 建筑物直線提取算法

本文結(jié)合直線擬合和角度先驗(yàn)思路對(duì)傳統(tǒng)Hough變換算法進(jìn)行改進(jìn)，提高了其建筑物的識(shí)別能力和效率。算法步驟如下：

(1) 影像輸入與預(yù)處理，使用5×5的高斯去噪聲模型去噪，使用canny算子檢測(cè)邊緣信息(檢測(cè)結(jié)果包含不需要的道路和綠化區(qū)信息)。

(2) 對(duì)邊緣信息通過角度θ和距離ρ進(jìn)行直線判斷，并對(duì)每條短直線構(gòu)造直線函數(shù)y=kix+bi。

(3) 以每條直線l0為基礎(chǔ)，搜尋其余直線，設(shè)置條件k0×ki=1，保存滿足條件的直線。若存在兩條或兩條以上與l0垂直的直線，且存在非l0的直線與剩下的直線相垂直，則提取4個(gè)相交端點(diǎn)并存儲(chǔ)。

(4) 通過對(duì)比四點(diǎn)矩形區(qū)域的面積，最大矩形區(qū)域即需要保留的建筑物區(qū)域。劃分區(qū)域并通過4個(gè)邊界點(diǎn)由公式α=arctan[(y1-y2)/(x1-x2)]計(jì)算建筑區(qū)的兩個(gè)方位角α1和α2(α1<α2)，即后續(xù)Hough變換的角度先驗(yàn)值。

(5) 以劃分區(qū)域圖上灰度值為255的點(diǎn)為中心、α1和α2為搜索方向進(jìn)行同類點(diǎn)搜尋，設(shè)置閾值為70，保留大于閾值的直線并提取。

(6) 以提取結(jié)果圖上灰度值為255的點(diǎn)為中心，沿著角度α1和α2方向進(jìn)行同類點(diǎn)搜索，設(shè)置α1方向的閾值為40，設(shè)置α2方向的閾值為20，保留滿足閾值的起始點(diǎn)，刪除已搜索過的中心點(diǎn)。

(7) 連接起始點(diǎn)，輸出直線提取結(jié)果。

3 成圖分析

為了更加直觀地判斷方法之間的優(yōu)劣，本文運(yùn)用MATLAB軟件對(duì)提取結(jié)果進(jìn)行圖形建立，并與原影像圖進(jìn)行對(duì)比。對(duì)比結(jié)果如圖3—圖5所示。

圖3 原影像圖

圖4 傳統(tǒng)Hough變換提取

通過對(duì)圖3—圖5的對(duì)比分析可以得出，傳統(tǒng)Hough變換對(duì)試驗(yàn)區(qū)域內(nèi)建筑物的識(shí)別能力很差，提取出大量斷裂的短直線和弧線，整體效果較差。本文算法排除了周邊道路和綠化區(qū)域，對(duì)建筑物有很強(qiáng)的識(shí)別能力，由于加入了直線擬合和角度先驗(yàn)值，有效地抗拒了建筑物周邊和內(nèi)部的陰影與綠化噪聲，較好地提取了建筑物的基本骨架線，線與線之間交點(diǎn)明顯，提取效果較好。

圖5 本文算法提取圖

4 精度檢驗(yàn)

通過算法識(shí)別的測(cè)區(qū)建筑物輪廓線，南北走向提取較為完整，東西走向出現(xiàn)多處漏線現(xiàn)象。對(duì)于提取直線精度，通過圖上量取和實(shí)地丈量對(duì)相應(yīng)的建筑物10條外圍輪廓線進(jìn)行對(duì)比檢查，見表1。

表1 直線精度對(duì)比 m

通過表1可以得出，10條提取直線的檢查結(jié)果最大誤差為0.488 m，最小誤差為0.137 m，平均誤差為0.295 m。按最大誤差，滿足1∶5000地形圖的線長精度要求；按平均誤差，滿足1∶2000地形圖線長精度要求。

5 結(jié) 論

本文主要進(jìn)行了Hough變換、直線擬合與角度先驗(yàn)、成圖分析和精度檢驗(yàn)的研究與分析，得到如下結(jié)論：

(1) 直線擬合的加入，有效地識(shí)別了建筑物區(qū)域，對(duì)周邊綠化、道路有很強(qiáng)的抗性。

(2) 角度先驗(yàn)方法的加入使得建筑物提取線更接近基本骨架線，有效地解決了建筑區(qū)域內(nèi)綠化和陰影噪聲的影響。

(3) 保留直線起始點(diǎn)的方法，使得提取的直線間連接點(diǎn)明顯，為后續(xù)的影像匹配及目標(biāo)跟蹤等圖像處理步驟提供了很大的便捷。

(4) 提取精度較高，按最大誤差，可達(dá)到1∶5000地形圖的線長精度要求；按平均誤差，可達(dá)到1∶2000地形圖的線長精度要求。

該算法仍存在較多改進(jìn)之處，如對(duì)建筑物上較小矩形提取效果不佳，在建筑物上的南北走向線提取較完整，東西走向線中出現(xiàn)多處漏線等現(xiàn)象。

[1] 張建偉，張啟衡.基于直線邊緣特征匹配的擴(kuò)展目標(biāo)跟蹤方法[J].光學(xué)學(xué)報(bào), 2009, 29(10):2826-2831.

[2] 孫顯，王宏琦，張正，等.基于遙感圖像的城市道路自動(dòng)測(cè)繪方法研究[J].光學(xué)學(xué)報(bào), 2009, 29(1):86-92.

[3] WU Wenyu，YIN Dong，ZHANG Rong，et al.Bridge Recognition of Median-resolution SAR Images Using Pun Histogram Entropy[J].Chin. Opt. Lett., 2009, 7(7):572-575.

[4] 王競(jìng)雪，朱慶，王偉璽，等.結(jié)合邊緣編組的Hough變換直線提取[J].遙感學(xué)報(bào), 2014, 18(10):378-389.

[5] 袁廣林，薛模根，李從利，等.基于相位編組圖像分塊的快速變換直線檢測(cè)[J].計(jì)算機(jī)測(cè)量與控制, 2005, 13(7):647-649.

[6] 潘大夫，汪渤.基于邊緣方向的直線提取算法[J].北京理工大學(xué)學(xué)報(bào), 2008, 28(6):513-516.

[7] GURU D S，SHEKAR B H，NAGABHUSHAN P.A Simple and Robustline Detection Algorithm Based on Small Eigenvalue Analysis[J].Pattern Recognition Lett., 2004, 25(1):1-13.

[8] 馬克，阿爾貝托.特征提取與圖像處理[M].2版.石英，楊高波，譯.北京：電子工業(yè)出版社，2010:313-324.

[9] 李增兵，張曉帆，李瑞華，等.直線擬合快速實(shí)現(xiàn)的一種新算法[J].計(jì)算機(jī)測(cè)量與控制, 2006, 11(32):1524-1525.

[10] 郭靜鈺, 齊宏亮, 袁媛，等.先驗(yàn)圖像約束的有限角度CT圖像重建算法[J].核電子學(xué)與探測(cè)技術(shù), 2014,12(3):1421-1424.

[11] 戴激光，宋偉東，李建軍.一種有效的直線提取算法[J].測(cè)繪通報(bào), 2012(11):59-61.

[12] 王耀貴.圖像高斯平滑濾波分析[J].計(jì)算機(jī)與信息技術(shù), 2008,8(30):79-81.

[13] 王植, 賀賽先.一種基于Canny理論的自適應(yīng)邊緣檢測(cè)方法[J].中國圖象圖形學(xué)報(bào), 2004,8(11):65-70.

[14] 劉春閣.基于自適應(yīng)中值濾波和改進(jìn)變換的直線提取[J].測(cè)繪科學(xué), 2011, 36(6):267-269.

[15] 巫兆聰, 萬茜婷, 梁靜，等.粒度Hough變換及其在遙感影像直線檢測(cè)中的應(yīng)用[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2007, 10(4):860-863.

[16] 李明磊, 李廣云, 王力，等.3D Hough Transform在激光點(diǎn)云特征提取中的應(yīng)用[J].測(cè)繪通報(bào), 2015(2):29-33.

A Linear Feature Extraction Algorithm for Buildings in Remote Sensing Images

ZHU Tianyi，F(xiàn)AN Qiang，DU Manfei

(School of Mapping and Geographical Science， Liaoning Technical University， Fuxin 123000， China)

With the continuous development of remote sensing technology, efficient and fast remote sensing image processing has become a hot research topic. As the basis of internal processing, people pay more and more attention to the feature extraction. Using the 1∶2000 orthophoto map of the university of science and technology of Liaoning according to the teaching area, building area and road greening characteristics mixed, this experiment puts forward a kind of automatic recognition and extraction of building linear feature algorithm. The algorithm combines straight line fitting and Hough transform, and adds a priori method. The experimental results show that the algorithm is superior to the traditional Hough transform in line extraction, and it can well adapt to the teaching area, and has a great improvement in the recognition ability and the linear extraction efficiency of the building.

Hough transform; linear fitting; teaching area; building; straight line extraction

2016-11-03

朱添翼(1992—)，男，碩士生，研究方向?yàn)橛跋裥畔⑻崛±碚撆c方法。E-mail：494765448@qq.com

P237

A

0494-0911(2017)06-0036-04