一種考慮通道噪聲的極化SAR定標(biāo)改進(jìn)算法

張海瀛，李真芳，李錦偉，房 超

(西安電子科技大學(xué) 雷達(dá)信號(hào)處理國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710071)

?

一種考慮通道噪聲的極化SAR定標(biāo)改進(jìn)算法

張海瀛，李真芳，李錦偉，房 超

(西安電子科技大學(xué) 雷達(dá)信號(hào)處理國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710071)

在利用極化合成孔徑雷達(dá)數(shù)據(jù)進(jìn)行各項(xiàng)研究時(shí)，首先需要對(duì)極化合成孔徑雷達(dá)數(shù)據(jù)進(jìn)行標(biāo)定．現(xiàn)有Ainsworth算法在利用分布目標(biāo)估計(jì)交叉極化通道不平衡參數(shù)時(shí)，沒(méi)有考慮交叉極化通道噪聲的影響，這將影響該參數(shù)的估計(jì)穩(wěn)健性．針對(duì)此種情況，提出了一種改進(jìn)的Ainsworth定標(biāo)算法．該改進(jìn)算法充分考慮了交叉極化通道噪聲對(duì)交叉極化通道不平衡參數(shù)估計(jì)的影響，在交叉極化通道噪聲不一致以及低信噪比的情況下具有良好的性能．基于仿真數(shù)據(jù)和實(shí)測(cè)數(shù)據(jù)的實(shí)驗(yàn)結(jié)果驗(yàn)證了該改進(jìn)算法的有效性．

極化合成孔徑雷達(dá);極化定標(biāo);Ainsworth算法

極化合成孔徑雷達(dá)(Polarimetric Synthetic Aperture Radar，PolSAR)遙感已經(jīng)在地物分類以及植被參數(shù)估計(jì)等方面得到廣泛的研究和應(yīng)用[1-2]．在實(shí)際應(yīng)用中，所有的極化分析與應(yīng)用均需建立在各極化通道間的相對(duì)幅度和相對(duì)相位已進(jìn)行精確標(biāo)定的前提下，因而在利用極化合成孔徑雷達(dá)(Synthetic Aperture Radar，SAR)數(shù)據(jù)進(jìn)行各應(yīng)用研究時(shí)，首要，需要校正由系統(tǒng)性能不理想引起的極化數(shù)據(jù)的畸變．目前，國(guó)內(nèi)外眾多學(xué)者已提出了大量的極化定標(biāo)算法及其改進(jìn)算法，主要包括使用點(diǎn)目標(biāo)、使用分布目標(biāo)[3]及混合使用點(diǎn)目標(biāo)和分布目標(biāo)[4-8]等3類定標(biāo)算法，其中，混合使用分布目標(biāo)和點(diǎn)目標(biāo)的定標(biāo)算法對(duì)極化系統(tǒng)和分布目標(biāo)進(jìn)行了一定的假設(shè)，簡(jiǎn)化了定標(biāo)過(guò)程，因而該類算法在實(shí)際中得到了廣泛應(yīng)用．

由于目前的極化SAR系統(tǒng)可獲得很高的極化隔離度，并能通過(guò)精確的內(nèi)定標(biāo)技術(shù)進(jìn)一步提高極化數(shù)據(jù)的精度，因而，文獻(xiàn)[6]提出了一種假設(shè)條件更加寬泛的定標(biāo)算法——Ainsworth算法，其僅要求目標(biāo)滿足散射互易性，從而大大擴(kuò)展了可用定標(biāo)區(qū)域的選取范圍，其已被成功應(yīng)用于F-SAR系統(tǒng)[9]和無(wú)人機(jī)合成孔徑雷達(dá)(Unmanned Aerial Vehicle Synthetic Aperture Radar，UAVSAR)系統(tǒng)[10]極化數(shù)據(jù)的標(biāo)定．在實(shí)際應(yīng)用中，現(xiàn)有Ainsworth算法在估計(jì)不平衡參數(shù)α的過(guò)程中并沒(méi)有考慮交叉極化通道噪聲的影響，因而當(dāng)交叉極化通道噪聲不可忽略時(shí)，參數(shù)α估計(jì)精度將會(huì)下降，從而導(dǎo)致極化數(shù)據(jù)的畸變無(wú)法被正確校正．為此，筆者提出了一種考慮交叉極化通道噪聲的極化SAR定標(biāo)改進(jìn)算法，可消除交叉極化通道噪聲對(duì)參數(shù)α估計(jì)的不良影響，提高參數(shù)α的估計(jì)精度．

1 極化定標(biāo)模型

PolSAR定標(biāo)算法的有效性與PolSAR系統(tǒng)模型的有效性緊密相關(guān)．觀測(cè)的極化散射矩陣O與目標(biāo)的真實(shí)散射矩陣S間的關(guān)系，可用PolSAR系統(tǒng)誤差模型進(jìn)行描述[5]，即

其中，R為接收失真矩陣，T為發(fā)射失真矩陣，N為系統(tǒng)加性噪聲．式(1)中的矩陣均為 2×2 矩陣．式(1)可詳細(xì)描述為

其中，下標(biāo)H和V分別表示水平和垂直極化狀態(tài)；SHV表示發(fā)射V極化、接收H極化時(shí)的回波，其余元素具有相似的含義．假設(shè)各極化通道噪聲與信號(hào)不相關(guān)且各極化通道噪聲間也不相關(guān)，則觀測(cè)協(xié)方差矩陣Co和真實(shí)協(xié)方差矩陣Cs間的關(guān)系可表示為

其中，Co=o oH，Cs=s sH，Cn=n nH，[·]H表示共軛轉(zhuǎn)置操作，·表示求集合平均，o為觀測(cè)極化散射矢量，s為真實(shí)散射矢量，n為系統(tǒng)噪聲矢量．各變量定義如下:

其中，k為接收通道不平衡參數(shù)，α為交叉極化通道不平衡參數(shù)，u、v、w和z為極化通道串?dāng)_參數(shù)，[·]T表示轉(zhuǎn)置操作．通常情況下，參數(shù)Y和k的求取與參數(shù)u、v、w、z、α的求取可獨(dú)立進(jìn)行．Y和k參數(shù)的求取需要散射矩陣精確已知的先驗(yàn)?zāi)繕?biāo)，詳細(xì)的描述可參考文獻(xiàn)[4，8]．

2 改進(jìn)的Ainsworth定標(biāo)算法

Ainsworth極化定標(biāo)方法僅僅要求場(chǎng)景中地物滿足互易性，即SHV=SVH，并不要求場(chǎng)景中地物滿足同極化通道和交叉極化通道間不相關(guān)的假設(shè)，即=0，i、j=h 或v，i≠j，從而擴(kuò)大了定標(biāo)區(qū)域的選取范圍．由式(3)定義可得

當(dāng)場(chǎng)景中的地物滿足互易性時(shí)，矩陣D具有以下形式:

在實(shí)際應(yīng)用中，Ainsworth算法通常忽略噪聲項(xiàng)Cn，利用

α=(C33C22)1/2

實(shí)現(xiàn)參數(shù)α的估計(jì)，其中，Cxy表示矩陣C的xy項(xiàng)，下文中所有類似符號(hào)的含義與此處相同，arg(·)表示取相角操作．式(8)并沒(méi)有考慮交叉通道噪聲對(duì)參數(shù)α估計(jì)的影響．當(dāng)交叉極化通道噪聲的影響不可忽略，但在實(shí)際處理中不予考慮時(shí)，參數(shù)α的估計(jì)精度將受到影響．文中加性噪聲被建模為零均值高斯白噪聲[5]，即

在實(shí)際中，各極化通道的噪聲功率可能并不相等，尤其對(duì)于H和V極化接收通道分置的系統(tǒng)，如相控陣型L波段合成孔徑雷達(dá)(Phased Array type L-band Synthetic Aperture Radar，PALSAR)系統(tǒng)等，因而文中假定VH通道和HV通道的噪聲功率比 m= N3/ N2．m值可利用符合一階Bragg散射模型的光滑表面進(jìn)行估計(jì)，一般平靜的水面是最佳目標(biāo)．為消除交叉極化通道噪聲對(duì)參數(shù)α估計(jì)精度的影響，文中將Anisworth算法基本步驟中參數(shù)α的估計(jì)式(8)修改為[12]

式(10)在估計(jì)α?xí)r，不僅考慮了交叉極化通道存在噪聲時(shí)的情況，而且還考慮了交叉極化通道噪聲不一致的情況．改進(jìn)的Anisworth定標(biāo)算法的基本處理流程如下所述:

(1) 根據(jù)式(10)估計(jì)參數(shù)α的值．

(2) 計(jì)算矩陣A=diag(α，1，α，1)，并利用式(7)計(jì)算矩陣L．

(3) 利用式(11)估計(jì)G和H:

(4) 利用式(12)～(13)估計(jì)串?dāng)_參數(shù)u、v、w、z，即

其中，κ=[u，v，z，w]T，且有

(5) 利用步驟(1)中估計(jì)得到的α值、步驟(4)中估計(jì)得到的串?dāng)_參數(shù)及式(7)，完成矩陣D的估計(jì)．

因此，在忽略參數(shù)u、v、w、z的高階項(xiàng)(二階及以上)的情況下，可得轉(zhuǎn)換后的噪聲矩陣P中的元素為

因而在此步驟中需要按照

更新m值．

利用上述步驟可進(jìn)行循環(huán)迭代，即將每次迭代估計(jì)得到的矩陣D作為下一次迭代的輸入．如此循環(huán)迭代，直到迭代次數(shù)達(dá)到設(shè)定的次數(shù)或迭代過(guò)程中所有串?dāng)_參數(shù)的中間值低于設(shè)定閾值．根據(jù)文獻(xiàn)[6，11]，在實(shí)驗(yàn)中，最大迭代次數(shù)設(shè)定為12，閾值設(shè)定為10-8．最終的參數(shù)估計(jì)值可由每次迭代的中間值(ui，vi，wi，zi，αi)給出，即

3 實(shí)驗(yàn)結(jié)果分析

為驗(yàn)證文中改進(jìn)定標(biāo)算法在精確估計(jì)參數(shù)α方面的有效性，分別利用仿真數(shù)據(jù)和德國(guó)Oberpfaffenhofen地區(qū)全極化數(shù)據(jù)對(duì)改進(jìn)算法進(jìn)行了實(shí)驗(yàn)分析．由于實(shí)際數(shù)據(jù)的失真矩陣不可能事先知道，因此，利用已標(biāo)定好的全極化數(shù)據(jù)，通過(guò)設(shè)定失真參數(shù)值來(lái)產(chǎn)生實(shí)驗(yàn)數(shù)據(jù)．

真實(shí)協(xié)方差矩陣Cs中的元素被設(shè)置為代表植被的典型值[11]．在文中，主要分析交叉極化通道噪聲對(duì)參數(shù)α的影響，以及改進(jìn)算法對(duì)于α估計(jì)精度的提升性能．由于目前的極化雷達(dá)系統(tǒng)大都具有較高的極化隔離度，這意味著串?dāng)_參數(shù)均較小(幅度值小于0.1)，所以串?dāng)_參數(shù)對(duì)參數(shù)α估計(jì)精度的影響可忽略．文中將串?dāng)_參數(shù)u、v、w和z設(shè)置為固定值，如表1所示．參數(shù)α的幅度和相位分別在0.9～1.1和0～2π范圍內(nèi)選取．構(gòu)造噪聲矩陣Cn，使其具有diag(N，N/m，N，N)的形式，其中，N的取值范圍為0.00～0.05，m的取值范圍為0.5～2.0．

表1 仿真時(shí)串?dāng)_參數(shù)設(shè)置值

為分析改進(jìn)算法估計(jì)參數(shù)α的性能隨m值的變化情況，設(shè)定3個(gè)具有代表性的α值，其幅度值分別為0.9、1.0和1.1，相位均為 -70.08°，并設(shè)定 N= 0.05．圖1表示的是m取不同值時(shí)參數(shù)α的估計(jì)值及估計(jì)誤差．圖1(a)為利用原始的Ainsworth算法估計(jì)得到的α幅度值，圖1(d)為利用改進(jìn)算法估計(jì)得到的幅度值．從圖1(a)和圖1(d)可以看出，原始Ainsworth算法的α幅度值的估計(jì)精度受交叉極化通道噪聲影響明顯，而改進(jìn)后的算法對(duì)噪聲表現(xiàn)了良好的魯棒性，α幅度值的估計(jì)精度大大提高．從圖1(c)可看出，當(dāng)交叉通道噪聲功率不一致時(shí)，原始算法對(duì)于α值的估計(jì)將會(huì)產(chǎn)生較大的誤差，達(dá)到10-1數(shù)量級(jí)，而由圖1(f)可看出，利用改進(jìn)算法可大大提高α的估計(jì)精度，誤差下降到10-3數(shù)量級(jí)．另外，從圖1(b)和圖1(e)可以看出，原始Ainsworth算法和改進(jìn)算法對(duì)于α相位值的估計(jì)精度差異不大，誤差均在 ±1° 之內(nèi)．由于原始Ainsworth算法和改進(jìn)算法對(duì)于α相位值的估計(jì)準(zhǔn)確度均很高，因此，下面將主要在α的幅度值估計(jì)精度方面對(duì)改進(jìn)定標(biāo)算法的估計(jì)性能展開(kāi)討論．

圖1 基于原始Ainsworth算法和改進(jìn)算法的α估計(jì)值比較

為分析交叉極化通道信噪比對(duì)參數(shù)α估計(jì)精度的影響，分別將α和m的值設(shè)置為固定值1.1和0.5，而N的值由0.00到0.05均勻變化．文中分別利用原始Ainsworth算法和改進(jìn)算法對(duì)參數(shù)α值進(jìn)行估計(jì)．由于在上述設(shè)置下，HV通道的信噪比要低于VH通道的信噪比(參考噪聲矩陣的形式)，因而僅畫(huà)出α估計(jì)值以及其估計(jì)誤差隨HV通道信噪比的變化趨勢(shì)，如圖2所示．圖2(a)為利用原始Ainsworth算法和改進(jìn)算法估計(jì)得到的α幅度值隨信噪比的變化，圖2(b)為其估計(jì)誤差隨信噪比的變化．從圖2可看出，當(dāng)交叉極化通道信噪比較高時(shí)(此時(shí)對(duì)應(yīng)低噪聲水平)，原始Ainsworth算法可獲得良好的α估計(jì)性能，但當(dāng)信噪比降低時(shí)，估計(jì)精度明顯降低．而改進(jìn)算法的估計(jì)性能對(duì)信噪比的變化表現(xiàn)了良好的魯棒性，低信噪比下估計(jì)誤差仍然很小，因而改進(jìn)算法對(duì)于參數(shù)α的估計(jì)性能要明顯優(yōu)于原始算法．

圖2 α幅度估計(jì)值及其估計(jì)誤差隨信噪比的變化

為進(jìn)一步驗(yàn)證改進(jìn)定標(biāo)算法的性能，利用德國(guó)Oberpfaffenhofen地區(qū)全極化數(shù)據(jù)，通過(guò)設(shè)定失真參數(shù)值和加性噪聲值產(chǎn)生了一組實(shí)驗(yàn)數(shù)據(jù)．該數(shù)據(jù)中，串?dāng)_參數(shù)設(shè)置值仍然采用表1中所示的值，α值的幅度和相位值分別為1.1和 -70.08°，VH通道和HV通道的所加噪聲的功率之比m取值為0.14．圖3為基于該數(shù)據(jù)利用原始Ainsworth算法和改進(jìn)算法估計(jì)得到的參數(shù)α的幅度和相位沿距離向的變化．從圖3可看出，基于原始Ainsworth算法和改進(jìn)算法的α的幅度和相位估計(jì)值隨距離向均有所波動(dòng)，這是由于實(shí)測(cè)數(shù)據(jù)中噪聲沿距離向不均勻造成的．從圖3(b)可看出，原始算法和改進(jìn)算法在α相位值估計(jì)方面的性能基本相同，這與上面的分析結(jié)果相同．但值得注意的是，基于改進(jìn)算法的α幅度估計(jì)值隨距離向的波動(dòng)明顯比基于原始Ainsworth算法的小得多，如圖3(a)所示，這說(shuō)明改進(jìn)算法的穩(wěn)健性要明顯優(yōu)于原始算法．

圖3 基于德國(guó)Oberpfaffenhofen地區(qū)全極化數(shù)據(jù)的α估計(jì)值隨距離的變化

為更直觀地說(shuō)明改進(jìn)算法的性能提升，這里給出了德國(guó)Oberpfaffenhofen地區(qū)全極化數(shù)據(jù)無(wú)失真時(shí)的合成圖以及分別基于原始算法和改進(jìn)算法失真參數(shù)估計(jì)值校正后的合成圖，如圖4所示．從圖4中可看出，相比于原始算法，基于改進(jìn)算法估計(jì)得到的失真參數(shù)校正后的合成圖更接近于無(wú)失真時(shí)的合成圖．

圖4 德國(guó)Oberpfaffenhofen地區(qū)全極化數(shù)據(jù)合成圖

4 結(jié) 束 語(yǔ)

現(xiàn)有Ainsworth算法在利用分布目標(biāo)估計(jì)參數(shù)α?xí)r沒(méi)有考慮交叉極化通道噪聲的影響，當(dāng)交叉極化通道噪聲不可忽略時(shí)，參數(shù)α的估計(jì)精度將會(huì)明顯降低．針對(duì)此種情況，文中提出了一種考慮交叉極化通道噪聲的極化SAR定標(biāo)改進(jìn)算法．實(shí)驗(yàn)結(jié)果表明，該改進(jìn)算法在交叉通道噪聲功率不一致以及信噪比較低的情況下仍能保持良好的估計(jì)性能，從而可有效保證參數(shù)α的估計(jì)精度，實(shí)現(xiàn)極化SAR數(shù)據(jù)的精確標(biāo)定，利于后續(xù)的極化分類等應(yīng)用的研究．

[1] 劉璐, 王蓉芳, 焦李成, 等. 采用聯(lián)合相似度譜聚類集成的極化SAR 分類[J]. 西安電子科技大學(xué)學(xué)報(bào), 2015, 42(3): 48-53.

LIU Lu, WANG Rongfang, JIAO Licheng, et al. Combined Similarity Based Spectral Clustering Ensemble for POLSAR Classification[J]. Journal of Xidian University, 2015, 42(3): 48-53.

[2]盧紅喜, 李東霖, 劉宏偉, 等. 雙基線極化干涉SAR 植被參數(shù)反演[J]. 西安電子科技大學(xué)學(xué)報(bào), 2015, 42(6): 23-29.

LU Hongxi, LI Donglin, LIU Hongwei, et al. Forest Parameters Retrieval with Dual-baseline Polarimetric SAR Interferometry Based on Clustering Analysis[J]. Journal of Xidian University, 2015, 42(6): 23-29.

[3]VILLA A, IANNINI L, GIUDICI D, et al. Calibration of SAR Polarimetric Images by Means of a Covariance Matching Approach[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(2): 674-686.

[4]VAN ZYL J J. Calibration of Polarimetric Radar Images Using Only Image Parameters and Trihedral Corner Reflector Responses[J]. IEEE Transactions on Geoscience and Remote Sensing, 1990, 28(3): 337-348.

[5]QUEGAN S. A Unified Algorithm for Phase and Cross-talk Calibration of Polarimetric Data-theory and Observations[J]. IEEE Transactions on Geoscience and Remote Sensing, 1994, 32(1): 89-99.

[6]AINSWORTH T L, FERRO-FAMIL L, LEE J S. Orientation Angle Preserving a Posteriori Polarimetric SAR Calibration[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(4): 994-1003.

[7]ZHANG H, LU W, ZHANG B, et al. Improvement of Polarimetric SAR Calibration Based on the Ainsworth Algorithm for Chinese Airborne PolSAR Data[J]. IEEE Geoscience and Remote Sensing Letters, 2013, 10(4): 898-902.

[8]SHI L, YANG J, LI P. Co-polarization Channel Imbalance Determination by the Use of Bare Soil[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2014, 95: 53-67.

[9]FISCHER J, BAUMGARTNER S V, REIGBER A, et al. Geometric, Radiometric, Polarimetric and Along-track Interferometric Calibration of the New F-SAR System of DLR in X-band [C]//2008 7th European Conference on Synthetic Aperture Radar. Frankfurt: VDE, 2008: 1-4.

[10]FORE A G, CHAPMAN B D, HAWKINS B P, et al. UAVSAR Polarimetric Calibration[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(6): 3481-3491.

[11]GOH A S, PREISS M, Gray D A, et al. Comparison of Parameter Estimation Accuracy of Distributed-target Polarimetric Calibration Techniques[C]//International Geoscience and Remote Sensing Symposium. Piscataway: IEEE, 2007: 4175-4178.

[12]KIMURA H, MIZUNO T, PAPATHANASSIOU K P, et al. Improvement of Polarimetric SAR Calibration Based on the Quegan Algorithm[C]//International Geoscience and Remote Sensing Symposium. Piscataway: IEEE, 2004: 184-187.

(編輯：齊淑娟)

Improved polarimetric SAR calibration algorithm considering channel noise

ZHANGHaiying，LIZhenfang，LIJinwei，F(xiàn)ANGChao

(National Key Lab. of Radar Signal Processing, Xidian Univ., Xi’an 710071, China)

It is necessary to calibrate polarimetric SAR(Synthetic Aperture Radar) data in order to use the data for scientific applications. The present Ainsworth algorithm does not consider the cross-polarization channel noise in estimating the cross-polarization imbalance parameter and therefore the estimation of this parameter is not robust. To solve the above problem, an improved Ainsworth algorithm is presented, which takes the cross-polarization channel noise into consideration. This algorithm can preserve a high performance under the situations of difference cross-polarization channel noise powers and a low signal noise ratio. Experimental results based on simulated and real data confirm the effectiveness of the improved algorithm.Key Words: polarimetric SAR;polarimetric calibration; Ainsworth algorithm

2015-10-25

時(shí)間：2016-04-01

國(guó)家自然科學(xué)基金資助項(xiàng)目(61471276，41371439)

張海瀛(1990-)，男，西安電子科技大學(xué)博士研究生，E-mail：zhanghaiying4155@163.com．

http://www.cnki.net/kcms/detail/61.1076.tn.20160401.1622.016.html

10.3969/j.issn.1001-2400.2016.06.008

TN957

A

1001-2400(2016)06-0045-06