基于组合条件随机场的极化SAR图像监督地物分类

邹焕新 罗天成 张月 周石琳

邹焕新, 罗天成, 张月, 周石琳. 基于组合条件随机场的极化SAR图像监督地物分类[J]. 雷达学报, 2017, 6(5): 541-553. doi: 10.12000/JR16109
引用本文: 邹焕新, 罗天成, 张月, 周石琳. 基于组合条件随机场的极化SAR图像监督地物分类[J]. 雷达学报, 2017, 6(5): 541-553. doi: 10.12000/JR16109
Zou Huanxin, Luo Tiancheng, Zhang Yue, Zhou Shilin. Combined Conditional Random Fields Model for Supervised PolSAR Images Classification[J]. Journal of Radars, 2017, 6(5): 541-553. doi: 10.12000/JR16109
Citation: Zou Huanxin, Luo Tiancheng, Zhang Yue, Zhou Shilin. Combined Conditional Random Fields Model for Supervised PolSAR Images Classification[J]. Journal of Radars, 2017, 6(5): 541-553. doi: 10.12000/JR16109

基于组合条件随机场的极化SAR图像监督地物分类

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

    邹焕新(1973–),男,广东人,现任国防科技大学电子科学与工程学院副教授,硕士生导师,主要研究方向为SAR图像解译、多源遥感信息融合等。E-mail: hxzou2008@163.com

    罗天成(1991–),男,湖南人,现为国防科技大学电子科学与工程学院硕士研究生,主要研究方向为极化SAR图像地物分类。E-mail: 15616045932@163.com

    张   月(1990–),女,河南人,现为国防科技大学电子科学与工程学院硕士研究生,主要研究方向为极化SAR图像地物分类。E-mail: yue1554415@163.com

    周石琳(1965–),男,湖南人,现任国防科技大学电子科学与工程学院教授,博士生导师,主要研究方向为计算机视觉与智能信息处理、多源遥感信息融合等。E-mail: slzhoumail@163.com

    通讯作者:

    邹焕新   hxzou2008@163.com

  • 中图分类号: TN957

Combined Conditional Random Fields Model for Supervised PolSAR Images Classification

Funds: The National Natural Science Foundation of China (61331015)
  • 摘要: 在复杂极化SAR (Synthetic Aperture Radar, SAR)图像的地物分类中,可以通过提取多种特征及利用上下文信息来提高分类精度。特征维度的增加会产生过拟合和特征干扰,从而导致分类器性能降低。针对这个问题,该文提出了一种基于组合条件随机场(Conditional Random Fields, CRF)的极化SAR图像监督地物分类方法。不同于CRF模型通过将多个特征矢量堆叠形成一个新的特征矢量以利用多种特征信息的传统方式,组合CRF模型首先将不同类型特征分成多个特征子集合分别训练同一个CRF模型得到多个有差异的子分类器,从而得到多个地物分类结果,然后以单个子分类器的归一化总体分类精度作为加权系数将多个分类结果进行融合以得到最终的分类结果。两组真实极化SAR图像分类实验结果表明,该文方法的分类精度比利用单个特征矢量单个子分类器的分类精度有明显提升。对于实验采用的两组数据,该文方法的分类精度比利用多个特征矢量堆叠的分类精度分别提高13.38%和11.55%,同时也比基于SVM (Support Vector Machines, SVM)的分类精度分别提高13.78%和14.75%。
  • 图  1  基于组合CRF模型的极化SAR图像分类流程图

    Figure  1.  Flowchart of the combined CRF model for PolSAR images classification

    图  2  本文采用的两组实测极化SAR图像数据

    Figure  2.  Two real-world PolSAR images for experiments

    图  3  基于堆叠特征的CRF分类器分类结果

    Figure  3.  Classification result of CRF model based on stacking features

    图  4  基于不同特征的子分类器分类结果

    Figure  4.  Classification results based on child classifiers using different features

    图  5  组合CRF分类结果

    Figure  5.  Classification result of combined CRF model

    图  6  基于l-1特征选择的CRF分类器分类结果

    Figure  6.  Classification result of CRF model based on l-1 features selection

    图  7  SVM分类器分类结果

    Figure  7.  Classification result of SVM model

    图  8  基于投票方法的组合CRF分类结果

    Figure  8.  Classification result of combined CRF model based on majority voting

    图  9  基于堆叠特征的CRF分类结果

    Figure  9.  Classification result of CRF model based on stacking features

    图  10  子分类器分类结果

    Figure  10.  Classification results based on different child classifiers

    图  11  组合CRF分类器分类结果

    Figure  11.  Classification result of combined CRF model

    图  12  基于l-1特征选择分类器分类结果

    Figure  12.  Classification result of CRF model based on l-1 features selection

    图  13  SVM分类器分类结果

    Figure  13.  Classification result of SVM model

    图  14  基于投票表决法的组合CRF分类结果

    Figure  14.  Classification result of CRF model based on majority voting

    表  1  本文方法使用的特征

    Table  1.   Features used in the proposed method

    特征 描述 特征内容 数量
    基于测量和目标分解的特征 协方差矩阵上三角 |C11|, |C12|, |C13|, |C22|, |C23|, |C33| 6
    VV/HH比值[25] 10log(|SVV|2$\big/$|SHH|2) 1
    HV/HH比值[25] 10log(|SHV|2$\big/$|SHH|2) 1
    HV/VV比值[25] 10log(|SHV|2$\big/$|SVV|2) 1
    Yamaguchi 四分量分解 表面散射(Ps)、二次散射(Pd)、体散射(Pv)、螺旋体散射(Pc) 4
    ${\rm{H}}/{\rm{\alpha }}/{\rm A}$分解 熵(H)、平均散射机制( ${\rm{\alpha }}$)、各向异性系数(A) 3
    Freeman分解 表面散射(Ps)、二次散射(Pd)、体散射(Pv) 3
    纹理特征 基于GLCM的纹理特征 相关性、对比度、熵、同质性 4
    形态学特征 腐蚀和膨胀[21] 重建开操作、重建闭操作 8
    下载: 导出CSV

    表  2  堆叠特征分类精度

    Table  2.   Classification accuracy of stacking features

    类别 豆类 森林 马铃薯 苜蓿 小麦 裸土 甜菜 油菜籽 豌豆 草地 水体 用户精度
    豆类 6486 66 30 0 0 0 48 0 153 66 0 0.9470
    森林 1210 15171 1333 46 74 11 284 11 343 4 62 0.8179
    马铃薯 340 220 11831 15 1292 44 1521 14 2736 92 65 0.6378
    苜蓿 523 42 110 7743 84 164 69 39 5 1266 476 0.7360
    小麦 13 63 94 258 15723 21 656 76 51 26 1073 0.8709
    裸土 6 2 5 2 0 3266 0 24 1 0 0 0.9879
    甜菜 33 0 2716 53 27 22 7309 30 42 172 6 0.7021
    油菜籽 5 219 63 25 1095 1757 404 10013 126 2 284 0.7156
    豌豆 3 117 888 12 332 26 37 45 8630 5 12 0.8539
    草地 392 12 808 481 48 2703 415 151 49 2312 100 0.3092
    水体 18 57 429 18 207 3449 30 4 0 2 6192 0.5950
    制图精度 0.9035 0.9500 0.6463 0.8948 0.8327 0.2859 0.6785 0.9621 0.7111 0.5858 0.7487
    总体分类精度:74.06%;Kappa系数:0.7113
    下载: 导出CSV

    表  3  组合CRF分类精度

    Table  3.   Classification accuracy of combined CRF model

    类别 豆类 森林 马铃薯 苜蓿 小麦 裸土 甜菜 油菜籽 豌豆 草地 水体 用户精度
    豆类 6467 93 36 60 12 0 64 0 1 116 0 0.9442
    森林 97 17095 1234 97 18 0 0 2 3 3 2 0.9216
    马铃薯 37 534 17155 41 14 2 180 46 74 86 1 0.9441
    苜蓿 263 1 15 8482 460 15 5 45 1 1234 0 0.8062
    小麦 14 0 8 0 16807 2 2 1145 76 0 0 0.9309
    裸土 0 0 0 0 0 3173 0 17 0 0 116 0.9598
    甜菜 14 11 857 8 102 9 8598 400 141 270 0 0.8259
    油菜籽 0 0 8 5 4134 658 2 8748 77 361 0 0.6252
    豌豆 0 0 131 3 289 4 0 283 9388 8 1 0.9289
    草地 172 4 138 568 73 255 1 33 0 6263 0 0.8376
    水体 4 0 0 1 1 1152 0 0 1 0 9247 0.8886
    制图精度 0.9150 0.9638 0.8761 0.9155 0.7671 0.6055 0.9713 0.8161 0.9617 0.7509 0.9874
    总体分类精度:87.44%;Kappa系数:0.8596
    下载: 导出CSV

    表  4  堆叠特征分类精度

    Table  4.   Classification accuracy of stacking features

    类别 道路 林地 城郊 农田 其他 用户精度
    道路 122208 1456 6531 9828 18352 0.7435
    林地 1913 91840 50760 106297 8531 0.3541
    城郊 25348 9464 204174 16170 40917 0.6896
    农田 2552 5649 6351 64194 30475 0.5880
    其他 109049 150 10698 24033 473749 0.7670
    制图精度 0.4877 0.8539 0.7384 0.2925 0.8282
    总体分类精度:66.78%;Kappa系数:0.5547
    下载: 导出CSV

    表  5  组合CRF分类精度

    Table  5.   Classification accuracy of combined CRF model

    类别 道路 林地 城郊 农田 其他 用户精度
    道路 93169 1970 7725 3805 57661 0.5668
    林地 434 201255 50737 4032 2883 0.7760
    城郊 955 10321 267059 9783 7955 0.9020
    农田 407 12851 18465 56485 20968 0.5174
    其他 22984 2355 23078 54041 515211 0.8341
    制图精度 0.7899 0.8798 0.7276 0.4406 0.8520
    总体分类精度:78.33%;Kappa系数:0.7013
    下载: 导出CSV
  • [1] Zhang La-mei, Wang Xiao, Sun Liang-jie, et al.. Contextual sparse representation and its application in polarimetric SAR image classification[C]. IET International Radar Conference, Hangzhou, China, 2015: 1–5.
    [2] Shang Fang and Hirose A. Use of Poincare sphere parameters for fast supervised PolSAR land classification[C]. Proceedings of 2013 IEEE International Geoscience and Remote Sensing Symposium, Melbourne, VIC, 2013: 3175–3178.
    [3] Xie Wen, Jiao Li-cheng, and Zhao Jin. PolSAR image classification via D-KSVD and NSCT-Domain features extraction[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(2): 227–231. doi:  10.1109/LGRS.2015.2506718
    [4] Kong J A, Swartz A A, Yueh H A, et al.. Identification of terrain cover using the optimum polarimetric classifier[J]. Journal of Electromagnetic Waves and Applications, 1988, 2(2): 171–194.
    [5] Silva W B, Freitas C C, Sant’Anna S J S, et al.. Classification of segments in PolSAR imagery by minimum stochastic distances between Wishart distributions[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2013, 6(3): 1263–1273. doi:  10.1109/JSTARS.2013.2248132
    [6] 孙即祥. 现代模式识别[M]. 长沙: 国防科技大学出版社, 2002: 385–447.Sun Ji-xiang. Modern Pattern Recognition[M]. Changsha: Press of National University of Defense Technology, 2002: 385–447.
    [7] Lee J S and Ainsworth T L. An overview of recent advances in Polarimetric SAR information extraction: Algorithms and applications[C]. Proceedings of 2010 IEEE International Geoscience and Remote Sensing Symposium, Honolulu, HI, 2010: 851–854.
    [8] 钟平. 面向图像标记的随机场模型研究[D]. [博士论文], 国防科学技术大学, 2008.Zhong Ping. Random fields model for image label[D]. [Ph.D. dissertation], National University of Defense Technology, 2008.
    [9] Zhong Ping and Wang Run-sheng. A multiple conditional random fields ensemble model for urban area detection in remote sensing optical images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(12): 3978–3988. doi:  10.1109/TGRS.2007.907109
    [10] Lafferty J, McCallum A, and Pereira F. Conditional random fields: Probabilistic models for segmenting and labeling sequence data[C]. Proceedings of the Eighteenth International Conference on Machine Learning, San Francisco, CA, USA, 2001: 282–289.
    [11] Kumar S and Hebert M. Discriminative random fields: A discriminative framework for contextual interaction in classification[C]. IEEE International Conference on Computer Vision, Piscataway, NJ, 2003: 1150–1157.
    [12] 吴立珍. 面向UAV战场感知的目标特征建模与应用研究[D]. [博士论文], 国防科学技术大学, 2012.Wu Li-zhen. Research on object feature modeling and applications for battlefield awareness of unmanned aerial vehicle[D]. [Ph.D. dissertation], National University of Defense Technology, 2012.
    [13] Pearl J. Probabilistic Reasoning in Intelligent Systems: Networks of Plausible Inference[M]. California: Morgan Kaufmann, 1988: 247–289.
    [14] Wang Lei-guang, Dai Qin-ling, and Huang Xin. Spatial regularization of pixel-based classification maps by a two-step MRF method[C]. Proceedings of IEEE International Geoscience and Remote Sensing Symposium, Beijing, 2016: 2407–2410.
    [15] Zhong Yan-fei, Zhao Ji, and Zhang Liang-pei. A hybrid object-oriented conditional random field classification framework for high spatial resolution remote sensing imagery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(11): 7023–7037. doi:  10.1109/TGRS.2014.2306692
    [16] Pieczynski W and Tebbache A N. Pairwise markov random fields and segmentation of textured images[J]. Machine Graphics and Vision, 2000, 9(3): 705–718.
    [17] Kumar S and Hebert M. A hierarchical field framework for unified context-based classification[C]. Proceedings of the Tenth IEEE International Conference on Computer Vision, Beijing, 2005: 1284–1291.
    [18] Jiang Wei, Chang S F, and Loui A C. Context-based concept fusion with boosted conditional random fields[C]. Proceedings of IEEE International Conference on Acoustics, Speech and Signal Processing, Honolulu, HI, 2007: I-949–I-952.
    [19] Lee C H, Schmidt M, Murtha A, et al.. Segmenting brain tumors with conditional random fields and support vector machines[C]. Proceedings of the First International Conference on Computer Vision for Biomedical Image Applications, Beijing, China, 2005: 469–478.
    [20] Do T M T and Artières T. Polynomial conditional random fields for signal processing[C]. Proceedings of the 2006 Conference on ECAI 2006: 17th European Conference on Artificial Intelligence, Riva del Garda, Italy, 2006: 797–798.
    [21] Du Pei-jun, Samat A, Waske B, et al.. Random forest and rotation forest for fully polarized SAR image classification using polarimetric and spatial features[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2015, 105: 38–53. doi:  10.1016/j.isprsjprs.2015.03.002
    [22] Zou Tong-yuan, Yang Wen, Dai Deng-xin, et al.. Polarimetric SAR image classification using multifeatures combination and extremely randomized clustering forests[J]. EURASIP Journal on Advances in Signal Processing, 2010, 2010: Article No. 4.
    [23] 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
    [24] Haralick R M, Shanmugam K, and Dinstein I. Textural features for image classification[J]. IEEE Transactions on Systems, Man, and Cybernetics, 1973, SMC-3(6): 610-621. doi:  10.1109/TSMC.1973.4309314
    [25] 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
  • [1] 陈诗强, 洪文.  基于极化轴比参数的圆极化波发射误差分析方法 . 雷达学报, 2020, 9(2): 343-353. doi: 10.12000/JR19063
    [2] 庞雷, 张风丽, 王国军, 刘娜, 邵芸, 张家萌, 赵钰川, 庞蕾.  建筑物Ku波段极化SAR成像仿真及损毁评估特征分析 . 雷达学报, 2020, 9(3): 578-587. doi: 10.12000/JR20061
    [3] 朱庆涛, 殷君君, 曾亮, 杨健.  基于邻域一致性的极化SAR图像仿射配准 . 雷达学报, 2020, 9(): 1-12. doi: 10.12000/JR20120
    [4] 秦先祥, 余旺盛, 王鹏, 陈天平, 邹焕新.  基于复值卷积神经网络样本精选的极化SAR图像弱监督分类方法 . 雷达学报, 2020, 9(3): 525-538. doi: 10.12000/JR20062
    [5] 邹焕新, 李美霖, 马倩, 孙嘉赤, 曹旭, 秦先祥.  一种基于张量积扩散的非监督极化SAR图像地物分类方法 . 雷达学报, 2019, 8(4): 436-447. doi: 10.12000/JR19057
    [6] 滑文强, 王爽, 郭岩河, 谢雯.  基于邻域最小生成树的半监督极化SAR图像分类方法 . 雷达学报, 2019, 8(4): 458-470. doi: 10.12000/JR18104
    [7] 张腊梅, 张思雨, 董洪伟, 朱厦.  基于Pinball损失函数支持向量机的极化SAR图像鲁棒分类 . 雷达学报, 2019, 8(4): 448-457. doi: 10.12000/JR19055
    [8] 张向荣, 于心源, 唐旭, 侯彪, 焦李成.  基于马尔科夫判别谱聚类的极化SAR图像分类方法 . 雷达学报, 2019, 8(4): 425-435. doi: 10.12000/JR19059
    [9] 胡涛, 李卫华, 秦先祥, 王鹏, 余旺盛, 李军.  基于深度卷积神经网络和条件随机场模型的PolSAR图像地物分类方法 . 雷达学报, 2019, 8(4): 471-478. doi: 10.12000/JR18065
    [10] 陶臣嵩, 陈思伟, 李永祯, 肖顺平.  结合旋转域极化特征的极化SAR地物分类 . 雷达学报, 2017, 6(5): 524-532. doi: 10.12000/JR16131
    [11] 折小强, 仇晓兰, 雷斌, 张薇, 卢晓军.  一种极化熵结合混合GEV模型的全极化SAR潮间带区域地物分类方法 . 雷达学报, 2017, 6(5): 554-563. doi: 10.12000/JR16149
    [12] 钟能, 杨文, 杨祥立, 郭威.  基于混合Wishart模型的极化SAR图像非监督分类 . 雷达学报, 2017, 6(5): 533-540. doi: 10.12000/JR16133
    [13] 杨文, 钟能, 严天恒, 杨祥立.  基于黎曼流形的极化SAR图像分类 . 雷达学报, 2017, 6(5): 433-441. doi: 10.12000/JR17031
    [14] 计科峰, 王海波, 冷祥光, 邢相薇, 康利鸿.  星载简缩极化SAR船舶目标检测技术研究 . 雷达学报, 2016, 5(6): 607-619. doi: 10.12000/JR16083
    [15] 张杰, 张晰, 范陈清, 孟俊敏.  极化SAR在海洋探测中的应用与探讨 . 雷达学报, 2016, 5(6): 596-606. doi: 10.12000/JR16124
    [16] 洪文.  基于混合极化架构的极化SAR:原理与应用(英文) . 雷达学报, 2016, 5(6): 559-595. doi: 10.12000/JR16074
    [17] 邢艳肖, 张毅, 李宁, 王宇, 胡桂香.  一种联合特征值信息的全极化SAR图像监督分类方法 . 雷达学报, 2016, 5(2): 217-227. doi: 10.12000/JR16019
    [18] 孙勋, 黄平平, 涂尚坦, 杨祥立.  利用多特征融合和集成学习的极化SAR图像分类 . 雷达学报, 2016, 5(6): 692-700. doi: 10.12000/JR15132
    [19] 余凡, 袁捷.  一种改进的随机粗糙地表的双尺度微波散射模型研究(英文) . 雷达学报, 2015, 4(5): 560-570. doi: 10.12000/JR15067
    [20] 滑文强, 王爽, 侯彪.  基于半监督学习的SVM-Wishart极化SAR图像分类方法 . 雷达学报, 2015, 4(1): 93-98. doi: 10.12000/JR14138
  • 加载中
图(14) / 表 (5)
计量
  • 文章访问数:  1132
  • HTML全文浏览量:  343
  • PDF下载量:  482
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-09-22
  • 修回日期:  2016-12-14
  • 网络出版日期:  2017-07-31
  • 刊出日期:  2017-10-28

基于组合条件随机场的极化SAR图像监督地物分类

doi: 10.12000/JR16109
    基金项目:  国家自然科学基金(61331015)
    作者简介:

    邹焕新(1973–),男,广东人,现任国防科技大学电子科学与工程学院副教授,硕士生导师,主要研究方向为SAR图像解译、多源遥感信息融合等。E-mail: hxzou2008@163.com

    罗天成(1991–),男,湖南人,现为国防科技大学电子科学与工程学院硕士研究生,主要研究方向为极化SAR图像地物分类。E-mail: 15616045932@163.com

    张   月(1990–),女,河南人,现为国防科技大学电子科学与工程学院硕士研究生,主要研究方向为极化SAR图像地物分类。E-mail: yue1554415@163.com

    周石琳(1965–),男,湖南人,现任国防科技大学电子科学与工程学院教授,博士生导师,主要研究方向为计算机视觉与智能信息处理、多源遥感信息融合等。E-mail: slzhoumail@163.com

    通讯作者: 邹焕新   hxzou2008@163.com
  • 中图分类号: TN957

摘要: 在复杂极化SAR (Synthetic Aperture Radar, SAR)图像的地物分类中,可以通过提取多种特征及利用上下文信息来提高分类精度。特征维度的增加会产生过拟合和特征干扰,从而导致分类器性能降低。针对这个问题,该文提出了一种基于组合条件随机场(Conditional Random Fields, CRF)的极化SAR图像监督地物分类方法。不同于CRF模型通过将多个特征矢量堆叠形成一个新的特征矢量以利用多种特征信息的传统方式,组合CRF模型首先将不同类型特征分成多个特征子集合分别训练同一个CRF模型得到多个有差异的子分类器,从而得到多个地物分类结果,然后以单个子分类器的归一化总体分类精度作为加权系数将多个分类结果进行融合以得到最终的分类结果。两组真实极化SAR图像分类实验结果表明,该文方法的分类精度比利用单个特征矢量单个子分类器的分类精度有明显提升。对于实验采用的两组数据,该文方法的分类精度比利用多个特征矢量堆叠的分类精度分别提高13.38%和11.55%,同时也比基于SVM (Support Vector Machines, SVM)的分类精度分别提高13.78%和14.75%。

English Abstract

邹焕新, 罗天成, 张月, 周石琳. 基于组合条件随机场的极化SAR图像监督地物分类[J]. 雷达学报, 2017, 6(5): 541-553. doi: 10.12000/JR16109
引用本文: 邹焕新, 罗天成, 张月, 周石琳. 基于组合条件随机场的极化SAR图像监督地物分类[J]. 雷达学报, 2017, 6(5): 541-553. doi: 10.12000/JR16109
Zou Huanxin, Luo Tiancheng, Zhang Yue, Zhou Shilin. Combined Conditional Random Fields Model for Supervised PolSAR Images Classification[J]. Journal of Radars, 2017, 6(5): 541-553. doi: 10.12000/JR16109
Citation: Zou Huanxin, Luo Tiancheng, Zhang Yue, Zhou Shilin. Combined Conditional Random Fields Model for Supervised PolSAR Images Classification[J]. Journal of Radars, 2017, 6(5): 541-553. doi: 10.12000/JR16109
    • 极化SAR (Synthetic Aperture Radar, SAR)图像地物分类在地质勘探、地形分析以及灾害监测等方面具有广泛的应用[1]。普通SAR图像由于含有较少的信息,使地物分类精度较低,而多极化SAR图像包含多种极化方式使图像能够获得目标的全部极化散射特性,从而能够取得较好的分类结果[2,3]。根据分类过程中是否需要加入人为因素的干预,可以将分类过程分为监督分类和无监督分类两大类。一般来说,无监督分类由于缺少先验知识的引导,因此其分类精度相对较低;而监督分类在选取的训练样本具有足够代表性的前提下,可以达到很高的分类精度。极化SAR图像地物分类方法在1988年由麻省理工学院的Jin-Au Kong领导的研究组首次提出,这是极化SAR图像Bayes有监督分类的最初形式[3,4]。经过20多年的发展,国内外众多学者针对极化SAR数据有监督地物分类提出了很多方法,包括线性(或非线性)判别函数法、支持矢量机、神经网络、Bayes决策等[57]

      近年来,随着获取的极化SAR图像越来越复杂,对地物分类提出了更多的挑战。例如,Bayes决策分类方法需要对极化SAR图像数据进行复杂的统计建模;传统的提取一种或少数几种特征的方式难以表达图像中的目标包含的本质信息[8,9]。因此,为了得到更高的极化SAR图像地物分类精度,对原有的分类模型进行改进或者提出新的分类模型势在必行。

      Lafferty教授[10]在2001年提出的一种最初应用于文本序列分割和标记的CRF (Conditional Random Fields, CRF)模型引起了分类领域学者的高度重视。Kumar和Hebert[11]将CRF模型扩展到2维领域,使CRF模型可以使用不同的势函数以满足机器视觉领域的各种具体应用要求,极大地增加了CRF模型的灵活性。CRF模型是一种对标记图像和观测图像直接建模的后验概率模型,可以避免因统计模型不准确而造成的分类精度下降问题。CRF模型不仅能够利用上下文信息,而且有很强的融合多特征的能力,因此对地物类型复杂的极化SAR图像有良好的分类能力。CRF模型对从图像中提取的多种特征的传统利用方法是将任一像素点的全部特征首尾相接而堆叠成一个特征矢量。理论上,提取的极化SAR图像特征越多,分类器的性能越好,得到的极化SAR图像分类精度越高[8]。但是,特征维度的增加将导致以下两方面的问题:一是多种特征之间存在干扰使得模型在训练时难以获取某些重要特征所表征的目标本质信息[8];二是同其它分类器一样,过多的特征可能会导致在训练分类模型时产生过拟合的问题,从而使模型的泛化能力下降。解决上述问题的第1种技术途径是特征稀疏化处理,包括PCA[8]l-1正则化法[8]等特征选择的方法;另一种途径是多分类器组合方法。

      本文针对基于堆叠特征的CRF模型容易引起特征干扰和过拟合的问题,提出了一种组合CRF模型用于极化SAR图像的地物分类问题。组合CRF模型首先将不同类型的特征分成多个特征子集合分别训练同一个CRF模型得到多个有差异的子分类器,从而得到多个地物分类结果,然后以单个子分类器的分类精度作为加权系数将多个分类结果进行融合以得到最终的分类结果。两组真实数据的实验结果表明,组合CRF模型相比于传统的特征堆叠利用多特征的方式,其分类精度有明显提高。

    • 概率图模型是图论和概率论相结合而发展起来的产物,能够描述机器视觉、图像处理领域日益复杂的数据关系[12]。概率图模型根据变量之间是否存在因果性的依赖关系,可以分为:有向图模型,如Bayes网络[13]等;无向图模型,如马尔科夫随机场[14](Markov Random Field, MRF)和CRF[15]等随机场模型。随机场模型的应用始于对MRF的研究,经过多年的发展,MRF已经成功地应用于图像分析的各个领域。在研究的过程中,MRF的一些劣势限制了其发展,比如MRF需要对输入特征的似然分布进行详细的建模以及对观测到的图像特征集作条件独立假设,而CRF是对后验概率直接建模,不存在这种限制,从而使其在极化SAR图像分类领域获得了广泛的应用[16]

    • Lafferty教授[10]对条件随机场进行了详细的定义:对于随机场xy,当给定随机场y时,x的分布满足马尔科夫性质[8]

      $$p({x_i}|{{y}},{x_j},j \ne i,j \in S) = p({x_i}|{{y}},{x_j},j \in {\eta _i})$$ (1)

      其中, ${\eta _i}$为位置i的邻域系统,S表示一个定义有邻域系统的有限集合,那么此时(x, y)称为一个条件随机场[8]

      x是关于特征y的标记,且 $C = \{ \{ {{{y}}_c},{x_c}\} \} $是邻域系统上的基团集,则关于标记x的CRF模型为:

      $$p(x|{{y}},\theta ) = \frac{1}{{Z\left( {{{y}},\theta } \right)}}\exp \left\{ {\sum\limits_{c \ \in C} {{\phi _c}\left( {{x_c},{{y}},\theta } \right)} } \right\}$$ (2)

      其中, $Z({{y}},\theta ) =\displaystyle \sum\nolimits_x {\exp \left\{\displaystyle \sum\nolimits_{c \ \in C} {{\phi _c}({x_c},{{y}},\theta )} \right\} } $为归一化函数; ${\phi _c}{\rm{(}}{x_c}{\rm{,}}{{y}}{\rm{,}}\theta {\rm{)}}$是定义在基团c上带有参数 $\theta $的势函数。

      假设用 $x = {({x_s})_{s \in S}}$表示标记图像, ${{y}} = {({{{y}}_s})_{s \in S}}$表示观测图像。当只考虑单位置基团和双位置基团时,式(2)给出的条件随机场数学表达式可以进一步表示为:

      $$p\left( {x\left| {\mathit{{y}},\theta } \right.} \right) = \frac{1}{{Z\left( {\mathit{{y}},\theta } \right)}}\exp \left\{ {\sum\limits_{i \in S} {{\phi _i}\left( {{x_i},\mathit{{y,w}}} \right)} + \sum\limits_{i \in S} {\sum\limits_{j \in {\eta _i}} {{\phi _{ij}}\left( {{x_i},{x_j},\mathit{{y,v}}} \right)} } } \right\}$$ (3)

      其中, ${\phi _i}( \cdot )$和 ${\phi _{ij}}( \cdot )$分别表示具有参数w和参数v的单位置势函数和双位置势函数; $\theta $表示参数集 $\theta = \{ {{w}},{{v}}\} $。

    • 针对不同的应用背景,需要采用不同的势函数。常见的势函数包括多元对数回归(Multinomial Logistic Regression, MLR)[17]、Boost[18]、支持向量机(Support Vector Machines, SVM)[19]、多项式函数[20]等。MLR分类器不需要观测数据满足正态分布等限制条件,可以适应更多的观测数据,因此本文选择MLR分类器作为单位置势函数[8]

      $${\phi _i}({x_i},{{y}},{{w}}) = \sum\limits_{l = 1}^L {\delta ({x_i} = l)} \log p({x_i} = l|{{y}},{{w}})$$ (4)

      其中, $\delta ( \cdot )$表示Kronecker delta函数,且:

      $$p({x_i} = l|{y},{w}) = \left\{ {\begin{array}{*{20}{c}}{\frac{{\exp \left({w}_l^{\rm{T}}{{y}_i}\right)}}{{1 +\displaystyle \sum\limits_{k = 1}^{L{\rm{ - 1}}} {\exp \left({w}_k^{\rm{T}}{{y}_i}\right)} }}{\rm{,\; if }}\; l < L}\\{\frac{1}{{1 +\displaystyle \sum\limits_{k = 1}^{L - 1} {{\rm{exp}}\left({w}_k^{\rm{T}}{{y}_i}\right)} }}{\rm{,\; if }}\; l = L}\end{array}} \right.$$ (5)

      其中, $l \in \{ 1,2, ·\!·· ,L\} $且L表示图像的类别数目;wk为对应第k类的参数向量 ${\left[ {{w_{k1}},{w_{k2}}, ·\!·\!· ,{w_{kD}}} \right]^{\rm{T}}}$; w表示所有的wk( $k = 1,2, ·\!·\!· ,L - 1$)首尾相接得到的参数向量,上标T表示转置。

      双位置势函数可以表达图像空间的上下文信息,本文选择如式(6)所示的广义Ising/Potts模型作为双位置势函数模型[8]

      $$\begin{aligned}{\phi _{ij}}({x_i},{x_j},{{y}},{{v}}) = & \sum\limits_{k,l \in \{ 1, 2, \cdots ,L\} } {{{{v}}_{kl}^{\rm{T}}}} \\&\cdot {{{g}}_{ij}}({{y}}) \delta ({x_i} = k)\delta ({x_j} = l)\end{aligned}$$ (6)

      其中,vkl为参数向量,v表示所有的vkl( $k,l = 1,$ $2, ·\!·\!· ,L$)首尾相接得到的参数向量, ${{{g}}_{ij}}\left( \cdot \right)$表示从整个观测数据提取的位置对(i, j)的特征向量, $\delta ( \cdot )$表示Kronecker delta函数,上标T表示转置。

    • CRF分类器利用多特征的方式为:把从图像中任一位置提取得到的多种特征首尾相接而堆叠形成特征矢量,并以这个特征矢量作为分类器的输入矢量。这种利用多特征的方式,不仅操作简单,而且理论上采用的特征越多,得到的分类精度应该越高,但实际中这种方式容易引起特征干扰和过拟合问题,从而影响最终的分类精度。为了更好地利用多种特征信息,同时避免出现特征干扰和过拟合问题,本文提出采用组合CRF模型来提高极化SAR图像的地物分类精度。

    • 组合模型是利用多个子分类器组合分类来获得比单个分类器更高的分类精度。组合模型结合决策树的思想在机器学习领域应用非常广泛,这种思想通常用来解决单个决策树判别能力不足以及容易产生过拟合的问题[21]。本文根据这种思想提出一种组合CRF模型。组合CRF模型主要分为两部分:(1)利用特定的策略得到多个有差异性的子分类器;(2)使用一种规则将多个子分类器融合。为得到有差异性的子分类器,本文将提取的特征分成多个子特征矢量分别训练同一个CRF模型。对于组合子分类器的规则,研究工作者给出了大量的组合方法,包括投票表决法[21]、Boosting方法[8]等。其中投票表决法是较早提出并取得广泛应用的方法,成为比较和评估新提出的融合规则的基准。为得到组合CRF模型,首先提取训练样本的K个子特征矢量,然后利用K个子特征矢量训练得到K个有差异的子分类器,最后使用融合公式将K个子分类器融合从而得到最终的组合CRF模型。

    • 本文首先提取了4种类型的特征,包括基于功率测量的特征[21]、基于目标分解的特征[22,23]、基于形态学的特征[21]、基于灰度共生矩阵(Gray Level Co-occurrence Matrix, GLCM)的纹理特征[24],然后将这些特征分成3组子特征矢量。表1为划分后的3组子特征矢量,分别用 $\{ f_i^k({{y}}),k = 1,2,3\} $表示。3组子特征矢量将作为输入用于训练CRF模型得到3个子CRF分类器。

      表 1  本文方法使用的特征

      Table 1.  Features used in the proposed method

      特征 描述 特征内容 数量
      基于测量和目标分解的特征 协方差矩阵上三角 |C11|, |C12|, |C13|, |C22|, |C23|, |C33| 6
      VV/HH比值[25] 10log(|SVV|2$\big/$|SHH|2) 1
      HV/HH比值[25] 10log(|SHV|2$\big/$|SHH|2) 1
      HV/VV比值[25] 10log(|SHV|2$\big/$|SVV|2) 1
      Yamaguchi 四分量分解 表面散射(Ps)、二次散射(Pd)、体散射(Pv)、螺旋体散射(Pc) 4
      ${\rm{H}}/{\rm{\alpha }}/{\rm A}$分解 熵(H)、平均散射机制( ${\rm{\alpha }}$)、各向异性系数(A) 3
      Freeman分解 表面散射(Ps)、二次散射(Pd)、体散射(Pv) 3
      纹理特征 基于GLCM的纹理特征 相关性、对比度、熵、同质性 4
      形态学特征 腐蚀和膨胀[21] 重建开操作、重建闭操作 8
    • 对极化SAR图像的任一位置i,提取的K种类型的子特征矢量可以表示为 $f_i^{1:K}({{y}}) = \{ f_i^1({{y}}),$ $f_i^2({{y}}), ·\!·\!· ,f_i^K({{y}})\} $,其中 $f_i^k({{y}}) = [f_{i1}^k({{y}}),f_{i2}^k({{y}}), ·\!·\!· ,$ $f_{iD}^k({{y}})]{\rm{, }}1 \le k \le K$表示图像位置i处的第k种类型特征提取到的D维特征,则第k种类型特征表示为 ${f^k}({{y}}) = \{ f_i^k({{y}}),i \in S\} $,其中S表示一个定义有邻域系统的有限集合。用 $f^{1:K}({{y}}) = \{ f^1({{y}}),f^2({{y}}),·\!·\!· , $ $f^K({{y}})\} $表示输入图像的特征y,首先将K个子特征矢量分别训练同一个CRF模型得到K个有差异性的子CRF分类器,然后使用式(7)的组合公式将子分类器融合。

      $$p\left(x|{f^{1:K}}({{y}})\right) = \sum\limits_{k = 1}^K {{\alpha _k}p\left(x|{f^k}({{y}})\right)} $$ (7)

      其中, $p\left(x|{f^k}({{y}})\right)$表示第k个子CRF分类器的分类结果; ${\alpha _k}$表示第k个子CRF分类器在组合公式中的自适应权值系数,取值范围为 ${\alpha _k} \in \left[ {0,1} \right]$。在本文的实验中,为了简便地得到每个子分类器对应的自适应权值系数,利用选定的测试样本测试子分类器并得到对应的分类精度,最后将所有子分类器得到的分类精度归一化得到 ${\alpha _k}(k = 1,2, ·\!·\!· ,K)$。

      为了清楚地说明本文的算法过程,图1给出了本文基于组合CRF模型的极化SAR图像分类流程图。

      图  1  基于组合CRF模型的极化SAR图像分类流程图

      Figure 1.  Flowchart of the combined CRF model for PolSAR images classification

    • 为了对本文提出算法的性能进行分析和评估,我们选择了两组实测极化SAR图像数据进行实验。图2(a)是由NASA/JPL提供的尺寸为750×1024像素的荷兰Flevoland地区极化AIRSAR图像数据。图2(c)是由German Aerospace Center提供的尺寸为1300×1200像素的德国Oberpfaffenhofen地区的E-SAR图像数据。图2(a)图2(c)中的红色方框表示本文方法在这两组实测数据中选择的训练样本,图2(b)图2(d)分别为两幅极化SAR图像对应的真值地物图,两幅真值地物图分别标记有11种和5种地物类别。

      图  2  本文采用的两组实测极化SAR图像数据

      Figure 2.  Two real-world PolSAR images for experiments

    • 本文首先将从Flevoland地区极化SAR图像中提取的3组特征矢量以首尾相接的堆叠方式组成共31维的特征矢量用于训练单个CRF模型得到CRF分类器。图3(a)图3(b)分别为基于堆叠特征的CRF分类器的分类结果以及分类结果与真实地物图相减得到的错分地物图,表2为分类结果的混淆矩阵。

      表 2  堆叠特征分类精度

      Table 2.  Classification accuracy of stacking features

      类别 豆类 森林 马铃薯 苜蓿 小麦 裸土 甜菜 油菜籽 豌豆 草地 水体 用户精度
      豆类 6486 66 30 0 0 0 48 0 153 66 0 0.9470
      森林 1210 15171 1333 46 74 11 284 11 343 4 62 0.8179
      马铃薯 340 220 11831 15 1292 44 1521 14 2736 92 65 0.6378
      苜蓿 523 42 110 7743 84 164 69 39 5 1266 476 0.7360
      小麦 13 63 94 258 15723 21 656 76 51 26 1073 0.8709
      裸土 6 2 5 2 0 3266 0 24 1 0 0 0.9879
      甜菜 33 0 2716 53 27 22 7309 30 42 172 6 0.7021
      油菜籽 5 219 63 25 1095 1757 404 10013 126 2 284 0.7156
      豌豆 3 117 888 12 332 26 37 45 8630 5 12 0.8539
      草地 392 12 808 481 48 2703 415 151 49 2312 100 0.3092
      水体 18 57 429 18 207 3449 30 4 0 2 6192 0.5950
      制图精度 0.9035 0.9500 0.6463 0.8948 0.8327 0.2859 0.6785 0.9621 0.7111 0.5858 0.7487
      总体分类精度:74.06%;Kappa系数:0.7113

      图  3  基于堆叠特征的CRF分类器分类结果

      Figure 3.  Classification result of CRF model based on stacking features

      图3(a)图3(b)表2可以看出,基于堆叠特征方式得到的分类结果较差,总体分类精度只有74.06%。尤其是草地和水体两类地物被大量地错分,分类精度分别只有30.92%和59.50%。这说明高维特征矢量虽然可以表征图像更多的信息,但是,由于存在特征之间的相互干扰以及分类器过拟合问题,从而导致分类精度不高。

      图4(a)图4(f)分别为基于测量数据的特征和分解特征、纹理特征、形态学特征3组特征矢量分别训练同一个CRF模型得到的子分类器的分类结果图以及错分地物图。从图4中可以看出,各个子分类器比基于堆叠特征的CRF分类器的分类效果差,其原因为:虽然单个子分类器采用了较少的特征,减少了特征之间的相互干扰,但是较少的特征只包含了较少的地物信息,导致子分类器对地物的分类能力较弱;同时,由于子分类器使用不同的特征组进行训练,因此,对同一地物具有不同的分类能力。如图4(c)中水体基本分类正确,而图4(a)图4(e)中与水体对应的像素点基本被错分。

      图  4  基于不同特征的子分类器分类结果

      Figure 4.  Classification results based on child classifiers using different features

      图5为组合CRF模型的分类结果图和错分地物图,表3为分类结果的混淆矩阵。在本文实验中,通过使用各个子分类器先对图2(a)对应的图像数据进行分类得到3个总体分类精度,然后对3个分类精度归一化得到最终的 ${\alpha _i}\left( {i = 1,2,3} \right)$,取值分别为0.4877, 0.2040和0.3083。对比观察图3图4图5的分类结果,可以看出组合CRF模型的分类结果有很大的改善。组合CRF分类器的分类结果优于基于堆叠特征的CRF分类器和各个子分类器的分类结果。同时,对比表2表3,组合CRF分类器总体分类精度从74.06%提高到87.44%,尤其是草地和水体两种地物在组合CRF分类器中分类精度分别达到了83.76%和88.86%,这也进一步说明了组合CRF分类器的有效性。

      图  5  组合CRF分类结果

      Figure 5.  Classification result of combined CRF model

      表 3  组合CRF分类精度

      Table 3.  Classification accuracy of combined CRF model

      类别 豆类 森林 马铃薯 苜蓿 小麦 裸土 甜菜 油菜籽 豌豆 草地 水体 用户精度
      豆类 6467 93 36 60 12 0 64 0 1 116 0 0.9442
      森林 97 17095 1234 97 18 0 0 2 3 3 2 0.9216
      马铃薯 37 534 17155 41 14 2 180 46 74 86 1 0.9441
      苜蓿 263 1 15 8482 460 15 5 45 1 1234 0 0.8062
      小麦 14 0 8 0 16807 2 2 1145 76 0 0 0.9309
      裸土 0 0 0 0 0 3173 0 17 0 0 116 0.9598
      甜菜 14 11 857 8 102 9 8598 400 141 270 0 0.8259
      油菜籽 0 0 8 5 4134 658 2 8748 77 361 0 0.6252
      豌豆 0 0 131 3 289 4 0 283 9388 8 1 0.9289
      草地 172 4 138 568 73 255 1 33 0 6263 0 0.8376
      水体 4 0 0 1 1 1152 0 0 1 0 9247 0.8886
      制图精度 0.9150 0.9638 0.8761 0.9155 0.7671 0.6055 0.9713 0.8161 0.9617 0.7509 0.9874
      总体分类精度:87.44%;Kappa系数:0.8596

      上述实验结果表明,由于组合CRF分类器解决了特征干扰以及过拟合问题,模型在训练的时候可以充分地获取图像特征表达的本质信息,而且组合分类器能够综合利用各个子分类器的优势,因此,组合CRF分类器可以获得比基于堆叠特征的分类器和各个子分类器更高的分类精度。

      由本文引言中的介绍可以知道,特征选择(也称为特征稀疏化处理)也可以用来解决特征干扰和过拟合问题,其中l-1正则化法是一种广泛使用并且效果较好的方法[8]。作为对比,图6给出了基于l-1正则化法的CRF分类器对Flevoland地区极化SAR图像的分类结果以及错分地物图。对比图5图6可以看出,组合CRF分类器与利用l-1正则化法进行特征选择得到的CRF分类器的分类效果基本一致,其中图6的总体分类精度为86.96%(限于篇幅,我们未给出相应的分类精度表),略小于图4的总体分类精度87.44%。这个实验表明,相比于基于l-1正则化法的CRF分类器,组合CRF分类器解决特征干扰和过拟合问题具有更好的效果。

      图  6  基于l-1特征选择的CRF分类器分类结果

      Figure 6.  Classification result of CRF model based on l-1 features selection

      为了更好地评估本文方法的分类性能,我们还利用目前机器学习领域非常流行的SVM分类器对Flevoland地区极化SAR图像进行了分类。观察图7中利用SVM分类器的分类结果(限于篇幅,我们未给出相应的分类精度表)并且与图5的分类结果进行对比,可以发现组合CRF模型的分类性能有较大的优势,图5图7的分类精度高13.78%。这说明组合CRF分类器在对极化SAR图像分类时,通过有效利用上下文信息并且减少特征干扰可以取得更好的分类性能。

      图  7  SVM分类器分类结果

      Figure 7.  Classification result of SVM model

      此外,为了验证本文提出的融合公式的有效性,图8给出了基于投票表决法的组合CRF模型分类结果和错分地物图,对比发现本文的组合CRF分类器模型分类结果优于投票表决法。基于投票表决法的组合CRF分类器取得了83.94%的分类正确率(限于篇幅,我们未给出相应的分类精度表),低于本文提出的组合CRF模型87.44%的总体分类精度,表明本文提出的融合公式对融合子分类器有较好的效果。

      图  8  基于投票方法的组合CRF分类结果

      Figure 8.  Classification result of combined CRF model based on majority voting

    • 为了进一步说明本文方法的适用性,我们利用本文方法对Oberpfaffenhofen地区极化SAR图像数据也开展了分类实验。图9(a)图9(b)为基于堆叠特征的CRF分类器分类结果,表4为分类结果的混淆矩阵。基于堆叠特征的CRF分类器的总体分类精度为66.78%,且分类结果中用红色标记的林地分类结果较差。

      图  9  基于堆叠特征的CRF分类结果

      Figure 9.  Classification result of CRF model based on stacking features

      表 4  堆叠特征分类精度

      Table 4.  Classification accuracy of stacking features

      类别 道路 林地 城郊 农田 其他 用户精度
      道路 122208 1456 6531 9828 18352 0.7435
      林地 1913 91840 50760 106297 8531 0.3541
      城郊 25348 9464 204174 16170 40917 0.6896
      农田 2552 5649 6351 64194 30475 0.5880
      其他 109049 150 10698 24033 473749 0.7670
      制图精度 0.4877 0.8539 0.7384 0.2925 0.8282
      总体分类精度:66.78%;Kappa系数:0.5547

      图  10  子分类器分类结果

      Figure 10.  Classification results based on different child classifiers

      图10为各个子分类器的分类结果及错分地物图。观察各个子分类器的分类结果,可以得出基于测量数据和分解特征的CRF分类器取得了较好的分类结果,但3个子分类器结果均比图9的分类结果差,表明测量数据特征和分解特征包含较多的表征Oberpfaffenhofen地区极化SAR图像的目标信息,但是各个子分类器均不能全面表征图像目标信息。同时观察发现,图10(a)中黑色方框标注的区域出现了明显的错分,但在图10(e)中此区域被基本正确分类。这是由于不同类型的特征包含了不同的目标信息,因此,各个子分类器对不同地物的分类能力也不同。

      图11为组合CRF分类器分类结果图及错分地物图,表5为组合CRF分类器的分类结果混淆矩阵。在本文实验中,通过使用各个子分类器先对图2(c)对应的图像数据进行分类得到3个总体分类精度,然后对3个分类精度归一化得到最终的 ${\alpha _i}\left( {i = 1,2,3} \right)$,取值分别为0.5783, 0.1330和0.2887。

      表 5  组合CRF分类精度

      Table 5.  Classification accuracy of combined CRF model

      类别 道路 林地 城郊 农田 其他 用户精度
      道路 93169 1970 7725 3805 57661 0.5668
      林地 434 201255 50737 4032 2883 0.7760
      城郊 955 10321 267059 9783 7955 0.9020
      农田 407 12851 18465 56485 20968 0.5174
      其他 22984 2355 23078 54041 515211 0.8341
      制图精度 0.7899 0.8798 0.7276 0.4406 0.8520
      总体分类精度:78.33%;Kappa系数:0.7013

      图12为基于l-1正则化法进行特征选择后的CRF分类器分类结果和错分地物图。观察对比发现其与组合CRF分类器的分类效果基本相同,进一步验证了组合分类器在解决特征干扰时的有效性。

      图13为SVM分类器对Oberpfaffenhofen地区极化SAR数据分类结果和错分地物图。与组合CRF分类器结果相比,SVM分类器得到的分类精度(63.58%,限于篇幅,我们未给出相应的分类精度表)与组合CRF分类器的分类精度(78.33%)存在较大的差距,这说明组合CRF分类器可以充分地利用特征信息和上下文信息来提高分类器性能。

      通过对比图11图13表4表5可以发现,组合CRF分类器的效果有较大的提升,而且组合分类器的分类精度达到78.33%,比基于堆叠特征的CRF分类器的分类精度(66.78%)提高了11.55%。这表明组合CRF分类器可以有效地解决基于堆叠特征的CRF分类器出现的特征干扰等问题,同时可以综合利用各个子分类器对不同地物分类性能的优势来提高分类精度。

      图  11  组合CRF分类器分类结果

      Figure 11.  Classification result of combined CRF model

      图  12  基于l-1特征选择分类器分类结果

      Figure 12.  Classification result of CRF model based on l-1 features selection

      图  13  SVM分类器分类结果

      Figure 13.  Classification result of SVM model

      图14为基于投票表决法的组合CRF分类器的分类结果和错分地物图,其总体分类精度为73.02% (限于篇幅,我们未给出相应的分类精度表)。上述实验结果表明,本文提出的组合CRF模型的分类精度优于投票表决法,同时也验证了本文提出的融合公式的有效性。

      图  14  基于投票表决法的组合CRF分类结果

      Figure 14.  Classification result of CRF model based on majority voting

    • 在极化SAR图像地物分类中,为了解决基于堆叠特征的CRF模型由于特征干扰和过拟合问题而导致分类器分类精度下降的问题,本文提出了一种组合CRF分类模型用于极化SAR图像地物分类。该模型首先将提取的极化SAR图像多种特征进行分组,然后利用各个特征组分别训练同一个CRF模型得到多个有差异性的子分类器,再通过融合公式融合各个子分类器的分类结果得到最终的组合CRF模型分类结果。利用两组实测极化SAR图像数据的实验结果表明,本文提出的组合CRF模型不仅能够有效地解决由于特征堆叠而引起的特征干扰和过拟合问题,而且能够充分地利用各个子分类器的性能优势,在对极化SAR图像分类时取得了较好的分类结果。

参考文献 (25)

目录

    /

    返回文章
    返回