基于改進變分水平集方法的遙感影像道路提取

胡秀云

（信陽師范學院華銳學院 社會科學系，河南 信陽 464000）

0 引言

長期以來，遙感影像道路提取被認為是遙感、計算機視覺、GIS等領域的研究熱點之一。在高分辨率遙感影像中，道路表現為具有一定長度和寬度，拓撲結構復雜的一類地物[1]，且受非目標地物（車輛、行人、樹木陰影等）影響較大。因此，在高分辨率遙感影像上進行道路提取一直被認為是一項很具難度的工作。

活動輪廓模型是一種自上而下的能夠有效融合上層知識和底層圖像特征的處理過程，特別適用于建模以及提取任意形狀的變形輪廓等。已有許多學者將活動輪廓模型應用于道路提取，并取得了一定的研究成果[1-2]，但是這些已有的方法有的僅適用于鄉村等背景簡單的場景[2]；有的人工干預較多，且時間效率較低[1]。

水平集表達模型作為活動輪廓模型研究的一個重要分支，當提取對象是圖像中的所有目標的邊界時，水平集表達模型是一種很有效的方法[3-7]。但是，如果僅需要提取感興趣的目標，特別是對于背景比較復雜的圖像，則利用該方法進行目標提取就變得很具難度。本文以城區高分辨遙感影像主干道路為研究對象，在現有研究成果的基礎上，提出一種適用于提取復雜背景下感興趣目標的變分水平集方法，并在第3節，利用多幅QuickBird衛星高分辨率遙感影像對該模型的有效性進行了驗證。實驗結果表明該模型能夠較精確地從復雜背景中提取出道路，具有較強的抗干擾能力。

1GACV模型

C-V模型[4]對弱邊界和離散邊界一定提取優勢，為進一步提高其在復雜背景下的分割能力，以及對局部邊界的定位精度，Chen等在文獻[6]中按照一種合理的方式將GAC（Geodesic Active Contour）模型[3]與 C-V（Chan-Vese）模型[4]相結合，提出一種 GACV(Geodesic-Aided C-V)模型。該模型關于彩色圖像的水平集演化方程為：

該算法綜合了GAC模型和C-V模型的優點，但仍存在一些不足，如：對于一些較弱的邊界，仍然存在邊界“泄露”現象；對于諸如高分辨率遙感影像的背景相當復雜的圖像，不能有效地分割出用戶感興趣的目標，提取結果不理想。

2 一種改進的彩色GACV模型

本節通過構造目標識別函數，根據Beltrami框架重新構造彩色梯度流，以及將文獻[7]中提出的懲罰能量函數引入GACV模型，提出一種改進的彩色GACV模型。

2.1 新約束條件構造

文獻[8]中，提出一種基于PCA的彩色區域生長算法。該方法首先分析選取的道路種子點R,G,B特征的主成分，并把整個圖像映射到主要特征分量上；然后根據種子點在主分量圖中的信息，利用區間估計方法確定道路特征的置信區間，此處記為[a，b]；最后，依次以這些道路種子點為生長起始點，以區間[a，b]為生長準則，進行區域生長，從而獲得如下初始道路圖 β（x，y）：

在GACV模型的邊緣檢測子中引入β（x，y）作為目標識別函數，并利用圖像顏色梯度信息重新構造如下邊緣檢測子：

式中Gσ表示方差為σ的Gauss函數；Λ是一種基于灰度圖像的梯度幅度對應于彩色圖像的推廣形式，根據Beltrami框架[5]，其可以按照如下過程進行計算：

把彩色圖像看作五維譜空間（x，y，R，G，B）中的二維平面（x，y，R（x，y），G（x，y），B（x，y）），則 Λ 為該譜空間的度量張量 gij的特征值的最大值，gij如式（4）所示：

2.2 改進的GACV彩色模型

其中，μ≥0，v≥0；λoi＞0，λbi＞0， 它們對應于每個顏色通道；＝），＝）是有待于在曲線演化過程中確定的向量；H（φ），δ（φ）分別為的規整化的 Heaviside函數和 Delta函數[4]。

與文獻[4]中C-V模型的求解過程類似，通過最小化能量函數Fε），可得到水平集體演化方程：

公式（10）中的第一項為懲罰項，第二項為“區域檢測子”，第三項為吸引項，同時引入目標識別函數β（x，y）。

3 遙感影像道路提取算法

本節我們將第2節提出的變分水平集方法應用于高分辨率遙感影像道路提取，并通過仿真實驗，在Matlab7.6環境下驗證其有效性。

3.1 算法描述

該道路提取算法的具體實現過程為：

首先，利用基于PCA的彩色區域生長算法對道路進行預分割，獲取大致的道路區域，輸出結果為一個二值圖像（1表示道路，0表示背景），我們定義該二值圖像為目標識別函數β（x，y）。

然后，按照2.1節的方法，根據Beltrami框架，利用目標識別函數β（x，y）和RGB顏色空間中各通道的顏色信息，構造彩色梯度流gcolor，并利用這些約束項，按照2.2節所提出的方法，構造水平集模型。

最后，利用預分割所獲得的二值道路圖初始化水平集函數φ0，進一步利用該改進的變分水平集方法進行道路邊界演化，以更精確的定位道路。

3.2 仿真實驗

該道路提取算法的仿真實驗結果如圖1所示。為了更清楚的看出改進模型的有效性，將其與基于GACV模型的道路提取結果（圖2）進行了對比。兩種方法的時間性能對比見表1。在實驗中，改進變分水平集方法的參數設置為：Δt＝1，η＝0.2 ／Δt，μ＝0.5，λoi＝λbi＝1，i＝1，2，3，v＝0，GACV 模型的參數設置為：μ＝1，v＝0，λoi＝λbi＝1，i＝1，2，3，τ＝1，Δt＝0.1。

圖1 基于改進變分水平集模型的道路提取效果

圖2 基于GACV模型的道路提取效果

表1 改進變分水平集方法與GACV模型性能對比

圖1（a）可見利用彩色區域生長算法能夠較有效地獲取大致的道路區域,只是局部邊界存在“過生長”和“鋸齒”現象。圖1（b）可見利用改進變分水平集方法能夠較精確、完整地定位道路邊界，具有較強的抗干擾能力，只有在道路附近噪聲較大時，存在局部定位不太準確的現象。此外，由表1可知，該算法經過較少次的迭代即能收斂到真實的道路邊界，時間效率較高。

圖2為GACV模型進行道路提取的結果，可見很多與道路光譜特性相近的地物也被分割出來，因此結果比較雜亂，且由表1可知，由于初始輪廓位置離真實道路邊界較遠，且在演化過程中需要不斷地對水平集函數進行分段初始化，基于該模型的道路提取算法時間效率較低。

4 結束語

本文在現有水平集方法的基礎上，提出一種應用于高分辨率遙感影像道路提取的變分水平集方法，從復雜背景中提取感興趣目標。在該模型的設計過程中，不僅考慮了圖像的邊緣、區域信息，還通過構造目標識別函數將目標的顏色信息引入到模型中，同時通過引入一個懲罰能量項完全消除了水平集演化過程中的重新初始化操作。通過驗證該模型能夠較精確地提取高分辨率遙感影像中的道路，時間效率較高，具有較強的抗干擾能力和實用性。

［1］馬振，吳紀桃，羅中華.基于水平集方法的遙感影像道路提取算法[C]//第十三屆全國圖象圖形學學術會議.2007:433-437.

［2］I.Laptev,H.Mayer,T.Lindeberg,W.Eckstein,C.Stergerand A.Baumgartner.Automatic extraction of roads from aerial images based on scale space and snakes[J].Machine Vision and Application,2000,12(1):23-31.

［3］CassellesV,KimmelR,Sap iro G.Geodesic active contours[C]//In:Proceedings of International Conference on ComputerVision.Boston,MA,USA,1995:694-699.

［4］Chan T,Vese L.Active contours without edges[J].IEEE Transactions on Image Processing,2001,10(2):266-277.

［5］R.Goldenberg,R.Kimmel,E.Rivlin,and M.Rudzsky.Fast geodesic active contours[J].IEEE Trans.Image Process,2001,10:1467-1475.

［6］L.Chen,Y.Zhou,Y.G.Wang,J.Yang.GACV:geodesic-aided C-V method[J].Pattern Recognition,2006,39(7):1391－1395.

［7］LI Chun-Ming,XU Chen-Yang and GUI Chang-Feng,et al.Level set evolution without re-initialization:a new variational formulation[C]//IEEE International Conference on Computer Vision and Pattern Recognition(CVPR).San Diego,2005,1:430-436.

［8］S.Guillermo,Color snakes,Comput.Vision Image Under standing[J].1997,68:247－253.