基于加權(quán)最大范數(shù)的SAR自聚焦方法

黃大榮 張 磊 邢孟道 周 峰 保 錚

?

基于加權(quán)最大范數(shù)的SAR自聚焦方法

黃大榮*張 磊 邢孟道 周 峰 保 錚

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

基于最大似然估計的特征向量分解自聚焦算法利用最大特征值對應(yīng)的特征向量實現(xiàn)對相位誤差的估計。該方法雖然具備精確和穩(wěn)健的性能，但需要對協(xié)方差矩陣進行特征分解，導致實際數(shù)據(jù)在處理中運算量巨大，對內(nèi)存要求也很高，難以在實時合成孔徑雷達(SAR)成像處理中應(yīng)用。該文提出一種基于加權(quán)最大范數(shù)的自聚焦方法，通過求解二范數(shù)最大化的優(yōu)化函數(shù)對目標特征向量進行估計，避免了特征值的分解過程，有效提升了運算效率；利用信噪比加權(quán)的思想，對不同距離單元賦予不同的權(quán)值，增強了優(yōu)質(zhì)特顯點樣本對相位誤差的估計貢獻，有效改善了自聚焦精度。通過實測SAR和ISAR數(shù)據(jù)處理驗證了算法的有效性。

逆合成孔徑雷達；自聚焦；加權(quán)信噪比；最大范

1 引言

經(jīng)典的自聚焦方法主要包括：特顯點法[4](Prominent Point Processing, PPP)和相位梯度自聚焦方法[5](Phase Gradient Autofocus, PGA)，這兩種方法沒有對誤差相位的模型做任何的假設(shè)，因此可以估計出大部分的相位誤差。但是，由于它們在聚束模式的假設(shè)條件下，僅利用相鄰兩次回波作為相位誤差的估計，因此在其它模式下聚焦精度受到一定的損失。加權(quán)相位估計自聚焦方法[6](Weig- hted Phase Estimation, WPE)對于相位誤差的階數(shù)無要求，具有很強的魯棒性。但該算法先提取各像素單元的相位信息，并逐個進行相位展開，巨大的運算量無法應(yīng)用于實際。最大似然估計的特征向量分解法[7]利用最大特征值對應(yīng)的特征向量估計相位誤差，估計的方差可以達到克拉美羅界，因此，可以比其它方法獲得更理想的聚焦效果。但由于該方法需要對協(xié)方差矩陣進行特征分解，龐大的運算量無法滿足實時處理要求，且對內(nèi)存要求很高。

為了解決特征向量分解自聚焦算法運算量大的問題，本文提出一種基于加權(quán)最大范數(shù)的自聚焦方法。通過求解2-范數(shù)最大化的優(yōu)化函數(shù)，對目標特征向量直接估計，避免了特征值分解過程，從而大大提升了運算效率；利用信噪比加權(quán)的思想，對不同距離單元賦予不同權(quán)值，增強了優(yōu)質(zhì)特顯點樣本對相位誤差的估計貢獻，有效改善了自聚焦精度。最后，通過實測數(shù)據(jù)處理，證明本文提出的方法具有比PGA方法更好的聚焦效果，在與特征向量分解法聚焦精度相當?shù)那闆r下，運算效率遠高于特征向量分解法。

本文第2節(jié)給出了含有相位誤差的信號模型；第3節(jié)論述了相位誤差的估計及補償方法，包括多普勒相位的補償，隨機初相的補償以及相位誤差的補償3部分；第4節(jié)對算法的運算復雜度進行分析；第5節(jié)給出了整個算法的流程圖；第6節(jié)給出了實測數(shù)據(jù)處理及結(jié)果分析，驗證了算法的有效性；第7節(jié)對全文進行了總結(jié)。

2 信號模型

3 相位誤差的估計與補償

3.1多普勒相位的補償

3.2 隨機初相的補償

圖1 特顯點樣本信號導向矢量

每個特顯點樣本信號信噪比不同，其對相位誤差估計精度的貢獻也不同。因此，在對隨機初相補償前，可通過對信噪比高的特顯點樣本信號加高權(quán)值，對質(zhì)量較差的特顯點樣本信號加低權(quán)值，最終優(yōu)化隨機初相的補償精度，權(quán)值的選擇[3,11]滿足：

3.3 相位誤差的補償

當完成多普勒相位補償和隨機初相補償后，特顯點樣本信號相位表示如下：

4 運算復雜度分析

圖 2 運算復雜度比值與采樣點M的關(guān)系

5 算法流程

綜上所述，整個算法的流程如圖3所示，具體描述如表1。

表1算法流程

算法流程步驟1 特顯點單元選擇。由式(4)篩選出高信噪比的距離單元，并按歸一化幅度方差從大到小排列，構(gòu)成特顯點樣本信號。步驟2 多普勒相位補償。首先對特顯點樣本信號進行傅里葉變換到圖像域，乘以式(6)將特顯點循環(huán)移位到圖像的方位中心位置，消除多普勒相位對相位誤差估計的影響。步驟3 加窗濾波。通過預(yù)先設(shè)定逐步減小的窗長，過濾大部分雜波能量，提高樣本的信雜比。步驟4 計算權(quán)值。由式(15)計算出各特顯點回波的權(quán)值，對經(jīng)過多普勒相位補償后的特顯點樣本信號加權(quán)，提高特顯點的信噪比。步驟5 隨機初相補償。利用式(14)補償隨機初相對相位誤差估計的影響。步驟6 相位誤差補償。利用式(18)估計并補償相位誤差，補償時通過線性擬合濾除估計中的線性分量以避免引入圖像偏移。步驟7 重復迭代步驟2至步驟6，直到估計達到收斂精度。步驟8 相位誤差的補償。對所有距離壓縮后的信號進行步驟6的操作。

6 實際數(shù)據(jù)處理與分析

為了驗證本文算法的有效性，本節(jié)對兩組實測數(shù)據(jù)進行對比分析。運算平臺為主頻3.2 GHz的個人計算機，運算代碼采用MATLAB.7.10.0平臺編寫。

6.1 ISAR實驗及性能分析

圖5表示每次迭代估計出的相位誤差。其中，實線表示由特顯點單元直接估計出的相位誤差，將估計出的相位誤差補償后，進行1次迭代，得到圖5中虛線所示的相位誤差，再經(jīng)過1次相位誤差補償，得到圖5中點線所示的相位誤差?？梢姡?jīng)過3次迭代后，估計出的相位誤差趨近于零，與3.3節(jié)分析的吻合，說明本文提出的相位誤差估計和補償方法是收斂的。

為了保證實驗的公平性，各算法均采用相同的迭代次數(shù)和相同的窗長，采用本文方法最終的成像結(jié)果如圖4(c)所示。作為比較，同時給出采用PGA方法的處理結(jié)果，如圖4(b)所示。由圖4(a)和圖4(b)可知，相比于不采用自聚焦處理的圖4(a)，采用PGA自聚焦方法后，圖像聚焦效果有很大的提高。但是，PGA方法未獲得理想的聚焦效果，圖4(b)中還是有許多散焦單元，采用本文提出的方法后大大改善了圖4(b)中散焦的狀況。

采用圖像的熵定量分析方法的聚焦效果。圖像熵[9]的定義為

圖6所示為特顯點單元的1維方位剖面圖。這里選取第33距離單元的包絡(luò)進行比較，如圖6(a)所示，圖6(b)是圖6(a)中標注部分的局部放大圖。其中，虛線表示PGA算法，實線表示本文方法，可以明顯看出本文方法相比于PGA方法，聚焦處理后信號最窄幅度最高，因而本文方法聚焦效果要優(yōu)于PGA算法。

圖4 成像結(jié)果分析

圖5 相位誤差估計

圖6 1維方位剖面圖

6.2 SAR數(shù)據(jù)實驗及性能分析

將圖7(a)中所標注的部分放大后，得到圖8(a)所示。畫出其中特顯點的1維方位向剖面圖，如圖8(b)所示。其中，實線表示原始數(shù)據(jù)，點劃線表示散焦的數(shù)據(jù)，點線表示特征向量分解法處理的結(jié)果，虛線表示本文方法處理的結(jié)果，可見，特征向量分解法和本文方法處理結(jié)果基本接近原始數(shù)據(jù)，都可以達到理想的聚焦效果。

為了說明本文方法對場景的適應(yīng)性，圖9給出了不同場景(曠野和城區(qū))下各算法的最終處理結(jié)果。郊區(qū)曠野場景成像中孤立散射中心很少，自聚焦處理中樣本選擇是較為困難的，而城區(qū)場景中具有分布較廣的建筑的強散射，對自聚焦處理而言，強建筑散射并不是理想的孤立散射，其具有較寬的聚焦響應(yīng)特性，也將對自聚焦形成較強的挑戰(zhàn)。圖9表明，對于不同的場景，本文方法具有較好的穩(wěn)健性。

7 結(jié)束語

為了解決實際應(yīng)用中特征向量分解自聚焦算法運算量大的缺點，本文提出一種加權(quán)最大范數(shù)的自聚焦方法。通過對目標函數(shù)的最優(yōu)化求解，避免了對矩陣進行特征值分解的操作，大大降低了運算量；利用信噪比加權(quán)的思想，對不同距離單元按信噪比賦予不同權(quán)值，提高了強特顯點樣本信號對相位誤差估計的貢獻，最終優(yōu)化相位誤差估計精度。實測數(shù)據(jù)處理結(jié)果表明，本文提出的方法具有比PGA方法更好的聚焦效果，在與特征向量分解自聚焦算法精度相當?shù)那闆r下，運算復雜度遠低于特征向量分解自聚焦算法。

圖7 Sandia實驗室無人機SAR圖像

圖8 聚焦精度分析

圖9 不同場景處理結(jié)果

[1] Jiang R, Zhu D, Shen M,.. Synthetic aperture radar autofocus based on projection approximation subspace tracking[J].&, 2012, 6(6): 465-471.

[2] Schulz T. Optimal sharpness function for SAR autofocus[J]., 2007, 14(1): 27-30.

[3] Zhang L, Qiao Z, Xing M,.. A robust motion compensation approach for UAV SAR imagery[J]., 2012, 50(8): 3202-3218.

[4] Carrara W, Goodman R, and Majewski R. Spotlight Synthetic Aperture Radar: Signal Processing Algorithms[M]. Norwood: Artech House, 1995: 268-282.

[5] Wahl D. Phase gradient autofocus-a robust tool for high resolution SAR phase correction[J]., 1994, 30(3): 827-834.

[6] Ye W, Yeo T, and Bao Z. Weighted least-square estimation of phase errors for SAR/ISAR autofocus[J]., 1999, 37(5): 2488-2494.

[7] Jakowatz C and Wahl D. Eigenvector method for maximum-likelihood estimation of phase errors in synthetic-aperture-radar imagery[J]., 1993, 10(12): 2539-2546.

[8] 王琦. 空間目標ISAR成像的研究[D]. [博士論文], 西安電子科技大學,2007.

Wang Q. Study of ISAR imaging for space targets[D]. [Ph.D. dissertation], Xidian University, 2007.

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

[10] Zhu D, Wang L, Yu Y,.. Robust ISAR range alignment via minimizing the entropy of the average range profile[J]., 2009, 6(2): 204-208.

[11] Li Y, Liu C, Wang Y,.. A robust motion error estimation method based on raw data[J]., 2012, 50(7): 2780-2790.

[12] Yang B. Projection approximation subspace tracking[J]., 1995, 43(1): 95-107.

黃大榮： 男，1986年生，博士生，研究方向為SAR成像.

張 磊： 男，1984年生，講師，研究方向為ISAR超分辨、稀疏微波信號處理.

邢孟道： 男，1975年生，教授，博士生導師，主要研究方向為雷達成像和目標識別等.

周 峰： 男，1980年生，副教授，主要研究方向為SAR成像、成像雷達電子對抗等.

保 錚： 男，1927年生，教授，中國科學院院士，主要研究方向為SAR成像、目標識別、天線超視距、空時二維自適應(yīng)等.

A Weighted Maximize Norm Method for SAR Autofocus

Huang Da-rong Zhang Lei Xing Meng-dao Zhou Feng Bao Zheng

(,,710071,)

The eigenvector method for maximum-likelihood estimation of phase error can obtain ideal performance of phase error estimation by using the eigenvector corresponding to its largest eigenvalue. Although the method is accurate and robust, it requires eigen-decomposition of the sample covariance matrix, which is computationally expensive and limits its real-time applications. In this paper, a Weighted Maximum Norm Method (WMNM) for phase error estimation is proposed. The eigenvector of the maximum eigenvalue can be obtained directly by solving the problem of maximizing L-2 norm, which avoids the eigen-decomposition of the sample covariance matrix and reduces the computational cost greatly. By adding different weights to each range bin, the contribution of the range cells with high SNR can be enhanced. Experimental results of the measured data by SAR and Inverse SAR (ISAR) verify the validity of the proposed algorithm.

Inverse Synthetic Aperture Radar (ISAR); Autofocus; Weight Signal-to-Noise Ratio (WSNR); Maximum norm

TN958

A

1009-5896(2014)01-0202-07

10.3724/SP.J.1146.2012.01699

2012-12-24收到，2013-09-26改回

國家自然科學基金(61301280)和中央高?；究蒲袠I(yè)務(wù)費(K5051302001, K5051302038)資助課題

黃大榮 rsphdr@163.com