基于藍噪聲理論的遙感圖像森林植被紋理測量

劉小丹， 楊燊

(遼寧師范大學計算機與信息技術學院,大連 116029)

基于藍噪聲理論的遙感圖像森林植被紋理測量

劉小丹， 楊燊

(遼寧師范大學計算機與信息技術學院,大連 116029)

遙感圖像分割中森林植被是重要的一類目標，有效確定森林植被的紋理尺度是紋理分割的重要問題。提出一種用藍噪聲理論描述遙感圖像森林植被紋理特征的方法，是一種新的植被紋理刻畫和紋理尺度計算方法。研究尺度與植被紋理形態的對應關系，對于選定的探測區域，迭代尋找藍噪聲特征。迭代過程包含通過幾何變換縮小區域的尺寸，用快速傅里葉變換獲取區域的頻譜響應，從頻譜響應中提取藍噪聲特征。對于具有藍噪聲特征的區域，計算森林植被紋理的灰度分布，根據當前區域尺寸計算紋理的尺寸。實驗表明，森林植被紋理單元的尺度和灰度分布測量結果準確，為進一步紋理分割提供了可靠的基礎。

遙感圖像；紋理分割；藍噪聲；森林植被

0 引言

森林植被是區域生態環境的標志之一，也是面積較大的地面目標。遙感圖像森林植被分割在生態系統分析、國土資源管理及國防建設等領域均有重要應用。利用統計分析、頻域分析和紋理結構模式分析等方法分割遙感圖像中的森林植被區域是目前主要研究方法，其中，紋理特征的提取和分析是基礎問題[1-2]。紋理尺度和其他屬性(如形狀和灰度分布等)是分割的重要條件。在同一幅圖像中不同地物具有不同的最優空間尺度[3]。對于森林植被的紋理分割來說，森林植被紋理空間結構性較強，準確測量紋理尺度和灰度分布等特征非常重要，它決定分割后的森林植被區域有恰當的邊界，既不過分割，也不欠分割。在森林植被紋理結構模式分析和數學形態學濾波等方法中，紋理尺度和灰度分布參數已經可以用于實現自適應的紋理結構模板和形態學結構元，從而得到準確的區域分割結果[4]。雖然目前已有多種組合的紋理特征提取方法的紋理模型，但普遍存在紋理特征表達確切度不高的問題，尋找有效的森林植被紋理特征測量手段，是實現森林植被分割的關鍵[5-6]。

本文根據藍噪聲理論，提出一種探測森林植被紋理尺度和灰度分布的有效方法。除冬季的落葉喬木外，遙感圖像森林植被的紋理在特定的尺度上呈高頻隨機信號狀態，具有藍噪聲的基本特性。使用森林植被基礎特征快速選擇圖像的探測區域，通過快速傅里葉變換探測區域在不同尺度下的藍噪聲特征，利用探測到的藍噪聲特征及對應尺度計算森林植被的紋理尺度和灰度分布。這些紋理特征可以為后續的統計分析、結構模式分析及形態學濾波等基于紋理的分割方法提供可靠的基礎數據。

1 森林植被紋理的信號特征

1.1 紋理信號分析

一般而言，紋理形態的描述分宏觀和微觀2個方面，宏觀屬性包含周期、方向和單元尺度等，微觀屬性包含單元幾何形狀和對比度等。在高分辨率的遙感圖像中識別森林植被紋理，首先應找到恰當的空間尺度，這種特定的尺度與森林植被的宏觀屬性相關[7]。在遙感圖像中，自然狀態生長的森林植被在色調、植株位置、樹冠尺寸、高光及陰影對比度等方面的特性構成區別于其他地物的主要特征，其中，植株位置包括位置的非周期隨機性和空間密集性。宏觀上，森林植被由單個的樹冠連接構成，呈非周期性隨機狀態，在陽光照射下有較弱的方向性；微觀上，紋理單元呈簇團狀結構，中心與周圍對比度明顯，結構形狀隨植物種類有一定變化。當遙感圖像縮小到特定尺度時，樹冠的高光區域和陰影區域成為單個像素，高頻性質明顯。非周期、高頻和隨機的森林植被紋理特點恰好具有藍噪聲的基本特性。

藍噪聲指具有最小的低頻分量并且頻譜中沒有明顯峰值出現的非周期性隨機信號，在有限帶寬內，其功率譜密度與頻率成正比[8]。在計算機圖形學中，藍噪聲用于生成半色調抖動圖像，產生適合印刷等再現方式的像素空間結構,這種結構更符合人類視覺特點[9]。遙感圖像森林植被紋理的藍噪聲特征正是人眼在遙感圖像中判斷森林植被區域的基礎依據。由此，可以形成森林植被紋理識別的新方法，其中紋理尺度探測的有效性是關鍵。

1.2 藍噪聲特征探測

遙感圖像森林植被紋理單元空間分布規律具有藍噪聲的特點，即樹冠紋理單元以單個像素點為尺度時所形成的宏觀信號特性。判斷遙感圖像區域在某個尺度上具有藍噪聲特征，是探測的核心問題。森林植被紋理的空間形態受植物種類和拍攝條件等因素的影響，在紋理單元尺度、幾何形狀、對比度及分布規律等方面有較大差異，但在同一幅遙感圖像中，這些差異并不明顯。通過圖像尺度的動態變換來尋找紋理宏觀的藍噪聲特征。對每個尺度下的區域圖像進行快速傅里葉變換，分析信號頻域特征，完成判斷。

首先，取遙感圖像中典型森林植被區域進行分析，全色圖像尺寸為512像素×512像素，全色分辨率為0.65 m，多光譜分辨率為2.76 m。圖1是區域對應的多光譜圖像、增強處理圖像及其頻譜響應圖。增強處理主要是自適應的對比度拉伸，目的是增大信號的振幅，便于提取頻率信息。樹冠作為紋理單元，尺寸較大，傅里葉變換后頻譜中低頻成分較多，不具備藍噪聲特征。當該圖像尺寸縮小到64像素×64像素時(圖1(b))，樹冠紋理單元的尺寸接近2，即高光部分和陰影部分各對應1個像素，頻譜中低頻成分明顯減少，高頻成分劇增，表現出典型的藍噪聲特征。

(a) 512像素×512像素區域(左)、增強處理圖像(中)及其頻譜響應(右)

(b) 縮小到64像素×64像素區域(左)、增強處理圖像(中)及其頻譜響應(右)

對于森林與其他地物并存的區域，如果森林占區域面積的比例不夠大，就稱為非典型區域。這類區域在尺度縮小的過程中不會出現典型的藍噪聲現象。在圖2中，3個區域都經過尺寸縮小和對比度增強處理。城市區域的頻譜響應中低頻信號很強，有明顯的方向性，這是因為城市道路和房屋是有方向性的。有波浪的水面區域的紋理對比度較小，對比度拉伸處理可導致大面積平滑區域出現，從而產生較強的低頻信號，并且由于水面的波浪紋理有一致的方向，對應的高頻信號較強但有明顯的方向性。草場區域的紋理對比度較小，且缺乏細節，在縮小和對比度拉伸處理后基本上成為平滑區域，頻譜響應中低頻信號很強。存在非森林地物具有藍噪聲信號特征的情況，如亂石灘等，可以通過色調等其他特征來區別。對于非自然生長的人造林區域，由于植株位置較規律，頻譜響應中會在相關方向上出現高頻區域的高峰，在低頻區域仍然具有較低的信號能量。當出現顯著的藍噪聲特征時，由當前區域尺寸與原區域尺寸之比，可以方便地計算出樹冠紋理的尺度。

(a) 縮小到64像素×64像素城市區域及其頻譜響應

(b) 縮小到64像素×64像素有波浪水面區域及其頻譜響應

(c) 縮小到64像素×64像素草場區域及其頻譜響應

2 紋理測量的區域選擇

探測藍噪聲特征，需要選擇合適的區域，以提高探測效率。將遙感圖像N等分，根據各區域的基礎特征快速篩選最合適的區域。本文選用區域的色調和灰度分布2項特征來實現典型森林區域判斷。

2.1 區域的平均色調

遙感圖像中大部分的森林植被以綠色調為主。冬季對于落葉喬木來說可能會有其他色調，而且樹冠的形態也會有很大變化，這有待其他方法解決[10]。由于HSV色彩空間的顏色與亮度無關，本文將圖像區域從RGB空間變換到HSV空間，然后對整個區域統計平均H值。根據HSV色彩模型，綠色調的H值以120為最優[11]。對區域的平均H值進行歸一化處理，作為基礎特征之一。

2.2 區域灰度正態分布

遙感圖像中典型森林植被呈成片分布。樹冠作為紋理單元微觀上呈簇團狀結構，中心與周圍對比度明顯，灰度過渡平滑，中間灰度像素居多。在典型森林區域中，樹冠紋理單元覆蓋整個區域，紋理單元的這種灰度分布重復出現，并且集中在低灰度范圍內。對典型森林區域統計灰度概率密度(即像元數)，其直方圖呈正態分布且灰度均值處于低灰度段，見圖3(a)。這種特性可作為選擇合適測量區域的基礎特征之一。非森林地物的紋理在灰度分布規律上明顯不同于森林紋理，混合區域的直方圖是多種灰度分布的合成，故容易與典型森林區域直方圖區分，如圖3(b)(c)所示。

(a) 森林區域及其直方圖

(b) 森林與城區混合區域及其直方圖

(c) 森林與水面混合區域及其直方圖

在圖3中城市道路區域在快速傅里葉變換后的頻譜圖上具有明顯的方向性，這是因為城市道路和房屋的建設是有規律性的，人為建筑具有一定的方向性。所以在頻譜中也能很明顯發現其高頻區域(即道路)也具有方向性。

正態分布的密度函數是關于均值對稱，在均值處達到最大值的一條位于X軸上方的鐘形曲線。Jarque-Bera檢驗是用來判斷一組樣本是否來自正態總體的方法[12]。設區域圖像灰度數據為向量X，數據容量為n，X的偏度(三階矩)為S，峰度(四階矩)為K，則Jarque-Bera檢驗(JB)公式為

(1)

在正態分布情況下，X的偏度S=0，峰度K=3，JB=0。如果JB較大，則不能認為樣本符合正態分布。

將JB做歸一化處理，再根據平均色調和正態分布2個基礎特征對區域進行選擇，為此作為藍噪聲特征檢測的區域圖像。

3 紋理特征測量

3.1 算法流程

N等分遙感圖像，計算各區域的色調和灰度分布，篩選最合適的區域，并對該區域進行自適應對比度增強處理。對于選定的探測區域，迭代尋找藍噪聲特征。迭代過程包含通過幾何變換縮小區域的尺寸，用快速傅里葉變換獲取區域的頻譜響應，從頻譜響應中提取藍噪聲特征。對于具有藍噪聲特征的區域，計算森林植被紋理的灰度分布，根據當前區域尺寸計算紋理尺寸。紋理特征測量的算法流程見圖4。

圖4 紋理特征測量的算法流程Fig.4 Texture measurement algorithm flow

在區域尺寸較大的階段，森林植被紋理的微觀結構特性明顯，紋理單元的尺寸較大且灰度過渡平滑，低頻信號較強，呈紅噪聲或粉紅噪聲特征。頻譜響應會隨著空間尺寸縮變的每一步發生變化，當森林植被紋理尺寸縮小到接近像素級時，由微觀結構形成的低頻信號消失，而突出高頻振蕩的隨機特性。此時能量譜密度函數呈藍噪聲特征，即低頻段接近零值，高頻段能量譜密度與頻率成正比。如果在達到最小尺寸之前，區域沒有這種藍噪聲特性出現，則認為本區域不是典型森林植被紋理區域。

3.2 區域藍噪聲特征提取

預處理包括對快速傅里葉變換后的區域頻譜進行移頻和平滑。

為判斷方向性，首先計算以原點為中心等角度間隔的N個方向的一維能量譜Ei1，Ei2，…，EiK′，再沿方向i(i=1,2,…,N)求能量之和Pi，即

(2)

式中：K′為單一方向上的能量譜個數之和；j為臨時變量。最后計算各方向能量方差σ,即

(3)

式中：pi為方向上隨機一點的能量值;pμ為pi的均值。若σ=0，則無方向性。

如方向性不強，將N個方向的一維能量譜合并為全方位均值一維能量譜E1,E2,…，EK′，其能量之和為PT。利用

(4)

衡量全方位一維能量譜的偏度S。式中:fx為頻率；fm為頻率的中心值；Ex為N個方向一維頻譜值之和。一維情況下，藍噪聲的能量譜密度分布偏右是主要特征之一，S>0時為右偏。藍噪聲的另一個主要特征是能量譜密度與頻率成正比，理想情況下，函數k(x)>0,且為常數，k(x)定義為

(5)

對k(x)計算方差，以方差接近零的程度評價符合能量譜密度與頻率成正比的程度。

方向能量方差σ接近0，能量譜密度分布偏度S>0，且能量譜密度與頻率之比k(x)的方差接近0的區域，可以判定為具有藍噪聲特征。

3.3 紋理尺度計算

當探測到藍噪聲特征時，根據當前區域尺寸和原區域尺寸，計算樹冠紋理的尺度，即樹冠直徑的像素個數。由于紋理單元的高光區域和陰影區域分別對應一個像素點，故原區域尺寸與當前區域尺寸之比再乘2就是紋理的尺度。當然，這是原區域中森林植被紋理的平均尺度。

3.4 紋理灰度分布計算

紋理的灰度分布是基于形態學濾波分割方法設計結構元的依據。當探測到藍噪聲特征時，區域圖像內容應以森林植被為主。樹冠紋理重復出現，概率密度函數在低灰度區呈正態分布。此時可以用原區域的灰度分布直方圖來計算樹冠紋理的平均灰度分布狀況，包括平均灰度、灰度對比度及峰度。灰度對比度C計算公式為

(6)

式中：相鄰像素間灰度差δ(i,j)=|i-j|；Pδ(i,j)為當δ=G時的像素分布概率；M為區域中像素相鄰關系的總個數。

4 實驗結果與分析

本文采用Matlab7.12.0為實驗環境。所用遙感圖像為QuickBird衛星(全色和多光譜)數據產品，全色分辨率為0.65 m，多光譜分辨率為2.76 m。本文算法對遙感圖像全色分辨率的要求是高于1.5 m，以保證樹冠紋理單元有足夠的像素用于尺寸縮小處理以及提取灰度分布處理，而對多光譜分辨率則無特殊要求。

實驗選取包含森林植被的遙感圖像6幅，如表1所示。圖像尺寸為2 048像素×2 048像素，但不同圖像中森林植被紋理單元具有不同的尺度。

表1 森林植被紋理測量結果Tab.1 Forest vegetation texture measurements

設N=16，即對原圖像做16等分用于選擇經色調運算與灰度正態分布統計后的最佳探測區域，故測量算法開始時原區域尺寸為512像素×512像素，每次水平和垂直方向各縮小25%。區域尺寸縮小處理中對最小尺寸設定為50像素×50像素。首先考慮遙感圖像中樹冠的尺度與圖像尺度的比例，其次考慮快速傅里葉變換精度對圖像尺度的要求。遙感圖像的拍攝條件以及地形等因素會嚴重影響成像的清晰度和圖像整體的灰度分布，并且同一片森林中經常存在不同樹種，樹冠的尺寸并不是很一致。這些情況反映在信號的藍噪聲特征上，往往是藍噪聲探測的3個主要參數不很典型，所以，本文在閾值選擇上做了適當放寬。

森林植被紋理灰度分布的3項數據計算都是常規的統計處理，其數據準確度主要依賴于參與計算的圖像區域是否典型。由藍噪聲探測得到的區域是整個圖像中相對典型的森林植被區域，所以這種統計的效果應該是最佳的。由表1可以看出，6幅圖像中探測區域的選擇是較準確的。人工測量的紋理尺度是探測區域中樹冠紋理抽樣的人工判讀平均尺度。對每個512像素×512像素的探測區域，使用圖像處理工具Photoshop沿中部水平方向和中部垂直方向抽取2個條帶，人工判讀樹冠并累計個數，計算出平均樹冠直徑的像素數。本文算法獲得的紋理尺度結果與人工測量的結果偏差不大。沒有選擇冬季的落葉喬木林，因為本文方法還不能完成對這種目標的識別。人工林的處理在本文方法基礎上適當調整方向能量方差判斷部分，應該可以完成，不過本文沒有做相關的實驗。對于全色分辨率<1.5 m的遙感圖像，本算法無效。

5 結論

1)紋理的特征提取是紋理圖像分割最重要的環節之一，關系到分割的準確度。遙感圖像中的森林植被受季節和拍攝條件的影響，在顏色和灰度等方面一致性較差，但紋理信息卻穩定。以藍噪聲特征探測為手段，在遙感圖像典型區域中測量森林植被紋理的尺度及灰度分布信息，是本文提出的一種新的森林植被紋理測量方法，也是一種新的紋理描述手段，為進一步紋理分割提供了可靠的基礎。

2)本文測量結果數據可以為進一步的森林植被紋理分割提供支持。例如，采用灰度形態學紋理濾波方法進行分割時，就會涉及與紋理尺度相關的結構元尺度、與紋理灰度分布相關的結構元灰度分布等。結構元設計必須與待處理圖像中分割目標的紋理相適應，這是分割效果的決定性因素。本文證明了遙感圖像植被紋理在特定空間尺度上具有藍噪聲特征，這對處理類似遙感圖像地物紋理是一種方法上的參考。

[1] Chubey M S,Franklin S E,Wulder M A.Object-based analysis of IKONOS-2 imagery for extraction of forest inventory parameters[J].Photogrammetry and Remote Sensing,2006,72(4):383-394.

[2] 劉雯,駱劍承,沈占鋒,等.多尺度植被信息提取模型研究[J].計算機應用研究,2009,26(6):2398-2400. Liu W,Luo J C,Shen Z F,et al.Research of multi-scale iterative vegetation extraction model[J].Application Research of Computers,2009,26(6):2398-2400.

[3] 張學良,肖鵬峰,馮學智.基于圖像內容層次表征的遙感圖像分割方法[J].中國圖象圖形學報,2012,17(1):142-149. Zhang X L,Xiao P F,Feng X Z.Remotely sensed imagery segmentation based on the hierarchical representation of image content[J].Journal of Image and Graphics,2012,17(1):142-149.

[4] 李杰,苗長云,武志剛,等.基于數學形態學的圖像邊緣檢測算法的研究[J].計算機科學,2012,39(6A):546-548. Li J,Miao C Y,Wu Z G,et al.Algorithm for edge detection of image based on mathematical morphology[J].Computer Science,2012,39(6A):546-548.

[5] Rodriguez-Galiano V F,Chica-Olmo M,Abarca-Hernandez F,et al.Random forest classification of Mediterranean land cover using multi-seasonal imagery and multi-seasonal texture[J].Remote Sensing of Environment,2012,121:93-107.

[6] Hu M X.Extraction of trees and grasses in urban green land based on texture character of high resolution remote sensing image[C]//2011 International Conference on Remote Sensing,Environment and Transportation Engineering.Nanjing:IEEE,2011:2178-2181.

[7] Lian L,Chen J F.Research on segmentation scale of multi-resources remote sensing data based on object-oriented[J].Procedia Earth and Planetary Science,2011,2:352-357.

[8] Ulichney R A.Dithering with blue noise[J].Proceedings of the IEEE,1988,76(1):56-79.

[9] Rodríguez J B,Arce G R,Lau D L.Blue-noise multitone dithering[J].IEEE Transactions on Image Processing,2008,17(8):1368-1382.

[10]李成,陳仁喜,王秋燕.改進的基于視覺認知特征的植被識別方法[J].國土資源遙感,2013,25(2):75-80.doi:10.6046/gtzyyg.2013.02.14. Li C,Chen R X,Wang Q Y.Improvement of the automatic recognition method based on vegetation visual characteristics[J].Remote Sensing for Land and Resources,2013,25(2):75-80.doi:10.6046/gtzyyg.2013.02.14.

[11]陳曉蘭.遙感圖像的色調統一研究[D].杭州:浙江大學,2007. Chen X L.Study on Color Harmonization of Remote Sensing Images[D].Hangzhou:Zhejiang University,2007.

[12]Jarque C M,Bera A K.A test for normality of observations and regression residuals[J].International Statistical Review,1987,55(2):163-172.

(責任編輯：邢宇)

Forest vegetation texture measurement of remote sensing images based on the blue noise theory

LIU Xiaodan,YANG Shen

(CollegeofComputerandInformationTechnology,LiaoningNormalUniversity,Dalian116029,China)

Forest vegetation remote sensing image segmentation is an important kind of target,and effective determination of the scale of forest vegetation texture segmentation is an important research topic. This paper presents a method in which the blue noise theory is used to describe the characteristics of remote sensing images for forest vegetation texture. This is a new method for vegetation texture characterization and texture scale calculation. The correspondence between the research scale morphology and vegetation textures can be used in the selected detection area to iteratively search for blue noise characteristics. Iteration consists of the reduction of the region size through the geometric transformation and the obtaining of a spectral response region by fast Fourier transform so as to extract the blue noise characteristics from the spectral response. For regions with blue noise characteristics,the intensity distribution of forest vegetation texture is computed,and the texture size is calculated based on the current size of the area. Experimental results show that the gray scale and the distribution of forest vegetation texture units can be accurately measured,which lays reliable foundation for further texture segmentation.

remote sensing image; texture segmentation；blue noise; forest vegetation

2014-03-26；

2014-05-14

國家自然科學基金項目(編號：41271422)和遼寧省教育廳自然科學基金項目(編號：L2012379)共同資助。

10.6046/gtzyyg.2015.02.10

劉小丹，楊燊.基于藍噪聲理論的遙感圖像森林植被紋理測量[J].國土資源遙感，2015，27(2)：63-68.(Liu X D,Yang S.Forest vegetation texture measurement of remote sensing images based on the blue noise theory[J].Remote Sensing for Land and Resources,2015,27(2)：63-68.)

TP 751.1

A

1001-070X(2015)02-0063-06

劉小丹(1957-)，男，教授，主要從事數字圖像處理和數字印刷技術方面的研究。Email：liuxiao8866@sina.com。

楊燊(1988-)，男，碩士研究生，主要從事數字圖像處理方面的研究。Email：yangshen310@126.com。