结合旋转域极化特征的极化SAR地物分类

陶臣嵩 陈思伟 李永祯 肖顺平

陶臣嵩, 陈思伟, 李永祯, 肖顺平. 结合旋转域极化特征的极化SAR地物分类[J]. 雷达学报, 2017, 6(5): 524-532. doi: 10.12000/JR16131
引用本文: 陶臣嵩, 陈思伟, 李永祯, 肖顺平. 结合旋转域极化特征的极化SAR地物分类[J]. 雷达学报, 2017, 6(5): 524-532. doi: 10.12000/JR16131
Tao Chensong, Chen Siwei, Li Yongzhen, Xiao Shunping. Polarimetric SAR Terrain Classification Using Polarimetric Features Derived from Rotation Domain[J]. Journal of Radars, 2017, 6(5): 524-532. doi: 10.12000/JR16131
Citation: Tao Chensong, Chen Siwei, Li Yongzhen, Xiao Shunping. Polarimetric SAR Terrain Classification Using Polarimetric Features Derived from Rotation Domain[J]. Journal of Radars, 2017, 6(5): 524-532. doi: 10.12000/JR16131

结合旋转域极化特征的极化SAR地物分类

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

    陶臣嵩(1993–),男,广西人,国防科学技术大学电子科学与工程学院在读硕士研究生,主要研究方向为成像雷达极化信息的解译与处理、极化特征的提取与优选等。E-mail: taochensongnudt@163.com

    陈思伟(1984–),男,四川人,博士,国防科学技术大学电子科学与工程学院讲师,主要研究方向包括雷达极化信息处理、成像雷达信息处理、目标散射建模与解译、微波遥感大数据处理与应用、微波遥感环境与灾害应用研究等。E-mail: chenswnudt@163.com

    李永祯(1977–),男,内蒙古人,博士后,国防科学技术大学电子科学与工程学院研究员,电子信息系统复杂电磁环境效应国家重点实验室副主任,主要研究方向为新体制雷达与电子对抗。E-mail: e0061@sina.com

    肖顺平(1964–),男,江西人,博士,国防科学技术大学电子科学与工程学院教授,博士生导师,电子信息系统复杂电磁环境效应国家重点实验室主任,国家高技术863专家,总装备部仿真专业组专家,中国电子学会高级会员,主要研究方向包括雷达极化信息处理、电子信息系统仿真评估技术、雷达目标识别等

    通讯作者:

    陈思伟   chenswnudt@163.com

  • 中图分类号: TN957.52

Polarimetric SAR Terrain Classification Using Polarimetric Features Derived from Rotation Domain

Funds: The National Natural Science Foundation of China (41301490, 61490692)
  • 摘要: 地物分类是极化合成孔径雷达(SAR)图像理解与解译的重要应用方向。利用H/A/α/SPAN等旋转不变特征参数的极化SAR地物分类是一种常用的分类方法。然而,目标的后向散射响应与其方位取向等姿态密切相关,极易引起散射机理的解译模糊,也限制了仅使用旋转不变特征参数作为分类特征集的极化SAR地物分类方法的精度。针对这一问题,有文献提出了在绕雷达视线的旋转域中解译目标散射特性的统一的极化矩阵旋转理论,并导出了一系列旋转域极化特征,刻画目标旋转域隐含信息。基于该理论,该文将旋转域极化特征用于极化SAR地物辨识与分类,并发展了一种结合旋转域极化特征与旋转不变特征H/A/α/SPAN的极化SAR地物分类方法。该方法将优选的旋转域极化特征参数和H/A/α/SPAN作为支持向量机(SVM)分类器的输入,利用两类特征对不同地物类别区分辨识能力的互补,以达到更好的分类性能。对AIRSAR和UAVSAR实测数据的对比实验表明,相较于仅以H/A/α/SPAN作为SVM分类器输入的传统方法,该方法得到更好的分类精度和稳健性。其中,对于AIRSAR数据十五类地物的分类,该方法总体分类精度达到92.3%,优于传统方法的91.1%。此外,对于多时相UAVSAR数据七类地物的分类,该方法平均总体分类精度达到95.72%,显著优于传统方法的87.80%,验证了该方法对多时相数据的稳健性。该文研究进一步证实了通过深入挖掘旋转域中目标极化散射信息能够为极化SAR图像的解译与应用提供新的可行途径。
  • 图  1  本文方法具体流程图

    Figure  1.  Flowchart of proposed method

    图  2  AIRSAR数据

    Figure  2.  AIRSAR data

    图  3  AIRSAR数据的分类结果

    Figure  3.  Classification results of AIRSAR data

    图  4  多时相UAVSAR数据滤波后Pauli RGB图

    Figure  4.  Filtered Pauli RGB images of multi-temporal UAVSAR data

    图  5  所用区域的真值图

    Figure  5.  Gound truth of the multi-temporal data

    图  6  传统方法对多时相UAVSAR数据分类结果

    Figure  6.  Classification results of multi-temporal UAVSAR data using conventional method

    图  7  本文方法对多时相UAVSAR数据分类结果

    Figure  7.  Classification results of multi-temporal UAVSAR data using proposed method

    表  1  旋转域极化特征参数[12]

    Table  1.   Polarimetric feature parameters derived from rotation domain[12]

    散射矩阵元素项 $A = \sqrt \bullet $ B $\omega $ ${\theta _0} = \frac{1}{\omega }{\rm{Angle}}\left\{ \bullet \right\}$
    ${\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$ ${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{13}}} \right]$ 0 2 ${\mathop{\rm Re}\nolimits} \left[ {{T_{13}}} \right] + j{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}} \right]$
    ${\mathop{\rm Re}\nolimits} \left[ {{T_{13}}\left( \theta \right)} \right]$ ${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{13}}} \right]$ 0 2 $ - {\mathop{\rm Re}\nolimits} \left[ {{T_{12}}} \right] + j{\mathop{\rm Re}\nolimits} \left[ {{T_{13}}} \right]$
    ${\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$ ${{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{13}}} \right]$ 0 2 ${\mathop{\rm Im}\nolimits} \left[ {{T_{13}}} \right] + j{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}} \right]$
    ${\mathop{\rm Im}\nolimits} \left[ {{T_{13}}\left( \theta \right)} \right]$ ${{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{13}}} \right]$ 0 2 $ - {\mathop{\rm Im}\nolimits} \left[ {{T_{12}}} \right] + j{\mathop{\rm Im}\nolimits} \left[ {{T_{13}}} \right]$
    ${\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$ $\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$ 0 4 $\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right) + j{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right]$
    ${T_{22}}\left( \theta \right)$ $\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$ $\frac{1}{2}\left( {{T_{22}} + {T_{33}}} \right)$ 4 ${\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] + j\frac{1}{2}\left( {{T_{22}} - {T_{33}}} \right)$
    ${T_{33}}\left( \theta \right)$ $\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$ $\frac{1}{2}\left( {{T_{22}} + {T_{33}}} \right)$ 4 $ - {\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] + j\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right)$
    ${\left| {{T_{12}}\left( \theta \right)} \right|^2}$ ${\rm R}{{\rm e}^2}\left[ {{T_{12}}T_{13}^*} \right] + \frac{1}{4}{\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)^2}$ $\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} + {{\left| {{T_{13}}} \right|}^2}} \right)$ 4 ${\mathop{\rm Re}\nolimits} \left[ {{T_{12}}T_{13}^*} \right] + j\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} - {{\left| {{T_{13}}} \right|}^2}} \right)$
    ${\left| {{T_{13}}\left( \theta \right)} \right|^2}$ ${\rm R}{{\rm e}^2}\left[ {{T_{12}}T_{13}^*} \right] + \frac{1}{4}{\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)^2}$ $\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} + {{\left| {{T_{13}}} \right|}^2}} \right)$ 4 $ - {\mathop{\rm Re}\nolimits} \left[ {{T_{12}}T_{13}^*} \right] + j\frac{1}{2}\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)$
    ${\left| {{T_{23}}\left( \theta \right)} \right|^2}$ $\frac{1}{4}{\left\{ {\frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]} \right\}^2}$ $\begin{array}{l}\frac{1}{2}\left\{ {\frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]} \right\}\\ + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{23}}} \right]\end{array}$ 8 $\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right){\rm Re}\left[ {{T_{23}}} \right] \\+ j\frac{1}{2}\left[ {{{{\mathop{\rm Re}\nolimits} }^2}\left[ {{T_{23}}} \right] - \frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2}} \right]$
    下载: 导出CSV

    表  2  针对不同极化SAR实测数据的特征优选结果

    Table  2.   Selected features for different PolSAR data

    实测数据 优选所得旋转域极化特征(相应地物类别对的个数)
    AIRSAR ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]\left( {18} \right)$, ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]\left( {15} \right)$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]\left( {71} \right)$, $B\_{T_{33}}\left( \theta \right)\left( {1} \right)$
    UAVSAR 6月17日 ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(12), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(4)
    6月22日 ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(14), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(2)
    7月03日 ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(3), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(18)
    7月17日 ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(7), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(9)
    下载: 导出CSV

    表  3  两种方法所得AIRSAR数据15类地物及总体的分类精度(%)

    Table  3.   Classification accuracy of different terrains in AIRSAR data using two methods (%)

    地物 传统方法 本文方法
    茎豆 97.2 98.0
    豌豆 93.7 96.9
    森林 92.6 93.7
    苜蓿 96.8 96.6
    小麦1 88.7 85.9
    甜菜 93.8 93.8
    土豆 92.6 93.3
    裸地 95.5 87.2
    草地 59.3 77.3
    油菜籽 83.9 88.0
    大麦 92.6 91.5
    小麦2 89.2 89.4
    小麦3 94.3 95.9
    水域 98.0 98.5
    建筑物 84.9 83.2
    总体精度 91.1 92.3
    下载: 导出CSV

    表  4  两种方法所得多时相UAVSAR数据7类地物及总体的分类精度 (%)

    Table  4.   The classification accuracy of different terrains in multi-temporal UAVSAR data using two methods (%)

    日期 方法 阔叶林 草料 大豆 玉米 小麦 油菜籽 燕麦 总体
    6月17日 传统 98.47 62.24 92.64 96.12 93.63 91.70 86.37 90.19
    本文 98.49 81.65 96.76 98.19 96.08 92.25 96.32 94.98
    6月22日 传统 98.05 61.38 94.14 97.30 97.89 93.82 77.29 90.75
    本文 97.96 72.60 96.86 98.18 97.07 96.84 95.13 95.12
    7月3日 传统 97.41 54.38 90.45 98.89 68.75 98.81 63.46 80.87
    本文 97.77 76.68 98.12 99.08 96.95 98.93 94.22 95.99
    7月17日 传统 96.86 64.51 97.38 99.78 84.76 92.19 82.98 89.39
    本文 97.27 93.15 99.31 99.58 94.73 99.71 92.16 96.78
    平均 传统 97.70 60.63 93.65 98.02 86.26 94.13 77.53 87.80
    本文 97.87 81.02 97.76 98.76 96.21 96.93 94.46 95.72
    下载: 导出CSV
  • [1] Lee J S and Pottier E. Polarimetric Radar Imaging: From Basics to Applications[M]. Boca Raton, FL, USA, CRC Press, 2009.
    [2] 吴永辉, 计科峰, 郁文贤. SVM全极化SAR图像分类中的特征选择[J]. 信号处理, 2007, 23(6): 877–881. http://www.cnki.com.cn/Article/CJFDTOTAL-XXCN200706021.htm

    Wu Yonghui, Ji Kefeng, and Yu Wenxian. A new feature selection algorithm for SVM-based fully polarimetric SAR image classification[J]. Signal Processing, 2007, 23(6): 877–881. http://www.cnki.com.cn/Article/CJFDTOTAL-XXCN200706021.htm
    [3] Cloude S R and Pottier E. A review of target decomposition theorems in radar polarimetry[J]. IEEE Transactions on Geoscience and Remote Sensing, 1996, 34(2): 498–518. doi:  10.1109/36.485127
    [4] Freeman A and Durden S L. A three-component scattering model for polarimetric SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(3): 963–973. doi:  10.1109/36.673687
    [5] Yamaguchi Y, Moriyama T, Ishido M, et al.. Four-component scattering model for polarimetric SAR image decomposition[J]. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(8): 1699–1706. doi:  10.1109/TGRS.2005.852084
    [6] Chen S W, Wang X S, Xiao S P, et al.. General polarimetric model-based decomposition for coherency matrix[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(3): 1843–1855. doi:  10.1109/TGRS.2013.2255615
    [7] Chen S W, Li Y Z, Wang X S, et al.. Modeling and interpretation of scattering mechanisms in polarimetric synthetic aperture radar: Advances and perspectives[J]. IEEE Signal Processing Magazine, 2014, 31(4): 79–89. doi:  10.1109/MSP.2014.2312099
    [8] Cloude S R and Pottier E. An entropy based classification scheme for land applications of polarimetric SARs[J]. IEEE Transactions on Geoscience and Remote Sensing, 1997, 35(1): 68–78. doi:  10.1109/36.551935
    [9] Lee J S, Grunes M R, Ainsworth T L, et al.. Unsupervised classification of polarimetric SAR images by applying target decomposition and complex Wishart distribution[J]. IEEE Transactions on Geoscience and Remote Sensing, 1999, 37(5): 2249–2258. doi:  10.1109/36.789621
    [10] 赵一博. 基于SVM和RBF神经网络的极化SAR图像分类方法研究[D]. [硕士论文], 国防科学技术大学, 2013.

    Zhao Yibo. Classification of polarimetric SAR image based on SVM and RBF neural network[D]. [Master dissertation], National University of Defense Technology, 2013.
    [11] 吴永辉, 计科峰, 郁文贤. 基于H-alpha和改进C-均值的全极化SAR图像非监督分类[J]. 电子与信息学报, 2007, 29(1): 30–34. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX200701007.htm

    Wu Yonghui, Ji Kefeng, and Yu Wenxian. Unsupervised classification of fully polarimetric SAR image using H-alpha decomposition and modified C-Mean algorithm[J]. Journal of Electronics & Information Technology, 2007, 29(1): 30–34. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX200701007.htm
    [12] Chen S W, Wang X S, and Sato M. Uniform polarimetric matrix rotation theory and its applications[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(8): 4756–4770. doi:  10.1109/TGRS.2013.2284359
    [13] Chen S W, Li Y Z, and Wang X S. Crop discrimination based on polarimetric correlation coefficients optimization for PolSAR data[J]. International Journal of Remote Sensing, 2015, 36(16): 4233–4249. doi:  10.1080/01431161.2015.1079345
    [14] Xiao S P, Chen S W, Chang Y L, et al.. Polarimetric coherence optimization and its application for manmade target extraction in PolSAR data[J]. IEICE Transactions on Electronics, 2014, 97(6): 566–574.
    [15] Chang C C and Lin C J. LIBSVM: Alibrary for support vector machines[J]. ACM Transactions on Intelligent Systems and Technology, 2011, 2(3): 389–396.
    [16] Chen S W, Wang X S, and Sato M. PolInSAR complex coherence estimation based on covariance matrix similarity test[J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(11): 4699–4709. doi:  10.1109/TGRS.2012.2192937
  • [1] 崔兴超, 粟毅, 陈思伟.  融合极化旋转域特征和超像素技术的极化SAR舰船检测 . 雷达学报, 2021, 10(1): 35-48. doi: 10.12000/JR20147
    [2] 王雪松, 陈思伟.  合成孔径雷达极化成像解译识别技术的进展与展望 . 雷达学报, 2020, 9(2): 259-276. doi: 10.12000/JR19109
    [3] 陈世超, 高鹤婷, 罗丰.  基于极化联合特征的海面目标检测方法 . 雷达学报, 2020, 9(4): 664-673. doi: 10.12000/JR20072
    [4] 邹焕新, 李美霖, 马倩, 孙嘉赤, 曹旭, 秦先祥.  一种基于张量积扩散的非监督极化SAR图像地物分类方法 . 雷达学报, 2019, 8(4): 436-447. doi: 10.12000/JR19057
    [5] 胡涛, 李卫华, 秦先祥, 王鹏, 余旺盛, 李军.  基于深度卷积神经网络和条件随机场模型的PolSAR图像地物分类方法 . 雷达学报, 2019, 8(4): 471-478. doi: 10.12000/JR18065
    [6] 肖东凌, 刘畅.  基于精调的膨胀编组-交叉CNN的PolSAR地物分类 . 雷达学报, 2019, 8(4): 479-489. doi: 10.12000/JR19039
    [7] 张腊梅, 张思雨, 董洪伟, 朱厦.  基于Pinball损失函数支持向量机的极化SAR图像鲁棒分类 . 雷达学报, 2019, 8(4): 448-457. doi: 10.12000/JR19055
    [8] 滑文强, 王爽, 郭岩河, 谢雯.  基于邻域最小生成树的半监督极化SAR图像分类方法 . 雷达学报, 2019, 8(4): 458-470. doi: 10.12000/JR18104
    [9] 折小强, 仇晓兰, 雷斌, 张薇, 卢晓军.  一种极化熵结合混合GEV模型的全极化SAR潮间带区域地物分类方法 . 雷达学报, 2017, 6(5): 554-563. doi: 10.12000/JR16149
    [10] 邹焕新, 罗天成, 张月, 周石琳.  基于组合条件随机场的极化SAR图像监督地物分类 . 雷达学报, 2017, 6(5): 541-553. doi: 10.12000/JR16109
    [11] 徐丰, 王海鹏, 金亚秋.  深度学习在SAR目标识别与地物分类中的应用 . 雷达学报, 2017, 6(2): 136-148. doi: 10.12000/JR16130
    [12] 胡丁晟, 仇晓兰, 雷斌, 徐丰.  极化串扰对基于Cloude分解的地物散射机制特征量影响分析 . 雷达学报, 2017, 6(2): 221-228. doi: 10.12000/JR16129
    [13] 陈思伟, 李永祯, 王雪松, 肖顺平.  极化SAR目标散射旋转域解译理论与应用 . 雷达学报, 2017, 6(5): 442-455. doi: 10.12000/JR17033
    [14] 田壮壮, 占荣辉, 胡杰民, 张军.  基于卷积神经网络的SAR图像目标识别研究 . 雷达学报, 2016, 5(3): 320-325. doi: 10.12000/JR16037
    [15] 黄晓菁, 杨祥立, 黄平平, 杨文.  基于原型理论的极化SAR图像特征表达 . 雷达学报, 2016, 5(2): 208-216. doi: 10.12000/JR15071
    [16] 孙勋, 黄平平, 涂尚坦, 杨祥立.  利用多特征融合和集成学习的极化SAR图像分类 . 雷达学报, 2016, 5(6): 692-700. doi: 10.12000/JR15132
    [17] 邢艳肖, 张毅, 李宁, 王宇, 胡桂香.  一种联合特征值信息的全极化SAR图像监督分类方法 . 雷达学报, 2016, 5(2): 217-227. doi: 10.12000/JR16019
    [18] 杨汝良, 戴博伟, 李海英.  极化合成孔径雷达极化层次和系统工作方式 . 雷达学报, 2016, 5(2): 132-142. doi: 10.12000/JR16013
    [19] 邵璐熠, 洪文.  基于二维极化特征的PolSAR图像决策分类 . 雷达学报, 2016, 5(6): 681-691. doi: 10.12000/JR16002
    [20] 滑文强, 王爽, 侯彪.  基于半监督学习的SVM-Wishart极化SAR图像分类方法 . 雷达学报, 2015, 4(1): 93-98. doi: 10.12000/JR14138
  • 加载中
图(7) / 表 (4)
计量
  • 文章访问数:  1141
  • HTML全文浏览量:  369
  • PDF下载量:  753
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-11-30
  • 修回日期:  2017-01-24
  • 网络出版日期:  2017-04-07
  • 刊出日期:  2017-10-28

结合旋转域极化特征的极化SAR地物分类

doi: 10.12000/JR16131
    基金项目:  国家自然科学基金(41301490, 61490692)
    作者简介:

    陶臣嵩(1993–),男,广西人,国防科学技术大学电子科学与工程学院在读硕士研究生,主要研究方向为成像雷达极化信息的解译与处理、极化特征的提取与优选等。E-mail: taochensongnudt@163.com

    陈思伟(1984–),男,四川人,博士,国防科学技术大学电子科学与工程学院讲师,主要研究方向包括雷达极化信息处理、成像雷达信息处理、目标散射建模与解译、微波遥感大数据处理与应用、微波遥感环境与灾害应用研究等。E-mail: chenswnudt@163.com

    李永祯(1977–),男,内蒙古人,博士后,国防科学技术大学电子科学与工程学院研究员,电子信息系统复杂电磁环境效应国家重点实验室副主任,主要研究方向为新体制雷达与电子对抗。E-mail: e0061@sina.com

    肖顺平(1964–),男,江西人,博士,国防科学技术大学电子科学与工程学院教授,博士生导师,电子信息系统复杂电磁环境效应国家重点实验室主任,国家高技术863专家,总装备部仿真专业组专家,中国电子学会高级会员,主要研究方向包括雷达极化信息处理、电子信息系统仿真评估技术、雷达目标识别等

    通讯作者: 陈思伟   chenswnudt@163.com
  • 中图分类号: TN957.52

摘要: 地物分类是极化合成孔径雷达(SAR)图像理解与解译的重要应用方向。利用H/A/α/SPAN等旋转不变特征参数的极化SAR地物分类是一种常用的分类方法。然而,目标的后向散射响应与其方位取向等姿态密切相关,极易引起散射机理的解译模糊,也限制了仅使用旋转不变特征参数作为分类特征集的极化SAR地物分类方法的精度。针对这一问题,有文献提出了在绕雷达视线的旋转域中解译目标散射特性的统一的极化矩阵旋转理论,并导出了一系列旋转域极化特征,刻画目标旋转域隐含信息。基于该理论,该文将旋转域极化特征用于极化SAR地物辨识与分类,并发展了一种结合旋转域极化特征与旋转不变特征H/A/α/SPAN的极化SAR地物分类方法。该方法将优选的旋转域极化特征参数和H/A/α/SPAN作为支持向量机(SVM)分类器的输入,利用两类特征对不同地物类别区分辨识能力的互补,以达到更好的分类性能。对AIRSAR和UAVSAR实测数据的对比实验表明,相较于仅以H/A/α/SPAN作为SVM分类器输入的传统方法,该方法得到更好的分类精度和稳健性。其中,对于AIRSAR数据十五类地物的分类,该方法总体分类精度达到92.3%,优于传统方法的91.1%。此外,对于多时相UAVSAR数据七类地物的分类,该方法平均总体分类精度达到95.72%,显著优于传统方法的87.80%,验证了该方法对多时相数据的稳健性。该文研究进一步证实了通过深入挖掘旋转域中目标极化散射信息能够为极化SAR图像的解译与应用提供新的可行途径。

English Abstract

陶臣嵩, 陈思伟, 李永祯, 肖顺平. 结合旋转域极化特征的极化SAR地物分类[J]. 雷达学报, 2017, 6(5): 524-532. doi: 10.12000/JR16131
引用本文: 陶臣嵩, 陈思伟, 李永祯, 肖顺平. 结合旋转域极化特征的极化SAR地物分类[J]. 雷达学报, 2017, 6(5): 524-532. doi: 10.12000/JR16131
Tao Chensong, Chen Siwei, Li Yongzhen, Xiao Shunping. Polarimetric SAR Terrain Classification Using Polarimetric Features Derived from Rotation Domain[J]. Journal of Radars, 2017, 6(5): 524-532. doi: 10.12000/JR16131
Citation: Tao Chensong, Chen Siwei, Li Yongzhen, Xiao Shunping. Polarimetric SAR Terrain Classification Using Polarimetric Features Derived from Rotation Domain[J]. Journal of Radars, 2017, 6(5): 524-532. doi: 10.12000/JR16131
    • 极化合成孔径雷达(Polarimetric Synthetic Aperture Radar, PolSAR)具有全天时和几乎全天候的工作能力,通过收发极化状态正交的电磁波以获取目标的全极化散射信息[1]。地物分类是农作物生长监控、农村与城市用地普查、环境监测等应用领域的共性基础问题,也是极化SAR图像理解与解译的重要应用方向。高精度的地物分类结果能够为上述应用领域提供可靠的信息支撑。

      通常,提高极化SAR地物分类精度主要有两种途径[2]。第1种途径专注于极化特征的挖掘与优选,通过精细化的极化散射机理建模与解译,从全极化信息中提取出对不同地物类别具有更强区分度的特征。常用的极化散射机理解译方法有基于特征值分解的方法和基于模型分解的方法。基于这些极化目标分解方法所得到的极化特征参数经常被用于极化SAR地物分类,例如Cloude-Pottier分解所得的极化熵/极化平均角/极化反熵(H/ $\alpha $/A)参数[3],Freeman-Durden分解[4]、Yamaguchi分解[5]和近年来提出的精细化极化目标分解[6]所得的各散射机理的散射能量参数(如奇次散射、偶次散射、体散射、螺旋散射等)[7]。第2种途径则从分类器入手,使用性能更好的分类器,以对现有的极化特征进行充分利用。常用的分类器包括C均值分类器、Wishart分类器、支持向量机(Support Vector Machine, SVM)分类器、随机森林分类器、神经网络分类器以及近来年在诸多领域取得成功应用的以卷积神经网络为代表的深度学习分类方法等[811]。当然,对特征和分类器同时进行优化和优选也是提高极化SAR地物分类精度的有效途径。

      在传统基于特征的极化SAR地物分类中,具有旋转不变特性的极化特征参数得到了广泛应用。例如,基于H/ $\alpha $/A和总散射能量SPAN的极化SAR地物分类就是一种常用的分类方法。然而,目标的极化响应与目标和SAR的相对几何关系密切相关。同一目标在不同方位取向下,其后向散射可以是显著不同的。同时,不同目标在某些特定方位取向下,其后向散射又是十分相似的。例如,具有不同方位取向的建筑物与森林等植被就是极化SAR图像解译的难点。这是诸多传统极化目标分解方法存在散射机理解译模糊的重要原因之一,同时也限制了基于旋转不变极化特征参数的传统分类方法所得精度的进一步提升。为避免这种解译模糊,一种思路是构建更精细化的目标散射模型和精细化的极化目标分解方法。而另一种思路则是挖掘利用目标方位取向与其后向散射机理之间的隐含关系。文献[12]提出的统一的极化矩阵旋转理论就是一种代表性的方法。该方法提出了在绕雷达视线的旋转域中理解目标散射特性的新思路,并导出了一系列旋转域极化特征。部分旋转域极化特征参数已经在农作物辨识[13]、目标对比增强[12]、人造目标提取[14]等领域获得了成功应用。

      由于这些旋转域极化特征包含有目标在旋转域中隐含的极化散射信息,且与其方位取向具有一定关系。若将它们与传统的旋转不变极化特征参数于H/ $\alpha $/A/SPAN联合作为地物分类特征集,则从极化特征挖掘的角度来看,两类不同的极化特征对于不同地物类别的区分能力势必会形成一定程度的互补,进而使分类精度得到进一步提升。基于这一思路,本文提出了一种结合旋转域极化特征与旋转不变特征H/ $\alpha $/A/SPAN的极化SAR地物分类方法。具体即基于不同地物类别样本集类间距最大的特征优选准则,以部分优选的旋转域极化特征参数与H/ $\alpha $/A/SPAN联合作为地物分类所用特征,并选用性能较为稳定的SVM[15]作为分类器进行分类处理。由于该分类方法额外使用了目标在方位取向方面的隐含信息,故相较于仅使用旋转不变特征H/ $\alpha $/A/SPAN作为输入的SVM分类器[10],其能够达到更优的分类性能表现。

      本文第2节简要介绍了统一的极化矩阵旋转理论及其所导出的旋转域极化特征参数;第3节提出结合旋转域极化特征的极化SAR地物分类方法;第4节基于AIRSAR和多时相UAVSAR实测数据开展了地物分类对比实验及分析;第5节总结本文方法并对后续研究工作进行展望。

    • 极化SAR获得的目标全极化信息可以通过极化相干矩阵T表示。满足互易性原理时,极化相干矩阵T可以表示为:

      $${{T}} = \left\langle {{{{k}}_{\rm{P}}}{{k}}_{\rm{P}}^{\rm{H}}} \right\rangle = \left[ {\begin{array}{*{20}{c}}{{T_{{\rm{11}}}}}&{{T_{{\rm{12}}}}}&{{T_{{\rm{13}}}}}\\{{T_{{\rm{21}}}}}&{{T_{{\rm{22}}}}}&{{T_{{\rm{23}}}}}\\{{T_{{\rm{31}}}}}&{{T_{{\rm{32}}}}}&{{T_{{\rm{33}}}}}\end{array}} \right]$$ (1)

      其中, ${{{k}}_{\rm{P}}} = \displaystyle\frac{1}{{\sqrt 2 }}{\left[ {\begin{array}{*{20}{c}}{{S_{{\rm{HH}}}} + {S_{{\rm{VV}}}}}&{{S_{{\rm{HH}}}} - {S_{{\rm{VV}}}}}&{2{S_{{\rm{HV}}}}}\end{array}} \right]^{\rm{T}}}$为Pauli散射矢量。 ${S_{{\rm{HV}}}}$为以垂直极化天线发射并以水平极化天线接收条件下的散射系数, ${{{k}}_{\rm{P}}}$中其它元素可类似定义。 $\left\langle { \ } \right\rangle $表示集合平均。 ${T_{ij}}$则表示极化相干矩阵 ${{T}}$中第i行第j列所对应的元素。

      将极化相干矩阵 ${{T}}$绕雷达视线进行旋转处理,则可得到旋转域中极化相干矩阵的表达式为:

      $${{T}}\left( \theta \right) = {{{k}}_{\rm{P}}}\left( \theta \right){{k}}_{\rm{P}}^{\rm{H}}\left( \theta \right) = {{{R}}_{\rm{3}}}\left( \theta \right){{TR}}_{\rm{3}}^{\rm{H}}\left( \theta \right)$$ (2)

      其中,旋转矩阵为:

      $${{{R}}_{\rm{3}}}\left( \theta \right) = \left[ {\begin{array}{*{20}{c}}1&0&0\\0&{\cos 2\theta }&{\sin 2\theta }\\0&{ - \sin 2\theta }&{\cos 2\theta }\end{array}} \right]$$ (3)

      在旋转域中极化相干矩阵 ${{T}}\left( \theta \right)$的每个元素经过相应的数学变换即可被统一地由一个正弦函数进行表征[12]

      $$f\left( \theta \right) = A\sin \left[ {\omega \left( {\theta + {\theta _{\rm{0}}}} \right)} \right] + B$$ (4)

      其中,A为振荡幅度,B为振荡中心, $\omega $为角频率, ${\theta _{\rm{0}}}$为初始角度。文献[12]将这4类极化特征参数 $\left\{ {A,B,\omega ,{\theta _0}} \right\}$称为振荡参数集,其完整表征极化相干矩阵的各元素在旋转域中的特性。这样就可以导出一系列旋转域极化特征参数,如表1所示。其中, ${\rm{Angle}}\left\{ a \right\}$表示复数a的相位,相应取值范围为 $\left[ { - {{π}} ,{{π}} } \right]$。

      表 1  旋转域极化特征参数[12]

      Table 1.  Polarimetric feature parameters derived from rotation domain[12]

      散射矩阵元素项 $A = \sqrt \bullet $ B $\omega $ ${\theta _0} = \frac{1}{\omega }{\rm{Angle}}\left\{ \bullet \right\}$
      ${\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$ ${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{13}}} \right]$ 0 2 ${\mathop{\rm Re}\nolimits} \left[ {{T_{13}}} \right] + j{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}} \right]$
      ${\mathop{\rm Re}\nolimits} \left[ {{T_{13}}\left( \theta \right)} \right]$ ${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{13}}} \right]$ 0 2 $ - {\mathop{\rm Re}\nolimits} \left[ {{T_{12}}} \right] + j{\mathop{\rm Re}\nolimits} \left[ {{T_{13}}} \right]$
      ${\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$ ${{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{13}}} \right]$ 0 2 ${\mathop{\rm Im}\nolimits} \left[ {{T_{13}}} \right] + j{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}} \right]$
      ${\mathop{\rm Im}\nolimits} \left[ {{T_{13}}\left( \theta \right)} \right]$ ${{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{13}}} \right]$ 0 2 $ - {\mathop{\rm Im}\nolimits} \left[ {{T_{12}}} \right] + j{\mathop{\rm Im}\nolimits} \left[ {{T_{13}}} \right]$
      ${\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$ $\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$ 0 4 $\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right) + j{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right]$
      ${T_{22}}\left( \theta \right)$ $\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$ $\frac{1}{2}\left( {{T_{22}} + {T_{33}}} \right)$ 4 ${\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] + j\frac{1}{2}\left( {{T_{22}} - {T_{33}}} \right)$
      ${T_{33}}\left( \theta \right)$ $\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$ $\frac{1}{2}\left( {{T_{22}} + {T_{33}}} \right)$ 4 $ - {\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] + j\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right)$
      ${\left| {{T_{12}}\left( \theta \right)} \right|^2}$ ${\rm R}{{\rm e}^2}\left[ {{T_{12}}T_{13}^*} \right] + \frac{1}{4}{\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)^2}$ $\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} + {{\left| {{T_{13}}} \right|}^2}} \right)$ 4 ${\mathop{\rm Re}\nolimits} \left[ {{T_{12}}T_{13}^*} \right] + j\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} - {{\left| {{T_{13}}} \right|}^2}} \right)$
      ${\left| {{T_{13}}\left( \theta \right)} \right|^2}$ ${\rm R}{{\rm e}^2}\left[ {{T_{12}}T_{13}^*} \right] + \frac{1}{4}{\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)^2}$ $\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} + {{\left| {{T_{13}}} \right|}^2}} \right)$ 4 $ - {\mathop{\rm Re}\nolimits} \left[ {{T_{12}}T_{13}^*} \right] + j\frac{1}{2}\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)$
      ${\left| {{T_{23}}\left( \theta \right)} \right|^2}$ $\frac{1}{4}{\left\{ {\frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]} \right\}^2}$ $\begin{array}{l}\frac{1}{2}\left\{ {\frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]} \right\}\\ + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{23}}} \right]\end{array}$ 8 $\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right){\rm Re}\left[ {{T_{23}}} \right] \\+ j\frac{1}{2}\left[ {{{{\mathop{\rm Re}\nolimits} }^2}\left[ {{T_{23}}} \right] - \frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2}} \right]$

      基于上述振荡参数集,文献[12]还导出了一系列的极化角参数集,如极化零角参数、极化最大化角参数以及极化最小化角参数等。其中,极化零角参数的定义为在绕雷达视线的旋转域中使极化相干矩阵某元素取值为零的旋转角,即:

      $$f\left( \theta \right) \!=\! A{\rm{sin}}\left[ {\omega \left( {{\theta _{{\rm{null}}}} \!+\! {\theta _0}} \right)} \right] \!+\! B \!=\!\! 0 \Rightarrow {\theta _{{\rm{null}}}} \!=\! - {\theta _0}$$ (5)

      其中, ${\theta _{{\rm{null}}}}$即极化零角参数。由于表1中相互独立的5个初始角度 ${\theta _0}$分别为 ${\theta _0}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _0}\_{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _0}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$, ${\theta _0}\_{\left| {{T_{12}}\left( \theta \right)} \right|^2}$和 ${\theta _0}\_{\left| {{T_{23}}\left( \theta \right)} \right|^2}$,故相应的极化零角参数有 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[\! {{T_{12}}\left( \theta \right)} \!\right]$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Im}\nolimits} \left[\! {{T_{12}}\left( \theta \right)} \!\right]$, ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ \!{{T_{23}}\left( \theta \right)} \!\right]$, ${\theta _{{\rm{null}}}}\_{\left| {{T_{12}}\left( \theta \right)} \right|^2}$和 ${\theta _{{\rm{null}}}}\_{\left| {{T_{23}}\left( \theta \right)} \right|^2}$。由文献[12]可知,各初始角度与其相应极化零角参数所包含的极化信息是相互等价的,且极化零角参数具有相对明确的物理意义,故在本文的后续部分均以极化零角参数代替相应的初始角度。

    • 文献[12]使用极化零角参数 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{12}}\left( \theta \right)} \right]$和 ${\theta _{{\rm{null}}}}\_{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$的组合能够成功辨识7类不同农作物,初步证实了极化零角参数集对于不同地物类别具有较好的区分能力。在此基础上,本文挖掘利用旋转域极化特征所蕴含目标在旋转域中的隐含信息,并将其应用于极化SAR地物分类。

      在此之前,需要基于地物分类的应用背景对众多的旋转域极化特征进行优选处理。在文献[12]所导出的一系列旋转域极化特征之中,以不同地物类别样本集相互之间的“类间距最大化”为准则,进行相应的旋转域极化特征优选。具体步骤为:首先对各旋转域极化特征参数进行归一化处理;然后将不同的地物类别两两组合形成若干的地物类别对;接着针对各地物类别对,以其中两地物类别之间的类间距为标准,优选出使其取值达到最大的旋转域极化特征,则每个地物类别对均对应于一个优选的旋转域极化特征;最后,将各地物类别对的优选结果进行“取并集”处理,进而得到最终的优选结果。

      文献[12]所导出相互独立的旋转域极化特征共有12个,分别为 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{23}}\left( \theta \right)} \right]$, ${\theta _{{\rm{null}}}}\_{\left| {{T_{12}}\left( \theta \right)} \right|^2}$, ${\theta _{{\rm{null}}}}\_{\left| {{T_{23}}\left( \theta \right)} \right|^2}$, $A\_{\rm Re}\left[ {{T_{12}}\left( \theta \right)} \right]$, $A\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$, $A\_{T_{12}}\left( \theta \right)$, A_ ${T_{23}}\left( \theta \right)$, $B\_{T_{12}}\left( \theta \right)$, $B\_{T_{33}}\left( \theta \right)$, $B\_{T_{23}}\left( \theta \right)$。针对之后实验部分所使用的AIRSAR数据(15类地物,两两组合形成105个地物类别对;其它说明见4.1节)以及多时相UAVSAR数据(7类地物,两两组合形成21个地物类别对;4个数据获取日期;其它说明见4.2节),上述特征优选流程所得结果如表2所示。

      表 2  针对不同极化SAR实测数据的特征优选结果

      Table 2.  Selected features for different PolSAR data

      实测数据 优选所得旋转域极化特征(相应地物类别对的个数)
      AIRSAR ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]\left( {18} \right)$, ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]\left( {15} \right)$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]\left( {71} \right)$, $B\_{T_{33}}\left( \theta \right)\left( {1} \right)$
      UAVSAR 6月17日 ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(12), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(4)
      6月22日 ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(14), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(2)
      7月03日 ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(3), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(18)
      7月17日 ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(7), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(9)

      综合考虑表2中的优选结果,并在追求较高地物分类精度的同时,将两组实测数据优选得到的旋转域极化特征进行统一,故本文优选部分的最终结果为3个极化零角参数,即 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$和 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{23}}\left( \theta \right)} \right]$。

    • 为了将目标在旋转域中的隐含信息充分利用在极化SAR地物分类中,同时又发挥传统的旋转不变极化特征参数H/A/ $\alpha $/SPAN在极化散射机理解译方面的优点,本文提出了一种结合旋转域极化特征的极化SAR地物分类方法,其流程图如图1所示,相应的具体操作如下:

      图  1  本文方法具体流程图

      Figure 1.  Flowchart of proposed method

      (1) 在进行Cloude-Pottier分解之前,需要对极化SAR数据进行相干斑滤波处理。本文采用新近提出的一种基于矩阵相似性检验的SimiTest自适应相干斑滤波方法[16]对极化SAR数据进行滤波预处理。

      (2) 基于滤波后的极化相干矩阵,计算总散射能量SPAN。

      (3) 同样地,基于滤波后的极化相干矩阵,进行Cloude-Pottier分解,得到极化特征量H/ $\alpha $/A

      (4) 同时,将滤波后的极化相干矩阵绕雷达视线旋转,计算上述优选部分所得的3个极化零角参数。

      (5) 对上述7个极化特征参数分别进行归一化处理,以作为地物分类特征集输入至SVM分类器。

      (6) 通过SVM相应的训练与测试过程,实现对不同地物类别的分类处理。

    • 为了验证新极化特征(即3个旋转域极化零角参数)的引入对于传统地物分类方法性能的提升作用,在对极化相干矩阵中全部极化信息进行利用的前提之下,将本文方法与仅使用旋转不变特征H/A/ $\alpha $/SPAN作为SVM分类器输入的传统方法进行对比。首先使用AIRSAR数据15类地物的分类验证本文方法的分类性能,再使用多时相UAVSAR数据7类地物的分类进一步验证本文方法对多时相数据的稳健性。在对此两组数据分别进行SimiTest相干斑滤波[16]时,所用滑窗大小均为15×15。对SVM分类器,各类地物样本的一半用于训练,另一半用于测试。

    • 本文首先使用NASA/JPL AIRSAR系统在荷兰Flevoland地区所获取的L波段全极化SAR数据进行地物分类实验。该数据方位向分辨率为12.1 m,距离向分辨率为6.6 m,所用区域大小为736×1010。SimiTest相干斑滤波后的Pauli RGB图如图2(a)所示。该区域的真值图如图2(b)所示,其中主要包含茎豆、豌豆、森林、苜蓿、小麦1、甜菜、土豆、裸地、草地、油菜籽、大麦、小麦2、小麦3、水域以及建筑物等15类地物。

      图  2  AIRSAR数据

      Figure 2.  AIRSAR data

      使用传统方法和本文方法分别对滤波后的数据进行分类处理,所得结果如图3所示。

      图  3  AIRSAR数据的分类结果

      Figure 3.  Classification results of AIRSAR data

      两种方法对AIRSAR数据15类地物分类处理所得精度如表3所示。通过比较可知,本文方法得到的总体分类精度为92.3%,优于传统方法91.1%的分类精度。且本文方法对草地77.3%的分类精度相较于传统方法的59.3%提升了18个百分点。另外,由于SVM分类器所用分类策略以总体分类精度的最大化为目标,无法保证单一地物类别的分类精度均达到最优。例如,本文方法在苜蓿、小麦1、裸地、大麦以及建筑物等5种地物类别区域所得分类精度均不及传统方法。针对其中分类精度差距最大(约8.3%)的裸地,由于其相应区域的主要散射机制为“面散射”,不同方位取向对其后向散射的影响较小,使用传统的旋转不变极化特征已经能较好地对其进行区分与辨识,本文方法额外引入的3个旋转域极化零角参数可能造成了分类信息的冗余,进而导致所得分类精度的较大幅度下降。

      表 3  两种方法所得AIRSAR数据15类地物及总体的分类精度(%)

      Table 3.  Classification accuracy of different terrains in AIRSAR data using two methods (%)

      地物 传统方法 本文方法
      茎豆 97.2 98.0
      豌豆 93.7 96.9
      森林 92.6 93.7
      苜蓿 96.8 96.6
      小麦1 88.7 85.9
      甜菜 93.8 93.8
      土豆 92.6 93.3
      裸地 95.5 87.2
      草地 59.3 77.3
      油菜籽 83.9 88.0
      大麦 92.6 91.5
      小麦2 89.2 89.4
      小麦3 94.3 95.9
      水域 98.0 98.5
      建筑物 84.9 83.2
      总体精度 91.1 92.3
    • 本文使用NASA/JPL UAVSAR系统在加拿大Manitoba地区所获取的多时相L波段全极化SAR数据进行地物分类实验。该数据方位向分辨率为7 m,距离向分辨率为5 m,所用区域大小为1325×1011。多时相极化SAR数据分别获取于6月17日、6月22日、7月3日以及7月17日。SimiTest相干斑滤波处理之后多时相极化SAR数据对应的Pauli RGB图如图4所示。该区域的主要地物类型是以谷物和油种产品为代表的混合型牧场农作物。相应的真值图如图5所示,其中主要包含阔叶林、草料、大豆、玉米、小麦、油菜籽以及燕麦等7类地物。

      图  4  多时相UAVSAR数据滤波后Pauli RGB图

      Figure 4.  Filtered Pauli RGB images of multi-temporal UAVSAR data

      图  5  所用区域的真值图

      Figure 5.  Gound truth of the multi-temporal data

      使用传统方法和本文方法分别对滤波后的多时相极化SAR数据进行相互独立的分类处理,所得结果分别如图6图7所示。

      图  6  传统方法对多时相UAVSAR数据分类结果

      Figure 6.  Classification results of multi-temporal UAVSAR data using conventional method

      图  7  本文方法对多时相UAVSAR数据分类结果

      Figure 7.  Classification results of multi-temporal UAVSAR data using proposed method

      图6(c)图7(c)所示,基于7月3日获取的数据,传统方法将红色圆框内小麦与燕麦的绝大部分错分为了大豆,而本文方法在该区域的分类性能相较于前者有显著提升。又如图6(d)图7(d)所示,基于7月17日获取的数据,传统方法将白色圆框内小麦的绝大部分错分为了大豆,而本文方法在该区域的分类精度相较于前者也有较大提升。

      两种方法对多时相UAVSAR数据7类地物分类处理所得精度如表4所示。通过比较可知,对不同日期获取的数据,本文方法所得各类地物及总体的分类精度均优于或相当于传统方法。其中,对6月17日、6月22日、7月3日以及7月17日4个不同日期所获取的数据,本文方法得到的总体分类精度分别为94.98%, 95.12%, 95.99%以及96.78%,而传统方法所得总体分类精度则波动于80.87%至90.75%之间,出现约10%的起伏。具体就小麦和燕麦而言,本文方法得到的分类精度均分别保持在94%和92%以上,而传统方法所得相应分类精度则分别出现了约30%和23%的波动起伏。另外,本文方法95.72%的平均总体分类精度相较于传统方法的87.80%提升了约8个百分点。故本文方法较好的分类性能对于同一系统的多时相数据更具稳健性。

      表 4  两种方法所得多时相UAVSAR数据7类地物及总体的分类精度 (%)

      Table 4.  The classification accuracy of different terrains in multi-temporal UAVSAR data using two methods (%)

      日期 方法 阔叶林 草料 大豆 玉米 小麦 油菜籽 燕麦 总体
      6月17日 传统 98.47 62.24 92.64 96.12 93.63 91.70 86.37 90.19
      本文 98.49 81.65 96.76 98.19 96.08 92.25 96.32 94.98
      6月22日 传统 98.05 61.38 94.14 97.30 97.89 93.82 77.29 90.75
      本文 97.96 72.60 96.86 98.18 97.07 96.84 95.13 95.12
      7月3日 传统 97.41 54.38 90.45 98.89 68.75 98.81 63.46 80.87
      本文 97.77 76.68 98.12 99.08 96.95 98.93 94.22 95.99
      7月17日 传统 96.86 64.51 97.38 99.78 84.76 92.19 82.98 89.39
      本文 97.27 93.15 99.31 99.58 94.73 99.71 92.16 96.78
      平均 传统 97.70 60.63 93.65 98.02 86.26 94.13 77.53 87.80
      本文 97.87 81.02 97.76 98.76 96.21 96.93 94.46 95.72

      另外,对于6月22日所获取数据中的阔叶林和小麦,以及7月17日所获取数据中的玉米,本文方法所得分类精度均略低于传统方法,且分类精度的差距均在1%以内。

      在上述两组相互独立的对比实验所得结果中,本文方法所得分类精度均优于传统方法。故本文方法所表现出的较好分类性能对于不同系统的数据也具有较强稳健性。

    • 目标方位取向对其后向散射响应的直接影响极易引起散射机理的解译模糊,进而限制仅使用旋转不变特征参数作为分类特征集的极化SAR地物分类所得精度。针对这一问题,本文将刻画目标旋转域隐含信息的旋转域极化特征用于极化SAR地物分类,并提出了一种结合旋转域极化特征和旋转不变特征H/A/ $\alpha $/SPAN的极化SAR地物分类方法,该方法将旋转域极化零角参数和H/A/ $\alpha $/SPAN联合作为分类特征集输入至SVM分类器。

      将本文方法与仅使用旋转不变特征H/A/ $\alpha $/SPAN作为SVM分类器输入的传统方法进行比较:对AIRSAR数据15类地物分类而言,本文方法总体分类精度达到92.3%,优于传统方法的91.1%。对多时相UAVSAR数据7类地物分类而言,本文方法平均总体分类精度达到95.72%,显著优于传统方法的87.80%,表明本文方法对同一系统的多时相数据更具稳健性。这两组对比实验也表明本文方法较好的分类性能对于不同系统的数据具有较强稳健性。

      通过对旋转域中目标极化散射信息的深入挖掘,能够为极化SAR图像的解译与应用提供一条新的可行途径。下一步将考虑旋转域极化特征与具有深度学习能力的卷积神经网络等分类器相结合,以实现更高的分类精度。另外,对极化特征参数更优的选择准则及相互融合也是我们未来将要深入研究讨论的内容。

参考文献 (16)

目录

    /

    返回文章
    返回