留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

烟花算法在SAR图像属性散射中心参数估计中的应用

周志洪 陈秀真 马进 夏正敏

周志洪, 陈秀真, 马进, 夏正敏. 烟花算法在SAR图像属性散射中心参数估计中的应用[J]. 红外与激光工程, 2022, 51(8): 20210581. doi: 10.3788/IRLA20210581
引用本文: 周志洪, 陈秀真, 马进, 夏正敏. 烟花算法在SAR图像属性散射中心参数估计中的应用[J]. 红外与激光工程, 2022, 51(8): 20210581. doi: 10.3788/IRLA20210581
Zhou Zhihong, Chen Xiuzhen, Ma Jin, Xia Zhengmin. Application of firework algorithm into parameter estimation of attributed scattering centers in SAR images[J]. Infrared and Laser Engineering, 2022, 51(8): 20210581. doi: 10.3788/IRLA20210581
Citation: Zhou Zhihong, Chen Xiuzhen, Ma Jin, Xia Zhengmin. Application of firework algorithm into parameter estimation of attributed scattering centers in SAR images[J]. Infrared and Laser Engineering, 2022, 51(8): 20210581. doi: 10.3788/IRLA20210581

烟花算法在SAR图像属性散射中心参数估计中的应用

doi: 10.3788/IRLA20210581
基金项目: 上海市工业强基专项(GYQJ-2018-3-03)
详细信息
    作者简介:

    周志洪,男,讲师,博士,主要研究方向为车联网安全、网络安全检评等

    通讯作者: 陈秀真,女,副教授,博士,主要研究方向为大数据与人工智能应用、网络空间安全综合管理新技术等
  • 中图分类号: TP753

Application of firework algorithm into parameter estimation of attributed scattering centers in SAR images

Funds: Shanghai Industry Strong Base Special Programm(GYQJ-2018-3-03)
  • 摘要: 针对合成孔径雷达(SAR)属性散射中心估计问题,提出基于烟花算法的方法。首先,在图像域对SAR图像中高能量区域进行分割解耦,获得单个独立散射中心在图像域的表现形式。在此基础上,以属性散射中心参数化模型为基础,构建优化问题,对分离出来的单个散射中心进行最优参数的搜索。在此阶段,引入烟花算法进行参数寻优。该算法具有强大的全局和局部搜索能力,在保证优化精度的条件下避免陷入局部最优,从而保证散射中心参数估计的可靠性。在原始图像中剔除求解后的单个散射中心,对残余图像进行高能量区域分割,序惯估计下一个散射中心的属性参数。最终,获取输入SAR图像上所有散射中心的参数集。实验中,首先基于MSTAR数据集中的SAR图像进行参数估计验证,通过参数估计结果与原始图像的对比以及基于估计参数集对原始图像进行重构,反映了提出算法的有效性。此外,实验还基于估计得到的属性参数进行SAR目标识别算法验证,通过与其他参数估计算法在相同条件下进行识别性能的对比,进一步体现了提出方法在属性散射中心参数估计上的性能优势。
  • 图  1  MSTAR目标的光学及SAR图像

    Figure  1.  Optical and SAR images of MSTAR targets

    图  2  估计得到的散射中心位置示意

    Figure  2.  Illustration of locations of estimated scattering centers

    图  3  基于估计散射中心的重构目标及残差

    Figure  3.  Reconstructed target and residuals based on the estimated scattering centers

    表  1  一幅MSTAR SAR图像的属性散射中心参数估计结果

    Table  1.   Parameter estimation results of attributed scattering centers in a MSTAR SAR image

    No.$ \left| A \right| $${x_p}/{\text{m} }$${y_p}/{\text{m} }$${L_p}/{\text{m} }$$ \alpha $$ {\gamma _p} $$ {\phi _p} $
    1 513.62 −1.97 7.53 0.00 0.00 0.00 0.00
    2 483.50 −8.08 12.28 2.57 0.00 0.00 0.17
    3 216.39 3.18 3.75 0.00 0.00 1.82 0.00
    4 131.06 10.60 −1.87 0.00 0.00 3.70 0.00
    5 164.36 −5.68 6.08 0.00 0.00 0.00 0.00
    6 138.23 −14.85 3.40 0.00 0.00 −0.18 0.00
    7 133.72 6.99 1.63 0.00 1.00 0.00 0.00
    8 523.03 −4.67 2.71 5.94 1.00 0.00 0.29
    9 120.12 −13.06 7.07 0.00 0.00 0.00 0.00
    10 148.23 6.69 −5.80 2.50 0.00 0.00 −0.02
    11 95.01 −10.32 5.00 0.00 0.00 0.00 0.00
    12 88.15 1.26 1.89 0.00 1.00 −0.19 0.00
    13 86.96 0.99 7.27 0.00 1.00 −0.19 0.00
    14 83.63 −2.94 10.30 0.00 1.00 0.00 0.00
    15 80.11 −8.16 6.75 0.00 0.00 −0.19 0.00
    下载: 导出CSV

    表  2  10类目标的训练和测试样本

    Table  2.   Training and test samples of the ten targets

    ClassBMP2BTR70T72T62BDRM2BTR60ZSU23/4D7ZIL1312S1
    Training (17°)232233231299298256299299299299
    Testing (15°)195196196273274195274274274274
    下载: 导出CSV

    表  3  散射中心匹配识别方法在不同提取算法下的平均识别率

    Table  3.   Average recognition rates of the scattering center matching method using different parameter estimation algorithms

    Parameter estimation algorithmAverage recognition rate
    Simulated annealing98.16
    Dynamic particle swarm optimization98.56
    Sparse representation98.74
    Proposed99.02
    下载: 导出CSV
  • [1] 郭琨毅, 殷红成, 盛新庆. 雷达目标的散射中心建模研究[J]. 电波科学学报, 2020, 35(1): 106–115

    Guo K Y, Yin H C, Sheng X Q. Research on scattering center modeling for radar target [J]. Chinese Journal of Radio Science, 2020, 35(1): 106-115. (in Chinese)
    [2] 文贡坚, 朱国强, 殷红成等. 基于三维电磁散射参数化模型的SAR目标识别方法[J]. 雷达学报, 2017, 6(2): 115–135

    Wen G J, Zhu G Q, Yin H C, et al. SAR ATR based on 3D parametric electromagnetic scattering model [J]. Journal of Radar, 2017, 6(2): 115-135. (in Chinese)
    [3] 丁柏圆, 钟金荣, 马聪慧, 等. 基于属性散射中心模型的SAR目标重构可视化增强方法[J]. 雷达学报, 2013, 2(4): 499–506

    Ding B Y, Zhong J R, Ma C H, et al. SAR target reconstruction visualization enhancement based on attributed scattering center model [J]. Journal of Radar, 2013, 2(4): 499-506. (in Chinese)
    [4] 刘阳. 基于属性散射中心多层次匹配的SAR目标识别方法[J]. 电子测量与仪器学报, 2019, 33(11): 192-198

    Liu Y. Target recognition of SAR images based on multi-level matching of attributed scattering centers [J]. Journal of Electronic Measurement and Instrumentation, 2019, 33(11): 192-198. (in Chinese)
    [5] 丁柏圆, 文贡坚, 余连生等. 属性散射中心匹配及其在SAR目标识别中的应用[J]. 雷达学报, 2017, 6(2): 157–166

    Ding B Y, Wen G J, Yu L S, et al. Matching of attributed scattering center and its application to synthetic aperture radar automatic target recognition [J]. Journal of Radar, 2017, 6(2): 157-166. (in Chinese)
    [6] 张婷, 蔡德饶. 基于属性散射中心匹配的噪声稳健SAR目标识别方法[J]. 中国电子科学研究院学报, 2019, 14(6): 557-563 doi:  10.3969/j.issn.1673-5692.2019.06.001

    Zhang T, Cao D R. SAR target recognition method based on attributed scattering center matching with robustness to noise corruption [J]. Journal of China Academy of Electronics and Information Technology, 2019, 14(6): 557-563. (in Chinese) doi:  10.3969/j.issn.1673-5692.2019.06.001
    [7] Ding B Y, Wen G J, Zhong J R, et al. A robust similarity measure for attributed scattering center sets with application to SAR ATR [J]. Neurocomputing, 2017, 219: 130-143. doi:  10.1016/j.neucom.2016.09.007
    [8] 张华, 张素丽, 何树吉. 基于幅相分离和动态粒子群算法的SAR图像属性散射中心参数估计[J]. 中国电子科学研究院学报, 2019, 14(9): 993-1000 doi:  10.3969/j.issn.1673-5692.2019.09.017

    Zhang H, Zhang S L, He S J. Parameter estimation of attributed scattering centers in SAR images based on amplitude-phase separation and DPSO [J]. Journal of China Academy of Electronics and Information Technology, 2019, 14(9): 993-1000. (in Chinese) doi:  10.3969/j.issn.1673-5692.2019.09.017
    [9] 蒋文, 李王哲. 基于幅相分离的属性散射中心参数估计新方法[J]. 雷达学报, 2019, 8(5): 606-615 doi:  10.12000/JR18097

    Jiang W, Li W Z. A new method for parameter estimation of attributed scattering centers based on amplitude-phase separation [J]. Journal of Radar, 2019, 8(5): 606-615. (in Chinese) doi:  10.12000/JR18097
    [10] 段佳, 张磊, 盛佳恋, 等. 独立属性散射中心参数降耦合估计方法[J]. 电子与信息学报, 2012, 34(8): 1853-1859

    Duan J, Zhang L, Sheng J L, et al. Parameter decouple and estimation of independent attributed scattering centers [J]. Journal of Electronics & Information Technology, 2012, 34(8): 1853-1859. (in Chinese)
    [11] 陈丛. 基于狼群算法的SAR 图像属性散射中心参数估计[J]. 红外与激光工程, 2021, 50(4): 20200409 doi:  10.3788/IRLA20200409

    Chen C. Parameter estimation of attributed scattering centers in SAR images based on wolf pack algorithm [J]. Infrared and Laser Engineering, 2021, 50(4): 20200409. (in Chinese) doi:  10.3788/IRLA20200409
    [12] 占荣辉, 胡杰民, 张军. 基于压缩感知的二维GTD模型参数估计方法[J]. 电子与信息学报, 2013, 35(2): 419-425

    Zhan R H, Hu J M, Zhang J. A novel method for parametric estimation of 2D geometrical theory of diffraction model based on compressed sensing [J]. Journal of Electronics & Information Technology, 2013, 35(2): 419-425. (in Chinese)
    [13] 李飞, 纠博, 刘宏伟, 等. 基于稀疏表示的SAR图像属性散射中心参数估计算法[J]. 电子与信息学报, 2014, 36(4): 931-937

    Li F, Jiu B, Liu H W, et al. Sparse representation based algorithm for estimation of attributed scattering center parameter on SAR imagery [J]. Journal of Electronics & Information Technology, 2014, 36(4): 931-937. (in Chinese)
    [14] 张肖, 闫恒庄, 周建江, 等. 基于组稀疏表示的二维全极化散射中心提取[J]. 电光与控制, 2016, 23(2): 26-30

    Zhnag X, Yan H Z, Zhou J J, et al. Extraction of 2D full polarization scattering centers based on group sparse representation [J]. Electronics Optics & Control, 2016, 23(2): 26-30. (in Chinese)
    [15] Tan Y, Zhu Y. Fireworks algorithm for optimization[C]// International Conference in Swarm Intelligence. Berlin: Springer, 2010: 355-364.
    [16] 范虹, 候存存, 朱艳春, 等. 烟花算法优化的软子空间MR图像聚类算法[J]. 软件学报, 2017, 28(11): 3080-3093

    Fan H, Hou C C, Zhu Y C, et al. Soft subspace algorithm for MR image clustering based on fireworks optimization algorithm [J]. Journal of Software, 2017, 28(11): 3080-3093. (in Chinese)
    [17] 薛俊杰, 王瑛, 孟祥飞, 等. 二进制反向学习烟花算法求解多维背包问题[J]. 系统工程与电子技术, 2017, 39(2): 451-458 doi:  10.3969/j.issn.1001-506X.2017.02.33

    Xue J J, Wang Y, Meng X F, et al. Binary opposite backward learning fireworks algorithm for multidimensional knapsack problem [J]. Systems Engineering and Electronics, 2017, 39(2): 451-458. (in Chinese) doi:  10.3969/j.issn.1001-506X.2017.02.33
    [18] 包晓晓, 叶春明, 黄霞. 烟花算法求解JSP问题的研究[J]. 计算机工程与应用, 2017, 53(3): 247-251 doi:  10.3778/j.issn.1002-8331.1505-0008

    Bao X X, Ye C M, Huang X. Solving job-shop scheduling problem by firework algorithm [J]. Computer Engineering and Applications, 2017, 53(3): 247-251. (in Chinese) doi:  10.3778/j.issn.1002-8331.1505-0008
  • [1] 李正伟, 黄孝斌, 胡尧.  基于二维随机投影特征典型相关分析融合的SAR ATR方法 . 红外与激光工程, 2022, 51(10): 20220029-1-20220029-8. doi: 10.3788/IRLA20220029
    [2] 陆建华.  融合CNN和SRC决策的SAR图像目标识别方法 . 红外与激光工程, 2022, 51(3): 20210421-1-20210421-7. doi: 10.3788/IRLA20210421
    [3] 游丽.  基于块稀疏贝叶斯学习的SAR图像目标方位角估计方法 . 红外与激光工程, 2022, 51(4): 20210282-1-20210282-6. doi: 10.3788/IRLA20210282
    [4] 吴剑波, 陆正武, 关玉蓉, 王庆东, 姜国松.  二维压缩感知多投影矩阵特征融合的SAR目标识别方法 . 红外与激光工程, 2021, 50(6): 20200531-1-20200531-7. doi: 10.3788/IRLA20200531
    [5] 李宁, 王军敏, 司文杰, 耿则勋.  基于最大熵准则的多视角SAR目标分类方法 . 红外与激光工程, 2021, 50(12): 20210233-1-20210233-7. doi: 10.3788/IRLA20210233
    [6] 李亚娟.  结合多决策准则稀疏表示的SAR图像目标识别方法 . 红外与激光工程, 2021, 50(8): 20210138-1-20210138-8. doi: 10.3788/IRLA20210138
    [7] 尚珊珊, 余子开, 范涛, 金利民.  高斯过程模型在SAR图像目标识别中的应用 . 红外与激光工程, 2021, 50(7): 20200337-1-20200337-7. doi: 10.3788/IRLA20200337
    [8] 伍友龙.  多元经验模态分解及在SAR图像目标识别中的应用 . 红外与激光工程, 2021, 50(4): 20200236-1-20200236-7. doi: 10.3788/IRLA20200236
    [9] 马丹丹.  图像分块匹配的SAR目标识别方法 . 红外与激光工程, 2021, 50(10): 20210120-1-20210120-8. doi: 10.3788/IRLA20210120
    [10] 刘志超, 屈百达.  复数二维经验模态分解在SAR目标识别中的应用 . 红外与激光工程, 2021, 50(5): 20200309-1-20200309-8. doi: 10.3788/IRLA20200309
    [11] 陈丛.  基于狼群算法的SAR图像属性散射中心参数估计 . 红外与激光工程, 2021, 50(4): 20200409-1-20200409-7. doi: 10.3788/IRLA20200409
    [12] 王保平, 马健钧, 张研, 方阳.  基于最小能量准则的参数估计圆周SAR三维成像方法 . 红外与激光工程, 2019, 48(7): 726001-0726001(7). doi: 10.3788/IRLA201948.0726001
    [13] 郭力仁, 胡以华, 王云鹏.  相干激光探测中微动参数估计的克拉美-罗界 . 红外与激光工程, 2017, 46(7): 706002-0706002(7). doi: 10.3788/IRLA201746.0706002
    [14] 杨新锋, 张金娜, 滕书华, 崔宇.  利用压缩感知实现随机变频雷达散射中心估计 . 红外与激光工程, 2016, 45(5): 526004-0526004(6). doi: 10.3788/IRLA201645.0526004
    [15] 蒋立辉, 章典, 陈星, 熊兴隆.  基于谱估计的Fabry-Perot标准具透过率曲线参数估计方法 . 红外与激光工程, 2016, 45(9): 906006-0906006(8). doi: 10.3788/IRLA201645.0906006
    [16] 刘小阳, 孙广通, 宋萍, 刘军, 李峰.  利用MODIS水汽数据进行ASAR干涉测量大气改正研究 . 红外与激光工程, 2014, 43(12): 4111-4116.
    [17] 张闻文, 刘婧婧, 陈钱, 顾国华.  电子倍增CCD噪声参数的估计方法 . 红外与激光工程, 2014, 43(12): 4146-4152.
    [18] 陈骁, 石志广, 杨卫平, 凌建国.  利用三维点云的圆锥状弹头目标参数估计 . 红外与激光工程, 2014, 43(10): 3480-3485.
    [19] 邹盼, 刘晖, 张闻文, 陈钱, 顾国华, 张连东.  利用期望最大化算法的EMCCD噪声分布模型的参数估计 . 红外与激光工程, 2013, 42(1): 268-272.
    [20] 何莲, 蔡敬菊, 张启衡.  改进的基于弦切变换的目标检测方法 . 红外与激光工程, 2012, 41(1): 239-247.
  • 加载中
图(3) / 表(3)
计量
  • 文章访问数:  187
  • HTML全文浏览量:  39
  • PDF下载量:  31
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-12-02
  • 修回日期:  2022-02-20
  • 刊出日期:  2022-08-31

烟花算法在SAR图像属性散射中心参数估计中的应用

doi: 10.3788/IRLA20210581
    作者简介:

    周志洪,男,讲师,博士,主要研究方向为车联网安全、网络安全检评等

    通讯作者: 陈秀真,女,副教授,博士,主要研究方向为大数据与人工智能应用、网络空间安全综合管理新技术等
基金项目:  上海市工业强基专项(GYQJ-2018-3-03)
  • 中图分类号: TP753

摘要: 针对合成孔径雷达(SAR)属性散射中心估计问题,提出基于烟花算法的方法。首先,在图像域对SAR图像中高能量区域进行分割解耦,获得单个独立散射中心在图像域的表现形式。在此基础上,以属性散射中心参数化模型为基础,构建优化问题,对分离出来的单个散射中心进行最优参数的搜索。在此阶段,引入烟花算法进行参数寻优。该算法具有强大的全局和局部搜索能力,在保证优化精度的条件下避免陷入局部最优,从而保证散射中心参数估计的可靠性。在原始图像中剔除求解后的单个散射中心,对残余图像进行高能量区域分割,序惯估计下一个散射中心的属性参数。最终,获取输入SAR图像上所有散射中心的参数集。实验中,首先基于MSTAR数据集中的SAR图像进行参数估计验证,通过参数估计结果与原始图像的对比以及基于估计参数集对原始图像进行重构,反映了提出算法的有效性。此外,实验还基于估计得到的属性参数进行SAR目标识别算法验证,通过与其他参数估计算法在相同条件下进行识别性能的对比,进一步体现了提出方法在属性散射中心参数估计上的性能优势。

English Abstract

    • 合成孔径雷达(Synthetic aperture radar, SAR)具备全天时、全天候的高分辨率成像能力,在军民领域都有着广泛运用。与清晰直观的光学图像不同,SAR图像反映目标的电磁散射特性,其中目标的外观特点不明显,给图像解译带来了障碍。为定量化描述SAR图像的分布特性,研究人员通过散射中心模型的方式对局部散射现象进行建模,最为代表性的是属性散射中心模型[1-3]。属性散射中心模型包含了与目标局部结构相关联的属性参数,如相对强度、位置、长度等。基于属性散射中心的SAR目标识别方法也充分验证了属性参数对于区分不同目标的有效性[4-7]。然而,属性散射中心模型形式复杂,参数估计十分复杂。早期,研究人员通过图像域分治策略,实现逐个散射中心的序惯估计。在对单个散射中心进行参数估计时,通常采用较为经典的优化算法,如模拟退火、粒子群算法、狼群算法等[8-11]。根据压缩感知理论,参考文献[12-14]设计了基于稀疏表示的属性散射中心估计算法,认为单幅SAR图像的属性散射中心在整个参数空间属于稀疏分布。稀疏表示方法避免了图像域的局部解耦过程,但其最终结果与图像域的一致性往往难以保证,具有较强的随机性。为此,从确保散射中心参数物理意义的角度,图像域的散射中心估计算法具有更强的可靠性。

      文中在传统图像域属性散射中心估计方法的基础上,引入烟花算法,实现对单个散射中心的参数优化。对于图像域解耦的单个散射中心图像区域,通过烟花算法获得最佳的属性参数。烟花算法启发于烟花燃放过程,由Tan和Zhu在2010年设计提出[15]。其中,每一个烟花被视为解空间的一个可行解,烟花爆炸产生一定数量火花的过程则反映邻域搜索的过程。每个烟花爆炸半径和产生的火花数不同,适应度较差的烟花爆炸半径更大,这让优化算法具有全局搜索能力;适应度较差的烟花爆炸半径更小,这让优化算法具有局部搜索能力。不同烟花按照各自适应度进行动态分配和交互使得优化算法在全局和局部搜索能力之间达到平衡。此外,烟花算法加入了变异火花增加了种群的多样性,避免算法陷入局部最优。目前,烟花算法已经在任务分配、图像处理等领域得到了广泛应用和验证,具有更高的精度和稳健性[16-18]。为此,文中将烟花算法引入到单个属性散射中心的参数估计中,获得更为可靠的估计结果。实验中,采用MSTAR数据集对提出方法进行验证,分别测试了方法对于SAR图像的估计结果以及对于基于属性散射中心匹配的识别方法的贡献。实验结果证明了提出方法的有效性和相比现有几类方法的优越性。

    • 作为描述高频区域雷达目标局部散射特性的参数化模型,属性散射中心模型相比先前的GTD模型具有更强的物理层和信号层描述能力。参考文献[4-7]基于属性散射中心设计SAR目标方法,验证了其在提高目标识别稳健性方面的显著优势。基于属性散射中心模型描述单个散射中心的散射特性如下:

      $$ \begin{split} {E_i}(f,\phi ;{\theta _i}) =& {A_i} \cdot {\Bigg(j\frac{f}{{{f_c}}}\Bigg)^{{\alpha _i}}} \cdot \exp \Bigg(\frac{{ - j4\pi f}}{c}({x_i}\cos \phi + {y_i}\sin \phi )\Bigg)\cdot \\& {\text{ }} {\text{sinc}}\Bigg(\frac{{2\pi f}}{c}{L_i}{\text{sin(}}\phi - \overline {{\phi _i}} {)\Bigg)} \cdot {{{\rm{exp}}( - 2}}\pi f{\gamma _i}{\text{sin}}\phi {\text{)}} \end{split} $$ (1)

      式中:$ f $代表入射波频率;$ \phi $为方位角;$ {\theta _i} = [{A_i},{\alpha _i},{x_i},{y_i},{L_i}, \overline {{\phi _i}} ,{\gamma _i}] $为属性参数集,$ {A_i} $代表不同散射中心的相对幅度,$ {x_i},{y_i} $记录散射中心位置;$ {\alpha _i} $为频率依赖因子,其取值与局部结构直接相关;对于分布式散射中心,$ {L_i} $$ \overline {{\phi _i}} $分别表示其长度和方向角;对于局部式散射中心,$ {\gamma _i} $指代方位依赖性参数。

      对于一幅包含目标的SAR图像,其中存在若干个散射中心。那么基于属性散射中心模型可以对目标整体的散射特性描述如下:

      $$ E(f,\phi ;\theta ) = \sum\limits_{i = 1}^P {{E_i}(f,\phi ;{\theta _i})} $$ (2)

      式中:$ \theta {\text{ = }}\left\{ {{\theta _i}} \right\}(i = 1,2, \cdots ,P) $为目标的属性散射中心集;P代表散射中心数目。

    • 尽管属性散射中心模型的描述精度更高,但随之而来的参数估计问题也愈加复杂。总体而言,待分析的SAR图像数据可以表示为:

      $$ D(f,\phi ){\text{ = }}E(f,\phi ;\theta ){\text{ + }}N(f,\phi ) $$ (3)

      式中:$ E(f,\phi ;\theta ) $为公式(2)描述的模型分量;$ N(f,\phi ) $图像中的噪声分量,一般建模为高斯分布。在极大似然估计的框架下,目标属性散射中心集合的联合估计过程描述如下:

      $$ {\hat \theta _{ML}} = \arg \mathop {\min }\limits_\theta ||D - E(\theta )|{|^2} $$ (4)

      公式(4)给出属性散射中心参数估计的理想形式。然而,由于模型自身的复杂度以及散射中心参数数目较多,这种联合估计的优化过程非凸,很容易陷入局部最优,导致估计结果不稳健、不准确。为此,研究人员提出了图像域分治策略,即根据有效散射中心在图像域的显著特点实现单个散射中心的图像区域分离[8-11]。在此基础上,可采用优化算法实现单个散射中心的参数估计,进而序惯获得所以散射中心的参数。根据现有文献报道,图像域解耦的属性散射中心参数估计算法流程一般可区分为图像域解耦、单个散射中心参数估计以及序惯估计三个步骤。首先,对输入SAR图像进行分析,采用图像分割算法获取其中能量最大的区域,作为单个散射中心的图像域表现形式。然后,对于解耦后的单个散射中心数据,在属性散射中心模型的约束下进行优化求解,获得最佳属性参数。此处,需要引入适当的优化算法,确保参数估计的正确性。最后,在原始图像中剔除已估计成分(通过原始图像减去散射中心重构结果),更新输入图像,获得下一个散射中心的参数。通过序惯迭代更新,在满足终止条件的情况下输出目标所有散射中心的属性参数集。

    • 烟花算法受启发于烟花爆炸过程,设计过程中有效地平衡了算法的局部搜索和全局搜索,具有较强的稳健性[15-18]。该算法根据实际求解问题,首先在解空间设置N个烟花,然后计算每个烟花的适应度值以及烟花的爆炸半径和火花数。对于烟花$ {x_i} $,采用公式(5)和公式(6)分别计算器爆炸半径$ {A_i} $与火花数$ {S_i} $

      $$ {A_i} = \hat A \times \frac{{f({x_i}) - {y_{\min }} + \varepsilon }}{{\displaystyle\sum\nolimits_{i = 1}^N {\left( {f({x_i}) - {y_{\min }}} \right) + \varepsilon } }} $$ (5)
      $$ {S_i} = M \times \frac{{{y_{\max }} - f({x_i}) + \varepsilon }}{{\displaystyle\sum\nolimits_{i = 1}^N {\left( {{y_{\max }} - f({x_i})} \right) + \varepsilon } }} $$ (6)

      式中:$ \hat A $$ M $均为常数,用于限制最大爆炸幅度和火花总数。

      对于适应值较小的烟花,它能在较小范围内产生较多火花,从而获得强大的局部搜索能力。反之,对于适应值较大的烟花,它能在较大返回产生较少火花,从而获得强大的全局搜索能力。同时,为增强种群多样性,烟花算法通过随机选定烟花进行变异操作(如,高斯变异)。此外,对于超过可行域的火花,则设置一定的映射关系映射到可行域。为保持种群中的优秀信息,在爆炸和变异操作后,在初始烟花及产生火花中选择N个个体构成下一代烟花,直至迭代结束。烟花算法的具体实施按照以下步骤:

      Step1:随机产生N个初始烟花;

      Step2:计算每个烟花在爆炸操作中的爆炸半径和火花个数;

      Step3:根据爆炸半径和火花个数生成爆炸火花和变异火花;

      Step4:在烟花、爆炸火花及变异火中选取当前最优值;

      Step5:若满足迭代终止条件则输出最优值并停止迭代,否则进入Step6;

      Step6:按照既定规则选取N个新烟花,并进入Step2。

      鉴于烟花算法所具有的特点和优势,文中将其引入到属性散射中心参数估计中。对于图像域解耦得到的单个散射中心数据,采用烟花算法求解其最佳属性参数。按照序惯估计的策略,获得SAR图像中所有属性散射中心的参数集。

    • 文中基于MSTAR数据集中的SAR图像样本对所提方法进行测试。该数据集中包含了图1所示的10类车辆目标SAR图像。这些图像的分辨率达到0.3 m,直观可见目标的局部散射现象明显。文中实验中,采用所提方法对MSTAR SAR图像进行散射中心参数估计。由于这些SAR图像的散射中心参数真值未知,只能通过对比估计参数与目标局部强散射的对应性(主要是位置参数)以及重构图像与原始图像的一致性对估计结果的有效性进行初步评估。除此之外,由于属性散射中心已在SAR目标识别方法中得到广泛运用,文中通过对比文中参数估计结果与其他现有散射中心参数估计对识别结果的贡献进一步评估所提方法的有效性。

      图  1  MSTAR目标的光学及SAR图像

      Figure 1.  Optical and SAR images of MSTAR targets

    • 采用所提方法对MSTAR数据集中的一幅SAR图像进行典例验证,选取15作为散射中心数目的上限,获得具体的参数估计结果如表1所示。图2所示为位置参数在被估计SAR图像上的投影结果,可见与目标局部强散射中心吻合性较强,说明估计结果中位置参数的精度具有较高的精度。图3(a)所示为基于估计散射中心重构的目标图像,采用经典的图像相关作为测量度,计算其与被估计SAR图像的相似度达到0.92,表明重构结果与原始图像具有很强的相关性,说明散射中心各类属性参数估计结果的有效性。图3(b)为被估计图像与重构图像对比(相减)之后的残差分量。直观可见,重构残差主要是被估计图像中的背景成分,由于杂波、噪声等因素噪声。通过对被估计SAR图像中目标散射中心的估计,并据此重构,有效剔除了原始图像中存在的背景因素,对于针对性分析目标特性具有有力支撑。上述结果均表明所提方法对于SAR图像属性散射中心参数估计的有效性。

      表 1  一幅MSTAR SAR图像的属性散射中心参数估计结果

      Table 1.  Parameter estimation results of attributed scattering centers in a MSTAR SAR image

      No.$ \left| A \right| $${x_p}/{\text{m} }$${y_p}/{\text{m} }$${L_p}/{\text{m} }$$ \alpha $$ {\gamma _p} $$ {\phi _p} $
      1 513.62 −1.97 7.53 0.00 0.00 0.00 0.00
      2 483.50 −8.08 12.28 2.57 0.00 0.00 0.17
      3 216.39 3.18 3.75 0.00 0.00 1.82 0.00
      4 131.06 10.60 −1.87 0.00 0.00 3.70 0.00
      5 164.36 −5.68 6.08 0.00 0.00 0.00 0.00
      6 138.23 −14.85 3.40 0.00 0.00 −0.18 0.00
      7 133.72 6.99 1.63 0.00 1.00 0.00 0.00
      8 523.03 −4.67 2.71 5.94 1.00 0.00 0.29
      9 120.12 −13.06 7.07 0.00 0.00 0.00 0.00
      10 148.23 6.69 −5.80 2.50 0.00 0.00 −0.02
      11 95.01 −10.32 5.00 0.00 0.00 0.00 0.00
      12 88.15 1.26 1.89 0.00 1.00 −0.19 0.00
      13 86.96 0.99 7.27 0.00 1.00 −0.19 0.00
      14 83.63 −2.94 10.30 0.00 1.00 0.00 0.00
      15 80.11 −8.16 6.75 0.00 0.00 −0.19 0.00

      图  2  估计得到的散射中心位置示意

      Figure 2.  Illustration of locations of estimated scattering centers

      图  3  基于估计散射中心的重构目标及残差

      Figure 3.  Reconstructed target and residuals based on the estimated scattering centers

    • 参考文献[4-7]均是基于属性散射中心设计SAR目标识别方法,其中一种典型场景就是在标准操作条件下对图1所示的10类目标进行分类,具体设置的训练和测试集参见表2。为进一步验证提出的散射中心参数估计方法的有效性,可在同一个基于属性散射中心的SAR目标识别方法上采用不同估计算法的结果,进而对比识别性能。按照这一思路,文中采用现有文献中多种不同属性散射中心估计算法进行属性散射中心提取,进而采用参考文献[5]中的散射中心匹配方法对表2中的10类目标的测试样本进行分类。其中,选取的散射中心估计算法包括参考文献[8]中基于模拟退火算法、参考文献[9]中基于动态粒子群算法、参考文献[13]中的稀疏表示估计算法。在不同算法下,散射中心匹配识别方法取得的平均识别率如表3所示。对比可见,基于所提方法估计得到的属性散射中心可以取得最高的平均识别率,表明其对于目标识别的贡献更大,反映了其估计精度相对更高。因此,通过不同散射中心提取算法对于目标识别的贡献可以间接反映了所提方法具有更强的有效性。

      表 2  10类目标的训练和测试样本

      Table 2.  Training and test samples of the ten targets

      ClassBMP2BTR70T72T62BDRM2BTR60ZSU23/4D7ZIL1312S1
      Training (17°)232233231299298256299299299299
      Testing (15°)195196196273274195274274274274

      表 3  散射中心匹配识别方法在不同提取算法下的平均识别率

      Table 3.  Average recognition rates of the scattering center matching method using different parameter estimation algorithms

      Parameter estimation algorithmAverage recognition rate
      Simulated annealing98.16
      Dynamic particle swarm optimization98.56
      Sparse representation98.74
      Proposed99.02
    • 文中针对SAR图像目标属性散射中心估计问题,提出基于烟花算法的参数优化方法。在传统图像域分治的基础上,在对每一个属性散射中心进行参数估计是引入烟花算法。烟花算法通过模拟烟花燃放过程在解空间进行全局和局部搜索,可以获得稳健性更强、精度更高的估计结果。实验中,基于MSTAR数据集对所提方法进行性能测试。首先,基于MSTAR SAR图像开展参数估计的典例演示,通过比较估计位置参数和图像强散射点的位置以及重构图像与原始图像的相似度验证了所提方法的有效性。然后,通过比较不同参数估计算法对于SAR目标识别的贡献进一步反映了所提方法的优越性能。

参考文献 (18)

目录

    /

    返回文章
    返回