基于方差系數(shù)統(tǒng)計的合成孔徑雷達(dá)相干斑強度估計

鳳宏曉 屈建社

(陜西黃河集團(tuán)有限公司 西安 710043)

0 引言

相干斑是SAR成像系統(tǒng)本身所固有的，在不同成像系統(tǒng)或參數(shù)下產(chǎn)生的相干斑的強度也是不同的。理想情況下，對同一場景的L幅單視圖像進(jìn)行平均時，就會獲得L視的SAR圖像，此時我們也將知道相應(yīng)的相干斑的方差值。但是，這需要假設(shè)多普勒帶寬分割得到的子帶寬之間是相互獨立的，而在實際中各個子帶寬之間往往存在一定的相關(guān)性［1-2］。同時，若采用空間平均的多視處理方法，由于SAR圖像內(nèi)部的像素之間存在很明顯的空間相關(guān)性，空間獨立的假設(shè)也是不能很好滿足的［1］。當(dāng)子帶寬或空間的獨立性不能滿足時，在完全發(fā)展相干斑的假設(shè)下，經(jīng)過多視處理所獲得名義上的視數(shù)要大于其有效的視數(shù)值［2］(Effective Number of Looks)。此外，SAR圖像的灰度范圍發(fā)生變化［3］、圖像的尺寸調(diào)整與裁剪、不同圖像格式之間轉(zhuǎn)換以及圖像經(jīng)過信道編碼傳輸?shù)榷紩ο喔砂叩膹姸犬a(chǎn)生影響。因此，對SAR圖像中相干斑的有效視數(shù)的估計顯得非常重要。

最早研究視數(shù)估計的學(xué)者是華人科學(xué)家Lee J.S.，在文獻(xiàn)［4］中，他提出一種有監(jiān)督的視數(shù)估計方法，但由于該方法涉及人工的操作使得估計結(jié)果存在比較大的不確定性。為了克服這個不足，Lee J.S.等［3］又提出了一種無監(jiān)督的視數(shù)估計方法，但該方法依然需要依靠經(jīng)驗人工地選擇較多個參數(shù)，如:扇形的角度、小區(qū)域的大小，并且參數(shù)值對最終的估計結(jié)果很敏感。Sathit I.等［5］也提出了一種無監(jiān)督的視數(shù)估計，但是該方法只能用于強度SAR圖像，并且對估計圖像廣義Gamma分布的假設(shè)使得估計非常復(fù)雜。

基于Lee J.S.的思想［3］，本文提出了一種新的無監(jiān)督相干斑強度估計方法。首先把SAR圖像劃分為多個圖像塊，接著對每個圖像塊計算其方差系數(shù)，最后統(tǒng)計已得到的方差系數(shù)的概率密度函數(shù)，并根據(jù)概率密度函數(shù)的最大值得到最終的相干斑強度估計。相比Lee的算法［3］，本文所提算法在參數(shù)選擇上更加容易，估計思想及過程非常簡便，大量模擬以及真實數(shù)據(jù)的測試證明了本文所提算法的有效性。

1 基于均值—標(biāo)準(zhǔn)差平面的相干斑強度估計

在已知SAR圖像視數(shù)信息的情況下，相干斑名義上的強度值可使用下面的式子計算獲得:

其中，Ln代表從數(shù)據(jù)源獲得的名義上的視數(shù)值。

這里先給出圖像方差系數(shù)γ(Variation Coefficient)的定義［4］:

其中，f代表某種信號。

若假設(shè)信號f等于相干斑Z，考慮到SAR圖像相干斑Z的均值為1，我們定義相干斑的方差系數(shù)如下:

可以看出，相干斑的方差系數(shù)實際上是和其標(biāo)準(zhǔn)差相等的。在一個同質(zhì)區(qū)域內(nèi)，真實的SAR圖像X是一個常數(shù)X0［4］，應(yīng)用SAR圖像乘性模型Y=X·Z(Y表示SAR圖像的觀測值)，我們可得到

(4)式說明，在一個同質(zhì)區(qū)域內(nèi)，其觀測圖像的方差系數(shù)與相干斑的標(biāo)準(zhǔn)差是相等的，換句話說，相干斑強度的估計可以轉(zhuǎn)化為同質(zhì)區(qū)域內(nèi)的方差系數(shù)的估計。

通常，為了計算γY，需要選擇許多同質(zhì)區(qū)域。接著計算每個區(qū)域的均值mY和標(biāo)準(zhǔn)差σY，并將每個塊的均值與標(biāo)準(zhǔn)差作為一個點(mY，σY)加入到一個二維平面上。實際中，這些點會聚成一類，找個一條經(jīng)過原點的直線，同時讓該直線通過聚類中心，此時直線的斜率就是所要計算的γY，這就是Lee J.S.等在文獻(xiàn)［4］中提到的估計方法，圖1給出了這種方法的一個估計示例。由于需要人工的選擇同質(zhì)區(qū)域，因此這是一種有監(jiān)督的估計方法。為了改進(jìn)人工選擇的缺陷，Lee J.S.等在文獻(xiàn)［3］中提出了一種估計方法，其基本思路是將圖像劃分為一個個互不重疊的小區(qū)域，然后采用兩種方法自動的估計。文獻(xiàn)［3］中的無監(jiān)督估計方法實際上是需要人工設(shè)定多個選擇參數(shù)的，并且每個參數(shù)都對估計結(jié)果敏感。

圖1 均值平方—方差平面估計相干斑強度示意圖

2 本文算法——基于方差系數(shù)統(tǒng)計的相干斑強度估計

針對Lee J.S.等［3］算法的不足，本文提出了一種基于方差系數(shù)統(tǒng)計的相干斑強度估計算法，具體的思路是:

先將整個SAR圖像劃分為一個個互不重疊的正方形小區(qū)域(假設(shè)尺寸為W×W)，再計算每個區(qū)域的均值和標(biāo)準(zhǔn)差并繪制均值平方—方差平面圖，最后使用一個估計方法獲得結(jié)果，其核心步驟是如何估計通過均值平方—方差平面中點數(shù)最多的直線的斜率，換種角度考慮就是出現(xiàn)頻率最多的方差系數(shù)。這個估計過程可從統(tǒng)計學(xué)角度求解，估計所有的方差系數(shù)的概率密度函數(shù)(Probability Density Function-PDF)，在PDF最大值出現(xiàn)的地方所對應(yīng)的方差系數(shù)值就是所要求的最優(yōu)估計:

其中，f(γY)代表γY的概率密度函數(shù)。

直方圖是一種應(yīng)用最為廣泛，并且簡單、有效的PDF估計方法［6］。假設(shè)x0為起始點，h為條形小區(qū)間的寬度，那么在這些小區(qū)間的直方圖可以定義在一系列區(qū)間［x0+mh，x0+(m+1)h］上，對于x點處的PDF可以估計為:

式中 m=0，1，．．．，M-1，其中，Nm表示區(qū)間［x0+mh，x0+(m+1)h］內(nèi)的樣本個數(shù)，并假設(shè)條形小區(qū)間的總數(shù)為M，樣本總數(shù)為n。

在方差系數(shù)的概率密度函數(shù)f(γY)的估計中，起始點γY0等于0(當(dāng)小區(qū)域的方差為0時)，那么區(qū)間就可以簡化為，那么此時的估計可以表示為:

其中，nr表示整幅SAR圖像劃分的小區(qū)域的個數(shù)。

本文所提出算法的具體執(zhí)行步驟是:

a.將一幅給定的SAR圖像劃分為數(shù)個互不重疊的W×W的圖像塊;

b.根據(jù)公式(2)計算劃分得到的所有圖像塊的方差系數(shù)γY;

c.根據(jù)公式(7)利用直方圖法計算這些方差系數(shù)的概率密度函數(shù);

3 實驗結(jié)果與分析

為了測試本文所提算法的性能，我們首先選擇了兩幅自然圖像:Cameraman(256×256，圖2(a))和Peppers(256×256，圖2(b))并分別加不同視數(shù)的幅度模擬相干斑。Peppers和Cameraman分別代表了包含較多同質(zhì)場景與異質(zhì)(Heterogeneous)的兩類圖像。在實驗中，所有圖像被劃分為4×4(W=4)的小塊，在方差系數(shù)的PDF的估計中，取條形小區(qū)域的寬度h=0．001。表1給出了在不同等級的模擬相干斑下，使用本文所提算法對相干斑強度的估計結(jié)果。可以看出，對不同等級和不同場景的圖像，所提算法都可以獲得很高的估計精度。對于異質(zhì)區(qū)域較多的Cameraman，其平均偏差為5.17%，相比較而言，同質(zhì)區(qū)域較多的Peppers，其平均估計偏差則只有2.03%，并且最低的偏差只有0.46%。此外，也可以發(fā)現(xiàn)隨著相干斑強度的降低，本文所提算法的估計精度逐漸提高;所提算法對同質(zhì)區(qū)域較多的Peppers在大多數(shù)情況下的估計精度要高于異質(zhì)區(qū)域較多的Cameraman。

圖2 不同等級和不同場景的模擬相干斑比較圖

一幅真實的單視幅度SAR圖像Howland Forest(401×386，圖2(c))被用來測試本文所提算法，以及對比所提算法與Lee J.S.等在文獻(xiàn)［3］中提出的扇形區(qū)間估計算法的性能。本文所提算法的對相干斑方差的估計為0.509，而Lee J.S.在文獻(xiàn)［3］中的估計結(jié)果為0.508，其理論值為0.5227。對于真實圖像而言，我們的估計結(jié)果與扇形區(qū)間估計法的結(jié)果非常接近。但是相比較而言，Lee J.S.的方法需要額外選擇扇形區(qū)間角度Δθ以及扇形區(qū)間重疊角度Od。由于一個小的Δθ可能會導(dǎo)致錯誤的結(jié)果，而一個大的Δθ則會使估計精度降低，因此Δθ的選擇需要更多的人工經(jīng)驗，而本文所提算法則是自動的不要人工選擇參數(shù)，實現(xiàn)了真正意義上的無監(jiān)督估計。

表1 本文所提算法在不同等級的模擬相干斑下的估計結(jié)果

4 結(jié)束語

本文提出一種基于方差系數(shù)統(tǒng)計無監(jiān)督的SAR相干斑強度估計算法，首先將整幅SAR圖像劃分為一個個互不重疊的小區(qū)域，接著計算每個小區(qū)域的方差系數(shù)。根據(jù)最優(yōu)的相干斑方差系數(shù)估計就是出現(xiàn)頻率最多的那個方差系數(shù)值這一基本思想，利用直方圖估計方差系數(shù)的概率密度函數(shù)以及區(qū)間掃描的方法求得概率密度函數(shù)的最大值和相應(yīng)的最優(yōu)相干斑方差估計值。相比于傳統(tǒng)估計方法，本文算法自動執(zhí)行、執(zhí)行過程簡單并且算法復(fù)雜度低，很適合工程應(yīng)用。

［1］Oliver C.，Quegan S.Understanding synthetic aperture radar images［M］.Boston，MA:Artech House，1998.

［2］Henri Maitre［法］編，孫洪等譯.合成孔徑雷達(dá)圖像處理［M］.北京:電子工業(yè)出版社，2005.

［3］Lee J.S.，Hoppel K.，and Mango S.A.Unsupervised estimation of speckle noise in radar images［J］.International Journal of Imaging Systems and Technology，1992，4:298-305.

［4］Lee J.S.，Hoppel K.Noise modeling and estimation of remotely-sensed images［C］.In Proc.IGARSS，1989，2:1005-1008.

［5］Sathit I.Speckle noise estimation with generalized Gamma distribution［C］.In Proceedings SICE-ICASE International Joint Conference，2006:1164-1167.

［6］Knuth D.The art of computer programming，volume 3:sorting and searching［M］.third ed.New Jersey:Addison-Wesley Publishing Company，1997:106-110.