基于半正定规划的压缩感知线阵三维SAR自聚焦成像算法

韦顺军 田博坤 张晓玲 师君

引用本文:
Citation:

基于半正定规划的压缩感知线阵三维SAR自聚焦成像算法

    作者简介: 韦顺军(1983–),男,广西柳州人,博士,2013年获电子科技大学工学博士学位,目前为电子科技大学信息与通信工程学院副教授,主要从事合成孔径雷达成像、阵列雷达3维成像等技术研究,已发表论文30余篇。E-mail: weishunjun@uestc.edu.cn;田博坤(1993–),男,河北沧州人,电子科技大学博士生,主要从事合成孔径雷达成像研究。E-mail: 3544348143@qq.com;张晓玲(1964–),女,四川成都人,博士,2002年获电子科技大学工学博士学位,目前为电子科技大学信息与通信工程学院教授,博士生导师,主要从事SAR成像技术、雷达探测技术研究,已发表论文50余篇。E-mail: xlzhang@uestc.edu.cn;师 君(1979–),男,河南南阳人,博士,2009年获电子科技大学工学博士学位,目前为电子科技大学信息与通信工程学院副教授,博士生导师,主要从事SAR成像技术、雷达信号处理研究,已发表论文50余篇。E-mail: shijun@uestc.edu.cn.
    通讯作者: 韦顺军, weishunjun@uestc.edu.cn
  • 基金项目:

    国家自然科学基金(61501098),博士后面上基金(2015M570778),高分青年基金项目(GFZX04061502),中央高校科研基本业务费(ZYGX2016KYQD107)

  • 中图分类号: TN957.52

Compressed Sensing Linear Array SAR Autofocusing Imaging via Semi-definite Programming

    Corresponding author: Wei Shunjun, weishunjun@uestc.edu.cn ;
  • Fund Project: The National Natural Science Foundation of China (61501098), The China Postdoctoral Science Foundation (2015M570778), The High Resolution Earth Observation Youth Foundation (GFZX04061502), The Fundamental Research Funds for the Central Universities (ZYGX2016KYQD107)

    CLC number: TN957.52

  • 摘要: 线阵合成孔径雷达(Linear Array Synthetic Aperture Radar, LASAR)3维成像技术是一种具有重要潜在应用价值的新体制成像雷达,压缩感知稀疏重构是近几年实现LASAR高分辨3维成像的热点研究之一。但相对于传统2维SAR,受线阵稀疏分布及阵列-平台2维联动,压缩感知LASAR成像面临回波数据欠采样、多维度高阶相位误差等问题,传统SAR自聚焦算法难以适用于压缩感知LASAR 3维稀疏自聚焦成像。为克服欠采样条件下多维度高阶相位误差对LASAR成像的影响,该文提出了一种基于半正定规划的压缩感知LASAR自聚焦成像算法。首先,结合压缩感知成像理论、图像最大锐度及最小均方误差准则,构造欠采样条件下稀疏目标的相位误差估计模型;其次,利用松弛半正定规划方法估计相位误差;最后,利用迭代逼近方法提高相位误差估计精度,实现压缩感知LASAR高精度稀疏自聚焦成像。另外,通过主散射目标区域提取,仅采用主散射区域进行相位误差估计,进一步提高自聚焦算法运算效率。仿真数据和实测数据验证了该文算法的有效性。
  • 图 1  LASAR正下视3维成像的几何模型

    Figure 1.  The geographic model of down-looking LASAR imaging

    图 3  沿航迹及切航迹随机分布相位误差下仿真点目标成像结果(上:全部回波数据;中:50%回波数据;下:25%回波数据)

    Figure 3.  The results of the point targets in the case of the along-track random phase errors and cross-track random phase errors(Top: all samples; Middle: 50% samples; Bottom: 25% samples)

    图 4  SDPSA算法MSE变化曲线

    Figure 4.  The MSE curve of SDPSA

    图 2  沿航迹二次项及切航迹随机分布相位误差下仿真点目标成像结果(上:全部回波数据;中:50%回波数据;下:25%回波数据)

    Figure 2.  The results of the point targets in the case of the along-track quadratic phase errors and cross-track random phase errors(Top: all samples; Middle: 50% samples; Bottom: 25% samples)

    图 5  地基等效LASAR成像实验

    Figure 5.  The ground-based LASAR experiment

    图 6  地基LASAR合成阵列平面及回波数据

    Figure 6.  The virtual array antenna and the echo of the ground-based LSAR

    图 8  路灯目标实测数据稀疏自聚焦成像结果(左:BP-PGA算法;中:IRLS-PGA算法;右:SDPSA算法)

    Figure 8.  The sparse autofocusing results of the experimental light data(Left: BP-PGA; Middle: IRLS-PGA; Right: SDPSA)

    表 1  SDPSA算法

    Table 1.  SDPSA algorithm

    输入:代价函数 ${J_{\rm{1}}}\left( {{{{f}}},{{{φ}}} } \right)$ ${J_{\rm{2}}}\left( {{{{f}}},{{{φ}}} } \right)$
    输出:估计值 $\left( {{{{\stackrel \frown{f}}}} ,{{\stackrel \frown{{φ}}} } } \right)$
    初始化:迭代次数 $i = 0$ ,相位误差 ${{{{φ}}} ^{\left( {\rm{0}} \right)}} = 0$ ,门限 ${{δ}} $ ,迭代总次数 $K$
    ${\left\| {{{{{{\stackrel \frown{f}}}} }^{\left( {i + 1} \right)}} - {{{{{\stackrel \frown{f}}}} }^{\left( i \right)}}} \right\|_2}\Bigr/{\left\| {{{{{{\stackrel \frown{f}}}} }^{\left( i \right)}}} \right\|_2} > {{δ}} $ $i \le K$ ,循环开始
     步骤1 构造稀疏目标重构的代价函数:
    ${J_{\rm{1}}}\left( {{{{f}}},{{{{{\stackrel \frown{φ}}} } }^{\left( i \right)}}} \right) = {{λ}} {\left\| {{{f}}} \right\|_1}{\rm{ + }} \left\| {{{{y}}} - {{{R}}}\left( {{{{{{\stackrel \frown{φ}}} } }^{\left( i \right)}}} \right){{{Af}}}} \right\|_2^{}$
     步骤2 利用IRLS算法进行稀疏目标重构:
    ${{{{\stackrel \frown{f}}}} ^{\left( {i{\rm{ + 1}}} \right)}} = \arg \mathop {\min }\limits_{f} {J_{\rm{1}}}\left( {{{{f}}},{{{{{\stackrel \frown{φ}}} } }^{\left( i \right)}}} \right)$
     步骤3 构造相位误差估计的代价函数:
    ${J_{\rm{2}}}\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i + 1} \right)}},{{{φ}}} } \right) \approx {{{γ} }^{\rm H}}{{{{Q}}}^{\left( {i + 1} \right)}}{{γ} }$
     步骤4 利用SDP方法求解最优问题:
    ${{{{\stackrel \frown{X}}}} _{\rm{opt}}} = \arg \mathop {\min }\limits_{X} {\rm tr}\left({{{{Q}}}^{\left( {i + 1} \right)}}{{{X}}}\right),\ {\rm s.t.} \ \ {{{X}}} \succeq 0,{{{{X}}}_{ii}} = {g_n},n = 1,2, ·\!·\!· ,N$
     步骤5 利用 ${{{{\stackrel \frown{X}}}} _{\rm {opt}}}$ 估计相位误差 ${{{{\stackrel \frown{φ}}} } ^{\left( {i + 1} \right)}}$
    $i \leftarrow i + 1$
    循环结束
    返回结果: ${{{\stackrel \frown{f}}}} \leftarrow {{{{\stackrel \frown{f}}}} ^{\left( {i + 1} \right)}}$ , $\stackrel \frown{{{{φ}}} } \leftarrow {\stackrel \frown{{{{φ}}} } ^{\left( {i + 1} \right)}}$
    下载: 导出CSV
  • [1] Du L, Wang Y P, Hong W, et al. A three-dimensional range migration algorithm for downward-looking 3D-SAR with single-transmitting and multiple-receiving linear array antennas[J]. EURASIP Journal on Advances in Signal Processing, 2010, 2010: 957916. DOI: 10.1155/2010/957916
    [2] Liao K F, Zhang X L, and Shi J. Plane-wave synthesis and RCS extraction via 3-D linear array SAR[J]. IEEE Antennas and Wireless Propagation Letters, 2015, 14: 994–997. DOI: 10.1109/LAWP.2015.2389264
    [3] Han K Y, Wang Y P, Tan W X, et al. Efficient pseudopolar format algorithm for down-looking linear-array SAR 3-D imaging[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(3): 572–576. DOI: 10.1109/LGRS.2014.2351792
    [4] Zhang S Q, Zhu Y T, and Kuang G Y. Imaging of downward-looking linear array three-dimensional SAR based on FFT-MUSIC[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(4): 885–889. DOI: 10.1109/LGRS.2014.2365611
    [5] Wei S J, Zhang X L, and Shi J. Linear array SAR imaging via compressed sensing[J]. Progress In Electromagnetics Research, 2011, 117: 299–319. DOI: 10.2528/PIER11033105
    [6] Zhang S Q, Zhu Y T, Dong G G, et al. Truncated SVD-based compressive sensing for downward-looking three-dimensional SAR imaging with uniform/nonuniform linear array[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(9): 1853–1857. DOI: 10.1109/LGRS.2015.2431254
    [7] Zhang S Q, Dong G G, Kuang G Y, et al. Superresolution downward-looking linear array three-dimensional SAR imaging based on two-dimensional compressive sensing[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(6): 2184–2196. DOI: 10.1109/JSTARS.2016.2549548
    [8] Peng X M, Tan W X, Hong W, et al. Airborne DLSLA 3-D SAR image reconstruction by combination of polar formatting and L1 regularization[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(1): 213–226. DOI: 10.1109/TGRS.2015.2453202
    [9] Tian J H, Sun J P, Han X, et al.. Motion compensation for compressive sensing SAR imaging with autofocus[C]. Proceedings of the 6th IEEE Conference on Industrial Electronics and Applications (ICIEA), Beijing, China, 2011: 1564–1567. DOI: 10.1109/ICIEA.2011.5975839.
    [10] Cetin M, Stojanovic I, Onhon O, et al. Sparsity-driven synthetic aperture radar imaging: Reconstruction, autofocusing, moving targets, and compressed sensing[J]. IEEE Signal Processing Magazine, 2014, 31(4): 27–40. DOI: 10.1109/MSP.2014.2312834
    [11] Onhon N Ö and Cetin M. A sparsity-driven approach for joint SAR imaging and phase error correction[J]. IEEE Transactions on Image Processing, 2012, 21(4): 2075–2088. DOI: 10.1109/TIP.2011.2179056
    [12] Zhe Z, Yao Z, Jiang C L, et al.. Autofocus of sparse microwave imaging radar based on phase recovery[C]. Proceedings of 2013 IEEE International Conference on Signal Processing, Communication and Computing (ICSPCC), Kunming, China, 2013: 1–5. DOI: 10.1109/ICSPCC.2013.6663989.
    [13] Chen Y C, Li G, Zhang Q, et al. Motion compensation for airborne SAR via parametric sparse representation[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(1): 551–562. DOI: 10.1109/TGRS.2016.2611522
    [14] Camlica S, Gurbuz A C, Arikan O, et al. Autofocused spotlight SAR image reconstruction of off-grid sparse scenes[J]. IEEE Transactions on Aerospace and Electronic Systems, 2017, 53(4): 1880–1892. DOI: 10.1109/TAES.2017.2675138
    [15] Uḡur S and Arıkan O. SAR image reconstruction and autofocus by compressed sensing[J]. Digital Signal Processing, 2012, 22(6): 923–932. DOI: 10.1016/j.dsp.2012.07.011
    [16] Kelly S, Yaghoobi M, and Davies M. Sparsity-based autofocus for undersampled synthetic aperture radar[J]. IEEE Transactions on Aerospace and Electronic Systems, 2014, 50(2): 972–986. DOI: 10.1109/TAES.2014.120502
    [17] Uḡur S, Arıkan O, and Gürbüz A C. Off-grid sparse SAR image reconstruction by EMMP algorithm[C]. Proceedings of 2013 IEEE Radar Conference (RADAR), Ottawa, ON, Canada, 2013: 1–4. DOI: 10.1109/RADAR.2013.6586034.
    [18] Wei S J and Zhang X L. Sparse autofocus recovery for under-sampled linear array SAR 3-D imaging[J]. Progress In Electromagnetics Research, 2013, 140: 43–62. DOI: 10.2528/PIER13020614
    [19] Wei S J, Zhang X L, and Shi J. Sparse autofocus via Bayesian learning iterative maximum and applied for LASAR 3-D imaging[C]. Proceedings of 2014 IEEE Radar Conference, Cincinnati, OH, USA, 2014: 666–669. DOI: 10.1109/RADAR.2014.6875674.
    [20] Donoho D L. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289–1306. DOI: 10.1109/TIT.2006.871582
    [21] Figueiredo M A T, Nowak R D, and Wright S J. Gradient projection for sparse reconstruction: Application to compressed sensing and other inverse problems[J]. IEEE Journal of Selected Topics in Signal Processing, 2007, 1(4): 586–597. DOI: 10.1109/JSTSP.2007.910281
    [22] Ji S H, Xue Y, and Carin L. Bayesian compressive sensing[J]. IEEE Transactions on Signal Processing, 2008, 56(6): 2346–2356. DOI: 10.1109/TSP.2007.914345
    [23] Grant M and Boyd S. CVX: Matlab software for disciplined convex programming, version 1.21[R]. CVX Research, Inc., 2010. Available from: URL: http://cvxr.com/cvx.
    [24] Toh K C, Todd M J, and Tütüncü R H. SDPT3—A Matlab software package for semidefinite programming, version 1.3[J]. Optimization Methods and Software, 1999, 11(1/4): 545–581. DOI: 10.1080/10556789908805762
    [25] Liu K H, Wiesel A, and Munson D C. Synthetic aperture radar autofocus via semidefinite relaxation[J]. IEEE Transactions on Image Processing, 2013, 22(6): 2317–2326. DOI: 10.1109/TIP.2013.2249084
  • [1] 王爱春向茂生 . 基于块压缩感知的SAR层析成像方法. 雷达学报, 2016, 5(1): 57-64. doi: 10.12000/JR16006
    [2] 顾福飞张群杨秋霍文俊王敏 . 基于NCS算子的大斜视SAR压缩感知成像方法. 雷达学报, 2016, 5(1): 16-24. doi: 10.12000/JR15035
    [3] 肖鹏吴有明于泽李春升 . 一种基于压缩感知恢复算法的SAR图像方位模糊抑制方法. 雷达学报, 2016, 5(1): 35-41. doi: 10.12000/JR16004
    [4] 刘向阳杨君刚孟进张晓牛德智 . 低信噪比下基于Hough变换的前视阵列SAR稀疏三维成像. 雷达学报, 2017, 6(3): 316-323. doi: 10.12000/JR17011
    [5] 田鹤李道京 . 稀疏重航过阵列SAR运动误差补偿和三维成像方法. 雷达学报, 2018, 7(6): 717-729. doi: 10.12000/JR18101
    [6] 闫敏韦顺军田博坤张晓玲师君 . 基于稀疏贝叶斯正则化的阵列SAR高分辨三维成像算法. 雷达学报, 2018, 7(6): 705-716. doi: 10.12000/JR18067
    [7] 张增辉郁文贤 . 稀疏微波SAR图像特征分析与目标检测研究. 雷达学报, 2016, 5(1): 42-56. doi: 10.12000/JR15097
    [8] 李杭梁兴东张福博吴一戎 . 基于高斯混合聚类的阵列干涉SAR三维成像. 雷达学报, 2017, 6(6): 630-639. doi: 10.12000/JR17020
    [9] 杨军张群罗迎邓冬虎 . 基于压缩感知的认知雷达多目标跟踪方法. 雷达学报, 2016, 5(1): 90-98. doi: 10.12000/JR14107
    [10] 吴一戎洪文张冰尘蒋成龙张柘赵曜 . 稀疏微波成像研究进展(科普类). 雷达学报, 2014, 3(4): 383-395. doi: 10.3724/SP.J.1300.2014.14105
    [11] 赵曜张冰尘洪文吴一戎 . 基于RIPless 理论的稀疏微波成像波形分析方法. 雷达学报, 2013, 2(3): 265-270. doi: 10.3724/SP.J.1300.2013.13032
    [12] 仲利华胡东辉丁赤飚张问一 . 一种稀疏孔径下大尺寸目标的ISAR 成像方法. 雷达学报, 2012, 1(3): 292-300. doi: 10.3724/SP.J.1300.2012.20033
    [13] 陈文峰李少东杨军马晓岩 . 基于线性Bregman迭代类的多量测向量ISAR成像算法研究. 雷达学报, 2016, 5(4): 389-401. doi: 10.12000/JR16057
    [14] 任笑真杨汝良 . 一种基于幅度和相位迭代重建的四维合成孔径雷达成像方法. 雷达学报, 2016, 5(1): 65-71. doi: 10.12000/JR15135
    [15] 刘振魏玺章黎湘 . 一种新的随机PRI脉冲多普勒雷达无模糊MTD算法. 雷达学报, 2012, 1(1): 28-35. doi: 10.3724/SP.J.1300.2013.10063
    [16] 赵娟白霞 . 一种适用于TDOMP算法的测量矩阵优化方法. 雷达学报, 2016, 5(1): 8-15. doi: 10.12000/JR15131
    [17] 廖明生魏恋欢汪紫芸TimoBalz张路 . 压缩感知在城区高分辨率SAR层析成像中的应用. 雷达学报, 2015, 4(2): 123-129. doi: 10.12000/JR15031
    [18] 李烈辰李道京黄平平 . 基于变换域稀疏压缩感知的艇载稀疏阵列天线雷达实孔径成像. 雷达学报, 2016, 5(1): 109-117. doi: 10.12000/JR14159
    [19] 赵逸超朱宇涛粟毅杨猛 . 用于线阵三维SAR成像的二维快速ESPRIT算法. 雷达学报, 2015, 4(5): 591-599. doi: 10.12000/JR15065
    [20] 钟金荣文贡坚 , . 基于块稀疏贝叶斯学习的雷达目标压缩感知(英文). 雷达学报, 2016, 5(1): 99-108. doi: 10.12000/JR15056
  • 加载中
图(8)表(1)
计量
  • 文章访问数:  425
  • HTML浏览量:  143
  • PDF下载量:  178
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-11-09
  • 录用日期:  2018-03-28
  • 网络出版日期:  2018-05-23
  • 刊出日期:  2018-12-28

基于半正定规划的压缩感知线阵三维SAR自聚焦成像算法

    通讯作者: 韦顺军, weishunjun@uestc.edu.cn
    作者简介: 韦顺军(1983–),男,广西柳州人,博士,2013年获电子科技大学工学博士学位,目前为电子科技大学信息与通信工程学院副教授,主要从事合成孔径雷达成像、阵列雷达3维成像等技术研究,已发表论文30余篇。E-mail: weishunjun@uestc.edu.cn;田博坤(1993–),男,河北沧州人,电子科技大学博士生,主要从事合成孔径雷达成像研究。E-mail: 3544348143@qq.com;张晓玲(1964–),女,四川成都人,博士,2002年获电子科技大学工学博士学位,目前为电子科技大学信息与通信工程学院教授,博士生导师,主要从事SAR成像技术、雷达探测技术研究,已发表论文50余篇。E-mail: xlzhang@uestc.edu.cn;师 君(1979–),男,河南南阳人,博士,2009年获电子科技大学工学博士学位,目前为电子科技大学信息与通信工程学院副教授,博士生导师,主要从事SAR成像技术、雷达信号处理研究,已发表论文50余篇。E-mail: shijun@uestc.edu.cn
  • 电子科技大学电子工程学院   成都   611731
基金项目:  国家自然科学基金(61501098),博士后面上基金(2015M570778),高分青年基金项目(GFZX04061502),中央高校科研基本业务费(ZYGX2016KYQD107)

摘要: 线阵合成孔径雷达(Linear Array Synthetic Aperture Radar, LASAR)3维成像技术是一种具有重要潜在应用价值的新体制成像雷达,压缩感知稀疏重构是近几年实现LASAR高分辨3维成像的热点研究之一。但相对于传统2维SAR,受线阵稀疏分布及阵列-平台2维联动,压缩感知LASAR成像面临回波数据欠采样、多维度高阶相位误差等问题,传统SAR自聚焦算法难以适用于压缩感知LASAR 3维稀疏自聚焦成像。为克服欠采样条件下多维度高阶相位误差对LASAR成像的影响,该文提出了一种基于半正定规划的压缩感知LASAR自聚焦成像算法。首先,结合压缩感知成像理论、图像最大锐度及最小均方误差准则,构造欠采样条件下稀疏目标的相位误差估计模型;其次,利用松弛半正定规划方法估计相位误差;最后,利用迭代逼近方法提高相位误差估计精度,实现压缩感知LASAR高精度稀疏自聚焦成像。另外,通过主散射目标区域提取,仅采用主散射区域进行相位误差估计,进一步提高自聚焦算法运算效率。仿真数据和实测数据验证了该文算法的有效性。

English Abstract

    • 合成孔径雷达(Synthetic Aperture Radar, SAR)作为一种具有全天时、全天候、高分辨力、大区域观测及信息量丰富的雷达成像技术,已成为当今遥感对地观测的重要手段。因具备3维成像能力,线阵SAR (Linear Array Synthetic Aperture Radar, LASAR)是近几年来被广泛关注的一种新体制SAR 3维成像技术[14]。本质上,LASAR主要利用线阵天线运动合成一个大的虚拟2维阵列实现观测目标的2维分辨,再利用距离脉冲压缩技术实现观测目标的3维成像。相对于其它常规SAR 3维成像体制(如圆周SAR、层析SAR等), LASAR具备成像模式多样、应用灵活等优势,克服了常规SAR 3维成像体制无法下视、前视成像的缺陷,在城市地形测绘、飞行器导航及盲降、目标精定位与识别、战场情报获取等民用和军用领域具有重要的研究价值和应用前景。

      目前,LASAR 3维成像处理主要采用基于匹配滤波理论的成像算法[1,2],如距离多普勒算法和后向投影算法。受载荷平台(如飞机、无人机、卫星等)空间及线阵天线长度限制,传统匹配滤波成像算法在LASAR线阵方向的分辨能力通常比其它维向的分辨力低,一定程度约束了LASAR高分辨成像应用。然而,由于地面、空中等观测目标在3维空间通常具有强稀疏性,故可利用目标稀疏先验信息提升LASAR 3维成像精度。近几年,基于压缩感知(Compressed Sensing, CS)理论的稀疏成像已经成为LASAR高分辨成像的研究热点之一,并提出了多种基于CS稀疏重构的LASAR 3维高分辨稀疏成像方法[58]。但是,CS稀疏重构对成像模型的精确性要求非常高,如在X波段雷达时位置精度通常需要达到毫米级,现有压缩感知LASAR 3维稀疏成像算法通常以成像模型精确已知为前提,并没有考虑实际平台运动误差的影响。当LASAR存在运动误差、阵列抖动时,则会导致成像模型存在偏差,CS稀疏重构出现散焦、畸变、虚假目标等成像质量恶化现象,甚至不能成像。另外,与传统SAR系统单个或少量天线工作不同,LASAR需要成百上千天线阵元同时工作,故仅利用雷达系统中单个位置的导航测量系统(如IMU (Inertial Measurement Unit)系统、GPS (Global Positioning System)或北斗系统)数据难以实现LASAR多个天线阵元运动误差的精确补偿。为了提高压缩感知LASAR 3维成像质量,在利用导航测量数据进行运动误差补偿后,还需结合回波数据和自聚焦成像技术实现多天线阵元的残余运动误差补偿。

      至今,针对传统SAR成像中的运动误差补偿,相关学者已经提出了多种基于不同准则的自聚焦成像算法,如最大似然估计(Maximum Likelihood Estimation, MLE)自聚焦、最小熵估计(Minimum Entropy Estimation, MEE)自聚焦、最大锐度(Maximum Sharpness Estimation, MSE)自聚焦等方法。相位梯度自聚焦(Phase Gradient Autofocus, PGA)算法作为一种在SAR实际成像中广泛应用的典型MLE自聚焦成像算法,已被成功用于压缩感知SAR稀疏成像中的残余相位误差补偿[9]。但是,当回波数据欠采样时,传统自聚焦算法不一定适用于压缩感知SAR稀疏自聚焦成像[10]。针对压缩感知SAR 稀疏自聚焦问题,近几年相关学者也已经提出一些解决方法。文献[11]提出了一种基于稀疏驱动的联合SAR稀疏成像及相位误差校正方法,采用迭代估计方法求解非二次正则化问题实现SAR稀疏自聚焦成像。文献[12]提出一种针对稀疏微波成像的自聚焦方法,采用贪婪算法实现相位误差的高精度估计。文献[13]提出了一种基于子孔径参数稀疏表示的机载SAR自聚焦成像算法,实测数据结果表明该方法相对于传统PGA算法具有更好的自聚焦性能。文献[14]针对运动误差造成的网格偏离现象,提出了一种基于稀疏重构的聚束式SAR自聚焦成像算法,采用正交匹配追踪及最速下降法实现稀疏自聚焦成像。另外,还有一些文献提出基于其它准则或估计模型的SAR 稀疏自聚焦算法[1517]。然而,与传统SAR运动误差形式不同,由于受线阵稀疏分布及平台颤动、线阵抖动等联合运动误差影响,LASAR成像模型中相位误差具有多维度、高阶项分布等更复杂的形式,现有 SAR稀疏自聚焦算法不一定适用于LASAR自聚焦成像。对此,本文作者也提出了基于稀疏贝叶斯学习的LASAR自聚焦成像算法,在欠采样数据条件下可有效估计相位误差[18,19],但是该方法需已知目标及相位的先验分布特性并且要求较高的回波信噪比条件。

      为了补偿欠采样条件下LASAR多维度高阶相位误差,本文提出了一种基于半正定规划(Semi-Definite Programming, SDP)的压缩感知LASAR自聚焦成像算法,即半正定规划自聚焦(Semi-Definite Programming Sparse Autofocus, SDPSA)算法,简称为SDPSA算法。该方法利用LASAR压缩感知成像原理、图像最大锐度及最小均方误差准则,构造了欠采样条件下稀疏目标相位误差估计模型,并通过松弛半正定规划及迭代逼近方法实现压缩感知LASAR的高精度稀疏自聚焦成像。为了进一步提高SDPSA算法自聚焦运算效率,利用主散射目标区域提取,仅采用主散射区域进行相位误差估计。仿真数据和实测数据验证了本文SDPSA算法的有效性,相对于PGA自聚焦算法,SDPSA算法在欠采样数据时具有更好的成像质量。本文的主要结构如下:第2节简单介绍了压缩感知LASAR成像模型和稀疏成像原理,阐述了其阵列多相位中心的相位误差影响;第3节建立LASAR相位误差估计的半正定规划模型,提出了基于半正定规划的压缩感知LASAR稀疏自聚焦成像方法,并介绍了该方法的基本步骤;第4节利用LASAR仿真及实测数据,分析SDPSA稀疏自聚焦成像算法性能,验证了算法的有效性;第5节给出了本文算法的主要结论。

    • LASAR 3维成像的典型几何模型如图1所示,其中X, YZ轴分别表示切航迹、沿航迹和高度向,线阵天线平行于X轴放置。假设LASAR系统工作于正下视成像模式, $ {{{{P}}}_A} =\{{{\cal Q }_n} = \left( {{x_n},{y_n},{z_n}} \right); $ $ n \in \left[ {1,2, ·\!·\!· ,{N_A}} \right] \} $ 表示合成孔径时间内LASAR合成的2维等效阵列天线相位中心(Antenna Phase Center, APC)位置集, ${N_A}$ 为2维APC总个数, ${x_n}$ , ${y_n}$ ${z_n}$ 分别为第 $n$ 个APC位置 ${\cal Q_n}$ X, YZ轴位置。在远场条件下,LASAR观测场景可近似为点散射模型,并对观测场景进行3维离散化。令 ${{{{P}}}_S} \!=\! \left\{ {{{\cal P}_m} \!=\! \left( {{x_m},{y_m},{z_m}} \right);m \in {{{Ω}}} } \right\}$ ${{{{σ}}} _S} \!=\! \left\{\! {{{{σ}} _m};m \!\!\in\!\! {{{Ω}}} } \right\}$ 分别表示LASAR离散观测场景的位置集和散射系数集,其中 ${{{Ω}}} = \left[ {1,2, ·\!·\!· ,{M_S}} \right]$ 为观测场景的散射单元序列, ${M_S}$ 为散射单元总数, ${x_m}$ , ${y_m}$ ${z_m}$ 分别为第 $m$ 个散射单元 ${{\cal P}_m}$ X,YZ轴位置, ${{{σ}} _m}$ ${{\cal P}_m}$ 的散射系数。

      图  1  LASAR正下视3维成像的几何模型

      Figure 1.  The geographic model of down-looking LASAR imaging

      假设LASAR系统发射线性调频信号,其原始回波经过距离压缩后,散射单元 ${{\cal P}_m}$ 回波信号可表示为:

      ${s_r}\left( {r,n,m} \right) = {{{σ}} _m}{{{χ}} _R}\left( {r - {R_{n,m}}} \right)\exp \left( { - {\rm j}2k{R_{n,m}}} \right)$

      其中, $r\ $ 表示距离域, $k$ 为波数, ${R_{n,m}} = {\left\| {{{\cal Q}_n} - {{\cal P}_m}} \right\|_2}$ 为散射单元 ${{\cal P}_m}$ 到第 $n$ 个APC位置 ${{\cal Q}_n}$ 的距离, ${{{χ}} _R}\left( \cdot \right)$ 为LASAR距离脉压模糊函数。式(1)第1项 ${{{χ}} _R}\left( {r - {R_{n,m}}} \right)$ 提供了目标距离向分辨力,第2项 $\exp \left( { - {\rm j}2k{R_{n,m}}} \right)$ 包含了目标沿航向及切航向信息。对于离散全观测场景,LASAR回波信号可表示为各个散射单元回波之和:

      ${s_r}\left( {r,n} \right) = \sum\limits_{m \in {{{Ω}}} } {{s_r}\left( {r,n,m} \right)} $

      此时,回波信号 ${s_r}\left( {r,n} \right)$ 可表示为散射系数向量 ${{{f}}} \in {\mathbb{C}^{{M_S} \times 1}}$ 和时延相位向量 ${{{ψ}}} \left( {{r_i},n} \right) \in {\mathbb{C}^{{M_S} \times 1}}$ 的乘积:

      $s\left( {{r_i},n} \right) \!=\! {{{ψ}}} {\left( {{r_i},n} \right)^{\rm T}}{{{f}}},\ i \!=\! 1, ·\!·\!· ,{N_R}, n \!=\! 1, ·\!·\!· ,{N_A}$

      其中, ${{{f}}} = {\rm{Vec}}\left[ {{{{σ}} _m}} \right],m = 1,·\!·\!·, {M_S}$ 为观测场景 ${{{Ω}}} $ 的散射系数向量, ${{{ψ}}} \left( {{r_i},n} \right) = {\rm{Vec}}\left[ {{{{χ}} _R}\left( {r - {R_{n,m}}} \right) } \right.$ $ \left. \cdot\exp \left( { - {\rm j}2k{R_{n,m}}} \right) \right]$ 为LASAR时延相位向量, ${N_R}$ 为距离采样点数, ${\rm{Vec}}\left[ \cdot \right]$ 表示向量化符号。因此,LASAR回波信号线性测量模型可以表示为:

      ${{{y}}} = {{{Af}}} + {{{n}}}$

      其中, ${{{y}}} = {\rm{Vec}}\left[ {{s_r}\left( {{r_i},n} \right)} \right] \in {\mathbb{C}^{{N_R}{N_A} \times 1}}$ 表示LASAR 脉压后的回波信号向量, ${{{A}}} = {\left[ {{{{ψ}}} \left( {{r_i},n} \right)} \right]^{\rm T}} \in {\mathbb{C}^{{N_R}{N_A} \times {M_S}}}$ 为回波 ${{{y}}}$ 的线性测量矩阵, ${{{n}}}$ 为回波中的噪声向量。

    • CS理论指出当信号可稀疏表示时,可利用远低于Nyquist采样率进行采样,将原始信号重构问题转化为线性约束最优化求解问题,实现稀疏信号的精确重构[20]。若LASAR观测场景散射系数 ${{{f}}}$ 具有稀疏特征,根据CS稀疏重构理论,散射系数 ${{{f}}}$ 可通过求解如下 ${\ell _1}$ 范数最优化问题进行重构:

      ${\stackrel \frown {{f}}} = \arg \mathop {\min }\limits_{f} {\left\| {{f}} \right\|_1},\ {\rm{s}}{\rm{.t.}} {\left\| {{{y}} - {{Af}}} \right\|_2} \le {{ε}} $

      其中, ${{ε}} \ge {\left\| {{{n}}} \right\|_2}$ 为回波信号 ${{{y}}}$ 中噪声门限。若测量矩阵 ${{{A}}}$ 满足低相关条件,式(5)可等效为 ${\ell _1}{{\tiny{-}}}{\ell _2}$ 联合范数最小化问题求解:

      ${\stackrel \frown {{f} } = \arg \mathop {\min }\limits_{f}{ \Bigr( {{{λ}} {{\left\| {{{f}}} \right\|}_1}{\rm{ + }} \left\| {{{{y}}} - {{{Af}}}} \right\|_2^{}} \Bigr) $

      其中, ${{λ}} \ge 0$ 为稀疏约束正则化参数,用于均衡重构误差及约束信号稀疏性。针对式(5)和式(6)中稀疏信号 ${{{f}}}$ 求解,至今学者们已提出多种经典稀疏重构算法,如正交匹配追踪(Orthogonal Matching Pursuit, OMP)算法、最小绝对值收缩和选择(Least Absolute Shrinkage Selection Operator, LASSO)算法、迭代加权 ${\ell _1}$ 范数最小二乘(Iteration Reweighted ${\ell _1}$ norm Least Square, IRLS)算法、贝叶斯压缩感知(Bayesian Compressed Sensing, BCS)算法等[21,22]

    • 由于受平台颤动、线阵抖动等联合运动误差影响,LASAR 2维APC存在位置误差,会导致 LASAR成像模型存在多维度、高阶项分布的复杂相位误差。令 ${{{φ}}} = {\left[ {{{{φ}} _1},{{{φ}} _2}, ·\!·\!· ,{{{φ}} _{{N_R}{N_A}}}} \right]^{\rm T}\ }$ 表示LASAR回波数据中2维APC位置不确定引起的相位误差向量,此时LASAR回波信号的线性测量模型表示为:

      ${{{y}}} = {{{R}}}\left( {{{φ}}} \right){{{Af}}} + {{{n}}}$

      其中, ${{{R}}}\left( {{{φ}}} \right) = \operatorname{diag} \left\{ {{{\rm e}^{{\rm j}{{{φ}} _1}}},{{\rm e}^{{\rm j}{{{φ}} _2}}}, ·\!·\!· ,{{\rm e}^{{\rm j}{{{φ}} _{{N_R}{N_A}}}}}} \right\}$ 表示相位误差矩阵。针对式(7),LASAR稀疏自聚焦成像目的是根据已知的回波 ${{{y}}}$ 及测量矩阵 ${{{A}}}$ ,估计未知散射系数向量 ${{{f}}}$ 和相位误差向量 ${{{φ}}} $ ,其可表示为以下 ${\ell _1} {\tiny{-}} {\ell _2}$ 组合范数的联合最优化问题:

      $ \begin{align} \left( {{{ {{\stackrel \frown{f}}}} },{{{\stackrel \frown {φ}}} } } \right) =& \arg \mathop {\min }\limits_{{{{f}}},{{{φ}}} } J\left( {{{f}},{{φ}} } \right) \\ = & \arg \mathop {\min }\limits_{{{f}},{{{φ}}} } \Bigr( {{{λ}} {{\left\| {{f}} \right\|}_1}{\rm{ + }} \left\| {{ {y}} - {{{R}}}\left( {{{φ}}} \right){{{Af}}}} \right\|_2^{}} \Bigr) \end{align} $

      若线阵天线阵元稀疏分布或回波信号欠采样,LASAR自聚焦稀疏本质上是一个病态求逆问题,直接求解 $J\left( {{{{f}}},{{{φ}}} } \right)$ 最优化问题十分困难。目前,稀疏自聚焦算法通常采用基于某个准则的迭代算法以估计 $\left( {{{{\stackrel \frown {f}}}} ,{{{\stackrel \frown {φ}}} } } \right)$ 。如,文献[11,13]采用稀疏重构联合最小MSE 准则实现SAR稀疏自聚焦,文献[19]采用最大后验概率准则实现稀疏自聚焦,而不同稀疏自聚焦准则具有不同的优劣及特点。

    • 为了实现欠采样数据条件下LASAR稀疏自聚焦,联合压缩感知LASAR成像原理、图像最大锐度及最小均方误差重构准则,本文提出一种基于半正定规划的LASAR稀疏自聚焦成像算法。该方法采用迭代估计方法,给定估计初值及参数后,主要由稀疏目标重构、相位误差估计两个独立过程循环实现,其具体理论推导及实现过程如下。

      (1) 稀疏目标重构。在算法第 $i + 1$ 次迭代中,根据式(8)及第 $i$ 次迭代的相位误差 ${{{\stackrel \frown{{φ}}} } ^{\left( i \right)}}$ ,构造目标重构代价函数 ${J_{\rm{1}}}\left( {{{{f}}},{{{{\stackrel \frown{{φ}}} } }^{\left( i \right)}}} \right)$ ,则稀疏目标重构过程可表示为:

      $ \begin{align} {{\stackrel \frown {f}} ^{\left( {i{\rm{ + 1}}} \right)}} =& \arg \mathop {\min }\limits_{{f}} J\left( { {f},{{{\stackrel \frown {φ}} }^{\left( i \right)}}} \right) \\ =& \arg \mathop {\min }\limits_{{f}} \left( {\lambda {{\left\| {{f}} \right\|}_1}\! +\! \left\| {{{y}} \!-\! {{R}}\left( {{{{\stackrel \frown {φ}} }^{\left( i \right)}}} \right){{Af}}} \right\|_2^{}} \right) \end{align} $

      该过程为标准CS稀疏重构,可采用现有稀疏重构方法求解 ${{{{\stackrel \frown {f}}}} ^{\left( {i{\rm{ + 1}}} \right)}}$ 。兼顾精度及稳健性,本文利用IRLS算法进行稀疏目标重构,此时第 $i + 1$ 次迭代得到的LASAR稀疏目标散射系数估计向量为:

      $ \begin{align} {{{{\stackrel \frown {f}}}} ^{\left( {i + 1} \right)}} =& {\left[ {{{{{A}}}^{\rm H}}{{{A}}} + {{λ}} {{Λ}}\left( {{{{{{\stackrel \frown {f}}}} }^{\left( {i + 1} \right)}}} \right)} \right]^{ - 1}}\\ & \cdot {{{{A}}}^{\rm H}}{{{{R}}}^{\rm H}}\left( {{\stackrel \frown { {{φ}} }^{\left( i \right)}}} \right){{{y}}} \end{align} $

      其中, ${{Λ}}\left( {{{f}}} \right)$ 为对角矩阵,且 ${{Λ}}\left( {f} \right) = {\rm diag}\left\{ {{{\rm{1}} \mathord{\left /\!\! {\vphantom {{\rm{1}} {\sqrt {{{\left| {{f_1}} \right|}^2} + \eta } }}} \right.} {\sqrt {{{\left| {{f_1}} \right|}^2} \!+\! \eta } }},} \right. $ $\left. {·\!·\!· ,{{\rm{1}} \mathord{\left/ {\vphantom {{\rm{1}} {\sqrt {{{\left| {{f_{{M_S}}}} \right|}^2} + \eta } }}} \right.} {\sqrt {{{\left| {{f_{{M_S}}}} \right|}^2} + \eta } }}} \right\}$ , ${{η}} $ 为小常数,一般选择为 ${{η}} = {10^{ - 6}}$

      (2) 相位误差估计。利用式(9)得到的 ${{{{\stackrel \frown {f}}}} ^{\left( {i{\rm{ + 1}}} \right)}}$ ,联合图像最大锐度及最小MSE准则构造相位误差估计的代价函数 ${J_{\rm{2}}}\left( {{{\stackrel \frown{{{{f}}}} }^{\left( {i{\rm{ + 1}}} \right)}},{{{φ}}} } \right)$

      对于最大锐度准则,采用传统匹配滤波重构结果进行判定,因 ${{{\stackrel \frown {f}}}} _{\rm {MF}}^{}\left( {{{φ}}} \right) = {{{A}}}_{}^{\rm H}{{{R}}}_{}^{\rm H}\left( {{{φ}}} \right){\cal {y}}$ ,则其图像锐度可近似表示为:

      $ \begin{align} {{{χ}} _{\rm {MF}}}\left( {{φ}} \right) \approx &{\stackrel \frown {f}} _{\rm {MF}}^{\rm H}\left( {{{φ}}} \right){{{\stackrel \frown{f}}}} _{\rm {MF}}^{}\left( {{{φ}}} \right) \\ = & { {y}}_{}^{\rm H}{{{R}}}\left( {{{φ}}} \right){{{AA}}}_{}^{\rm H}{{{R}}}_{}^{\rm H}\left({{{φ}}} \right) { {y}} \end{align} $

      对于最小MSE准则,根据式(9),回波数据的重构MSE表示为:

      $ E\left( {{{{{{\stackrel \frown {f}}}} }^{\left( {i{\rm{ + 1}}} \right)}}} \right) = \left\| {{\cal {y}} - {{{R}}}\left( {{{φ}}} \right){{{A}}}{{{{{\stackrel \frown {f}}}} }^{\left( {i{\rm{ + 1}}} \right)}}} \right\|_2^{} $

      相位误差估计的代价函数 ${J_{\rm{2}}}\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i{\rm{ + 1}}} \right)}},{{{φ}}} } \right)$ 构造为:

      ${J_{\rm{2}}}\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i{\rm{ + 1}}} \right)}},{{{φ}}} } \right) = E\left( {{{{{\stackrel \frown{{f}}}} }^{\left( {i{\rm{ + 1}}} \right)}}} \right) - {{β}} {{{χ}} _{\rm {MF}}}\left( {{{φ}}} \right)$

      其中, ${{β}} > 0$ 为加权因子,其控制图像锐度准则的权重,本文选择 ${{β}} = {\rm{1}}$ 。可知, ${{β}} \to 0$ 时,式(13)相位误差估计等效为仅利用最小均方误差(Mean Square Error, MSE)准则,其估计过程与文献[11,13]相似; ${{β}} \to \infty $ 时,相位误差估计等效为用图像最大锐度准则。

      因此,第 $i + 1$ 次迭代中SDPSA算法求解代价函数 ${J_{\rm{2}}}\left( {{{{{{\stackrel \frown {f}}}} }^{\left( {i{\rm{ + 1}}} \right)}},{{{φ}}} } \right)$ 最小化问题估计LASAR阵列APC的相位误差 ${{{{\stackrel \frown{φ}}} } ^{\left( {i{\rm{ + 1}}} \right)}}$ ,其表示为:

      $ \begin{align} {{{{\stackrel \frown {φ}}} } ^{\left( {i{\rm{ + 1}}} \right)}} =& \arg \mathop {\min }\limits_{{{φ}}} {J_{\rm{2}}}\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i + 1} \right)}},{{{φ}}} } \right) \\=& \arg \mathop {\min }\limits_{{{φ}}} \left\{ {E\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i{\rm{ + 1}}} \right)}}} \right) - {{β}} {{{χ}} _{\rm {MF}}}\left( {{{φ}}} \right)} \right\} \end{align} $

      若令相位误差向量 ${{γ} } = {{\rm e}^{ - {\rm j}{{φ}} }}$ ,矩阵 ${{{Y}}} = \operatorname{diag} \left\{ {{y}} \right\}$ ,因 ${{{{R}}}^{\rm H}}\left( {{{φ}}} \right){{{R}}}\left( {{φ}} \right) = {{I}}$ ,推导可得:

      $ \begin{align} {\rm{}}& E\left( {{{{{{\stackrel \frown {f}}}} }^{\left( {i{\rm{ + 1}}} \right)}}} \right) \\ {\rm{}}& \!=\! \left\| {{ y} \!-\! {{R}}\left( {{{φ}}} \right)\!\!{{{A}}}{{\left[ {{{{{A}}}^{\rm H}}{{{A}}} \!+\! {{λ}} {{Λ}}\left( {{{{{{\stackrel \frown {f}}}} }^{\left( {i + 1} \right)}}} \right)} \right]}^{ - 1}}\!\!\!\!{{{{A}}}^{\rm H}}{{{{R}}}^{\rm H}}\!\!\left( {{φ}} \right)\!{y}} \right\|_2 \\ {\rm{}} &\!=\! \! \left\| {{{{R}}}\left( \!{{φ}} \right)\!\!\left\{ {{{{I}}} \!-\! {{{A}}}{{\left[ {{{{{A}}}^{\rm H}}{{{A}}} \!\!+\!\! {{λ}} {{Λ}}\!\!\left( {{{{{{\stackrel \frown {f}}}} }^{\left( {i + 1} \right)}}} \right)} \right]}^{ - 1}}\!\!\!\!{{{{A}}}^{\rm H}}}\!\! \right\}\!\!{{{{R}}}^{\rm H}}\!\!\left( {{φ}} \right)\!{y}} \right\|_2^{} \\ {\rm{}} & \!=\!\! \left\| {{{{R}}}\left( {{φ}} \right)} \right\|_2^2\!\left\|\! {\left\{\!\! {{{{I}}} \!\!-\!\! {{{A}}}\!\!{{\left[\! {{{{{A}}}^{\rm H}}\!\!{{{A}}} \!\!+\!\! {{λ}} {{Λ}}\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i + 1} \right)}}} \right)} \!\right]}^{ - 1}}\!\!\!\!{{{{A}}}^{\rm H}}} \right\}\!\!{{{{R}}}^{\rm H}}\!\!\left(\! {{φ}} \!\right)\!{y}} \right\|_2 \\ {\rm{}}& \!=\!\! \left\| {{{{R}}}\left( {{φ}} \right)} \right\|_2^{}\left\|\!\! {\left\{ {{{{I}}}\! -\! {{{A}}}\!\!{{\left[ {{{{{A}}}^{\rm H}}{{{A}}} \!+\! {{λ}} {{Λ}}\left( {{{\stackrel \frown{{{{f}}}} }^{\left( {i + 1} \right)}}} \right)} \right]}^{ - 1}}\!\!\!{{{{A}}}^{\rm H}}}\!\! \right\}\!\!{{{Y}}}{{γ} }} \right\|_2 \\ {\rm{}} & \!=\! \left\| {{{{R}}}\left( {{φ}} \right)} \right\|_2^{}{{{γ} }^{\rm H}}{{{Φ} }^{\left( {i + 1} \right)}}{{γ} } \end{align} $

      ${{{χ}} _{\rm {MF}}}\left( {{{φ}}} \right) \approx {{y}}_{}^{\rm H}{{{R}}}\left( {{φ}} \right){{AA}}_{}^{\rm H}{{{R}}}_{}^{\rm H}\left( {{φ}} \right){y} = {{{{γ}}} ^{\rm H}}{{{Θ}}} {{{γ}}} $

      ${J_{\rm{2}}}\left( {{{\stackrel \frown{{{{f}}}} }^{\left( {i + 1} \right)}},{{{φ}}} } \right) \approx {{{γ} }^{\rm H}}{{{Φ} }^{\left( {i + 1} \right)}}{{γ} } - {{β} }{{{γ} }^{\rm H}}{{Θ} {{γ}} } $

      其中,矩阵 ${{{Φ} }^{\left( {i + 1} \right)}} \!=\!\! {{{{C}}}^{\rm H}}{{{C}}}$ 且矩阵 $ {C} \!=\! \Biggl\{ {{I} \!-\! {A}} \left[ {{{A}^{\rm H}}{A} \! } \right.$ $\left. +{{{\left. {\lambda {{Λ} }\left( {{{{\stackrel \frown {f}} }^{\left( {i + 1} \right)}}} \right)} \right]}^{{ - 1}}}\!\!\!{{A}^{\rm{H}}}} \right\}\!\!{{Y}\ $ ,矩阵 ${{Θ} } \!=\!\! {{{{D}}}^\rm{H}}{{{D}}}$ ${D} \!=\! {{{{A}}}^{\rm H}}{{{Y}}}$

      由于 $\left\| {{{{R}}}\left( {{{φ}}} \right)} \right\|_2^{}{\rm{ = 1/}}{N_A}{N_R}$ 为常数,令矩阵 ${{{{Q}}}^{\left( {i + 1} \right)}} = {{{Φ} }^{\left( {i + 1} \right)}} - {{β}} {{Θ} }$ ,相位误差估计代价函数近似为 ${J_{\rm{2}}}\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i + 1} \right)}},{{{φ}}} } \right) \approx {{{γ} }^{\rm H}}{{Q}^{\left( {i + 1} \right)}}{{γ} }$ ,则式(14)最优化问题变为:

      $ \begin{align} {\rm{}}&{{{{\stackrel \frown{φ}}} } ^{\left( {i{\rm{ + 1}}} \right)}} = \arg \mathop {\min }\limits_{{{φ}}} {J_{\rm{2}}}\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i + 1} \right)}},{{{φ}}} } \right) \\ {\rm{}}&\quad\quad\ \, = \arg \mathop {\min }\limits_{{γ}} \left( {{{{γ} }^{\rm H}}{{ Q}^{\left( {i + 1} \right)}}{{γ} }} \right), \\ {\rm{}}&{\rm {s.t.}}\ \left| {{{{γ}} _n}} \right| = 1,n = 1,2, ·\!·\!·,N \end{align} $

      显然,式(18)最优化问题可采用恒定常数二次规划(Constant Modulus Quadratic Programming, CMQP)方法求解。相位误差向量 ${{γ} }$ 的最优化估计可等效为在 ${{γ} }$ 幅度值为1约束条件下的线性最小化求解。但是,由于式(18)是非确定性多项式困难问题,当向量 ${{γ} }$ 维数大时,CMQP等式的最优化求解非常困难。因此,一般采用近似估计算法求解CMQP最优解。

      近年来,半正定规划(SDP)方法成为求解CMQP等式最优化问题的热门方法。对于CMQP最优化问题求解,SDP过程的典型表达式为:

      $ \begin{align} {\rm{}}& {{{{\stackrel \frown {X}}}} } = \arg \mathop {\min }\limits_{{X}} {\rm {tr}}\left({{{{Q}}}^{\left( {i + 1} \right)}}{{{X}}}\right),\\ {\rm{}}& {\rm s.t.}\ \ {{{X}}} \succeq 0,{\rm tr}({ {B}_n}{{{X}}}) = {g_n},n = 1,2, ·\!·\!· ,N \end{align} $

      其中, ${{{X}}} = {{γ} }{{{γ} }^{\rm H}}$ 表示LASAR相位误差矩阵, ${\rm {tr}}\left( \cdot \right)$ 表示矩阵对角元素之和。因 $ {{{γ} }^{\rm H}}{{{{Q}}}^{\left( {i + 1} \right)}}{{γ} } = $ ${\rm {tr}}\left( {{{γ} }{{{{Q}}}^{\left( {i + 1} \right)}}{{{γ} }^{\rm H}}} \right) = {\rm {tr}}\left( {{{{{Q}}}^{\left( {i + 1} \right)}}{{γ} }{{{γ} }^{\rm H}}} \right)$ ,式(19)的CMQP最优化求解可近似转换为以下半正定矩阵的最优化估计问题。

      $ \begin{align} {\rm{}}& \stackrel \frown{ X} _{\rm {opt}}^{\left( {i + 1} \right)} = \arg \mathop {\min }\limits_{{X}} {\rm {tr}}\left({{{Q}}^{\left( {i + 1} \right)}} X\right),\\ {\rm{}}& {\rm {s.t.}}\ \ {{X}} \succeq 0,{{{{X}}}_{ii}} = {g_n},n = 1,2, ·\!·\!· ,N \end{align} $

      为了求解式(20)中的SDP问题,可利用现有的凸优化算法估计求解,如共轭梯度方法、内点法等。另外,目前许多学者也开发了各种SDP求解软件工具箱,如CVX工具箱[23]、SDPT3工具箱[24]等,可用于式(20) SDP问题的快速求解。获得矩阵 ${{{\stackrel \frown{X}}}} _{\rm {opt}}^{\left( {i + 1} \right)}$ 后,最优向量 ${{{{\stackrel \frown{γ}}} } }_{\rm {opt}}^{\left( {i + 1} \right)}$ 可通过矩阵 ${{{\stackrel \frown{X}}}} _{\rm {opt}}^{\left( {i + 1} \right)}$ 分解得到,具体实现方法可借鉴文献[25],本文不再详述。若 ${{{{\stackrel \frown {γ}}} } }_{\rm {opt}}^{\left( {i + 1} \right)}$ 已估,则第 $i + 1$ 次迭代中压缩感知LASAR成像模型的相位误差估计 ${{{{\stackrel \frown{φ}}} } ^{\left( {i + 1} \right)}}$ 可表示为:

      ${{{{{\stackrel \frown{φ}}}} } ^{\left( {i + 1} \right)}} = - \angle {{{{\stackrel \frown{{γ}}}} } }_{\rm {opt}}^{\left( {i + 1} \right)}$

      获取相位误差估计 ${{{{{\stackrel \frown{φ}}}} } ^{\left( {i + 1} \right)}}$ 后,利用 ${{{{{\stackrel \frown{φ}}}} } ^{\left( {i + 1} \right)}}$ 对压缩感知LASAR成像模型进行补偿,采用迭代逼近方法循环处理,直至SDPSA算法结果满足精度要求。

    • 为了提高相位误差 ${\rm{{{φ}}}} $ 的估计精度并且减少估计过程的运算量,本文SDPSA算法中仅采用主散射区域构造式(14)中代价函数。主散射点区域的估计采用以下方法:(1) 在距离压缩后回波数据中根据图像强度选择主散射目标距离单元;(2)根据所选主散射目标距离单元图像估计目标主瓣宽度;(3)根据目标主瓣宽度通过阈值选择主散射目标区域。获取主散射目标区域后,令 ${{{{Ω}} }_0} = \left[ {1,2, ·\!·\!· ,{M_0}} \right]$ 表示这些主散射目标的序列集合, ${M_0}$ 是主散射目标单元的总个数。此时,主散射目标区域 ${{{{Ω}} }_0}$ 所对应的LASAR测量矩阵变为 ${{{A}}_{{{{{Ω}} }_0}}} \subset {{A}}$

      针对LASAR稀疏自聚焦成像,根据半正定规划自聚焦成像原理,本文SDPSA算法的主要步骤如表1所示,其主要步骤是迭代估计过程。

      输入:代价函数 ${J_{\rm{1}}}\left( {{{{f}}},{{{φ}}} } \right)$和 ${J_{\rm{2}}}\left( {{{{f}}},{{{φ}}} } \right)$
      输出:估计值 $\left( {{{{\stackrel \frown{f}}}} ,{{\stackrel \frown{{φ}}} } } \right)$
      初始化:迭代次数 $i = 0$,相位误差 ${{{{φ}}} ^{\left( {\rm{0}} \right)}} = 0$,门限 ${{δ}} $,迭代总次数 $K$
      若 ${\left\| {{{{{{\stackrel \frown{f}}}} }^{\left( {i + 1} \right)}} - {{{{{\stackrel \frown{f}}}} }^{\left( i \right)}}} \right\|_2}\Bigr/{\left\| {{{{{{\stackrel \frown{f}}}} }^{\left( i \right)}}} \right\|_2} > {{δ}} $和 $i \le K$,循环开始
       步骤1 构造稀疏目标重构的代价函数:
      ${J_{\rm{1}}}\left( {{{{f}}},{{{{{\stackrel \frown{φ}}} } }^{\left( i \right)}}} \right) = {{λ}} {\left\| {{{f}}} \right\|_1}{\rm{ + }} \left\| {{{{y}}} - {{{R}}}\left( {{{{{{\stackrel \frown{φ}}} } }^{\left( i \right)}}} \right){{{Af}}}} \right\|_2^{}$
       步骤2 利用IRLS算法进行稀疏目标重构:
      ${{{{\stackrel \frown{f}}}} ^{\left( {i{\rm{ + 1}}} \right)}} = \arg \mathop {\min }\limits_{f} {J_{\rm{1}}}\left( {{{{f}}},{{{{{\stackrel \frown{φ}}} } }^{\left( i \right)}}} \right)$
       步骤3 构造相位误差估计的代价函数:
      ${J_{\rm{2}}}\left( {{{{{{\stackrel \frown{f}}}} }^{\left( {i + 1} \right)}},{{{φ}}} } \right) \approx {{{γ} }^{\rm H}}{{{{Q}}}^{\left( {i + 1} \right)}}{{γ} }$
       步骤4 利用SDP方法求解最优问题:
      ${{{{\stackrel \frown{X}}}} _{\rm{opt}}} = \arg \mathop {\min }\limits_{X} {\rm tr}\left({{{{Q}}}^{\left( {i + 1} \right)}}{{{X}}}\right),\ {\rm s.t.} \ \ {{{X}}} \succeq 0,{{{{X}}}_{ii}} = {g_n},n = 1,2, ·\!·\!· ,N$
       步骤5 利用 ${{{{\stackrel \frown{X}}}} _{\rm {opt}}}$估计相位误差 ${{{{\stackrel \frown{φ}}} } ^{\left( {i + 1} \right)}}$;
      $i \leftarrow i + 1$
      循环结束
      返回结果: ${{{\stackrel \frown{f}}}} \leftarrow {{{{\stackrel \frown{f}}}} ^{\left( {i + 1} \right)}}$, $\stackrel \frown{{{{φ}}} } \leftarrow {\stackrel \frown{{{{φ}}} } ^{\left( {i + 1} \right)}}$

      表 1  SDPSA算法

      Table 1.  SDPSA algorithm

      一般情况下,SDPSA算法经过20次迭代左右可实现收敛。另外,SDPSA算法估计得到的相位误差可能会存在线性偏移项,该误差仅导致SDPSA自聚焦成像结果与原始图像存在固定偏移,此时可利用定标技术或其它线性相位误差估计方法进行补偿校正,本文不再详述。

    • 为了验证SDPSA稀疏自聚焦成像方法的性能,本节先利用LASAR仿真数据进行成像验证,并对比基于传统PGA相位误差估计补偿的后向投影(Back Projection, BP)算法及迭代最小二乘(IRLS)稀疏成像算法结果进行分析,分别简称为BP-PGA算法和IRLS-PGA算法。其中,IRLS-PGA算法是将BP-PGA算法得到的相位误差对回波数据进行补偿后,再利用IRLS算法进行稀疏成像。为了简便,该仿真实验中只对LASAR单个等距离单元切面进行成像,实际中对全部等距离单元切面进行成像后进行合并即可得LASAR观测场景3维成像。

      仿真中假设LASAR系统采用正下视工作模式,主要仿真参数如下:雷达中心频率 ${f_{\rm c}}$ 为30 GHz,发射信号带宽 ${B_{\rm r}}$ 为1 GHz,信号采样率 ${f_{\rm s}}$ 为1.2 GHz,脉冲重复频率PRF为1200,载荷平台飞行速度 ${V_{\rm S}}$ $\left[ {0,50,0} \right]{\rm{m/s}}$ ,飞行平台高度 $H$ 为1000 m,线阵天线长度 ${L_{\rm A}}$ 为5 m,线阵天线阵元为等间隔均匀分布,线阵天线阵元数 ${N_{\rm C}}$ 为128,系统沿航向采样点数 ${N_{\rm S}}$ 为128,距离向采样点数 ${N_{\rm R}}$ 为1024。值得注意的是,为了便于CS测量矩阵 ${{A}}$ 构造及计算,仿真中线阵天线阵元间隔约为4个波长,只适用于小观测场景成像仿真,但大场景仿真时应要求阵元间隔小于半个波长以避免栅瓣影响。仿真中原始场景存在4个单点目标,其散射系数相同。在成像处理中,等距离切面成像空间与原始场景相同,平面空间范围为64 m×64 m,被均匀离散化成64×64个分辨单元。为了分析SDPSA算法在不同类型相位误差情况下稀疏自聚焦性能,在仿真回波数据中分别加入缓变和陡变的相位误差,并且在方位向-切航向平面随机抽取50%和25%回波数据进行成像对比分析。

      首先,在回波数据中加入沿航向缓变、切航迹陡变的相位误差,其中沿航迹相位误差在 $\left[ {{\rm{0,8}}{{π}} } \right]$ 为二次项函数、切航迹相位误差在 $\left[ { - {{π}} ,{{π}} } \right]$ 区间上服从均匀随机分布,脉压后回波数据信噪比(Signal to Noise Ratio, SNR)为25 dB。图2给出了存在相位误差时全部回波数据、50%方位-切航随机抽取回波数据和25%方位-切航随机抽取回波数据情况下BP未自聚焦、BP-PGA、BP-SDPSA(即利用SDPSA算法估计得到的相位误差对BP算法进行补偿)、IRLS未自聚焦、IRLS-PGA和SDPSA算法的仿真点目标成像,其中图像显示经过归一化处理且幅度取对数,图像显示门限为–30 dB。从图2可知,在此相位误差情况下BP算法和IRLS算法未自聚焦时成像失效,点目标散焦严重;经过PGA算法相位误差补偿后,BP和IRLS沿航迹向成像质量得到一定提升,但切航迹向成像提升不明显,BP-PGA和IRLS-PGA结果非目标区域仍然存在大量的旁瓣或虚假目标,并且随着回波数据利用率减少自聚焦成像性能下降。相对于PGA算法相位误差补偿,SDPSA算法在全部回波数据、50%回波数据和25%回波数据条件下均得到点目标的良好聚焦成像结果,而且非目标区域几乎没有旁瓣和虚假目标。另外,经过SDPSA算法相位误差补偿后,在不同回波数据情况下BP-SDPSA成像结果较BP-PGA也抑制了沿航迹向和切航迹向的旁瓣水平,提高了自聚焦成像质量。

      其次,在回波数据沿航向和切航迹向均加入随机相位误差,其在 $\left[ { - {{π}} ,{{π}} } \right]$ 区间上服从均匀随机分布,脉压后回波数据SNR为25 dB。图3给出了存在相位误差时全部回波数据、50%回波数据和25%回波数据情况下BP未自聚焦、BP-PGA、BP-SDPSA、IRLS未自聚焦、IRLS-PGA和SDPSA算法的仿真点目标成像,其图像显示同图2。从图3可知,在此随机相位误差情况下PGA算法性能提升很小,尤其是非全部回波数据时BP-PGA和IRLS-PGA目标散射严重。然而, SDPSA算法在全部回波数据、50%回波数据和25%回波数据时都实现了点目标良好聚集成像,原始场景中4个单点目标被精确恢复,并且BP-SDPSA成像结果旁瓣水平较BP-PGA也大大降低。实验结果验证了SDPSA算法在压缩感知LASAR稀疏自聚焦成像的有效性,说明了该算法可对LASAR稀疏回波数据条件下缓变及陡变相位进行估计和补偿。

      图  3  沿航迹及切航迹随机分布相位误差下仿真点目标成像结果(上:全部回波数据;中:50%回波数据;下:25%回波数据)

      Figure 3.  The results of the point targets in the case of the along-track random phase errors and cross-track random phase errors(Top: all samples; Middle: 50% samples; Bottom: 25% samples)

      最后,为了分析SDPSA算法在不同SNR及迭代次数情况下自聚焦性能,将图2中SNR变为0至30 dB,图4(a)给出了全部、50%和25%回波数据量时SDPSA算法成像结果与原始仿真场景的MSE。可知,SDPSA算法的MSE随着SNR增大而变小,并且在SNR为0 dB至5 dB时MSE大于0.1,但SNR大于10 dB后MSE均小于0.1且变化很小。因此,当脉压后回波数据SNR大于10 dB时,SDPSA算法可良好实现稀疏自聚焦成像。对于LASAR原始回波数据,通常其脉冲压缩后回波数据SNR可优于10 dB,满足SDPSA算法稀疏自聚焦成像要求。图4(b)给出了回波数据SNR为20 dB,全部、50%和25%回波数据量时SDPSA算法MSE随迭代次数变化曲线,可知全部数据成像时5次迭代后基本收敛,而在50%和25%回波数据时10次迭代后基本收敛,此时算法MSE均小于0.1。

      图  4  SDPSA算法MSE变化曲线

      Figure 4.  The MSE curve of SDPSA

      图  2  沿航迹二次项及切航迹随机分布相位误差下仿真点目标成像结果(上:全部回波数据;中:50%回波数据;下:25%回波数据)

      Figure 2.  The results of the point targets in the case of the along-track quadratic phase errors and cross-track random phase errors(Top: all samples; Middle: 50% samples; Bottom: 25% samples)

    • 为了进一步验证本文SDPSA算法的有效性,利用本课题组地基等效LASAR实验系统获取的实测数据进行稀疏自聚焦成像分析。地基等效LASAR实验系统实物图如图5(a)所示,系统主要参数如下:雷达中心频率 ${f_{\rm c}}$ 为9.62 GHz,发射信号带宽 ${B_{\rm r}}$ 为80 MHz,信号采样率 ${f_{\rm s}}$ 为120 MHz,线阵长度 ${L_{\rm A}}$ 为1.25 m。实验场景为地面足球场并布置3个参考球目标,其光学图像如图5(b)图5(c)所示,包含了两种目标类型:简单的球目标和较复杂的路灯铁栅栏,该观测场景中心到实验平台的距离大约为100 m。

      图  5  地基等效LASAR成像实验

      Figure 5.  The ground-based LASAR experiment

      图6(a)给出了实验系统天线扫描的2维轨迹,可等效为一个虚拟的2维天线阵列,此虚拟阵列大小为1.25 m×1.25 m,阵列的阵元个数为10416,此时观测区域成像不会出现栅瓣。图6(b)为图6(a)对应的脉压后回波数据,其中球1回波对应于第18个距离单元数据、球2和球3对应于第29个距离单元数据,路灯对应于第76个距离单元数据。为了分析测量阵元样本数对压缩感知LASAR稀疏自聚焦成像的影响,本节图6(b)选择 $N = {\rm{10}}000$ 个阵元的脉压后回波数据,并以此数据作为参考的全部回波数据,然后分别从该10000样本阵元数中随机选择5000和2000个阵元的回波数据作为稀疏2维阵列回波。在回波数据加入2维阵列平面相位误差,该相位误差在 $\left[ { - {{π}} /2,{{π}} /2} \right]$ 区间上服从均匀随机分布,并对比传统PGA算法,分析压缩感知LASAR自聚焦成像性能。

      图  6  地基LASAR合成阵列平面及回波数据

      Figure 6.  The virtual array antenna and the echo of the ground-based LSAR

      图7给出了10000, 5000和2000阵元回波数据条件下采用BP-PGA, IRLS-PGA及本文SDPSA算法获得的实验球目标自聚焦成像结果,其中图像已归一化且显示门限为最大值–25 dB。从图7成像结果看出,在10000个阵元回波数据时3种算法对3个参考球目标均能良好自聚焦成像,且IRLS-PGA算法和本文SDPSA算法结果相似,说明PGA算法和本文算法在10000个阵元回波数据时自聚焦性能相当。但是,在5000和2000阵元回波数据成像时,BP-PGA算法出现严重旁瓣干扰,3个参考球目标在X-Z平面存在很高的旁瓣串扰;IRLS-PGA算法在球体目标周围出现虚假目标,且2000阵元样本时结果更加恶化,说明PGA自聚焦算法在数据欠采样时相位误差估计精度下降,并且随着样本数减少自聚焦性能变差。然而,SDPSA算法在5000和2000稀疏阵元样本时,成像结果与10000阵元基本一致,较IRLS-PGA算法结果大大抑制了虚假目标,说明SDPSA算法的相位误差估计性能优于PGA算法,在稀疏欠采样数据时仍可良好实现LASAR相位误差估计和校正。

      图8给出了10000, 5000和2000阵元回波数据条件下采用BP-PGA, IRLS-PGA及本文SDPSA算法获得的路灯铁栅栏目标自聚焦成像结果,其中图像已归一化且显示门限为最大值–20 dB。从图8成像结果看出,相对传统BP-PGA和IRLS-PGA算法,SDPSA算法在10000, 5000和2000阵元样本时对路灯和铁栅栏目标成像均有明显提升,同样验证了SDPSA算法稀疏自聚焦性能优于传统PGA算法,在稀疏阵元采样数据时可实现LASAR相位误差估计和校正。综上,地基LASAR实测数据实验结果验证了本文SDPSA算法的有效性。

      图  8  路灯目标实测数据稀疏自聚焦成像结果(左:BP-PGA算法;中:IRLS-PGA算法;右:SDPSA算法)

      Figure 8.  The sparse autofocusing results of the experimental light data(Left: BP-PGA; Middle: IRLS-PGA; Right: SDPSA)

    • 本文提出了一种基于半正定规划的压缩感知LASAR自聚焦3维成像算法。该SDPSA方法结合LASAR压缩感知成像原理及最小均方误差准则,利用松弛半正定规划方法及迭代最优方法估计相位误差,实现压缩感知LASAR的高精度自聚焦成像,并且为了进一步提高自聚焦算法运算效率,通过主散射目标区域进行相位误差估计。仿真数据和实测数据验证了本文SDPSA算法的有效性,实验结果说明相对于传统PGA自聚焦算法,SDPSA算法在稀疏采样数据条件下具备更好的稀疏自聚焦成像结果,可为LASAR 3维成像技术实际应用提供技术途径。在未来工作中,可开展算法优化研究,进一步提升SDPSA算法的运算效率。

参考文献 (25)

目录

    /

    返回文章
    返回