傳統分水嶺變換具有簡單、速度快、能檢測出弱邊緣物件及能獲得物件完整邊界等諸多優點,被廣泛的用於各種領域。但因其一般是在梯度影像進行變換,受暗雜訊等因素的影響,存在大量偽區域性極小值區域,會出現過分割現象[1]。高解析度遙感影像在大氣傳遞、成像過程容易產生雜訊,特別是地物內部紋理豐富,利用傳統分水嶺方法分割時更容易出現過分割現象。
同時,利用傳統分水嶺演算法分割高解析度影像區域時,在強邊界附近易於造成交叉,使分割效果變差[2,3]。
基於上述問題,本文利用方向分水嶺變換演算法增強弧和初始方向邊緣強度的訊號之間的一致性,反映邊緣的強弱,並將其作為區域合併條件,完成影像多尺度分割。
1 方向分水嶺變換原理
方向分水嶺變換演算法的區域提取方法的主要思想:由傳統分水嶺變換對上述得到的邊緣影像進行區域提取,通過邊緣方向的擬合併修正邊緣的強弱,獲取區域分割結果。作為方向分水嶺變換演算法資料的輸入,要先獲得影像的邊緣的方向和強度資訊——邊緣影像。
為進一步提高邊緣的定位精度,將影像的光譜資訊特徵和紋理特徵相結合,提取不同尺度下影像的邊緣。因此,根據影像的亮度和紋理特徵,利用直方圖差分方法分別得到影像多尺度的分割結果,通過對它們在不同方向下設定權重,獲得不同尺度下影像的方向邊緣強度[4]。
(1)首先,在三個尺度下,根據灰度值和紋理基元特徵分別計算八個方向的直方圖的方向梯度,檢測出灰度影像邊緣和紋理區域邊緣。
(2)找出八個方向的響應極大值,作為邊界。
(3)將上述得到的多特徵、多尺度得到的影像,設定權重以下公式形式組合起來,得到邊緣強度mpb。
2 方向分水嶺變換影像多尺度分割演算法
經過方向分水嶺變換分割演算法能夠得到影像的初始分割結果,但是通常是過分割的結果,本文利用區域邊緣的強度,逐步合併弱邊緣的區域,獲得影像的層次結構。在合併過程中,由於需要知道區域和邊界的拓撲關係,因此多層次合併的方向分水嶺演算法分為以下兩步。
第一步,利用分水嶺方法得到過分割的影像,作為最細粒度的分割影像,並獲取像元、邊緣、區域和區域邊界的相關資訊。該步驟分為以下幾個方面內容。
(1)計算影像所有畫素點八個方向下的全域性梯度的最大值。根據邊界響應的極大值作為邊緣,得到過分割影象。
(2)對上一步結果進行分水嶺變換。根據的極小值得到種子位置,為分水嶺進行同質區分割做準備,接著,運用分水嶺演算法得到分割區域。
(3)提取邊界和鄰接區域資訊。內容包括:提取邊界兩側的區域、建立頂點-邊緣對映圖、提取頂點座標、提取邊緣的端點、提取邊緣上各個像元的座標、提取封閉邊緣的標識以及邊界擬合。
因此,得到相應折線段上的像元。
第二步,將區域之間邊界的pb均值作為區域的不相似性。
(1)逐弧段計算權重,即鄰接區域間的合併條件。
(2)搜尋權重最小的弧段,合併其鄰接區域,並將弧段融入合併後的區域中。
(3)區域性調整拓撲鄰接關係,即調集成並後區域的邊界弧段資訊以及邊界弧段的臨界區域資訊。
(4)區域性調整弧段權重,修改合併後區域邊界弧段的權重。
(5)判斷是否滿足合併終止條件,若不滿足,則返回(2);若滿足,則合併終止。通過建立拓撲鄰接關係,可快速查詢弧段與區域、區域與區域之間空間關係,及時區域性調整以應對合併產生的變化。
3 實驗與分析
為驗證本文提出方法的有效性,選用航空影像進行分割實驗。圖1(a)中航空影像的大小為509×368,影像上主要包含居民地、道路、樹木、水體及耕地等。圖1(b～d)是圖1(a)的層次分割結果,圖1(b)是最底層的分割,獲得區域數最多,得到影像的過分割結果。
k是區域之間的邊界強度,k的取值範圍為[0,255]。k=0的情況對應圖1(b),當k=255時,區域個數為1,即整幅影像。
圖1中,(b)上明顯存在著過分割現象。隨著區域的合併,某些區域的過分割現象逐漸消失。如k=50時耕地被分割為完整一塊,finest、k=30上裸地也被分為完整的一塊,finest、k=30上房屋過分割現象得以消除。
但是,隨著區域合併的增加,有些地方也出現了欠分割的現象。如k=50上,右下角的地塊與相鄰的小路被劃分為一塊。隨著k的取值由小變大,較小區域之間逐漸合併,過分割現象慢慢減少;當分割尺度達到某一值時,某些地方會出現欠分割現象。
上述實驗結果,隨著區域合併的增加,分割結果中斑塊的尺寸也在增大,而且大尺度上斑塊是小尺度上斑塊的集合。可知,區域合併初期,最底層影像的區域物件內部相似性較強。但是,隨著區域合併,區域變大,區域之間的異質性增強。
因此,基於方向分水嶺變換方法可以獲得同質性最強的分割結果,這為以區域物件為節點的尺度綜合提供了較理想的分割結果。
4 結語
本文提出了一種方向分水嶺高解析度遙感影像分割方法,先由傳統分水嶺變換對上述得到的邊緣影像進行區域提取,通過邊緣方向的擬合併修正邊緣的強弱,獲取區域分割結果。接著,提出對方向分水嶺高解析度遙感影像分割結果進行區域之間的合併,根據區域邊界的強弱作為區域合併條件,從而獲得多尺度的分割結果。實驗表明不同層次的分割結果能夠較好地反映地物的尺度特性。
參考文獻
[1]章毓晉.影象分割[m].北京:科學出版社,2001:10～17.
[2]梅天燦,李德仁,秦前清.分水嶺變換在遙感影像現狀特徵提取中的應用[j].武漢大學學報(資訊科學版),2004,29(4):338～341.
[3]宮鵬,黎夏,徐冰.高解析度影像解譯理論與應用方法中的一些研究問題[j].遙感學報,2006,10(1):1～5.
[4]程杰.一種基於直方圖的分割方法[j],華中理工大學學報,1999,27(1):84～86.
影象處理分水嶺演算法實驗報告
一 實驗目的 1.通過matlab實驗的具體操作來掌握空間影象處理演算法 2.通過實驗來驗證所學知識，達到學以致用 3.通過實驗來了解分水嶺演算法的優缺點。二 實驗原理 分水嶺分割方法，是一種基於拓撲理論的數學形態學的分割方法，其基本思想是把影象看作是測地學上的拓撲地貌，影象中每一點畫素的灰度值表示...
自然水系 或分水嶺 自動提取的matlab實現方案
水系提取演算法 d8 的matlab實現 d8演算法是當今非常成熟的提取水系 或分水嶺 的計算機實現方法。這裡筆者結合自身的程式設計實踐介紹d8演算法的matlab實現過程，謹與大家分享。d8演算法分為兩部分實現 每個柵格點水流方向的計算 每個柵格點匯水面積的計算。1 水流方向的計算 在水流方向的計...
分水嶺完小紅領巾心向黨活動總結
為進一步增強廣大少年兒童對黨和社會主義祖國的樸素感情，我校結合實際，開展了豐富多彩的活動，深入推進了 紅領巾心向黨 主題教育活動的開展，為迎接黨的十八大勝利召開,表達 紅領巾心向黨 的美好心願。並取得較好的效果，現將活動總結如下 一 舉行 唱紅歌 誦紅詩 贊祖國 為主題的活動 通過此次活動，加強同學...