一種基于脈沖壓縮的機(jī)載條帶SAR重疊子孔徑實(shí)時(shí)成像算法

詹學(xué)麗王巖飛 王 超 劉碧丹

(中國(guó)科學(xué)院電子學(xué)研究所 北京 100190)

一種基于脈沖壓縮的機(jī)載條帶SAR重疊子孔徑實(shí)時(shí)成像算法

詹學(xué)麗*王巖飛 王 超 劉碧丹

(中國(guó)科學(xué)院電子學(xué)研究所 北京 100190)

該文提出了一種基于距離向脈沖壓縮體制的機(jī)載條帶SAR重疊子孔徑(OSA)實(shí)時(shí)成像算法，可直接應(yīng)用于未采用去斜接收技術(shù)的條帶機(jī)載SAR實(shí)時(shí)成像處理系統(tǒng)。該算法在距離向采用脈沖壓縮技術(shù)、方位向采用OSA信號(hào)處理方法，能夠在成像過(guò)程中完成距離徙動(dòng)校正并消除空變相位誤差，得到條帶模式下的高分辨率圖像。首先分析了條帶模式SAR幾何關(guān)系，然后在脈沖壓縮體制下建立了回波模型并對(duì)成像流程進(jìn)行了詳細(xì)的數(shù)學(xué)推導(dǎo)，最后對(duì)運(yùn)算量、存儲(chǔ)量與成像范圍進(jìn)行了計(jì)算分析，對(duì)點(diǎn)目標(biāo)數(shù)據(jù)進(jìn)行了仿真，利用實(shí)時(shí)處理平臺(tái)對(duì)原始數(shù)據(jù)進(jìn)行了成像，驗(yàn)證了算法的高時(shí)效性。

合成孔徑雷達(dá)(SAR)；實(shí)時(shí)成像；脈沖壓縮；重疊子孔徑算法(OSA)

1 引言

重疊子孔徑算法[1](Overlapped Subaperture Algorithm, OSA)是一種高分辨率合成孔徑雷達(dá)(Synthetic Aperture Radar, SAR)成像算法。OSA在時(shí)域?qū)Ψ轿幌蚧夭ㄐ盘?hào)進(jìn)行重疊子孔徑劃分，然后在子孔徑內(nèi)與子孔徑間的多次處理過(guò)程中引入誤差補(bǔ)償因子，完成方位向的相位誤差補(bǔ)償及距離向與方位向的耦合誤差補(bǔ)償，生成高分辨率雷達(dá)圖像。OSA成像中對(duì)子孔徑進(jìn)行壓縮處理后，可以得到每個(gè)子孔徑的方位向粗分辨率信息，這時(shí)子孔徑內(nèi)處于圖像域、子孔徑間處于時(shí)域，這種既在圖像域又在時(shí)域是OSA獨(dú)有的特點(diǎn)，利用這一特點(diǎn)可以采用圖像域的粗分辨率信息，在成像過(guò)程中補(bǔ)償子孔徑間時(shí)域的相位誤差。OSA與距離多普勒(Range Doppler, RD)算法、線性調(diào)頻變標(biāo) (ChirpScaling, CS)算法相比，通過(guò)在方位向的分塊處理能夠降低對(duì)存儲(chǔ)量的需求、采用方位向的子孔徑內(nèi)與子孔徑之間的短長(zhǎng)度快速傅里葉變換(FFT)有效降低了運(yùn)算量，更加適合于機(jī)載SAR實(shí)時(shí)成像。

文獻(xiàn)[2,3-5]針對(duì)去斜接收的聚束SAR系統(tǒng)提出了能夠校正空變相位誤差的OSA算法，可根據(jù)雷達(dá)系統(tǒng)參數(shù)與姿態(tài)/位置測(cè)量參數(shù)實(shí)時(shí)計(jì)算相位補(bǔ)償因子，有效解決了圖像的幾何失真與圖像空變散焦問(wèn)題，并實(shí)際應(yīng)用到Lynx型無(wú)人機(jī)載SAR雷達(dá)實(shí)時(shí)成像處理器中[6]。文獻(xiàn)[7-11]建立了極坐標(biāo)下的聚束SAR幾何模型，對(duì)聚束工作模式下的OSA成像算法進(jìn)行了研究。

OSA算法大部分應(yīng)用于去斜接收的聚束SAR系統(tǒng)，但是去斜技術(shù)限制了SAR系統(tǒng)的條帶寬度，無(wú)法獲取地面的大面積測(cè)繪范圍；另外條帶SAR是成像系統(tǒng)中常用的一種工作模式，其應(yīng)用更簡(jiǎn)便直觀，因此是SAR系統(tǒng)的基本工作模式。

本文重點(diǎn)提出了基于距離向脈沖壓縮體制的機(jī)載條帶OSA實(shí)時(shí)成像算法，該算法可應(yīng)用于未采用去斜接收技術(shù)的條帶SAR系統(tǒng)，可高效地得到條帶模式圖像。本文首先分析了條帶模式SAR幾何關(guān)系，然后給出了接收端去斜與非去斜時(shí)的SAR回波模型，針對(duì)非去斜接收系統(tǒng)對(duì)距離向脈沖壓縮體制的OSA成像流程進(jìn)行了詳細(xì)的數(shù)學(xué)推導(dǎo)與說(shuō)明，給出了算法的運(yùn)算量、存儲(chǔ)量與有效成像范圍，最后通過(guò)對(duì)點(diǎn)目標(biāo)仿真數(shù)據(jù)與實(shí)際回波數(shù)據(jù)的處理對(duì)算法時(shí)效性進(jìn)行了驗(yàn)證。

2 脈沖壓縮體制條帶SAR的原始數(shù)據(jù)回波模型

基于去斜接收的聚束式SAR OSA成像算法在文獻(xiàn)[3]中有詳細(xì)介紹，其簡(jiǎn)要描述見(jiàn)附錄。本文基于脈沖壓縮的機(jī)載條帶SAR系統(tǒng)提出了OSA實(shí)時(shí)成像算法，本節(jié)主要描述脈沖壓縮體制條帶SAR的原始數(shù)據(jù)回波模型。

圖1為SAR幾何關(guān)系示意圖。在圖1中，rb為天線波束中心到地面的垂直距離；rs為平臺(tái)到點(diǎn)目標(biāo)的瞬時(shí)距離；天線的波束中心指向(Xc,rb)處；在(sa,rb+sr)處有一個(gè)點(diǎn)目標(biāo)；飛機(jī)位置為x。當(dāng)Xc=0時(shí)，是理想的正側(cè)視條帶SAR模式，雷達(dá)天線波束中心保持與飛行航跡垂直，均勻照射地面。本文主要討論正側(cè)視情況。

雷達(dá)發(fā)射的線性調(diào)頻信號(hào)遇到點(diǎn)目標(biāo)后反射回雷達(dá)，SAR接收到的回波信號(hào)為：

圖1 SAR幾何關(guān)系示意圖Fig. 1 SAR geometry

其中：t為時(shí)間變量，c為光速，rs為點(diǎn)目標(biāo)與雷達(dá)之間的斜距，fn：第n個(gè)雷達(dá)脈沖發(fā)射信號(hào)的中心頻率，一般情況下每個(gè)脈沖發(fā)射的信號(hào)中心頻率不變均為f0,kn：第n個(gè)雷達(dá)脈沖發(fā)射信號(hào)的調(diào)頻斜率，一般情況下為定值k。

接收端的本振信號(hào)為：

其中2rb/c為回波起始采樣時(shí)間，在SAR系統(tǒng)中為一個(gè)固定值，一般取為測(cè)繪帶起始點(diǎn)或中心點(diǎn)。

將回波信號(hào)與參考函數(shù)混頻后得到：

若令TsFi=t-(2rb/c) ，其中-F/2+1≤Fi≤F/2,Ts為采樣間隔，F(xiàn)i為距離向采樣點(diǎn)序列編號(hào)，F(xiàn)為距離向采樣總點(diǎn)數(shù)，則SAR原始回波為：

3 脈沖壓縮體制條帶SAR的OSA成像方法

脈沖壓縮方式有兩種，分別為相關(guān)法與去斜法。在去斜體制下，根據(jù)SAR回波模型，利用去斜法得到距離向信息。而在非去斜接收的機(jī)載SAR系統(tǒng)中，發(fā)射較寬線性調(diào)頻信號(hào)以獲取高分辨率，若在視頻數(shù)字域采用去斜法會(huì)在頻域產(chǎn)生混疊，因此一般采用相關(guān)法進(jìn)行脈沖壓縮。采用相關(guān)法實(shí)現(xiàn)脈沖壓縮時(shí)，需要把信號(hào)FFT變換到頻域，將信號(hào)頻域與含有2次共軛相位的頻域?yàn)V波器進(jìn)行復(fù)乘，再IFFT變換到時(shí)域。將脈沖壓縮與方位向OSA處理相結(jié)合，SAR成像流程具體可以表示為圖2。

圖2 OSA算法應(yīng)用于條帶SAR距離向脈沖壓縮的流程示意圖Fig. 2 Processing flow of OSA for range pulse compression strip SAR

OSA應(yīng)用于距離向脈沖壓縮的條帶SAR的具體流程為：首先將距離向信號(hào)變換為頻域，在頻域與距離向參考函數(shù)、方位向去斜函數(shù)復(fù)乘，接著在距離向進(jìn)行插值，然后將去斜后的方位向數(shù)據(jù)劃分為多個(gè)重疊的子孔徑，再對(duì)每個(gè)子孔徑的方位向進(jìn)行FFT變換即方位向第1級(jí)濾波粗處理。子孔內(nèi)的方位向處理完畢后，可以得到方位向的粗分辨率信息，利用該信息計(jì)算得到每個(gè)子孔徑的距離方位耦合補(bǔ)償因子，對(duì)每個(gè)子孔徑進(jìn)行距離方位耦合補(bǔ)償，補(bǔ)償后完成距離向IFFT，由此完成了每個(gè)子孔徑的距離向脈沖壓縮與方位向粗分辨率成像，恢復(fù)了每個(gè)子孔徑內(nèi)的距離方位信息，得到每個(gè)子孔徑的圖像。利用子孔徑的圖像信息計(jì)算得到方位向精補(bǔ)償因子，補(bǔ)償方位向冗余的1次、2次空變相位誤差，得到方位向精分辨率圖像。最后對(duì)方位向數(shù)據(jù)進(jìn)行輸出選擇，去掉重疊部分的圖像并進(jìn)行幅度調(diào)制調(diào)整。以下對(duì)具體實(shí)現(xiàn)步驟進(jìn)行詳細(xì)說(shuō)明。

3.1 距離向FFT與方位向去斜

對(duì)條帶SAR回波進(jìn)行距離向FFT，變換到距離向頻域，其頻域表達(dá)式為：

3.2 子孔徑劃分與方位向FFT

將方位向位置x寫成離散化形式，按照?qǐng)D3在方位向進(jìn)行子孔徑劃分。假設(shè)全孔徑總長(zhǎng)度為N點(diǎn)，劃分為S個(gè)子孔徑，每個(gè)子孔徑長(zhǎng)度為M點(diǎn)，相鄰子孔徑間重疊點(diǎn)數(shù)為 (M-Δ)個(gè)，則有：

圖3 子孔徑劃分示意圖Fig. 3 Aperture divided into subapertures

子孔徑劃分大小的因素與子孔徑重疊率、第1級(jí)粗分辨率有關(guān)。子孔徑的重疊率即M/Δ，為抑制信號(hào)旁瓣，一般取為2～3[1,2]。子孔徑長(zhǎng)度的選取還受到第1級(jí)粗分辨率的限制，具體詳見(jiàn)式(16)和式(29)。通過(guò)子孔徑劃分，數(shù)據(jù)變成了N個(gè)2維矩陣，構(gòu)成了一個(gè)關(guān)于i,si,m的3維數(shù)據(jù)矩陣，具體表示為：

其中si為子孔徑內(nèi)采樣點(diǎn)的索引，－S/2+1＜si＜S/2,m為子孔徑序列的索引，－M/2+1＜m＜M/2。

式(12)的第4項(xiàng)指數(shù)項(xiàng)是關(guān)于子孔徑內(nèi)序列m的函數(shù)，若2次相位小于具體表示為：

則可以忽略關(guān)于m的2次項(xiàng)，對(duì)變量m做FFT變換，F(xiàn)FT變換后得到子孔內(nèi)方位向粗分辨率圖像，具體表示為：

由式(16)可以看到，每個(gè)子孔徑的壓縮峰值位置在隨子孔徑si，距離向i遷移，為了能夠得到精聚焦圖像，需要對(duì)遷移量進(jìn)行限制，一般將這個(gè)遷移量限制在半個(gè)粗分辨率單元以內(nèi)，具體可以表示為：

3.3 距離向IFFT

根據(jù)方位向粗處理結(jié)果，可以求出關(guān)于sa的第1次估計(jì)值由此可以求出的第1次估計(jì)值通過(guò)這些估計(jì)值可以算得距離方位的耦合補(bǔ)償向量：

補(bǔ)償后對(duì)變量i做IFFT變換得到：

由此恢復(fù)出距離向信息，得到了每個(gè)子孔徑的圖像，其中距離向?yàn)楦叻直媛?、方位向?yàn)榇址直媛?。由?20)可以看到距離向的位置隨si發(fā)生遷移，若遷移量能夠滿足半個(gè)距離分辨單元，則遷移可忽略，具體表示為：

3.4 方位向精處理

根據(jù)每個(gè)子孔徑圖像，可以求出sr的估計(jì)值進(jìn)一步求出的第2次估計(jì)值利用第2次估計(jì)值算得第2次方位向?yàn)V波處理前所需的誤差補(bǔ)償向量，表示為：

式(25)為正側(cè)視條帶SAR模式下的圖像重建，其中方位向輸出為：

圖像的距離向分辨率由發(fā)射信號(hào)帶寬決定，圖像的方位向分辨率為：

3.5 方位向數(shù)據(jù)輸出選擇與幅度校正

式(25)是一個(gè)關(guān)于3維向量(u1,u2,v)的數(shù)學(xué)表達(dá)式，而雷達(dá)圖像實(shí)際上是個(gè)2維圖像，因此需要將式(24)的方位向進(jìn)行數(shù)據(jù)輸出選擇，即將方位向數(shù)據(jù)進(jìn)行重排。v對(duì)應(yīng)的是距離向，(u1,u2)對(duì)應(yīng)的是方位向，因此數(shù)據(jù)輸出選擇是針對(duì)(u1,u2)進(jìn)行的數(shù)據(jù)重排，數(shù)據(jù)重排的原則與子孔徑劃分方法一致。另外在方位向經(jīng)過(guò)了兩次FFT運(yùn)算，第1次FFT運(yùn)算即粗處理過(guò)程中引入了幅度調(diào)制，影響了第2次FFT運(yùn)算即精處理輸出的幅度，導(dǎo)致了精處理輸出后的數(shù)據(jù)有明暗“條紋”。為保證最終圖像幅度一致性，在方位向精處理完畢后需要對(duì)幅度進(jìn)行校正，幅度校正因子與sinc函數(shù)幅度一致。

4 算法性能分析

4.1 運(yùn)算量與存儲(chǔ)量分析

以一個(gè)大小為M×N點(diǎn)的SAR原始數(shù)據(jù)塊為例，其中距離向點(diǎn)數(shù)為M，取值為4096，方位向點(diǎn)數(shù)為N，取值為4096點(diǎn)。采用本方法進(jìn)行成像的運(yùn)算量如表1所示。其中長(zhǎng)度為N的FFT或IFFT的浮點(diǎn)運(yùn)算(Floating Point Operations, FLOP)為5Nlog2(N)，一次復(fù)乘需要6個(gè)FLOP。OSA的總運(yùn)算量與RD, CS的5.61 G, 4.05 G[12]相比，分別下降了30%與5%。在實(shí)際系統(tǒng)處理時(shí)，一般脈沖重復(fù)頻率比是方位向多普勒帶寬的2倍以上，這樣在第1級(jí)子孔徑處理完畢后冗余點(diǎn)會(huì)達(dá)到子孔徑的一半，去掉冗余點(diǎn)后，OSA運(yùn)算量還會(huì)進(jìn)一步降低。

由于OSA將方位向數(shù)據(jù)劃分多個(gè)小數(shù)據(jù)塊，每次僅針對(duì)小數(shù)據(jù)塊進(jìn)行操作，因此OSA應(yīng)用于高速實(shí)時(shí)并行處理時(shí)對(duì)數(shù)據(jù)存儲(chǔ)量的要求方面具有很大的優(yōu)勢(shì)。對(duì)一個(gè)16384×16384點(diǎn)的數(shù)據(jù)塊進(jìn)行實(shí)時(shí)成像時(shí)，降4倍降采樣濾波，成像中距離向有效點(diǎn)輸出為12000點(diǎn)，圖像采用16 bit量化。方位向劃分S=256個(gè)子孔徑，每個(gè)子孔徑長(zhǎng)度為M=32,Δ=16點(diǎn)，第1級(jí)子孔徑處理完畢后去掉冗余點(diǎn)保留16點(diǎn)。分別采用RD, CS, OSA成像時(shí)需要的存儲(chǔ)量如表2所示。

由表2中可以看到，相對(duì)比CS方法，OSA對(duì)轉(zhuǎn)置存儲(chǔ)容量要求降低一倍以上；相比RD算法，轉(zhuǎn)置存儲(chǔ)量基本相當(dāng)，但是從運(yùn)算量、性能上要優(yōu)于RD。由于進(jìn)行了子孔徑劃分，第1次和第2次轉(zhuǎn)置存儲(chǔ)量很小，在實(shí)時(shí)處理時(shí)，利用高速數(shù)字信號(hào)處理芯片內(nèi)部的存儲(chǔ)空間就能夠?qū)崿F(xiàn)轉(zhuǎn)置，降低了對(duì)硬件資源的要求，還便于運(yùn)算與數(shù)據(jù)讀取的流水設(shè)計(jì)。因此在高分辨率處理時(shí)，從處理性能、運(yùn)算量和存儲(chǔ)量幾個(gè)方面考慮，OSA是一種比較優(yōu)化的實(shí)時(shí)處理方法。

表1 脈沖壓縮條帶SAR的OSA運(yùn)算量估計(jì)Tab. 1 Computation load of OSA for pulse compression strip SAR

表2 不同實(shí)時(shí)成像算法的存儲(chǔ)量對(duì)比Tab. 2 Data storage of different real-time imaging algorithms

4.2 成像范圍分析

式(13)，式(18)，式(21)與式(24)為本文算法對(duì)成像范圍的限制，根據(jù)計(jì)算比較分析，得到最終對(duì)方位向Dx，距離向成像范圍Dy與方位向第1級(jí)粗分辨率單元的成像限制，分別為：

圖4為本文算法中成像范圍與圖像分辨率的關(guān)系，在Ku波段時(shí)，若要得到0.3 m分辨率，方位向成像尺寸為1.2 km，距離向尺寸大于8 km，已經(jīng)能夠滿足大部分實(shí)際機(jī)載SAR系統(tǒng)實(shí)時(shí)處理的要求。通過(guò)成像范圍與方位向第1級(jí)粗分辨率限制的分析，可以指導(dǎo)子孔徑的具體劃分。

4.3 點(diǎn)目標(biāo)與原始數(shù)據(jù)成像分析

仿真中共放置21個(gè)點(diǎn)目標(biāo)，在場(chǎng)景中心有1個(gè)點(diǎn)目標(biāo)，其余的20個(gè)點(diǎn)目標(biāo)圍繞中心點(diǎn)目標(biāo)組成一個(gè)矩形場(chǎng)景，相鄰點(diǎn)方位向間距為100 m，距離向間距為400 m，水平方向?yàn)榉轿幌?，垂直方向?yàn)榫嚯x向。用本文算法重建圖像后，不僅圖像中心處的點(diǎn)聚焦良好圖像沒(méi)有形變，而且圖像邊緣處的點(diǎn)聚焦效果也滿足分辨率要求，具體如圖5(a)，圖5(b)所示。SAR系統(tǒng)參數(shù)與成像參數(shù)見(jiàn)表3。方位向子孔徑劃分時(shí)，重疊率選取為2，根據(jù)式(17)，式(30)計(jì)算得到子孔徑長(zhǎng)度需要大于12，在實(shí)際的SAR實(shí)時(shí)處理系統(tǒng)中，一般FFT的長(zhǎng)度取為2n，結(jié)合實(shí)時(shí)處理效率考慮，子孔徑長(zhǎng)度選取為64，子孔徑個(gè)數(shù)選為512個(gè)。選取了場(chǎng)景中心與邊緣共5個(gè)點(diǎn)進(jìn)行了成像評(píng)估，經(jīng)測(cè)量所有點(diǎn)的3 dB寬度、積分旁瓣比與峰值旁瓣比均達(dá)到指標(biāo)要求。其中點(diǎn)目標(biāo)1，點(diǎn)目標(biāo)2，點(diǎn)目標(biāo)4與點(diǎn)目標(biāo)5為圖像最邊緣處的4個(gè)點(diǎn)，分別位于圖像左上、右上、左下與右下角，點(diǎn)目標(biāo)3為場(chǎng)景中心處的點(diǎn)目標(biāo)。

圖4 OSA算法的成像范圍Fig. 4 Scene diameter limits for OSA

圖5 點(diǎn)目標(biāo)仿真結(jié)果Fig. 5 Point-targets simulation

圖6 是某無(wú)人機(jī)載SAR原始數(shù)據(jù)利用高速實(shí)時(shí)專用信號(hào)處理平臺(tái)的實(shí)時(shí)成像結(jié)果。SAR系統(tǒng)參數(shù)與成像參數(shù)與點(diǎn)目標(biāo)仿真參數(shù)相同，具體如表3所示。通過(guò)對(duì)實(shí)測(cè)數(shù)據(jù)的成像結(jié)果表明，實(shí)時(shí)成像處理速度滿足系統(tǒng)實(shí)時(shí)性要求，成像處理結(jié)果滿足系統(tǒng)對(duì)圖像質(zhì)量的要求。本文描述的適用于距離向脈沖壓縮條帶SAR的OSA算法能夠應(yīng)用于中小型無(wú)人機(jī)載SAR系統(tǒng)的實(shí)時(shí)成像處理器。

本文算法在時(shí)域?qū)⒎轿幌驍?shù)據(jù)劃分為多個(gè)孔徑，更適用于無(wú)人機(jī)載SAR系統(tǒng)合成孔徑時(shí)間長(zhǎng)的使用條件；本文算法對(duì)每個(gè)子孔徑數(shù)據(jù)分別處理，有利于并行實(shí)現(xiàn)、降低了數(shù)據(jù)存儲(chǔ)與實(shí)時(shí)信號(hào)處理硬件的設(shè)備要求同時(shí)還便于運(yùn)動(dòng)補(bǔ)償[13]，有助于無(wú)人機(jī)載SAR系統(tǒng)的小型化與輕量化；本文算法中相位誤差補(bǔ)償與成像處理同時(shí)完成，直接生成高分辨率圖像，具有較高的計(jì)算效率；本文算法在脈沖壓縮體制下工作，無(wú)需在接收端對(duì)回波信號(hào)進(jìn)行去斜處理，降低接收端設(shè)備的復(fù)雜度。但是本文算法也具有一定的局限性，在一定的成像范圍內(nèi)才有效，尤其是在0.1 m及以上高分辨率模式時(shí)，對(duì)成像范圍要求較為苛刻，因此為提高高分辨率SAR的成像范圍可以進(jìn)一步考慮多級(jí)子孔徑的成像方法。

表3 SAR系統(tǒng)與成像參數(shù)Tab. 3 Parameters of SAR system and imaging

圖6 某無(wú)人機(jī)載SAR數(shù)據(jù)實(shí)時(shí)成像結(jié)果Fig. 6 UAV SAR real-time image

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

本文提出了一種面向距離向脈沖壓縮的機(jī)載條帶SAR重疊子孔徑實(shí)時(shí)成像算法，該算法具有以下特點(diǎn)：(1)雷達(dá)數(shù)據(jù)經(jīng)過(guò)子孔徑劃分后，在成像過(guò)程中可以通過(guò)在多普勒域復(fù)乘校正因子消除每個(gè)子孔徑的距離向與方位向耦合，降低了運(yùn)算復(fù)雜度；(2)通過(guò)擴(kuò)展SAR回波模型的泰勒級(jí)數(shù)展開項(xiàng)數(shù)，通過(guò)子孔徑劃分、距離遷移校正因子、精聚焦因子等消除方位向高階相位誤差。與常規(guī)接收端去斜聚束SAR系統(tǒng)中的極坐標(biāo)OSA算法相比，本文算法有如下優(yōu)勢(shì)：(1)可適用于脈沖壓縮體制的SAR系統(tǒng)，獲取更高測(cè)繪帶寬的雷達(dá)圖像。(2)可適用于條帶式SAR工作模式。(3) 可降低實(shí)時(shí)成像處理中數(shù)據(jù)存儲(chǔ)量與運(yùn)算量要求，提高了系統(tǒng)硬件并行性。通過(guò)對(duì)模型分析以及仿真和實(shí)測(cè)數(shù)據(jù)的處理驗(yàn)證了該算法的有效性，并在中小型無(wú)人機(jī)載高分辨率SAR實(shí)時(shí)成像領(lǐng)域具有廣泛的應(yīng)用前景。

附錄 接收端去斜的SAR回波模型與OSA成像方法流程

SAR發(fā)射一個(gè)線性調(diào)頻信號(hào)，接收到的回波信號(hào)為：

其中t為時(shí)間變量，c為光速，fn：第n個(gè)雷達(dá)脈沖發(fā)射信號(hào)的中心頻率，kn：第n個(gè)雷達(dá)脈沖發(fā)射信號(hào)的調(diào)頻斜率。

在接收端的本振信號(hào)為：

則采樣后的雷達(dá)回波視頻信號(hào)最終形式為：

對(duì)式(A-9)中的距離向進(jìn)行FFT，方位向進(jìn)行進(jìn)行OSA處理，就能夠得到雷達(dá)圖像。

接收端采用一個(gè)與發(fā)射信號(hào)符號(hào)相反的本振、調(diào)頻斜率相反的混頻信號(hào)與SAR回波信號(hào)進(jìn)行混頻，在高頻段完成去斜處理。在Lynx系統(tǒng)中不僅要發(fā)射中心頻率、調(diào)頻斜率按規(guī)律變化的線性調(diào)頻信號(hào)，還要在接收端利用混頻信號(hào)發(fā)生設(shè)備，生成中心頻率、調(diào)頻斜率實(shí)時(shí)變化的本振信號(hào)，系統(tǒng)的起始采樣時(shí)間也實(shí)時(shí)發(fā)生變化。利用這些措施，在接收端完成了距離向、方位向的2維去斜，完成了距離向插值操作，有利于OSA的實(shí)時(shí)實(shí)現(xiàn)、簡(jiǎn)化了后端實(shí)時(shí)處理。圖7為基于接收端去斜的OSA成像方法流程圖。

圖7 基于接收端去斜的OSA成像方法流程圖Fig. 7 Processing flow of OSA for dechirp-on-reciever SAR

[1]Burns B L and Cordaro J T. SAR image formation algorithm that compensates for the spatially variant effects of antenna motion[J].Proceedings of SPIE, 1994, 2230(4): 14-24.

[2]Burns B L and Cordaro J T. Imaging synthetic aperture radar[P]. United States Patent, 1997, No.5608404.

[3]Doerry A W. Synthetic aperture radar processing with tiered subapertures[J].Communications and Radar, 1994, 30(4): 1125-1129.

[4]Doerry A W. Synthetic aperture radar processing with polar formatted subapertures[C]. Proceedings of 28th Asilomar Conferrence Signals System Computer, Pacific Grove, CA, 1994: 1210-1215.

[5]Doerry A W. Wavefront curvature limitations and compensation to polar format processing for synthetic aperture radar images[R]. Sandia National Labs, New Mexico, CA, Technical Report, SAND2007-0046, 2007.

[6]Walker B, Sander G, Thompson M,et al.. A highresolution, four-band SAR Testbed with real-time image formation[C]. International Geoscience and Remote Sensing Symposium(IGARSS) 1996, Lincoln, NE, 1996(3): 1881-1885.

[7]謝冬冬, 禹衛(wèi)東, 徐峰, 等. 利用OSA算法處理?xiàng)l帶SAR數(shù)據(jù)[J]. 系統(tǒng)工程與電子技術(shù), 2005, 27(6): 1003-1006. Xie Dong-dong, Yu Wei-dong, Xu Feng,et al.. OSA algorithm for the processing of stripmap SAR data[J].Systems Engineering and Electronics, 2005, 27(6): 1003-1006.

[8]Tang Yu, Zhang Bo, Xing Meng-dao,et al.. Azimuth overlapped subaperture algorithm in frequency domain for highly squinted synthetic aperture radar[J].IEEE Geoscience and Remote Sensing Letters, 2013, 10(4): 692-696.

[9]唐禹, 邢孟道, 保錚, 等. 基于重疊子孔徑極坐標(biāo)算法的波前彎曲效應(yīng)的補(bǔ)償[J]. 電子學(xué)報(bào), 2008, 36(6): 1108-1113. Tang Yu, Xing Meng-dao, Bao Zheng,et al.. Wavefront curvature compensation based on overlapped subaperture polar format algorithm[J].Acta Electronica Sinica, 2008, 36(6): 1108-1113.

[10]毛新華, 曹海洋, 朱岱寅, 等. 基于先驗(yàn)知識(shí)的SAR兩維自聚焦算法[J]. 電子學(xué)報(bào), 2013, 41(6): 1041-1047. Mao Xin-hua, Cao Hai-yang, Zhu Dai-yin,et al.. Prior knowledge aided two dimensional autofocus approach for synthetic aperture radar[J].Acta Electronica Sinica, 2013, 41(6): 1041-1047.

[11]Mao Xin-hua, Zhu Dai-yin, and Zhu Zhao-da. Polar formatalgorithm wavefront curvature compensation under arbitrary radar flight path[J].IEEE Geoscience and Remote Sensing Letters, 2012, 9(3): 526-530.

[12]Cumming I G and Wong F H. Digital Processing of Synthetic Aperture Radar Data[M]. Norwood, MA: Artech House, Inc., 2005, Chapter 11.

[13]田雪, 梁興東, 李焱磊, 等. 基于子孔徑包絡(luò)誤差校正的SAR高精度運(yùn)動(dòng)補(bǔ)償方法[J]. 雷達(dá)學(xué)報(bào), 2014, 3(4): 583-590.Tian Xue, Liang Xing-dong, Li Yan-lei,et al.. High precision motion compensation method based on the subaperture envelope error correction for SAR[J].Journal ofRadars, 2014, 3(4): 583-590.

Research on Overlapped Subaperture Real-time Imaing Algorithm for Pulse Compression Airborne Strip SAR System

Zhan Xue-li Wang Yan-fei Wang Chao Liu Bi-dan
(Institute of Electronics, Chinese Academy of Sciences,Beijing100190,China)

A high-resolution real-time subaperture imaing formation for range direction pulse compression airborne strip Synthetic Aperture Radar (SAR) system is presented. It can be used in no dechirp strip SAR system. By pulse compression in range direction and Overlapped Subaperture Algorithm (OSA) in azimuth direction this algorithm can compensate the range-azimuth cross error and the space variant phase error. In this study, first strip SAR geometry is analyzed and the strip SAR model is derived, and then the processing flow of OSA for range pulse compression strip SAR system is described in detail. Computation load, data storage and limitations of patch are analyzed then. The point-target simulations and live data processing results show the proposed approach is feasible and effective.

Synthetic Aperture Radar (SAR); Real-time imaging; Pulse compression; Overlapped Subaperture Algorithm (OSA)

TN958

A

2095-283X(2015)02-0199-10

10.12000/JR14126

詹學(xué)麗, 王巖飛, 王超, 等. 一種基于脈沖壓縮的機(jī)載條帶SAR重疊子孔徑實(shí)時(shí)成像算法[J]. 雷達(dá)學(xué)報(bào), 2015, 4(2): 199-208. http://dx.doi.org/10.12000/JR14126.

Reference format:Zhan Xue-li, Wang Yan-fei, Wang Chao,et al.. Research on overlapped subaperture realtime imaing algorithm for pulse compression airborne strip SAR system[J].Journal of Radars, 2015, 4(2): 199-208. http://dx.doi.org/10.12000/JR14126.

詹學(xué)麗(1976-)，女，中國(guó)科學(xué)院電子學(xué)研究所，博士，主要研究方向?yàn)镾AR成像理論及其信號(hào)處理等.

王巖飛(1963-)，男，中國(guó)科學(xué)院電子學(xué)研究所，博士生導(dǎo)師，研究員，主要研究方向?yàn)槲⒉ǔ上窭走_(dá)系統(tǒng)及其理論、數(shù)字信號(hào)處理等.

王 超(1983-)，男，中國(guó)科學(xué)院電子學(xué)研究所，博士，主要研究方向?yàn)镾AR成像理論及其信號(hào)處理等.

劉碧丹(1980-)，女，中國(guó)科學(xué)院電子學(xué)研究所，助理研究員，主要研究方向?yàn)镾AR成像理論及其信號(hào)處理等.

2014-11-05收到，2015-01-19改回；2015-04-07網(wǎng)絡(luò)優(yōu)先出版

國(guó)家自然科學(xué)基金 (61101201)資助課題

*通信作者： 詹學(xué)麗 16087415@qq.com