一种基于城市水体指数与分形几何算法的OLI遥感影像水体提取方法

城市河流是城市泄洪排涝、饮用水源或生态景观的重要组成部分。随着城市的扩张与人口的陡增，城市河流的管理与监测面临着严峻的压力[1]。遥感具有面积大、周期短、实时性、信息量大等优势，为水体信息提取提供了新的监测技术。高分辨率卫星影像被广泛用于水体信息的提取[2-4]，并已成为水资源调查、洪水灾害预测评估、水利规划、环境监测等工作的重要技术手段之一[5-6]。

目前，利用遥感手段提取水体信息的主要方法有：单波段阈值法、水体指数法[7-12]、面向对象法[13-16]等。阈值法的特点是简单易行，但提取精度不高，尤其是对城市中常见的细小且严重污染的水体提取效果不稳定[17]；水体指数法是通过特征波段间的比值运算，以区分水体与其他地物，但运用于城市水体提取时，易受到道路、建筑物等影响而出现斑点[18]；面向对象法将图像分割为由若干个同质像元组成的对象，将地物的光谱信息与空间结构信息、纹理信息结合，适合复杂环境中的地物提取[19-20]，是目前遥感水体提取的研究热点之一。

针对严重污染的城市水体与道路、建筑物、阴影等易于混分，以及遥感水体提取结果不连续、存在斑点问题，本研究以广州北部、东部的两条主要河流的OLI遥感影像为研究对象，基于城市水体指数，结合八邻域像元分析对水体进行标注，通过设置形状、面积等特征阈值，实现城市复杂环境下的城市河流信息的自动提取。

1 数据介绍及预处理

1.1 数据简介

Landsat 8携带的两个主要载荷之一的陆地成像仪(operational land imager，OLI)包含覆盖了可见光、近红外、热红外的9个通道，多光谱通道空间分辨率为30 m，辐射分辨率为16 bit，成像宽幅为185 km×185 km，重访周期为16 d。较高的空间分辨率、辐射分辨率、幅宽、重访周期以及水色通道能够满足城市河流监测与信息提取需求，是目前水体信息监测的重要数据资源。

广州市是我国一线城市，水系发达。流溪河水系广州段属于山区地形，部分河段及支流受工业污水、生活废水影响水体污染严重[21]；东江水系广州段位于广州与东莞交界处，水体受到工业的严重污染，周边建筑物、道路网密集[22]。因此，本文选择广州市北部流溪河水系和东部东江水系作为试验区。收集了流溪河(2016年12月7日，2017年4月30日)与东江(2016年2月7日，2017年10月23日)共4景少云的OLI影像用于本文的研究(详细信息见表1)。

在行业内就业的37位毕业生中，在煤炭行业内就业的为17人，占毕业生总数的45.95%；在金属和非金属矿产行业内就业的为12人，占32.43%；在其他地质行业内就业的为8人，占21.62%(图2)。

1.2 数据预处理

[13] 陈晨，陈静欣，李向前,等.面向对象的高分辨率影像水体信息提取研究[J].人民黄河，2013，35(9):68-70.

支墩为钢管结构，在顶部设置落梁，为保证整体稳定性，支墩间采用专门的连接件进行连接。承重梁为贝雷梁，共27排，贝雷梁间通过连接件实现连接，腹板以下的间距按照30cm控制，其他则按60cm控制。在贝雷梁的上部设置方木，其间距为25cm。同时在方木的基础上设置竹胶板，将其作为底部模板。梁体外侧方木避免设置钢轨，为龙门吊的行走提供轨道。由木模制成外模，梁底部与顶部分别设置拉紧予以加固。通过拌和站对混凝土进行集中搅拌，采用罐车将搅拌好的混凝土运输到施工现场，再由泵车进行连续浇筑[1]。

表1 试验中采用的OLI卫星影像

序号研究区日期云量/(%)轨道号行号列号12流溪河2016-12-0712.152017-04-3018.931224434东江2016-02-07 2.44 2017-10-230.0512144

2 研究方法

本文研究思路如图1所示，首先针对城市污染水体提出采取城市水体指数(CWI)实现城市水体的提取。其次对影像进行预处理，采用改进的归一化差异水体指数(MNDWI)判断影像中是否存在云：若存在云，则采用MNDWI先去除云阴影，然后采用CWI指数提取城市水体信息；若不存在云，则直接采用CWI指数提取城市水体信息，然后采用短波红外通道消除部分建筑物等噪声。最后采用分形几何算法消除其他噪声影响，基于形态学闭运算与开运算提取连续的河流信息。

2.1 城市水体指数

由于水体在短波红外通道强吸收使得该通道反射率接近于0，植被与建筑物在短波红外远大于0，因此许多学者采用短波红外单通道提取水体信息[4-6]。随着水体中悬浮物的增加，水体在短波红外的反射率远大于0，学者们基于绿光通道与短波红外通道提出了MNDWI指数法[23]，该算法对于混浊的水体达到良好效果。然而城市水体受到生活污水、工业废水等严重污染，以及复杂外界条件的影响，严重污染的城市水体在绿光通道反射率降低，导致污染水体MNDWI值降低，与道路易于混分。相关研究表明污染水体对绿光通道比较敏感[24]，因此笔者基于前人的研究基础提出了CWI指数为

(1)

式中， ρgreen与ρswir分别对应OLI影像的第3波段与第6波段。试验表明，阈值设为1.2可较好地识别城市污染水体与正常水体，且能避免部分建筑物、道路与桥梁等影响。

图1 总体流程

2.2 基于分形几何算法等的噪声消除

城市建筑物复杂多样，部分建筑物的CWI值与污染水体的CWI值相近，通过试验发现城市污染水体在短波红外通道的反射率值低于0.15，而混分的部分建筑物的反射率值高于0.15，采用短波红外单通道阈值法可剔除易于混分的建筑物噪声。

云的阴影与污染水体光谱近似，CWI指数难以区分云阴影与污染水体。笔者通过试验发现，云阴影的MNDWI值小于0.2，而水体的MNDWI值大于0.2，采用MNDWI指数能够识别云阴影。因此用CWI指数提取水体信息前，判断影像中是否存在云，若存在云，则采用MNDWI先去除云阴影，然后采用CWI指数提取城市水体信息；若不存在云，则直接采用CWI指数提取城市水体信息。

[5] 高佳华，王文惠，韩鑫,等.基于遥感影像的喀什地区水体提取[J].西部资源，2017(5)：131-133.

朱涛等[32]通过室内试验分析了宽级配砾质土的压实特性，结果表明，宽级配砾质土击实时干密度和含水率关系曲线呈上凸的抛物线型，具有一个最优含水率。击实后的最大干密度随着砾石含量的增加先迅速增大，后又开始逐渐减小。随着砾石含量的增加，其最优含水率先逐渐减小，后基本趋于稳定状态。

图2 连接分量示意图

2.3 数学形态学闭运算

利用水体指数可以快速准确实现大面积河流的提取，然而，对于易受到建筑物遮挡的细小支流提取时会出现河流间断。笔者通过对比研究了形态学的膨胀与腐蚀等相关计算，发现基于数学形态学的闭运算可解决城市细小支流出现的部分间断问题，同时又较好地保持了河流形状信息。

A·B=(A⊕B)ΘB

(2)

目前，纸病在线检测系统普遍采用“电荷耦合器件(Charge Coupled Device，CCD)相机+计算机”的硬件结构模式，通过CCD高速相机实时采集纸张图像数据并传送至检测系统PC端，系统的核心处理器利用图像处理技术对纸张图像进行处理和辨识，实现对黑斑、孔洞、划痕、裂边等缺陷的识别和分类[1- 3]。但是随着纸机车速的提升、纸幅的加宽以及对纸病辨识精度的提高，CCD相机采集到的海量图像数据，给系统的快速性甚至稳定性带来了巨大冲击，如何解决打破纸病辨识快速性这一瓶颈已成为当前研究纸病检测技术的关键问题。

3 结果与分析

3.1 CWI指数敏感性分析

CWI采用OLI影像的第3通道与第6通道计算，第3通道为绿光通道受大气影响严重，第6通道为短波红外通道受大气影响较弱，由于大气气溶胶等复杂多变，不同日期、不同区域的影像大气纠正前MNDWI差异显著，本文对影像进行大气纠正，大气纠正后CWI均小于或等于1.2。试验表明阈值设为1.2可较好地识别水体，尤其是城市污染水体，可以避免绝大部分建筑物、植被等影响。

PALOMA-1 [4]是一项Ⅱ期、开放的临床研究，入组的165例HR＋和HER2-绝经后不适合手术的转移或局部复发乳腺癌患者按1∶1的比例随机分为来曲唑单药（2.5 mg/d，持续给药）和来曲唑联合palbociclib（125 mg/d，连续3周停1周，4周为1个周期）2组。结果显示，联合组对比单药组的中位无进展生存期（median progressionfree survival，mPFS）为20.2个月和10.2个月[风险比（hazard ratio，HR）＝0.488，95%可信区间（confidence interval，CI）：0.319～0.748，P＜0.001]。

根据设计要求，履带式多吸盘行走机构对工作平面的适应性强、行进速度较快，可越过一般的障碍物。这种运动方案不仅可以使光伏清洁机器人安全可靠地在太阳能电池板表面上行走，并能顺利地克服排布间隙对吸附系统的影响。所以采用履带行走方式配合多个真空吸盘吸附的运动方案。

3.2 水体提取结果及分析

六是产业市场化。电气自动化企业在集中精力研发新技术、攻关新课题的同时，还将以市场化的形式采取外包的方式带动配套产业的市场化发展，有计划地对重大技术装备进行开发与研究，实现自动化产业的市场化。

表2 大气纠正前后CWI值对比

研究区影像日期大气纠正前CWI大气纠正后CWI水体植被建筑物水体植被建筑物流溪河2016-12-07>1.17<0.67<1.17>=1.20<1.20<1.202017-04-30>1.01<0.70<1.01>=1.20<1.20<1.20东江2016-02-07>1.33<0.86<1.33>=1.20<1.20<1.202017-10-23>1.29<0.83<1.29>=1.20<1.20<1.20

图3(a)与图4(a)为大气纠正后的流溪河流域影像，流溪河提取结果如图3(b)—(f)与图4(b)—(f)所示。由于山体阴影、坑塘等的影响，MNDWI(如图3(c)与图4(c)所示)算法与SVM算法提取的河流较连续，但是存在大量斑点；单通道阈值法与光谱角度算法提取的河流存在少量斑点，但是提取的河流出现间断；本文算法提取的河流连续且无噪声，图4(f)河流出现部分间断原因是受云的影响。结果表明，本文提出的算法能够实现河道宽度大于30 m的严重污染水体提取，并能克服山体阴影、坑塘以及桥梁等影响，适用于山区城市河流的提取。

本文通过监测CSP是否履行了SLA来作为评估CSP信任的基础。并且在监测过程中利用马尔科夫模型的预测功能，在一定程度上允许CSP提前调整SLA规定的性能指标，最大限度地保证了CSP的履约性。

使用该模型开发一款水光面膜配方，筛选19名消费者（排除条件同1.1），年龄为22～41周岁，其中干性肌肤和中性肌肤9名，混合性和油性肌肤10名。采用半脸随机单盲测试（盲测针对专家视觉评估人员和仪器测试人员），按1.6模型计算产品使用前后拟合水光指数提升值、拟合水光指数提升百分比，并将水光面膜产品使用组和空白组进行对照。

图5(a)与图6(a)为大气纠正后的东江流域影像，东江水体及其支流体提取结果如图5(b)—(f)与图6(b)—(f)所示。由于城市建筑物、道路、坑塘等影响，MNDWI算法与SVM算法能够较好地提取东江及其细小支流，但是提取结果存在大量斑点；单通道阈值法与光谱角度算法对大面积水体提取效果较好，但是细小的支流出现间断，另外也存在少量斑点。本文算法提取的东江及其细小支流连续且无斑点噪声。结果表明，本文提出的算法能够克服建筑物、道路坑塘等影响，可适用于城市建筑密集条件下的河流提取。

图3 2016年12月7日流溪河水体提取结果

图4 2017年4月30日流溪河水体提取结果

4 结 论

本文以广州市流溪河与东江水系为研究对象，基于2016年与2017年OLI遥感影像，采用CWI水体指数法，同时结合分形几何算法，通过设置形状面积等特征，实现城市复杂环境下水体信息的自动提取。并与单通道算法、MNDWI水体指数算法、SVM算法与光谱角度算法的水体提取结果进行了对比分析。具体结论如下：

将单通道阈值、MNDWI算法、SVM(支持向量机)算法、光谱角度算法以及本文改进的算法分别用于2016年12月7日、2017年4月30日大气纠正后的流溪河流域OLI影像以及2016年2月7日、2017年10月23日大气纠正后的东江流域OLI影像。

图5 2016年2月7日东江水体提取结果

图6 2017年10月23日东江水体提取结果

(1) CWI受大气影响显著，本文对影像进行大气纠正，大气纠正后CWI均大于或等于1.2，试验表明阈值设为1.2可较好地识别城市水体。

当1≥e1>e2>0时，对任意的w1，定理3结合定理1，有q*(e1,w1)>q*(e2,w1)，如图5所示。

(2) 将单通道阈值、MNDWI算法、SVM算法、光谱角度算法以及本文改进的算法分别用于2016年12月7日、2017年4月30日大气纠正后的流溪河流域OLI影像以及2016年2月7日、2017年10月23日大气纠正后的东江流域OLI影像。结果表明本文提出的算法能够较好地提取出大面积的河流及细小支流，提取结果连续且无斑点噪声，能够克服建筑物、道路坑塘等的影响，完全适用于城市建筑密集条件下的污染水体与一般水体提取。

本文研究采用30 m分辨率OLI影像进行河流提取，由于受到空间分辨率限制，难以提取细小河流，下一步计划将本文算法用于高空间分辨率影像。本文提出的算法能够实现城市污染水体以及一般水体连续、准确提取，为城市水资源的管理与保护提供基础数据支撑。

参考文献：

城市水体污染较严重，部分河流出现黑臭而且由于污染物浓度过高，导致水体在绿光通道有较高的反射，CWI指数中水体与部分建筑物难以区分。但是由于水体在红短波红外通道强吸收，在短波红外通道水体反射率小于部分建筑物反射率，而部分建筑物反射率小于植被反射率，可以设置一定的阈值将城市严重污染的水体与道路分开提取。

[1] 杜云艳，周成虎.水体的遥感信息自动提取方法[J].遥感学报，1998，2(4):364-369.

[2] 周蕾，崔云霞，韦玉春,等.空间分辨率对水体遥感信息提取的影响[J].南京师范大学学报(工程技术版)，2015，15(4):66-74.

[3] 程滔，刘若梅，周旭.基于高分辨率遥感影像的地理国情普查水体信息提取方法[J].测绘通报，2014(4):86-89.

对照组:应用生物敷料进行创面修复,敷料包含微量元素、类黏蛋白、胶原纤维等,首先对创面进行清理和止血,顺延创面边缘将无生物活性或是焦痂的坏死组织进行清除,经电凝进行止血,应用经稀释后的碘伏进行清洗,用无菌纱布包扎,每隔2天更换一次纱布。

[4] 杨果果，高原，周蕾.水体信息遥感自动提取与演变分析——以南京城区为例[J].测绘与空间地理信息，2016(10):35-39.

在初步提取的结果中，1代表水体(前景)，0代表其他地物(背景)，由于城市水体还受到其他噪声等影响，存在大量孤立分量，是需要消除的。本文将各连接分量采用空间邻域法进行标注(如图2所示)，将前景影像分割成不同的子集，并用数值进行标注，然后采用直方图统计面积，通过交互式方法设置阈值面积，有针对性提取水体。

[6] 陈婧，邓一兵.基于ETM+遥感影像水体信息提取方法探讨[J].测绘与空间地理信息，2011，34(1):177-180.

传统工艺的七个本质特征与其价值是一一对应的:实用的品格赋予它使用价值和交换价值,也就是经济和民生价值;审美的品格使传统工艺具有很高的艺术价值和美育价值,传统工艺不仅给人们以美的享受,还起着陶冶情操、培育国民人文素养等重要作用;理性的品格赋予传统工艺丰富多样的科学技术基因和重要的学术价值。只要承认传统工艺具有这些品格,传统工艺就有现代价值,这是符合内在逻辑的。传统工艺的品格具有永恒性,因此它是与人类共存的，它的价值不会被湮没。而且,人的本性是愿意亲近手工艺的。

[7] 骆剑承，盛永伟，沈占锋,等.分步迭代的多光谱遥感水体信息高精度自动提取[J].遥感学报，2009，13(4):604-615.

[8] 张伟，赵理君，郑柯,等.一种改进光谱角匹配的水体信息提取方法[J].测绘通报，2017(10):34-38.

[9] 钱立辉，于成龙，徐进,等.基于Canny算子与决策树模型的水体信息提取研究[J].测绘与空间地理信息，2016(10):178-181.

[10] 何海清，杜敬，陈婷,等.结合水体指数与卷积神经网络的遥感水体提取[J].遥感信息，2017(5):82-86.

式中，A为待处理图像；B为计算模板。闭运算先用模板B对待处理图像进行膨胀，然后再对处理后的图像进行腐蚀。闭运算可有效实现间断线的连接。

[11] 李小涛，黄诗峰，郭怀轩.基于纹理特征的SPOT 5影像水体提取方法研究[J].人民黄河，2010，32(12):5-6.

[12] 童李霞，燕琴，骆成凤,等.基于NDWI分割与面向对象的水体信息提取[J].地理空间信息，2017，15(5):57-59.

对OLI影像的预处理主要包括精确大气校正和边界切割处理。卫星影像可见光通道受大气影响严重，传感器接收到的地表反射信息，两次经过大气，大气中气体分子、气溶胶、云雾水滴、冰晶等对太阳光线的作用使得传感器接收到的地表信息发生变化；且水体反射信息微弱，大气信息在总信号中占到90%以上。因此，遥感影像上精确提取水体前需对遥感影像进行精确的大气纠正。当前大气纠正可以采用的方法有基于统计分析的大气校正方法和基于辐射传输模型的大气纠正方法。基于统计分析的大气校正方法需要有同步的地面实测数据做支持，该方法不具有较好的普适性；基于辐射传输模型的大气纠正方法具有较好的物理意义和普适性，因此研究中采用了基于辐射传输模型的Flaash模型进行大气纠正。由于缺乏同步的大气信息，在大气纠正中，采用迭代求解大气气溶胶光学厚度。

[14] 李士进，王声特.基于混合特征空间MRF(Markov Random Filed)模型的高分辨率遥感影像水体提取[J].南京师大学报(自然科学版)，2017，40(1):13-19.

[15] 殷亚秋，李家国，余涛,等.基于高分辨率遥感影像的面向对象水体提取方法研究[J].测绘通报，2015(1):81-85.

[16] 杜斌.基于面向对象的高分辨率遥感影像水体信息提取优势研究[D].昆明：云南师范大学，2014.

[17] 李艳华，丁建丽，闫人华.基于国产GF-1遥感影像的山区细小水体提取方法研究[J].资源科学，2015，37(2):408-416.

这句话不仅攻击征兵宣传片嗜杀黩武，还暗示中国的青年人对此十分买账，以至于宣传部门如果想要得到青年人注意，就得用“straight outta Beiijng”这种露骨的宣战字眼。这无疑同西方叙事中中国民族主义泛滥的观点相关联，符合西方观众对新闻事件主体（中国和解放军）的预期，从而造就了新闻的负面性和一致性价值。

[18] 张红梅，王大卫，高杨,等.基于OLI数据与决策树法的去山体阴影水体信息提取研究[J].测绘工程，2017，26(11):45-48.

[19] 李长春，张光胜，慎利.基于AdaBoost算法的遥感影像水体信息提取[J].测绘科学，2013，38(2)：104-106.

[20] 宋英强，杨联安，许婧婷,等.基于Landsat-8卫星OLI遥感影像和AdaBoost算法的水体信息提取[J].测绘地理信息，2017，42(3):44-47.

[21] 高强.珠三角感潮河网水体黑臭评价方法初探[J].生态与环境工程，2011(16)：200-202.

[22] 杨宇，董天明，王丽萍,等.广州市重点整治河涌水质指数WQI的研究思路与应用[J].广东化工，2017，44(13):171-172.

[23] 徐涵秋.利用改进的归一化差异水体指数(MNDWI)提取水体信息的研究[J].遥感学报，2005，9(5):589-595.

[24] 童小华，谢欢，仇雁翎,等.黄浦江上游水域的多光谱遥感水质监测与反演模型[J].武汉大学学报(信息科学版),2006，31(10):851-854.