余小龍，張 炎，陶小輝，姜力暉，彭國良，曹 銳

（1.中國電子科技集團公司第三十八研究所，安徽合肥 230088；2.孔徑陣列與空間探測安徽省實驗室，安徽合肥 230088）

0 引 言

采用SAR 技術(shù)，將收發(fā)天線沿著長線陣的軌跡等速移動并輻射相參信號，從而實現(xiàn)等同于長線陣的方位向高分辨率［1］。線性合成孔徑陣列可以通過條帶模式進行數(shù)據(jù)采集［2］，這種模式下隨著雷達中心頻率的增加至毫米級，SAR 系統(tǒng)獲得的方位向分辨率能力也隨之更高。近年來，基于CMOS 的調(diào)頻連續(xù)波雷達能夠?qū)⒛M單元與射頻單元集成到一個緊湊芯片中［3］，這使超寬帶毫米波雷達系統(tǒng)小型化并低成本化，被廣泛用于健康監(jiān)測、武器藏匿檢查、建筑結(jié)構(gòu)檢查等視覺失靈領(lǐng)域的高分辨率近場成像［4-9］。傳統(tǒng)的SAR成像系統(tǒng)采用由單個陣元構(gòu)成的密集陣列，而MIMO陣列則由多個陣元構(gòu)成的子陣進行信號收發(fā)，實現(xiàn)了數(shù)據(jù)采集時間的減少。

對于單個陣元，現(xiàn)階段長期使用的SAR 全息成像方法是基于RMA 的波數(shù)域方法［10］，為了提高計算速度常常需要在時間-空間域、頻率-波數(shù)域通過Stolt 插值法得到均勻采樣的后向散射信號數(shù)據(jù)后進行快速傅里葉變換（Fast Fourier Transform,FFT）［11］，然而這種插值過程繁瑣，成像效率較低。本文討論了引入NUFFT 代替Stolt 插值的方法，極大提高了RMA 重建圖像速度［12-16］，這對于相比單個陣元復(fù)雜度更高的MIMO 的采用有著重要意義［17］。

在實際應(yīng)用中，MIMO 陣列不同陣元拾取的后向散射信號之間存在著相位或者距離不匹配［18］，同時多站雷達的收發(fā)天線之間存在耦合的情況［19］，未經(jīng)過校準(zhǔn)的SAR 成像存在嚴(yán)重的偽影。本文提出了通過對平行于平面MIMO-SAR 陣列的金屬平面板照射后拾取的散射信號進行分析，實現(xiàn)192個虛擬子陣列陣元的同時校準(zhǔn)。

1 MIMO-SAR成像

1.1 虛擬陣列拓?fù)浣Y(jié)構(gòu)

如圖1所示，將射頻板與數(shù)據(jù)采集板的級聯(lián)板安裝在由兩個絲桿步進電機組成的平面掃描平臺上。通過絲桿步進電機的帶動級聯(lián)板橫向掃描和斜向復(fù)位，實現(xiàn)了收發(fā)天線陣列在整個平面內(nèi)的位置移動。

圖1 級聯(lián)板掃描機械運動結(jié)構(gòu)

虛擬陣元位置是由TX 天線和RX 天線在空間上位置卷積產(chǎn)生［20］。在近場條件下，若天線空間位置rT和rR詳細表示為(xt,yt,zt)和(xr,yr,zr)，目標(biāo)位置表示為(x′,y′,z′) ，收發(fā)天線空間位置中點到目標(biāo)點距離為R，則我們可以將新的虛擬陣元與目標(biāo)之間的往返距離表示為［21］

如圖2所示，發(fā)射天線TX1和接收天線RX1可以構(gòu)成虛擬天線陣列的一個陣元CH1，其位置由TX1 和RX1 的空間位置卷積產(chǎn)生，所有的2 個TX天線和4個RX天線可以構(gòu)成8個虛擬陣元位置。

圖2 MIMO拓?fù)浣Y(jié)構(gòu)

如圖3所示，如果對天線陣列水平和豎直位移進行合理的調(diào)整，則可以獲得一個更大的虛擬陣列。

圖3 子陣列構(gòu)成200 mm×200 mm掃描孔徑陣列

1.2 基于NUFFT改良的RMA

經(jīng)典的SAR 數(shù)據(jù)成像主要有基于時域的后向投影算法［22］（Back-Projection Algorithm, BPA）和基于波數(shù)域的RMA［23］，兩者在SAR 成像聚焦原理上具有一致性，但在成像過程中RMA 的速度要明顯快于BPA［24］，因此本文選擇以RMA為基礎(chǔ)，進行圖像復(fù)現(xiàn)。

可以進行如下總結(jié)：

基于NUFFT的波數(shù)域三維全息圖像重建算法：

1）對回波信號建立波數(shù)域上的模型s(u,v,k,h)；

2）對u和v在波數(shù)域上進行3D-FFT，得到S(ku,kv,k,kz)；

4）在ky通過Stolt 插值使得后面的FFT 運算得以實現(xiàn)/在方位向上使用NUFFT代替插值；

5）通過對(ku,kv)維上的降維積累，實現(xiàn)4D 數(shù)據(jù)降維為3D；

7）實現(xiàn)圖像重建。

2 基于反射面的MIMO陣列校準(zhǔn)

MIMO 陣列成像系統(tǒng)的實際使用中，各個通道之間存在時間延時、通道響應(yīng)、系統(tǒng)噪聲等導(dǎo)致的誤差。本文對各個誤差的產(chǎn)生原因不作更為細致的分析，對客觀存在的不同通道之間的系統(tǒng)誤差值進行計算，從而實現(xiàn)MIMO陣列校準(zhǔn)［28-29］。

同時，對于多站雷達采集信號進行range-FFT分析，結(jié)果如圖4所示。可以更清晰地看到除了在理論位置附近獲得了一個幅值尖峰，在離天線較近位置處也有一個由收發(fā)天線耦合造成的幅值尖峰，本文在下述校準(zhǔn)過程中對互耦的影響也進行了補償。

圖4 點目標(biāo)回波信號的range-FFT結(jié)果

2.1 基于金屬反射面的MIMO陣列校準(zhǔn)方法

本文使用平面金屬反射面作為成像校準(zhǔn)目標(biāo)進行實驗，其特點在于當(dāng)虛擬陣列平面和收發(fā)天線陣列相平行且金屬平面足夠大時，虛擬陣元均落于金屬反射面投影內(nèi)，可以認(rèn)為每個陣元與目標(biāo)始終保持相等的距離。此時通道之間的差異全部來自于系統(tǒng)誤差。

如圖5所示，目標(biāo)物體平面金屬反射面平行虛擬陣列平面放置，距離為z0，則對于每一個虛擬陣元距離目標(biāo)的距離均為z0。本文假定虛擬陣列平面上有一組額外的收發(fā)天線對構(gòu)成的虛擬陣元，理想條件下，忽略信號振幅和殘余視頻相位項，采集到的后向散射信號為

圖5 利用平行金屬平面進行采樣

式中，fb=Kτ為拍頻信號的頻率，?= 2πf0τ為拍頻信號的相位，可以看出均與信號往返時延τ有關(guān)，式中τ= 2z0/c。假設(shè)對拍頻信號在時域上均勻采樣，采樣周期為Ts，采樣頻率為N，則可以由s(t)連續(xù)信號得到離散信號s[n]=s(nTs)，式 中n=[0,1,2,…,N- 1]。通過對采樣拍頻信號進行N點離散傅里葉變換（DTFT）得到

式中，ωb= 2πfb= 4πKz0/c為拍頻信號角頻率。由于拍頻信號是在時間上的均勻采樣，上式的采樣版本可以通過range-FFT獲得。

將上述生成的理想后向散射信號s(t)作為參考信號，將實際采集到的信號與其進行對比，從而得到對應(yīng)通道的校準(zhǔn)參數(shù)。對于某個MIMO 虛擬陣元采集到的未校準(zhǔn)信號，存在收發(fā)增益al、系統(tǒng)時延（包括儀器帶來的系統(tǒng)延遲τi、收發(fā)天線耦合導(dǎo)致的時延τ）l，忽略熱噪聲和殘余視頻相位項，則未校準(zhǔn)信號可表示為

式中，s*(t)為s(t)的復(fù)共軛。從而將估算ηi與fi轉(zhuǎn)化為估算噪聲中的單頻

式中Wl(f)表示為

如果ωl(t)是在時域上的均勻采樣，則可以通過FFT 得到Wl(f)的離散型。最后對于復(fù)增益系數(shù)ηi可以通過對中插入估算得到的來計算。

2.2 192個虛擬通道的MIMO陣列校準(zhǔn)

在本文實驗中，采用雷達天線的橫向掃描構(gòu)建虛擬陣列，對每個通道的后向散射信號而言，信號相位變化有

式中，k為波數(shù)，D為虛擬陣元與目標(biāo)的往返距離，Δd=D/2 為單程距離，λ為波長，f為頻率，c為光速，考慮到對于平行金屬反射面照射，在各個陣元位置上Δd值不變，因此信號相位應(yīng)該為常數(shù)φ=4πfΔd/c。對一次單行橫向掃描采樣結(jié)果進行相位分析，選取其中1 個通道進行相位展開，得到相位信息如圖6 所示，x軸為采樣幀序列，y軸為采集信號的相位值?？梢郧逦乜吹?，在整個掃描過程中，相位值相當(dāng)穩(wěn)定，整體在2 rad 以內(nèi)，這種情況下對于金屬平面的掃描結(jié)果可以等效于每一個通道始終對于距離不變的目標(biāo)進行掃描。

圖6 單對天線在不同陣元采集信號相位角變化

經(jīng)過對本文裝置的分析，認(rèn)為圖像最左段與最右段出現(xiàn)的相位抖動情況，來自于掃描天線在加速移動中帶來的掃描架振動；而在圖像中段出現(xiàn)絕對值小于1.5π 的相位變化，來自于實驗金屬平面本身在加工中不夠光滑導(dǎo)致。這種誤差在對于使用點目標(biāo)頂點校準(zhǔn)時也同時產(chǎn)生影響。

同時，如果對于點目標(biāo)進行照射，對于理想的信號相位角應(yīng)當(dāng)是一條拋物線，但是實際可以看到信號相位角變化過程中有明顯的振動，如圖7所示，這正是掃描架機械振動帶來的影響，因此對于這種情況需要進行函數(shù)擬合使得曲線更為平滑。

圖7 機械振動帶來的信號相位的微小變化

如圖8所示，本文對圖4中的實驗抽取了20個通道的結(jié)果進行了擬合，擬合結(jié)果和預(yù)想的平面結(jié)果基本吻合。

圖8 部分陣元散射信號相位采集

根據(jù)對金屬平面的掃描得到每個通道采集信號的距離和相位信息進行分析，得到了192個通道不匹配的具體數(shù)值。這192 個通道中的144 個負(fù)責(zé)距離向成像，48 個通道負(fù)責(zé)方位向成像。根據(jù)2.1節(jié)對于不同通道之間偏差的計算可以得到如圖9 所示的距離偏差和相位偏差結(jié)果，通過乘以加權(quán)值矩陣補償偏差至同一基線，即實現(xiàn)了多通道相位信息校準(zhǔn)。

圖9 144通道在計算上的偏差

2.3 在金屬平面傾斜情況下的距離估算

在2.1 節(jié)中進行多通道校準(zhǔn)時，需要金屬平面和虛擬陣元二維平面平行，而在實際操作放置實驗裝置時這一理想條件很難實現(xiàn)。本節(jié)就金屬平面傾斜放置的情況下，對多通道校準(zhǔn)系數(shù)產(chǎn)生的影響進行討論。

如圖10 所示，在笛卡兒坐標(biāo)系中，平面MIMO陣列位于xOy面上，目標(biāo)平面相對MIMO 陣列傾斜放置。假設(shè)平面無限大，平面方程為aX+bY+cZ+d= 0；選取虛擬陣列上任一陣元，位置為(m,n,0)，則經(jīng)過該虛擬陣元且與目標(biāo)平面垂直的直線方程為；假設(shè)虛擬陣元在目標(biāo)平面上的投影為(m′,n′,l′)，則知(m′,n′,l′)滿足

圖10 利用傾斜金屬平面進行采樣

求解得到

則此時陣元距離反射平面距離從2.1節(jié)的z0更新為

3 MIMO-SAR成像系統(tǒng)搭建

本文構(gòu)建了一套基于W 波段的毫米波調(diào)頻連續(xù)波雷達成像實驗平臺，以驗證本文所提成像算法的重構(gòu)效果與校準(zhǔn)算法的有效性。

3.1 系統(tǒng)硬件組成

本文采用基于平面陣列的SAR 成像，考慮到數(shù)據(jù)采樣后處理的簡便性，雷達平臺在橫縱方向的運動均采用逐行光柵掃描，即在橫向完成一行的掃描后移回初始位置，同時縱移一定距離，然后在新的縱向位置上完成橫向單行掃描。雷達運動掃描平臺的控制、射頻板雷達信號的發(fā)送、數(shù)據(jù)采集板的后向散射數(shù)據(jù)接收，均由上位機控制。設(shè)計成像平臺硬件架構(gòu)如圖11所示。

圖11 設(shè)計成像平臺硬件架構(gòu)

本文中上位機對射頻板和數(shù)據(jù)采集板級聯(lián)而成的信號收發(fā)模塊發(fā)出指令，采集到的雷達數(shù)據(jù)通過以太網(wǎng)傳輸回上位機以進行后處理。采用兩根絲桿步進電機垂直交叉構(gòu)成一個二維運動結(jié)構(gòu)，并將水平軸步進電機安裝在豎直軸電機上。安裝完成的成像平臺實物圖如圖12所示。

圖12 硬件平臺系統(tǒng)

3.2 系統(tǒng)軟件組成

考慮到雷達運動平臺和信號收發(fā)模塊在時間上的同步性，因此不能直接使用毫米波雷達成像的配套工具。毫米波MIMO-SAR 成像系統(tǒng)處理流程如圖13所示，包括雷達運動與數(shù)據(jù)采集模塊、校準(zhǔn)模塊和圖像重構(gòu)模塊。本文設(shè)計了一個可視化用戶操作界面，如圖14所示，集合了雷達運動和數(shù)據(jù)采集、保存的功能，實現(xiàn)了輸入運動參數(shù)和信號參數(shù)后對整個成像平臺的快速工作配置，并解決了運動模塊和信號收發(fā)模塊的在時間同步性上的問題，完成了信號采集與數(shù)據(jù)保存的一鍵化實現(xiàn)。

圖13 MIMO SAR成像系統(tǒng)處理流程

圖14 軟件操作界面

在可視化用戶操作界面中，可以直接修改平臺運動參數(shù)和信號參數(shù)，包括了雷達的水平和垂直初始位置、掃描模式選擇、橫向掃描距離、縱向移動距離、幀周期、數(shù)據(jù)保存路徑等配置。

4 實驗成像與分析

4.1 實驗參數(shù)設(shè)計

根據(jù)方位向理論分辨率與距離的關(guān)系以及Nyquist 采樣定律的要求，設(shè)置成像實驗采用的雷達參數(shù)如表1 所示，合成孔徑參數(shù)配置如表2 所示。上位機發(fā)出指令對MIMO-SAR 成像平臺基本參數(shù)進行配置，將雷達采集板中的原始中頻數(shù)據(jù)通過千兆網(wǎng)口傳輸回上位機存儲，再對實際采集的數(shù)據(jù)進行后處理得到可直接用于后續(xù)算法處理的數(shù)據(jù)立方體。

表1 成像實驗雷達參數(shù)

表2 合成孔徑參數(shù) mm

4.2 點目標(biāo)成像分析

對于已有成像系統(tǒng)實際應(yīng)用中的成像分辨率測算，可以利用點目標(biāo)擴展函數(shù)（Point Spread Function, PSF）對某一點目標(biāo)的解析能力進行描述。通過對某一理想點目標(biāo)進行成像，分析其PSF函數(shù)-3 dB 寬度，以此作為對成像分辨率的評價指標(biāo)。

本文實驗中，角反目標(biāo)與掃描平面距離為100 mm，掃描孔徑為200 mm×200 mm，橫向、縱向步進掃描間隔設(shè)置均為1 mm，具體仿真參數(shù)采用4.1 節(jié)參數(shù)。利用采集到的信號對點目標(biāo)進行圖像重建，得到如圖15 所示的成像結(jié)果，可以觀察到清晰明亮的單個點目標(biāo)，其平面坐標(biāo)位置與真實目標(biāo)位置相一致。以重建點目標(biāo)圖像為原點，沿水平方向和垂直方向分別有多個旁瓣，如圖16所示，對于能量較低的旁瓣，可以通過設(shè)置動態(tài)范圍將其濾除。理想點目標(biāo)成像結(jié)果的亮斑大小可定性說明成像系統(tǒng)的分辨率性能。對于理想點目標(biāo)，其空間體積可以忽略，根據(jù)調(diào)頻連續(xù)波雷達分辨率理論

圖15 點目標(biāo)成像結(jié)果圖

圖16 PSF的二維展示

得到計算理論分辨率為10 mm。而對于圖15 的成像結(jié)果，可以采用PSF對實際成像分辨率進行大致估算，如圖16 所示。將成像結(jié)果分別沿著水平軸方向和豎直軸方向做切面，對所得結(jié)果的幅值量進行歸一化處理并取dB 計算，如圖17 所示。觀察分析主瓣-3 dB 寬度范圍可知，本文雷達成像平臺的成像分辨率約為9.5 mm，這與前述所說的理論值基本符合。與理論計算值基本一致。

圖17 PSF切面展示

進一步觀察，可以看到水平方向相比垂直方向整體上分辨率一致，但波形明顯更為光滑，分析原因為水平方向運動掃描間隔更小。

4.3 MIMO-SAR成像

4.3.1 金屬尺成像

選擇金屬直尺作為成像目標(biāo)。鋼制直尺長度為120 mm，寬度為25 mm，豎直平行立于距離雷達發(fā)射平面前，目標(biāo)直尺上方有一個直徑7 mm 的圓形孔洞，方便驗證成像的準(zhǔn)確性。目標(biāo)中心位置與雷達掃描平面中心位置基本對齊，距離雷達掃描平面中心距離約100 mm，成像目標(biāo)在x軸方向的長度約為25 mm，y軸方向長度約為120 mm。設(shè)置水平和垂直兩個方向運動平臺掃描孔徑范圍大小均為200 mm。x軸方向雷達運動步進間隔為1 mm，y軸方向步進間隔為1 mm，采樣步進值滿足空間采樣定律理論值。掃描目標(biāo)如圖18所示。

圖18 鋼尺目標(biāo)實際圖像與尺寸

對于實驗所得數(shù)據(jù)分別直接利用1.2 節(jié)中改進的RMA 進行圖像重建和基于2.1 節(jié)進行多通道校準(zhǔn)后再進行圖像重建，兩次實驗的成像結(jié)果如圖19所示。

圖19 未校準(zhǔn)圖像與校準(zhǔn)后圖像對比

可以看出，在校準(zhǔn)前，MIMO-SAR 成像帶有嚴(yán)重的偽影和虛焦的情況，經(jīng)過校準(zhǔn)，兩種情況都得到了改善，并且鋼尺上方的圓形空洞有較清晰的成像結(jié)果。

值得注意的是，實驗如果采用0.9 mm 孔徑的單發(fā)單收毫米波雷達實現(xiàn)陣列孔徑為200 mm×200 mm SAR 成像，數(shù)據(jù)采集時間約為8 080 s；本文采用了MIMO-SAR，可以實現(xiàn)了86 個通道的同時掃描，則實現(xiàn)相同陣列孔徑SAR 成像理論機械掃描時間最低可降為原來的1/86，而對于200 mm×200 mm SAR 成像需要采集時間則為120 s，這顯著提高數(shù)據(jù)采集的速度。

4.3.2 鏤空鋁箔成像

將印有“CETC”和“KLAASA”鏤空字樣的鋁箔附著在泡沫板上，在距離射頻板約200 mm 進行掃描采樣。采用獲得的校準(zhǔn)數(shù)據(jù)對采樣結(jié)果進行校準(zhǔn)，得到的成圖效果如圖20所示。

圖20 對小尺寸鏤空圖像進行成像處理

成像結(jié)果相對理想，可以根據(jù)成像目標(biāo)的尺寸可以看到在對于6 mm 的小尺寸目標(biāo)成像有良好的效果，這驗證了在整個實驗裝置不移動的情況下，通過一次校準(zhǔn)數(shù)據(jù)獲得后，可以穩(wěn)定用于之后成像校準(zhǔn)計算中。

5 結(jié)束語

本文主要討論了在近場MIMO-SAR 成像中采用RMA 實現(xiàn)圖像重現(xiàn)并通過引入NUFFT 對其運算速度進行優(yōu)化；對于采用MIMO系統(tǒng)討論了利用照射金屬反射面實現(xiàn)多通道的快速校準(zhǔn)，有效減少了多基地雷達帶來的成像偽影。本文搭建了200 mm×200 mm 孔徑大小的SAR 平臺，通過校準(zhǔn)后的圖像良好地重現(xiàn)了原目標(biāo)尺寸，成像分辨率達到了2.5 mm。相較于原有的單通道SAR 成像，本文實現(xiàn)了86 個虛擬通道的SAR 成像，使得機械掃描時間理論最高可以縮短至原來的1/86，這對需要實現(xiàn)快速SAR成像的應(yīng)用場景有重要的意義。