基于稀疏重构的全极化SAR联合多维重建

孙豆 路东伟 邢世其 杨潇 李永祯 王雪松

孙豆, 路东伟, 邢世其, 等. 基于稀疏重构的全极化SAR联合多维重建[J]. 雷达学报, 2020, 9(5): 865–877. doi:  10.12000/JR20092
引用本文: 孙豆, 路东伟, 邢世其, 等. 基于稀疏重构的全极化SAR联合多维重建[J]. 雷达学报, 2020, 9(5): 865–877. doi:  10.12000/JR20092
SUN Dou, LU Dongwei, XING Shiqi, et al. Full-polarization SAR joint multidimensional reconstruction based on sparse reconstruction[J]. Journal of Radars, 2020, 9(5): 865–877. doi:  10.12000/JR20092
Citation: SUN Dou, LU Dongwei, XING Shiqi, et al. Full-polarization SAR joint multidimensional reconstruction based on sparse reconstruction [J]. Journal of Radars, 2020, 9(5): 865–877. doi:  10.12000/JR20092

基于稀疏重构的全极化SAR联合多维重建

doi: 10.12000/JR20092
基金项目: 国家自然科学基金(61971429, 61901499)
详细信息
    作者简介:

    孙 豆(1992–),女,博士研究生,主要研究方向为极化雷达成像和雷达信号处理。E-mail: sundou14@nudt.edu.cn

    路东伟(1992–),男,博士研究生,主要研究方向为合成孔径雷达对抗和雷达目标识别。E-mail: bookwormldw@qq.com

    邢世其(1984–),男,副研究员,主要研究方向为极化雷达成像、雷达信号处理以及合成孔径雷达对抗。E-mail: xingshiqi_paper@163.com

    杨 潇(1983–),男,助教,主要研究方向为雷达信号处理和雷达目标识别。E-mail: 297414430@qq.com

    李永祯(1977–),男,研究员,博士生导师,主要研究方向为极化雷达与电子对抗。E-mail: e0061@sina.com

    王雪松(1972–),男,教授,博士生导师,主要研究方向为极化雷达、目标识别与电子对抗。E-mail: wxs1019@vip.sina.com

    通讯作者:

    邢世其 xingshiqi_paper@163.com

  • 责任主编:仇晓兰 Corresponding Editor: QIU Xiaolan
  • 中图分类号: TN95

Full-polarization SAR Joint Multidimensional Reconstruction Based on Sparse Reconstruction

Funds: The National Natural Science Foundation of China (61971429, 61901499)
More Information
  • 摘要: 各极化通道独立处理和三维分步成像会忽视数据之间的关联性,造成散射中心的失配以及极化散射矩阵获取的不准确。鉴于此,该文提出一种基于稀疏重构的全极化联合多维重建方法。该方法通过设置联合稀疏约束对所有极化通道及所有维度进行联合,将全极化多维重建建模为多通道联合稀疏重构问题。通过数据插值对模型简化后,结合三维快速傅里叶变换、共轭梯度法和牛顿迭代法给出一种高效的模型求解方法,可以同时得到极化散射矩阵和目标三维信息。该文方法保证了不同极化通道、不同维度的稀疏支撑集一致,且充分利用了数据之间的关联性带来的额外信息。基于仿真数据和电磁计算数据的实验结果表明,该方法的性能不受目标类型影响,具有一定的抗噪性,能有效地获取目标的多维重建结果,得到的三维成像结果分辨率高且极化散射矩阵估计精度高。
  • 图  1  各个极化通道的仿真目标成像结果

    Figure  1.  Each polarization channel’s imaging results of simulated targets

    图  2  仿真目标全极化联合多维重建结果

    Figure  2.  Full polarization joint multi-dimensional reconstruction results of simulated targets

    图  3  Slicy的CAD模型

    Figure  3.  CAD model of Slicy

    图  4  各个极化通道的Slicy成像结果

    Figure  4.  Each polarization channel’s imaging results of Slicy

    图  5  Slicy全极化联合多维重建结果

    Figure  5.  Full polarization joint multi-dimensional reconstruction results of Slicy

    图  6  卫星的CAD模型

    Figure  6.  CAD model of satellite

    图  7  各个极化通道的卫星成像结果

    Figure  7.  Each polarization channel’s imaging results of satellite

    图  8  卫星全极化联合多维重建结果

    Figure  8.  Full polarization joint multi-dimensional reconstruction results of satellite

    表  1  全极化联合多维重建方法的步骤

    Table  1.   Steps of full polarization joint multi-dimensional reconstruction method

     (1) 设定初值:${{\tilde{ \beta }}^n} = {{0}}$;
     (2) 对每个极化通道,
       (a) 使用3-D NUFFT对${G_l}(k_x^{},k_y^{},k_z^{})$进行插值,得到
         ${\hat G_l}(\hat k_x^{},\hat k_y^{},\hat k_z^{})$,
       (b) 对${\hat G_l}(\hat k_x^{},\hat k_y^{},\hat k_z^{})$进行3-D IFFT,并向量化结果得到
         ${{\hat{ A}}^{\rm{H}}}{{\hat{ b}}_l}$;
     (3) 根据所有极化通道的${{\hat{ A}}^{\rm{H}}}{{\hat{ b}}_l}$,得到${{\hat{ A}}^{\rm{H}}}{\hat{ b}}$;
     (4) 结合3-D FFT,3-D IFFT和共轭梯度法计算
       ${(2{{\hat{ A}}^{\rm{H}}}{\hat{ A}} + \mu p{{D}}({{\tilde{ \beta }}^n}))^{ - 1}}$;
     (5) 根据${{\hat{ A}}^{\rm{H}}}{\hat{ b}}$,计算${(2{{\hat{ A}}^{\rm{H}}}{\hat{ A}} + \mu p{{D}}({{\tilde{ \beta }}^n}))^{ - 1}}2{{\hat{ A}}^{\rm{H}}}{\hat{ b}}$;
     (6) 按式(17)迭代计算${{\tilde{ \beta }}^{n + 1}}$,当${ {\left\| { { { {\tilde{ \beta } } }^{n + 1} } - { { {\tilde{ \beta } } }^n} } \right\|_2^2} \Bigr/ {\left\| { { { {\tilde{ \beta } } }^n} } \right\|_2^2} } < \tau$时
       得到解${\tilde{ \beta }} = {{\tilde{ \beta }}^{n + 1}}$;
     (7) 对${\tilde{ \beta }}$进行Cameron分解,得到全极化联合多维重建结果。
    下载: 导出CSV

    表  2  仿真目标信息

    Table  2.   Information of simulated targets

    类型散射矩阵幅度位置
    三面角$\left[ {\begin{array}{*{20}{c}} 1&0 \\ 0&1 \end{array}} \right]$1$x = 1.0\,{\rm{m} },y = - 0.5\,{\rm{m} },z = 0.7\,{\rm{m} }$
    偶极子$\left[ {\begin{array}{*{20}{c}} 1&0 \\ 0&0 \end{array}} \right]$1$x = - 1.0\,{\rm{m} },y = 0.5\,{\rm{m} },z = - 0.7\,{\rm{m} }$
    30°二面角$\left[ {\begin{array}{*{20}{c}} {0.5}&{0.866} \\ {0.866}&{ - 0.5} \end{array}} \right]$1$x = - 0.5\,{\rm{m} },y = - 1.0\,{\rm{m} },z = 0.7\,{\rm{m} }$
    45°二面角$\left[ {\begin{array}{*{20}{c}} 0&1 \\ 1&0 \end{array}} \right]$1$x = 0.5\,{\rm{m} },y = 1.0\,{\rm{m} },z = - 0.7\,{\rm{m} }$
    下载: 导出CSV

    表  3  目标的仿真参数

    Table  3.   Simulation parameters of simulated targets

    雷达扫描参数
    频率范围[8 GHz, 12 GHz]
    频率采样间隔20 MHz
    方位角范围[–4°, 6°]
    方位角采样间隔1/14°
    俯仰角范围[18°, 42°]
    俯仰角采样间隔1/14°
    极化方式HH, HV, VH, VV
    下载: 导出CSV

    表  4  仿真目标的极化散射矩阵估计结果

    Table  4.   Polarization scattering matrix estimation results of simulated targets

    方法目标类型三面角偶极子30°二面角45°二面角
    联合多维重建变型前$\left[ {\begin{array}{*{20}{c}} {0.4\operatorname{j} }&0 \\ 0&{0.4{\rm{j}}} \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} { - 0.38{\rm{j}}}&0 \\ 0&0 \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} { - 0.21}&{ - 0.37} \\ { - 0.37}&{0.21} \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} 0&{ - 0.42} \\ { - 0.42}&0 \end{array}} \right]$
    变型后${\color{Blue}{0.4} }{ {\rm{e} }^{ {\rm{j} }{ {90}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } 1&0 \\ 0&1 \end{array} } \right] } }$${\color{Blue}{0.38} }{ {\rm{e} }^{ {\rm{j} }( - { {90}^ \circ })} }{\color{red}{\left[ {\begin{array}{*{20}{c} } 1&0 \\ 0&0 \end{array} } \right] } }$${\color{Blue}{0.42} }{ {\rm{e} }^{ {\rm{j} }{ {180}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } {0.5}&{0.88} \\ {0.88}&{ - 0.5} \end{array} } \right]} }$${\color{Blue}{0.42} }{ {\rm{e} }^{ {\rm{j} }{ {180}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } 0&1 \\ 1&0 \end{array} } \right]}}$
    独立多维重建变型前$\left[ {\begin{array}{*{20}{c}} {0.38{\rm{j}}}&0 \\ 0&{0.38{\rm{j}}} \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} { - 0.38{\rm{j}}}&0 \\ 0&0 \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} { - 0.13 + 0.02{\rm{j}}}&{ - 0.34 + 0.06{\rm{j}}} \\ { - 0.34 + 0.06{\rm{j}}}&{0.13 - 0.02{\rm{j}}} \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} 0&{ - 0.4 - 0.07{\rm{j}}} \\ { - 0.4 - 0.07{\rm{j}}}&0 \end{array}} \right]$
    变型后${\color{Blue}{0.38} }{ {\rm{e} }^{ {\rm{j} }{ {89}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } 1&0 \\ 0&1 \end{array} } \right]}}$${\color{Blue}{0.38} }{ {\rm{e} }^{ {\rm{j} }( - { {89}^ \circ })} }{\color{red}{\left[ {\begin{array}{*{20}{c} } 1&0 \\ 0&0 \end{array} } \right]}}$${\color{Blue}{0.26} }{ {\rm{e} }^{ {\rm{j} }{ {170}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } {0.5}&{1.35} \\ {1.35}&{ - 0.5} \end{array} } \right]}}$${\color{Blue}{0.41} }{ {\rm{e} }^{ {\rm{j} }( - { {170}^ \circ })} }{\color{red}{\left[ {\begin{array}{*{20}{c} } 0&1 \\ 1&0 \end{array} } \right]}}$
    下载: 导出CSV

    表  5  不同SNR下仿真目标的极化散射矩阵估计结果

    Table  5.   Polarization scattering matrix estimation results of simulated targets under different SNR

    目标类型SNR=13 dBSNR=18 dBSNR=23 dB


    变型
    $\left[\!\! {\begin{array}{*{20}{c}} {0.01 + 0.39{\rm{j}}} \!\!\!&\!\!\! {0.02} \\ {0.01 + 0.01{\rm{j}}} \!\!\!&\!\!\! { - 0.02 + 0.39{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.02 + 0.4{\rm{j}}} \!\!\!&\!\!\! {0.01 - 0.02{\rm{j}}} \\ {0.03{\rm{j}}} \!\!\!&\!\!\! {0.01 + 0.4{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} {0.01 + 0.41{\rm{j}}} \!\!\!&\!\!\! {0.02{\rm{j}}} \\ { - 0.01} \!\!\!&\!\!\! {0.01 + 0.40{\rm{j}}} \end{array}} \!\!\right]$
    变型
    ${\color{blue}{0.39}}{ {\rm{e} }^{ {\rm{j8} }{ {8.7}^ \circ } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! {0.02 - 0.06{\rm{j} } } \\ {0.03 - 0.04{\rm{j} } } \!\!\!&\!\!\! {1.03 + 0.07{\rm{j} } } \end{array} } \!\!\right]}}$${\color{blue}{0.4}}{ {\rm{e} }^{ {\rm{j92} }{\rm{.} }{ {\rm{2} }^{\rm{o} } } }}{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! { - 0.04 - 0.03{\rm{j} } } \\ {0.09 - 0.01{\rm{j} } } \!\!\!&\!\!\! {1.01 - 0.06{\rm{j} } } \end{array} } \!\!\right] }}$${\color{blue}{0.41} }{ {\rm{e} }^{ {\rm{j} }{ {88.6}^ \circ } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! {0.04} \\ {0.02 + 0.03{\rm{j} } } \!\!\!&\!\!\! {0.97 - 0.01{\rm{j} } } \end{array} } \!\!\right]} }$


    变型
    $\left[\!\! {\begin{array}{*{20}{c}} {0.02 - 0.38{\rm{j}}} \!\!\!&\!\!\! {0.01{\rm{j}}} \\ { - 0.02 + 0.02{\rm{j}}} \!\!\!&\!\!\! { - 0.03{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} {0.02 - 0.41{\rm{j}}} \!\!\!&\!\!\! {0.02 + 0.01{\rm{j}}} \\ { - 0.01 + 0.01{\rm{j}}} \!\!\!&\!\!\! {0.01} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.4{\rm{j}}} \!\!\!&\!\!\! {0.02{\rm{j}}} \\ 0 \!\!\!&\!\!\! 0 \end{array}} \!\!\right]$
    变型
    ${\color{blue}{0.3} }{\rm{8} }{ {\rm{e} }^{ {\rm{j( - 86} }{\rm{.} }{ {\rm{8} }^{\rm{o} } }{\rm{)} } } }\!{\color{red}{\left[\!\!\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! { - 0.03} \\ { - 0.05 \!-\! 0.05{\rm{j} } } \!\!\!&\!\!\! {0.08 \!-\! 0.13{\rm{j} } } \end{array} } \!\!\!\right]} }$$\;{\color{blue}{0.41} }{ {\rm{e} }^{ {\rm{j( - 87} }{\rm{.} }{ {\rm{3} }^{\rm{o} } }{\rm{)} } } }\!{\color{red}{\left[\!\!\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! { - 0.02 \!+\! 0.06{\rm{j} } } \\ { - 0.02 \!-\! 0.01{\rm{j} } } \!\!\!&\!\!\! { - 0.01 \!+\! 0.01{\rm{j} } } \end{array} } \!\!\!\!\right]} }$${\color{blue}{0.4} }{ {\rm{e} }^{ {\rm{j - 9} }{ {\rm{0} }^{\rm{o} } } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! { - 0.05} \\ 0 \!\!\!&\!\!\! 0 \end{array} } \!\!\right]} }$
    30°
    二面
    变型
    $\left[\!\! {\begin{array}{*{20}{c}} { - 0.16 + 0.02{\rm{j}}} \!\!\!&\!\!\! { - 0.28 + 0.04{\rm{j}}} \\ { - 0.39 + 0.07{\rm{j}}} \!\!\!&\!\!\! {0.2 - 0.04{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.22 + 0.03{\rm{j}}} \!\!\!&\!\!\! { - 0.37 + 0.06{\rm{j}}} \\ { - 0.39 + 0.06{\rm{j}}} \!\!\!&\!\!\! {0.24 - 0.04{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.22 + 0.04{\rm{j}}} \!\!\!&\!\!\! { - 0.39 + 0.06{\rm{j}}} \\ { - 0.37 + 0.05{\rm{j}}} \!\!\!&\!\!\! {0.22 - 0.04{\rm{j}}} \end{array}} \!\!\right]$
    变型
    ${\color{blue}{0.33} }{ {\rm{e} }^{ {\rm{j172} }{\rm{.} }{ {\rm{4} }^{\rm{o} } } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } {0.5} \!\!\!&\!\!\! {0.84} \\ {1.18 - 0.07{\rm{j} } } \!\!\!&\!\!\! { - 0.6 + 0.04{\rm{j} } } \end{array} } \!\!\right]} }$${\color{blue}{0.45} }{ {\rm{e} }^{ {\rm{j172} }{\rm{.} }{ {\rm{3} }^{\rm{o} } } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } {0.5} \!\!\!&\!\!\! {0.85 - 0.03{\rm{j} } } \\ {0.89 - 0.02{\rm{j} } } \!\!\!&\!\!\! { - 0.53 + 0.03{\rm{j} } } \end{array} } \!\!\right]} }$${\color{blue}{0.45} }{ {\rm{e} }^{ {\rm{j169} }{\rm{.} }{ {\rm{4} }^{\rm{o} } } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } {0.5} \!\!\!&\!\!\! {0.87 + 0.02{\rm{j} } } \\ {0.84 + 0.04{\rm{j} } } \!\!\!&\!\!\! { - 0.5} \end{array} } \!\!\right]} }$
    45°
    二面
    变型
    $\left[\!\! {\begin{array}{*{20}{c}} {0.02} \!\!\!&\!\!\! { - 0.41 - 0.06{\rm{j}}} \\ { - 0.45 - 0.06{\rm{j}}} \!\!\!&\!\!\! { - 0.02} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.01 + 0.01{\rm{j}}} \!\!\!&\!\!\! { - 0.46 - 0.06{\rm{j}}} \\ { - 0.44 - 0.08{\rm{j}}} \!\!\!&\!\!\! { - 0.01{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} 0 \!\!\!&\!\!\! { - 0.43 - 0.01{\rm{j}}} \\ { - 0.42 - 0.08{\rm{j}}} \!\!\!&\!\!\! {0.02{\rm{j}}} \end{array}} \!\!\right]$
    变型
    ${\color{blue}{0.42} }{ {\rm{e} }^{ {\rm{j( - 171} }{\rm{.} }{ {\rm{8} }^{\rm{o} } }{\rm{)} } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } { - 0.05} \!\!\!&\!\!\! 1 \\ {1.01} \!\!\!&\!\!\! {0.04} \end{array} } \!\!\right]} }$${\color{blue}{0.46} }{ {\rm{e} }^{ {\rm{j( - 172} }{\rm{.} }{ {\rm{4} }^{\rm{o} } }{\rm{)} } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } {0.02 \!-\! 0.01{\rm{j} } } \!\! & \!\! 1 \\ {0.96 \!+\! 0.05{\rm{j} } } \!\! & \!\! {0.01 \!+\! 0.02{\rm{j} } } \end{array} } \!\!\right]} }$${\color{blue}{0.44} }{ {\rm{e} }^{ {\rm{j( - 170} }{\rm{.} }{ {\rm{7} }^{\rm{o} } }{\rm{)} } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 0 \!\!\!&\!\!\! 1 \\ {0.98 \!+\! 0.01{\rm{j} } } \!\!\!&\!\!\! { - 0.02 \!-\! 0.05{\rm{j} } } \end{array} } \!\!\right]} }$
    下载: 导出CSV

    表  6  Slicy的仿真参数

    Table  6.   Simulation parameters of Slicy

    参数类型参数取值
    雷达扫描参数频率范围[8 GHz, 12 GHz]
    频率采样间隔20 MHz
    方位角范围[–4°, 6°]
    方位角采样间隔1/14°
    俯仰角范围[18°, 42°]
    俯仰角采样间隔1/14°
    极化方式HH, HV, VH, VV
    场景参数方位角0°沿x 轴正方向
    俯仰角0°沿z 轴正方向
    x 轴方向场景范围[–0.6 m, 0.6 m]
    y 轴方向场景范围[–0.9 m, 0.9 m]
    z 轴方向场景范围[0 m, 0.75 m]
    下载: 导出CSV

    表  7  卫星的仿真参数

    Table  7.   Simulation parameters of satellite

    参数类型参数取值
    雷达扫描参数频率范围[9 GHz, 11 GHz]
    频率采样间隔20 MHz
    方位角范围[–20°, 20°]
    方位角采样间隔0.1°
    俯仰角范围[45°, 65°]
    俯仰角采样间隔0.2°
    极化方式HH, HV, VH, VV
    场景参数方位角0°沿x 轴正方向
    俯仰角0°沿z 轴正方向
    x 轴方向场景范围[–0.5 m, 0.5 m]
    y 轴方向场景范围[–4.105 m, 4.105 m]
    z 轴方向场景范围[–1.775 m, 1.775 m]
    下载: 导出CSV
  • [1] 保铮, 邢孟道, 王彤. 雷达成像技术[M]. 北京: 电子工业出版社, 2005.BAO Zheng, XING Mengdao, and WANG Tong. Radar Imaging[M] Beijing: Publishing House of Electronics Industry, 2005.
    [2] CUMMING I G and WONG F H. Digital Processing of Synthetic Aperture Radar Data: Algorithm and Implementation[M]. Boston: Artech House, 2005.
    [3] LEE J S and POTTIER E. Polarimetric Radar Imaging: From Basics to Applications[M]. Boca Raton, FL: CRC Press, 2009.
    [4] 庄钊文, 肖顺平, 王雪松. 雷达极化信息处理及其应用[M]. 北京: 国防工业出版社, 1999.ZHUANG Zhaowen, XIAO Shunping, and WANG Xuesong. Radar Polarization Information Processing and Application[M] Beijing: National Defense Industry Press, 1999.
    [5] FREY O and MEIER E. Analyzing tomographic SAR data of a forest with respect to frequency, polarization, and focusing technique[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(10): 3648–3659. doi:  10.1109/TGRS.2011.2125972
    [6] GUILLASO S, FERRO-FAMIL L, REIGBER A, et al. Building characterization using L-band polarimetric interferometric SAR data[J]. IEEE Geoscience and Remote Sensing Letters, 2005, 2(3): 347–351. doi:  10.1109/LGRS.2005.851543
    [7] PONCE O, PRATS-IRAOLA P, SCHEIBER R, et al. First airborne demonstration of holographic SAR tomography with fully polarimetric multicircular acquisitions at L-band[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(10): 6170–6196. doi:  10.1109/TGRS.2016.2582959
    [8] 丁赤飚, 仇晓兰, 吴一戎. 全息合成孔径雷达的概念、体制和方法[J]. 雷达学报, 2020, 9(3): 399–408. doi:  10.12000/JR20063DING Chibiao, QIU Xiaolan, and WU Yirong. Concept, system, and method of holographic synthetic aperture radar[J]. Journal of Radars, 2020, 9(3): 399–408. doi:  10.12000/JR20063
    [9] FENG Dong, AN Daoxiang, HUANG Xiaotao, et al. A phase calibration method based on phase gradient autofocus for airborne holographic SAR imaging[J]. IEEE Geoscience and Remote Sensing Letters, 2019, 16(12): 1864–1868. doi:  10.1109/LGRS.2019.2911932
    [10] SUN Dou, XING Shiqi, LI Yongzhen, et al. Sub-aperture partitioning method for three-dimensional wide-angle synthetic aperture radar imaging with non-uniform sampling[J]. Electronics, 2019, 8(6): 629. doi:  10.3390/electronics8060629
    [11] XING Shiqi, LI Yongzhen, DAI Dahai, et al. Three-dimensional reconstruction of man-made objects using polarimetric tomographic SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(6): 3694–3705. doi:  10.1109/TGRS.2012.2220145
    [12] AUSTIN C D, ERTIN E, and MOSES R L. Sparse signal methods for 3-D radar imaging[J]. IEEE Journal of Selected Topics in Signal Processing, 2011, 5(3): 408–423. doi:  10.1109/JSTSP.2010.2090128
    [13] SUN Dou, PANG Bo, XING Shiqi, et al. Direct 3-D sparse imaging using non-uniform samples without data interpolation[J]. Electronics, 2020, 9(2): 321. doi:  10.3390/electronics9020321
    [14] ERTIN E, MOSES R L, and POTTER L C. Interferometric methods for three-dimensional target reconstruction with multipass circular SAR[J]. IET Radar, Sonar & Navigation, 2010, 4(3): 464–473.
    [15] HU Xiaowei, TONG Ningning, GUO Yiduo, et al. MIMO radar 3-D imaging based on multi-dimensional sparse recovery and signal support prior information[J]. IEEE Sensors Journal, 2018, 18(8): 3152–3162. doi:  10.1109/JSEN.2018.2810705
    [16] NANNINI M, SCHEIBER R, HORN R, et al. First 3-D reconstructions of targets hidden beneath foliage by means of polarimetric SAR tomography[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(1): 60–64. doi:  10.1109/LGRS.2011.2160329
    [17] NANNINI M, SCHEIBER R, and HORN R. Imaging of targets beneath foliage with SAR tomography[C]. The 7th European Conference on Synthetic Aperture Radar, Friedrichshafen, Germany, 2008: 1–4.
    [18] SAUER S, FERRO-FAMIL L, REIGBER A, et al. Three-dimensional imaging and scattering mechanism estimation over urban scenes using dual-baseline polarimetric InSAR observations at L-band[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(11): 4616–4629. doi:  10.1109/TGRS.2011.2147321
    [19] NGUYEN N H, BERRY P, and TRAN H T. Compressive sensing for tomographic imaging of a target with a narrowband bistatic radar[J]. Sensors, 2019, 19(24): 5515. doi:  10.3390/s19245515
    [20] ZHU Xiaoxiang and BAMLER R. Super-resolution power and robustness of compressive sensing for spectral estimation with application to spaceborne tomographic SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(1): 247–258. doi:  10.1109/TGRS.2011.2160183
    [21] CETIN M and KARL W C. Feature-enhanced synthetic aperture radar image formation based on nonquadratic regularization[J]. IEEE Transactions on Image Processing, 2001, 10(4): 623–631. doi:  10.1109/83.913596
    [22] CAMERON W L and LEUNG L K. Feature motivated polarization scattering matrix decomposition[C]. Proceedings of IEEE International Conference on Radar, Arlington, USA, 1990: 549–557.
    [23] 阙肖峰, 聂在平, 胡俊. 混合场积分方程结合MLFMA分析导体介质复合目标电磁散射问题[J]. 电子学报, 2007, 35(11): 2062–2066. doi:  10.3321/j.issn:0372-2112.2007.11.006QUE Xiaofeng, NIE Zaiping, and HU Jun. Analysis of EM scattering by composite conducting and dielectric object using combined field integral equation with MLFMA[J]. Acta Electronica Sinica, 2007, 35(11): 2062–2066. doi:  10.3321/j.issn:0372-2112.2007.11.006
  • [1] 毕辉, 张冰尘, 洪文, 吴一戎.  基于复图像的稀疏SAR成像方法在高分三号数据上的验证 . 雷达学报, 2020, 9(1): 123-130. doi: 10.12000/JR19092
    [2] 李晓峰, 张彪, 杨晓峰.  星载合成孔径雷达遥感海洋风场波浪场 . 雷达学报, 2020, 9(3): 425-443. doi: 10.12000/JR20079
    [3] 李永祯, 黄大通, 邢世其, 王雪松.  合成孔径雷达干扰技术研究综述 . 雷达学报, 2020, 9(5): 753-764. doi: 10.12000/JR20087
    [4] 黄岩, 赵博, 陶明亮, 陈展野, 洪伟.  合成孔径雷达抗干扰技术综述 . 雷达学报, 2020, 9(1): 86-106. doi: 10.12000/JR19113
    [5] 卫扬铠, 曾涛, 陈新亮, 丁泽刚, 范宇杰, 温育涵.  典型线面目标合成孔径雷达参数化成像 . 雷达学报, 2020, 9(1): 143-153. doi: 10.12000/JR19077
    [6] 田彪, 刘洋, 呼鹏江, 吴文振, 徐世友, 陈曾平.  宽带逆合成孔径雷达高分辨成像技术综述 . 雷达学报, 2020, 9(5): 765-802. doi: 10.12000/JR20060
    [7] 赵婉婉, 王鹏波, 门志荣, 李春升.  一种基于二维信号稀疏重构的互质采样星载SAR成像处理方法 . 雷达学报, 2020, 9(1): 131-142. doi: 10.12000/JR19086
    [8] 郑通, 蒋李兵, 王壮.  基于多快拍图像联合的MIMO雷达三维成像方法 . 雷达学报, 2020, 9(4): 739-752. doi: 10.12000/JR19069
    [9] 邢孟道, 林浩, 陈溅来, 孙光才, 严棒棒.  多平台合成孔径雷达成像算法综述 . 雷达学报, 2019, 8(6): 732-757. doi: 10.12000/JR19102
    [10] 高敬坤, 邓彬, 秦玉亮, 王宏强, 黎湘.  扫描MIMO阵列近场三维成像技术 . 雷达学报, 2018, 7(6): 676-684. doi: 10.12000/JR18102
    [11] 闫敏, 韦顺军, 田博坤, 张晓玲, 师君.  基于稀疏贝叶斯正则化的阵列SAR高分辨三维成像算法 . 雷达学报, 2018, 7(6): 705-716. doi: 10.12000/JR18067
    [12] 孙翔, 宋红军, 王宇, 李宁.  基于高分辨率全极化SAR图像的取向角校正方法 . 雷达学报, 2018, 7(4): 465-474. doi: 10.12000/JR18026
    [13] 金添, 宋勇平.  超宽带雷达建筑物结构稀疏成像 . 雷达学报, 2018, 7(3): 275-284. doi: 10.12000/JR18031
    [14] 祝晓静, 李飞, 王宇, 王伟, 孙翔.  基于改进方位相位编码的全极化SAR距离模糊抑制方法 . 雷达学报, 2017, 6(4): 420-431. doi: 10.12000/JR17015
    [15] 赵团, 邓云凯, 王宇, 李宁, 王翔宇.  基于扇贝效应校正的改进滑动Mosaic全孔径成像算法 . 雷达学报, 2016, 5(5): 548-557. doi: 10.12000/JR16014
    [16] 艾小锋, 曾勇虎, 高磊, 王肖洋, 汪连栋.  飞机目标全极化双基地散射特性研究 . 雷达学报, 2016, 5(6): 639-646. doi: 10.12000/JR16070
    [17] 钟金荣文贡坚, .  基于块稀疏贝叶斯学习的雷达目标压缩感知(英文) . 雷达学报, 2016, 5(1): 99-108. doi: 10.12000/JR15056
    [18] 任笑真, 杨汝良.  一种基于幅度和相位迭代重建的四维合成孔径雷达成像方法 . 雷达学报, 2016, 5(1): 65-71. doi: 10.12000/JR15135
    [19] 金添.  叶簇穿透合成孔径雷达增强成像方法 . 雷达学报, 2015, 4(5): 503-508. doi: 10.12000/JR15114
    [20] 王杰, 沈明威, 吴迪, 朱岱寅.  基于主瓣杂波高效配准的机载非正侧视阵雷达STAP 算法研究 . 雷达学报, 2014, 3(2): 235-240. doi: 10.3724/SP.J.1300.2014.13122
  • 加载中
图(8) / 表 (7)
计量
  • 文章访问数:  98
  • HTML全文浏览量:  127
  • PDF下载量:  58
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-07-06
  • 修回日期:  2020-09-24
  • 网络出版日期:  2020-10-15
  • 刊出日期:  2020-10-28

基于稀疏重构的全极化SAR联合多维重建

doi: 10.12000/JR20092
    基金项目:  国家自然科学基金(61971429, 61901499)
    作者简介:

    孙 豆(1992–),女,博士研究生,主要研究方向为极化雷达成像和雷达信号处理。E-mail: sundou14@nudt.edu.cn

    路东伟(1992–),男,博士研究生,主要研究方向为合成孔径雷达对抗和雷达目标识别。E-mail: bookwormldw@qq.com

    邢世其(1984–),男,副研究员,主要研究方向为极化雷达成像、雷达信号处理以及合成孔径雷达对抗。E-mail: xingshiqi_paper@163.com

    杨 潇(1983–),男,助教,主要研究方向为雷达信号处理和雷达目标识别。E-mail: 297414430@qq.com

    李永祯(1977–),男,研究员,博士生导师,主要研究方向为极化雷达与电子对抗。E-mail: e0061@sina.com

    王雪松(1972–),男,教授,博士生导师,主要研究方向为极化雷达、目标识别与电子对抗。E-mail: wxs1019@vip.sina.com

    通讯作者: 邢世其 xingshiqi_paper@163.com
  • 责任主编:仇晓兰 Corresponding Editor: QIU Xiaolan
  • 中图分类号: TN95

摘要: 各极化通道独立处理和三维分步成像会忽视数据之间的关联性,造成散射中心的失配以及极化散射矩阵获取的不准确。鉴于此,该文提出一种基于稀疏重构的全极化联合多维重建方法。该方法通过设置联合稀疏约束对所有极化通道及所有维度进行联合,将全极化多维重建建模为多通道联合稀疏重构问题。通过数据插值对模型简化后,结合三维快速傅里叶变换、共轭梯度法和牛顿迭代法给出一种高效的模型求解方法,可以同时得到极化散射矩阵和目标三维信息。该文方法保证了不同极化通道、不同维度的稀疏支撑集一致,且充分利用了数据之间的关联性带来的额外信息。基于仿真数据和电磁计算数据的实验结果表明,该方法的性能不受目标类型影响,具有一定的抗噪性,能有效地获取目标的多维重建结果,得到的三维成像结果分辨率高且极化散射矩阵估计精度高。

注释:
1)  责任主编:仇晓兰 Corresponding Editor: QIU Xiaolan

English Abstract

孙豆, 路东伟, 邢世其, 等. 基于稀疏重构的全极化SAR联合多维重建[J]. 雷达学报, 2020, 9(5): 865–877. doi:  10.12000/JR20092
引用本文: 孙豆, 路东伟, 邢世其, 等. 基于稀疏重构的全极化SAR联合多维重建[J]. 雷达学报, 2020, 9(5): 865–877. doi:  10.12000/JR20092
SUN Dou, LU Dongwei, XING Shiqi, et al. Full-polarization SAR joint multidimensional reconstruction based on sparse reconstruction[J]. Journal of Radars, 2020, 9(5): 865–877. doi:  10.12000/JR20092
Citation: SUN Dou, LU Dongwei, XING Shiqi, et al. Full-polarization SAR joint multidimensional reconstruction based on sparse reconstruction [J]. Journal of Radars, 2020, 9(5): 865–877. doi:  10.12000/JR20092
    • 合成孔径雷达(Synthetic Aperture Radar, SAR)具有全天时全天候对目标场景进行观测的优势,日益成为目标识别的关键传感器[1,2]。极化是电磁波的一个重要本质属性,深刻反映目标形状、结构和姿态等物理特征,对于提升目标识别能力具有重要作用[3,4]。近年来,SAR越来越朝着全极化和三维的方向发展[5-8]。三维信息使得雷达对目标的描述越来越精细,极化使得雷达对目标散射特性的描述更加完整,将全极化信息和三维散射信息相结合可以得到人造目标更为全面的重建结果,有利于提高目标分类识别精度。

      梳理现有的文献资料,将极化三维重建相关的研究工作进行分类。按照成像方法的不同,相关研究可以分为两类。一类是基于傅里叶变换的成像方法[9,10],这类方法是传统二维SAR成像方法的三维扩展,具有快捷高效的优势,但是这类方法只能单通道独立成像,无法进行多通道联合同时成像。另一类是基于稀疏重构的成像方法[11-13],由于稀疏成像本质是一种基于模型匹配的方法,因此这类方法的成像结果分辨率高,且对旁瓣和噪声的抑制效果好,此外,稀疏重构类的方法可以进行多通道联合处理。按照成像方式的不同,相关研究可以分为两类。一类是先进行二维成像再获取高度维信息[9,11,14],这类方法可以使用现有的二维成像方法,只需进一步对高度信息进行估计,其优势在于计算复杂度低,效率高,但是成像结果受二维图像配准质量的影响。另一类是直接三维成像[10,13,15],即同时得到目标的三维信息,这类方法的优势在于可以保证数据之间的关联性,不需要进行额外的图像配准,但是直接三维成像会面临巨大的计算压力。按照极化信息的提取方式的不同,相关研究也可以分为两类。一类是各极化通道独立处理[16,17],即先每个通道分别成像,再提取极化散射矩阵。这种方式的优势在于快捷,适用性广,但是由于忽略了各个极化通道数据之间的关联性,不能保证提取的极化散射矩阵准确。另一类是所有极化通道联合处理[11,18],即所有通道联合成像,同时得到极化散射矩阵,这种联合的方式保证了极化信息提取的准确度,可以有效克服独立处理中散射中心失配的问题,但是联合处理的计算复杂度较大。

      现有文献所采用的方法都是对上述3种分类的不同组合,但是还没有方法将稀疏成像、直接三维成像和极化联合处理这三者相结合。实际上,相比于基于傅里叶变换的成像,稀疏成像在分辨率、抗噪性等方面有着更多的优势[19,20]。相比于直接三维成像,先二维后一维的成像方法很可能破坏散射的一致性,特别是每个二维成像结果中散射点的数目和位置估计可能不同,这种不一致性会使得高度信息的提取不准确。此外,由于每个极化通道观测相同的区域,各通道数据高度相关,且之间存在着丰富的互补信息,三维信息和极化信息同时联合获取不破坏数据之间的关联性,才能充分利用这些互补信息来提高估计精度和分辨率。因此,考虑到稀疏成像、直接三维成像和极化联合处理这三种方式各自的优势,将它们结合进行极化多维重建可以得到质量、精度更高的重建结果。

      为了充分利用数据之间的关联性带来的额外信息,实现极化散射矩阵和目标三维信息的同时获取,本文提出基于稀疏重构的全极化联合多维重建方法。首先建立各极化通道的直接三维稀疏重建模型,然后对所有通道设置联合稀疏约束,将全极化三维重建建模为多通道联合稀疏重构问题。考虑到上述稀疏重建问题规模大、计算压力大,本文通过数据插值对模型进行简化,并结合三维快速傅里叶变换、共轭梯度法和牛顿迭代法给出了一种高效的模型求解方法。最后,对三维极化散射矩阵进行极化分解,得到包含目标三维位置信息以及散射类型的多维重建结果。本文提出的全极化联合多维重建方法保证了所有极化通道中散射中心的数目和三维位置一致,有效地克服了散射中心失配的问题,确保了极化特征提取的准确性。此外,该方法大大降低了多维联合处理面临的计算压力,能够在不改变结果精度的情况下有效地获取目标多维重建结果。基于仿真数据和电磁计算数据的实验验证了本文方法的有效性。

    • 假设雷达距离场景足够远,位于相对于场景中心方位角$\phi $和俯仰角$\theta $的位置,并发射中心频率为${f_{\rm{c}}}$,带宽为${\rm{BW}}$的观测信号,使用平面波模型,接收信号可以表示为

      $$ r(\tilde x;\phi ,\theta ) = \left[ {\iint {g\left(\tilde x = \frac{{{\rm{c}}t}}{2},\tilde y,\tilde z;\phi ,\theta \right)}{\rm{d}}\tilde y{\rm{d}}\tilde z} \right] * s(t) $$ (1)

      其中,${\rm{c}}$是光速,$t$表示时间,$s(t)$是已知的观测信号,$ * $表示卷积。$g\left(\tilde x = {{({\rm{c}}t)}}/{2},\tilde y,\tilde z;\phi ,\theta \right)$为场景的反射率函数,$\tilde x,\tilde y,\tilde z$表示目标在场景中的位置。

      式(1)对$\tilde x$做傅里叶变换,并匹配滤波,得到

      $$ R(k;\phi ,\theta ) \!=\!\! \int\! {\left[ {\iint {g\left(\tilde x = \frac{{{\rm{c}}t}}{2},\tilde y,\tilde z;\phi ,\theta \right)}{\rm{d}}\tilde y{\rm{d}}\tilde z} \right] \!\cdot\! {{\rm{e}}^{ - {\rm{j}}\tilde xk}}} {\rm{d}}\tilde x $$ (2)

      其中,波数$k = \dfrac{{4{\rm{\pi }}f}}{{\rm{c}}}$, $f$表示观测信号的频率,且$f \in \left[ {{f_{\rm{c}}} - \dfrac{{{\rm{BW}}}}{2},{f_{\rm{c}}} + \dfrac{{{\rm{BW}}}}{2}} \right]$

      式(2)可以理解为$g\left(\tilde x = \dfrac{{{\rm{c}}t}}{2},\tilde y,\tilde z;\phi ,\theta \right)$$\tilde x$方向上投影$\displaystyle\iint {g\left(\tilde x = \frac{{{\rm{c}}t}}{2},\tilde y,\tilde z;\phi ,\theta \right)}{\rm{d}}\tilde y{\rm{d}}\tilde z$的傅里叶变换,因此根据投影切片定理[12],可得

      $$ \begin{split} & \int \left[ {\iint {g\left(\tilde x = \frac{{{\rm{c}}t}}{2},\tilde y,\tilde z;\phi ,\theta \right)}{\rm{d}}\tilde y{\rm{d}}\tilde z} \right] \\ & \qquad \cdot {{\rm{e}}^{ - {\rm{j}}\tilde xk}} {\rm{d}}\tilde x = G\left( {k,0,0} \right) \end{split} $$ (3)

      其中,$G\left( {k,0,0} \right)$为反射率函数的三维傅里叶变换在$(k,0,0)$处的值。

      以角度$(\phi ,\theta )$$(\tilde x,\tilde y,\tilde z)$$\left( {k,0,0} \right)$旋转,得到旋转后的$(x,y,z)$$({k_x},{k_y},{k_z})$。根据傅里叶变换的旋转特性,式(3)变为

      $$ \begin{split} G({k_x},{k_y},{k_z}) =\,& \int {g(x,y,z;\phi ,\theta )}\\ & \cdot {{\rm{e}}^{ - {\rm{j}}({k_x}x + {k_y}y + {k_z}z)}}{\rm{d}}x{\rm{d}}y{\rm{d}}z \end{split} $$ (4)

      其中

      $$ \left. \begin{aligned} & {k_x^{} = k\cos \theta \cos \phi } \\ & {k_y^{} = k\cos \theta \sin \phi } \\ & {k_z^{} = k\sin \theta } \end{aligned} \right\} $$ (5)

      由式(4)可知,反射率函数是波数域采样$G({k_x}, {k_y},{k_z})$的三维傅里叶逆变换。

    • 对于每个极化通道,需要被重建的成像场景可以通过式(6)获得

      $$ \begin{split} {G_l}(k_x^{},k_y^{},k_z^{}) =\,& \int {{g_l}(x,y,z;\phi ,\theta )}\\ & \cdot {{\rm{e}}^{ - {\rm{j}}(k_x^{}x + k_y^{}y + k_z^{}z)}}{\rm{d}}x{\rm{d}}y{\rm{d}}z \end{split} $$ (6)

      其中,${g_l}(x,y,z;\phi ,\theta )$为第$l\left( {l = 1,2, ··· ,L} \right)$个极化通道的场景反射率函数,${G_l}(k_x^{},k_y^{},k_z^{})$是第$l$个极化通道的波数域采样。

      在三维重建空间中,定义$N$个位置为候选散射中心

      $${{C}} = \{ ({x_n},{y_n},{z_n})\} _{n = 1}^N$$ (7)

      通常,这些位置是从均匀的矩形网格中选择的。基于这些位置,定义$M \times N$的字典矩阵

      $$ {{A}} = {\left[ {{{\rm{e}}^{ - {\rm{j}}(k_{x,m}^{}{x_n} + k_{y,m}^{}{y_n} + k_{z,m}^{}{z_n})}}} \right]_{m,n}} $$ (8)

      其中,$m$表示波数域采样的索引,$n$表示${{C}}$$N$个位置的索引。

      根据字典矩阵${{A}}$,将式(6)写成矩阵形式,可得

      $$ {{{b}}_l} = {{A}}{{{\beta}} _l} $$ (9)

      其中,$N \times 1$的向量${{{\beta}} _l}$是要重建的第$l$个极化通道的三维场景,$M \times 1$的向量${{{b}}_l}$表示第$l$个极化通道的波数域采样。

      每个极化通道单独三维重建就是要解决下面的直接三维稀疏成像模型

      $$ {\tilde {{\beta}} _l} = \mathop {{\rm{argmin}}}\limits_{{{{\beta}} _l}} {\left\| {{{{\beta}} _l}} \right\|_0}{\rm{, s}}{\rm{.t}}{\rm{. }}\;\;{{{b}}_l} = {{A}}{{{\beta}} _l} $$ (10)

      对式(10)中的直接三维稀疏成像模型进行求解则可以同时得到目标的三维信息。相比于先二维稀疏成像再第三维稀疏重建的方法,直接三维稀疏成像避免了二维图像稀疏支撑集的不一致。

      对于每个极化通道,由于其观测场景是相同的,且采样方式不变,因此,各极化通道共享相同的字典矩阵${{A}}$。基于此,建立联合成像模型

      $$ \tilde {{\beta}} = \mathop {{\rm{argmin}}}\limits_{{\beta }} {\left\| {{\beta }} \right\|_0},\;{\rm{ s}}{\rm{.t}}{\rm{. }}\;\;\;{{{b}}_l} = {{A}}{{{\beta }}_l}\left( {l = 1,2 ,···, L} \right)\!\!\! $$ (11)

      其中${{\beta }} = \left[ {{{{\beta }}_1}}\ ··· \ {{{{\beta }}_l}}\ ···\ {{{{\beta }}_L}} \right]$是联合稀疏解矩阵,${\left\| {{\beta }} \right\|_0}{\rm{ = }} {\left\| {\left| {{{{\beta }}_1}} \right| + ··· + \left| {{{{\beta }}_l}} \right| + ··· + \left| {{{{\beta }}_L}} \right|} \right\|_0}$表示稀疏解矩阵${{\beta }}$中不为零的行数。

      在联合成像模型中,每个极化通道中${{{\beta }}_l}$的支撑集是相同的,只是不同极化通道的${{{\beta }}_l}$的幅度不同。联合稀疏约束项${\left\| {{\beta }} \right\|_0}{\rm{ = }}{\left\| {\left| {{{{\beta }}_1}} \right| + ···+ \left| {{{{\beta }}_l}} \right| + ··· + \left| {{{{\beta }}_L}} \right|} \right\|_0}$保证了不同极化通道的重建结果共享相同数量和位置的散射中心。联合稀疏为散射系数的估计提供了额外的约束条件,使得在重建过程中可以保留更多有价值的信息。

    • 式(11)中的数学模型不是一个凸优化问题,而是一个NP难问题。为了用数值方法求解该模型,首先对目标函数进行放松,定义混合规范

      $${\left( {{{\left\| {{\beta }} \right\|}_{2,p}}} \right)^p} = \sum\limits_i {\left( {{{\left( {\sum\limits_l {\left( {{{\left| {{{{\beta }}_l}(i)} \right|}^2}} \right)} } \right)}^{p/2}}} \right)} $$ (12)

      其中${{{\beta }}_l}(i)$表示${{{\beta }}_l}$中的第$i$个元素。混合范数${\ell _{2,p}}$沿不同极化通道的方向计算${\ell _2}$范数,即计算${{\beta }}$的张量,并沿散射位置的方向计算${\ell _p}$范数。

      结合式(12)中的混合范数,式(11)中的模型被松弛为式(13)的优化问题

      $${\tilde{ \beta }} = \arg \mathop {\min }\limits_{{\beta }} \left( {\sum\limits_l {\left\| {{{{b}}_l} - {{A}}{{{\beta }}_l}} \right\|_2^2} {\rm{ + }}\mu {{\left( {{{\left\| {{\beta }} \right\|}_{2,p}}} \right)}^p}} \right)$$ (13)

      其中$\mu $是控制稀疏性的正则化参数。联合多维重建就是对式(13)中的模型进行求解。

      为了求解式(13),本文在Cetin等人[21]处理思想的基础上,提出了一种快速迭代方法。定义代价函数

      $$J = \sum\limits_l {\left\| {{{{b}}_l} - {{A}}{{{\beta }}_l}} \right\|_2^2} {\rm{ + }}\mu {\left( {{{\left\| {{\beta }} \right\|}_{2,p}}} \right)^p}$$ (14)

      当代价函数$J$取最小值时,可以得到稀疏重建的结果。计算$J$${{\beta }}$的偏导数,得到

      $$ {\nabla _{{\beta }}}J = \Bigr(2{{{A}}^{\rm{H}}}{{A}} + \mu p{{D}}({{\beta }})\Bigr){{\beta }} - 2{{{A}}^{\rm{H}}}{{b}} $$ (15)

      其中${{D}}({{\beta }}) = {\rm{diag}}\left( {{{\left( {\displaystyle\sum\nolimits_l {{{\left| {{{{\beta }}_l}\left( i \right)} \right|}^2}} } \right)}^{p/2 - 1}}} \right)$$N \times N$的对角矩阵,${{b}} = \left[ {{{{b}}_1}} \ ··· \ {{{{b}}_l}}\ ··· \ {{{{b}}_L}} \right]$

      ${\nabla _{{\beta }}}J = {{0}}$时,则得到了散射系数的估计。考虑到式(15)中的$2{{{A}}^{\rm{H}}}{{A}} + \mu p{{D}}({{\beta }})$可近似看作Hessian矩阵,因此可以采用近似高斯迭代法求解${\nabla _{{\beta }}}J = {{0}}$,得到迭代式:

      $$ {{\tilde{ \beta }}^{n + 1}} \!=\! {{\tilde{ \beta }}^n} - {\varDelta ^{n + 1}}\Bigr({{\tilde{ \beta }}^n} \!-\! {\Bigr(2{{{A}}^{\rm{H}}}{{A}} \!+\! \mu p{{D}}({{\tilde{ \beta }}^n})\Bigr)^{ - 1}}2{{{A}}^{\rm{H}}}{{b}}\Bigr) $$ (16)

      其中${\varDelta ^{n + 1}}$是在迭代中的步长,为了使算法快速收敛,在迭代过程中采用变步长${\varDelta ^{n + 1}} = {\left( {{\varDelta ^n}} \right)^{0.9}}$。当${{\left\| {{{{\tilde{ \beta }}}^{n + 1}} - {{{\tilde{ \beta }}}^n}} \right\|_2^2} \Bigr/ {\left\| {{{{\tilde{ \beta }}}^n}} \right\|_2^2}}$小于预设阈值时,算法退出,此时就得到了所有极化通道的三维重建结果。为了进一步提取目标极化特性,还可对式(16)得到的极化散射矩阵进行Cameron分解[22],最终就得到了包含极化信息和目标三维信息的多维重建结果。

      实际上,式(16)的迭代计算还存在一些问题。由于本文方法进行三维直接重建,同时获取目标的全极化三维信息,不可避免地面临计算量巨大的问题,即式(16)中的字典矩阵维度很大,无法直接存储计算。具体地,字典矩阵${{A}}$$M \times N$的,$N$与场景的大小和分辨率有关, $M$与波数域的采样点数有关。假设要重建的场景大小为$5\;{\rm{m}} \times 5\;{\rm{m}} \times 5\;{\rm{m}}$,每个维度的分辨率是0.05 m,那么三维重建结果的维度是$100 \times 100 \times 100$,则$N = 1 \times {10^6}$。假设频率的采样点数为100,方位向的采样点数为100,俯仰向的采样点数为100,则$M = 1 \times {10^6}$。在这种情况下,由于$M \times N = 1 \times {10^{12}}$,字典矩阵${{A}}$需要的存储空间大小约为3GB。除了要存储巨大的字典矩阵之外,每次迭代都不可避免地要进行字典矩阵的乘法运算,如计算${{A}}{{{\beta }}_l}$,需要$M \times N = 1 \times {10^{12}}$次乘法和加法计算,因此直接对式(16)进行迭代的计算复杂度很高。为了解决上述计算规模大的问题,本文借鉴文献[12]中的思路,结合三维非均匀快速傅里叶变换(3-D Non-Uniform Fast Fourier Transform, 3-D NUFFT)和共轭梯度法求解式(16)中的${(2{{{A}}^{\rm{H}}}{{A}} + \mu p{{D}}({{\tilde{ \beta }}^n}))^{ - 1}}2{{{A}}^{\rm{H}}}{{b}}$部分。

      由式(5)可知,${k_x},{k_y},{k_z}$$f,\theta ,\phi $之间的关系不是线性的,因此${k_x},{k_y},{k_z}$处于非均匀网格上。对于字典矩阵${{A}} = {\left[ {{{\rm{e}}^{ - {\rm{j}}(k_{x,m}^{}{x_n} + k_{y,m}^{}{y_n} + k_{z,m}^{}{z_n})}}} \right]_{m,n}}$,虽然$x,y,z$处于均匀的网格上,但由于${k_x},{k_y},{k_z}$是非均匀的,因此字典矩阵不可以被看作三维傅里叶变换矩阵。为了可以使用三维快速傅里叶变换(3-D Fast Fourier Transform, 3-D FFT)代替字典矩阵,就需要在稀疏重构前对每个极化通道的波数域采样进行预处理,使得波数域采样处于均匀的网格上。

      这里使用3-D NUFFT对非均匀网格${k_x},{k_y},{k_z}$上的每个极化通道的波数域采样${G_l}(k_x^{},k_y^{},k_z^{})$进行插值,插值后的波数域采样${\hat G_l}(\hat k_x^{},\hat k_y^{},\hat k_z^{})$处于均匀网格$\hat k_x^{},\hat k_y^{},\hat k_z^{}$上。对应于插值后的波数域采样,字典矩阵则更新为${\hat{ A}} = {\left[ {{{\rm{e}}^{ - {\rm{j}}(\hat k_{x,m}^{}{x_n} + \hat k_{y,m}^{}{y_n} + \hat k_{z,m}^{}{z_n})}}} \right]_{m,n}}$,式(16)更新为

      $$ {{\tilde{ \beta }}^{n + 1}} \!=\! {{\tilde{ \beta }}^n} - {\varDelta ^{n + 1}}\Bigr({{\tilde{ \beta }}^n} - {\Bigr(2{{\hat{ A}}^{\rm{H}}}{\hat{ A}} + \mu p{{D}}({{\tilde{ \beta }}^n})\Bigr)^{ - 1}}\!2{{\hat{ A}}^{\rm{H}}}{\hat{ b}}\Bigr) $$ (17)

      其中${\hat{ b}} = \left[ {{{{\hat{ b}}}_1}} \ ··· \ {{{{\hat{ b}}}_l}} \ ··· \ {{{{\hat{ b}}}_L}} \right]$, ${{\hat{ b}}_l}$根据${\hat G_l}(\hat k_x^{},\hat k_y^{},\hat k_z^{})$得到。由于$\hat k_x^{},\hat k_y^{},\hat k_z^{}$是均匀网格,更新后的字典矩阵${\hat{ A}}$和三维傅里叶变换矩阵是等价的。因此,可以在不存储字典矩阵的情况下,使用三维快速傅里叶逆变换(3-D Fast Inverse Fourier Transform, 3-D IFFT)计算${{\hat{ A}}^{\rm{H}}}{\hat{ b}}$。根据得到的${{\hat{ A}}^{\rm{H}}}{\hat{ b}}$,可以用共轭梯度法求解${(2{{\hat{ A}}^{\rm{H}}}{\hat{ A}} + \mu p{{D}}({{\tilde{ \beta }}^n}))^{ - 1}}2{{\hat{ A}}^{\rm{H}}}{\hat{ b}}$,其中${\hat{ A}}$${{\hat{ A}}^{\rm{H}}}$使用3-D FFT和3-D IFFT代替,需要的计算量约为${D^3}{\rm{ln}} \left( {{D^3}} \right)$,其中$D$是成像场景中各维采样点数最大值,这里$D = 100$。相比较于直接进行字典矩阵的乘法计算需要$1 \times {10^{12}}$次计算,使用3-D FFT或3-D IFFT只需要$1.38 \times {10^7}$次计算,节约了$7.2 \times {10^4}$倍,这就解决了重建中计算规模大的困难。

      表1总结了本文所提出的全极化联合多维重建方法的具体步骤。

      表 1  全极化联合多维重建方法的步骤

      Table 1.  Steps of full polarization joint multi-dimensional reconstruction method

       (1) 设定初值:${{\tilde{ \beta }}^n} = {{0}}$;
       (2) 对每个极化通道,
         (a) 使用3-D NUFFT对${G_l}(k_x^{},k_y^{},k_z^{})$进行插值,得到
           ${\hat G_l}(\hat k_x^{},\hat k_y^{},\hat k_z^{})$,
         (b) 对${\hat G_l}(\hat k_x^{},\hat k_y^{},\hat k_z^{})$进行3-D IFFT,并向量化结果得到
           ${{\hat{ A}}^{\rm{H}}}{{\hat{ b}}_l}$;
       (3) 根据所有极化通道的${{\hat{ A}}^{\rm{H}}}{{\hat{ b}}_l}$,得到${{\hat{ A}}^{\rm{H}}}{\hat{ b}}$;
       (4) 结合3-D FFT,3-D IFFT和共轭梯度法计算
         ${(2{{\hat{ A}}^{\rm{H}}}{\hat{ A}} + \mu p{{D}}({{\tilde{ \beta }}^n}))^{ - 1}}$;
       (5) 根据${{\hat{ A}}^{\rm{H}}}{\hat{ b}}$,计算${(2{{\hat{ A}}^{\rm{H}}}{\hat{ A}} + \mu p{{D}}({{\tilde{ \beta }}^n}))^{ - 1}}2{{\hat{ A}}^{\rm{H}}}{\hat{ b}}$;
       (6) 按式(17)迭代计算${{\tilde{ \beta }}^{n + 1}}$,当${ {\left\| { { { {\tilde{ \beta } } }^{n + 1} } - { { {\tilde{ \beta } } }^n} } \right\|_2^2} \Bigr/ {\left\| { { { {\tilde{ \beta } } }^n} } \right\|_2^2} } < \tau$时
         得到解${\tilde{ \beta }} = {{\tilde{ \beta }}^{n + 1}}$;
       (7) 对${\tilde{ \beta }}$进行Cameron分解,得到全极化联合多维重建结果。

      该联合多维重建方法实现了极化散射矩阵和目标三维信息的同时获取,充分利用了数据之间的关联性带来的额外信息。实际上,本文方法的关键在于反映了联合稀疏性的${{D}}({{\tilde{ \beta }}^n})$矩阵包含了所有极化通道和全部三维的散射系数。联合稀疏性为散射系数的估计提供了额外的约束条件,保证了不同极化通道、不同维度的稀疏支撑集是一致的。不管散射中心是不是在所有极化通道都有散射信息,联合稀疏对4个极化通道进行联合,同时估计散射中心在各个极化通道的散射信息,有效地克服了独立处理和三维分步成像中散射中心失配的问题,进而保证了获取的极化散射矩阵的准确性。此外,同时获取极化和三维信息不可避免地要面对计算规模大的问题。本文方法结合3-D FFT和共轭梯度法对模型进行简化,避免了巨大字典矩阵的存储,大大降低了计算压力,使联合多维重建可以在不改变解的精度的前提下高效地进行。

    • 本节基于仿真数据和电磁计算数据开展实验,验证本文提出方法的有效性。实验在一台配备Intel Core I5-6500 CPU和12 GB RAM的计算机上的MATLAB R2016b中进行。

    • 选取三面角,偶极子,30°二面角和45°二面角作为仿真目标,进行全极化联合三维重建,并从目标类型依赖性、极化散射矩阵估计精度、噪声敏感度这3个方面分析本文方法的性能。表2给出了4个仿真目标的信息,其中三面角,偶极子和45°二面角的稀疏支撑集不一致,30°二面角的稀疏支撑集一致。回波数据由MATLAB生成,使用理想点散射模型,表3记录了仿真实验时的各项参数。

      表 2  仿真目标信息

      Table 2.  Information of simulated targets

      类型散射矩阵幅度位置
      三面角$\left[ {\begin{array}{*{20}{c}} 1&0 \\ 0&1 \end{array}} \right]$1$x = 1.0\,{\rm{m} },y = - 0.5\,{\rm{m} },z = 0.7\,{\rm{m} }$
      偶极子$\left[ {\begin{array}{*{20}{c}} 1&0 \\ 0&0 \end{array}} \right]$1$x = - 1.0\,{\rm{m} },y = 0.5\,{\rm{m} },z = - 0.7\,{\rm{m} }$
      30°二面角$\left[ {\begin{array}{*{20}{c}} {0.5}&{0.866} \\ {0.866}&{ - 0.5} \end{array}} \right]$1$x = - 0.5\,{\rm{m} },y = - 1.0\,{\rm{m} },z = 0.7\,{\rm{m} }$
      45°二面角$\left[ {\begin{array}{*{20}{c}} 0&1 \\ 1&0 \end{array}} \right]$1$x = 0.5\,{\rm{m} },y = 1.0\,{\rm{m} },z = - 0.7\,{\rm{m} }$

      表 3  目标的仿真参数

      Table 3.  Simulation parameters of simulated targets

      雷达扫描参数
      频率范围[8 GHz, 12 GHz]
      频率采样间隔20 MHz
      方位角范围[–4°, 6°]
      方位角采样间隔1/14°
      俯仰角范围[18°, 42°]
      俯仰角采样间隔1/14°
      极化方式HH, HV, VH, VV

      图1给出了仿真目标的4个极化通道的成像结果,其中左列是HH, HV, VH和VV极化通道下的二维投影结果,右列是HH, HV, VH和VV极化通道下的三维结果。以HH极化通道成像结果的峰值为0 dB,图1中所有子图显示成像结果的幅度阈值为–20 dB。观察HH通道的成像结果有3个点目标,对应三面角、偶极子、30°二面角;HV通道和VH通道的成像结果有两个点目标,对应30°二面角和45°二面角;VV通道的成像结果有两个点目标,对应三面角和30°二面角。各个极化通道的成像结果和目标的极化散射矩阵可以对应上,目标极化散射矩阵分量是0的通道的成像结果没有该目标。

      图  1  各个极化通道的仿真目标成像结果

      Figure 1.  Each polarization channel’s imaging results of simulated targets

      图2给出了仿真目标的多维重建结果,其中图2(a)为全极化三维成像结果,图2(b)为Cameron分解结果。对比图2(a)的全极化三维成像结果和图1中的各个极化通道的三维成像结果,可以看出图2(a)融合了所有极化通道的结果,包含了全部的四个目标,对应三面角、偶极子、30°二面角和45°二面角。图2(b)中不同的图标表示由Cameron分解得到的不同的散射类别。对比图2(b)中的结果和表2中的目标信息可知,各散射类型判别结果与实际情况相符合,验证了本文方法对极化散射矩阵提取的有效性。

      图  2  仿真目标全极化联合多维重建结果

      Figure 2.  Full polarization joint multi-dimensional reconstruction results of simulated targets

      为了进一步分析本文提出的联合多维重建方法对极化散射矩阵的估计精度,表4给出了仿真目标的极化散射矩阵估计结果,并同时列出了各个极化通道分别独立多维重建的估计结果进行对比。变型前的结果是多维重建得到的极化散射矩阵原始值。为了和目标真实的极化散射矩阵进行对比分析,我们对估计得到的每个散射矩阵进行归一化处理,将幅度和相位提到矩阵外,得到变形后的极化散射矩阵估计值。对于变形后的极化散射矩阵估计值,蓝色部分表示各个目标的幅度估计结果,红色部分表示各个目标的极化散射矩阵估计结果,指数部分表示与当前场景有关的相对相位,不影响结果的分析。从表4中联合多维重建方法得到的变型后的结果可以看出,各个目标的散射幅度相近,这与表2中目标真实的幅度情况一致。由于NUFFT插值和3-D FFT会带来增益的变化,因此幅度的估计结果不是1。与表2中目标真实的极化散射矩阵对比,可以看出对于三面角、偶极子和45°二面角,联合多维重建方法得到的极化散射矩阵的估计值和真实值完全一致,对于30°二面角,联合多维重建方法得到的极化散射矩阵估计结果基本上和真实值一致,交叉极化通道的估计结果略微大于真实值。观察表4中独立多维重建法得到的变型后的结果,可以看出极化散射矩阵的估计结果并不准确。具体地,30°二面角的散射幅度和其他目标的散射幅度差异很大,这与表2中目标真实的幅度情况不一致。此外,对于30°二面角,独立多维重建方法得到的极化散射矩阵估计结果与真实值差异很大,交叉极化通道的估计结果远大于真实值。因此,各个极化通道分别独立多维重建不能保证极化散射矩阵提取的准确性,而联合多维重建方法则可以得到较准确的极化散射矩阵的估计结果。

      表 4  仿真目标的极化散射矩阵估计结果

      Table 4.  Polarization scattering matrix estimation results of simulated targets

      方法目标类型三面角偶极子30°二面角45°二面角
      联合多维重建变型前$\left[ {\begin{array}{*{20}{c}} {0.4\operatorname{j} }&0 \\ 0&{0.4{\rm{j}}} \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} { - 0.38{\rm{j}}}&0 \\ 0&0 \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} { - 0.21}&{ - 0.37} \\ { - 0.37}&{0.21} \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} 0&{ - 0.42} \\ { - 0.42}&0 \end{array}} \right]$
      变型后${\color{Blue}{0.4} }{ {\rm{e} }^{ {\rm{j} }{ {90}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } 1&0 \\ 0&1 \end{array} } \right] } }$${\color{Blue}{0.38} }{ {\rm{e} }^{ {\rm{j} }( - { {90}^ \circ })} }{\color{red}{\left[ {\begin{array}{*{20}{c} } 1&0 \\ 0&0 \end{array} } \right] } }$${\color{Blue}{0.42} }{ {\rm{e} }^{ {\rm{j} }{ {180}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } {0.5}&{0.88} \\ {0.88}&{ - 0.5} \end{array} } \right]} }$${\color{Blue}{0.42} }{ {\rm{e} }^{ {\rm{j} }{ {180}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } 0&1 \\ 1&0 \end{array} } \right]}}$
      独立多维重建变型前$\left[ {\begin{array}{*{20}{c}} {0.38{\rm{j}}}&0 \\ 0&{0.38{\rm{j}}} \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} { - 0.38{\rm{j}}}&0 \\ 0&0 \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} { - 0.13 + 0.02{\rm{j}}}&{ - 0.34 + 0.06{\rm{j}}} \\ { - 0.34 + 0.06{\rm{j}}}&{0.13 - 0.02{\rm{j}}} \end{array}} \right]$$\left[ {\begin{array}{*{20}{c}} 0&{ - 0.4 - 0.07{\rm{j}}} \\ { - 0.4 - 0.07{\rm{j}}}&0 \end{array}} \right]$
      变型后${\color{Blue}{0.38} }{ {\rm{e} }^{ {\rm{j} }{ {89}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } 1&0 \\ 0&1 \end{array} } \right]}}$${\color{Blue}{0.38} }{ {\rm{e} }^{ {\rm{j} }( - { {89}^ \circ })} }{\color{red}{\left[ {\begin{array}{*{20}{c} } 1&0 \\ 0&0 \end{array} } \right]}}$${\color{Blue}{0.26} }{ {\rm{e} }^{ {\rm{j} }{ {170}^ \circ } } }{\color{red}{\left[ {\begin{array}{*{20}{c} } {0.5}&{1.35} \\ {1.35}&{ - 0.5} \end{array} } \right]}}$${\color{Blue}{0.41} }{ {\rm{e} }^{ {\rm{j} }( - { {170}^ \circ })} }{\color{red}{\left[ {\begin{array}{*{20}{c} } 0&1 \\ 1&0 \end{array} } \right]}}$

      为分析本文方法对噪声的敏感度,仿真中给波数域采样$G({k_x},{k_y},{k_z})$增加噪声,在图像域信噪比(Signal to Noise Ratio, SNR)为13 dB, 18 dB, 23 dB时分别进行三维联合重建,并记录极化散射矩阵估计结果如表5所示。从表5中变型后的结果可以看出,当SNR为18 dB和23 dB时,蓝色部分表示的各个目标的散射幅度相近,这与表2中目标真实的幅度情况一致。此外,红色部分表示的各个目标的极化散射矩阵估计结果也与表2中目标真实的极化散射矩阵基本一致。当SNR=13 dB时,可以看出各个目标的散射幅度估计结果差异略大,特别是30°二面角和45°二面角的散射幅度估计结果相差0.09。此外,与表2中目标真实的极化散射矩阵对比,可以看出对于三面角、偶极子和45°二面角,其极化散射矩阵的估计值和真实值基本一致,但对于30°二面角,受噪声影响,其极化散射矩阵的估计值和真实值略有差异。

      表 5  不同SNR下仿真目标的极化散射矩阵估计结果

      Table 5.  Polarization scattering matrix estimation results of simulated targets under different SNR

      目标类型SNR=13 dBSNR=18 dBSNR=23 dB


      变型
      $\left[\!\! {\begin{array}{*{20}{c}} {0.01 + 0.39{\rm{j}}} \!\!\!&\!\!\! {0.02} \\ {0.01 + 0.01{\rm{j}}} \!\!\!&\!\!\! { - 0.02 + 0.39{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.02 + 0.4{\rm{j}}} \!\!\!&\!\!\! {0.01 - 0.02{\rm{j}}} \\ {0.03{\rm{j}}} \!\!\!&\!\!\! {0.01 + 0.4{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} {0.01 + 0.41{\rm{j}}} \!\!\!&\!\!\! {0.02{\rm{j}}} \\ { - 0.01} \!\!\!&\!\!\! {0.01 + 0.40{\rm{j}}} \end{array}} \!\!\right]$
      变型
      ${\color{blue}{0.39}}{ {\rm{e} }^{ {\rm{j8} }{ {8.7}^ \circ } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! {0.02 - 0.06{\rm{j} } } \\ {0.03 - 0.04{\rm{j} } } \!\!\!&\!\!\! {1.03 + 0.07{\rm{j} } } \end{array} } \!\!\right]}}$${\color{blue}{0.4}}{ {\rm{e} }^{ {\rm{j92} }{\rm{.} }{ {\rm{2} }^{\rm{o} } } }}{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! { - 0.04 - 0.03{\rm{j} } } \\ {0.09 - 0.01{\rm{j} } } \!\!\!&\!\!\! {1.01 - 0.06{\rm{j} } } \end{array} } \!\!\right] }}$${\color{blue}{0.41} }{ {\rm{e} }^{ {\rm{j} }{ {88.6}^ \circ } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! {0.04} \\ {0.02 + 0.03{\rm{j} } } \!\!\!&\!\!\! {0.97 - 0.01{\rm{j} } } \end{array} } \!\!\right]} }$


      变型
      $\left[\!\! {\begin{array}{*{20}{c}} {0.02 - 0.38{\rm{j}}} \!\!\!&\!\!\! {0.01{\rm{j}}} \\ { - 0.02 + 0.02{\rm{j}}} \!\!\!&\!\!\! { - 0.03{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} {0.02 - 0.41{\rm{j}}} \!\!\!&\!\!\! {0.02 + 0.01{\rm{j}}} \\ { - 0.01 + 0.01{\rm{j}}} \!\!\!&\!\!\! {0.01} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.4{\rm{j}}} \!\!\!&\!\!\! {0.02{\rm{j}}} \\ 0 \!\!\!&\!\!\! 0 \end{array}} \!\!\right]$
      变型
      ${\color{blue}{0.3} }{\rm{8} }{ {\rm{e} }^{ {\rm{j( - 86} }{\rm{.} }{ {\rm{8} }^{\rm{o} } }{\rm{)} } } }\!{\color{red}{\left[\!\!\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! { - 0.03} \\ { - 0.05 \!-\! 0.05{\rm{j} } } \!\!\!&\!\!\! {0.08 \!-\! 0.13{\rm{j} } } \end{array} } \!\!\!\right]} }$$\;{\color{blue}{0.41} }{ {\rm{e} }^{ {\rm{j( - 87} }{\rm{.} }{ {\rm{3} }^{\rm{o} } }{\rm{)} } } }\!{\color{red}{\left[\!\!\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! { - 0.02 \!+\! 0.06{\rm{j} } } \\ { - 0.02 \!-\! 0.01{\rm{j} } } \!\!\!&\!\!\! { - 0.01 \!+\! 0.01{\rm{j} } } \end{array} } \!\!\!\!\right]} }$${\color{blue}{0.4} }{ {\rm{e} }^{ {\rm{j - 9} }{ {\rm{0} }^{\rm{o} } } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 1 \!\!\!&\!\!\! { - 0.05} \\ 0 \!\!\!&\!\!\! 0 \end{array} } \!\!\right]} }$
      30°
      二面
      变型
      $\left[\!\! {\begin{array}{*{20}{c}} { - 0.16 + 0.02{\rm{j}}} \!\!\!&\!\!\! { - 0.28 + 0.04{\rm{j}}} \\ { - 0.39 + 0.07{\rm{j}}} \!\!\!&\!\!\! {0.2 - 0.04{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.22 + 0.03{\rm{j}}} \!\!\!&\!\!\! { - 0.37 + 0.06{\rm{j}}} \\ { - 0.39 + 0.06{\rm{j}}} \!\!\!&\!\!\! {0.24 - 0.04{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.22 + 0.04{\rm{j}}} \!\!\!&\!\!\! { - 0.39 + 0.06{\rm{j}}} \\ { - 0.37 + 0.05{\rm{j}}} \!\!\!&\!\!\! {0.22 - 0.04{\rm{j}}} \end{array}} \!\!\right]$
      变型
      ${\color{blue}{0.33} }{ {\rm{e} }^{ {\rm{j172} }{\rm{.} }{ {\rm{4} }^{\rm{o} } } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } {0.5} \!\!\!&\!\!\! {0.84} \\ {1.18 - 0.07{\rm{j} } } \!\!\!&\!\!\! { - 0.6 + 0.04{\rm{j} } } \end{array} } \!\!\right]} }$${\color{blue}{0.45} }{ {\rm{e} }^{ {\rm{j172} }{\rm{.} }{ {\rm{3} }^{\rm{o} } } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } {0.5} \!\!\!&\!\!\! {0.85 - 0.03{\rm{j} } } \\ {0.89 - 0.02{\rm{j} } } \!\!\!&\!\!\! { - 0.53 + 0.03{\rm{j} } } \end{array} } \!\!\right]} }$${\color{blue}{0.45} }{ {\rm{e} }^{ {\rm{j169} }{\rm{.} }{ {\rm{4} }^{\rm{o} } } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } {0.5} \!\!\!&\!\!\! {0.87 + 0.02{\rm{j} } } \\ {0.84 + 0.04{\rm{j} } } \!\!\!&\!\!\! { - 0.5} \end{array} } \!\!\right]} }$
      45°
      二面
      变型
      $\left[\!\! {\begin{array}{*{20}{c}} {0.02} \!\!\!&\!\!\! { - 0.41 - 0.06{\rm{j}}} \\ { - 0.45 - 0.06{\rm{j}}} \!\!\!&\!\!\! { - 0.02} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} { - 0.01 + 0.01{\rm{j}}} \!\!\!&\!\!\! { - 0.46 - 0.06{\rm{j}}} \\ { - 0.44 - 0.08{\rm{j}}} \!\!\!&\!\!\! { - 0.01{\rm{j}}} \end{array}} \!\!\right]$$\left[\!\! {\begin{array}{*{20}{c}} 0 \!\!\!&\!\!\! { - 0.43 - 0.01{\rm{j}}} \\ { - 0.42 - 0.08{\rm{j}}} \!\!\!&\!\!\! {0.02{\rm{j}}} \end{array}} \!\!\right]$
      变型
      ${\color{blue}{0.42} }{ {\rm{e} }^{ {\rm{j( - 171} }{\rm{.} }{ {\rm{8} }^{\rm{o} } }{\rm{)} } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } { - 0.05} \!\!\!&\!\!\! 1 \\ {1.01} \!\!\!&\!\!\! {0.04} \end{array} } \!\!\right]} }$${\color{blue}{0.46} }{ {\rm{e} }^{ {\rm{j( - 172} }{\rm{.} }{ {\rm{4} }^{\rm{o} } }{\rm{)} } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } {0.02 \!-\! 0.01{\rm{j} } } \!\! & \!\! 1 \\ {0.96 \!+\! 0.05{\rm{j} } } \!\! & \!\! {0.01 \!+\! 0.02{\rm{j} } } \end{array} } \!\!\right]} }$${\color{blue}{0.44} }{ {\rm{e} }^{ {\rm{j( - 170} }{\rm{.} }{ {\rm{7} }^{\rm{o} } }{\rm{)} } } }{\color{red}{\left[\!\! {\begin{array}{*{20}{c} } 0 \!\!\!&\!\!\! 1 \\ {0.98 \!+\! 0.01{\rm{j} } } \!\!\!&\!\!\! { - 0.02 \!-\! 0.05{\rm{j} } } \end{array} } \!\!\right]} }$

      根据上述仿真分析可得,本文方法的性能不受目标类型影响,不论目标的稀疏支撑集是否一致,本文方法都可以较准确地估计出极化散射矩阵。相比较于各个极化通道分别独立多维重建不能保证极化散射矩阵提取的准确性,本文的联合多维重建方法对极化散射矩阵的提取准确且估计精度高。此外,本文方法具有一定的抗噪性,在SNR=13 dB时,极化散射矩阵的估计结果开始变差。

    • 使用Slicy模型的电磁计算数据进行仿真实验。该模型由三面角、二面角、圆柱等组成,图3给出了其三维CAD模型及尺寸信息。回波数据由CST微波工作室生成,采用弹跳射线法。表6记录了Slicy仿真实验中的各项参数。

      图  3  Slicy的CAD模型

      Figure 3.  CAD model of Slicy

      表 6  Slicy的仿真参数

      Table 6.  Simulation parameters of Slicy

      参数类型参数取值
      雷达扫描参数频率范围[8 GHz, 12 GHz]
      频率采样间隔20 MHz
      方位角范围[–4°, 6°]
      方位角采样间隔1/14°
      俯仰角范围[18°, 42°]
      俯仰角采样间隔1/14°
      极化方式HH, HV, VH, VV
      场景参数方位角0°沿x 轴正方向
      俯仰角0°沿z 轴正方向
      x 轴方向场景范围[–0.6 m, 0.6 m]
      y 轴方向场景范围[–0.9 m, 0.9 m]
      z 轴方向场景范围[0 m, 0.75 m]

      图4给出了Slicy模型的4个极化通道的成像结果,其中左列是HH, HV, VH和VV极化通道下的二维投影结果,右列是HH, HV, VH和VV极化通道下的三维结果。以HH极化通道成像结果的峰值为0 dB,图4中所有子图显示成像结果的幅度阈值为–30 dB,可以看出成像结果稀疏,旁瓣基本都被去除掉,验证了本文基于稀疏重构的重建方法在旁瓣抑制上的有效性。观察两个主极化通道,其成像结果相似,都可清晰地看到Slicy模型的外形轮廓,Slicy模型中的两个帽型结构,一个三面角和两个二面角都被重建出来,成像结果的位置和尺寸与图3中的模型吻合。观察两个交叉极化通道,目标的散射相对较弱,轮廓不太清晰,有些部分没有被重建出来,且两个交叉极化通道的结果略有不同,Slicy模型的一个帽型结构的强度有些差异。

      图  4  各个极化通道的Slicy成像结果

      Figure 4.  Each polarization channel’s imaging results of Slicy

      图5给出了Slicy模型的多维重建结果,其中图5(a)为全极化三维成像结果,图5(b)图为Cameron分解结果。对比图5(a)的全极化三维成像结果和图4中的各个极化通道的三维成像结果,可以看出图5(a)图4中的主极化通道的结果基本一致,不存在散射中心失配的问题,这也就验证了本文的联合重建方法的有效性,即保证了不同极化通道、不同维度的稀疏支撑集一致。图5(b)中不同的图标表示由Cameron分解得到的不同的散射类别。对比图5(b)中的结果和图3中的模型可知,各散射类型判别结果基本与实际情况相符合,验证了本文方法对极化散射矩阵提取的有效性。但由于Cameron分解自身的局限性,散射类型判别存在一定的不确定性,Slicy模型中的二面角和帽型结构没法分辨。

      图  5  Slicy全极化联合多维重建结果

      Figure 5.  Full polarization joint multi-dimensional reconstruction results of Slicy

      使用卫星模型的电磁计算数据进行仿真实验,图6给出了其三维CAD模型,卫星中心位于坐标原点,太阳能两帆板外侧的距离为8.21 m,星体最长尺寸3.55 m。回波数据由电磁计算软件[23]生成,采用积分方程快速计算方法。表7记录了卫星仿真实验中的各项参数。

      图  6  卫星的CAD模型

      Figure 6.  CAD model of satellite

      表 7  卫星的仿真参数

      Table 7.  Simulation parameters of satellite

      参数类型参数取值
      雷达扫描参数频率范围[9 GHz, 11 GHz]
      频率采样间隔20 MHz
      方位角范围[–20°, 20°]
      方位角采样间隔0.1°
      俯仰角范围[45°, 65°]
      俯仰角采样间隔0.2°
      极化方式HH, HV, VH, VV
      场景参数方位角0°沿x 轴正方向
      俯仰角0°沿z 轴正方向
      x 轴方向场景范围[–0.5 m, 0.5 m]
      y 轴方向场景范围[–4.105 m, 4.105 m]
      z 轴方向场景范围[–1.775 m, 1.775 m]

      图7给出了卫星的4个极化通道的成像结果,其中左列是HH, HV, VH和VV极化通道下的二维投影结果,右列是HH, HV, VH和VV极化通道下的三维结果。以HH极化通道成像结果的峰值为0 dB,图7中所有子图显示成像结果的幅度阈值为–30 dB,可以看出成像结果稀疏,卫星的轮廓很清晰,旁瓣基本都被去除掉了,进一步验证了本文基于稀疏重构的重建方法在旁瓣抑制上的有效性。观察两个主极化通道,都可清晰地看到卫星模型的外形轮廓,卫星主体的柱状结构和卫星两侧的太阳能帆板都被重建出来,与图6中的模型吻合。不过两个主极化通道的成像结果略有不同,HH极化通道中卫星两个帆板的散射强度明显大于VV极化通道。观察两个交叉极化通道,其成像结果相似,目标的散射相对较弱,轮廓不太清晰,卫星两侧的太阳能帆板没有被重建出来,只有卫星主体柱状结构被重建出来。

      图  7  各个极化通道的卫星成像结果

      Figure 7.  Each polarization channel’s imaging results of satellite

      图8给出了卫星的多维重建结果,其中图8(a)为全极化三维成像结果,图8(b)为Cameron分解结果。对比图8(a)的全极化三维成像结果和图7中的各个极化通道的三维成像结果,可以看出图8(a)图7中的HH极化通道的结果基本一致,因此不存在散射中心失配的问题,进一步验证了本文联合重建方法的有效性。图8(b)中不同的图标表示由Cameron分解得到的不同的散射类别,可以看出卫星的两个帆板上的散射类型为圆柱体,两个帆板下沿以及帆板与卫星主体连接处的散射类型为平板/三面角,卫星主体上部的散射类型为圆柱体,卫星主体中部比较复杂,散射类型较多,卫星主体下部的散射类型为平板/三面角和帽型/二面角的组合。由于Cameron分解自身的局限性,二面角和帽型结构,三面角和平板没法分辨。对比图8(b)中的结果和图6中的模型,各散射类型判别结果基本与实际情况相符合,进一步验证了本文方法对极化散射矩阵提取的有效性。

      图  8  卫星全极化联合多维重建结果

      Figure 8.  Full polarization joint multi-dimensional reconstruction results of satellite

    • 本文提出一种基于稀疏重构的全极化联合多维重建方法。该方法对所有极化通道和所有维度的散射系数进行了联合稀疏约束,不仅保证了不同极化通道、不同维度的稀疏支撑集一致,并且充分利用了数据之间的关联性带来的额外信息。此外,该方法解决了同时获取极化和三维信息面临的计算规模大的问题,使联合多维重建可以在不改变解的精度的前提下高效地进行。实验结果表明本文方法的性能不受目标类型影响,具有一定的抗噪性,可以实现极化散射矩阵和目标三维信息同时获取,有效地克服了散射中心失配的问题,得到了准确的目标多维重建结果。此外,该方法保持了稀疏重构的优势,有效地抑制了旁瓣,得到的重建结果分辨率高。

参考文献 (23)

目录

    /

    返回文章
    返回