一种基于高分辨率卫星影像的入海排口疑似排污识别方法

技术领域

本发明属于海洋环境保护技术领域，尤其是涉及一种基于高分辨率卫星影像的入海排口疑似排污识别方法。

背景技术

入海排口是指通过溪闸、管道、沟渠等直接向近岸海域排放内陆水体的口门，呈现出面积小、分布广的特点。特殊的地理位置使其受到陆源污染、河流输入、泥沙再悬浮、邻接效应等多重影响，水体性质及污染来源复杂，是生态环境治理的重要节点。由于陆源污染物通过入海排口直排入海会影响近岸海域生态环境质量，体现在氨氮、总磷、总氮、化学需氧量、生化需氧量及其他水质指标超地表水环境质量标准或海水水质标准，因此对入海排口排水状况的监测与管理十分重要。

近年来，各级海洋环境监测部门进行了详细的入海排口排查工作，但是传统的人工现场巡查费时费力，并且难以快速且全面地摸清入海排口处排水状况，导致排查成果不完善且效率低。目前，部分重点入海排污口已经安装水质自动监测装置，但由于入海排口数量巨大，大部分排口目前缺乏在线监测设备，是当前海洋环境监管的“痛点”。卫星遥感具有高空间分辨率、高时间分辨率观测的特点，可以实现对沿岸入海排口的快速响应与动态监测，有利于全面摸清入海排口数量、分布、动态变化等特征。

现有技术中已经出现了利用卫星遥感进行水体监测的技术，例如Zhao, Y等人在Recognition of Water Colour Anomaly by Using Hue Angle and Sentinel 2 Image（Remote Sensing 12, 716）一文中利用色度角的方法来判断是否为水色异常水体，当色度角超过阈值则认为是水色异常水体。但是，该方法只适用于光学深水，对光学浅水不适用，且不适用于富营养化水体的识别。目前，对于入海排口的污染监管，有待提供一种适用场合更广、通用性更强、识别更加准确的卫星监测方法。

发明内容

针对现有技术中存在的问题，本发明提供了一种基于高分辨率卫星影像的入海排口疑似排污识别方法，该方法主要包括以下步骤：

S1：通过高分辨率卫星获得多种水色类型、多时相的遥感反射率光谱，并对遥感反射率进行归一化处理，突出谱形特征，构成由归一化的遥感反射率光谱形成的训练样本集；其中，获取的多种水色类型除了异常水色也包含其他类型的水色，如高悬浮泥沙水体、离岸较远的清洁水体、水库水体、河道水体、养殖池塘水体、邻近岸的水体等，从而建立覆盖各个区域、季节、水色类型的遥感反射率光谱库。

S2：依据光谱之间的相似性特征，对训练样本集中的归一化的遥感反射率光谱进行聚类处理，得到多种水体光学类型及各自对应的中心光谱；根据所述中心光谱的光谱特征，确定各类水体光学类型所对应的水体状态；所述水体状态包括清洁水体状态、浑浊水体状态、富营养化水体状态、泥沙含量极高的水体状态或是容易被误分的滩涂，以及疑似污染水体状态。

S3：根据各类水体光学类型所对应的水体状态及光谱特征，提取确定对不同水体光学类型或水体状态具有敏感度的多个光谱参数，并将所述训练样本集中归一化的遥感反射率光谱的光谱参数、对应的水体光学类型作为真值，对机器学习模型进行迭代训练，得到水体光学类型分类模型；

S4：通过高分辨率卫星获取包括入海排口的目标区域的遥感反射率光谱，将提取确定的光谱参数输入所述水体光学类型分类模型，并以输出的水体光学类型作为识别结果。

具体的，所述步骤S1在进行归一化处理前，对高分辨率卫星获得的遥感影像进行了陆地掩膜、云掩膜处理，仅保留晴空水体像元；然后对晴空水体像元进行大气校正，提取得到晴空水体像元的遥感反射率光谱。优选的，基于通用瑞利散射查找表，采用适用于近岸区域的短波红外大气校正方法进行大气校正。

具体的，所述步骤S2中采用非监督聚类，优选使用Kmeans++聚类。通过聚类处理，至少能够获得五种以上的水体光学类型及各自对应的中心光谱，通过分析中心光谱的光谱特征，例如特征峰、特征谷、变化趋势等，可以确定各类水体光学类型所对应的水体状态，有利于疑似排污的遥感识别。

具体的，所述步骤S3中采用的机器学习模型优选为随机森林分类模型。

具体的，所述步骤S4中的所述目标区域包括入海排口在内的大气校正后的遥感影像，所述步骤S4逐个提取确定遥感影像上各个像元的光谱参数并输入所述水体光学类型分类模型，得到各个水体像元的水体光学类别，根据水体光学类型对应的水体状态，对是否存在疑似排污进行判断。

进一步的，所述方法还包括步骤S5：根据多期历史遥感影像获得所述目标区域的历史识别结果，追溯其是否长期排污，分析异常排水的过程及演变。

进一步的，所述方法还包括步骤S6：基于新的入海排口排污事件和对应的遥感影像，提取遥感影像上污染像元的遥感反射率光谱，对所述训练样本集进行动态更新，进一步完善水体光学类型分类模型，提高入海排口疑似排污识别准确率。

本发明提供的方法基于高分辨率卫星影像的光谱特征，建立入海排口疑似排污识别模型，可以快速筛查疑似污染排口，对沿岸入海排口水色状况进行动态全面监测，迅速、详细掌握污染产生、演变过程，有利于从面上对入海排口进行监督管理，节省人力、物力等资源，严格监督陆源污染物入海，改善近岸海域生态环境质量。

附图说明

通过结合附图对本公开示例性实施例进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施例中，相同的参考标号通常代表相同部件。

图1示出了具体实施方式中识别方法的流程图；

图2示出了具体实施方式中归一化的遥感反射率光谱的聚类结果；

图3示出了具体实施方式中的入海排口疑似排污识别示意图。

具体实施方式

下面将参照附图更详细地描述本公开的实施例。虽然附图中显示了本公开的实施例，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些例举性的实施例是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。

在本文中使用的术语“包括”及其变形表示开放性包括，即“包括但不限于”。除非特别申明，术语“或”表示“和/或”。术语“基于”表示“至少部分地基于”。术语“连接”、“连通”表示直接或通过其他组件间接的连接或连通。术语“第一”、“第二”等可以指代不同的或相同的对象，但并不直接表示先后顺序或重要程度的不同。下文还可能包括其他明确的和隐含的定义。

本具体实施方式以Sentinel-2的遥感影像为例，提供了一种基于高分辨率卫星影像的入海排口疑似排污识别方法，如图1所示，具体包括以下步骤：

（1）获取多时相的Sentinel-2高分辨率多光谱遥感影像，空间分辨率为10米，双星组合，重访周期为5天，各波段波长及波段号如下表所示：

表1. Sentinel-2各波段波长及对应波段号简称。

（2）对获取的多景遥感影像进行陆地掩膜，获取水体区域。基于近红外波段和红光波段两个波段的大气顶反射率，具体为865nm的大气顶反射率R(865nm)以及665nm的大气顶反射率R(665nm)，计算归一化植被指数（NDVI），采用阈值法对陆地进行掩膜，NDVI通过下式计算得出。

（3）由于光学卫星影像会受云的干扰，需要剔除云覆盖的像元，获取晴空的水体像元。采用2190nm波段瑞利校正反射率，根据设定的阈值对云进行掩膜，最终仅保留晴空水体像元；

（4）对晴空水体像元进行大气校正，消除大气分子散射、气溶胶散射和水面反射的影响，获取水体遥感反射率。由于入海排口主要沿岸线分布，水体通常较为浑浊，因此采用适用于浑浊水体的短波红外大气校正方法。具体来讲，遥感器接收到的大气顶反射率如下式：

式中，ρ

式中，ρ

式中，τ

（5）选取大量异常水色及正常水色的像元，作为训练样本点，提取各波段的遥感反射率光谱R

式中，

（1）由于缺乏对逐像元水体光学类别的先验知识，对建立的训练样本集进行非监督聚类，使用Kmeans++聚类，根据Davies-Bouldin指数及轮廓系数(silhouettecoefficient)确定最优聚类数目，从而确定每一条光谱对应的分类类别。在本实施例中，共确定了14类水体光学类型，如附图1所示；

（2）根据对光谱形状的先验知识，清洁水体的峰值通常出现在490nm或560nm，富营养化水体通常在705nm附近出现峰值，悬浮泥沙含量较高的水体会在665nm出现峰值，并且在近红外波段的遥感反射率会抬升，而污染水体会呈现黑色、红棕色等水色，对应于遥感反射率的谱形上体现为整体反射率偏低，但仍在红光/近红外波段有明显的峰/谷值。因此，根据聚类得到的中心光谱的光谱特征，对特征峰、特征谷及变化趋势等进行分析，确定每一类水体光学类型对应的水体状态；在本实施例中，如图2所示，水体光学类型1~4对应的水体状态均是清洁水体状态，清洁度逐渐下降，5~10对应的是浑浊水体状态，且泥沙含量逐渐增加，11~12对应的是富营养化水体状态，13对应的是泥沙含量极高的水体状态或是容易被误分的滩涂，14为疑似污染水体状态。

（1）除了样本点的前9个波段的归一化遥感反射率，基于不同水体光学类型的光谱特征，即对每一类的遥感反射率光谱，判断峰值和谷值出现的波段，特别是与其他类别相比，该类特有的峰、谷值位置或者变化趋势，确定了6个额外的光谱指数(OCI1~OCI6)作为输入变量，旨在强调突出不同水体状态的光谱特征。新增的6个额外的光谱指数分别经过了多景影像阈值法的测试，呈现出对黑臭水体、富营养化水体、高悬浮泥沙水体敏感性强的特点，可以用于疑似排污和正常水体的区分。分别为：

式(7)~式(12)中，b1~b8a分别为Sentinel-2前9个波段的波段号对应的遥感反射率的值的简称,具体波长见表1。OCI1~OCI6实际代表的含义分别为归一化黑臭水体差值指数、490nm波峰基线高、665nm波峰基线高、红边指数、近红外指数和四个波段波峰基线高的平均值。

（2）基于步骤S1建立的训练样本集、上述15个光谱参数，以及作为输出的步骤S2中聚类分析确定的每一条光谱的水体光学类型，对随机森林分类模型进行训练。模型参数包含树的数目、深度等，通过设置这些参数的变化范围，进行组合与遍历，最后利用混淆矩阵，即预测的类和实际的类之间的比较，来评估分类模型的准确性，以此来确定最优的模型参数。在最优的模型参数的基础上，最后确定了基于15个光谱参数的水体光学类型分类模型。

（1）获取准实时更新的高分辨率卫星影像，对包括入海排口的目标区域进行准实时动态监测。通过与步骤S1类似的方法进行陆地掩膜、云掩膜及大气校正预处理，仅保留晴空水体像元，提取得到目标区域晴空水体像元的遥感反射率光谱；

（2）逐像元计算15个光谱参数，输入到步骤S3建立的水体光学类型分类模型，得到逐像元的水体光学类型分布图，识别出水色异常或疑似排污的入海排口。如图3所示，本实施例识别出在某入海排口附近出现明显的富营养化水体状态，闸后水塘内也存在部分疑似污染的水体状态。

进一步的，识别出异常状态后，可以根据多期历史遥感影像获得所述目标区域的历史识别结果，追溯其是否长期排污，并且根据准实时更新的影像分析异常排水的过程及演变。基于新的遥感影像和新的入海排口排污事件，可以进一步提取遥感影像上污染像元的遥感反射率光谱，对训练样本集进行动态更新，进一步完善水体光学类型分类模型。

通过本发明提供的识别方法，相关生态环境监测部门根据识别结果及时进行现场核查，并对可能的排污主体或单位进行跟踪调查与处理。

以上已经描述了本公开的各实施例，本领域技术人员应当理解上述说明仅仅是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对现有技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。