2.793

2018影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

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

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

基于局部空间信息的可变类模糊阈值光学遥感图像分割

杨蕴 李玉 赵泉华

杨蕴, 李玉, 赵泉华. 基于局部空间信息的可变类模糊阈值光学遥感图像分割. 自动化学报, 2020, 46(x): 1−12. doi: 10.16383/j.aas.c190412
引用本文: 杨蕴, 李玉, 赵泉华. 基于局部空间信息的可变类模糊阈值光学遥感图像分割. 自动化学报, 2020, 46(x): 1−12. doi: 10.16383/j.aas.c190412
Yang Yun, Li Yu, Zhao Quan-Hua. Fuzzy threshold optical remote sensing image segmentation with variable class number based on local spatial information. Acta Automatica Sinica, 2020, 46(x): 1−12. doi: 10.16383/j.aas.c190412
Citation: Yang Yun, Li Yu, Zhao Quan-Hua. Fuzzy threshold optical remote sensing image segmentation with variable class number based on local spatial information. Acta Automatica Sinica, 2020, 46(x): 1−12. doi: 10.16383/j.aas.c190412

基于局部空间信息的可变类模糊阈值光学遥感图像分割


DOI: 10.16383/j.aas.c190412
详细信息
    作者简介:

    辽宁工程技术大学测绘与地理科学学院博士研究生. 主要研究方向高分辨遥感图像的地物目标几何以及特征提取.E-mail: m13147945981@163.com

    辽宁工程技术大学测绘与地理科学学院教授. 主要研究方向为遥感数据处理理论与应用基础研究. 本文通信作者.E-mail: liyu@lntu.edu.cn

    辽宁工程技术大学测绘与地理科学学院教授. 主要研究方向为遥感图像建模与分析随机几何在遥感图像处理中的应用.E-mail:zhaoquanhua@lntu.edu.cn

  • 基金项目:  国家自然科学基金 (41301479, 41271435)资助

Fuzzy Threshold Optical Remote Sensing Image Segmentation with Variable Class Number Based on Local Spatial Information

More Information
  • Fund Project:  Supported by National Natural Science Foundation of China (41301479, 41271435)
  • 摘要: 阈值法分割在光学遥感图像分析中被得到广泛的应用, 然而传统阈值法也存在诸多局限性, 如对噪声敏感, 需人为设定类别数, 计算复杂度高等. 针对传统阈值法的局限性, 提出一种基于局部空间信息的可变类模糊阈值光学遥感图像分割方法. 首先, 以图像光谱的一阶矩为初始类中心, 利用二分法原理和区域间最大相似度准则来快速确定类别数及其中心. 然后, 通过岭形模糊隶属函数计算各像素点对不同类的隶属程度, 同时考虑到像素点的隶属度局部空间信息, 在隶属度域中定义一个模糊加权滤波器对各类的隶属度矩阵进行滤波, 以滤波后的隶属度集合为依据, 按照最大隶属原则确定图像的标号场. 最后, 对标号场中的局部异常标号进行替换, 将修正后的标号场由对应的类中心赋色得到分割图像. 视觉和统计分析评价结果表明, 与传统阈值法相比, 该方法能在减少计算时间的同时获得更好的分割结果, 可适用于光学遥感图像的多阈值分割.
  • 图  1  多级岭形隶属函数

    Fig.  1  Multilevel ridge membership function

    图  2  算法流程图

    Fig.  2  Algorithm flow chart

    图  3  模拟图像

    Fig.  3  Simulated images

    图  4  模拟图像分割结果

    Fig.  4  Simulated image segmentation results

    图  5  全色遥感图像和分割结果

    Fig.  5  Panchromatic remote sensing images and segmentation results

    图  6  多光谱遥感图像和分割结果

    Fig.  6  Multispectral remote sensing images and segmentation results

    表  1  各同质区域的高斯分布参数

    Table  1  Gaussian distribution parameters of homogeneous regions

    模拟图像 参数
    图3(b1) 均值 70 90 130 180 160
    方差 6 2 7 4 8
    图3(b2) 均值 20/120/40 70/80/200 120/160/80 150/60/160 200/200/110
    方差 5/7/4 7/5/3 4/2/7 3/4/5 5/6/2
    下载: 导出CSV

    表  2  模拟图像分割的定量评价结果

    Table  2  Quantitative evaluation results of simulated image segmentation

    图像 指标 区域 Kmeans FCM 文献[19] 本文方法
    图3(b1) 用户精度/% 69.8 72.1 96.9 99.9
    92.2 90.1 98.6 99.9
    44.7 81.7 98.1 99.9
    90.4 80.4 97.0 99.8
    58.8 69.7 94.5 99.3
    产品精度/% 56.3 66.5 99.6 99.7
    88.5 86.1 78.3 100
    33.9 75.6 98.5 99.5
    86.5 86.5 97.1 100
    63.3 71.7 90.4 99.6
    总精度/% 55.4 75.2 88.6 99.4
    Kappa系数/% 53.9 74.6 85.3 99.7
    图3(b2) 用户精度/% 42.4 96.5 97.5 99.5
    38.9 84.5 70.6 96.0
    63.2 96.1 96.6 99.5
    85.6 59.6 95.9 99.0
    88.4 86.7 88.3 97.2
    产品精度/% 55.3 69.7 90.4 98.7
    55.3 90.5 76.8 99.9
    48.7 81.7 88.4 93.7
    78.3 71.4 69.5 98.5
    90.1 73.8 90.3 99.2
    总精度/% 53.2 81.2 89.7 98.3
    Kappa系数/% 48.9 80.0 88.1 98.6
    下载: 导出CSV

    表  3  全色遥感图像分割质量评价指标

    Table  3  Quality evaluation of panchromatic remote sensing image segmentation

    图像 MV JM
    Kmeans FCM 文献[19] 本文方法 Kmeans FCM 文献[19] 本文方法
    图5(a1) 1.626 1.335 0.973 0.632 0.887 0.742 0.712 0.633
    图5(b1) 2.344 1.698 1.335 0.966 0.831 0.787 0.737 0.596
    图5(c1) 1.886 1.475 1.203 1.079 0.759 0.703 0.663 0.645
    图5(d1) 0.982 0.875 0.619 0.512 0.692 0.640 0.598 0.582
    下载: 导出CSV

    表  4  计算复杂度对比

    Table  4  Computational complexity comparison

    方法 计算复杂度
    Kmeans $O((K+M \times N / K) \times t)$
    FCM $O((M \times N \times K \times t)$
    文献[19] $O\left((1+M \times N) \times K \times t \times \omega^{2}\right)$
    本文方法 $O\left(M \times N \times K+2 \times M \times N \times \omega^{2}\right)$
    下载: 导出CSV

    表  5  全色图像分割时间对比(s)

    Table  5  Panchromatic images segmentation time comparison (s)

    方法 图5(a1) 图5(b1) 图5(c1) 图5(d1)
    Kmeans 1.49 1.62 4.17 2.38
    FCM 2.63 2.87 7.46 4.28
    文献[19] 38.04 41.66 110.31 60.48
    本文方法 1.66 1.80 3.87 2.29
    下载: 导出CSV

    表  6  多光谱遥感图像分割质量评价

    Table  6  Quality evaluation of multispectral remote sensing image segmentation

    指标 方法 图6(a1) 图6(b1) 图6(c1) 图6(d1)
    MV Kmeans 1.971 1.613 2.316 2.146
    FCM 1.813 1.404 1.833 1.799
    文献[19] 1.570 1.071 1.279 1.344
    本文方法 1.376 0.796 0.941 1.001
    JM Kmeans 0.832 0.797 0.774 0.808
    FCM 0.748 0.624 0.647 0.734
    文献[19] 0.662 0.588 0.541 0.631
    本文方法 0.575 0.534 0.532 0.565
    E Kmeans 0.671 0.572 0.607 0.632
    FCM 0.524 0.466 0.573 0.597
    文献[19] 0.456 0.403 0.434 0.463
    本文方法 0.347 0.332 0.293 0.306
    下载: 导出CSV

    表  7  多光谱图像分割计算时间对比(s)

    Table  7  Multispectral images segmentation time comparison (s)

    方法 图6(a1) 图6(b1) 图6(c1) 图6(d1)
    Kmeans 4.15 2.75 3.01 5.39
    FCM 7.33 4.99 5.34 9.54
    文献[19] 106.02 66.65 77.49 138.51
    本文方法 3.56 3.15 3.37 4.68
    下载: 导出CSV
  • [1] Min Wang, Jiru Huang, Dongping Ming. Region-line association constraints for high-resolution image segmentation. IEEE Journal of Selected Topics in Applied Earth Observations & Remote Sensing, 2016, 10(2): 1−10
    [2] Andrés Troya-Galvis, Gancarski P, Berti-Equille L. Remote sensing image analysis by aggregation of segmentation-classification collaborative agents. Pattern Recognition, 2018, 73: 259−274 doi:  10.1016/j.patcog.2017.08.030
    [3] 李修霞, 荆林海, 李慧, 唐韵玮, 戈文艳. 参考1维光谱差异的区域生长种子点选取方法. 中国图象图形学报, 2016, 21(9): 1256−1264 doi:  10.11834/jig.20160915

    Li Xiu-Xia, Jing Lin-Hai, Li Hui, Tang Yun-Wei, Ge Wen-Yan. Seed extraction method for seeded region growing based on one-dimensional spectral differences. Journal of Image and Graphics, 2016, 21(9): 1256−1264 doi:  10.11834/jig.20160915
    [4] 游江, 唐力伟, 邓士杰, 苏续军. 完全基于边缘信息的目标靶快速分割算法. 激光与红外, 2017, 47(3): 372−378 doi:  10.3969/j.issn.1001-5078.2017.03.023

    You Jiang, Tang Li-Wei, Deng Shi-Jie, Su Xue-jun. Fast target segmentation algorithm fully based on edge information. Laser & Infrared, 2017, 47(3): 372−378 doi:  10.3969/j.issn.1001-5078.2017.03.023
    [5] Funke J, Tschopp F D, Grisaitis W, Sheridan A, Singh C, Saalfeld Se. Large scale image segmentation with structured loss based deep learning for connectome reconstruction. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2018, 41(7): 1669−1680
    [6] 李擎, 唐欢, 迟健男, 邢永跃, 李华通. 基于改进最大类间方差法的手势分割方法研究. 自动化学报, 2017, 43(4): 528−537

    Li Qing, Tang Huan, Chi Jian-Nan, Xing Yong-Yue, Li Hua-Tong. Gesture segmentation with improved maximum between-cluster variance algorithm. Acta Automatica Sinica, 2017, 43(4): 528−537
    [7] Zeggada A, Melgani, Farid, Ba zi, Yakoub. A deep learning approach to UAV image multilabeling. IEEE Geoscience & Remote Sensing Letters, 2017, 14(5): 694−698
    [8] 田娟秀, 刘国才, 谷珊珊, 鞠忠建, 刘劲光, 顾冬冬. 医学图像分析深度学习方法研究与挑战. 自动化学报, 2018, 44(3): 401−424

    Tian Juan-Xiu, Liu Guo-Cai, Gu Shan-Shan, Ju Zhong-Jian, Liu Jin-Guang, GU Dong-Dong. Deep learning in medical image analysis and its challenges. Acta Automatica Sinica, 2018, 44(3): 401−424
    [9] Fredo A R J, Abilash R S, Kumar C S. Segmentation and analysis of damages in composite images using multi-level threshold methods and geometrical features. Measurement, 2017, 100(100): 270−278
    [10] Yamini B, Sabitha R. Image steganalysis: adaptive color image segmentation using Otsu’s method. Journal of Computational & Theoretical Nanoscience, 2017, 14(9): 4502−4507
    [11] 黄扬, 郭立君, 张荣. 融合全局和局部相关熵的图像分割. 中国图象图形学报, 2018, 20(12): 1619−1628

    Huang Yang, Guo Li-Jun, Zhang Rong. Integration of global and local correntropy image segmentation algorithm. Journal of Image and Graphics, 2018, 20(12): 1619−1628
    [12] Choy S K, Shu Y L, Yu K W, Lee W Y, Leung K T. Fuzzy model-based clustering and its application in image segmentation. Pattern Recognition, 2017, 100(68): 141−157
    [13] Kurt B, Nabiyev V V, Turhan K. A novel automatic suspicious mass regions identification using Havrda & Charvat entropy and Otsu's N thresholding. Computer Methods & Programs in Biomedicine, 2014, 114(3): 349−360
    [14] Mousavirad S J, Ebrahimpour-Komleh H. Multilevel image thresholding using entropy of histogram and recently developed population-based metaheuristic algorithms. Evolutionary Intelligence, 2017, 10(1): 45−75
    [15] 产思贤, 周小龙, 张卓, 陈胜勇. 一种基于超像素的肿瘤自动攻击交互式分割算法. 自动化学报, 2017, 43(10): 1829−1840

    Chan Si-Xian, Zhou Xiao-Long, Zhang Zhuo, Chen Sheng-Yong. Interactive multi-label image segmentation with multi-layer tumors automata. Acta Automatica Sinica, 2017, 43(10): 1829−1840
    [16] 姜枫, 顾庆, 郝慧珍, 李娜, 郭延文, 陈道蓄. 基于内容的图像分割方法综述. 软件学报, 2017, 28(1): 160−183

    Jiang Feng, Gu Qing, Hao Hui-Zhen, Li Na, Guo Yan-Wen, Chen Dao-Xu. Survey on content-based image segmentation methods. Journal of Software, 2017, 28(1): 160−183
    [17] 肖满生, 文志诚, 张居武, 汪新凡. 一种改进隶属度函数的FCM聚类算法. 控制与决策, 2015, 30(12): 2270−2274

    Xiao Man-Sheng, Wen Zhi-Cheng, Zhang Ju-Wu, Wan Xin-Fan. An FCM clustering algorithm with improved membership function. Control and Decision, 2015, 30(12): 2270−2274
    [18] Meena Prakash R, Shantha S K R. Fuzzy C means integrated with spatial information and contrast enhancement for segmentation of MR brain images. International Journal of Imaging Systems & Technology, 2016, 26(2): 116−123
    [19] 赵雪梅, 李玉, 赵泉华. 参数自适应的可变类FLICM灰度图像分割算法. 控制与决策, 2017, 32(2): 262−268

    Zhao Xue-Mei, Li Yu, Zhao Quan-Hua. Self-adaptive FLICM algorithm for gray image segmentation with unknown number of clusters. Control and Decision, 2017, 32(2): 262−268
    [20] Ning J, Zhang L, Zhang D, Chengke W. Interactive image segmentation by maximal similarity based region merging. Pattern Recognition, 2010, 43(2): 445−456 doi:  10.1016/j.patcog.2009.03.004
    [21] 杜茂康, 王忠思, 宋强. 基于Bhattacharyya系数的改进相似度度量方法. 重庆邮电大学学报(自然科学版), 2018, 30(05): 115−120

    Du Mao-Kang, Wang Zhong-Si, Song Qiang. Research of improving similarity measure based on Bhattacharyya coefficient. Journal of Chongqing University of Posts and Telecommunications (Natural Science Edition), 2018, 30(05): 115−120
    [22] Choi S H, Jung H Y, Kim H. Ridge fuzzy regression model. International Journal of Fuzzy Systems, 2019, 21(7): 2077−2090 doi:  10.1007/s40815-019-00692-0
    [23] Mansoori G, Zolghadri J, Katebi D. A weighting function for improving fuzzy classification systems performance. Fuzzy Sets & Systems, 2007, 158(5): 583−591
    [24] Wang Y, Qi Q, Liu Y. Unsupervised segmentation evaluation using area-weighted variance and Jeffries-Matusita distance for remote sensing images. Remote Sensing, 2018, 10(8): 1193−1217 doi:  10.3390/rs10081193
    [25] Chen H C, Wang S J. Visible colour difference-based quantitative evaluation of colour segmentation. IEE Proceedings-Vision, Image and Signal Processing, 2006, 153(5): 598−609 doi:  10.1049/ip-vis:20045221
    [26] 李玉, 徐艳, 赵雪梅, 赵泉华. 利用高斯混合模型的多光谱图像模糊聚类分割. 光学精密工程, 2017, 25(2): 509−518 doi:  10.3788/OPE.20172502.0509

    Li Yu, Xu Yan, Zhao Xue-Mei, Zhao Quan-Hua. Multispectral image segmentation by fuzzy clustering algorithm used Gaussian mixture model. Optics and Precision Engineering, 2017, 25(2): 509−518 doi:  10.3788/OPE.20172502.0509
    [27] Yellamraju T, Boutin M. Clusterability and clustering of images and other “real” high-dimensional data. IEEE Transactions on Image Processing, 2018, 27(4): 1927−1938 doi:  10.1109/TIP.2017.2789327
  • [1] 孙文燕, 董恩清, 曹祝楼, 郑强. 一种基于模糊主动轮廓的鲁棒局部分割方法[J]. 自动化学报, doi: 10.16383/j.aas.2017.c160260
    [2] 许允喜, 陈方. 特征联合和旋转不变空间分割联合的局部图像描述符[J]. 自动化学报, doi: 10.16383/j.aas.2016.c150206
    [3] 邵昊阳, 张英涛, 鲜敏, 李致勋, 唐降龙. 基于多域先验的乳腺超声图像协同分割[J]. 自动化学报, doi: 10.16383/j.aas.2016.c150199
    [4] 刘丁, 张新雨, 陈亚军. 基于多目标人工鱼群算法的硅单晶直径检测图像阈值分割方法[J]. 自动化学报, doi: 10.16383/j.aas.2016.c150587
    [5] 李文平, 杨静, 印桂生, 张健沛. 数据场典型相关分析及其在图像分割中的应用[J]. 自动化学报, doi: 10.16383/j.aas.2015.c130896
    [6] 张瑞, 冯象初, 王斯琪, 常莉红. 基于稀疏梯度场的非局部图像去噪算法[J]. 自动化学报, doi: 10.16383/j.aas.2015.c140903
    [7] 谭乐怡, 王守觉. 基于双重超像素集的快速路径相似度图像分割算法[J]. 自动化学报, doi: doi{10.3724/SP.J.1004.2013.01653
    [8] 龙建武, 申铉京, 陈海鹏. 自适应最小误差阈值分割算法[J]. 自动化学报, doi: 10.3724/SP.J.1004.2012.01134
    [9] 王骏, 王士同, 邓赵红, 应文豪. 面向小目标图像的快速核密度估计图像阈值分割算法[J]. 自动化学报, doi: 10.3724/SP.J.1004.2012.01679
    [10] 刘松涛, 王慧丽, 殷福亮. 基于图割和模糊连接度的交互式舰船红外图像分割方法[J]. 自动化学报, doi: 10.3724/SP.J.1004.2012.01735
    [11] 魏巍, 申铉京, 千庆姬. 工业检测图像灰度波动变换自适应阈值分割算法[J]. 自动化学报, doi: 10.3724/SP.J.1004.2011.00944
    [12] 刘博, 黄剑华, 唐降龙, 刘家锋, 张英涛. 结合全局概率密度差异与局部灰度拟合的超声图像分割[J]. 自动化学报, doi: 10.3724/SP.J.1004.2010.00951
    [13] 侯彪, 徐婧, 刘凤, 焦李成. 基于第二代Bandelet域隐马尔科夫树模型的图像分割[J]. 自动化学报, doi: 10.3724/SP.J.1004.2009.00498
    [14] 范九伦, 雷博. 灰度图像最小误差阈值分割法的二维推广[J]. 自动化学报, doi: 10.3724/SP.J.1004.2009.00386
    [15] 康文雄, 邓飞其. 基于方向场分布率的静脉图像分割方法[J]. 自动化学报, doi: 10.3724/SP.J.1004.2009.01496
    [16] 李启翮, 罗予频, 萧德云. 基于流向标量场与快速扫掠法的图像分割[J]. 自动化学报, doi: 10.3724/SP.J.1004.2008.00993
    [17] 谢衍涛, 桑农, 张天序. 基于自适应隶属度函数的特征选择[J]. 自动化学报
    [18] 陈果, 左洪福. 图像的自适应模糊阈值分割法[J]. 自动化学报
    [19] 张恩勤, 施颂椒, 翁正新. 采用三角型隶属度函数的模糊系统的插值特性[J]. 自动化学报
    [20] 刘健庄, 栗文青. 灰度图像的二维Otsu自动阈值分割法[J]. 自动化学报
  • 加载中
计量
  • 文章访问数:  6
  • HTML全文浏览量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-27
  • 录用日期:  2019-12-02

基于局部空间信息的可变类模糊阈值光学遥感图像分割

doi: 10.16383/j.aas.c190412
    基金项目:  国家自然科学基金 (41301479, 41271435)资助
    作者简介:

    辽宁工程技术大学测绘与地理科学学院博士研究生. 主要研究方向高分辨遥感图像的地物目标几何以及特征提取.E-mail: m13147945981@163.com

    辽宁工程技术大学测绘与地理科学学院教授. 主要研究方向为遥感数据处理理论与应用基础研究. 本文通信作者.E-mail: liyu@lntu.edu.cn

    辽宁工程技术大学测绘与地理科学学院教授. 主要研究方向为遥感图像建模与分析随机几何在遥感图像处理中的应用.E-mail:zhaoquanhua@lntu.edu.cn

摘要: 阈值法分割在光学遥感图像分析中被得到广泛的应用, 然而传统阈值法也存在诸多局限性, 如对噪声敏感, 需人为设定类别数, 计算复杂度高等. 针对传统阈值法的局限性, 提出一种基于局部空间信息的可变类模糊阈值光学遥感图像分割方法. 首先, 以图像光谱的一阶矩为初始类中心, 利用二分法原理和区域间最大相似度准则来快速确定类别数及其中心. 然后, 通过岭形模糊隶属函数计算各像素点对不同类的隶属程度, 同时考虑到像素点的隶属度局部空间信息, 在隶属度域中定义一个模糊加权滤波器对各类的隶属度矩阵进行滤波, 以滤波后的隶属度集合为依据, 按照最大隶属原则确定图像的标号场. 最后, 对标号场中的局部异常标号进行替换, 将修正后的标号场由对应的类中心赋色得到分割图像. 视觉和统计分析评价结果表明, 与传统阈值法相比, 该方法能在减少计算时间的同时获得更好的分割结果, 可适用于光学遥感图像的多阈值分割.

English Abstract

杨蕴, 李玉, 赵泉华. 基于局部空间信息的可变类模糊阈值光学遥感图像分割. 自动化学报, 2020, 46(x): 1−12. doi: 10.16383/j.aas.c190412
引用本文: 杨蕴, 李玉, 赵泉华. 基于局部空间信息的可变类模糊阈值光学遥感图像分割. 自动化学报, 2020, 46(x): 1−12. doi: 10.16383/j.aas.c190412
Yang Yun, Li Yu, Zhao Quan-Hua. Fuzzy threshold optical remote sensing image segmentation with variable class number based on local spatial information. Acta Automatica Sinica, 2020, 46(x): 1−12. doi: 10.16383/j.aas.c190412
Citation: Yang Yun, Li Yu, Zhao Quan-Hua. Fuzzy threshold optical remote sensing image segmentation with variable class number based on local spatial information. Acta Automatica Sinica, 2020, 46(x): 1−12. doi: 10.16383/j.aas.c190412
  • 随着遥感图像分辨率的提高, 其在人们生活中的使用比例越来越大, 如何快速、精确地从其中提取目标信息已成为遥感领域的一个重要研究课题[1]. 分辨率的提高在带来信息丰富的同时也给目标分析带来巨大困难, 而图像分割技术通过给每个像素分配标签, 使同一标签的像素具有相同的视觉特性, 可提高图像分析的效率[2]. 因此执行图像分割用于更好描述图像是由图像处理到图像分析的关键步骤, 对于遥感图像目标的高效分析具有重要意义.

    现有图像分割方法大致可分为四类: 基于区域[3]、边缘[4]、学习分类[5]和阈值[6] 的方法. 其中, 区域法的优点是对噪声不敏感, 但需要制定合适且高效的相似性准则, 易造成图像过分割; 而单纯基于边缘检测的方法往往不能提供较好的区域结构; 学习分类法利用大量的训练样本能较好地解决图像中的噪声和均匀问题, 但目前存在网络模型难以确定, 易出现欠学习和局部最优等问题[7-8]; 阈值法具有简单易行、性能稳定等优点[9], 目前流行的阈值方法有Otsu法[10]、最大熵法[11]和聚类法[12]等, 这些方法基于不同准则来选取阈值, 都具有各自的优点, 但也有自身的缺点. Otsu和最大熵法分别以分割后图像的类间方差和熵最大为目标函数, 穷举选取使目标函数获得最大值的阈值, 但其仅考虑像素本身对分割结果的影响, 对噪声和异常值敏感[13]. 为此, 人们引入邻域像素进行二维阈值化, 在一定程度上减弱了噪声和异常值对分割结果的影响, 但存在信息损失严重、耗时过长等缺点, 特别是在多阈值图像分割中更为突出[14]. 聚类法是一种不依赖于分类的先验知识, 遵循同质性准则对像素划分的方法, 能够有效分离具有全局特征的数据[15], 其中最具代表是Kmeans和模糊C 均值(Fuzzy C-Means, FCM)算法. Kmeans算法以欧式距离作为相似度测度, 以误差平方和为聚类准则函数, 迭代计算使得准则函数收敛为止[16]. 它具有简单、计算快速等优点, 但其强制图像的各个像素属于特定的类别, 使得对具有复杂场景的遥感图像的分割结果很差. FCM 方法利用模糊集理论中的隶属函数对每个像素点进行划分, 考虑了模糊性和不确定性, 可对纹理和背景简单的图像取得良好的效果[17], 但FCM仅考虑像素光谱信息, 对噪声较敏感. 为此常将图像局部空间信息引入目标函数, 如文献[18-19] 中, 通过迭代计算像素邻域和聚类中心之间的距离来提高FCM算法对图像分割的鲁棒性. 然而, 局部空间信息的引入会导致算法复杂度升高, 降低了该算法的实际应用价值.

    为了克服传统阈值法需人为设定类别数、对噪声敏感、计算复杂度高等局限性, 提出了一种基于局部空间信息的可变类模糊阈值遥感图像分割方法. 首先, 以图像中待定类像素点的光谱一阶矩和二阶矩为初始类中心和自适应阈值, 利用二分法原理来初始化类别数目. 对于任意待定类别, 当相邻两次迭代得到类别的像素光谱均值差小于规定值时, 该类查找结束, 并以此均值作为该类的中心. 并在确定各类的迭代过程中, 计算与其二阶邻域系统间的相似度, 以最大相似度准则进行类别区域合并. 由确定的类别数和类中心为输入, 通过岭形模糊隶属函数计算图像中各像素点对于不同类别的隶属程度. 考虑图像中属于同一类别的像素具有一定程度的连通性, 即各像素与其邻域之间关系紧密, 不可能存在孤立的像素, 而隶属度和标号场是各像素属性的映射, 因此也满足此局部空间关系. 为了利用此局部空间信息同时降低计算复杂度, 在隶属度域中定义一个加权模糊滤波器对每一类别中像素点的隶属度进行滤波, 以滤波后的隶属度集合为依据, 按照最大隶属原则确定图像的标号场. 最后, 对标号场中局部异常标号进行替换, 将修正后的标号场由对应的类中心赋色得到分割图像.

    • ${{I}} \!=\! \{{{I}}_{p}(i, j):(i, j) \in {{S}}, i \!=\! 1,2,\cdots, M, j = 1,$ $ 2, \cdots, N\} $ 为定义在图像空域 $ {{S}} $ 、大小为 $ M \times N $ 的待分割遥感影像, 其中, ${I}_p(i, j) \!=\! (I^1_ p(i, j), I^2 _p(i, j), \cdots,$ $ I^l _p(i, j))^{\rm T} $ 为像素光谱测度矢量, $ p $ 为像素索引, $ (i, j) $ 表示像素 $ p $ 的格点位置, $ l $ 为光谱测度矢量的维度. 对于全色遥感图像, $ l = 1; $ 对于多光谱遥感图像, $ l \ge 3 .$

    • 利用图像的全局信息快速地对分割图像的类别数及其中心粗略估计, 其包含以下步骤: 类中心初始化、划分阈值选取、待定类识别、类别区域融合和查找结束.

    • 为了避免随机选择初始值而降低类别查找的动态性, 本文以图像光谱的一阶矩作为初始类中心, 即

      $$ \begin{array}{l} {{{C}}}_s^0 = \dfrac{{\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {{I^s}\left( {i,j} \right)} } }}{{M \times N}} \end{array} $$ (1)

      其中, $ s\in\{1, 2,\cdots, l\} $ 为波段索引, $ {{C}}^0_ s $ 为波段 $ s $ 的初始类中心.

    • 以图像中待定类各像素点与类中心的二阶中心矩的算术平方根为自适应阈值, 较之各待定类像素点与类中心的距离从而确定类属划分, 波段 $ s $ 中第 $ k $ 个类别的自适应阈值为

      $$ \begin{array}{l} T_s^k = \sqrt {\dfrac{{\sum\limits_{i = 1}^{M'} {\sum\limits_{j = 1}^{N'} {{{\left( {{I^s}\left( {i,j} \right) - {{{C}}}_s^k} \right)}^2}} } }}{{M' \times N'}}} \end{array} $$ (2)

      其中, $ s \in \{1, 2, \cdots, l\} , k\in \{1, 2, \cdots, K\} $ 为类别索引, $ K $ 为确定的类别数目, $ M' \times N' $ 为确定第 $ k $ 类中心时图像中待定类像素点个数, 则第 $ k $ 个类别的自适应阈值为 $ { T}^k = [T^k_1, T^k _2, \cdots, T^k _l]. $

    • 当查找第 $ k $ 类时, 以欧式距离为准则, 计算波段 $ s $ 中各待定类像素点与类中心数值特征的绝对差异,

      $$ \begin{array}{l} D_s^k\left( {i,j} \right) = \sqrt {{{\left( {{I^s}\left( {i,j} \right) - C_s^k\left( {i,j} \right)} \right)}^2}} \end{array} $$ (3)

      其中, $i \!\!=\!\! 1, 2, \cdots, M', j \!\!=\!\! 1, 2, \cdots, N' , {{D}}^k \!\!=\!\! [D^k_1, D^k _2,\cdots,$ $D^k _l].$ 标记 $ {{D}}^k $ 小于 $ {{T}}^k $ 的像素点, 以这些像素点光谱测度均值更新第 $ k $ 类的聚类中心 $ {{C}}^k _s ,$ 当相邻两次迭代得到类别的光谱均值差小于0.87 (以各波段坐标均相差0.5为限)时, 该类查找结束, 并以此均值作为该类的中心.

    • 考虑到类别查找时, 由于遥感图像的局部非均匀性, 在区域上会出现局部最优现象, 为此通过计算相邻区域的相似度[20]来指导合并过程, 确定类别数. 其将图像每一维度均匀量化为16个等级, 然后在特征空间中计算每个类别区域的归一化直方图, 使用Bhattacharyya系数[21] 来度量相邻类别区域间的相似性, 可表示为

      $$ \begin{array}{l} \rho \left( {{{{R}}},{{{Q}}}} \right) = \displaystyle\sum\limits_{u = 1}^{{{16}^l}} {\sqrt {{\rm{Hist}}_{{{R}}}^u \cdot {\rm{Hist}}_{{{Q}}}^u} } \end{array} $$ (4)

      其中, $ l $ 为光谱测度矢量的维度, $ {\rm{Hist}}_{ R} $ $ {\rm{Hist}}_{ Q} $ 分别是类别区域 $ {{R}} $ $ {{Q}} $ 的归一化直方图, 上标 $ u $ 代表它们的第 $ u $ 个元素. 设以 $ {{R}} $ 为中心的二阶邻域系统内所有类别区域为 $ {{Q}}_k , k \in \{1, 2, \cdots, 8\} ,$ 计算相似度 $ \rho({{R}}, {{Q}}_k), $ 取其最大值对应的类别区域索引, 可示为

      $$ \begin{array}{l} {\rm{km}} = \mathop {\arg \max }\limits_k \left( {\rho \left( {{{{R}}},{{{{Q}}}_k}} \right)} \right) \end{array} $$ (5)

      其中, $ {\rm km}\in\{1, 2, \cdots, 8\}. $ $ \rho({{R}}, {{Q}}_{\rm km}) $ 大于阈值, 将 $ {{R}} $ $ {{Q}}_{\rm km} $ 合并, 并更新类别及其中心, 其中阈值越大, 合并条件越严格, 反之亦然, 本文取0.85.

    • 当第 $ k $ 类确定后, 在待定类的像素点中重复执行上面步骤, 直到所有像素完成标记, 类别查找结束.

    • 定义一个模糊隶属函数 $ \mu_k( \cdot ) $ , 由初始化的类别数和类中心使各像素与每个类相关联, 得到各像素所属类别的隶属度. 岭形分布隶属函数(见图1)采用正弦函数作为核心, 具有主值区间宽、过渡平稳的特点[22], 能很好反映各像素与类别隶属程度的模糊关系, 定义为

      图  1  多级岭形隶属函数

      Figure 1.  Multilevel ridge membership function

      $$ \begin{split} &{\mu _k}\left( {I_p^s\left( {i,j} \right)} \right) = \\ &\left\{\!\!\begin{array}{l} \begin{array}{*{20}{c}} {}&{}&{\begin{array}{*{20}{c}} {}&{}&{\begin{array}{*{20}{c}} 0&{\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {}&{} \end{array}}&{}&, \end{array}}&{I_p^s\left( {i,j} \right) \le a} \end{array}} \end{array}} \end{array}\\ \dfrac{1}{2} \!+\! \dfrac{1}{2}\sin \!\!\left( {\dfrac{\pi }{{b \!-\! a}}\left( {I_p^s\left( {i,j} \right) \!-\! \dfrac{{a \!+\! b}}{2}} \right)} \right),a \!<\! I_p^s\!\left( {i,j} \right) \!\le\! b\\ \dfrac{1}{2} \!-\! \dfrac{1}{2}\sin\!\! \left( {\dfrac{\pi }{{c \!-\! b}}\left( {I_p^s\left( {i,j} \right) \!-\! \dfrac{{b \!+\! c}}{2}} \right)} \right),b \!<\! I_p^s\left( {i,j} \right) \!\le\! c\\ \begin{array}{*{20}{c}} {}&{}&{\begin{array}{*{20}{c}} {}&{}&{\begin{array}{*{20}{c}} 0&{\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {}&{} \end{array}}&{}&, \end{array}}&{I_p^s\left( {i,j} \right) > c} \end{array}} \end{array}} \end{array} \end{array} \right. \end{split} $$ (6)

      其中, $ a , b, c $ 为模糊参数, 将其与类中心和图像相关联, 当 $ k = K = 1 $ 时, $ a = \min(I^s) ,b = {{C}}_K; $ $ k = K $ 时, $ a = {{C}}_k - 1, b = {{C}}_k,c = \max(I^s); $ $ 1 < k < K $ 时, $ a = {{C}}_k - 1, b = {{C}}_k , c = {{C}}_k + 1 $ .

      $ \mu_{k}\left(I_{p}^s(i, j)\right) $ 可得图像 $ I^s $ 中像素 $ p $ 属于第 $ k $ 类的隶属度, 其输出是一个隶属度向量,

      $$ \begin{split} {{\mu}}\left(I_{p}^{s}(i, j)\right) =& [\mu_{1}\left(I_{p}^{s}(i, j)\right) \mu_{2}\left(I_{p}^{s}(i, j)\right) \cdots \\ &\mu_{K}\left(I_{p}^{s}(i, j)\right)] \end{split} $$ (7)

      其中, $ K $ 为类别数. 图像中所有像素的隶属度向量组成一个隶属度集 $ {\mu} ({I}). $

    • 同一类别区域内部在分割过程由于随机噪声产生的孤立像素可认为是误分割, 在隶属度矩阵中表现为异常的隶属度值. 为了利用此局部空间信息同时降低计算复杂度, 将模糊加权函数[23] 应用于均值滤波器移动窗口内的隶属度值, 在隶属度域中定义一个模糊加权均值滤波器来过滤随机噪声, 以 $ \mu_{k}(i, j),k \in \{1,2, \cdots, K\} $ 为二维模糊滤波器的输入, 其输出为

      $$ \begin{array}{l} \mu_{k}^{F}(i, j) \!=\! \dfrac{\sum\limits_{(m, n) \in \eta_{p}} W_{A}\left[\mu_{k}(i\!+\!m, j\!+\!n)\right] \!\cdot\! \mu_{k}(i\!+\!m, j\!+\!n)}{\sum\limits_{(m, n) \in \eta_{p}} W_{A}\left[\mu_{k}(i\!+\!m, j\!+\!n)\right]} \end{array} $$ (8)
      $$ \begin{array}{l} {W_A}\left[ {\mu \left( {i + m,j + n} \right)} \right] = \left\{ \begin{array}{l} 1 - \dfrac{{{\mu _{\rm mav}}\left( {i,j} \right) - \mu \left( {i + m,j + n} \right)}}{{{\mu _{{\rm{mav}}}}\left( {i,j} \right) - {\mu _{{\rm{min}}}}\left( {i,j} \right)}},{\mu _{{\rm{min}}}}\left( {i,j} \right) \le \mu \left( {i + m,j + n} \right) \le {\mu _{{\rm{mav}}}}\left( {i,j} \right)\\ 1 - \dfrac{{\mu \left( {i + m,j + n} \right) - {\mu _{{\rm{mav}}}}\left( {i,j} \right)}}{{{\mu _{{\rm{max}}}}\left( {i,j} \right) - {\mu _{{\rm{mav}}}}\left( {i,j} \right)}},{\mu _{{\rm{mav}}}}\left( {i,j} \right) \le \mu \left( {i + m,j + n} \right) \le {\mu _{{\rm{max}}}}\left( {i,j} \right)\\ 1,{{\mu _{{\rm{mav}}}}\left( {i,j} \right) - {\mu _{{\rm{min}}}}\left( {i,j} \right) = 0/{\mu _{{\rm{max}}}}\left( {i,j} \right) - {\mu _{{\rm{mav}}}}\left( {i,j} \right) = 0} \end{array} \right. \end{array} $$ (9)

      其中, $ \mu_{k}^{F}(i, j) $ 为第 $ k $ 类隶属度空间格点位置坐标 $ (i,j) $ 滤波后的输出值; $ \eta_{p} $ 是以 $ p $ 为中心的窗口邻域, 本文 $ \eta_{p} = 25 $ ; $ W_{A}(\cdot) $ 是以移动窗口内均值构建的窗口加权函数, 定义为

      其中, $ (m, n) \in \eta_{p} , \mu_{\max}(i,j) $ $ \mu_{\min}(i,j) $ $ \mu_{\rm mav}(i,j) $ 分别是 $ \eta_{p} $ 内输入值的最大、最小和平均值.

    • 由滤波后的隶属函数, 按最大隶属原则得到标号场, 如式(10)所示

      $$ \begin{array}{l} {{{{L}}}_s} = \mathop {\arg \max }\limits_k \left\{ {\mu _k^F\left( {{I^s}} \right)} \right\} \end{array} \tag{10}$$ (10)

      其中, $k \!\in\!\{1,2, \cdots, K\}, { L}_{s} \!=\! \{L_{s}(i, j) \!=\! k:i \!=\! 1,2,\cdots,$ $M, j = 1,2, \cdots, N\}.$ 为了消除标号场 $ L_s $ 中由于误判的类别标号(假设表现为脉冲噪声)并保持其边缘, 以移动窗口内的中值构建窗口加权函数 $ W_M(\cdot) $ 来定义模糊加权中值滤波器, 如式(11)所示

      其中, $ \mu_{\rm med}(i,j) $ $ \eta_p $ 内所有输入值的中值. 经滤波后得到修正的标号场 $ {L}^F_s $ ,

      $$ \begin{split} &{{{L}}}_s^F\left( {i,j} \right) =\\ &\dfrac{{\sum\limits_{\left( {m,n} \right) \in {\eta _p}} {{W_M}\left[ {{{{{L}}}_s}\left( {i \!+\! m,j \!+\! n} \right)} \right] \!\cdot\! {{{{L}}}_s}\left( {i \!+\! m,j \!+\! n} \right)} }}{{\sum\limits_{\left( {m,n} \right) \in {\eta _p}} {{W_M}\left[ {{{{{L}}}_s}\left( {i + m,j + n} \right)} \right]} }} \end{split} \tag{12}$$ (12)

      其中, $ i = 1,2, \cdots, M ,j = 1,2, \cdots, N. $ $ {C}_s $ 的颜色值赋予 $ {L}^F_s $ 相应的标号位置得到分割图像.

    • 算法流程图如图2所示.

      图  2  算法流程图

      Figure 2.  Algorithm flow chart

      综上所述, 可对本文方法的流程总结如下.

      S1 利用二分法原理和区域间最大相似度准则对图像的类别数及其中心初始化;

      S2 将类中心和各像素与岭形分布隶属函数相关联, 得到各像素所属类别的隶属度;

      S3 在隶属度域中定义一个模糊加权均值滤波器来过滤由于噪声引起的异常隶属度值;

      S4 对滤波后的隶属度集按照最大隶属原则确定图像的标号场, 并由模糊加权中值滤波器来消除其中由于类别初始化时误判的类别标号;

      S5 对修正后的标号场由对应的类中心赋色得到分割图像.

    • 选择模拟图像、全色遥感图像和多光谱遥感图像对本文方法进行验证, 其中研究区的真彩色图像是由WorldView-2 的红、绿、蓝三波段(波段5 = 红, 波段3 = 绿, 波段2 = 蓝)合成.

    • 实验是在具有Intel(R) Core (TM) CPU, i7-6700, 3.4GHz, 16GB RAM的DELL计算机上, 应用Windows 8 环境下的MATLAB R2017a实现的. 为了验证提出方法的优越性, 将所提出方法与一些著名的阈值算法的性能进行比较, 选择了Kmeans、FCM和文献[19]的方法作为对比算法, 其中Kmeans 和FCM算法可在文献[16] 和文献[17]中得到详细描述, 对其分割结果从目视分析、监督评价、非监督评价和运行时间来进行比较分析. 对于模拟图像, 由于已知参考图像, 选择监督评价中应用最广泛的混淆矩阵来计算用户精度、产品精度、总精度和Kappa系数对分割结果进行定量评价. 对于全色遥感图像, 利用文献[24] 提出的面积加权方差(Area-Weighted Variance, WV)和Jeffries-Matusita(JM) 距离对分割结果进行非监督评价. 对多光谱遥感图像, 其光谱信息增加, 为了较全面评价分割质量, 除了利用WV和JM外, 增加文献[25]提出的可见色差E 对分割结果进行非监督评价. WV用于测量全局分割同质区域内的优度, 并按每个同质区域的面积进行加权. 其定义如下

      $$ \begin{array}{l} {W_M}\left[ {L\left( {i + m,j + n} \right)} \right] = \left\{ \begin{array}{l} 1 - \dfrac{{{L_{{\rm{med}}}}\left( {i,j} \right) - L\left( {i + m,j + n} \right)}}{{{L_{{\rm{med}}}}\left( {i,j} \right) - {L_{{\rm{min}}}}\left( {i,j} \right)}},{L_{{\rm{min}}}}\left( {i,j} \right) \le L\left( {i + m,j + n} \right) \le {L_{{\rm{med}}}}\left( {i,j} \right)\\ 1 - \dfrac{{L\left( {i + m,j + n} \right) - {L_{{\rm{med}}}}\left( {i,j} \right)}}{{{L_{{\rm{max}}}}\left( {i,j} \right) - {L_{{\rm{med}}}}\left( {i,j} \right)}},{L_{{\rm{med}}}}\left( {i,j} \right) \le L\left( {i + m,j + n} \right) \le {L_{{\rm{max}}}}\left( {i,j} \right)\\ 1,{{L_{{\rm{med}}}}\left( {i,j} \right) - {L_{{\rm{min}}}}\left( {i,j} \right) = 0/{L_{{\rm{max}}}}\left( {i,j} \right) - {L_{{\rm{med}}}}\left( {i,j} \right){\rm{ = }}0} \end{array} \right. \end{array}\tag{11} $$ (11)
      $$ \begin{array}{l} {\mathop{\rm WV}\nolimits} = \dfrac{{\sum\limits_{i = 1}^n {{a_i} \times {v_i}} }}{{\sum\limits_{i = 1}^n {{a_i}} }} \end{array} $$ (13)
      $$ \begin{array}{l} {v_i} = \dfrac{{\sum\limits_{b = 1}^m {{v_{ib}}} }}{m} \end{array} $$ (14)

      其中, $ i\in\{1,2,\cdots,n\} $ 为同质区域索引, $ n $ 是同质区域数, $ a_i $ 为区域 $ i $ 的面积, $ v_i $ 是区域 $ i $ 所有波段的平均方差, $ l $ 是遥感图像的波段数, $ v_{ib} $ $ b $ 波段中第 $ i $ 个区域的方差. 较低的WV表明分割结果具有较高的区域内同质性. JM用于测量分割同质区域间的整体优度, 并按每个同质区域的面积进行加权, 其定义如下

      $$ \begin{array}{l} {\mathop{\rm JM}\nolimits} = \dfrac{{\sum\limits_{i = 1}^n {{a_i} \times \left( {{ {1 \over l}}\sum\limits_{b = 1}^l {{J_i}\left( b \right)} } \right)} }}{{\sum\limits_{i = 1}^n {{a_i}} }} \end{array} $$ (15)

      其中, $ J_i(b) $ 是波段 $ b $ 中区域 $ i $ 的JM距离, 可表示为

      $$ \begin{array}{l} {J_i} = \sum\limits_{k = 1}^\Omega {\dfrac{{{L_{ik}}}}{{{L_i}}}} \times \left( {2\left( {1 - {e^{ - {B_{ik}}}}} \right)} \right) \end{array} $$ (16)

      其中, $ L_i $ 为区域 $ i $ 的边界长度, $ L_{ik} $ 为区域 $ i $ 及其相邻区域 $ k\in\{1, 2,\cdots,\Omega\} $ 的共同边界长度, $ B_{ik} $ 是根据区域邻接图计算区域 $ i $ $ k $ 的Bhattacharyya距离, 可表示为

      $$ \begin{array}{l} {B_{ik}} \!=\! \dfrac{1}{8}{\left( {{\mu _i} \!-\! {\mu _k}} \right)^2} \!\times\! \dfrac{2}{{v_i^2 \!+\! v_k^2}} \!+\! \dfrac{1}{2}\ln \left( {\dfrac{{v_i^2 \!+\! v_k^2}}{{2{v_i}{v_k}}}} \right) \end{array} $$ (17)

      其中, $ \mu_i , \mu_k $ $ v_i , v_k $ 分别是相邻区域 $ i $ $ k $ 的均值和方差. 较低的JM表明分割结果具有较高的区域间异质性.

      在可见色差的情况下, E结合区域内和区域间视觉误差来评价分割结果的质量, 定义为

      $$ \begin{array}{l} {\mathop{\rm E}\nolimits} = {{\mathop{\rm E}\nolimits} _{{\rm{intra}}}} + {{\mathop{\rm E}\nolimits} _{{\rm{inter}}}} \end{array} $$ (18)

      其中, $ \rm E_{ intra} $ 表示区域内视觉误差, 其旨在评估欠分割的程度, 可表示为

      $$ \begin{array}{l} {{\mathop{\rm E}\nolimits} _{\rm intra}} = \dfrac{{\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {u\left( {{{\left\| {I\left( {i,j} \right) - I'\left( {i,j} \right)} \right\|}_{{\mathop{\rm Lab}\nolimits} }} - th} \right)} } }}{{M \times N}} \end{array} $$ (19)

      其中, $ M\times N $ 为图像大小, $ I $ $ I' $ 分别表示原始和分割图像, $ \left\|\cdot \right\|_{\rm Lab} $ 表示CIELab颜色空间中的色差, $ th $ 表示可见色差的阈值, $ u(\cdot) $ 表示阶跃函数,

      $$ \begin{array}{l} u\left( t \right) = \left\{ \begin{array}{l} 1,\begin{array}{*{20}{c}} {}&{t > 0} \end{array}\\ 0,\begin{array}{*{20}{c}} {}&{{\mathop{\rm otherwise}\nolimits} } \end{array} \end{array} \right. \end{array} $$ (20)

      其中, $ th $ 为可调阈值, 本文 $ th = 6 $ . $ \rm E_{inter} $ 为区域间视觉误差, 用来评估过分割程度, 定义为

      $$ \begin{array}{l} \mathrm{E}_{\text {inter }} = \dfrac{\sum\limits_{i = 1}^{K} \sum\limits_{j = 1 \atop j \neq i}^{K} w_{i j} \times u\left(t h-\left\|I_{i}^{\prime}-I_{j}^{\prime}\right\|_{\rm Lab}\right)}{c \times M \times N} \end{array} $$ (21)

      其中, $ K $ 表示分割的同质区域数量, $ w_{ij} $ 表示区域 $ i $ 和区域 $ j $ 之间的连接长度, 如果区域 $ i $ $ j $ 未连接, 则等于零, $ c $ 表示归一化因子, 本文选择 $ c = 1/6. $ 由于 $ \rm E_{intra} $ $ \rm E_{inter} $ 之间存在着相互关系, 因此当E 达到较小的值时, 被视为实现较好的分割.

    • 图3(a1)(a2)是尺寸为 $ 128\times 128 $ 像素的模板, 用于生成标记为I-V的五个均匀区域的模拟图像(见图3(b1)图3(b2)), 其各波段中均匀区域的像素强度服从表1所示的高斯分布.

      图  3  模拟图像

      Figure 3.  Simulated images

      表 1  各同质区域的高斯分布参数

      Table 1.  Gaussian distribution parameters of homogeneous regions

      模拟图像 参数
      图3(b1) 均值 70 90 130 180 160
      方差 6 2 7 4 8
      图3(b2) 均值 20/120/40 70/80/200 120/160/80 150/60/160 200/200/110
      方差 5/7/4 7/5/3 4/2/7 3/4/5 5/6/2

      分别利用对比算法和本文方法对图3(b1)3(b2)进行分割实验, 分割结果如图4所示, 其中第一行和第二行分别为图3(b1)图3(b2)的分割结果. 由图4可以看出, Kmeans的分割结果中含有大量分割噪声, 这是因为其仅考虑像素的测量, 对光谱的变化很敏感, 因此很难区分不同的均匀区域. 对于FCM, 它利用隶属度对像素模糊划分, 隶属度提供的自由度使得比传统硬阈值划分具有更好的不确定性处理能力, 从而减弱噪声和异常值对分割结果的影响, 它可以成功地区分区域I、III和IV, 但区域II和V的分割结果比较模糊. 文献[19] 通过定义中心像素与邻域像素隶属度间的制约关系改善了FCM得到的分割结果, 但对部分区域不能有效分割, 如图4(c2)的区域II. 本文方法在正确识别类别数的情况下, 考虑了隶属度空间中各像素隶属度之间的空间信息, 从而加强了邻域系统中像素之间的连接, 此外在标号域中建立的模糊中值滤波剔除了局部异常标号, 取得了更好的分割结果.

      图  4  模拟图像分割结果

      Figure 4.  Simulated image segmentation results

      为了进行定量精度评价, 计算不同同质区域的误差矩阵, 从而计算出产品精度、用户精度、总体精度和Kappa系数[26] 来描述分割结果的准确性(见表2). 从表2可以看出, 本文方法各同质区域的用户精度和产品精度均在96% 以上, 总体精度和Kappa系数平均为98.8%和99.1%. 在对比算法中, Kmeans的总体精度和Kappa系数平均分别为54.3%和51.4%; FCM的总体精度和Kappa系数平均分别为78.1%和77.3%, 其优于Kmeans; 文献[19]的总体精度和Kappa系数分别为89.1%和87.5%, 优于Kmeans和FCM, 但精度仍低于本文方法. 定量结果表明, 本文所提出的方法可取得良好的分割效果, 具有广阔的应用前景.

      表 2  模拟图像分割的定量评价结果

      Table 2.  Quantitative evaluation results of simulated image segmentation

      图像 指标 区域 Kmeans FCM 文献[19] 本文方法
      图3(b1) 用户精度/% 69.8 72.1 96.9 99.9
      92.2 90.1 98.6 99.9
      44.7 81.7 98.1 99.9
      90.4 80.4 97.0 99.8
      58.8 69.7 94.5 99.3
      产品精度/% 56.3 66.5 99.6 99.7
      88.5 86.1 78.3 100
      33.9 75.6 98.5 99.5
      86.5 86.5 97.1 100
      63.3 71.7 90.4 99.6
      总精度/% 55.4 75.2 88.6 99.4
      Kappa系数/% 53.9 74.6 85.3 99.7
      图3(b2) 用户精度/% 42.4 96.5 97.5 99.5
      38.9 84.5 70.6 96.0
      63.2 96.1 96.6 99.5
      85.6 59.6 95.9 99.0
      88.4 86.7 88.3 97.2
      产品精度/% 55.3 69.7 90.4 98.7
      55.3 90.5 76.8 99.9
      48.7 81.7 88.4 93.7
      78.3 71.4 69.5 98.5
      90.1 73.8 90.3 99.2
      总精度/% 53.2 81.2 89.7 98.3
      Kappa系数/% 48.9 80.0 88.1 98.6
    • 图5(a1)-(d1)为待分割的5幅全色遥感图像, 其中, 图5(a1)是分辨率为0.5 m, 尺度为 $ 322\times 288 $ 像素的Pleiades-1 机场图像, 人为判读类别数为4; 图5(b1)是分辨率为1 m, 尺度为 $ 270 \times 350 $ 像素的IKONOS郊区图像, 人为判读类别数为5; 图5(c1)(d1)是分辨率为0.7 m 的EROS-B图像, 其中, 图5(c1) $ 512\times 512 $ 像素的码头图像, 人为判读类别数为6, 图5(d1) $ 256\times 256 $ 像素的房屋图像, 人为判读其类别数为7. 采用本文算法和对比算法对全色遥感图像进行分割实验, 其中图5(a2)-(d2) 为采用Kmeans方法进行分割实验得到的结果; 图5(a3)-(d3)为采用FCM方法进行分割实验得到的结果; 图5(a4)-(d4)为采用文献[19]方法进行分割实验得到的结果; 5(a5)-(d5)为采用本文方法进行分割实验得到的结果.

      图  5  全色遥感图像和分割结果

      Figure 5.  Panchromatic remote sensing images and segmentation results

      图5可以看出, 通过各方法的分割, 图像表现为不同的同质区域, 简化了图像表示, 有利于图像的后续分析. 对Kmeans 法, 各同质区域叠加在一起, 使其含有大量的分割噪声, 分割结果较差; FCM法分割精度有所提高, 但仍有一些明显的噪声; 文献[19]的方法由于引入局部信息, 具有一定的抗噪性. 由于在隶属度域和标号域进行了模糊加权滤波, 本文方法在视觉上提供了比对比算法更好的分割结果.

      表3为利用WV和JM对图5(a1)-(d1)进行定量评价的结果. 从表中可以看出, 本文方法均具有较小的WV和JM 值, 说明本文方法相比对比算法, 分割结果的同质区域内部具有较低的方差, 且各同质区域面积较为完整, 相邻同质区域间具有较大公共边界, 满足图像分割结果对空间连续且光谱均匀的要求.

      表 3  全色遥感图像分割质量评价指标

      Table 3.  Quality evaluation of panchromatic remote sensing image segmentation

      图像 MV JM
      Kmeans FCM 文献[19] 本文方法 Kmeans FCM 文献[19] 本文方法
      图5(a1) 1.626 1.335 0.973 0.632 0.887 0.742 0.712 0.633
      图5(b1) 2.344 1.698 1.335 0.966 0.831 0.787 0.737 0.596
      图5(c1) 1.886 1.475 1.203 1.079 0.759 0.703 0.663 0.645
      图5(d1) 0.982 0.875 0.619 0.512 0.692 0.640 0.598 0.582

      由于场景的不同, 属于每一类的像素不尽相同, 存在稀疏类、密集类和均匀类[27]. 为了便于分析, 本文假设各场景为均匀类, 其像素在各类别中均匀分布. 基于此, 在表4中给出了对比算法和本文方法的计算复杂度, 其中 $ M\times N $ 是图像的像素数, $ K $ 是类别的数量, $ t $ 是迭代数, $ \omega $ 是滤波窗口的大小, $ M\times N/K $ 为每一类包含的像素数.

      表 4  计算复杂度对比

      Table 4.  Computational complexity comparison

      方法 计算复杂度
      Kmeans $O((K+M \times N / K) \times t)$
      FCM $O((M \times N \times K \times t)$
      文献[19] $O\left((1+M \times N) \times K \times t \times \omega^{2}\right)$
      本文方法 $O\left(M \times N \times K+2 \times M \times N \times \omega^{2}\right)$

      根据表4, 由于每个类包含 $ M\times N/K $ 个元素, Kmeans在每次迭代中最多进行K次比较就可以确定适当的聚类, 并进行 $ M\times N/K $ 比较以确定适当的元素, 具有低计算复杂度, 文献[19]的方法与迭代次数成线性关系且每次迭代时都需要计算像素邻域值, 具有较高的计算 复杂度, 本文方法考虑空间邻域信息时利用在隶属度和标号域滤波, 降低了计算复杂度. 此外, 为了估计算法的实用性, 我们比较了这些算法的运行时间, 如表5所示(以秒为单位).

      表 5  全色图像分割时间对比(s)

      Table 5.  Panchromatic images segmentation time comparison (s)

      方法 图5(a1) 图5(b1) 图5(c1) 图5(d1)
      Kmeans 1.49 1.62 4.17 2.38
      FCM 2.63 2.87 7.46 4.28
      文献[19] 38.04 41.66 110.31 60.48
      本文方法 1.66 1.80 3.87 2.29

      表5可以看出, 由于Kmeans和FCM仅考虑像素光谱测度对分割结果的影响, 忽略了局部空间信息, 能够在较短的时间内对图像进行分割. 文献[19]在聚类的迭代过程中利用空间邻域信息, 需要在每次迭代中计算像素邻域与聚类中心之间的距离, 使其具有较高的计算复杂度, 耗时较长且随图像尺度增大呈现非线性增长, 相比之下, 本文方法利用图像的全局信息初始化类中心, 耗时短, 而在考虑局部空间信息时, 本文在隶属度域和标号域使用滤波方法, 降低了计算复杂度, 使本文方法整体计算时间较短.

    • 为了体现本文方法的适用性, 对多光谱遥感图像进行了分割. 图6(a1)-(d1)显示了0.5米分辨率的WorldView-2图像, 其含有不同土地利用和覆盖类型. 其中, 图6(a1)-(d1)的尺度分别为 $ 476\times 432 $ $ 397\times 356 $ $ 445\times 342 $ $ 529 \times 435 $ 像素, 目视判读分别含有5、6、6和7 个同质区域. 利用本文方法和对比算法对其进行分割实验, 其中图6(a2)-(d2)为利用Kmeans方法所得的结果; 图6(a3)-(d3)为利用FCM 方法所得的结果; 图6(a4)-(d4)为利用文献[19]方法 所得的结果; 图6(a5)-(d5)为利用本文方法所得的结果.

      图  6  多光谱遥感图像和分割结果

      Figure 6.  Multispectral remote sensing images and segmentation results

      图6可以看出, 由于多光谱遥感图像光谱信息丰富, 具有更多的地物色彩, 增加了分割的难度, Kmeans 算法不能对同质区域内部的异质性进行有效分离. FCM由于考虑了模糊隶属关系而能区分部分内部区域, 但没有克服其对噪声的敏感性, 分割仍然不能令人满意. 文献[19]利用空间邻域信息, 计算每个迭代中像素邻域与聚类中心之间的距离, 提高了图像分割效果, 但固定的空间距离对不同局部信息的图像不具有鲁棒性, 使得如图6(b4)的耕地内部出现误分割现象. 本文方法所采用的隶属度空间滤波和相应的标号域滤波之间的相互作用可以增强同一邻域系统中像素之间的联系, 可以在视觉上更好地捕获同质区域, 因此分割结果比文献[19]方法更令人满意.

      表6为利用MV、JM和E对图6中的4幅图像进行评价的结果. 从表6可以看出, 本文方法的指标均是最优的, 说明本文方法可有效分割出同质区域, 且在各区域内部具有较高的纯度, 满足人眼对各区域分离的色差.

      表 6  多光谱遥感图像分割质量评价

      Table 6.  Quality evaluation of multispectral remote sensing image segmentation

      指标 方法 图6(a1) 图6(b1) 图6(c1) 图6(d1)
      MV Kmeans 1.971 1.613 2.316 2.146
      FCM 1.813 1.404 1.833 1.799
      文献[19] 1.570 1.071 1.279 1.344
      本文方法 1.376 0.796 0.941 1.001
      JM Kmeans 0.832 0.797 0.774 0.808
      FCM 0.748 0.624 0.647 0.734
      文献[19] 0.662 0.588 0.541 0.631
      本文方法 0.575 0.534 0.532 0.565
      E Kmeans 0.671 0.572 0.607 0.632
      FCM 0.524 0.466 0.573 0.597
      文献[19] 0.456 0.403 0.434 0.463
      本文方法 0.347 0.332 0.293 0.306

      表7为本文方法与对比算法对图6中的4幅图像的分割时间, 并将其与表5的全色遥感图像运行时间进行对比.

      表 7  多光谱图像分割计算时间对比(s)

      Table 7.  Multispectral images segmentation time comparison (s)

      方法 图6(a1) 图6(b1) 图6(c1) 图6(d1)
      Kmeans 4.15 2.75 3.01 5.39
      FCM 7.33 4.99 5.34 9.54
      文献[19] 106.02 66.65 77.49 138.51
      本文方法 3.56 3.15 3.37 4.68

      与全色遥感图像分割的运行时间对比可以看出, 由于多光谱遥感图像数据维数和各均匀区域中的内部异质性的增加, 使得计算复杂度增大, 各方法对同尺度图像的分割时间增加, 较之不考虑局部空间信息的Kmeans和FCM, 文献[19]的耗时增加更明显, 而本文方法是在隶属度域和标号域考虑局部空间信息, 因此耗时短且鲁棒性强, 具有实用性.

    • (1)提出一种利用图像全局信息进行类别及其中心初始化的方法, 以图像中待定类像素点的光谱一阶矩和二阶矩为初始类中心和自适应阈值, 利用二分法原理和区域相似度准则对类别数及其中心粗略估计.

      (2)以估计的类别数和类中心为输入, 通过定义的岭形模糊隶属函数快速得到图像中各像素点对于不同类别的隶属程度, 将图像域转化为隶属度域.

      (3)考虑到隶属度和标号场是各像素属性的映射, 首先在隶属度域中定义一个模糊加权均值滤波器来过滤由于噪声引起的异常隶属度值, 对滤波后的隶属度集按照最大隶属原则确定图像的标号场. 然后由移动窗口内的中值构建模糊加权中值滤波器, 以消除标号场中由于误判的类别标号.

      (4)在隶属度域和标号域滤波时, 由于原图像分辨率、场景复杂度等的不同, 造成滤波窗口邻域大小难以自动获取, 其产生不同的平滑度, 对最终分割结果中同质区域的边界有较大影响.

      综上, 本文利用图像全局信息进行类别及其中心初始化可避免由于随机选择初始值而降低类别查找的动态性, 后经岭形模糊隶属函数可快速对像素类属划分, 同时由隶属度域和标号域滤波, 充分利用了图像的局部空间信息以对类别数进行精确划定, 为可见光遥感图像的阈值分割问题提供了一个比较实用的解决方案.

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回