向 聰 王亞軍 張開生

(西安電子工程研究所 西安 710100)

機載SAR方位空變補償算法

向 聰 王亞軍 張開生

(西安電子工程研究所 西安 710100)

當載機的理想航跡與實際航跡存在偏差時，方位上分布不同的目標對應的運動誤差歷程具有方位空變性。結(jié)合傳統(tǒng)子孔徑頻域補償?shù)乃枷?，提出了一種新的方位空變誤差精確補償方法。該方法采用時域卷積替代頻域點乘，在回波方位的慢時間域逐點對不同時刻、不同方位角目標的運動誤差進行精確補償。經(jīng)仿真和實測數(shù)據(jù)驗證表明，本文所提方法可在精確補償方位空變誤差的同時，能有效克服傳統(tǒng)補償方法存在的方位柵瓣問題，并具有運算復雜度低、利于工程實現(xiàn)的優(yōu)點。

運動誤差；方位空變誤差；子孔徑；方位柵瓣

Abstract: When there is deviation between ideal track and practical track of aircraft, the motion error course correspondent to different targets distributed in azimuth possesses azimuth-variant. A new azimuth-variant error precise compensation method is presented by combining idea of conventional subaperture compensation in frequency domain. This method employs time domain convolution to substitute frequency domain dot product to perform precise compensation for motion error of target in different azimuth at different time point by point in slow time domain of echo azimuth. The simulated and real measured data verifies that the presented method can be used to overcome azimuth grating lobe issue existed in conventional compensation method and provide advantageous like low operation complexity and benefit for implementation in project while compensating azimuth-variant error precisely.

Keywords：motion error; azimuth-variant error; subaperature; azimuth grating lobe

0 引言

合成孔徑雷達(SyntheticAperture Radar，SAR)可在各種惡劣戰(zhàn)場環(huán)境下全天時、全天候地獲取大面積的距離-方位二維高分辨地表圖像，其獲取距離高分辨的手段是通過發(fā)射大帶寬線性調(diào)頻信號，而方位高分辨則依賴于平臺與目標之間的相對運動以形成大虛擬陣列，從而突破天線尺寸對方位分辨率的限制[1]。

近年來，隨著SAR成像技術(shù)的不斷發(fā)展，其圖像質(zhì)量可與光學圖像相比擬，分辨率已能達到厘米級。分辨率的提高對自動目標檢測、特征提取、識別等能力的提升有很大幫助，極大地擴展其功能和應用。因而提高分辨率始終是雷達成像領域不懈追求的目標，但高分辨率同時意味著需要更高的運動補償精度[2-5]。

眾所周知，機載SAR成像要求載機平臺作勻速直線運動，然而受氣流或其他因素的影響，實際中勻直運動很難得到保證，不可避免的會出現(xiàn)偏離航線及橫滾、縱搖和偏流等非理想運動狀態(tài)，從而導致沿雷達視線(LOS)方向存在距離誤差、沿航向方向采樣不均勻等問題[3]，進而引起SAR圖像質(zhì)量下降，如分辨率降低、圖像對比度差，甚至無法成像。目前常用的運動誤差補償方法大多只能對波束中心的運動誤差進行補償，忽略了誤差在方位向的空變效應。此類方法僅適合于窄波束條件下，運動誤差量級較小或分辨率較低的情況。但對高分辨成像而言，相應的方位積累角較大，運動誤差的時變性將會導致目標方位譜產(chǎn)生畸變和移位，從而影響圖像的局部聚焦效果。

針對上述問題，本文結(jié)合傳統(tǒng)子孔徑頻域補償?shù)乃枷?，提出了一種新的方位空變誤差精確補償方法。該方法采用時域卷積替代頻域點乘，在回波方位的慢時間域逐點對不同時刻、不同方位角目標的運動誤差進行精確補償。經(jīng)仿真和實測數(shù)據(jù)驗證表明，本文所提方法可在精確補償方位空變誤差的同時，能有效克服傳統(tǒng)補償方法存在的方位柵瓣問題，并具有運算復雜度低、利于工程實現(xiàn)的優(yōu)點。

1 運動誤差方位空變性

SAR成像過程中不僅有方位空變誤差，還有距離空變誤差，限于篇幅本文不討論距離空變誤差補償。為了簡明的展現(xiàn)方位空變誤差，下面給出圖1所示的成像斜平面示意圖(波束內(nèi)任意點目標與航跡均能形成一個斜平面，可旋轉(zhuǎn)至同一平面)。

設平臺在成像孔徑時間內(nèi)從A點運動到B點，以A、B兩點的直接連線作為平臺參考航跡，平臺在某時刻實際位置為P點，此時在參考航跡上的等效位置為Q，且PQ=Δy(n)，則對波束內(nèi)任意目標F(方位角θ、斜距R)，運動誤差補償?shù)倪^程是將平臺從實際位置P補償?shù)降刃恢肣，則補償?shù)恼`差為：

(1)

方位波束中心處目標E(方位角θ0、斜距R) 的誤差為ΔR(θ0)，為波束內(nèi)共性誤差。根據(jù)空間幾何關系，不同方位角θ的目標，所需的誤差補償與波束中心誤差不同，即稱為方位空變誤差，可表示為：

Ψ=ΔR(θ)-ΔR(θ0)

(2)

當雷達波束較窄或運動誤差較小時，方位空變誤差可忽略，此時僅補償波束中心處目標的誤差即可(共性誤差)?？紤]到實際飛行過程中平臺存在飛行方向的非均勻采樣及高度向的運動誤差，因此還需增加平臺從實際位置到等效位置間還存在X向(飛行航向)誤差Δx(n)及Z向(垂直于XY平面)誤差Δz(n)。針對方位角為θ、作用距離為R的目標，將平臺從真實位置補償?shù)降刃恢脮r需補償?shù)目臻g距離差為：

(3)

忽略二次及二次以上項，上式可簡化為：

ΔR(θ)≈-Δx(n)sinθ-Δy(n)cosθ

(4)

由此可見，雷達波束較寬或運動誤差量級較大時，方位空變誤差的影響不可忽略。為不影響相參積累，上述距離差對相位的影響應小于pi/4，對包絡的影響應小于距離分辨單元的一半。實際中，方位空變誤差對包絡的影響往往可以忽略，僅考慮相位誤差即可。

2 傳統(tǒng)子孔徑精確補償算法

2.1 方位空變誤差精確補償原理

下面以正側(cè)視為例，在無運動誤差的情況下，分析方位信號多普勒譜與成像平面中方位波束角的對應關系。假設系統(tǒng)的脈沖重復頻率為fr，雷達波長為λ，平臺速度為va，回波信號方位脈沖個數(shù)為N。對方位信號進行FFT變換，根據(jù)多普勒譜形成原理，在成像平面內(nèi)與航跡法線夾角為θ的波束對應的多普勒頻率為：

fa=2vasinθ/λ

(5)

則相應的波束指向角與多普勒頻率的對應關系有：

(6)

顯然對于方位多普勒譜而言，每個頻點可以對應波束的局部范圍。

存在運動誤差時，將第n個脈沖時刻平臺實際位置等效到對應參考位置，針對不同頻率fa(n)需要補償?shù)目臻g差為ΔR(n，fa(n))，才能實現(xiàn)對各方位角的目標實現(xiàn)精確補償，從而不僅避免造成方位譜畸變，同時消除了誤差空變性對方位相位和距離向徙動造成的影響。

假設對距離頻域-方位時域的回波信號進行方位FFT變換，可表示為：

S(fr，fa)=fft(s(fr，tm))

(7)

(8)

上式相位中，第一項為包絡補償，第二項為相位補償。以上即為方位空變誤差精確補償?shù)幕驹怼?/p>

忽略方位空變誤差對包絡的影響，上述誤差補償因子可以簡化為：

(9)

由此，方位空變誤差的補償僅需在距離時域-方位多普勒域進行，其補償過程可表示為：

(10)

2.2 子孔徑頻域補償算法

合成孔徑雷達單個回波脈沖對應了波束范圍內(nèi)所有目標，平臺的運動誤差對在波束范圍內(nèi)的目標產(chǎn)生了誤差差異，但在單個脈沖內(nèi)，無法通過頻域變換實現(xiàn)不同方位角目標的補償。因此，需要將方位全孔徑數(shù)據(jù)劃分成若干子孔徑處理，并假設子孔徑內(nèi)的誤差恒定，采用子孔徑交疊滑窗處理的方式進行補償，處理示意圖如圖2 所示。

圖2 分子孔徑補償示意圖

從方位空變誤差補償?shù)男Ч紤]，子孔徑長度N值越大，意味著波束劃分的越細，空變誤差補償?shù)母?，但考慮運動誤差的時變性，子孔徑內(nèi)誤差恒定的假設越偏離實際。根據(jù)相干積累原理，當相位誤差小于pi/4時基本不影響脈沖壓縮特性，則對于N的選取條件是使得相鄰子波束間的空變誤差小于λ/16，即N點FFT變換時有

ΔR(fa(n))-ΔR(fa(n+1))<λ/16

(11)

此外，子段間的步進量M選擇，決定了算法實現(xiàn)的運算量。M選擇過大，雖然運算量較小，但子孔徑間的誤差量階躍變化越明顯，同時回波信號的相位也將出現(xiàn)階躍變化，從而導致方位脈壓后，SAR圖像中的強點出現(xiàn)周期性柵瓣的問題。而M選擇太小，其實時運算量太大，往往不能被接受。

理論上，當上述基于子孔徑補償方法的步進量M為1時，其補償效果最好，但其運算量極為龐大，不利于工程實現(xiàn)。針對上述問題，下面提出了一種新的方位空變誤差精確補償方法。該方法采用時域卷積替代頻域點乘，在回波方位的慢時間域逐點對不同時刻、不同方位角目標的運動誤差進行精確補償，可在不增加運算量的基礎上，有效克服傳統(tǒng)補償方法存在的方位柵瓣問題。

3 改進的方位空變誤差精確補償算法

3.1 改進算法

首先，介紹一下循環(huán)卷積與傅里葉變化的關系。

對任意兩個有限長序列x1(n)和x2(n)，長度均為N。設x1(n)和x2(n)的N點FFT變換分別為：X1(k)=FFT[x1(n)]和X2(k)=FFT[x2(n)]。如果X(k)=X1(k)·X2(k)，則有：

x(n)=IFFT[X(k)]=x1(n)?x2(n)

(12)

其中?表示循環(huán)卷積。

由此，上述(10)式采用循環(huán)卷積的原理可以等價表示如下：

(13)

實時補償時，同樣需要將方位全孔徑數(shù)據(jù)劃分成若干子孔徑處理，假設第k個子孔徑內(nèi)，方位空變誤差補償前與補償后的回波可以寫成如下向量的形式：

Sk=[s-N/2，s-N/2+1，…，s0，…，sN/2-1]

(14)

(15)

(16)

(17)

令子孔徑間的步進量M為1，每次只針對該孔徑中間時刻t0處的回波進行補償。上式中的循環(huán)卷積并不需要全部計算，只需計算中間時刻對應的補償系數(shù)即可，其處理過程如圖3所示。補償后中間時刻對應回波數(shù)據(jù)為：

(18)

圖3 改進算法的處理流程圖

3.2 計算量對比

假設全孔徑長度的方位總脈沖數(shù)為Na，距離像采樣點為M。采用傳統(tǒng)方法與本文所提方法時，子孔徑長度均為N，而傳統(tǒng)方法子孔徑間的步進量為K，本文所提方法的步進量為1，兩種算法涉及的FFT、IFFT和點乘運算量對比如下表所示：

(19)

本文所提方法總的運算復雜度為：

(20)

兩者運算量之比為：

(21)

4 仿真實驗

為了驗證本文方位空變誤差補償方法的實際有效性，選取一塊運動誤差較大的Ku波段實測回波數(shù)據(jù)進行處理。主要系統(tǒng)參數(shù)如下表所示。

表2 雷達及載機主要參數(shù)

由慣導計算，成像斜平面內(nèi)的徑向和沿航向距離誤差分別為圖4(a)和(b)所示。

圖4和圖5分別給出了徑向誤差最大時(-140m)和航向誤差最大時(-15.8m)，波束內(nèi)不同角度對應的包絡和相位空變誤差。

由此可以看出，徑向誤差導致的合成孔徑時間內(nèi)包絡影響0.07m左右，可以忽略，但相位影響需考慮。由航向誤差導致的合成孔徑時間內(nèi)包絡影響0.6m左右，1m以上分辨率可以忽略，但相位影響需考慮。

圖4 載機運動誤差量級

圖5 最大徑向誤差對包絡和相位的影響

圖6 最大航向誤差對包絡和相位的影響

圖7為傳統(tǒng)頻域補償方法進行方位空變誤差補償后的成像結(jié)果，其中子孔徑長度為256，子孔徑間的步進量為16。圖8為采用本文所提方法進行方位空變誤差補償后的成像結(jié)果。從兩圖對比中可以明顯發(fā)現(xiàn)，在聚焦效果相同的情況下，本文所提方法有效消除了方位柵瓣的影響。

5 結(jié)論

本文基于時域循環(huán)卷積結(jié)合運動誤差補償?shù)姆椒?，從原理上很好的解決了傳統(tǒng)頻域補償方法成像過程中的方位柵瓣問題，通過理論分析和實測機載〗數(shù)據(jù)的驗證本文所提方法的有效性。

圖7 傳統(tǒng)頻域補償方法的成像結(jié)果

圖8 本文所提方法的成像結(jié)果

[1] 于春銳，張永勝，董臻，梁甸農(nóng).基于特征分解的SAR射頻干擾抑制方法[J]．信號處理，2011，27 ( 11 )：1696-1700．

[2] 李建陽，常文革，黎向陽．二階運動補償對低頻超寬帶SAR視線誤差補償?shù)挠绊慬J]．電子與信息學報，2011，33 (4)：792-797．

[3] 安道祥，黃曉濤，李欣，周智敏，齊星．機載超寬帶SAR運動補償方法[J]．信號處理，2011，27 (1)：73-80．

[4] 李楊寰，宋千，周智敏．一種適合低空機載平臺SAR成像運動補償方法[J]．信號處理，2010，26 ( 4 )：545-551．

[5] 張雙喜，孫光才，周峰，劉艷陽，邢孟道．一種基于子孔徑自聚焦的高頻運動誤差估計和補償方法[J].電子與信息學報，2010，32 (12)：3013-3017．

AirborneSARAzimuth-variantCompensationAlgorithm

Xiang Cong， Wang Yajun， Zhang Kaisheng

(Xi’an Electronic Engineering Research Institute， Xi’an 710100)

TN959

A

1008-8652(2017)02-042-06

2017-04-26

向聰(1985-) ，男，博士研究生。研究方向為陣列信號處理，機載空時二維處理，SAR成像技術(shù)等。