基于改进SLIC算法的Landsat 8 OLI遥感影像水体提取.pdf

收稿日期：基金项目：安徽省重点研究与开发计划项目（）；国家自然科学基金青年项目（）；安徽省高校省级自然科学研究项目（）作者简介：方悦（），女，安徽六安人，硕士研究生，研究方向为遥感图像处理 通信作者：颜普（），男，安徽亳州人，讲师，博士，主要从事计算机视觉方面的研究工作：【水利信息化】基于改进 算法的 遥感影像水体提取方 悦，颜 普，陈 杰，（安徽建筑大学 电子与信息工程学院，安徽 合肥；安徽省国际古建筑智能感知与高维建模联合研究中心，安徽 合肥）摘 要：洪涝灾害时常发生，基于遥感影像绘制水体分布图观察水体变化很有必要。针对简单线性迭代聚类（）算法在遥感影像水体分割时的易误分割问题，提出一种改进的新超像素分割算法（）。首先基于预处理后的 遥感影像计算新水体指数，突出水体特征；再结合线性特征增强和局部三值模式（）纹理特征获得新纹理特征（），纳入 算法对 图像进行超像素分割处理；最后对水体超像素进行区域合并，得到完整水体轮廓，通过提取同一区域不同时间的水体轮廓有效观察水体变化。结果表明，与 算法相比，算法的边界召回率更高、欠分割错误率更低，解决了小水体难分割的问题。关键词：遥感影像；水体提取；超像素分割；算法；局部三值模式中图分类号：；文献标志码：引用格式：方悦，颜普，陈杰基于改进 算法的 遥感影像水体提取人民黄河，（）：，（，；，）：，（）（），（）（），：；在所有自然灾害中洪涝灾害发生频率很高，通过提取水体信息、精确绘制水体图来观察水体变化，对灾害救援和评估具有重要意义。与实地勘测相比，基于遥感影像提取水体信息可以节约大量劳动力成本、提高水体信息提取精度。近年来，随着卫星技术快速发展，遥感影像数据种类越来越多，其中 遥感影像数据可免费获得，其具有中等空间分辨率，可满足水体提取要求，因此被广泛用于水体提取。在遥感影像水体提取方法中最常用的是水体指数法，但该方法十分依赖阈值选择，需要重复大量实验才可确定最佳阈值。简单线性迭代聚类（）算法是一种高效的超像素分割算法，其将图像中的像素点聚类到超像素中，具有优异的分割性能。已有多数研究者利用 算法在遥感影像中提取水体，湛南渝等采用 算法在哨兵 号雷达数据中提取水体，并对 算法提取信息与 阈值法提取信息进行比较分析，发现 算法的提取效率与精度更高。但采用 算法进行像素相似性度量时只考虑影像颜色特征和空间特征，这对复杂遥感影像分割是不利的。遥感影像中往往有很多颜色接近的地物，如第 卷第 期 人 民 黄 河 ，年 月 ，水体一般呈现深色调，其与黑色建筑物表面颜色很接近，因此需要加入纹理特征来精确分割此类地物。等将灰度共生矩阵纹理特征纳入 算法，针对高分二号全色图像、多光谱图像和融合图像验证改进算法的有效性，结果表明，改进算法的分割完整性和正确性都高于 算法的。此外，算法常出现过分割问题，需要利用区域合并方法来解决此问题，区域邻接图（）、近邻图（）和层次聚类法是 种主要的区域合并方法。相较于区域邻接图和近邻图，层次聚类法的输入参数较少。等将 纹理特征纳入层次聚类法，验证了该方法具有良好的区域合并性能。因此，笔者也选择将纹理特征纳入层次聚类方法进行区域合并。为将水体与相似颜色地物分割开，笔者结合线性特征增强方法与局部三值模式（）纹理特征提出一种新的纹理特征（），将 纹理特征纳入 算法，得到新超像素分割方法（），之后对水体超像素进行区域合并得到完整水体，通过提取同一区域不同时间的水体，分析水体变化来检测洪涝灾害，以期为洪涝灾害的监测和预警提供更加可靠的手段。研究区域选取和数据预处理 研究区域选取选择 个区域进行研究，个区域涉及黄河及其 条主要支流（分别为位于河南省巩义市的伊洛河、陕西省渭南市的渭河、甘肃省兰州市的湟水），个区域的伪彩色图像见图。各图像都包含具有线性特征的大水体和小水体，小水体较为狭窄，增大了水体分割的难度。图 伪彩色图像（波段组合）遥感影像和 图像选用的 遥感影像数据源自地理空间数据云网站（：），数据已经过几何校正和地形校正，具体信息见表。遥感影像数据往往包含云层，为了避免云阴影导致的水体分割误差增大，裁剪出无云阴影图像作为研究区域。表 研究区域 遥感影像数据信息研究区域条带 行影像获取日期云量 巩义渭南兰州 把高分辨率 图像作为辅助信息，对 遥感影像数据的“真实”水体边界进行人工数字化标定，以获得 图像用来评估分割性能。数据预处理仅仅经过几何校正和地形校正的 遥感影像数据是数字图像，在计算水体指数之前需要先对影像数据进行辐射定标，将影像呈现的原始像元值转换为对应像元的地表反射率。因大气吸收和散射机制等因素会影响地表反射率的精确性，需对辐射定标后的数据再进行大气校正处理，大气校正处理模型是 模型。研究方法水体提取流程见图，主要环节有：计算水体指数，利用 超像素分割算法对 图像进行分割，基于 纹理特征和颜色特征对超像素进行区域合并得到完整水体。图 水体提取流程人 民 黄 河 年第 期 水体指数 归一化差异水体指数（）和修正归一化差异水体指数（）等传统水体指数对小水体的增强效果较差，因此通过（）变换、变换、自动水体提取指数（）和 变换提出一个新水体指数（）。变换可产生绿度、亮度、湿度 个分量。变换可产生地物总辐射水平（）、地物可见光红外光辐射平衡（）和辐射变化矢量（）个分量。与其他地物相比，水体的红外辐射最弱，采用 分量可以突出水体与其他地物的差异，显著反映水体特征。可以有效突出深色建筑物与水体的差异，减小误分类概率。自动水体提取指数计算公式为 （）（）（）式中：为自动水体提取指数，、分别为 遥感影像数据的波段（波段）、波段（波段）、波段（波段）、波段（波段）的反射率。将湿度分量、分量和 分别作为 图像的红色通道、绿色通道和蓝色通道，然后对 图像进行 变换，变换得到的 分量即为 图像。兰州区域、水体指数图像见图，从图上红色矩形框中可看出，相比其他 个传统水体指数，小水体与相邻地物的对比度更高，更有利于后续水体分割。图 种水体指数图像 超像素分割 纹理特征提取在 超像素分割聚类过程中度量周边邻域像素与聚类中心像素的相似性时，只考虑颜色特征和空间特征，这对复杂遥感影像的水体分割很不利，因此可在相似性度量中加入纹理特征距离，得到高度相似的超像素块。容易受噪声影响，是由 扩展的三值模式，为了提高描述符的稳健性，选择 算子来提取 图像的纹理特征。计算公式如下：（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）正 （）（）负 （）（）式中：、分别为 邻域的中心像素、周边像素，为阈值，（）为 编码 个位置的值，为 编码，（）、（）分别为正模式编码、负模式编码个位置的值，正、负分别为正模式编码、负模式编码。经过实验发现正模式编码得到的纹理特征更适合用来区分水体和非水体，因此本文选用正模式编码提取 纹理特征。纹理特征为了使 图像的纹理特征突出线性水体与非水体的差异，需要对 图像进行线性特征增强，再提取 纹理特征。线性特征增强需先使用、和 这 个方向的算子对图像进行滤波。滤波后图像的每个像素在 个方向的最大值所组成的图像就是线性特征增强图像，见图。对比兰州区域 图像的 与 纹理特征，从图上红色矩形框内可明显看出，纹理特征的水体与相邻地物的对比度较高，有利于后续水体分割。图 兰州区域 图像的两种纹理特征 算法复杂遥感影像中有很多与水体颜色相似的地物，为了解决水体与相邻地物难分割的难题，结合 纹理特征和 算法对遥感影像进行超像素分割。算法有 个输入参数，即预期的超像素个数（初始聚类中心）和超像素紧密度。在对灰度图像进行 超像素分割时，从初始化步骤开始，个初始聚类中心分别分布在间隔 个像素的规则网格上，网格间隔 ，其中 为图像的总像素数。为了避人 民 黄 河 年第 期免聚类中心处于图像边缘或者噪声上，将聚类中心移到其 邻域最低梯度位置上，之后在聚类中心周围 区域中搜索相似像素。算法中原像素相似性度量公式为（）（）（）（）（）（）（）（）（）（）（）式中：为灰度距离；为聚类中心 的灰度值；为像素 的灰度值；（，）是聚类中心 的坐标；（，）是像素 的坐标；为空间距离；为像素 与聚类中心 之间基于灰度特征和颜色特征的距离；为类内最大灰度距离，设定为输入参数紧密度；为类内最大空间距离，；为 对应的最终的相似性度量。在 算法中提出一种新相似性度量，在原来只考虑颜色特征和空间特征的基础上，加入纹理特征距离，纹理特征距离计算公式为 （）（）式中：为聚类中心 的 纹理特征，为像素 的 纹理特征，为 纹理特征距离。使用类内最大纹理特征距离将纹理特征与颜色、空间特征合并为一个度量值，新相似性度量 对应的最终相似性度量 计算公式为（）（）（）（）（）（）（）（）式中：为类内最大纹理特征距离，为像素 与聚类中心 之间基于灰度、颜色和纹理特征的距离。每个像素都与距离最小的聚类中心聚类，一遍聚类之后，将当前超像素的平均矢量更新为新的聚类中心，以此类推，重复迭代 次。通过后处理把很小的超像素分配给相邻超像素以加强连通性。总体来说，算法相比于传统 算法不同的是，算法结合 纹理特征改进了聚类相似性度量。基于 和灰度特征的区域合并 超像素分割算法得到的分割结果并不是完整水体提取结果，此分割结果属于过分割，因此需要通过区域合并将相似度高的超像素迭代合并为一个超像素来解决过分割问题。在基于层次聚类的区域合并算法中，相似性度量一般为颜色距离。但应用于遥感影像时，水体颜色与相邻地物颜色过于相似，仅仅使用颜色距离会导致误合并，因此选择将 纹理特征纳入相似性度量以优化合并性能。选用余弦相似度计算超像素与相邻超像素之间的灰度特征和纹理特征距离，即把图像表示为矢量，矢量间的余弦距离即为图像间的相似度。首先将两个超像素缩放至大小一致，获得超像素灰度或纹理特征的直方图分布，再将直方图划分为 个区，每个区包括连续 个灰度等级，对每个区的 个值进行求和得到 个数据，最终得到的 个数据为超像素灰度或纹理特征对应的矢量，这样就把两个超像素的灰度值或纹理特征转化成矢量 和矢量，见式（）。按照式（）计算 个矢量的余弦距离，余弦值（，）越接近，表明夹角越接近，个矢量越相似，即 个超像素越相似。（，），（，）（）（，）（）（，）（）（，）（）（）式中：为矢量 中的 个数据；为矢量 中的 个数据；为超像素与相邻超像素的灰度特征余弦距离；、分别为 个超像素灰度特征对应的矢量；、分别为 个超像素 纹理特征对应的矢量；为超像素与相邻超像素的 纹理特征余弦距离；为相似性度量，同时考虑灰度特征和 纹理特征距离。超像素只与相邻的余弦相似度最高且达到阈值的超像素合并，通过不断迭代，直到所有余弦相似度较高的超像素合并为一个超像素。把 图像中过分割的水体超像素合并之后就可获得完整的水体提取结果，再通过比较同一区域不同时间的水体变化来监测该区域是否发生洪涝灾害。算法评价指标采用边界召回率、欠分割错误率验证 算法的有效性和实用性。边界召回率越大，表示算法分割的边界越准确。边界召回率 计算公式为 （）式中：为算法分割边界像素出现在 边界周围两个像素区域内的像素总数，为 边界像素总数。欠分割错误率 是算法分割得到的超像素相比 分割结果的“溢出”。值越小，表示算人 民 黄 河 年第 期法分割结果相比 分割结果的“溢出”越少，分割结果越准确。计算公式为 （）式中：为算法分割结果的超像素；为 分割结果的超像素；为取与 交集最大的。结果与分析 种算法对遥感影像分割结果的比较为了验证 算法对 遥感影像的分割性能，采用 算法、融合 纹理特征的 算法、算法分别在巩义区域、渭南区域和兰州区域进行超像素分割处理。这 个研究区域的最佳分割超像素块数分别为 、，为了分割出呈窄线性的小水体，紧密度 取值都较大，分别为、。采用不同算法在最佳分割超像素块数和最佳紧密度下得到的超像素分割图像见图 图（分割图像上方图为对应颜色方框的局部放大图，绿框、蓝框和黑框从左至右均依次对应 图像、图像、图像的放大图），种算法的输入参数统一。在 个研究区域内，算法性能差异主要体现在小水体上，从绿框和蓝框局部放大图可看出，算法和算法的水体边界分割性能都明显比 算法的好，算法存在较多的欠分割情况。从黑框局部放大图可看出，算法的分割性能比 算法的稍好一些，算法也存在欠分割情况。总的来说，算法的分割性能最好，产生的水体超像素很规整。图 采用不同算法得到的巩义区域超像素分割图像图 采用不同算法得到的渭南区域超像素分割图像图 采用不同算法得到的兰州区域超像素分割图像人 民 黄 河 年第 期 、三种算法的边界召回率、欠分割错误率见图 图，算法与 算法的分割性能较接近，均优于 算法。算法的边界召回率比 算法提高了 ，算法的欠分割错误率比 算法降低了，超像素块数越少，算法的两项指标值与其他两种算法的差异越明显。图 巩义区域不同算法的边界召回率和欠分割错误率图 渭南区域不同算法的边界召回率和欠分割错误率图 兰州区域不同算法的边界召回率和欠分割错误率 水体变化监测采用 算法分割的各区域水体超像素区域合并后的水体提取结果见图，水体轮廓图箭头左边图像为对应颜色方框的局部放大图。可以看出，提取水体轮廓较完整，因此可用 算法提取同一区域不同时间的水体来比较水体变化。图 区域合并后的水体轮廓提取的不同时间的水体轮廓见图，水体轮廓图箭头左边图像为对应颜色方框的局部放大图。图（）中红色轮廓表示提取巩义区域 年 月 日的水体轮廓，绿色轮廓表示提取巩义区域 年 月 日的水体轮廓；图（）中红色轮廓表示提取渭南区域 年 月 日的水体轮廓，绿色轮廓表示提人 民 黄 河 年第 期取渭南区域 年 月 日的水体轮廓；图（）中红色轮廓表示提取兰州区域 年 月 日的水体轮廓，绿色轮廓表示提取兰州区域 年 月 日的水体轮廓。可以看出，个研究区域的水体随时间变化都发生不同程度的变化，其中巩义区域和渭南区域的黄河水体变化较明显，对于这些水体发生明显变化的地区可做标记，在下次遇到连续暴雨天气时进行洪涝灾害防控，以降低对附近居民的伤害。图 不同时间水体轮廓对比 结论本文将 纹理特征纳入 算法的相似性度量，提出一种新超像素分割方法（），对 遥感影像进行超像素分割。经与 算法和 算法对比，算法在边缘召回率和欠分割错误率上都得到了优化，能够很好地将水体和背景分割开，并解决了小水体难以分割的问题。此外，对水体超像素进行区域合并，解决了超像素过分割问题，得到完整水体轮廓。通过提取同一区域不同时间的水体轮廓，可有效观察水体变化，提前对水体变化明显的区域做标记，为洪涝防控提供参考。参考文献：，（）：王冬梅，陈琳，冯峰面向对象的 影像水体信息提取研究人民黄河，（）：，种丹，张世强，李浩杰，等基于光学遥感的湖泊水体信息提取研究进展人民黄河，（）：，赵亚杰，王立辉，孔祥兵，等基于 和 数据融合的土地利用分类研究福建农林大学学报（自然科学版），（）：白翠，向洋，邱春霞，等基于多源遥感数据的水体提取方法研究人民黄河，（）：孙根云，邵宝婕，丁孙金衍，等基于 平台的黄河流域水体指数研究人民黄河，（）：宋熙煜，周利莉，李中国，等图像分割中的超像素方法研究综述中国图象图形学报，（）：，（）：湛南渝，李小涛，路京选，等基于 的哨兵 号雷达数据水体信息提取人民长江，（）：，：，：，（）：，（）：，：，：，（）：，（）：宋晓宇，王纪华，刘良云，等基于高光谱遥感影像的大气纠正：用 数据评价大气纠正模块 遥感技术与应用，（）：赵宇晴，陈广胜，景维鹏对简单线性迭代聚类算法改进的遥感影像超像素分割方法东北林业大学学报，（）：，：，（）：【责任编辑 栗 铭】人 民 黄 河 年第 期