摘   要: 為了快速有效地從遙感圖像中提取目標建筑物，采用小波分析和Markov隨機場(MRF)相結(jié)合的方法進行遙感圖像建筑物目標分割。首先將小波分解得到圖像的多尺度序列作為各個尺度特征場的觀測值，采用高斯混合模型建模特征場，以MRF模型作為標記場先驗概率分布模型,通過EM算法迭代使得參數(shù)估計和圖像分割交替進行，最后采用模板匹配檢測建筑物目標的位置。選擇多幅遙感圖片進行實驗，結(jié)果表明,采用分解層數(shù)為3的Haar小波，類別數(shù)為2，MLL模型勢函數(shù)β為1時，該方法能夠完成復雜背景下的建筑物目標分割并能對規(guī)則目標進行檢測。
關(guān)鍵詞: 遙感影像; 圖像分割; 馬爾可夫隨機場; 建筑物檢測; 模板匹配

    建筑物是城市遙感影像的典型地物之一，建筑物等人造目標的分割提取是遙感圖像分析和處理的關(guān)鍵。對建筑物提取的方法有很多:(1)利用建筑物的顏色紋理信息[1-2]，如參考文獻[1]采用圖像子塊灰度的標準差和直方圖的熵作為模糊C-均值聚類的特征矢量用于建筑群的粗略分割，參考文獻[2]將Gabor紋理塊的遙感影像對象模型方法應(yīng)用于遙感影像城鎮(zhèn)建筑物的提取；(2)從建筑物形狀特征出發(fā)[3-4]，如田昊等將多個先驗形狀競爭模型引入水平集方法中，在標記函數(shù)的指導下，利用先驗形狀能量來約束曲線的演化，對圖像分割的同時完成建筑物的檢測和提取;(3)基于建筑物幾何直線特征的提取方法[5-6]。
    本文介紹無監(jiān)督環(huán)境下基于馬爾可夫隨機場MRF(Markov Random Field)的可見光遙感圖像中建筑物目標分割。利用建筑物顏色信息，從統(tǒng)計學角度出發(fā),將影像各個像素顏色值視為具有一定概率分布的隨機變量， MRF表示影像像素間的相互關(guān)系，結(jié)合影像的多尺度特性，采用貝葉斯規(guī)則,通過數(shù)學模型進行影像分割。對分割結(jié)果進行數(shù)學形態(tài)學處理后，利用相關(guān)性和常見建筑物形狀模板進行匹配，可以檢測出建筑物的位置。
1 小波域影像建模
1.1小波變換
    單尺度的基于MRF的圖像分割方法中都有一個基本假設(shè)，即圖像是一個平穩(wěn)的隨機過程，這種假設(shè)過于嚴格而與實際的圖像特性出入較大。遙感數(shù)字圖像通常不是平穩(wěn)的隨機信號，因此將MRF模型與小波分解結(jié)合起來進行圖像分析，可以更好地描述圖像的統(tǒng)計特點，從而很好地呈現(xiàn)其非平穩(wěn)特性。
    對圖像進行離散小波變換,經(jīng)過一次小波分解，獲取4個頻帶(LL、LH、HL和HH)的小波系數(shù)，低頻小波系數(shù)為原始圖像的一個近似，高頻小波系數(shù)則在不同方向上捕獲了圖像的奇異性，對應(yīng)位置的小波系數(shù)構(gòu)成系數(shù)向量,表示該尺度、該位置處圖像的觀測特征。經(jīng)過L層小波分解后，遙感數(shù)字圖像分解為:

1.3 多尺度圖像特征場統(tǒng)計模型
    小波分解后的各尺度特征場是一個向量場，觀測到的小波系數(shù)向量圖像是該特征場的一個實現(xiàn)。假設(shè)圖像的各個位置的像素是獨立同分布的，即
 
1.5 影像分割的MAP估計
       由于隨機場X⁽ⁿ⁾是MRF,具有正概率性和Markov性，不同分辨率尺度的標記場之間也滿足Markov性，標記場和特征場的聯(lián)合分布可以表示為:
 

3 實驗結(jié)果與分析
　選擇512×512像素大小的不同的實驗遙感圖片，使用Haar小波基進行小波分解，層數(shù)L=3，類別數(shù)M=2，MLL模型勢函數(shù)β=1。各個分辨率的初始權(quán)重分別設(shè)為5、10、15、20。分割結(jié)果共得到包括原始分辨率影像在內(nèi)的4個分辨率影像數(shù)據(jù)，這里僅給出尺度n=0(原始分辨率)的分割結(jié)果，如圖2所示。

    實驗結(jié)果證明，圖2(b)、圖2(d)的分割效果較好， 圖2(f)并未達到理想分割。圖2(a)中建筑物目標較少，從圖2(b)可看出所有建筑物被分割出來，但建筑物低矮部分受到陰影遮擋，建筑物目標并不完整。圖2(c)中建筑物目標密集而排列無規(guī)則，從圖2(d)可看出大多數(shù)建筑物被分割出來，建筑物和道路有明顯的分界。 圖2(e)中建筑物目標不是很復雜，圖2(f)是其分割結(jié)果，從中可看出分割出的部分建筑物被“擴大化”，中間位置的建筑物周圍部分不是建筑物而被當成建筑物分割出來,這是因為被誤分的建筑物和真正建筑物是同質(zhì)性區(qū)域，反射性質(zhì)相似，造成了混淆。
    對圖2(b)進行形態(tài)學運算，將其和圖3(b)所示的模板進行模板匹配，結(jié)果如圖4所示，其中模板和目標相關(guān)性最大處為最佳匹配位置(圖中最亮處)。圖3(a)中存在與前兩個模板形狀相似的目標，圖4(a)、圖4(b)中最亮部分是這些目標的匹配點位置。圖3(a)中沒有第3個模板圖形狀的目標，圖4(c)中最亮點不明確，該模板不能檢測出圖中的建筑物。第4個模板和圖3(a)中左下角建筑物形狀相似，圖4(d)下方有一個最亮點，目標得到正確的識別。

圖4 模板與目標匹配結(jié)果

    將Markov隨機場模型用于遙感影像建筑物目標分割，可以實現(xiàn)單幅影像中多目標建筑物分割。分割結(jié)果經(jīng)形態(tài)學處理后，與常見建筑物形狀模板進行相關(guān)匹配，能夠較為準確地檢測出目標建筑物位置。在目標分割階段，MRF通過建立鄰域系統(tǒng),考慮到相鄰像素間的相互作用,與單純的K-均值聚類相比，分割效果有很大的提高。但由于該方法只利用了像素的顏色值信息，因此對顏色值比較依賴。建筑物受到陰影遮擋時目標分割不完整的‘欠分割’，或是遇到反射性質(zhì)與其相近區(qū)域造成分割被擴大的‘過分割’狀況，都是這一原因造成的。在目標檢測階段，利用形狀模板和目標建筑物的相關(guān)性進行匹配，能把圖像區(qū)域和模板形狀一致的目標識別出來。但只有在模板和目標相似度較高的情況下才能得到準確的匹配位置，需要進一步改進以降低幾何失真對匹配性能的影響。
參考文獻
[1] 謝興燦，裴繼紅，黃建軍，等.城市航空影像中建筑群的FCM聚類分割方法[J].計算機仿真，2004,21(4)：43-46.
[2] 王世偉，方濤. 基于對象建模的遙感影像建筑物提取方 法[J].計算機仿真，2010,27(12)：254-257.
[3] 田昊，楊劍，汪彥明，等.基于先驗形狀約束水平集模型的建筑物提取方法[J].自動化學報，2010,36（11）：1502-1511.
[4] 李志軍，牛照東，陳曾平.基于模板的前視紅外圖像建筑物分割方法[J].光電工程，2012,39(5)：109-114.
[5] YANG M, WEI Z Q, CHU J, et al. A framework for automatic recovery of camera rotation and parallelepiped contour extraction in simple urban scene based on morphological sieves[J]. IEEE Transactions on Image Processing,2008,10(9):1332-1345.
[6] TRINH H H, KIM D N, JO K H. Urban building detection by visual and geometrical features[C]. International Conference on Control, Automation and Systems, 2007：1779-1784.