張榆紅 邢孟道 徐 剛

1 引言

逆合成孔徑雷達（Inverse Synthetic Aperture Radar, ISAR）的觀測對象一般為非合作性目標(biāo)[1]，如空中目標(biāo)和空間目標(biāo)等，通過傳統(tǒng)的距離-多普勒成像算法可以獲得目標(biāo)在兩維平面上的2維圖像。然而對機動性很強的目標(biāo)而言，在不同時間不同姿態(tài)下目標(biāo)2維ISAR圖像變化很大，增加了目標(biāo)分類和識別的難度，這也促使對獲取目標(biāo)3維圖像的方法進行研究。

現(xiàn)有文獻提出了一些3維ISAR成像的方法[25]-。文獻[4]提出 InISAR（Interferometric ISAR, InISAR）成像技術(shù)，利用3個天線構(gòu)成垂直基線構(gòu)型，通過多通道ISAR圖像間的干涉相位進行干涉處理來實現(xiàn)對目標(biāo)形狀的估計[4]。然而為了得到高質(zhì)量的干涉相位，獲取聚焦良好的2維ISAR圖像或者提取精確的目標(biāo)散射中心是必要的前提，這在實際中具有較大難度。由于 ISAR/InISAR系統(tǒng)針對的非合作性目標(biāo)相對雷達的轉(zhuǎn)動軸心和速度具有時變性，機動性很強，會在回波信號中引入時變和空變的多普勒調(diào)制，然而傳統(tǒng)的自聚焦方法無法實現(xiàn)對該時變和空變的多普勒補償，使得成像的復(fù)雜性增加。同時，多功能 ISAR系統(tǒng)為了同時對多個目標(biāo)實現(xiàn)檢測、跟蹤和成像，雷達波束會在不同的波束來回進行切換，對每個目標(biāo)而言只能獲得稀疏孔徑觀測，這對傳統(tǒng)的ISAR成像提出了挑戰(zhàn)。針對上述問題，文獻[6]通過構(gòu)建 chirp-傅里葉字典，并利用稀疏化表征技術(shù)實現(xiàn)了對機動目標(biāo)稀疏孔徑下的高分辨成像。然而，包括文獻[6]等的現(xiàn)有算法都是基于單通道形式的，在信噪比比較低時，ISAR圖像的副瓣、柵瓣和噪聲水平會提高，不利于后續(xù)的3維成像。

本文將文獻[6]基于稀疏表征的機動目標(biāo) ISAR成像推廣到 InISAR多通道聯(lián)合稀疏表征的情況。根據(jù)多通道圖像的結(jié)構(gòu)稀疏特性將 InISAR成像問題轉(zhuǎn)換為聯(lián)合稀疏約束的求解問題。對于勻加速轉(zhuǎn)動的目標(biāo)，回波的多普勒調(diào)制可以建模成線性調(diào)頻的形式，并用 chirp-傅里葉字典來描述其機動性，并用改進的正交匹配追蹤（Orthogonal Matching Pursuit, OMP）算法對聯(lián)合稀疏約束最優(yōu)化問題進行求解，實現(xiàn)多通道 ISAR成像。接著利用各個通道估計的ISAR圖像和調(diào)頻參數(shù)實現(xiàn)對目標(biāo)的3維幾何重構(gòu)。由于基于多通道數(shù)據(jù)進行聯(lián)合處理，能夠?qū)崿F(xiàn)高質(zhì)量的目標(biāo)3維成像。最后，通過實測數(shù)據(jù)驗證了本文提出方法的有效性。

2 信號模型

InISAR系統(tǒng)一般采用3個天線構(gòu)成垂直基線構(gòu)型[7]，本文以XOY平面的干涉基線為例。如圖1所示，天線A和天線B構(gòu)成水平方向的干涉對，且基線長度為b，在觀測時間內(nèi)目標(biāo)從位置 1運動到位置2, φ （ tm）為目標(biāo)的瞬時方位角。假設(shè)天線A發(fā)射線性調(diào)頻（Linear Frequency Modulation, LFM）信號，分別對各個天線的回波信號做相干檢波和脈沖壓縮后各通道的信號模型為

圖 1 XOY平面干涉模型

將式（2）代入式（1）中，則可以得到機動目標(biāo)的全孔徑InISAR信號模型為

式（3）中， AA,p= Ap?ex p -j4 π fc（ R0+ yp）/c, AB,p=AA,p?ex p [ - j2π fc?b?φ0/c]。同理可對垂直方向干涉天線進行處理，得到機動目標(biāo)在通道C的全孔徑信號模型。

對式（3）的信號進行離散化，并根據(jù)壓縮感知理論[10]，考慮隨機缺失采樣和塊缺失采樣兩種稀疏孔徑采樣方式，如圖2所示，將其推廣到稀疏孔徑的情況：

圖 2 稀疏孔徑幾何模型

3 聯(lián)合稀疏約束InISAR 3維成像

3.1 多通道聯(lián)合稀疏約束成像模型

如式（3）所示，對每個距離單元的機動目標(biāo)回波信號模型可以建模為多分量chirp信號的形式。采用文獻[6]所提出的方法構(gòu)建chirp-傅里葉字典[11,12]，可以將式（5）進一步表示為

其中， xi表示第i個通道的 ISAR圖像， C Fsa表示chirp-傅里葉字典CF的稀疏孔徑采樣形式,

其中 W = e-j（2π/M）, κ =κ 和l=l分別表示第p

M

pp個散射點的離散頻率和調(diào)頻率。

為了避免獨立對各個通道進行處理時引起的干涉相位質(zhì)量的下降，確保各通道之間回波信號有較高的相干性，需要對多通道 chirp-傅里葉字典進行統(tǒng)一構(gòu)造。在本文中，對于空間相鄰分布的天線陣元，假設(shè)各個通道目標(biāo)散射系數(shù)相同[12]，將式（6）推廣到多通道聯(lián)合chirp-傅里葉字典的形式：

其中，X2,1表示所求矩陣X的聯(lián)合l2,1范數(shù)，ε表示噪聲水平。如文獻[13]所述，相比于l1范數(shù)約束，通過利用最小聯(lián)合l2,1范數(shù)約束，可以有效提高各個通道的稀疏表征能力，從而提高信號重構(gòu)的精度[14,15]。至此，多通道聯(lián)合ISAR成像問題轉(zhuǎn)化為對式（8）的最優(yōu)化求解問題。

3.2 改進OMP算法2維成像求解

本文提出了一種基于多通道數(shù)據(jù)聯(lián)合處理的OMP算法來解決聯(lián)合稀疏表征的問題，實現(xiàn)對式（8）的求解。具體步驟如下：

步驟 1 聯(lián)合多通道的調(diào)頻參數(shù)估計。該步驟的主要任務(wù)是估計 Ssa中的最大分量chirp信號的調(diào)頻參數(shù)。結(jié)合式（8）的聯(lián)合稀疏表征，進行聯(lián)合多通道調(diào)頻參數(shù)估計為其中， 0 ≤ κ1≤ M -1和 - M /2 + 1 ≤ l1≤ M /2。式（9）可以看作是一種聯(lián)合多通道數(shù)據(jù)的基追蹤求解過程。為了實現(xiàn)（CFsa）H的運算，首先對稀疏孔徑數(shù)據(jù)空缺的部分進行補零，然后基于補零后的數(shù)據(jù)利用 FFT可以實現(xiàn)改進的離散調(diào)頻傅里葉變換（Modified Discrete Chirp Fourier Transform,MDCFT）操作。需要注意的是，在該步驟中，為了提高調(diào)頻參數(shù)估計的精度，進行參數(shù)粗估計后，可以以該參數(shù)為中心截取部分圖像進行插值處理后再進行參數(shù)估計。最后利用估計的調(diào)頻參數(shù)進行聯(lián)合chirp-傅里葉字典基矩陣的構(gòu)造 Φ =CFl1sa,κ1，其中κ1和l1分別表示對最大 chirp分量信號估計的頻率和離散調(diào)頻率。

步驟 2 聯(lián)合多通道散射中心提取和觀測信號殘余分量估計。在這一步中，主要的目的是估計X的主要分量。由前面的假設(shè)可知，目標(biāo)的散射系數(shù)強度在所有的通道具有一致性，則圖像幅度在所有的通道相同，將ISAR復(fù)圖像分解為幅度和相位兩項，則式（7）可以轉(zhuǎn)換為

其中， F FT[·]表示 FFT操作。至此，多通道聯(lián)合成像算法的求解已經(jīng)完成，可以看出聯(lián)合多通道數(shù)據(jù)基于聯(lián)合稀疏表征的方式進行參數(shù)估計，能夠有效抑制稀疏孔徑觀測以及噪聲和雜波的影響，實現(xiàn)更高質(zhì)量的ISAR成像。

3.3 3維目標(biāo)幾何模型重構(gòu)

根據(jù)所得的多通道 InISAR圖像和調(diào)頻參數(shù)以本文提出的方法實現(xiàn)目標(biāo)3維幾何模型和轉(zhuǎn)動參數(shù)估計。將估計所得到的調(diào)頻參數(shù)（κp和lp）, InISAR圖像（ISAR 圖像 R D1, R D2和 R D3）和式（2）的目標(biāo)轉(zhuǎn)動信號模型作為該算法的輸入。首先，利用InISAR圖像的干涉相位對目標(biāo)的3維幾何位置進行粗略估計。接著，選擇強度較高的散射點進行樣本選擇，并利用選擇樣本的參數(shù)采用加權(quán)最小二乘估計（Weighted Least Square Errors, WLSE）的估計算法進行轉(zhuǎn)速估計：

本文提出的算法具體流程如圖3所示。

4 實驗驗證

下面對錄取的民用飛機目標(biāo)的實測數(shù)據(jù)進行InISAR 3維成像實驗。雷達以150 Hz的脈沖重復(fù)頻率發(fā)射中心頻率8.1 GHz，帶寬800 MHz的LFM信號。由系統(tǒng)測量數(shù)據(jù)可知，雷達與目標(biāo)之間的距離大約為29 km，基線近似為16 m。圖4（a）為該民航飛機的光學(xué)圖像。首先，對3個通道的數(shù)據(jù)進行統(tǒng)一的平動補償，包括包絡(luò)對齊和初相校正。接著，進行 InISAR圖像之間的像素配準(zhǔn)操作，剩余的由目標(biāo)轉(zhuǎn)動引入的距離維越距離單元徙動（Migration Through Resolution Cells, MTRC）利用文獻[16]中提出的算法進行校正。在全孔徑下，對距離MTRC校正后的信號利用傳統(tǒng)距離-多普勒算法進行2維成像處理[16]，其結(jié)果如圖 4（b）所示。由圖 4（b）可見，由于目標(biāo)的機動性，距離-多普勒算法的成像結(jié)果在方位向存在模糊。然后，考慮稀疏孔徑的情況，圖5（a）和圖 5（b）分別對應(yīng)隨機缺失采樣方式和塊缺失采樣方式下3/4稀疏孔徑采樣數(shù)據(jù)量?；谙∈杩讖綌?shù)據(jù)，利用本文提出的算法分別進行2維和3維成像處理。其中，圖6（a）和圖6（b）分別表示信噪比為5 dB和10 dB時的2維成像結(jié)果，且在圖6（a）和圖6（b）中，第1列和第2列均分別對應(yīng)傳統(tǒng)單通道算法（通道A）和本文聯(lián)合多通道算法的成像結(jié)果，第1行和第2行均分別對應(yīng)隨機缺失采樣和塊缺失采樣兩種方式下的成像結(jié)果，可見本文算法可以有效抑制稀疏孔徑的影響得到良好的聚焦成像，且在相同稀疏孔徑采樣數(shù)據(jù)量下，隨機缺失采樣方式比塊缺失采樣方式具有更好的成像結(jié)果。同時，比較圖 6（a）和圖 6（b）可以得出，隨著信噪比的降低，成像性能明顯下降。圖7為3維重構(gòu)結(jié)果，從其結(jié)果能夠看出完整的飛機外形，且機頭和機翼能夠明顯地區(qū)分開，對應(yīng)較好的3維重構(gòu)結(jié)果。需要指出的是，在相同的采樣數(shù)據(jù)量下，塊缺失采樣方式比隨機缺失采樣方式的成像性能低，在目標(biāo)3維重構(gòu)結(jié)果中“局外點”會增加。同時，利用本文的算法分別在全孔徑和3/4稀疏孔徑采樣數(shù)據(jù)量下進行轉(zhuǎn)動參數(shù)估計，并將全孔徑下的估計結(jié)果作為參考值。表1為轉(zhuǎn)動參數(shù)估計結(jié)果?？梢娫谙∈璨蓸訑?shù)據(jù)量下的估計結(jié)果與全孔徑下非常接近，驗證了本文算法的有效性。

圖 3 機動目標(biāo)聯(lián)合稀疏約束InISAR成像算法流程圖

圖4 民航飛機實測數(shù)據(jù)

圖5 3/4數(shù)據(jù)量的稀疏孔徑數(shù)據(jù)（通道A）

表1 轉(zhuǎn)動參數(shù)估計結(jié)果

圖6 本文算法的InISAR 2維成像結(jié)果（3/4數(shù)據(jù)量）

圖7 本文算法的InISAR 3維重構(gòu)結(jié)果（3/4數(shù)據(jù)量）

5 結(jié)束語

本文主要針對機動目標(biāo)研究了稀疏孔徑InISAR 3維成像算法。對于勻加速轉(zhuǎn)動的目標(biāo)，本文可以用 chirp-傅里葉字典來表示其回波多普勒調(diào)制。同時，考慮到 InISAR通道獨立成像處理不利于保持通道之間的相干性，本文提出了利用InISAR多通道圖像的結(jié)構(gòu)稀疏性進行聯(lián)合稀疏表征處理的多通道聯(lián)合處理2維成像算法。在稀疏孔徑觀測下，由于2維成像結(jié)果存在較多的“局外點”，本文采用對3維重構(gòu)和轉(zhuǎn)速估計進行聯(lián)合優(yōu)化求解的方式有效地實現(xiàn)了干涉相位濾波，得到了精確的目標(biāo)3維幾何重構(gòu)結(jié)果。最后，基于實測數(shù)據(jù)實驗驗證了本文算法的有效性。

[1] 保錚, 邢孟道, 王彤. 雷達成像技術(shù)[M]. 北京: 電子工業(yè)出版社, 2005: 229-276.

[2] 李麗亞, 劉宏偉, 曹向海, 等. 基于 InISAR 像的目標(biāo)識別方法[J]. 電子與信息學(xué)報, 2008, 30（9）: 2089-2093.Li Li-ya, Liu Hong-wei, Cao Xiang-hai, et al.. Radar automatic target recognition based on InISAR images[J]. Journal of Electronics & Information Technology, 2008, 30（9）: 2089-2093.

[3] 李飛, 糾博, 劉宏偉. 基于隨機霍夫變換的干涉 ISAR橫向定標(biāo)算法[J]. 電子與信息學(xué)報, 2013, 35（1）: 49-55.Li Fei, Jiu Bo, and Liu Hong-wei. A novel method of cross-range scaling interferometric ISAR based on randomized hough transform[J]. Journal of Electronics &Information Technology, 2013, 35（1）: 49-55.

[4] Wang Gen-yuan, Xia Xiang-gen, and Chen Victor C.Three-dimensional ISAR imaging of maneuvering target using three receivers[J]. IEEE Transactions on Image Processing, 2001, 10（3）: 436-447.

[5] Zhang Qun, Yeo T S, Du Gan, et al.. Estimation of three-dimensional motion parameters in interferometric ISAR imaging[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42（2）: 292-300.

[6] Xu Gang, Xing Meng-dao, Zhang Lei, et al.. Sparse-apertures ISAR imaging and scaling for maneuvering targets[J]. IEEE Journal of Selected Topics in Applied Earth Observation and Remote Sensing, 2014, 7（7）: 2942-2956.

[7] Liu Ya-bo, Li N, and Wang R. Achieving high-quality three-dimensional InISAR imageries of maneuvering target via super-resolution ISAR imaging by exploiting sparseness[J]. IEEE Geoscience and Remote Sensing Letters, 2014, 11（4）: 828-832.

[8] 徐剛, 楊磊, 張磊, 等. 一種加權(quán)最小熵的 ISAR 自聚焦算法[J]. 電子與信息學(xué)報, 2011, 33（8）: 1809-1815.Xu Gang, Yang Lei, Zhang Lei, et al.. Weighted minimum entropy autofocus algorithm for ISAR imaging[J]. Journal of Electronics & Information Technology, 2011, 33（8）: 1809-1815.

[9] Xu Gang, Xing Meng-dao, and Zhang Lei. Bayesian inverse synthetic aperture radar imaging[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8（6）: 1150-1154.

[10] 吳敏, 邢孟道, 張磊. 基于壓縮感知的二維聯(lián)合超分辨 ISAR成像算法[J]. 電子與信息學(xué)報, 2014, 36（1）: 187-193.Wu Min, Xing Meng-dao, and Zhang Lei. Two dimensional Joint super-resolution ISAR imaging algorithm based on compressive sensing[J]. Journal of Electronics & Information Technology, 2014, 36（1）: 187-193.

[11] Xia Xiang-gen. Discrete chirp-Fourier transform and its application to chirp rate estimation[J]. IEEE Transactions on Signal Processing, 2000, 48（11）: 3122-3133.

[12] Guo Xin, Sun H B, and Wang S L. Comments on “discrete chirp-Fourier transform and its application to chirp rate estimation”[J]. IEEE Transactions on Signal Processing, 2002,50（12）: 3115-3116.

[13] Ramakrishnan N, Ertin E, and Moses R L. Enhancement of coupled multi-channel images using sparsity constraints[J]. IEEE Transactions on Image Processing, 2010, 19（8）: 2115-2126.

[14] Eldar Y and Rauhut H. Average case analysis of multichannel sparse recovery using convex relaxation[J]. IEEE Transactions on Information Theory, 2010, 56（1）: 505-519.

[15] Baraniuk R G, Cevher V, Duarte M F, et al.. Model-based compressive sensing[J]. IEEE Transactions on Information Theory, 2010, 56（4）: 1982-2001.

[16] 陳倩倩, 徐剛, 李亞超, 等. 短孔徑 ISAR 方位定標(biāo)[J]. 電子與信息學(xué)報, 2013, 35（8）: 1854-1861.Chen Qian-qian, Xu Gang, Li Ya-chao, et al.. Cross-range scaling for ISAR with short aperture data[J]. Journal of Electronics &Information Technology, 2013, 35（8）: 1854-1861.