吳玉峰 孫光才 楊 軍 邢孟道 周 峰

(西安電子科技大學雷達信號處理國家重點實驗室 西安 710071)

1 引言

相對機載SAR而言，星載SAR要提供全球覆蓋的觀測區(qū)域，為了縮短對場景的再訪時間以及獲得場景的細節(jié)特征，對地面進行高分辨寬測繪帶的成像顯得尤為重要。星載SAR由于受最小天線面積的限制[1]，高分辨與寬測繪帶構成一對矛盾，現(xiàn)有的不同SAR成像模式僅在某一方面取其所長。采用聚束模式和發(fā)射寬頻帶信號可以獲得方位和距離高分辨率，但是其觀測帶不連續(xù)；ScanSAR利用不同脈沖時間控制波束指向照射不同子測繪帶以此獲得超寬地面觀測帶，然而它是以犧牲方位分辨率為代價的；條帶SAR則在兩者之間折中只能得到中等的分辨率和測繪帶寬度。

為了解決高分辨與寬測繪帶之間的矛盾，國內(nèi)外學者提出了很多方法。文獻[2]提出高PRF采樣回波信號，利用距離維多個子孔徑進行波束形成解距離模糊，然而該方法得到的寬測繪帶存在盲區(qū)。文獻[3]利用低PRF采樣，方位向多個子孔徑滿足偏置相位中心(DPC)條件實現(xiàn)高分辨寬測繪帶成像，然而其工作條件要求苛刻，條件不滿足時系統(tǒng)性能損失嚴重。在 DPC條件不滿足的情況下，文獻[4]通過空域濾波方法解多普勒模糊實現(xiàn)高分辨寬測繪帶成像。文獻[5-7]考慮了寬測繪帶回波信噪比問題，通過多維波形編碼及數(shù)字波束形成技術控制各子測繪帶子脈沖的時間寬度和帶寬，從而得到靈活的能量分配和多距離分辨的特性。文獻[8-10]利用頻域子頻帶合成提高距離向分辨率，并通過方位多波束實現(xiàn)寬測繪帶成像。以上方法都是基于多通道的高分辨寬測繪帶成像方法。然而，多通道SAR系統(tǒng)本身面臨著一些問題：(1)系統(tǒng)復雜度，多通道 SAR系統(tǒng)通過Q個通道最多可解Q次模糊，當模糊數(shù)增加時就需要更多的天線通道，這將大大增加系統(tǒng)的復雜度；(2)通道間距，常規(guī)多通道 SAR系統(tǒng)若采用分布式小衛(wèi)星則其間距是不可控的，即使采用精確的測量技術其誤差也不容忽略，若采用天線陣列則其天線系統(tǒng)將顯得很復雜，實現(xiàn)較難；(3)通道一致性，由于模擬器件及其電路的種種原因導致各個通道頻率響應特性不一致，需要進行額外的通道均衡方法來抑制其不利影響。

相對多通道SAR系統(tǒng)而言，單通道SAR系統(tǒng)就不存在這些問題，但是目前對星載單通道高分辨寬測繪帶成像方法的研究還比較少。文獻[11]論述了一種基于多調制方式的寬測繪帶方法，但其數(shù)據(jù)處理復雜且未考慮方位模糊的問題。文獻[12]通過連續(xù)變化的PRF使距離盲區(qū)出現(xiàn)在不同的距離單元內(nèi)，實現(xiàn)單通道寬測繪帶成像，然而其方位分辨率受限并且也未考慮方位模糊的問題。

本文提出了一種利用方位周期性非均勻采樣實現(xiàn)單通道高分辨寬測繪帶成像的方法。通過非均勻采樣序列的優(yōu)化設計，避開距離盲區(qū)重疊的現(xiàn)象，使得相同距離單元的盲區(qū)最多出現(xiàn)在一個采樣通道內(nèi)；然后利用方位采樣的周期性，構造出等效的多通道數(shù)據(jù)，利用多通道解模糊的方法實現(xiàn)方位頻譜恢復。文章的結構如下：第2節(jié)簡要分析星載SAR的信號模型以及采用單通道實現(xiàn)高分辨寬測繪帶成像時存在的問題；第3節(jié)針對上述問題給出解決方法；第4節(jié)詳細論述非均勻采樣序列的優(yōu)化設計方法，并給出信號的處理流程；第5節(jié)給出仿真結果，驗證本文方法的可行性；第6節(jié)對全文進行了總結。

2 信號分析

圖1 星載SAR工作幾何關系

如圖1所示為星載SAR工作幾何示意圖，衛(wèi)星平臺到場景中心的距離為RS，平臺等效速度為V，高度為H，以平臺運動方向為X軸建立如圖所示的坐標系。星載平臺收發(fā)不在同一位置，可以利用等效相位中心原理表示雷達的采樣位置。

假設雷達采用脈沖體制，斜距平面內(nèi)存在任意點目標P(RB,Xn),RB為目標到雷達航跡的最近距離，Xn為方位位置。則雷達接收到的回波信號經(jīng)脈沖壓縮后在距離頻域可以表示為

其中fr為距離頻率，B為發(fā)射信號帶寬，c為光速，fc為雷達載頻，tm為方位慢時間，R(tm;RB)為點目標P到雷達的瞬時斜距，rect(fr/B)為發(fā)射信號在頻域的窗函數(shù)表示形式。本文的重點集中于采用單通道實現(xiàn)高分辨寬測繪帶成像，并介紹相應的信號處理方法。對于單個通道而言，如果采用常規(guī)的方位均勻采樣方式錄取回波信號，不可避免會出現(xiàn)距離或方位模糊。

傳統(tǒng)的單通道星載 SAR如果要實現(xiàn)寬測繪帶成像，必須要求采用較小的PRF以避免距離出現(xiàn)盲區(qū)和模糊；而若要實現(xiàn)高分辨成像，又必須采用較大的PRF以保證多普勒不模糊。因此，采用常規(guī)的單通道 SAR系統(tǒng)實現(xiàn)高分辨寬測繪帶成像必然會出現(xiàn)距離模糊或方位模糊的問題，或兩者同時存在。方位模糊的雷達回波信號可以表示為

其中S(fr,fa)為信號在兩維頻域的表示形式，fa∈[-PRF/2:PRF/2]為方位頻率，M為方位模糊數(shù)。式(2)以方位奇數(shù)次模糊為例。

3 解模糊方法

3.1 雷達工作體制

本文雷達采用常規(guī)的單通道體制，不同的是本文采用方位周期性非均勻采樣的方式錄取數(shù)據(jù)。如圖2所示為雷達方位周期性非均勻采樣的示意圖。其中同種線型的‘×’表示重復頻率為PRF的發(fā)射位置，其間隔為VPRF 。一個VPRF間隔內(nèi)，不同種線型的‘×’表示非均勻的采樣位置，圖中以 4個非均勻采樣點為例。4種線型的‘×’相對于實線型‘×’的距離記為di(i= 0,1,2,3)，其中d0=0。雷達在‘×’時刻發(fā)射設計好的寬帶信號，在不發(fā)射信號時始終處于接收狀態(tài)。將接收的回波信號進行重新排列，不同發(fā)射點對應的回波構成的兩維數(shù)據(jù)記為Si，這些數(shù)據(jù)可以認為是4個通道的數(shù)據(jù)。

圖2 周期性非均勻采樣示意圖

3.2 距離解模糊

對于寬測繪帶成像，本文采用較小的脈沖重復頻率 PRF，保證距離向無盲區(qū)和模糊。然而 PRF間隔內(nèi)天線還需要發(fā)射信號，其相應位置會因為發(fā)射遮擋出現(xiàn)盲區(qū)，距離盲區(qū)的形成如圖3所示。

圖3 距離盲區(qū)形成示意圖

從圖中可以看出，實線型‘×’發(fā)射點的回波在接收過程中，天線必須發(fā)射信號，對于同一個收發(fā)天線而言，天線處于發(fā)射狀態(tài)時不能接收回波信號，因此數(shù)據(jù)錄取時會出現(xiàn)距離盲區(qū)問題，圖中陰影部分即為實線型‘×’發(fā)射脈沖回波包絡內(nèi)由于發(fā)射遮擋導致的3個距離盲區(qū)。通過非均勻采樣序列的設計，可以保證同一距離單元的盲區(qū)只出現(xiàn)在一個通道中，具體的設計將在后面作細致分析。同時還可以看出，不同線型發(fā)射脈沖的回波將重疊在一起，對于這一問題，可以采用發(fā)射正交多波形的方法解決。

3.3 方位解模糊[13]

對于不同的通道而言，其數(shù)據(jù)方位向是模糊的，然而 PRF內(nèi)的非均勻采樣提供了足夠的信息來消除模糊。根據(jù)文獻[4,13]的討論，如果通道數(shù)N滿足N≥M，那么可以利用方位的多個偏置通道解方位模糊。由于本體系中的多個通道實際上是單個傳感器通過方位非均勻采樣來實現(xiàn)的，它構成的多通道系統(tǒng)比常規(guī)的多通道系統(tǒng)存在很多的優(yōu)勢，比如通道一致性好、通道間距可控、測繪帶寬度不會對硬件產(chǎn)生很大的要求等等，因此該體制將為系統(tǒng)設計以及信號恢復提供很大的便利。

假設N個通道的數(shù)據(jù)可以排列成向量的形式

其中R為方位多普勒數(shù)據(jù)構成的協(xié)方差矩陣。式(4)中第1個式子體現(xiàn)了采樣系統(tǒng)對多普勒分量fa+m0× P RF的增益為 1，第 2個式子表示系統(tǒng)對其它模糊分量具有抑制作用，第3個式子是指在空間自由度足夠滿足上面兩個式子后，對副瓣的抑制作用。

對某個距離單元，若第n0個通道存在距離盲區(qū)，那么去掉該通道的導向矢量，ANM改寫為AN′M，其中N′=N- 1，同樣可以利用上面的式子恢復信號。

4 應用分析

4.1 重復頻率設計

方位重復頻率PRF要滿足距離無盲區(qū)的要求，并且保證發(fā)射信號不會對自身采樣通道的回波構成距離盲區(qū)，可以表示為

其中Wr為場景的斜距幅寬，mr為臨近的自身采樣通道發(fā)射脈沖個數(shù)。

4.2 采樣序列設計

在 PRF滿足式(5)的基礎上，采樣點間距di的設計對于避開距離盲區(qū)的重疊和信號占空比都非常重要，可以用圖4說明。假設方位向有3組非均勻采樣點，圖中第1行表示發(fā)射信號時序，第2行為雷達接收到的回波，由于在不發(fā)射信號時雷達始終處于接收狀態(tài)，因此除了由發(fā)射遮擋導致的距離盲區(qū)(圖中陰影部分所示)外，其他時間都存在回波信號。第3至5行分別表示等效通道1至通道3的接收回波，其中粗實線對應同一PRF?？梢钥闯?，距離盲區(qū)的存在將原始回波數(shù)據(jù)分成了3段。為了使各通道回波數(shù)據(jù)起始采樣時刻對齊后距離盲區(qū)不重疊，必須設計好非均勻采樣點的位置。圖4中的第3至5行可以看作3個線段的循環(huán)移位，因此采樣點的設計實質上可以看作是3個線段長度的設計。每一行中只有兩個距離盲區(qū)，其位置可以通過累加操作求得。

圖4 采樣通道與距離盲區(qū)關系

因此對于一般的非均勻采樣而言，存在序列向量l= [n1,n2,…,nN]T較優(yōu)的組合，滿足

圖5 距離盲區(qū)疊加

由于序列移位累加點對線段進行均分，因此其累加點必須充滿整個區(qū)域。從而可以得到如下規(guī)律：數(shù)1必須包含在內(nèi)以保證第1段的存在；數(shù)2必須包含在內(nèi)以保證第2段的存在；如果數(shù)2與數(shù)1相連，則第3段被它們組合，那么數(shù)4肯定存在，如果數(shù)2與數(shù)1不相連，數(shù)3存在，等等……。上面的規(guī)律帶有遞歸性質，可以通過遞歸結構來實現(xiàn)。

通過遞歸結構可以得到N為不同值時的序列組，如表1所示。

表1 非均勻采樣序列舉例

以表 1中序列{1,2,5,4,6,13}為例，通過循環(huán)移位及累加，可以得到如圖6所示的結果，可以看出，6×5+1=31被均勻地充滿了，且間距為1。此時可以得到發(fā)射信號脈沖寬度的上限為1/(2PRF)?[1/(N(N-1)+1)]，就可以避開距離盲區(qū)的重疊。

實際系統(tǒng)中由于存在發(fā)射信號的延遲抖動，可能會導致序列不能嚴格按照表1所給的比例，此時是否會造成距離盲區(qū)重疊呢？答案是否定的，考慮到延遲抖動，假設上述序列變?yōu)閧1.08,1.84,5.13,4.26, 5.87,12.81}，此時序列循環(huán)移位及累加結果如圖 7，雖然不是均勻分布，但各通道的盲區(qū)依然相互分開，只要采用合適的發(fā)射信號寬度就可以避免距離盲區(qū)的重疊。由于發(fā)射信號的延遲是通過數(shù)字計數(shù)器實現(xiàn)的，其延遲抖動也已知，因此不會對方位解模糊產(chǎn)生影響。

4.3 數(shù)據(jù)處理流程

通過上述非均勻采樣序列的優(yōu)化設計，可以避免各個通道距離盲區(qū)的重疊。將回波數(shù)據(jù)按照距離門延遲及發(fā)射脈沖延遲重新排列，得到等效的多通道數(shù)據(jù)后，就可以進行后續(xù)的處理。

圖8為本文方法的數(shù)據(jù)處理流程，圖中假設方位向有4組非均勻采樣點。接收到的回波數(shù)據(jù)重新排列后可以形成4個等效通道的數(shù)據(jù)。各通道數(shù)據(jù)按各自的發(fā)射波形完成距離向匹配濾波后，變換到方位頻率域，此時方位頻譜存在模糊。根據(jù)發(fā)射脈沖時序可以計算出距離盲區(qū)的位置，若某距離單元全部通道都不存在盲區(qū)，則采用全部通道進行方位解模糊；若某一通道存在距離盲區(qū)，則去掉該通道的數(shù)據(jù)，采用剩余的其他通道解方位模糊。方位解模糊完成后，將所有距離向數(shù)據(jù)合成可以得到完整的且方位不模糊的數(shù)據(jù)。此時距離向補償回一個與發(fā)射信號相同時寬和帶寬的線性調頻信號，恢復距離向匹配濾波前的數(shù)據(jù)，就可以采用 CS算法實現(xiàn)距離脈沖壓縮、二次距離脈沖壓縮及距離單元徙動(RCM)校正。最后通過方位脈沖壓縮就可以得到無模糊的SAR圖像。

5 仿真結果

為了驗證本文提出的高分辨寬測繪帶成像方法，本節(jié)給出了仿真數(shù)據(jù)處理結果。雷達工作方式采用本文提出的周期性非均勻采樣方式，并且通過發(fā)射正交多波形、方位解模糊失效及方位脈壓失配解決距離模糊的問題。仿真時考慮方位模糊數(shù)為3，為了解決距離盲區(qū)內(nèi)的模糊問題，非均勻采樣點至少為4個。考慮到效率，仿真中采用的雷達發(fā)射波形為如圖9所示的4種近似正交波形。實際應用中可以選擇其他更加優(yōu)化的波形，如正交的離散頻率編碼波形(DFCW)[14]作為發(fā)射波形。經(jīng)過分析可知，圖中所示信號的自相關峰值旁瓣約為-13.3 dB，而互相關系數(shù)約為-25 dB，由于方位脈沖壓縮失配，將導致距離模糊信號的幅度進一步降低，達到約-30 dB，可以認為不會對成像結果造成影響。

圖6 序列循環(huán)移位及累加結果

圖7 考慮延遲抖動后移位累加結果

圖8 數(shù)據(jù)處理流程

圖9 近似正交波形圖

表2 仿真參數(shù)

仿真所用參數(shù)如表2所示。根據(jù)參數(shù)的分析，方位信號帶寬為3000 Hz，模糊了3次，因此非均勻采樣序列采用表 1中的{1,7,2,3}，對應的采樣點位置為di={0,1,8,10}/13 ×V/PRF ，該采樣點序列可以有效避免距離盲區(qū)的重疊，使得相同距離單元的盲區(qū)最多出現(xiàn)在一個通道內(nèi)。

仿真時距離向場景寬度設置為96 km，均勻布置了65個點，點的間距為1.5 km。為了更加有效地驗證本文方法對方位點的有效性，方位坐標值隨機分布在(-1000,1000)內(nèi)。圖10為仿真產(chǎn)生的原始數(shù)據(jù)圖(為了便于顯示，距離向進行了8抽1處理)，可以看出，每個采樣通道都存在距離盲區(qū)，但是相同距離單元的盲區(qū)只出現(xiàn)在一個通道內(nèi)，這樣使得我們可以最大限度的利用方位采樣點的資源。

圖10 仿真原始數(shù)據(jù)圖(距離向8抽1)

完成方位解模糊及成像處理后，部分點目標成像結果如圖11所示，其中水平方向為距離向，垂直方向為方位向。從圖中可以看出，場景中的點目標聚焦良好。

進一步取各個模糊區(qū)域內(nèi)點目標成像結果進行分析，如圖 12所示，其中圖 12(a)～12(e)分別對應圖10距離帶1～5內(nèi)點目標成像結果的幅度圖及方位剖面圖。從圖中可以看出，無論有無盲區(qū)及哪個采樣通道距離盲區(qū)內(nèi)的點目標通過方位解模糊之后，模糊信號基本都在-40 dB以下，因此不會對成像結果產(chǎn)生影響。表3給出了上述5個點目標方位向的成像性能分析，其中PSLR為峰值旁瓣比，ISLR為積峰旁瓣比?？梢钥闯龇轿幌騊SLR均在-13.01 dB以下，ISLR均在-9.63 dB以下，成像分辨率也與理論值接近，因此成像結果較為理想。

表3 各點目標方位向成像性能分析

6 結論

本文針對傳統(tǒng)星載單通道 SAR系統(tǒng)高分辨與寬測繪帶之間的矛盾，提出了一種方位周期性非均勻采樣的解決方法。該方法通過對非均勻采樣序列的優(yōu)化設計，避免了距離盲區(qū)的重疊，使得相同距離單元的盲區(qū)只出現(xiàn)在一個采樣通道內(nèi)。利用方位采樣的周期性，構造出等效的多通道數(shù)據(jù)，通過多通道解模糊的方法實現(xiàn)方位信號頻譜恢復。該方法可以實現(xiàn)星載單通道 SAR系統(tǒng)高分辨寬測繪帶成像，且不存在距離盲區(qū)問題。仿真實驗驗證了本文方法的有效性。

圖11 場景點目標成像結果(部分)

圖12 不同通道盲區(qū)內(nèi)點目標成像結果及方位剖面圖

[1]Freeman A, Johnson W T K, Huneycutt B,et al.. The“Myth” of the minimum SAR antenna area constraint [J].IEEE Transactions on Geoscience and Remote Sensing, 2000,38(1): 320-324.

[2]Callaghan G D and Longstaff I D. Wide-swath space-borne SAR and range ambiguity [C]. Radar 97, Edinburgh, UK,October 14-16, 1997: 248-252.

[3]Currie A and Brown M A. Wide-swath SAR [J].IEE Proceedings-Radar,Sonar,Navigation, 1992, 139(2):122-135.

[4]邢孟道, 李真芳, 保錚, 等. 分布式小衛(wèi)星雷達空時頻成像方法研究[J]. 宇航學報, 2005, 26(增刊1): 70-76.Xing Meng-dao, Li Zhen-fang, Bao Zheng,et al.. Study of distributed microsatellites radar space-time-frequency imaging method [J].Journal of Astronautics, 2005,26(Suppl.1): 70-76.

[5]Krieger G, Gebert N, and Moreira A. Multidimensional waveform encoding: a new digital beamforming technique for synthetic aperture radar remote sensing [J].IEEE Transactions on Geoscience and Remote Sensing, 2008, 46(1):31-46.

[6]Krieger G, Younis M, Gebert N,et al.. Advanced concepts for high-resolution wide-swath SAR imaging [C]. EUSAR 2010,Aachen, Germany, June 7-10, 2010: 524-527.

[7]齊維孔, 禹衛(wèi)東, 祁海明. 星載MIMO-SAR與距離向DBF相結合系統(tǒng)研究[J]. 電子學報, 2010, 38(10): 2251-2257.Qi Wei-kong, Yu Wei-dong, and Qi Hai-ming. Study of the system combining spaceborne MIMO-SAR and elevation DBF [J].Acta Electronica Sinica, 2010, 38(10): 2251-2257.

[8]Zaugg E, Long D, Edwards M,et al.. Coherent multi-frequency-band resolution enhancement for synthetic aperture radar [C]. IGARSS 2009, Cape Town, South Africa,July 12-17, 2009: I-56-I-59.

[9]鄧云凱, 陳倩, 祁海明, 等. 一種基于頻域子帶合成的多發(fā)多收高分辨率SAR 成像算法[J]. 電子與信息學報, 2011, 33(5):1082-1087.Deng Yun-kai, Chen Qian, Qi Hai-ming,et al.. A highresolution imaging algorithm for MIMO SAR based on the sub-band synthesis in frequency domain[J].Journal of Electronics&Information Technology, 2011, 33(5):1082-1087.

[10]黃平平, 鄧云凱, 徐偉, 等. 基于頻域合成方法的多發(fā)多收SAR技術研究[J]. 電子與信息學報, 2011, 33(2): 401-406.Huang Ping-ping, Deng Yun-kai, Xu Wei,et al.. The research of multiple-input and multiple-output SAR based on frequency synthetic [J].Journal of Electronics&Information Technology, 2011, 33(2): 401-406.

[11]王小青, 朱敏慧. 一種寬測繪帶新方法的探討[J]. 電子與信息學報, 2003, 25(10): 1425-1429.Wang Xiao-qing and Zhu Min-hui. A discussion on a new method of wide swath SAR [J].Journal of Electronics&Information Technology, 2003, 25(10): 1425-1429.

[12]Gebert N and Krieger G. Ultra-wide swath SAR imaging with continuous PRF variation [C]. EUSAR 2010, Aachen,Germany, June 7-10, 2010: 966-969.

[13]Li Zhen-fang, Wang Hong-yang, Su Tao,et al.. Generation of wide-swath and high-resolution SAR images from multichannel small spaceborne SAR systems [J].IEEE Geoscience and Remote Sensing Letters, 2005, 2(1): 82-86.

[14]Deng Hai. Discrete frequency-coding waveform design for netted radar systems [J].IEEE Signal Processing Letters,2004, 11(2): 179-182.