第13卷第2期 智能系统学报 Vol.13 No.2 2018年4月 CAAI Transactions on Intelligent Systems Apr.2018 D0:10.11992/tis.201609012 网络出版地址：http:/kns.cnki.net/cms/detail/23.1538.TP.20170626.1740.016.html 直线截距直方图城区遥感图像多阈值分割 吴诗婚，吴一全25，周建江 (1.南京航空航天大学电子信息工程学院，江苏南京211106,2.城市空间信息工程北京市重点实验室，北京100038： 3.江西省数宇国土重点实验室，江西南昌330013,4.江苏省大数据分析技术重点实验室，江苏南京210044：5.浙江 省信号处理重点实验室，浙江杭州310023) 摘要：阈值分割简单有效，但现有的单阈值方法对城区图像分割效果不佳，难以取得令人满意的结果。为了快速准 确地对城区遥感图像进行分割，本文提出了基于直线截距直方图倒数灰度嫡和人工蜂群优化(artificial bee colony op- timization,.ABC)的多阈值分割方法。首先，给出直线截距直方图的定义并建立城区遥感图像的直线截距直方图：然 后，计算该直方图倒数灰度嫡的大小，推导出其单阈值选取公式：最后，将其推广到多阈值选取，并利用人工蜂群优化 算法，对多个阈值进行快速精确地寻优，以此最终实现城区遥感图像的多阈值分割。实验结果表明，该方法所分割的 图像中多目标的形状、边缘更为准确，纹理及细节特征更加清晰，且所需运行时间仅为同类多阈值分割方法的25% 是一种行之有效的城区遥感图像分割方法。 关键词：城区提取；遥感图像；图像分割：阈值化：多阈值选取；直线截距直方图；倒数灰度嫡；人工蜂群优化 中图分类号：TP751.1:P237文献标志码：A文章编号：1673-4785(2018)02-0227-09 中文引用格式：吴诗婳，吴一全，周建江.直线截距直方图城区遥感图像多阈值分割.智能系统学报，2018.13(2)：227-235 英文引用格式：WU Shihua,WU Yiquan,ZHOU Jianjiang.Multi--level thresholding for remote sensing image of urban area based on line intercept histogramJl.CAAI transactions on intelligent systems,2018,13(2):227-235. Multi-level thresholding for remote sensing image of urban area based on line intercept histogram WU Shihua',WU Yiquan245,ZHOU Jianjiang' (1.College of Electronic and Information Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 211106,China; 2.Beijing Key Laboratory of Urban Spatial Information Engineering,Beijing 100038,China;3.Jiangxi Province Key Laboratory for Digital Land,Nanchang 330013,China;4.Jiangsu Key Laboratory of Big Data Analysis Technology,Nanjing 210044,China;5.Zhejiang Province Key Laboratory for Signal Processing,Hangzhou 310023,China) Abstract:Threshold segmentation is a kind of simple and effective method,however,the existing single-threshold method is hard to realize satisfactory effect in segmenting the images of urban area.In order to segment the remote sens- ing images of urban area quickly and accurately,a multi-threshold segmentation method based on straight-line intercept histogram,reciprocal grayscale entropy and Artificial Bee Colony(ABC)Optimization was proposed in the paper. Firstly,the straight-line intercept histogram was defined and the straight-line intercept histogram of the urban remote sensing image was established;then the value of the reciprocal grayscale entropy of the histogram was calculated and the single-threshold selection formula was deduced;finally,the application was popularized to multi-threshold selection, ABC Optimization algorithm was utilized for precise optimization of many thresholds,so as to finally realize the multi- threshold segmentation of urban remote sensing images.A large number of experiments show that,the multi-object shape and edge in the images segmented by the method are more accurate,the textures and details are more explicit,in addition,its running time is only 25%of other similar multi-threshold segmentation methods.This is a kind of effective method for segmenting the remote sensing images of urban area. Keywords:extraction of urban area;remote sensing image;image segmentation;thresholding,multi-level threshold se- lection;straight-line intercept histogram;reciprocal grayscale entropy;optimization of artificial bee colony 收稿日期：2016-09-28.网络出版日期：2017-06-26. 随着遥感技术的飞速发展，利用卫星遥感和飞 基金项目：国家自然科学基金项目(61573183)：城市空间信息工程 机遥感等方式实时获取的地物图像质量也越来越 北京市重点实验室开放基金项目(2014203)：江西省数 字国土重点实验室开放基金项目(DLLJ201412):江苏 高。从人工地物的遥感图像中提取信息，可以避免 省大数据分析技术重点实验室开放基金项目 (KXK1403):浙江省信号处理重点实验室开放基金项目 传统的实地勘测，大大提高工作效率。城区作为遥 (ZJKL6SP.OP2014-02):江苏高校优势学科建设工程 感图像中一类重要的人工地物目标，其自动提取在 项目(2012). 通信作者：吴一全.E-mail:nuaaimage@163.com. 城市规划、地理信息系统更新、数字化城市以及军
DOI: 10.11992/tis.201609012 网络出版地址: http://kns.cnki.net/kcms/detail/23.1538.TP.20170626.1740.016.html 直线截距直方图城区遥感图像多阈值分割 吴诗婳1 ，吴一全1,2,3,4,5，周建江1 （1. 南京航空航天大学 电子信息工程学院，江苏 南京 211106; 2. 城市空间信息工程北京市重点实验室，北京 100038; 3. 江西省数字国土重点实验室，江西 南昌 330013; 4. 江苏省大数据分析技术重点实验室, 江苏 南京 210044; 5. 浙江 省信号处理重点实验室, 浙江 杭州 310023） 摘 要：阈值分割简单有效，但现有的单阈值方法对城区图像分割效果不佳，难以取得令人满意的结果。为了快速准 确地对城区遥感图像进行分割，本文提出了基于直线截距直方图倒数灰度熵和人工蜂群优化 (artificial bee colony optimization, ABC) 的多阈值分割方法。首先，给出直线截距直方图的定义并建立城区遥感图像的直线截距直方图；然 后，计算该直方图倒数灰度熵的大小，推导出其单阈值选取公式；最后，将其推广到多阈值选取，并利用人工蜂群优化 算法，对多个阈值进行快速精确地寻优，以此最终实现城区遥感图像的多阈值分割。实验结果表明，该方法所分割的 图像中多目标的形状、边缘更为准确，纹理及细节特征更加清晰，且所需运行时间仅为同类多阈值分割方法的 25%， 是一种行之有效的城区遥感图像分割方法。 关键词：城区提取；遥感图像；图像分割；阈值化；多阈值选取；直线截距直方图；倒数灰度熵；人工蜂群优化 中图分类号：TP751.1；P237 文献标志码：A 文章编号：1673−4785(2018)02−0227−09 中文引用格式：吴诗婳, 吴一全, 周建江. 直线截距直方图城区遥感图像多阈值分割[J]. 智能系统学报, 2018, 13(2): 227–235. 英文引用格式：WU Shihua, WU Yiquan, ZHOU Jianjiang. Multi-level thresholding for remote sensing image of urban area based on line intercept histogram[J]. CAAI transactions on intelligent systems, 2018, 13(2): 227–235. Multi-level thresholding for remote sensing image of urban area based on line intercept histogram WU Shihua1 ，WU Yiquan1,2,3,4,5 ，ZHOU Jianjiang1 (1. College of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China; 2. Beijing Key Laboratory of Urban Spatial Information Engineering, Beijing 100038, China; 3. Jiangxi Province Key Laboratory for Digital Land, Nanchang 330013, China; 4. Jiangsu Key Laboratory of Big Data Analysis Technology, Nanjing 210044, China; 5. Zhejiang Province Key Laboratory for Signal Processing, Hangzhou 310023, China) Abstract: Threshold segmentation is a kind of simple and effective method, however, the existing single-threshold method is hard to realize satisfactory effect in segmenting the images of urban area. In order to segment the remote sensing images of urban area quickly and accurately, a multi-threshold segmentation method based on straight-line intercept histogram, reciprocal grayscale entropy and Artificial Bee Colony (ABC) Optimization was proposed in the paper. Firstly, the straight-line intercept histogram was defined and the straight-line intercept histogram of the urban remote sensing image was established; then the value of the reciprocal grayscale entropy of the histogram was calculated and the single-threshold selection formula was deduced; finally, the application was popularized to multi-threshold selection, ABC Optimization algorithm was utilized for precise optimization of many thresholds, so as to finally realize the multithreshold segmentation of urban remote sensing images. A large number of experiments show that, the multi-object shape and edge in the images segmented by the method are more accurate, the textures and details are more explicit, in addition, its running time is only 25% of other similar multi-threshold segmentation methods. This is a kind of effective method for segmenting the remote sensing images of urban area. Keywords: extraction of urban area; remote sensing image; image segmentation; thresholding; multi-level threshold selection; straight-line intercept histogram; reciprocal grayscale entropy; optimization of artificial bee colony 随着遥感技术的飞速发展，利用卫星遥感和飞 机遥感等方式实时获取的地物图像质量也越来越 高。从人工地物的遥感图像中提取信息，可以避免 传统的实地勘测，大大提高工作效率。城区作为遥 感图像中一类重要的人工地物目标，其自动提取在 城市规划、地理信息系统更新、数字化城市以及军 收稿日期：2016−09−28. 网络出版日期：2017−06−26. 基金项目：国家自然科学基金项目 (61573183)；城市空间信息工程 北京市重点实验室开放基金项目 (2014203)；江西省数 字国土重点实验室开放基金项目 (DLLJ201412)；江苏 省大数据分析技术重点实验室开放基金项目 (KXK1403)；浙江省信号处理重点实验室开放基金项目 (ZJKL_6_SP-OP2014-02)；江苏高校优势学科建设工程 项目 (2012). 通信作者：吴一全. E-mail：nuaaimage@163.com. 第 13 卷第 2 期 智 能 系 统 学 报 Vol.13 No.2 2018 年 4 月 CAAI Transactions on Intelligent Systems Apr. 2018
·228· 智能系统学报 第13卷 事侦察等实际领域均有重要的应用)。遥感图像 ization,PSO)操作简单、搜索速度较快，但在搜索过 分割是遥感图像处理中最为基础和关键的环节之 程中易陷入局部极值，且在进化后期阶段收敛精度 一，是对图像进行分析和识别的前提，对图像描述 较低、速度较慢。而近年来提出的人工蜂群优化 和特征测量有重要的影响。 算法(artificial bee colony optimization,ABC)具有计 阈值分割是一类常用的遥感图像分割方法，因 算简单、收敛速度快、需要调整参数少等优点，有 其简单快速和易于实现，成为研究热点56。该方法 望进一步提高阈值搜索的准确性和实时性。 主要是依据直方图的概率分布计算相应的准则函 基于上述分析，本文提出了一种基于直线截距 数，同时运用智能优化算法加快阈值搜索速度，从 直方图倒数灰度嫡和人工蜂群优化的城区遥感图像 而选取合适的阈值进行分割。其中基于嫡的方法最 多阈值分割方法。该方法依据遥感图像中各个像素 受关注，主要包含最大嫡法m、最大指数熵法和最 点的灰度级一邻域平均灰度级联合信息，建立该 小交叉嫡法例等。城区遥感图像通常由密集的建筑 图像的直线截距直方图，并将倒数灰度熵的单阈值 物群、草坪树丛、湖泊河流以及贯穿其中的道路交 选取准则运用到此直方图中，将此推广得到基于直 通网所构成。对含有灰度值位于明显不同区间的多 线截距直方图的多阈值选取公式，再利用人工蜂群 类目标的图像进行分割时，采用单阈值分割方法无 优化算法，搜索最优多阈值，以此对城区遥感图像 法取得令人满意的效果。而模糊聚类法一般是根 进行分割。针对大量城区遥感图像进行实验，并与 据特定的相似性度量方式和隶属度准则函数，实现 近年来提出的改进核模糊C均值聚类(kernel fuzzy 对多类目标的分割。但该算法针对样本量较大、 C means clustering,.KFCM分割法u、基于粒子群 内容复杂的图像，尤其是对高空间分辨率遥感图 优化的指数熵单、多阈值分割法以及本文提出的 像，所需的运行时间远远超过阈值分割法。因此， 单阈值分割法进行了对比，给出了相应的评价，验 为了提高城区遥感图像的分割精度和速度，可以考 证了本文方法的优越性。 虑将基于嫡的单阈值选取推广到多阈值的情况，以 此对城区遥感图像进行多阈值分割。 1直线截距直方图 在单阈值选取中，文献[12]提出的直方图最大 嫡阈值分割方法简单、易操作、效果较好，一经提出 对于一幅大小为M×N、灰度级数为L的图像， 就受到关注，并被推广到多阈值分割1。然而在该 fm,n)表示图像中像素点(m,n)的灰度级，f(m,m)∈ 类方法中采用的对数熵在零点处存在无定义值的问 [0,1,…,L-1]:g(m,n)表示其邻域平均灰度级，g(0m,m)∈ 题，且其计算仅仅依赖于直方图的分布情况，没有 [0,1,…,L-1]e 顾及图像中目标和背景区域内灰度的均匀性。为了 现有的二维直方图划分方式通常分为直分和斜 克服这一缺点，文献14]定义了倒数灰度嫡，并由此 分两种，如图1所示。 提出了倒数灰度熵阈值分割方法，弥补了零点处无 其中直分法依据相互垂直的两条直线，将二维 定义值及不能反映类内灰度均匀性的缺陷，能获得 直方图分为4个区域。一般情况下，假设边缘和噪 更好的分割效果。但目前基于熵的多阈值分割方法 声区域（图1(a)中阴影部分）的概率分布为零，则目 大多只利用了图像的一维灰度级信息，抗噪性能 标、背景对应图1(a)中的区域O、区域B。但因忽 差。为此，二维直方图被应用到多阈值分割上，增 略了阴影部分，分割结果的准确性有所下降。鉴于 强了方法的抗噪能力。但由于利用了图像的灰度 此，为了使图像分割更加准确，斜分法随后被提出叫 级一一邻域平均灰度级直方图，增加了算法的复杂 由于其对阈值的搜索空间仍为二维，运行速度有待 度，导致占用的内存空间剧增，运行速度大大降低， 进一步的提升。该方法，如图1(b)所示，通过在二 使得基于二维熵的多阈值分割方法无法满足实时性 维直方图中确定一条与主对角线垂直的以T为截距 的要求，适用范围受到限制。另外，若采用传统穷 的阈值直线，将直方图区域分为目标(O)和背景 举算法搜索最优多阈值，其运算量会随所需分割的 (B)两个部分。从图1中可以看出，二维直方图中 处于不同灰度区间的目标种类数量增多呈指数形式 共有2L-1条阈值直线。实际上由于斜分法划分直 增长。为了保证方法的实时性，条件迭代61、遗 方图时，阈值直线的斜率是固定的22)，只要确定该 传)、粒子群等优化算子被引入，以此对多阈值 直线的截距，该直线即被唯一确定，则可将整个二 选取方法进行加速，一定程度上缩短了所需的运行 维直方图区域向其对角线方向做投影，以此得到该 时间。其中，粒子群优化算法(particle swarm optim- 图像的直线截距直方图
事侦察等实际领域均有重要的应用[1-3]。遥感图像 分割是遥感图像处理中最为基础和关键的环节之 一，是对图像进行分析和识别的前提，对图像描述 和特征测量有重要的影响[4]。 阈值分割是一类常用的遥感图像分割方法，因 其简单快速和易于实现，成为研究热点[5-6]。该方法 主要是依据直方图的概率分布计算相应的准则函 数，同时运用智能优化算法加快阈值搜索速度，从 而选取合适的阈值进行分割。其中基于熵的方法最 受关注，主要包含最大熵法[7] 、最大指数熵法[8]和最 小交叉熵法[9]等。城区遥感图像通常由密集的建筑 物群、草坪树丛、湖泊河流以及贯穿其中的道路交 通网所构成。对含有灰度值位于明显不同区间的多 类目标的图像进行分割时，采用单阈值分割方法无 法取得令人满意的效果[10]。而模糊聚类法一般是根 据特定的相似性度量方式和隶属度准则函数，实现 对多类目标的分割[11]。但该算法针对样本量较大、 内容复杂的图像，尤其是对高空间分辨率遥感图 像，所需的运行时间远远超过阈值分割法。因此， 为了提高城区遥感图像的分割精度和速度，可以考 虑将基于熵的单阈值选取推广到多阈值的情况，以 此对城区遥感图像进行多阈值分割。 在单阈值选取中，文献[12]提出的直方图最大 熵阈值分割方法简单、易操作、效果较好，一经提出 就受到关注，并被推广到多阈值分割[13]。然而在该 类方法中采用的对数熵在零点处存在无定义值的问 题，且其计算仅仅依赖于直方图的分布情况，没有 顾及图像中目标和背景区域内灰度的均匀性。为了 克服这一缺点，文献[14]定义了倒数灰度熵，并由此 提出了倒数灰度熵阈值分割方法，弥补了零点处无 定义值及不能反映类内灰度均匀性的缺陷，能获得 更好的分割效果。但目前基于熵的多阈值分割方法 大多只利用了图像的一维灰度级信息，抗噪性能 差。为此，二维直方图被应用到多阈值分割上，增 强了方法的抗噪能力[15]。但由于利用了图像的灰度 级——邻域平均灰度级直方图，增加了算法的复杂 度，导致占用的内存空间剧增，运行速度大大降低， 使得基于二维熵的多阈值分割方法无法满足实时性 的要求，适用范围受到限制。另外，若采用传统穷 举算法搜索最优多阈值，其运算量会随所需分割的 处于不同灰度区间的目标种类数量增多呈指数形式 增长。为了保证方法的实时性，条件迭代[ 1 6 ] 、遗 传 [17] 、粒子群[18]等优化算子被引入，以此对多阈值 选取方法进行加速，一定程度上缩短了所需的运行 时间。其中，粒子群优化算法 (particle swarm optimization, PSO) 操作简单、搜索速度较快，但在搜索过 程中易陷入局部极值，且在进化后期阶段收敛精度 较低、速度较慢[19]。而近年来提出的人工蜂群优化 算法 (artificial bee colony optimization, ABC) 具有计 算简单、收敛速度快、需要调整参数少等优点[20] ，有 望进一步提高阈值搜索的准确性和实时性。 基于上述分析，本文提出了一种基于直线截距 直方图倒数灰度熵和人工蜂群优化的城区遥感图像 多阈值分割方法。该方法依据遥感图像中各个像素 点的灰度级——邻域平均灰度级联合信息，建立该 图像的直线截距直方图，并将倒数灰度熵的单阈值 选取准则运用到此直方图中，将此推广得到基于直 线截距直方图的多阈值选取公式，再利用人工蜂群 优化算法，搜索最优多阈值，以此对城区遥感图像 进行分割。针对大量城区遥感图像进行实验，并与 近年来提出的改进核模糊 C 均值聚类 (kernel fuzzy C means clustering, KFCM) 分割法[11] 、基于粒子群 优化的指数熵单、多阈值分割法[18]以及本文提出的 单阈值分割法进行了对比，给出了相应的评价，验 证了本文方法的优越性。 1 直线截距直方图 M ×N L f(m,n) (m,n) f (m,n) ∈ [0,1,···,L−1] g(m,n) g(m,n) ∈ [0,1,···,L−1] 对于一幅大小为 、灰度级数为 的图像， 表示图像中像素点 的灰度级， ； 表示其邻域平均灰度级， 。 现有的二维直方图划分方式通常分为直分和斜 分两种，如图 1 所示。 T 其中直分法依据相互垂直的两条直线，将二维 直方图分为 4 个区域。一般情况下，假设边缘和噪 声区域 (图 1(a) 中阴影部分) 的概率分布为零，则目 标、背景对应图 1(a) 中的区域 O、区域 B。但因忽 略了阴影部分，分割结果的准确性有所下降。鉴于 此，为了使图像分割更加准确，斜分法随后被提出[21] ， 由于其对阈值的搜索空间仍为二维，运行速度有待 进一步的提升。该方法，如图 1(b) 所示，通过在二 维直方图中确定一条与主对角线垂直的以 为截距 的阈值直线，将直方图区域分为目标 (O′) 和背景 (B′) 两个部分。从图 1 中可以看出，二维直方图中 共有 2L–1 条阈值直线。实际上由于斜分法划分直 方图时，阈值直线的斜率是固定的[22-23] ，只要确定该 直线的截距，该直线即被唯一确定，则可将整个二 维直方图区域向其对角线方向做投影，以此得到该 图像的直线截距直方图。 ·228· 智 能 系 统 学 报 第 13 卷
第2期 吴诗婳，等：直线截距直方图城区遥感图像多阈值分割 ·229· g (m,n) g (m.n) 0,L-1) (0,L-1) B (0.D B 0,s) 0 (0.0) G0)(L-1,0) 了m,m (0.0) (T,0)L-1,0)fm,m (a)二维直方图直分 ()二维直方图斜分 图1二维直方图直分与斜分 Fig.1 Vertical and oblique segmentation of two-dimensional histogram 令f(m,n)+g(m,n)=F(m,n),现以一幅城区遥感 图像（图2(a),1111像素×741像素）为例，建立其直 线截距直方图的基本思路如下：1)将原始城区遥 感图像进行灰度化处理，得到灰度图像f(m,n)(图 2(b);2)以fm,m)为依据求遥感图像中像素点的邻 域平均灰度级大小，得到其邻域平均灰度图像 g(m,n)(图2(d):3)对灰度图像和邻域平均灰度图像 分别进行直方图统计，得到其灰度级直方图（图 (d邻域平均灰度图像 ×10 2(c)和邻域平均灰度级直方图（图2(e):(4)将灰度 420 图像f(m,m)与其邻域平均灰度图像g(m,m)叠加，则 叠加后的F(m,n)可以看作一幅大小M×N、灰度级数 6 42 2L的灰度+邻域平均灰度图像。据此，对F(m,n)进行 0 灰度直方图统计，得到城区遥感图像的直线截距直 50 100150200250300 邻域平均灰度级 方图（图2（①）。 (e)邻域平均灰度级直方图 6*10 100200300400500600 灰度级+邻域平均灰度级 ()直线截距直方图 图2建立直线截距直方图的流程 (a)城区遥感图像 Fig.2 Process of constructing the line intercept histogram 基于上述分析可知，利用斜分法获取最优阈值 的过程可转化为在直线截距直方图的基础上，选取 合适的阈值选取准则，求解最优截距阈值的过程。 以此有望在保证分割结果准确性的同时，将阈值搜 索空间由二维转化为一维，进一步提高斜分法的运 行速度。 (b)灰度图像 2基于直线截距直方图倒数灰度熵和 ×10 人工蜂群优化的多阈值选取 0 50 100150200250300 灰度级 2.1基于直线截距直方图倒数灰度熵的多阈值选取 (c)灰度级直方图 设一幅城区遥感图像中像素(m,n)的灰度级为
(0, 0) (0, L−1) g (m, n) g (m, n) (L−1, 0) f (m, n) O B (t, 0) (0, s) (0, 0) (0, L−1) (T, 0) (L−1, 0) f (m, n) (0, T) B′ O′ (a) θ㐠ⰠపⰠܲ (b) θ㐠Ⱐప᫈ܲ 图 1 二维直方图直分与斜分 Fig. 1 Vertical and oblique segmentation of two-dimensional histogram f (m,n)+g(m,n) = F (m,n) f (m,n) f (m,n) g(m,n) f (m,n) g(m,n) F (m,n) M ×N 2L F (m,n) 令 ，现以一幅城区遥感 图像 (图 2(a), 1 111 像素×741 像素) 为例，建立其直 线截距直方图的基本思路如下：1) 将原始城区遥 感图像进行灰度化处理，得到灰度图像 (图 2(b))；2) 以 为依据求遥感图像中像素点的邻 域平均灰度级大小，得到其邻域平均灰度图像 (图 2(d))；3) 对灰度图像和邻域平均灰度图像 分别进行直方图统计，得到其灰度级直方图 (图 2(c)) 和邻域平均灰度级直方图 (图 2(e))；(4) 将灰度 图像 与其邻域平均灰度图像 叠加，则 叠加后的 可以看作一幅大小 、灰度级数 的灰度+邻域平均灰度图像。据此，对 进行 灰度直方图统计，得到城区遥感图像的直线截距直 方图 (图 2(f))。 (c) ▜Ꮢ㏓Ⱐప 0 50 100 150 200 250 300 1 2 ×104 ▜Ꮢ㏓ 䶽 䖑ᙋపࡦ (a( (b) ▜Ꮢప (d) 䗧ഋ▜Ꮢప (e) 䗧ഋ▜Ꮢ㏓Ⱐప 0 50 100 150 200 250 300 䗧ഋ▜Ꮢ㏓ 2 4 6 8 10 12 14 䶽 ×103 (f) Ⱐ㏫䌉Ⱐప 0 100 200 300 400 500 600 ▜Ꮢ㏓+䗧ഋ▜Ꮢ㏓ 2 4 6 䶽 ×103 图 2 建立直线截距直方图的流程 Fig. 2 Process of constructing the line intercept histogram 基于上述分析可知，利用斜分法获取最优阈值 的过程可转化为在直线截距直方图的基础上，选取 合适的阈值选取准则，求解最优截距阈值的过程。 以此有望在保证分割结果准确性的同时，将阈值搜 索空间由二维转化为一维，进一步提高斜分法的运 行速度。 2 基于直线截距直方图倒数灰度熵和 人工蜂群优化的多阈值选取 2.1 基于直线截距直方图倒数灰度熵的多阈值选取 设一幅城区遥感图像中像素 的灰度级为 (m,n) 第 2 期 吴诗婳，等：直线截距直方图城区遥感图像多阈值分割 ·229·
·230· 智能系统学报 第13卷 fm,n),其邻域平均灰度级为g(m,n),h(k为该图像 T2,…,TN)为 中满足F(m,m)=f(m,n)+gm,n)=k的像素点频数。 Ew(T,T2,…,TN)= 设图1(b)中阈值直线左下方的区域为目标区 h(k)k h(k)k 域，阈值直线右上方的区域为背景区域。则可以用 40,T)+k+ 14(T1+1,T)+k++ 阈值T将该城区遥感图像按灰度级F(m,n)分为目标 h(k)k h(k)k 类2.和背景类2 (Tx-i+1.T) 2Tw+1,2L-0+ (5) 2。={(m,n)F(m,n)=0,1,,T}, 2。={m,n)lF(m,n)=T+1,T+2,,2L-1) EM(T1,T2,…,Tw)的值越大，意味着城区遥感图 像中同一类区域内像素点灰度级的差值越小。当总 并设 倒数灰度熵Ew(T,T2,…,Tw)最大时，此时的W个阈值 F(m,n) (m,n)∈2。 即为最优多阈值，即基于直线截距直方图的倒数灰 F(x,y) (r.y)EQ 度熵多阈值选取公式为 F(m,n) (1) (m,n)E (TI',T2',....TN)=argmax (EM(Ti,T2..TN)!(6) 0T0 (7) E(T)取最大时，表明目标区域内像素灰度级大小差 式中：X(i=1,2,,S)表示第i个引领蜂的位置，即 异与背景区域内像素灰度级东西差异之和最小，即 对应的可能解。在评估当前食物源收益情况的同 图像分割的效果最理想，此时对应的灰度级阈值 时，每只引领蜂会观察周围的食物，并在附近随机 T即为最优阈值。 搜索另一个食物源 T'=ag8E(T)引 (3) Z=X+(X,-X) (8) 可根据此最优阈值T对城区遥感图像进行单阈 式中：ε是范围在[-1,1]的随机数，X表示第1≠)个 值分割： 引领蜂的位置。在当前食物源和随机搜索的食物源 0,f(m,n)+g(m,n)≤T 之间，引领蜂依据适应度值择优选取。 f'(m.n)= L-1,f(m,n)+g(m,n)>T' (4) 2)观察。引领蜂经过对食物源的初步探索，向 若图像中只含有单类目标或多类目标区域与其 观察蜂发出信号，信号强弱由引领蜂所在食物源的 他区域的灰度级差异较明显时，采用单阈值选取方 花蜜收益度确定，观察蜂依据信号强弱比例以概率 法是有效的。但是实际上，城区遥感图像中可能含 P,选取所要跟随的引领蜂 有植被、道路、湖泊河流、建筑区域等多类目标，且 F(X) P:= (9) 各目标的灰度值也有差异。因此，现将式(2)推广， F(X) 以N个阈值T,T2,…,Tx(0≤T1<T2…<Tw<2L-1) 将城区遥感图像分成N+1个灰度级区间。 式中：X表示第j个引领蜂的位置。 设40，T)=∑hkk,uT-+l,T,)=之hkk,q= 在所跟随引领蜂的周围，观察蜂随机搜索另一 个食物源，搜索方式同样依照式（⑦）。然后将该食物 j=T- 2,3.,N,uTv+1,2L-1)=方hk,则基于直线在 源信息传达给引领蜂，引领蜂再次依据适应度值飞 =Tv+1 到较优的那个位置。 距直方图的城区遥感图像的总倒数灰度嫡EM(T, 3)侦查。侦查蜂是引领蜂的变种，当引领蜂陷
f(m,n) g(m,n) h(k) F (m,n) = f (m,n)+g(m,n) = k ，其邻域平均灰度级为 ， 为该图像 中满足 的像素点频数。 T F (m,n) Ωo Ωb 设图 1(b) 中阈值直线左下方的区域为目标区 域，阈值直线右上方的区域为背景区域。则可以用 阈值 将该城区遥感图像按灰度级 分为目标 类 和背景类 Ωo = {(m,n)|F (m,n) = 0,1,···,T} Ωb = {(m,n)|F (m,n) = T +1,T +2,···,2L−1} ， 并设 Pm,n = ∑ F(m,n) (x,y)∈Ωo F(x, y) , (m,n) ∈ Ωo ∑ F(m,n) (x,y)∈Ωb F(x, y) , (m,n) ∈ Ωb (1) E(T) 则基于直线截距直方图的城区遥感图像的总倒数灰 度熵 为 E(T) = Eo + Eb = ∑ (m,n)∈Ωo Pm,n 1 1+ Pm,n + ∑ (m,n)∈Ωo Pm,n 1 1+ Pm,n = ∑T k=0 h(k) k k+uo(T) + 2∑L−1 k=T+1 h(k) k k+ub(T) (2) Eo Eb uo(T) = ∑T k=0 h(k)k ub(T) = 2∑L−1 k=T+1 h(k)k E(T) T ∗ 式中： 和 分别表示目标、背景区域的倒数灰度 熵， ， ，其可以反映目 标、背景区域的类内灰度均匀性。当总倒数灰度熵 取最大时，表明目标区域内像素灰度级大小差 异与背景区域内像素灰度级东西差异之和最小，即 图像分割的效果最理想，此时对应的灰度级阈值 即为最优阈值。 T ∗ = arg max 0⩽T⩽2L−1 {E(T)} (3) T 可根据此最优阈值 对城区遥感图像进行单阈 ∗ 值分割： f ′ (m,n) = { 0, f (m,n)+g(m,n) ⩽ T ∗ L−1, f (m,n)+g(m,n) > T ∗ (4) N T1,T2,···,TN(0 ⩽ T1 0 (7) 式中： Xi(i = 1,2,···,S L) 表示第 i 个引领蜂的位置，即 对应的可能解。在评估当前食物源收益情况的同 时，每只引领蜂会观察周围的食物，并在附近随机 搜索另一个食物源 Zi = Xi +ε(Xi − Xl) (8) ε [−1, 1] Xl 式中: 是范围在 的随机数， 表示第 个 l(l , i) 引领蜂的位置。在当前食物源和随机搜索的食物源 之间，引领蜂依据适应度值择优选取。 2) 观察。引领蜂经过对食物源的初步探索，向 观察蜂发出信号，信号强弱由引领蜂所在食物源的 花蜜收益度确定，观察蜂依据信号强弱比例以概率 Pi 选取所要跟随的引领蜂 Pi = F(Xi) ∑S L j=1 F(Xj) (9) Xj 式中： 表示第 个引领蜂的位置。 j 在所跟随引领蜂的周围，观察蜂随机搜索另一 个食物源，搜索方式同样依照式 (7)。然后将该食物 源信息传达给引领蜂，引领蜂再次依据适应度值飞 到较优的那个位置。 3) 侦查。侦查蜂是引领蜂的变种，当引领蜂陷 ·230· 智 能 系 统 学 报 第 13 卷
第2期 吴诗婳，等：直线截距直方图城区遥感图像多阈值分割 ·231· 入局部极值时，该引领蜂将会变为侦查蜂，对新的 探查一个新的食物源。再比较两者的适应度值，将 位置进行搜索，以跳出该局部最优解。 适应度值高的位置赋值给X。 通过上述3个部分的循环迭代，可搜索到全局 4)如果引领蜂在C,次循环中食物源位置不变， 最优解，即基于直线截距直方图的倒数灰度熵最优 则该引领蜂变为侦查蜂，搜索新的解取代X。 多阈值。 5)在一次循环迭代结束时，该循环的最优解将 会被记录下来，循环次数加1。 3方法步骤与流程图 6)当循环次数达到CM时，运算停止，所得解即 本文提出的基于直线截距直方图倒数灰度熵和 为本文方法的最优多阈值，以此对城区遥感图像进 人工蜂群优化的城区遥感图像多阈值分割方法的基 行分割，否则跳到3)： 本思路为：根据图像像素点的二维信息，建立一维 4实验结果与分析 直线截距直方图，并通过人工蜂群算法对最大倒数 灰度嫡多阈值选取进行优化，从而完成城区遥感图 利用本文提出的基于直线截距直方图的倒数灰 像的多阈值分割。本文方法的流程图如图3所示。 度嫡和人工蜂群优化的多阈值分割方法对大量城区 遥感图像进行了实验，并与文献[11]提出的改进 待分割城区遥感图像 KFCM聚类分割法、文献[18]提出的基于PSO的指 数熵单阈值分割法、基于PSO的指数熵多阈值分割 按式(7初始化引领蜂位置 法、本文提出的单阈值分割法在分割效果和运行时 间等方面进行了比较。实验中，KFCM聚类数目为 按式(5)计算每个引领蜂 位置的适应度函数 3:PS0中的粒子个数为30，最大迭代次数为100： 该引领蜂变身为 ABC中的蜂群种群规模为10，最大迭代次数为10， 按式(8)更新每个引领 侦查蜂 局部循环最大次数为3。实验运行环境为Intel(R) 蜂位置 Core(TM)i5CPU2.0GHz、2 GB RAM、MATLAB 观察蜂在引领蜂周围搜索， R2013a。因篇幅限制，现以其中大小不一、特点不 并更新引领蜂位置 同的3幅城区遥感图像为例说明本文方法的性能， 分别为城区遥感图像1(219像素×221像素)、城区 引领蜂位置 <不变化的选代次数 遥感图像2(249像素×305像素)、城区遥感图像3 小于C (2499像素×2393像素)。图4~6中的(a)和（①）为 分别为这3幅城区遥感图像的原始图像及其直线截 距直方图，(b)~(h)图分别为采用KFCM、基于 总迭代次数达到个 Cu PSO的指数熵单阈值、基于PSO的指数熵多阈值 (双阈值、三阈值)、本文方法单阈值、本文方法多阈 输出最优多阀值，并以此 值（双阈值、三阈值）等5种方法所得到的3幅遥感 分制城区遥感图像 图像的分割结果。 图3本文方法流程图 城区遥感图像1主要含有建筑物、河流、道路 Fig.3 Flowchart of proposed method 等多类目标，基于嫡的单阈值分割法中，本文提出 1)输入待分割城区遥感图像，设置蜂群的种群 的单阈值法的效果较好，然而利用单一的阈值难以 规模为10（引领蜂和观察蜂分别占据一半）。ABC 有效地将多类目标区分出来。基于P$O的指数嫡 算法的最大迭代次数C,为10，判断局部极值的条件 单阈值分割法所得结果的误分割率很高，其中的建 为当前循环次数C≥3。将引领蜂位置X,(i=1,2,…, 筑物、道路连成一片，完全丢失各自的特征，且存在 5)初始化，并依据式(7)计算每个位置的适应度函 严重的欠分割，目标背景均无法辨识。KFCM聚类 数值。 分割法和基于PSO的指数嫡多阈值分割法在一定 2)引领蜂按照式(8)更新位置Z,若Z的适应度 程度上改善了分割效果，但是相比本文提出的多阈 优于当前X的适应度，则将Z替代X。 值分割法，其结果中目标边缘粗糙，河流区域灰度 3)引领蜂发出的信号，观察蜂通过信号择优选 不均匀，且含有部分噪声斑。 择部分引领蜂进行跟随，并在所选引领蜂周边随机 城区遥感图像2反映的是建筑物群呈块状分布
入局部极值时，该引领蜂将会变为侦查蜂，对新的 位置进行搜索，以跳出该局部最优解。 通过上述 3 个部分的循环迭代，可搜索到全局 最优解，即基于直线截距直方图的倒数灰度熵最优 多阈值。 3 方法步骤与流程图 本文提出的基于直线截距直方图倒数灰度熵和 人工蜂群优化的城区遥感图像多阈值分割方法的基 本思路为：根据图像像素点的二维信息，建立一维 直线截距直方图，并通过人工蜂群算法对最大倒数 灰度熵多阈值选取进行优化，从而完成城区遥感图 像的多阈值分割。本文方法的流程图如图 3 所示。 Y 䖑ᙋపࡦޞܲᒱ 䒿ܦᰬфๆ䬴ը, ᎢБₐ 䖑ᙋపࡦޞܲ ᠵᐻ(7)݉ࡂᑁ䶲㰮ѹ㒚 ᕧ䔙Џ⁍䓪ݜ CM ᑁ䶲㰮ѹ㒚 ̹ऄࡂ⮰䔙Џ⁍ ᄻκCL ᠵᐻ(8)ᰠ͖ᑁ䶲 㰮ѹ㒚 䄑ᑁ䶲㰮ऄ䏗ͦ Ӓᴑ㰮 ᠵᐻ(5)䃍ッ͖ᑁ䶲㰮 ѹ㒚⮰䔮ᏀᏒܩ 㻮ᄋ㰮ᑁ䶲㰮ঔఠ᥈㉎, Ꭲᰠᑁ䶲㰮ѹ㒚 Y N N 图 3 本文方法流程图 Fig. 3 Flowchart of proposed method CM CL ⩾ 3 Xi(i = 1,2,···, 5) 1) 输入待分割城区遥感图像，设置蜂群的种群 规模为 10 (引领蜂和观察蜂分别占据一半)。ABC 算法的最大迭代次数 为 10，判断局部极值的条件 为当前循环次数 。将引领蜂位置 初始化，并依据式 (7) 计算每个位置的适应度函 数值。 Zi Zi Xi Zi Xi 2) 引领蜂按照式 (8) 更新位置 ，若 的适应度 优于当前 的适应度，则将 替代 。 3) 引领蜂发出的信号，观察蜂通过信号择优选 择部分引领蜂进行跟随，并在所选引领蜂周边随机 Xi 探查一个新的食物源。再比较两者的适应度值，将 适应度值高的位置赋值给 。 CL Xi 4) 如果引领蜂在 次循环中食物源位置不变， 则该引领蜂变为侦查蜂，搜索新的解取代 。 5) 在一次循环迭代结束时，该循环的最优解将 会被记录下来，循环次数加 1。 6) 当循环次数达到 时，运算停止，所得解即 CM 为本文方法的最优多阈值，以此对城区遥感图像进 行分割，否则跳到 3)。 4 实验结果与分析 利用本文提出的基于直线截距直方图的倒数灰 度熵和人工蜂群优化的多阈值分割方法对大量城区 遥感图像进行了实验，并与文献[11]提出的改进 KFCM 聚类分割法、文献[18]提出的基于 PSO 的指 数熵单阈值分割法、基于 PSO 的指数熵多阈值分割 法、本文提出的单阈值分割法在分割效果和运行时 间等方面进行了比较。实验中，KFCM 聚类数目为 3；PSO 中的粒子个数为 30，最大迭代次数为 100； ABC 中的蜂群种群规模为 10，最大迭代次数为 10， 局部循环最大次数为 3。实验运行环境为 Intel(R) Core(TM) i5 CPU 2.0 GHz、2GB RAM、MATLAB R2013a。因篇幅限制，现以其中大小不一、特点不 同的 3 幅城区遥感图像为例说明本文方法的性能， 分别为城区遥感图像 1 (219 像素×221 像素)、城区 遥感图像 2 (249 像素×305 像素)、城区遥感图像 3 (2 499 像素×2 393 像素)。图 4～6 中的 (a) 和 (i) 为 分别为这 3 幅城区遥感图像的原始图像及其直线截 距直方图，(b)～(h) 图分别为采用 KFCM、基于 PSO 的指数熵单阈值、基于 PSO 的指数熵多阈值 (双阈值、三阈值)、本文方法单阈值、本文方法多阈 值 (双阈值、三阈值) 等 5 种方法所得到的 3 幅遥感 图像的分割结果。 城区遥感图像 1 主要含有建筑物、河流、道路 等多类目标，基于熵的单阈值分割法中，本文提出 的单阈值法的效果较好，然而利用单一的阈值难以 有效地将多类目标区分出来。基于 PSO 的指数熵 单阈值分割法所得结果的误分割率很高，其中的建 筑物、道路连成一片，完全丢失各自的特征，且存在 严重的欠分割，目标背景均无法辨识。KFCM 聚类 分割法和基于 PSO 的指数熵多阈值分割法在一定 程度上改善了分割效果，但是相比本文提出的多阈 值分割法，其结果中目标边缘粗糙，河流区域灰度 不均匀，且含有部分噪声斑。 城区遥感图像 2 反映的是建筑物群呈块状分布 第 2 期 吴诗婳，等：直线截距直方图城区遥感图像多阈值分割 ·231·
·232· 智能系统学报 第13卷 且相对集中，其余部分较为空旷的城市区域，道路、 分割法时，背景地面的边界轮廓不清晰，对位于图 草坪等多类目标的灰度值较为相近。基于PSO的 像中间的树丛区域分割不准确，含有大量阴影。而 指数嫡单阈值和双阈值分割法，能较为准确地将建 本文提出的基于直线截距直方图倒数灰度嫡和人工 筑物与地面区分开，但却无法提取道路、草坪等其 蜂群优化的三阈值分割法则能较为准确地提取建筑 他多类目标，分割后图像的大量细节信息丢失，这 物的边界形状，草坪、道路网和树丛的轮廓均清晰 势必影响后续的地物目标检测。本文提出的单阈值 可辨，且树丛的纹理、边缘和细节特征更为丰富。 和多阈值选取方法获得了很好的分割效果，其中， 这是由于倒数灰度嫡同时考虑了图像中目标区域和 相比于单阈值方法，本文提出的多阈值选取方法的 背景区域内的灰度均匀性，对纹理丰富的区域分割 分割结果中，多类目标边缘完整、特征鲜明。KFCM 效果更优。 聚类分割法和基于PSO的指数熵三阈值分割法的 分割结果受噪声干扰，有少许虚警点，且存在目标 模糊和边缘残缺的现象，如图中右下角建筑物的阴 影信息被湮没。 (a)原始图像 (b)KFCM (a)原始图像 (b)KFCM (c)指数嫡单國 (d)指数嫡双阈 (e)指数嫡三阈 值+PSO 值+PSO 值+PSO (©)指数嫡单阈 (d指数嫡双國 (e)指数嫡三 值+PSO 值+PSO 值+PSO (①本文方法 (g)本文方法 h)本文方法 单阈值 双圆值 三阈值 (本文方法 (g)本文方法 (h)本文方法 800 单阈值 双阈值 三阈值 600 600 兹400 500 200 400 300 050100150200250300350400450500 200 灰度级 100 ()直线截距直方图 0 50100150200250300350400450500 灰度级 图5城区遥感图像2、分割结果及直线截距直方图 (①直线截距直方图 Fig.5 Remote sensing image of urban area 2,its segmenta- tion results and its line intercept histogram 图4城区遥感图像1、分割结果及直线截距直方图 Fig.4 Remote sensing image of urban area 1,its segmenta- 下面再根据正确分割率对5种方法的分割效果 tion results and its line intercept histogram 进行定量评价。以人工分割结果为标准，这里正确 城区遥感图像3中建筑物密集，道路网交织， 分割率定义为正确分割像素数目与像素总数的比 且有多处草坪、树丛覆盖，具有较高的复杂度。由 值，正确分割像素数目是指分割后的图像与经人工 分割结果可以看出，对于此类图像而言，基于嫡的 精确分割后的图像相比存在差异的像素个数。正确 单阈值分割法虽能较为准确地将建筑物与背景分 分割率越大，表示分割结果越准确。表1所示为上 离，但却无法准确地提取其他目标的边界形状。采 述5种方法的正确分割率比较。从表中可以看出， 用KFCM聚类分割法和基于PSO的指数熵多阈值 本文方法的正确分割率最大，即分割效果最优
且相对集中，其余部分较为空旷的城市区域，道路、 草坪等多类目标的灰度值较为相近。基于 PSO 的 指数熵单阈值和双阈值分割法，能较为准确地将建 筑物与地面区分开，但却无法提取道路、草坪等其 他多类目标，分割后图像的大量细节信息丢失，这 势必影响后续的地物目标检测。本文提出的单阈值 和多阈值选取方法获得了很好的分割效果，其中， 相比于单阈值方法，本文提出的多阈值选取方法的 分割结果中，多类目标边缘完整、特征鲜明。KFCM 聚类分割法和基于 PSO 的指数熵三阈值分割法的 分割结果受噪声干扰，有少许虚警点，且存在目标 模糊和边缘残缺的现象，如图中右下角建筑物的阴 影信息被湮没。 (f) ᱘᪳∁ 䬴ըࢁ (g) ᱘᪳∁ 䬴ըࣸ (h) ᱘᪳∁ ̵䬴ը (a) ࣋ప) b) KFCM (c) ᠳ⛡ࢁ䬴 ը+PSO (d) ᠳ⛡ࣸ䬴 ը+PSO (e) ᠳ⛡̵䬴 ը+PSO 0 50 100 150 200 250 300 350 400 450 500 100 200 300 400 500 600 (i) Ⱐ㏫䌉Ⱐప 䶽 ▜Ꮢ㏓ 图 4 城区遥感图像 1、分割结果及直线截距直方图 Fig. 4 Remote sensing image of urban area 1, its segmentation results and its line intercept histogram 城区遥感图像 3 中建筑物密集，道路网交织， 且有多处草坪、树丛覆盖，具有较高的复杂度。由 分割结果可以看出，对于此类图像而言，基于熵的 单阈值分割法虽能较为准确地将建筑物与背景分 离，但却无法准确地提取其他目标的边界形状。采 用 KFCM 聚类分割法和基于 PSO 的指数熵多阈值 分割法时，背景地面的边界轮廓不清晰，对位于图 像中间的树丛区域分割不准确，含有大量阴影。而 本文提出的基于直线截距直方图倒数灰度熵和人工 蜂群优化的三阈值分割法则能较为准确地提取建筑 物的边界形状，草坪、道路网和树丛的轮廓均清晰 可辨，且树丛的纹理、边缘和细节特征更为丰富。 这是由于倒数灰度熵同时考虑了图像中目标区域和 背景区域内的灰度均匀性，对纹理丰富的区域分割 效果更优。 (c) ᠳ⛡ࢁ䬴 ը+PSO (d) ᠳ⛡ࣸ䬴 ը+PSO (e) ᠳ⛡̵䬴 ը+PSO (f) ᱘᪳∁ 䬴ըࢁ (g) ᱘᪳∁ 䬴ըࣸ (h) ᱘᪳∁ ̵䬴ը (a) ࣋ప) b) KFCM (i) Ⱐ㏫䌉Ⱐప 0 50 100 150 200 250 300 350 400 450 500 200 400 600 800 䶽 ▜Ꮢ㏓ 图 5 城区遥感图像 2、分割结果及直线截距直方图 Fig. 5 Remote sensing image of urban area 2, its segmentation results and its line intercept histogram 下面再根据正确分割率对 5 种方法的分割效果 进行定量评价。以人工分割结果为标准，这里正确 分割率定义为正确分割像素数目与像素总数的比 值，正确分割像素数目是指分割后的图像与经人工 精确分割后的图像相比存在差异的像素个数。正确 分割率越大，表示分割结果越准确。表 1 所示为上 述 5 种方法的正确分割率比较。从表中可以看出， 本文方法的正确分割率最大，即分割效果最优。 ·232· 智 能 系 统 学 报 第 13 卷
第2期 吴诗婳，等：直线截距直方图城区遥感图像多阈值分割 ·233· 分割方法与本文提出的基于直线截距直方图倒数灰 度嫡和人工蜂群优化的单阈值、多阈值分割方法的 最优阈值比较，可以看出基于PSO的指数熵单阈 (a)原始图像 (b)KFCM 值、多阈值分割方法是建立于一维灰度级直方图的 基础上，最优阈值的范围为(0,255)。而本文方法是 建立于直线截距直方图的基础上，同时考虑了灰度 级和邻域平均灰度级信息，最优阈值的范围为(0,511)。 表24种方法的最优分割阈值比较 (c)指数熵单阈 (d指数嫡双國 (e)指数嫡三阈 Table 2 Comparison of four methods in optimal thresholds 值+PSO 值+PSO 值+PSO 分割 城区遥 城区遥 城区遥 阈值 方法 感图像1 感图像2 感图像3 单阈值 59 136 110 指数 双阈值 (154.83) (189,123) (172,91) (①本文方法 (g)本文方法 )本文方法 嫡+PSO 单阈值 双阈值 三阈值 三阈值(171,116,56)(191,130,68)(194,133,72) 4.0×10 单阈值 278 155 221 3.0 本文 20 (141.200) (152.282) 方法 双阈值(210,335) 1.0 三阈值(277,364,454)(127,64,310)(131,236,318) 0 50100150200250300350400450500 表3为5种分割方法所需的运行时间对比。从 灰度级 ①直线截距直方图 表中可以看出，KFCM聚类分割法的运行时间受图 图6城区遥感图像3、分割结果及直线截距直方图 像大小影响严重，针对数据量大的图像，分割时间 Fig.6 Remote sensing image of urban area 3,its segmenta- 过长，无法满足实时性的要求。而本文提出的单阈 tion results and its line intercept histogram 值方法耗时最少，但从上述对分割效果的分析可知 表15种方法的正确分割率比较 其分割效果欠佳。本文提出的多阈值分割法的耗时 Table 1 Comparison of five methods in correct segmentation 与基于P$O的指数嫡单阈值分割法相当，而与基 rate % 于PSO的指数嫡多阈值分割法相比，仅为其运行时 城区遥 城区遥 城区遥 分割方法 阈值 间的25%。这主要是由于本文方法是基于灰度级 感图像1感图像2感图像3 邻域平均灰度级联合信息所建立的一维直线截距直 KFCM 三类 97.67 98.08 94.98 方图，将原有的阈值搜索空间由二维转化为一维， 单阈值 39.65 45.92 93.70 减小了运算量，缩短了所需运行时间，且采用 指数嫡+PSO 双阈值 89.42 41.46 92.81 ABC优化算法加速了最优多阈值的搜索过程。综 三阈值 上所述，本文提出的多阈值分割方法在进一步改善 87.22 90.51 93.33 分割效果的同时，大幅度提高了方法的运行速度。 单阈值 90.74 88.58 93.99 表35种方法的运行时间比较 本文方法 双阈值 97.25 94.13 97.08 Table 3 Comparison of five methods in running time 三阈值 98.46 98.37 98.21 城区遥 城区遥 城区遥 分割方法 阈值 上述主观视觉效果分析和客观定量评价指标两 感图像1 感图像2 感图像3 方面均表明，无论是针对灰度层次较简单的图像 KFCM 三类 1.2745 2.1660 129.7655 (如图5(a),还是针对纹理细节丰富、数据量较大的 单阈值 0.0780 0.0229 1.5354 图像（如图6(a),与KFCM聚类分割法、基于熵的 指数熵+PSO双阈值 0.3450 0.2837 1.6478 单阈值分割法、基于PSO的指数嫡多阈值分割法相 三阈值 0.4283 0.3634 1.8345 比，本文提出的基于直线截距直方图倒数灰度熵和 单阈值 0.0194 0.0169 0.3480 人工蜂群优化的多阈值分割方法在分割效果上均有 本文方法 双阈值 0.0706 0.072 0.4074 明显的优势。 三阈值 0.0823 0.0843 表2给出了基于PSO的指数嫡单阈值、多阈值 0.4738
(a) ࣋ప) b) KFCM (c) ᠳ⛡ࢁ䬴 ը+PSO (d) ᠳ⛡ࣸ䬴 ը+PSO (e) ᠳ⛡̵䬴 ը+PSO (f) ᱘᪳∁ 䬴ըࢁ (g) ᱘᪳∁ 䬴ըࣸ (h) ᱘᪳∁ ̵䬴ը ×104 (i) Ⱐ㏫䌉Ⱐప 0 50 100 150 200 250 300 350 400 450 500 1.0 2.0 3.0 4.0 䶽 ▜Ꮢ㏓ 图 6 城区遥感图像 3、分割结果及直线截距直方图 Fig. 6 Remote sensing image of urban area 3, its segmentation results and its line intercept histogram 表 1 5 种方法的正确分割率比较 Table 1 Comparison of five methods in correct segmentation rate % 分割方法 阈值 城区遥 感图像 1 城区遥 感图像 2 城区遥 感图像 3 KFCM 三类 97.67 98.08 94.98 指数熵+PSO 单阈值 39.65 45.92 93.70 双阈值 89.42 41.46 92.81 三阈值 87.22 90.51 93.33 本文方法 单阈值 90.74 88.58 93.99 双阈值 97.25 94.13 97.08 三阈值 98.46 98.37 98.21 上述主观视觉效果分析和客观定量评价指标两 方面均表明，无论是针对灰度层次较简单的图像 (如图 5(a))，还是针对纹理细节丰富、数据量较大的 图像 (如图 6(a))，与 KFCM 聚类分割法、基于熵的 单阈值分割法、基于 PSO 的指数熵多阈值分割法相 比，本文提出的基于直线截距直方图倒数灰度熵和 人工蜂群优化的多阈值分割方法在分割效果上均有 明显的优势。 表 2 给出了基于 PSO 的指数熵单阈值、多阈值 分割方法与本文提出的基于直线截距直方图倒数灰 度熵和人工蜂群优化的单阈值、多阈值分割方法的 最优阈值比较，可以看出基于 PSO 的指数熵单阈 值、多阈值分割方法是建立于一维灰度级直方图的 基础上，最优阈值的范围为 (0, 255)。而本文方法是 建立于直线截距直方图的基础上，同时考虑了灰度 级和邻域平均灰度级信息，最优阈值的范围为 (0, 511)。 表 2 4 种方法的最优分割阈值比较 Table 2 Comparison of four methods in optimal thresholds 分割 方法 阈值 城区遥 感图像 1 城区遥 感图像 2 城区遥 感图像 3 指数 熵+PSO 单阈值 59 136 110 双阈值 (154, 83) (189, 123) (172, 91) 三阈值 (171, 116, 56) (191, 130, 68) (194, 133, 72) 本文 方法 单阈值 278 155 221 双阈值 (210, 335) (141, 200) (152, 282) 三阈值 (277, 364, 454) (127, 64, 310) (131, 236, 318) 表 3 为 5 种分割方法所需的运行时间对比。从 表中可以看出，KFCM 聚类分割法的运行时间受图 像大小影响严重，针对数据量大的图像，分割时间 过长，无法满足实时性的要求。而本文提出的单阈 值方法耗时最少，但从上述对分割效果的分析可知， 其分割效果欠佳。本文提出的多阈值分割法的耗时 与基于 PSO 的指数熵单阈值分割法相当，而与基 于 PSO 的指数熵多阈值分割法相比，仅为其运行时 间的 25%。这主要是由于本文方法是基于灰度级、 邻域平均灰度级联合信息所建立的一维直线截距直 方图，将原有的阈值搜索空间由二维转化为一维， 减小了运算量，缩短了所需运行时间，且采用 ABC 优化算法加速了最优多阈值的搜索过程。综 上所述，本文提出的多阈值分割方法在进一步改善 分割效果的同时，大幅度提高了方法的运行速度。 表 3 5 种方法的运行时间比较 Table 3 Comparison of five methods in running time s 分割方法 阈值 城区遥 感图像 1 城区遥 感图像 2 城区遥 感图像 3 KFCM 三类 1.274 5 2.166 0 129.765 5 指数熵+PSO 单阈值 0.078 0 0.022 9 1.535 4 双阈值 0.345 0 0.283 7 1.647 8 三阈值 0.428 3 0.363 4 1.834 5 本文方法 单阈值 0.019 4 0.016 9 0.348 0 双阈值 0.070 6 0.072 0.407 4 三阈值 0.082 3 0.084 3 0.473 8 第 2 期 吴诗婳，等：直线截距直方图城区遥感图像多阈值分割 ·233·
·234· 智能系统学报 第13卷 5结束语 mensional Otsu image segmentation method and fast recurs- ive realization[J].Journal of electronics and information 本文提出的基于直线截距直方图倒数灰度熵 technology,2010,32(5):1100-1104. 和ABC的多阈值分割方法，以城区遥感图像像素 [7]RANJANI JJ,THIRUVENGADAM S J.Fast threshold se- 的二维联合信息为基础建立直线截距直方图，以此 lection algorithm for segmentation of synthetic aperture 将阈值搜索空间转化为一维。另一方面，以倒数灰 radar images[J].IET radar,sonar and navigation,2012,6(8): 度熵作为城区遥感图像的阈值选取准则函数，并采 788-795. 用ABC算法对最优阈值的搜寻进行优化，大大减 [8]张金矿，吴一全.基于Tent映射CPSO的二维斜分指数嫡 少了方法的运行时间。本文方法是对基于嫡理论的 阈值分割.信号处理，2010,26(5)：703-708. 城区遥感图像分割技术的进一步扩展和补充。实验 ZHANG Jinkuang,WU Yiquan.Image thresholding based 结果表明，与KFCM聚类分割法、基于PSO的指数 on 2-D oblique exponent entropy method and Tent map chaotic particle swarm algorithm[J].Signal processing, 嫡阈值分割法等进行对比，在主观视觉效果和客观 2010,26(5):703-708 评价指标两个方面，证明了本文方法的准确性和实 [9]SARKAR S,DAS S,CHAUDHURI S S.A multilevel col- 时性。本文提出的方法，已应用于地物识别中的城 or image thresholding scheme based on minimum cross en- 区、建筑物分割，取得了极佳的分割效果。 tropy and differential evolution[J].Pattern recognition let- 参考文献： ters,2015,54:27-35 [10]MALYSZKO D,STEPANIUK J.Adaptive multilevel [1]SIRMACEK B,UNSALAN C.Urban area detection using rough entropy evolutionary thresholding[J].Information local feature points and spatial voting[J].IEEE geoscience sciences,2010,180(7):1138-1158 and remote sensing letters,2010,7(1):146-150 [11]NIAZMARDI S,NAEINI AA,HOMAYOUNI S,et al. [2]SIRMACEK B,UNSALAN C.Using local features to Particle swarm optimization of kernel-based fuzzy C- measure land development in urban regions[J].Pattern re- means for hyperspectral data clustering[J].Journal of ap- cognition letters,2010,31(10):1155-1159. plied remote sensing,2012,6(1):063601. [3]朱江洪，李江风，叶菁.利用决策树工具的土地利用类型 [12]KAPUR J N,SAHOO P K,WONG A K C.A new method 遥感识别方法研究).武汉大学学报：信息科学版，2011， for gray-level picture thresholding using the entropy of the 36(3):301-305. histogram[J].Computer vision,graphics,and image pro- ZHU Jianghong,LI Jiangfeng,YE Jing.Land use informa- cessing,.1985,29(3:273-285. tion extraction from remote sensing data based on decision [13]CAO L,SHI Z,CHENG E K W.Fast automatic multilevel tree tool[J].Geomatics and information science of Wuhan thresholding method[J].Electronics letters,2002,38(16): university,2011,36(3):301-305. 868-870 「4]陈洪，陶超，邹峥嵘，等.一种新的高分辨率遥感影像城区 [14]吴一全，孟天亮，吴诗婳，等.基于二维倒数灰度嫡的河 提取方法[J】.武汉大学学报：信息科学版，2013,38(9)： 流遥感图像分割[.华中科技大学学报：自然科学版， 1063-1067 2014,42(12):70-74,80 CHEN Hong,TAO Chao,ZOU Zhengrong,et al.Automat- WU Yiquan,MENG Tianliang,WU Shihua,et al.Remote ic urban area extraction using a Gabor filter and high-resolu- sensing images segmentation of rivers based on two-di- tion remote sensing imagery[].Geomatics and information mensional reciprocal gray entropy[J].Journal of Huazhong science of Wuhan university,2013,38(9):1063-1067. university of science and technology:nature science,2014, [5]李丽，柴文婷，梅树立.基于自适应全局阈值融合标记的 42(12):70-74,80. 遥感图像建筑群分割[J].农业机械学报，2013,44(7)： [15]陈恺，陈芳，戴敏，等.基于萤火虫算法的二维嫡多阈值 222-228 快速图像分割).光学精密工程，2014,22(2：517-523. LILi,CHAI Wenting,MEI Shuli.Segmentation of remote CHEN Kai,CHEN Fang,DAI Min,et al.Fast image seg- sensing images based on adaptive global threshold and fused mentation with multilevel threshold of two-dimensional en- markers[].Transactions of the Chinese society for agricul- tropy based on firefly algorithm[J].Optics and precision tural machinery,2013.44(7):222-228. engineering,2014,22(2:517-523. [6陈琪，熊博莅，陆军，等.改进的二维Ots图像分割方法 「16]罗希平，田捷.用最大嫡原则作多阈值选择的条件迭代 及其快速实现[J.电子与信息学报，2010,32(5)： 算法.软件学报，2000,11(3：379-385. 1100-1104 LUO Xiping,TIAN Jie.The ICM algorithm for multi-level CHEN Qi,XIONG Boli,LU Jun,et al.Improved Two-Di- threshold selection by maximum entropy criterion[J]
5 结束语 本文提出的基于直线截距直方图倒数灰度熵 和 ABC 的多阈值分割方法，以城区遥感图像像素 的二维联合信息为基础建立直线截距直方图，以此 将阈值搜索空间转化为一维。另一方面，以倒数灰 度熵作为城区遥感图像的阈值选取准则函数，并采 用 ABC 算法对最优阈值的搜寻进行优化，大大减 少了方法的运行时间。本文方法是对基于熵理论的 城区遥感图像分割技术的进一步扩展和补充。实验 结果表明，与 KFCM 聚类分割法、基于 PSO 的指数 熵阈值分割法等进行对比，在主观视觉效果和客观 评价指标两个方面，证明了本文方法的准确性和实 时性。本文提出的方法，已应用于地物识别中的城 区、建筑物分割，取得了极佳的分割效果。 参考文献： SIRMACEK B, UNSALAN C. Urban area detection using local feature points and spatial voting[J]. IEEE geoscience and remote sensing letters, 2010, 7(1): 146–150. [1] SIRMAÇEK B, ÜNSALAN C. Using local features to measure land development in urban regions[J]. Pattern recognition letters, 2010, 31(10): 1155–1159. [2] 朱江洪, 李江风, 叶菁. 利用决策树工具的土地利用类型 遥感识别方法研究[J]. 武汉大学学报: 信息科学版, 2011, 36(3): 301–305. ZHU Jianghong, LI Jiangfeng, YE Jing. Land use information extraction from remote sensing data based on decision tree tool[J]. Geomatics and information science of Wuhan university, 2011, 36(3): 301–305. [3] 陈洪, 陶超, 邹峥嵘, 等. 一种新的高分辨率遥感影像城区 提取方法[J]. 武汉大学学报:信息科学版, 2013, 38(9): 1063–1067. CHEN Hong, TAO Chao, ZOU Zhengrong, et al. Automatic urban area extraction using a Gabor filter and high-resolution remote sensing imagery[J]. Geomatics and information science of Wuhan university, 2013, 38(9): 1063–1067. [4] 李丽, 柴文婷, 梅树立. 基于自适应全局阈值融合标记的 遥感图像建筑群分割[J]. 农业机械学报, 2013, 44(7): 222–228. LI Li, CHAI Wenting, MEI Shuli. Segmentation of remote sensing images based on adaptive global threshold and fused markers[J]. Transactions of the Chinese society for agricultural machinery, 2013, 44(7): 222–228. [5] 陈琪, 熊博莅, 陆军, 等. 改进的二维 Otsu 图像分割方法 及其快速实现[J]. 电子与信息学报, 2010, 32(5): 1100–1104. CHEN Qi, XIONG Boli, LU Jun, et al. Improved Two-Di- [6] mensional Otsu image segmentation method and fast recursive realization[J]. Journal of electronics and information technology, 2010, 32(5): 1100–1104. RANJANI J J, THIRUVENGADAM S J. Fast threshold selection algorithm for segmentation of synthetic aperture radar images[J]. IET radar, sonar and navigation, 2012, 6(8): 788–795. [7] 张金矿, 吴一全. 基于 Tent 映射 CPSO 的二维斜分指数熵 阈值分割[J]. 信号处理, 2010, 26(5): 703–708. ZHANG Jinkuang, WU Yiquan. Image thresholding based on 2-D oblique exponent entropy method and Tent map chaotic particle swarm algorithm[J]. Signal processing, 2010, 26(5): 703–708. [8] SARKAR S, DAS S, CHAUDHURI S S. A multilevel color image thresholding scheme based on minimum cross entropy and differential evolution[J]. Pattern recognition letters, 2015, 54: 27–35. [9] MAŁYSZKO D, STEPANIUK J. Adaptive multilevel rough entropy evolutionary thresholding[J]. Information sciences, 2010, 180(7): 1138–1158. [10] NIAZMARDI S, NAEINI A A, HOMAYOUNI S, et al. Particle swarm optimization of kernel-based fuzzy Cmeans for hyperspectral data clustering[J]. Journal of applied remote sensing, 2012, 6(1): 063601. [11] KAPUR J N, SAHOO P K, WONG A K C. A new method for gray-level picture thresholding using the entropy of the histogram[J]. Computer vision, graphics, and image processing, 1985, 29(3): 273–285. [12] CAO L, SHI Z, CHENG E K W. Fast automatic multilevel thresholding method[J]. Electronics letters, 2002, 38(16): 868–870. [13] 吴一全, 孟天亮, 吴诗婳, 等. 基于二维倒数灰度熵的河 流遥感图像分割[J]. 华中科技大学学报: 自然科学版, 2014, 42(12): 70–74, 80. WU Yiquan, MENG Tianliang, WU Shihua, et al. Remote sensing images segmentation of rivers based on two-dimensional reciprocal gray entropy[J]. Journal of Huazhong university of science and technology: nature science, 2014, 42(12): 70–74, 80. [14] 陈恺, 陈芳, 戴敏, 等. 基于萤火虫算法的二维熵多阈值 快速图像分割[J]. 光学精密工程, 2014, 22(2): 517–523. CHEN Kai, CHEN Fang, DAI Min, et al. Fast image segmentation with multilevel threshold of two-dimensional entropy based on firefly algorithm[J]. Optics and precision engineering, 2014, 22(2): 517–523. [15] 罗希平, 田捷. 用最大熵原则作多阈值选择的条件迭代 算法[J]. 软件学报, 2000, 11(3): 379–385. LUO Xiping, TIAN Jie. The ICM algorithm for multi-level threshold selection by maximum entropy criterion[J]. [16] ·234· 智 能 系 统 学 报 第 13 卷
第2期 吴诗婳，等：直线截距直方图城区遥感图像多阈值分割 ·235· Journal of software,2000,11(3):379-385 WU Shihua,WU Yiquan,ZHOU Jianjiang,et al.Segment- [1刀郑毅，郑苹.结合模糊嫡和遗传算法的双阈值图像分割 ation method based on line intercept histogram reciprocal ).应用科学学报，2014,32(4：427-433. cross entropy for medical image[J].Journal of data acquisi- ZHENG Yi,ZHENG Ping.Dual thresholding method us- tion and processing,2015,30(5):982-992. ing fuzzy entropy and genetic algorithm[J].Journal of ap- [23]吴一全，龙云淋.基于直线截距直方图的Arimoto嫡或 plied sciences,2014,32(4):427-433. Arimoto灰度嫡的食品图像分割[J].现代食品科技， [18)王树亮，赵合计.基于改进粒子群算法的多阈值灰度图 2016.32(1):164-169. 像分割计算机应用，2012,32(S2):147-150. WU Yiquan,LONG Yunlin.Food image segmentation WANG Shuliang,ZHAO Heji.Multilevel thresholding based on line intercept histogram Arimoto entropy or Ar- gray-scale image segmentation based on improved particle imoto gray entropy[J].Modern food science and techno- swarm optimization[J].Journal of computer applications, 1ogy,2016,32(1)164-169. 2012,32(S2):147-150. [19吴一全，张晓杰，吴诗婳，等.利用高速收敛PS0或分解 作者简介： 进行二维灰度嫡图像分割U,武汉大学学报：信息科学 吴诗婳，女，硕土研究生，主要研 版.2011,36(9：1059-1063. 究方向为图像处理。发表学术论文多 WU Yiquan,ZHANG Xiaojie,WU Shihua,et al.Two-di- 篇。 mensional gray entropy image thresholding based on particle swarm optimization with high speed convergence or decomposition[.Geomatics and information science of Wuhan university,2011,36(9):1059-1063. [20]HORNG M H.A multilevel image thresholding using the 吴一全，男，教授，博士生导师，博 士，主要研究方向为图像处理与分析 honey bee mating optimization[J].Applied mathematics 日标检测与识别、智能信息处理。发 and computation,2010,215(9):3302-3310. 表学术论文280余篇。 [2]吴一全，纪守新，吴诗婳，等.基于二维直分与斜分灰度 嫡的图像阈值选取[J.天津大学学报，2011,44(12)： 1043-1049 WU Yiquan,JI Shouxin,WU Shihua,et al.Gray entropy 周建江，男，教授，博土生导师，博 image thresholding based on 2-dimensional histogram ver- 士，主要研究方向雷达目标特性分析、 tical and oblique segmentation[J].Journal of Tianjin uni- 特征控制与目标识别、机载电子信息 versity,2011,44(12)1043-1049. 系统、DSP技术。 [22]吴诗婳，吴一全，周建江，等.面向医学图像分割的直线 截距直方图倒数交叉嫡方法[).数据采集与处理，2015， 30(5:982-992
Journal of software, 2000, 11(3): 379–385. 郑毅, 郑苹. 结合模糊熵和遗传算法的双阈值图像分割 [J]. 应用科学学报, 2014, 32(4): 427–433. ZHENG Yi, ZHENG Ping. Dual thresholding method using fuzzy entropy and genetic algorithm[J]. Journal of applied sciences, 2014, 32(4): 427–433. [17] 王树亮, 赵合计. 基于改进粒子群算法的多阈值灰度图 像分割[J]. 计算机应用, 2012, 32(S2): 147–150. WANG Shuliang, ZHAO Heji. Multilevel thresholding gray-scale image segmentation based on improved particle swarm optimization[J]. Journal of computer applications, 2012, 32(S2): 147–150. [18] 吴一全, 张晓杰, 吴诗婳, 等. 利用高速收敛 PSO 或分解 进行二维灰度熵图像分割[J]. 武汉大学学报：信息科学 版, 2011, 36(9): 1059–1063. WU Yiquan, ZHANG Xiaojie, WU Shihua, et al. Two-dimensional gray entropy image thresholding based on particle swarm optimization with high speed convergence or decomposition[J]. Geomatics and information science of Wuhan university, 2011, 36(9): 1059–1063. [19] HORNG M H. A multilevel image thresholding using the honey bee mating optimization[J]. Applied mathematics and computation, 2010, 215(9): 3302–3310. [20] 吴一全, 纪守新, 吴诗婳, 等. 基于二维直分与斜分灰度 熵的图像阈值选取[J]. 天津大学学报, 2011, 44(12): 1043–1049. WU Yiquan, JI Shouxin, WU Shihua, et al. Gray entropy image thresholding based on 2-dimensional histogram vertical and oblique segmentation[J]. Journal of Tianjin university, 2011, 44(12): 1043–1049. [21] 吴诗婳, 吴一全, 周建江, 等. 面向医学图像分割的直线 截距直方图倒数交叉熵方法[J]. 数据采集与处理, 2015, 30(5): 982–992. [22] WU Shihua, WU Yiquan, ZHOU Jianjiang, et al. Segmentation method based on line intercept histogram reciprocal cross entropy for medical image[J]. Journal of data acquisition and processing, 2015, 30(5): 982–992. 吴一全, 龙云淋. 基于直线截距直方图的 Arimoto 熵或 Arimoto 灰度熵的食品图像分割[J]. 现代食品科技, 2016, 32(1): 164–169. WU Yiquan, LONG Yunlin. Food image segmentation based on line intercept histogram Arimoto entropy or Arimoto gray entropy[J]. Modern food science and technology, 2016, 32(1): 164–169. [23] 作者简介： 吴诗婳，女，硕士研究生，主要研 究方向为图像处理。发表学术论文多 篇。 吴一全，男，教授，博士生导师，博 士，主要研究方向为图像处理与分析、 目标检测与识别、智能信息处理。发 表学术论文 280 余篇。 周建江，男，教授，博士生导师，博 士，主要研究方向雷达目标特性分析、 特征控制与目标识别、机载电子信息 系统、DSP 技术。 第 2 期 吴诗婳，等：直线截距直方图城区遥感图像多阈值分割 ·235·