基于灰度數學形態學的遙感影像薄云處理方法

曹 爽,岳建平,馬 文

(1.河海大學地球科學與工程學院測繪科學與工程系,江蘇南京 210098;2.南京信息工程大學遙感學院測繪工程系,江蘇南京 210044)

基于灰度數學形態學的遙感影像薄云處理方法

曹 爽1,2,岳建平1,馬 文1,2

(1.河海大學地球科學與工程學院測繪科學與工程系,江蘇南京 210098;2.南京信息工程大學遙感學院測繪工程系,江蘇南京 210044)

遙感成像過程中薄云是最常見的一種噪聲。在利用數學形態學進行薄云處理的基礎上,探討利用廣義多結構元素的數學形態學進行薄云處理的方法,并對 SPOT圖像進行試驗。結果表明,該方法對薄云去除的處理效果有一定的改善和提高。

數學形態學;廣義多結構元素;遙感影像;薄云去除

一、引 言

從遙感影像中獲取信息是學術界一直關注的課題之一。遙感成像過程中,云霧 (主要是薄云、薄霧)是最常見的一種噪聲,它使獲得的遙感圖像模糊,分辨率下降,有用遙感信息相對貧乏[1]。去云處理不僅是對遙感圖像進行準確解譯的基礎,也是增強遙感數據有效性、可用性的重要途徑。它是遙感圖像預處理中的一個重要環節,也是圖像處理、圖像分析的一大難點。目前,除合成孔徑雷達傳感器能穿透云層獲取地表信息外,其他傳感器均未能徹底解決影像數據的云覆蓋問題。從 20世紀 90年代開始,遙感影像的云去除就成為一個活躍的研究領域。因此,運用圖像處理技術,研究如何有效去除云霧的影響、提高遙感數據利用率具有重要意義。本文在研究數學形態學進行薄云處理的基礎上,引入廣義結構元素,旨在討論一種新的薄云處理方法,更有效地去除薄云。

二、現有薄云處理方法

云覆蓋類型不同,消除云的方法也不同。薄云霧覆蓋區域,地面景物的光譜特征仍在一定程度上得以保留,對此可以通過大氣糾正或數字圖像處理的相關技術有效去除云霧影響,恢復地面景物信息。國內外對遙感圖像云層干擾的去除作了大量的研究,提出了多種算法。薄云霧覆蓋區信息恢復的方法與途徑主要有以下幾種：

1)基于光譜信息的薄云處理方法。這種方法對研究條件要求比較高,如果利用多光譜影像中某些波段對薄云的較強敏感性來提取薄云信息,有時會使遙感影像的某些波段缺失,且處理后的波段與原波段對應關系不明確,如纓帽變換法[2]。

2)基于影像融合的薄云處理方法。該方法操作簡便,但是去薄云效果不是很理想,局限性很大[3]。

3)基于同態濾波的薄云處理方法。該方法利用同態濾波以及改進的同態濾波方法[4],根據薄云霧覆蓋信息在頻率域上通常占據低頻信息的特點,將圖像通過傅里葉變換轉換到頻率域,然后使用高通濾波器對影像進行空間濾波,除去云霧等低頻信息。該方法忽略了影像的低頻成分除薄云霧以外,還包含大量的其他背景信息的事實,因此,在去除低頻云霧成分的同時,該方法也會削弱影像的其他背景信息,從而改變了一些區域的影像特性。

4)基于混合像元分解的薄云處理方法[5]。該方法充分考慮水汽環境的特點,把水體像元光譜看作水、污染物和氣溶膠(薄云)等光譜的混合。混合像元模型能有效地消除薄云和氣溶膠的影響,可實現更精確的水質信息提取。不過此方法主要針對水域影像薄云去除,應用受到獲取數據等的限制。

5)基于直方圖匹配的薄云處理方法[6]。該方法首先假設薄云覆蓋區去除云影響后的圖像特征與無云覆蓋區圖像特征相同,因此可以利用同一幅圖像中非云覆蓋區的圖像特征作為參照,使云覆蓋區的圖像直方圖與之匹配,達到消除云影響的目的。該方法不足之處在于假設條件有時并不滿足,精確的云覆蓋區范圍的確定也有一定困難。

6)基于數學形態學的薄云處理方法。該方法利用數學形態學中的開、閉運算相結合可以分離高低頻的特點,構造結構元素,進行薄云處理[7]。

三、基于灰度數學形態學的薄云處理

1.灰度數學形態學薄云處理的基本方法

數學形態學 (mathematical morphology)誕生于1964年,它應用最廣泛的一個方向就是濾波作用,如將開、閉運算結合起來構成形態學噪聲濾波器,而且數學形態學的并行運算對處理數字圖像有很重要的意義,可以大大提高處理的速度。數學形態學中的基本運算是基于二值圖像進行定義的,但實際遙感圖像處理是針對灰度圖像 (多值圖像)進行的,因此,處理灰度圖像時采用灰度數學形態學。利用灰度數學形態學的濾波算法來處理遙感影像中處于低頻的薄云的基本步驟如下：

1)利用結構元素對含有薄云影像進行先閉后開運算,分離影像的高低頻成分,即

式中,x、y為影像坐標;F(x,y)為薄云影像的灰度函數;L(x,y)為影像的低頻部分;fc為閉運算;fo為開運算;g為結構元素,一般取 3×3或 5×5的矩陣。

2)L(x,y)處理后獲得的是薄云影像的低頻部分。其中不但有云信息也有地物信息,只要選取系數 K將L(x,y)中的云分量求出,再由原影像將這部分減去,就可以得到去薄云后的影像。即

式中,F′(x,y)為去薄云后的影像;K為云所占低頻的百分數,根據影像的薄云情況經驗選取。

根據試驗比較可知,灰度數學形態學中的先閉后開運算應用到遙感影像的去薄云處理中,可以達到去云的目的,云覆蓋區域下的地物紋理更清楚,影像細節保持更好,目視效果與傳統去云效果相比較也很理想。

2.數學形態學中的廣義多結構元素

在數學形態學中,結構元素是最重要、最基本的概念,結構元素在形態變換中的作用相當于信號處理中的“濾波窗口”,就如卷積運算中的模板。不同的場合結構元素的選擇及其相應的處理算法是不一樣的,對不同的目標圖像須設計不同的結構元素和不同的處理算法。結構元素的大小 、形態選擇合適與否,直接影響去云的效果。

開運算可以濾除圖像中的峰值(正脈沖)噪聲,閉運算可以濾除圖像中的低谷 (負脈沖)噪聲。傳統的開運算和閉運算只采用了一種結構元素,因此濾波輸出的圖像中只保留了與這一結構元素相匹配的幾何信息,而其他的幾何信息與噪聲一起被濾掉,所以不利于圖像幾何特征的保持。為了提高數學形態學濾波的有效性,采用多個結構元素并行處理的數學形態學濾波在圖像濾波和特征提取方面獲得應用。趙春暉[8]采用了多結構元素,構造了一類多結構元素并行復合形態濾波器。這類濾波器不僅可以有效地抑制圖像中的噪聲,而且較好地保持了圖像的幾何結構特征。

在多結構元素形態學的基礎上,還可以采用多結構元素構造一類適于圖像平滑并行復合形態的濾波器,分別稱為廣義形態開-閉最大濾波器和閉-開最小濾波器：

通過研究證明,傳統的閉-開濾波模糊了圖像的細節特征,濾波效果不好;普通的多結構元素雖然保持了細節,但不能完全達到濾波效果;而這種廣義的形態閉-開最小濾波器不僅可以有效地抑制圖像中的噪聲,而且較好地保持了圖像的幾何結構特征[8]。

3.廣義多結構元素在薄云處理中的應用

多結構元素的數學形態學應用到薄云處理中的基本公式為

式中,F(x,y)、L(x,y)、fc、fo含義同式 (1);g、b為不同形態的結構元素,且相對來說 g的尺寸比 b的尺寸小。

按多結構元素數學形態學薄云處理的基本公式,將廣義形態閉-開最小濾波器應用于去除薄云中,相當于用兩組結構元素組做不同結構元素不同尺寸的多結構元素薄云處理,之后再將兩組薄云處理后的圖在每一像元點進行比較,取最小值,得到的就是按廣義形態閉-開最小濾波器薄云處理后的圖像。例如:選用兩組結構元素集,每個結構元素集中有兩個結構元素,其中, B11為尺寸是 9像素的菱形結構元素,B12為尺寸是 9像素的方形結構元素,第二組元素組中的結構元素是比第一組結構元素尺寸大 5像素的一組結構元素,B21是方形結構元素,B22為菱形結構元素。分別用 B11與 B21、B12與 B22按式 (5)做閉-開運算,求出兩組薄云處理后的圖像,再對兩組圖像求最小值。

四、試驗結果與分析

試驗采用如圖 1所示的 SPOT影像,接收時間為 10月份,該影像受到薄云以及厚云遮擋的影響,薄云最大的遮擋區域達到 166像素 ×200像素。選取薄云區域,利用Matlab編程,為了對比分析,對薄云區域分別用同態濾波(截止頻率為 0.3),多結構元素(先尺寸為 5的菱形后尺寸為 9的方形結構元素)及廣義結構元素(按第三章第 3節中尺寸)的數學形態學方法進行處理,處理結果如圖2～圖 4所示。

圖1 原始圖

圖2 同態濾波去薄云

圖3 多結構元素去薄云

圖4 廣義結構元素去薄云

從目視角度可以看出：①同態濾波方法使得薄云區域邊緣銳化;灰度數學形態學中的先閉后開運算應用到遙感影像的薄云處理中,可以達到去云作用,云覆蓋區域下的地物紋理清楚,影像細節保持更好,目視效果與傳統去云效果相比較也很理想。②廣義形態閉-開最小濾波器去薄云后的效果優于不同形態不同尺寸的多結構元素薄云處理的效果,在細節保持上更好,圖4整體清晰度要比圖3好,薄云處理后圖4受結構元素形狀的影響比圖 3要小。通過比較可知,采用廣義結構元素去云效果是最好的。

在遙感影像的空間域中,含薄云區域具有局部能量大、區域灰度平均值高、標準差小的特點,因此,本文對處理后的云影像進行統計分析,主要采用的評價指標是影像的均值、標準差和熵[13-14]。設影像為 f(x, y),大小為M×N,則各指標可如下定義

云影像均值要大于薄云處理后影像的均值,因此可用平均灰度評價去除薄云是否有效。標準差反映影像的信息含量,在遙感影像薄云處理中可以反映圖像的細節信息變化情況。熵反映圖像的信息含量,一般認為圖像的熵越高越好。理想的薄云處理結果是,去薄云后的影像平均灰度降低,標準差和熵的增大表明所含信息量豐富。處理后的薄云影像數據分別按均值、標準差、熵進行統計,如表1所示。

表 1 薄云處理前后影像統計指標對比

從表 1可以看出各種薄云處理方法都使去除薄云后的影像的平均灰度值下降,說明各方法都達到了薄云去除的目的,本文所提出的數學形態學各種薄云去除方法也可以進行薄云處理。為達到理想的目視效果,各方法在處理時都選擇無云參考區域進行直方圖規定化,因此各方法處理后影像的平均灰度值比較接近。各種方法處理后影像的標準差都變大,說明處理后影像都變清晰。從表 1中的薄云影像及幾種薄云處理后影像的熵統計值對比可以看出,幾種方法處理后的圖像的熵值與原影像的熵值相比較都變小,說明各種方法處理后影像的標準差都變大,熵值與原影像的熵值相比較都變小,說明各種方法在一定程度上都會去掉影像的細節信息。其中,數字形態學去薄云后的影像的熵值比同態濾波去薄云后影像的熵值大,說明數學形態學去薄云方法比同態濾波方法更好地保持了影像的細節信息,而且廣義結構元素的數學形態學的熵值最大,說明這種方法保留的細節最豐富。

五、結束語

通過目視評價及統計分析可見,數學形態學方法進行去薄云處理是可行的,并且在選取合適的結構元素,特別是構造廣義結構元素后利用數學形態學處理薄云,可以更好地保持處理區域的細節,效果更佳。但本文只根據試驗影像的特點對結構元素進行了選取,而這種結構元素并不適合所有影像,在實踐中應針對特定的影像選取合適的結構元素進行薄云處理,才能取得理想效果。因此,利用本文提出的方法進行薄云處理過程中,結構元素的選取仍在探討之中,需要進一步深入研究。

[1] 朱錫芳,吳峰,孫燕濱.基于Mallat算法遙感圖像去云霧處理的改進方法 [J].遙感學報,2007,11(2)：241-246.

[2] GAO B C,KAUFMAN Y J,HAN W,et al.Correction of Thin Cirrus Path Radiance in the 0.4-1.0μm Spectral Region Using the Sensitive 1.375μm Cirrus Detecting Channel[J].Journal of Geophysical Research,1998,103 (D24)：32169-32176.

[3] ZΗANG Y.A NewMergingMethod and Its Spectral and Spatial Effects[J].IntJ Remote Sensing,1999,20(10)：2003-2014.

[4] 趙忠明,朱重光.遙感圖像中薄云的去除方法[J].環境遙感,1996,11(3)：195-199.

[5] 吳傳慶,王橋,楊志峰.基于混合像元分解的水體遙感圖像去云法[J].遙感學報,2006,10(2)：34-41.

[6] 李月臣,陳晉,劉春霞,等.一種除去遙感影像薄云霧的方法[J].成都理工大學學報：自然科學版,2006, 33(1)：58-63.

[7] 曹爽,李浩,馬文.基于數學形態學的遙感影像薄云處理方法[J].地理與地理信息科學,2009,25(4)：30-33.

[8] 趙春暉,李兵,惠俊英.一種適于圖像噪聲平滑的廣義形態濾波器[J].哈爾濱工程大學學報,2000,21(2)：55-59.

[9] 曹爽.高分辨率遙感影像去云處理研究[D].南京：河海大學,2006.

[10] GONZALEZ R C,WOODS R E,EDD INS S L.數字圖像處理：Matlab版[M].阮秋琦,等譯.北京：電子工業出版社,2005.

Removing Thin Cloud Arithmetic Based on GrayMathematic M orphology for Remote Sensing Images

CAO Shuang,YUE Jianping,MA Wen

0494-0911(2010)12-0054-04

P237

B

2010-08-23

曹 爽(1977—),女,遼寧彰武人,博士生,講師,主要研究方向為攝影測量和遙感數字圖像處理。