極化SAR圖像自適應(yīng)增強(qiáng)Lee濾波算法

郎豐鎧,楊 杰,李德仁

武漢大學(xué)測繪遙感信息工程國家重點(diǎn)實(shí)驗(yàn)室,湖北武漢 430079

極化SAR圖像自適應(yīng)增強(qiáng)Lee濾波算法

郎豐鎧,楊 杰,李德仁

武漢大學(xué)測繪遙感信息工程國家重點(diǎn)實(shí)驗(yàn)室,湖北武漢 430079

精制極化Lee濾波算法在其應(yīng)用中存在明顯的缺陷:扇貝效應(yīng)和虛假細(xì)線。對此,本文提出增加一組均勻窗口及一組線性方向窗口,并采用大小自適應(yīng)窗口機(jī)制,在同質(zhì)度高的區(qū)域用大窗口濾波,在同質(zhì)度低的區(qū)域用小窗口濾波,從而使得濾波窗口在形狀和尺寸上都能自動適應(yīng)實(shí)際場景。利用機(jī)載和星載全極化SAR數(shù)據(jù)進(jìn)行的濾波試驗(yàn)結(jié)果表明,本文提出的自適應(yīng)增強(qiáng)Lee濾波算法在同質(zhì)區(qū)域的噪聲抑制能力明顯優(yōu)于精制極化Lee濾波算法及改進(jìn)的Sigma濾波算法,同時(shí)在保持點(diǎn)、線等細(xì)節(jié)信息方面也優(yōu)于精制極化Lee濾波算法,并且能很好地保持圖像中地物的極化散射信息。

極化SAR;自適應(yīng);相干斑;噪聲;濾波

1 引 言

由于SAR系統(tǒng)采用相干成像,SAR圖像中不可避免地存在相干斑噪聲。相干斑噪聲降低了數(shù)據(jù)的信噪比,令SAR圖像解譯、參數(shù)反演等工作變得更加困難。為了減小相干斑噪聲的影響,通常在進(jìn)一步處理之前首先對SAR圖像進(jìn)行相干斑噪聲濾波,以盡可能地抑制噪聲,同時(shí)盡量保持圖像的細(xì)節(jié)信息。因此,自從SAR系統(tǒng)產(chǎn)生以來,SAR圖像濾波算法的研究就從未停止過,并且到目前為止仍然是一個(gè)熱點(diǎn)[1-9]。

相對于單極化SAR,全極化SAR系統(tǒng)可以獲取HH、HV、VH和VV共4個(gè)通道的數(shù)據(jù),從而提供更加豐富的全極化信息。因此全極化SAR系統(tǒng)正在逐漸代替單極化SAR在地表分類、目標(biāo)識別、變化監(jiān)測、參數(shù)反演等方面扮演越來越重要的角色[20-24]。但同時(shí),對于極化SAR來說,相干斑濾波不僅要分別考慮4個(gè)通道的濾波,還要考慮4個(gè)通道之間的相關(guān)性,因此相比于單極化SAR,全極化SAR圖像的相干斑抑制問題更加復(fù)雜[4]。盡管如此,很多經(jīng)典的極化SAR濾波算法已經(jīng)被提出。文獻(xiàn)[1]提出極化白化濾波器,通過將各極化通道最優(yōu)加權(quán)組合,獲得單通道的降噪圖像,損失了極化信息。文獻(xiàn)[2]提出最優(yōu)加權(quán)濾波器,文獻(xiàn)[3]提出極化矢量濾波器,均利用了乘性噪聲模型和最小均方差(MMSE)思想,但是輸入輸出均為多通道的強(qiáng)度或幅度圖像,不能對協(xié)方差矩陣的非對角線元素進(jìn)行濾波。針對以上問題,文獻(xiàn)[4]提出精制極化Lee濾波算法,采用非方形的邊緣方向窗口達(dá)到保持邊緣信息的目的,并且對協(xié)方差矩陣的所有元素利用相同的參數(shù)進(jìn)行濾波,以防止通道間串?dāng)_。這種防止通道間串?dāng)_的思想對極化SAR濾波算法的研究影響深遠(yuǎn),成為極化SAR濾波的一個(gè)基本原則[20]。雖然之后人們又提出很多優(yōu)秀的極化SAR濾波算法[5-9],但是,精制極化Lee濾波仍然是目前應(yīng)用較廣泛的濾波算法[6,21,23]。然而,此濾波算法存在明顯的缺陷:扇貝效應(yīng)和虛假細(xì)線[8-9]。因此對于有些應(yīng)用如道路提取、圖像分割等會造成較大影響。目前已經(jīng)有一些改進(jìn)的Lee濾波算法[13,16],這些算法的改進(jìn)思路主要有兩個(gè):①增加同質(zhì)像素選擇步驟[14-16];②令濾波窗口大小自適應(yīng)[13]。第1種思路在非方形窗口中進(jìn)一步篩選同質(zhì)像素,會降低噪聲抑制能力,增加算法復(fù)雜度;第2種思路僅能改善虛假細(xì)線問題,但不能從根本上克服這兩種缺陷。本文針對這兩種缺陷提出一種新的自適應(yīng)增強(qiáng)Lee濾波算法,旨在保持原始精制極化Lee濾波簡單、高效、健壯等優(yōu)點(diǎn)的同時(shí),克服扇貝效應(yīng)和虛假細(xì)線的問題。

2 精制極化Lee濾波算法

精制極化Lee濾波算法主要包括兩步:模板窗口選擇和局部統(tǒng)計(jì)濾波。具體步驟為:

(1)模板窗口選擇。將濾波窗口分解為3×3 的9個(gè)子窗格,計(jì)算各子窗格內(nèi)的總功率均值,然后采用如圖1所示4個(gè)邊緣梯度算子進(jìn)行檢測,每個(gè)邊緣梯度算子對應(yīng)圖1(e)—(l)中每列的2個(gè)邊緣方向窗口。選擇梯度最大的算子所對應(yīng)的兩個(gè)邊緣方向窗口中均值與中心像素均值比較接近的那個(gè)。

(2)局部統(tǒng)計(jì)濾波。在邊緣方向窗口中,利用線性最小均方差(LMMSE)濾波器進(jìn)行濾波。LMMSE的數(shù)學(xué)表達(dá)式為

圖1 邊緣梯度算子與3×3邊緣方向窗口Fig．1 Edge gradient operators and 3×3 edge-aligned windows

對于全極化SAR數(shù)據(jù),首先用Span圖像計(jì)算權(quán)重b,然后對相干矩陣或協(xié)方差矩陣的所有元素用相同的權(quán)重及濾波窗口進(jìn)行濾波。

3 自適應(yīng)增強(qiáng)Lee濾波算法

3．1 新模板設(shè)計(jì)

經(jīng)過試驗(yàn)分析發(fā)現(xiàn),由于原始的精制極化Lee濾波僅定義了一組邊緣方向窗口,在濾波時(shí)將均勻區(qū)域與非均勻區(qū)域都在邊緣方向窗口內(nèi)進(jìn)行濾波。在邊緣區(qū)域?yàn)V波效果很好,但是在均勻區(qū)域,由于濾波窗口的原因使得濾波后反而變得不再均勻,從而導(dǎo)致扇貝效應(yīng);另一方面,當(dāng)濾波窗口大于圖像中的線性目標(biāo)寬度時(shí),由于邊緣方向窗口不能正確擬合線性目標(biāo),無論最終選擇哪個(gè)邊緣方向窗口,窗口內(nèi)總存在異質(zhì)像素,從而導(dǎo)致虛假細(xì)線的產(chǎn)生。針對以上第1個(gè)問題,可以考慮增加一個(gè)方形窗口;針對第2個(gè)問題,則需要增加一組線性方向窗口。本文所采用的線性方向窗口如圖2(a)—(d)所示,它們又分別對應(yīng)圖2(e)—(l)中每列的兩個(gè)線性梯度算子。在計(jì)算時(shí),將兩個(gè)算子所得梯度中較大的那個(gè)作為最終線性方向窗口的梯度。

3．2 自適應(yīng)窗口

增加線性方向窗口后,精制極化Lee濾波效果有所改善,但是當(dāng)選擇較大窗口時(shí),無論是邊緣還是細(xì)線都將不能再用直線近似擬合,此時(shí)方向窗口的適應(yīng)性大大下降,從而導(dǎo)致濾波效果的下降。為避免此問題,最好采用小窗口進(jìn)行濾波。但是另一方面,對于同質(zhì)區(qū)域,窗口越大,相干斑抑制能力越強(qiáng),濾波效果越好。為了克服此矛盾,本文采用大小自適應(yīng)的窗口:在同質(zhì)度較低的區(qū)域采用小窗口進(jìn)行濾波,在同質(zhì)度較高的區(qū)域采用大窗口進(jìn)行濾波。從而在保持圖像中細(xì)節(jié)信息的同時(shí),盡可能地提高相干斑抑制能力。

圖2 3×3線性方向窗口與線性梯度算子Fig．2 3×3 linear windows and linear gradient operators

從以上分析中可以看出,自適應(yīng)窗口技術(shù)的關(guān)鍵在于同質(zhì)度量度的選擇。目前提出的同質(zhì)度量度主要是變差系數(shù)(Co V)[17-19],其計(jì)算公式為

但是,Co V理論閾值的設(shè)置要用到等效視數(shù)(ENL),而一幅SAR圖像的ENL往往是未知的,需要事先估計(jì)。而且由于SAR圖像中每個(gè)像素對應(yīng)的ENL是不同的[25-26],而需要估計(jì)的ENL只是一個(gè)全局值,并不能代表所有像素,因此會導(dǎo)致偏差。另外,由于Co V的物理意義不夠直觀,使得手動設(shè)置閾值也比較困難。基于以上原因,本文不采用Co V作為同質(zhì)度度量,而提出用局部ENL作為同質(zhì)度度量。強(qiáng)度及幅度SAR圖像ENL的計(jì)算公式為

式中,Cv1表示單視SAR圖像噪聲的理論Co V值,對于強(qiáng)度圖,Cv1＝1,對于幅度圖,Cv1＝0．522 7;Cy表示SAR圖像局部Co V。從式(4)中可以看出, ENL是Co V的函數(shù),因此兩者在同質(zhì)度度量方面具有同等的效力。但是,由于ENL的物理意義比Co V更加直觀,因此在閾值設(shè)置時(shí)更加方便。

自適應(yīng)窗口技術(shù)的另一個(gè)關(guān)鍵問題是窗口變化順序,即窗口是由小到大變化還是由大到小變化。文獻(xiàn)[13]采用的是由小到大的變化策略,而本文采用由大到小的變化策略。理由如下,如果采用由小到大的變化策略,當(dāng)圖像噪聲較嚴(yán)重時(shí),小窗口計(jì)算的ENL是不可信的,而且往往是偏低的,此時(shí)會導(dǎo)致窗口停止增長,反映在濾波結(jié)果上,圖像中會出現(xiàn)很多明顯的椒鹽狀斑點(diǎn)。而反過來,如果采用由大到小的變化策略,先計(jì)算大窗口的ENL,則在同質(zhì)區(qū),即使噪聲比較大,也能保證計(jì)算出的ENL的準(zhǔn)確度,因此與小窗口計(jì)算的ENL相比會比較大,從而會選擇大窗口進(jìn)行濾波,避免了椒鹽斑點(diǎn)的出現(xiàn)。

3．3 新的濾波算法

基于前面的分析,本文提出一種新的自適應(yīng)增強(qiáng)Lee濾波算法,其流程如圖3所示。具體計(jì)算步驟如下:

(1)設(shè)置最小濾波窗口Wmin和最大濾波窗口Wmax,以及ENL閾值。由于不同大小的窗口所估計(jì)的ENL范圍不同,因此其對應(yīng)的閾值也應(yīng)該不同,所以ENL閾值實(shí)際上是一系列閾值的集合,每個(gè)閾值對應(yīng)一個(gè)窗口。令當(dāng)前窗口大小W＝Wmax。

(2)利用式(4)計(jì)算W×W矩形窗口內(nèi)的局部ENL,并判斷是否小于閾值。如果是,則表明此窗口內(nèi)像素同質(zhì)度較低,轉(zhuǎn)到第(3)步;如果否,則表明此窗口內(nèi)像素同質(zhì)度較高,因此選擇方形窗口,轉(zhuǎn)到第(6)步。

(3)判斷是否當(dāng)前窗口大小W＝Wmin。如果是,則選擇非方形窗口進(jìn)行濾波,轉(zhuǎn)到第(4)步;如果否,令W＝W－2,轉(zhuǎn)到第(2)步。

(4)將窗口分解為3×3的9個(gè)子窗口,分解方法采用精制極化Lee濾波算法中的方法,然后利用圖1和圖2中共12個(gè)梯度算子計(jì)算梯度,并求出最大梯度所對應(yīng)的梯度算子代號ID。

(5)判斷是否ID∈{a,b,c,d}。如果是,則說明窗口位于邊緣,此時(shí)選擇梯度算子對應(yīng)的兩個(gè)邊緣方向窗口中均值與中心像素均值較接近的那個(gè)。如果否,則說明窗口位于線上,需要進(jìn)一步判斷:如果ID∈{e,g,i,k},則選擇線性方向窗口(圖2(a)、圖2(c))中均值與中心像素均值較接近的那個(gè);如果ID∈{f,h,j,l},則選擇線性方向窗口(圖2(b)、圖2(d))中均值與中心像素均值較接近的那個(gè)。

(6)在所選擇的窗口內(nèi)利用LMMSE濾波器對極化SAR數(shù)據(jù)進(jìn)行濾波。

圖3 算法流程圖Fig．3 Flowchart of the proposed algorithm

4 試驗(yàn)及結(jié)果分析

為了驗(yàn)證本文提出的濾波算法的有效性,本文分別采用機(jī)載和星載全極化SAR數(shù)據(jù)進(jìn)行試驗(yàn),對比分析Boxcar濾波器、精制極化Lee濾波器、改進(jìn)的Sigma濾波器[8]及自適應(yīng)增強(qiáng)Lee濾波器間的優(yōu)劣。

4．1 機(jī)載數(shù)據(jù)

首先采用機(jī)載全極化SAR數(shù)據(jù)進(jìn)行試驗(yàn)。本文所使用的數(shù)據(jù)為DLR ESAR系統(tǒng)在德國Oberpfaffenhofen試驗(yàn)區(qū)獲取的L波段全極化SAR數(shù)據(jù)。原始數(shù)據(jù)大小為2861像素×1540像素,為了方便展示,本文首先對其在方位向進(jìn)行了2視處理,然后截取其中540像素×888像素大小的一塊區(qū)域進(jìn)行濾波試驗(yàn)。如圖4所示,為所截取區(qū)域的Pauli-RGB合成圖像。從圖中可以看出,試驗(yàn)區(qū)地物覆蓋類型非常豐富,點(diǎn)、線、邊緣等細(xì)節(jié)信息較多,并且圖像中有較強(qiáng)的相干斑噪聲。這些因素給SAR圖像處理解譯帶來很大困難,反過來,又有利于評價(jià)濾波算法在噪聲抑制、細(xì)節(jié)信息保持等方面的優(yōu)劣。

由于精制極化Lee濾波器采用邊緣方向窗口,因此其濾波窗口中所包含的實(shí)際像素?cái)?shù)比同尺度窗口的Boxcar濾波器所包含的像素?cái)?shù)要少。為達(dá)到相當(dāng)?shù)南喔砂咭种菩Ч?精制極化Lee濾波器的窗口需要設(shè)置的比Boxcar濾波器的窗口略大。并且在選擇窗口大小時(shí)為了保持細(xì)節(jié)信息,往往設(shè)置為與圖像中點(diǎn)、線等目標(biāo)的尺寸相當(dāng)。對于改進(jìn)的Sigma濾波器和自適應(yīng)增強(qiáng)Lee濾波器,由于其實(shí)際參與濾波的像素?cái)?shù)是自適應(yīng)的,因此可以選擇較大的濾波窗口,不會導(dǎo)致細(xì)節(jié)信息的嚴(yán)重丟失。本小節(jié)中考慮試驗(yàn)數(shù)據(jù)中點(diǎn)、線等細(xì)節(jié)信息的尺寸,分別用5×5的Boxcar濾波器、7×7的精制極化Lee濾波器、11×11的改進(jìn)的Sigma濾波器及最小窗口為3×3、最大窗口為11×11的自適應(yīng)增強(qiáng)Lee濾波器對試驗(yàn)數(shù)據(jù)進(jìn)行濾波處理,所得部分結(jié)果如圖5所示。

從圖中可以看出,精制極化Lee濾波結(jié)果中有明顯的扇貝效應(yīng)(圖5(a)中區(qū)域1),并且在暗線中心有虛假亮線(圖5(a)中區(qū)域2)。而自適應(yīng)增強(qiáng)Lee濾波結(jié)果在同質(zhì)區(qū)域?yàn)V波效果較平滑,并且能很好地保持點(diǎn)、線等細(xì)節(jié)信息(圖5(b))。為了定量評價(jià)濾波器的邊緣保持能力,根據(jù)以下邊緣保持指數(shù)(EPI)[27]公式對幾種濾波結(jié)果進(jìn)行比較

表1 濾波前后邊緣保持指數(shù)Tab．1 EPI of the original and filtered images

利用Span圖計(jì)算得到的結(jié)果如表1所示。從表1可以看出,自適應(yīng)增強(qiáng)Lee濾波的EPI比Boxcar和精制極化Lee濾波都高。但是改進(jìn)的Sigma濾波EPI卻遠(yuǎn)遠(yuǎn)高于其他濾波器。其中部分原因是因?yàn)楦倪M(jìn)的Sigma濾波對散射回波較高的像素不進(jìn)行濾波處理,因此與原始圖像是一樣的。從式(5)中可以看出,如果濾波前后像素值不變,則EPI為1,從而這些未濾波的像素使EPI估計(jì)值偏高。另外,從表2中可以看出,改進(jìn)的Sigma濾波噪聲抑制能力比自適應(yīng)增強(qiáng)Lee濾波弱,因此其高EPI是以弱噪聲抑制力為代價(jià)的。

為了進(jìn)一步評價(jià)濾波器在噪聲抑制方面的能力,選擇如圖4所示3個(gè)矩形區(qū)域,利用Span圖統(tǒng)計(jì)其濾波前后的Co V值,統(tǒng)計(jì)結(jié)果見表2。從表2中可以看出,精制Lee濾波與改進(jìn)的Sigma濾波的噪聲抑制效果大致相當(dāng),略優(yōu)于Boxcar濾波,而本文提出的自適應(yīng)增強(qiáng)Lee濾波所得Co V值明顯低于其他濾波結(jié)果,具有最優(yōu)的噪聲抑制能力。

表2 濾波前后CoV對比表Tab．2 CoV comparison of the original and filtered images

4．2 星載數(shù)據(jù)

本小節(jié)利用ALOS PALSAR系統(tǒng)對武漢市獲取的全極化SAR數(shù)據(jù)進(jìn)行進(jìn)一步試驗(yàn)。原始數(shù)據(jù)大小為18 432像素×1248像素,為了方便展示,本文首先對其在方位向進(jìn)行了8視處理,然后截取其中460像素×640像素大小的一塊區(qū)域進(jìn)行濾波試驗(yàn)。如圖6所示,為所截取區(qū)域的Pauli-RGB合成圖像。從圖中可以看到,試驗(yàn)區(qū)地物覆蓋情況非常復(fù)雜,河流、湖泊、道路縱橫交錯(cuò),線性目標(biāo)非常豐富,并且圖像中同樣伴有明顯的相干斑噪聲。

考慮試驗(yàn)數(shù)據(jù)中點(diǎn)、線等細(xì)節(jié)信息的尺寸,分別用3×3的Boxcar濾波器、5×5的精制極化Lee濾波器、7×7的改進(jìn)的Sigma濾波器及最小窗口為3×3、最大窗口為7×7的自適應(yīng)增強(qiáng)Lee濾波器對試驗(yàn)數(shù)據(jù)進(jìn)行濾波處理,所得部分結(jié)果如圖7所示。從圖中可以看出,精制極化Lee濾波結(jié)果在暗線中心有明顯的虛假亮線(圖7(a)中區(qū)域1),并且較小的水體(可視為暗的點(diǎn)目標(biāo))及較細(xì)的亮線會被模糊掉(圖7(a)中區(qū)域2)。而自適應(yīng)增強(qiáng)Lee濾波能很好地保持暗線、亮線、暗的點(diǎn)目標(biāo)等細(xì)節(jié)信息,并且在同質(zhì)區(qū)域?yàn)V波效果比精制極化Lee濾波更加平滑(圖7(b))。

利用Span圖計(jì)算濾波前后EPI,得到的結(jié)果如表3所示。從表3可以看出,自適應(yīng)增強(qiáng)Lee濾波的EPI遠(yuǎn)高于精制極化Lee濾波,但低于改進(jìn)的Sigma濾波。

表3 濾波前后邊緣保持指數(shù)Tab．3 EPI of the original and filtered images

為了進(jìn)一步評價(jià)自適應(yīng)增強(qiáng)Lee濾波對極化散射特性的保持效果,選取圖6中所示的一小塊同質(zhì)區(qū)域,繪制其濾波前后平均相干矩陣的極化響應(yīng)圖,結(jié)果如圖8所示。從圖中可以看出,濾波前后共極化與交叉極化的極化響應(yīng)圖均有高度的一致性,說明自適應(yīng)增強(qiáng)Lee濾波器能很好地保持圖像中地物的極化散射特性。

圖4 ESAR L波段全極化SAR數(shù)據(jù)Pauli-RGB合成圖Fig．4 Pauli-RGB image of ESAR L-band PolSAR data

圖5 濾波結(jié)果對比圖Fig．5 Comparison of different filters

圖6 PALSAR L波段全極化SAR數(shù)據(jù)Pauli-RGB合成圖Fig．6 Pauli-RGB image of PALSAR L-band PolSAR data

圖7 濾波結(jié)果對比圖Fig．7 Comparison of different filters

圖8 濾波前后共極化與交叉極化的極化響應(yīng)圖Fig．8 The original and filtered co-polarization and crosspolarization signatures

5 總 結(jié)

本文針對精制極化Lee濾波結(jié)果中存在扇貝效應(yīng)和虛假細(xì)線的問題,提出增加一組均勻窗口及一組線性方向窗口,并采用大小自適應(yīng)窗口機(jī)制,在同質(zhì)度高的區(qū)域用大窗口濾波,在同質(zhì)度低的區(qū)域用小窗口濾波,從而使得濾波窗口在形狀和尺寸上都能自動適應(yīng)實(shí)際場景。利用機(jī)載和星載全極化SAR數(shù)據(jù)進(jìn)行的濾波試驗(yàn)結(jié)果表明,本文提出的自適應(yīng)增強(qiáng)Lee濾波算法在同質(zhì)區(qū)域的噪聲抑制能力明顯強(qiáng)于精制極化Lee濾波算法及改進(jìn)的Sigma濾波算法,同時(shí)在保持較暗的與較亮的點(diǎn)、線等細(xì)節(jié)信息方面也優(yōu)于精制極化Lee濾波算法,并且能很好地保持圖像中地物的極化散射信息。

然而,由于窗口自適應(yīng)機(jī)制需要設(shè)置更多的參數(shù),尤其是同質(zhì)度閾值參數(shù)的設(shè)置,使得自適應(yīng)增強(qiáng)Lee濾波算法在易用性上面不如精制極化Lee濾波。后續(xù)將研究同質(zhì)度參數(shù)閾值的自適應(yīng)設(shè)置方法,提高算法的易用性。

[1] NOVAK L M,BURL M C．Optimal Speckle Reduction in Polarimetric SAR Imagery[J]．IEEE Transactions on Aerospace and Electronic Systems,1990,26(2): 293-305．

[2] LEE J,GRUNES M R,MANGO S A．Speckle Reduction in Multipolarization,Multifrequency SAR Imagery[J]．IEEE Transactions on Geoscience and Remote Sensing, 1991,29(4):535-544．

[3] GOZE S,LOPES A．A MMSE Speckle Filter for Full Resolution SAR Polarimetric Data[J]．Journal of Electromagnetic Waves and Applications,1993,7(5):717-737．

[4] LEE J S,GRUNES M R,DE GRANDI G．Polarimetric SAR Speckle Filtering and Its Implication for Classification[J]．IEEE Transactions on Geoscience and Remote Sensing,1999,37(5):2363-2373．

[5] VASILE G,TROUVE E,BUZULOIU V．Intensity-driven Adaptive-neighborhood Technique for Polarimetric and Interferometric SAR Parameters Estimation[J]．IEEE Transactions on Geoscience and Remote Sensing,2006, 44(6):1609-1621．

[6] LEE J,GRUNES M R,SCHULER D L,et al．Scatteringmodel-based Speckle Filtering of Polarimetric SAR Data [J]．IEEE Transactions on Geoscience and Remote Sensing,2006,44(1):176-187．

[7] LOPEZ-M ARTINEZ C,FABREGAS X．Model-based Polarimetric SAR Speckle Filter[J]．IEEE Transactions on Geoscience and Remote Sensing,2008,46(11): 3894-3907．

[8] LEE J,WEN J-H,AINSWORTH T L,et al．Improved Sigma Filter for Speckle Filtering of SAR Imagery[J]．IEEE Transactions on Geoscience and Remote Sensing, 2009,47(1):202-213．

[9] CHEN J,CHEN Y,AN W,et al．Nonlocal Filtering for Polarimetric SAR Data:A Pretest Approach[J]．IEEE Transactions on Geoscience and Remote Sensing,2011, 49(5):1744-1754．

[10] H UANG Shiqi,LIU Daizhi．Research on Method and Application of Speckle Noise Reduction of SAR Image[J]．Acta Geodaetica et Cartographica Sinica,2006,35(3): 245-250．(黃世奇,劉代志．SAR圖像斑點(diǎn)噪聲抑制方法與應(yīng)用研究[J]．測繪學(xué)報(bào),2006,35(3):245-250．)

[11] YANG Shenbin,LI Bingbai,SHEN Shuanghe,et al．Structure Retaining Linear Multi-channel SAR Image Speckle Filter [J]．Acta Geodaetica et Cartographica Sinica,2006,35 (4):364-370．(楊沈斌,李秉柏,申雙和,等．基于特征保持的線性多通道最優(yōu)求和SAR圖像濾波算法[J]．測繪學(xué)報(bào), 2006,35(4):364-370．)

[12] ZHANG Zhongshan,YU Jie,YAN Qin,et al．Research on Polarimetric SAR Image Speckle Reduction Using Kernel Independent Component Analysis[J]．Acta Geodaetica et Cartographica Sinica,2011,40(3):289-295．(張中山,余潔,燕琴,等．基于核獨(dú)立成分分析的極化SAR圖像相干斑抑制[J]．測繪學(xué)報(bào),2011,40(3): 289-295．)

[13] JIANG Yong,ZHANG Xiaoling,SHI Jun．Speckle Reduction for Polarimetric SAR Images by Improved Lee Filter[J]．Journal of University of Electronic Science and Technology of China,2009,38(1):5-8．(江勇,張曉玲,師君．極化SAR改進(jìn)Lee濾波相干斑抑制研究[J]．電子科技大學(xué)學(xué)報(bào),2009,38(1):5-8．)

[14] YANG Dahai,MA Debao．Refined Polarimetric SAR Speckle Lee Filtering Algorithm Based on Relativity of Polarization Vector[J]．Journal of Information Engineering University,2010,11(6):737-740．(楊大海,馬德寶．基于極化矢量相似系數(shù)的極化Lee濾波改進(jìn)算法[J]．信息工程大學(xué)學(xué)報(bào),2010,11(6):737-740．)

[15] OUYANG Qundong,WU Zhaocong,PENG Jiangui．An Improved Method of Lee Refined Polarized Filter[J]．Science of Surveying and Mapping,2011,36(5):136-138．(歐陽群東,巫兆聰,彭檢貴．一種改進(jìn)的精制極化Lee濾波算法[J]．測繪科學(xué),2011,36(5):136-138．)

[16] LIU Rong,LOU Xiaoguang．An Improved Speckle Reducing Algorithm for Lee Filtering Based on Edge Direction[J]．Science Technology and Engineering,2011,11(11): 2497-2501．(劉蓉,婁曉光．基于邊緣特性的極化Lee濾波改進(jìn)算法[J]．科學(xué)技術(shù)與工程,2011,11(11): 2497-2501．)

[17] NEZRY E,LOPES A,TOUZI R．Detection of Structural and Textural Features for SAR Images Filtering[C]∥IGARSS’91:Proceedings of IEEE 1991 International Geoscience and Remote Sensing Symposium．Espoo: IEEE,1991:2169-2172．

[18] NICOLASJ M,TUPIN F,MAITRE H．Smoothing Speckled SAR Images by Using Maximum Homogeneous Region Filters:An Improved Approach[C]∥IGARSS’01: Proceedings of IEEE 2001 International Geoscience and Remote Sensing Symposium．Sydney:IEEE,2001: 1503-1505．

[19] TOUZI R．A Review of Speckle Filtering in the Context ofEstimation Theory[J]．IEEE Transactions on Geoscience and Remote Sensing,2002,40(11):2392-2404．

[20] LEE J,POTTIER E．Polarimetric Radar Imaging:From Basic to Application[M]．New York:CRC Press,2009．

[21] YANG Jie,LANG Fengkai,LI Deren．An Unsupervised Wishart Classification for Fully Polarimetric SAR Image Based on Cloude-Pottier Decomposition and Polarimetric Whitening Filter[J]．Geomatics and Information Science of Wuhan University．2011,36(1):104-107．(楊杰,郎豐鎧,李德仁．一種利用Cloude-Pottier分解和極化白化濾波的全極化SAR圖像分類算法[J]．武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2011,36(1):104-107．)

[22] LANG Fengkai,YANG Jie,ZHAO Lingli,et al．Polarimetric SAR Data Classification with Freeman Entropy and Anisotropy Analysis[J]．Acta Geodaetica et Cartographica Sinica,2012,41(4):556-562．(郎豐鎧,楊杰,趙伶俐,等．基于Freeman散射熵和各向異性度的極化SAR影像分類算法研究[J]．測繪學(xué)報(bào),2012,41(4):556-562．)

[23] YANG Jie,ZHAO Lingli,SHI Lei,et al．Interpretation of Oblique Buildings Based on Optimal Polarimetric Coherence Coefficient[J]．Acta Geodaetica et Cartographica Sinica,2012,41(4):577-583．(楊杰,趙伶俐,史磊,等．基于最優(yōu)極化相干系數(shù)的傾斜建筑物解譯研究[J]．測繪學(xué)報(bào),2012,41(4):577-583．)

[24] TAN Lulu,YANG Libo,YANG Ruliang．Investigation of Tree Height Retrieval with Polarimetric SAR Interferometry Based on ESPRIT Algorithm[J]．Acta Geodaetica et Cartographica Sinica,2011,40(3):296-300．(談璐璐,楊立波,楊汝良．基于ESPRIT算法的極化干涉SAR植被高度反演研究[J]．測繪學(xué)報(bào),2011,40(3):296-300．)

[25] ANFINSEN S N,DOULGERIS A P,ELTOFT T．Estimation of the Equivalent Number of Looks in Polarimetric Synthetic Aperture Radar Imagery[J]．IEEE Transactions on Geoscience and Remote Sensing,2009,47(11): 3795-3809．

[26] CUI Y,ZHOU G,YANG J．Unsupervised Estimation of the Equivalent Number of Looks in SAR Images[J]．IEEE Geoscience and Remote Sensing Letters,2011,8(4): 710-714．

[27] HAN Chunming,GUO Huadong,WANG Changlin,et al．Edge-preserving Filter for SAR Images[J]．High Technology Letters,2003,13(7):11-15．(韓春明,郭華東,王長林,等．保持邊緣的SAR圖像濾波方法[J]．高技術(shù)通訊,2003, 13(7):11-15．)

(責(zé)任編輯:叢樹平)

An Adaptive Enhanced Lee Speckle Filter for Polarimetric SARImage

LANG Fengkai,YANG Jie,LI Deren
State Key Laboratory of Information Engineering in Surveying,Mapping and Remote Sensing,Wuhan University, Wuhan 430079,China

In its application,the refined Lee filter has two severe defects:scallop effect and false-line effect．In this paper,it is extended and improved by adding a square window and a group of linear windows．To better suppress the speckle noise as well as preserve the detail information,an adaptive window mechanism is adopted:big window is used in homogeneous areas and small window is used in heterogeneous areas．So the proposed filter can adapt actual scenes in both shape and size aspects．Both airborne and spaceborne PolSAR data are used to demonstrate its overall speckle filtering characteristics by comparing with other filters．The results show that the adaptive enhanced Lee filter has better performance in speckle suppressing and detail preserving than the refined Lee filter,and can preserve the polarimetric information very well．

polarimetric SAR;adaptive window;speckle;noise;filter

LANG Fengkai(1987—),male,PhD candidate,majors in PolSAR image processing．

P237

A

1001-1595(2014)07-0690-08

2013-03-07

郎豐鎧(1987—),男,博士生,研究方向?yàn)闃O化SAR圖像處理。

E-mail:lfkupc＠126．com

LANG Fengkai,YANG Jie,LI Deren．An Adaptive Enhanced Lee Speckle Filter for Polarimetric SAR Image[J]．Acta Geodaetica et Cartographica Sinica,2014,43(7):690-697．(郎豐鎧,楊杰,李德仁．極化SAR圖像自適應(yīng)增強(qiáng)Lee濾波算法[J]．測繪學(xué)報(bào), 2014,43(7):690-697．)

10．13485/j．cnki．11-2089．2014．0112

國家自然科學(xué)基金(60890074);國家863計(jì)劃(2011AA120404);中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(201161902020003)

修回日期:2013-07-09