洪 文楊士林*③李 洋③張晶晶③郭勝龍③

①(微波成像技術(shù)國家級重點實驗室 北京 100190)

②(中國科學院電子學研究所 北京 100190)

③(中國科學院大學 北京 100190)

分布式目標的極化SAR距離模糊計算方法研究

洪 文①②楊士林*①②③李 洋①②③張晶晶①②③郭勝龍①②③

①(微波成像技術(shù)國家級重點實驗室 北京 100190)

②(中國科學院電子學研究所 北京 100190)

③(中國科學院大學 北京 100190)

該文提出了一種運用協(xié)方差矩陣計算星載全極化SAR距離模糊的方法。首先分析了極化SAR距離模糊的成因及同極化和交叉極化通道間距離模糊能量的差異?，F(xiàn)有的距離模糊計算方法以散射矩陣為主，受相干斑影響，單一像素的散射強度不能表征分布目標的反射率，因而在距離向相鄰的像素無法確定其距離模糊比。針對上述問題該文推導了基于分布目標的極化SAR距離模糊的計算方法。結(jié)合Radarsat-2實測數(shù)據(jù)對這種新方法進行了驗證。實驗表明，運用協(xié)方差矩陣能夠很好地表征分布式目標，由它計算得出的距離模糊比是平穩(wěn)、確定的。相比已有算法，所提方法能夠有效計算極化SAR分布式目標的距離模糊。

合成孔徑雷達；分布目標；極化SAR；距離模糊

1 引言

合成孔徑雷達(SAR)是一種全天時、全天候的高分辨率成像雷達。雷達天線主瓣照射的區(qū)域為測繪帶。距離模糊是由于天線主瓣以外的旁瓣回波信號與所需測繪帶內(nèi)信號回波同時到達天線所造成的。距離模糊對于寬測繪帶SAR系統(tǒng)的參數(shù)設(shè)計起著非常重要的作用，它是影響SAR成像的一個重要因素[1]。目前單極化SAR系統(tǒng)的距離模糊分析已有很多詳細的討論[2?4]，而多極化SAR距離模糊的計算研究相對較少，該文主要分析了星載全極化SAR系統(tǒng)的距離模糊及其詳細計算方法。目前的極化SAR系統(tǒng)距離模糊計算中[5?7]，用表征地面散射特性的極化散射矩陣來直接計算距離模糊并不妥當。這是因為：星載SAR目標是分布式目標，它們隨環(huán)境和時間而動態(tài)變化。受相干斑影響[8]，相鄰像素間的信號強度是不連續(xù)的，如果從散射矩陣出發(fā)分析距離模糊比，那么計算出的交叉極化功率仍在擾動，就無法判斷同極化和交叉極化功率的相對大小，從而出現(xiàn)相鄰像素距離模糊比不能確定的問題。另外，交叉極化能量較同極化能量低6～10 dB[9]，但不能直接體現(xiàn)在散射矩陣上。對于星載SAR目標來說，極化數(shù)據(jù)的二階統(tǒng)計量協(xié)方差矩陣是大量分辨單元平均后的結(jié)果，它是穩(wěn)定的。用協(xié)方差矩陣來參與計算距離模糊是合理的，并且可以表征交叉極化能量和同極化能量的關(guān)系。該文首先分析了極化SAR距離模糊的形成原因，并且分析了同極化和交叉極化通道模糊能量的差異；詳細推導了基于分布式目標的極化SAR距離模糊的計算方法。接著結(jié)合數(shù)據(jù)對這種新方法進行了仿真驗證。最后與現(xiàn)有的計算距離模糊的方法進行了對比分析。

2 極化SAR距離模糊分析與計算

2.1 距離模糊分析

極化SAR系統(tǒng)的工作模式為交替發(fā)射水平(H)、垂直(V)極化波。在發(fā)射脈沖間隔內(nèi)同時接收發(fā)射脈沖對應的H, V散射回波。如圖1所示，發(fā)射H極化波時，水平通道接收HH散射波，垂直通道接收HV散射波；發(fā)射V極化波時，水平通道接收VH 散射波，垂直通道接收VV散射波。工作時序如圖1所示。

圖1 極化SAR工作時序圖(交叉極化回波比同極化回波低6～10 dB)

在距離向，假設(shè)某一區(qū)域的回波延時和觀測帶內(nèi)目標的回波延時正好相差整數(shù)倍個脈沖周期，這就形成了距離模糊，產(chǎn)生模糊信號的那個區(qū)域稱為距離模糊區(qū)，或者與觀測處的斜距相差整數(shù)個脈沖數(shù)的區(qū)域即為模糊區(qū)[10?12]。對于這種交替發(fā)射，同時接收的極化SAR工作模式，某個回波接收時刻，接收信號除了感興趣的信號之外，也有正好相差整數(shù)j個脈沖重復周期(IPP)的模糊回波信號。當j為奇數(shù)時，模糊回波信號對應的發(fā)射脈沖和有用信號對應的發(fā)射脈沖極化方式相反。我們把對應的模糊區(qū)為稱奇數(shù)模糊區(qū)(odd ambiguity zones)。當j為偶數(shù)時，模糊回波信號對應的發(fā)射脈沖和有用信號對應的發(fā)射脈沖極化方式相同。我們把對應的模糊區(qū)稱作偶數(shù)模糊區(qū)(even ambiguity zones)。圖2所示是極化SAR距離模糊示意圖。與t=0時刻相差Tp(脈沖重復周期)，?Tp的發(fā)射脈沖與t=0時刻的發(fā)射脈沖極化方式相反，對應的模糊區(qū)域為奇數(shù)模糊區(qū)。與t=0時刻相差2Tp,?2Tp的發(fā)射脈沖與t=0時刻的發(fā)射脈沖極化方式相同，對應的模糊區(qū)域為偶數(shù)模糊區(qū)。

基于以上分析，由文獻[9]可知，極化SAR觀測散射矩陣為

式中，RAR表示距離模糊比。等號右邊第1項表示測繪帶(swath)內(nèi)感興趣目標的散射矩陣，剩余兩項為距離模糊區(qū)域散射矩陣，j=od 為奇數(shù)模糊區(qū)域，j=ev為偶數(shù)模糊區(qū)域。

圖2 極化SAR距離模糊示意圖

由式(1)看出，有用信號同極化通道模糊能量分別來自奇數(shù)模糊區(qū)交叉極化通道和偶數(shù)模糊區(qū)同極化通道；有用信號交叉極化通道模糊能量分別來自奇數(shù)模糊區(qū)同極化通道和偶數(shù)模糊區(qū)交叉極化通道。由于同極化通道能量比交叉極化通道能量高6～10 dB，因此交叉極化通道受來自奇數(shù)模糊區(qū)同極化模糊能量的主導，模糊比值比較大[13]。

2.2 分布式目標距離模糊計算推導

下面推導距離模糊的計算公式。雷達信號接收電壓表達式為

其中?=?2kR ,R是散射目標到天線的斜距， η0=μ0和ε0分別表示磁導率和電導率，η是本地入射角，Pt是天線發(fā)射功率，G表示天線增益，λ是波長。

那么根據(jù)式(1)可以得到極化SAR的接收電壓為

其中Ssig,Sev,Sod分別表示有用信號，偶數(shù)模糊區(qū)，奇數(shù)模糊區(qū)的散射矩陣。

由于星載SAR目標是隨環(huán)境動態(tài)變化的，因此可以假設(shè)有用信號與模糊信號不相關(guān)，模糊信號之間也不相關(guān)(上標H表示共軛轉(zhuǎn)置)，即

那么接收信號功率可以表示為

其中，C=SSH是協(xié)方差矩陣。j是模糊區(qū)標號，其中j=0表示測繪帶內(nèi)感興趣的脈沖。i是測繪帶內(nèi)采樣點標號，Ri0表示觀測帶內(nèi)的第i個采樣點對應的斜距。那么Rij表示給i點造成模糊的所有點的斜距值。由式(6)就可以計算得出整數(shù)j的值。

Rn和Rf分別表示測繪帶內(nèi)的近點斜距和遠點斜距值；Re和h分別表示地球半徑和軌道高度。圖3表示斜距幾何關(guān)系。

圖3中，SAR到地球中心的距離定義為Rs=Re+h,γij為對應于入射角ηij的天線視軸角。

Gij表示在給定的γij下的天線方向圖。對于均勻照明孔徑，遠場方向圖由式(9)定義：

其中離軸高度角計算式為

式中，γ0表示天線下視角(look angle), r表示橫滾角。

根據(jù)距離模糊比(RASR)的定義[14]，距離模糊比是模糊信號能量與有用信號能量的比值，即

圖3 斜距幾何關(guān)系圖

根據(jù)式(1)和式(11)，傳統(tǒng)極化HH通道距離模糊比為

同理，傳統(tǒng)極化HV通道距離模糊比為

其中Cmn表示C協(xié)方差矩陣第m行，第n列的元素。由式(1)可知，極化SAR系統(tǒng)中C項是同極化通道，C和C是模糊能量中給同極化通道造成模糊的項。

至此，已詳細給出由全極化協(xié)方差矩陣C推導出距離模糊的計算表達式。下面用Radarsat-2實測數(shù)據(jù)進行驗證與分析。

3 實驗結(jié)果

3.1 實驗數(shù)據(jù)介紹

實驗所用實測數(shù)據(jù)來自于Radarsat-2系統(tǒng)，系統(tǒng)主要參數(shù)如表1所示。

文章分別對舊金山地區(qū)分布式目標海洋和森林運用本文的方法進行了計算。

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

根據(jù)文獻[6]現(xiàn)有的計算全極化SAR距離模糊的方法，海洋和森林的距離模糊計算結(jié)果如圖4～圖7所示。

由圖4～圖7可以看出，用散射矩陣描述分布式觀測對象，無論是海洋還是森林，每條距離模糊比曲線變化很不穩(wěn)定，曲線中相鄰兩點的距離模糊比跳變很大，而且如果存在強散射點，距離模糊比甚至有超過0 dB的可能。這是由于在距離模糊計算中，受相干斑因素的影響，相鄰像素的信號強度是不連續(xù)的，單一像素的極化散射矩陣不能表征分布式目標的散射率，由散射矩陣計算得出極化功率時刻在隨時間和空間而擾動[15]。極化數(shù)據(jù)的二階統(tǒng)計量協(xié)方差矩陣是穩(wěn)定的，是大量分辨單元平均后的結(jié)果，因此本文的運用散射數(shù)據(jù)的二階量協(xié)方差矩陣計算距離模糊是合理的。圖8～圖11結(jié)果是依據(jù)本文方法對海洋和森林目標計算得到的。

表1 Radarsat-2系統(tǒng)參數(shù)

圖8～圖11是采用本文的方法計算距離模糊比。由圖可以看出距離模糊比曲線比較平緩光滑，計算結(jié)果是比較穩(wěn)定的，所以可以根據(jù)當前某點的距離模糊比值而推理得到與它相鄰點的距離模糊比值以及變化趨勢。

此外，圖8所示是同極化通道的距離模糊比值，圖中測繪帶內(nèi)大部分區(qū)域距離模糊比都低于-20 dB(一般要求距離模糊比低于-20 dB)。圖9所示的交叉極化通道低于-20 dB的測繪帶寬不到20 km。這是因為由式(1)可知HV, VH有用回波信號的模糊能量受同極化信號的主導，因為同極化回波能量要比交叉極化高6～10 dB，所以交叉極化通道的距離模糊嚴重一些。HH, VV同極化波在奇數(shù)模糊區(qū)受交叉極化波的影響，因此它們所受的影響較小。偶數(shù)模糊區(qū)域和單極化SAR系統(tǒng)的模糊一樣。

圖4 海洋同極化通道距離模糊比

圖5 海洋交叉極化通道距離模糊比

圖6 森林同極化通道距離模糊比

圖7 森林交叉極化通道距離模糊比

圖8 海洋同極化通道距離模糊比

圖9 海洋交叉極化通道距離模糊比

圖10 森林同極化通道距離模糊比

圖11 森林交叉極化通道距離模糊比

4 結(jié)束語

本文詳細地分析了極化SAR系統(tǒng)距離模糊形成的原因，針對現(xiàn)有的星載極化SAR距離模糊計算方法中用散射矩陣無法描述分布式觀測對象極化散射特征的問題，本文提出了一種運用協(xié)方差矩陣計算星載全極化SAR距離模糊的方法。詳細推導出極化SAR系統(tǒng)HH, HV通道距離模糊的計算公式。最后，通過Radarsat-2舊金山地區(qū)海洋和森林分布式目標距離模糊計算，驗證了運用協(xié)方差矩陣計算極化SAR系統(tǒng)距離模糊的合理性與精確性，為極化SAR距離模糊分析與計算提供一種有價值的參考。

[1] 韓曉東, 宋紅軍, 徐偉. 基于星載SAR實時波束形成的星下點及距離模糊抑制方法[J]. 電子與信息學報, 2013, 35(12): 2821-2828.

Han Xiao-dong, Song Hong-jun, and Xu Wei. Nadir and range ambiguity suppression based on real-time digital beam froming in spaceborne SAR[J]. Journal of Electronics & Information Technology, 2013, 35(12): 2821-2828.

[2] Wang Wen-qin. Mitigating range ambiguities in high-PRF SAR with OFDM waveform diversity[J]. IEEE Geoscience and Remote Sensing Letters, 2013, 10(1): 101-105.

[3] Yang Jun, Sun Guang-cai, Wu Yu-feng, et al.. Range ambiguity suppression by azimuth phase coding in multichannel SAR systems[C]. IET International Radar Conference 2013, Xi'an, China, 2013: 1-5.

[4] Davis Michael S and Cook Daniel A. Quantifying the impact of range and azimuth ambiguities on multichannel SAR and SAS image quality[C]. EUSAR 2014: 10th European Conference on Synthetic Aperture Radar, Berlin, Germany, 2014: 1-4.

[5] Ralph C and Pierluigi M. Range ambiguities for a polarimetric spaceborne SAR[C]. IGARSS 92: Proceedings of the 12th Annual International Geoscience and Remote Sensing Symposium, Houston, TX, US, 1992: 637-639.

[6] 李延, 丁澤剛, 郭佳佳. 多極化星載SAR模糊比分析與計算[J]. 航天器工程, 2011, 20(1): 45-49.

Li Yan, Ding Ze-gang, and Guo Jia-jia. Ambiguity analysis and calculation of polarimetric SAR[J]. Spacecraft Engineering, 2011, 20(1): 45-49.

[7] 戴博偉, 楊汝良. 多極化星載合成孔徑雷達系統(tǒng)新工作方式及模糊分析[J]. 電子與信息學報, 2001, 23(12): 1256-1263.

Dai Bo-wei and Yang Ru-liang. New operation modes and ambiguity analysis for spaceborne polarimetric SAR system [J]. Journal of Electronics & Information Technology, 2001, 23(12): 1256-1263.

[8] Foucher S and Lopez-Martinez C. Analysis, evaluation, and comparison of polarimetric SAR speckle filtering techniques [J]. IEEE Transactions on Image Processing, 2014, 23(4): 1751-1764.

[9] Raney R K, Freeman A, and Jordan R L. Improved range ambiguity performance in quad-pol SAR[J]. IEEETransactions on Geoscience and Remote Sensing, 2012, 50(2): 349-356.

[10] 劉軍鷹, 蔡竟業(yè), 王文欽. 星載距離向多波束SAR的系統(tǒng)模糊特性[J]. 信息與電子工程, 2012, 10(2): 138-142.

Liu Jun-ying, Cai Jing-ye, and Wang Wen-qin. Ambiguity characteristics of spaceborne range multi-beam SAR systems[J]. Information and Electronic Engineering, 2012, 10(2): 138-142.

[11] 席龍梅, 盛磊, 葛家龍. 星載滑動聚束SAR模糊特性分析與仿真[J]. 上海航天, 2011, 28(1): 1-6.

Xi Long-mei, Sheng Lei, and Ge Jia-long. Analysis and simulation of ambiguities on spaceborne sliding spotlight SAR[J]. Aerospace Shanghai, 2011, 28(1): 1-6.

[12] 荊麟角. 星載SAR距離模糊分布規(guī)律及其改進設(shè)計[J]. 電子科學學刊, 1994, 16(5): 490-496.

Jing lin-jiao. Range ambiguity distribution characteristic of spaceborne SAR and its design consideration[J]. Journal of Electronics, 1994, 16(5): 490-496.

[13] 陳祺, 黃海風, 余安喜, 等. 一種改善交叉極化距離模糊度的新工作模式[J]. 電波科學學報, 2013, 28(6): 1082-1087.

Chen Qi, Huang Hai-feng, Yu An-xi, et al.. A new working mode to improve the cross-polarization range ambiguity to signal ratio[J]. Chinese Journal of Radio Science, 2013, 28(6): 1082-1087.

[14] Curlander J C and McDonough R N. Synthetic Aperture Radar Systems and Signal Processing[M]. Wiley Series in Remote Sensing, Hoboken, NJ: Wiley, 1991: 303-306.

[15] Lee J S and Pottier E. Polarimetric Radar Imaging from Basic to Application[M]. Boca Raton, FL, Taylor & Francis Group, 2009: 74-75.

洪 文： 女，1968年生，研究員，博士生導師，研究方向為SAR成像與系統(tǒng)及其應用、極化/極化干涉SAR數(shù)據(jù)處理及應用、3維微波成像新概念新體制新方法等.

楊士林： 男，1989年生，碩士生，研究方向為混合極化SAR系統(tǒng)仿真.

李 洋： 男，1983年生，博士生，研究方向為極化SAR、簡縮極化SAR、相干層析理論研究與應用.

張晶晶： 男，1986年生，博士生，研究方向為極化SAR、簡縮極化SAR、極化定標理論研究.

郭勝龍： 男，1987年生，博士生，研究方向為極化SAR/簡縮極化SAR理論研究與應用.

Study on Polarimetric SAR Range Ambiguity Computation for Distributed Targets

Hong Wen①②Yang Shi-lin①②③Li Yang①②③Zhang Jing-jing①②③Guo Sheng-long①②③

①(Science and Technology on Microwave Imaging Laboratory, Beijing 100190, China)

②(Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China)

③(University of Chinese Academy of Sciences, Beijing 100190, China)

A new approach is proposed to compute spaceborne SAR range ambiguity using covariance matrix. The causation of range ambiguity of spaceborne SAR and the ambiguous energy difference between like-polarized channel and cross-polarized channel are analyzed. In current study, it is common to use polarimetric scattering matrix to compute range ambiguity. Because of the speckle, the scattering matrix of a pixel can not describe the distributed targets precisely. Therefore, it is hard to decide the range ambiguity signal ratio of adjacent pixels. In this paper, the formulation of computing range ambiguity is derived for distributed targets first. Second, this new method is tested by using Radarsat-2 data. The results show that covariance matrix can effectively illustrate distributed targets. The results are smooth and confirmed. The new method is rational to compute range ambiguity of distributed targets.

SAR; Distributed targets; Polarimetric SAR; Range ambiguity

TN958

： A

：1009-5896(2015)06-1437-06

10.11999/JEIT141234

2014-09-23收到，2015-02-13改回

*通信作者：楊士林 slyangcau@163.com