一种基于多光谱影像的废水污染监测方法和装置

摘要

本发明提供了一种基于多光谱影像的废水污染监测方法和装置，涉及环境监测的技术领域，包括：获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；对废水污染指数时间序列进行时空模式分析，得到分析结果，其中，分析结果用于表征待监测水域的水质变化趋势，解决了现有的废水污染遥感监测算法的适用性较差的技术问题。

权利要求

1.一种基于多光谱影像的废水污染监测方法，其特征在于，包括：获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；

基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；

基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；

对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势；

其中，所述光学特性数据包括：叶绿素‑a，有色溶解有机物和非藻类颗粒的吸收系数，悬浮颗粒物的后向散射系数和漫衰减系数；

基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数，包括：基于所述漫衰减系数，计算出多光谱遥感影像数据的漫衰减系数的离散和值；

基于主成分分析法对目标数据进行处理，得到多光谱遥感影像数据中每个栅格的特征向量，其中，所述目标数据包括：所述叶绿素‑a、所述有色溶解有机物和非藻类颗粒的吸收系数、所述悬浮颗粒物的后向散射系数和漫所述衰减系数的离散和值；

基于所述特征向量和谱角映射器算法，计算出多光谱遥感影像数据的废水污染指数；

其中，所述废水污染指数的计算公式为 ， 为所述废水污染指数，为中心处理和进行线性组合后的像素值， 为 像素值随时间变化的平均值，σ是随时间变化的标准差。

2.根据权利要求1所述的方法，其特征在于，基于所述特征向量和谱角映射器算法，计算出多光谱遥感影像数据的废水污染指数，包括：确定出多光谱遥感影像数据中每个栅格的特征向量，确定出每个栅格对应的目标栅格；

基于每个栅格的特征向量和每个栅格对应的目标栅格的特征向量，计算出多光谱遥感影像数据的废水污染指数。

3.根据权利要求1所述的方法，其特征在于，确定出多光谱遥感影像数据中每个栅格的特征向量，确定出每个栅格对应的目标栅格，包括：计算初始栅格的特征向量与其它栅格的特征向量之间的角度差，得到角度差集合，其中，所述初始栅格为所述多光谱遥感影像数据中的任意一个栅格，所述其它栅格为所述多光谱遥感影像数据中除所述目标栅格以外的栅格；

将所述角度差集合中的最小值对应的其它栅格，确定为所述初始栅格对应的目标栅格。

4.根据权利要求1所述的方法，其特征在于，对所述废水污染指数时间序列进行时空模式分析，得到分析结果，包括：基于主成分分析算法，对所述废水污染指数时间序列进行处理，得到所述废水污染指数时间序列的第一主成分；

基于所述第一主成分和预设置信区间，对所述废水污染指数时间序列进行最小二乘拟合，得到拟合结果；

基于所述拟合结果的线性趋势，确定出所述分析结果。

5.一种基于多光谱影像的废水污染监测装置，其特征在于，包括：获取单元，用于获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；

计算单元，用于基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；

构建单元，用于基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；

分析单元，用于对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势；

其中，所述光学特性数据包括：叶绿素‑a，有色溶解有机物和非藻类颗粒的吸收系数，悬浮颗粒物的后向散射系数和漫衰减系数；

所述计算单元，用于：

基于所述漫衰减系数，计算出多光谱遥感影像数据的漫衰减系数的离散和值；

基于主成分分析法对目标数据进行处理，得到多光谱遥感影像数据中每个栅格的特征向量，其中，所述目标数据包括：所述叶绿素‑a、所述有色溶解有机物和非藻类颗粒的吸收系数、所述悬浮颗粒物的后向散射系数和漫所述衰减系数的离散和值；

基于所述特征向量和谱角映射器算法，计算出多光谱遥感影像数据的废水污染指数；

其中，所述废水污染指数的计算公式为 ， 为所述废水污染指数，为中心处理和进行线性组合后的像素值， 为 像素值随时间变化的平均值，σ是随时间变化的标准差。

6.根据权利要求5所述的装置，其特征在于，所述计算单元，用于：确定出多光谱遥感影像数据中每个栅格的特征向量，确定出每个栅格对应的目标栅格；

基于每个栅格的特征向量和每个栅格对应的目标栅格的特征向量，计算出多光谱遥感影像数据的废水污染指数。

7.一种电子设备，其特征在于，包括存储器以及处理器，所述存储器用于存储支持处理器执行权利要求1至4任一项所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。

8.一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，其特征在于，计算机程序被处理器运行时执行上述权利要求1至4任一项所述方法的步骤。

说明书

一种基于多光谱影像的废水污染监测方法和装置

技术领域

[0001] 本发明涉及环境监测的技术领域，尤其是涉及一种基于多光谱影像的废水污染监测方法和装置。

背景技术

[0002] 目前的废水污染遥感监测技术，以多光谱影像数据为主，技术路线为基于机器学习，利用光谱特征信息和地面站点观测值数据作为训练数据，构建模型进行监测识别。

[0003] 在现有技术中，利用遥感数据和站点实测数据构建机器学习模型，该方法是数据驱动的，具有机器学习方法自身的缺陷，由于训练数据集的固有限制，导致最终模型的泛化能力和鲁棒性并不好，只适用于部分场景和数据。

[0004] 针对上述问题，还未提出有效的解决方案。

发明内容

[0005] 有鉴于此，本发明的目的在于提供一种基于多光谱影像的废水污染监测方法和装置，以缓解了现有废水污染遥感监测算法的适用性较差的技术问题。

[0006] 第一方面，本发明实施例提供了一种基于多光谱影像的废水污染监测方法，包括：获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势。

[0007] 进一步地，所述光学特性数据包括：叶绿素‑a，有色溶解有机物和非藻类颗粒的吸收系数，悬浮颗粒物的后向散射系数和漫衰减系数；基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数，包括：基于所述漫衰减系数，计算出多光谱遥感影像数据的漫衰减系数的离散和值；基于主成分分析法对目标数据进行处理，得到多光谱遥感影像数据中每个栅格的特征向量，其中，所述目标数据包括：所述叶绿素‑a、所述有色溶解有机物和非藻类颗粒的吸收系数、所述悬浮颗粒物的后向散射系数和漫所述衰减系数的离散和值；基于所述特征向量和谱角映射器算法，计算出多光谱遥感影像数据的废水污染指数。

[0008] 进一步地，基于所述特征向量和谱角映射器算法，计算出多光谱遥感影像数据的废水污染指数，包括：确定出多光谱遥感影像数据中每个栅格的特征向量，确定出每个栅格对应的目标栅格；基于每个栅格的特征向量和每个栅格对应的目标栅格的特征向量，计算出多光谱遥感影像数据的废水污染指数。

[0009] 进一步地，确定出多光谱遥感影像数据中每个栅格的特征向量，确定出每个栅格对应的目标栅格，包括：计算初始栅格的特征向量与其它栅格的特征向量之间的角度差，得到角度差集合，其中，所述初始栅格为所述多光谱遥感影像数据中的任意一个栅格，所述其它栅格为所述多光谱遥感影像数据中除所述目标栅格以外的栅格；将所述角度差集合中的最小值对应的其它栅格，确定为所述初始栅格对应的目标栅格。

[0010] 进一步地，对所述废水污染指数时间序列进行时空模式分析，得到分析结果 ，包括：基于主成分分析算法，对所述废水污染指数时间序列进行处理，得到所述废水污染指数时间序列的第一主成分；基于所述第一主成分和预设置信区间，对所述废水污染指数时间序列进行最小二乘拟合，得到拟合结果；基于所述拟合结果的线性趋势，确定出所述分析结果。

[0011] 第二方面，本发明实施例还提供了一种基于多光谱影像的废水污染监测装置，包括：获取单元，用于获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；计算单元，用于基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；构建单元，用于基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；分析单元，用于对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势。

[0012] 进一步地，所述光学特性数据包括：叶绿素‑a，有色溶解有机物和非藻类颗粒的吸收系数，悬浮颗粒物的后向散射系数和漫衰减系数；所述计算单元，用于：基于所述漫衰减系数，计算出多光谱遥感影像数据的漫衰减系数的离散和值；基于主成分分析法对目标数据进行处理，得到多光谱遥感影像数据中每个栅格的特征向量，其中，所述目标数据包括：所述叶绿素‑a、所述有色溶解有机物和非藻类颗粒的吸收系数、所述悬浮颗粒物的后向散射系数和漫所述衰减系数的离散和值；基于所述特征向量和谱角映射器算法，计算出多光谱遥感影像数据的废水污染指数。

[0013] 进一步地，所述计算单元，用于：确定出多光谱遥感影像数据中每个栅格的特征向量，确定出每个栅格对应的目标栅格；基于每个栅格的特征向量和每个栅格对应的目标栅格的特征向量，计算出多光谱遥感影像数据的废水污染指数。

[0014] 第三方面，本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述第一方面中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。

[0015] 第四方面，本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序。

[0016] 在本发明实施例中，通过获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势，达到了无需机器学习模型就可以对废水污染进行监测的目的，进而解决了现有的废水污染遥感监测算法的适用性较差的技术问题，从而实现了提高废水污染遥感监测算法适用性的技术效果。

[0017] 本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。

[0018] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。

附图说明

[0019] 为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

[0020] 图1为本发明实施例提供的一种基于多光谱影像的废水污染监测方法的流程图；

[0021] 图2为本发明实施例提供的一种基于多光谱影像的废水污染监测装置的示意图；

[0022] 图3为本发明实施例提供的一种电子设备的示意图。

具体实施方式

[0023] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0024] 实施例一：

[0025] 根据本发明实施例，提供了一种基于多光谱影像的废水污染监测方法的实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。

[0026] 图1是根据本发明实施例的一种基于多光谱影像的废水污染监测方法的流程图，如图1所示，该方法包括如下步骤：

[0027] 步骤S102，获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；

[0028] 在本发明实施例中，多光谱遥感影像数据为Sentinel‑2遥感影像数据。

[0029] 步骤S104，基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；

[0030] 步骤S106，基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；

[0031] 步骤S108，对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势。

[0032] 在本发明实施例中，通过获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势，达到了无需机器学习模型就可以对废水污染进行监测的目的，进而解决了现有的废水污染遥感监测算法的适用性较差的技术问题，从而实现了提高废水污染遥感监测算法适用性的技术效果。

[0033] 在本发明实施例中，利用2SeaColor解析模型，得到叶绿素‑a、有色溶解有机物和非藻类颗粒的吸收系数、悬浮颗粒物的后向散射系数。然后结合区域内的时序多光谱数据（以Sentinel‑2为例）和区域实测的遥感反射率数据，得到漫衰减系数。

[0034] 其中，2SeaColor模型中的变量参数如下所示：

[0035] 地表遥感反射率 ，其中， 辐亮度和辐照度的换算系数， 是水的折射率， 和 分别为空气对水和水对空气的透过率，
为水对空气的反射率。

[0036] 地表下辐照度反射率 ，其中，是漫反射光的比例。

[0037] 半无限介质的双半球反射率 。

[0038] 半无限介质的定向半球面反射率 ，其中， 是水面下太阳天顶角的余弦值。

[0039] 后向散射吸收比 。

[0040] 总吸收系数 ，其中， 、 、 下标分别对应的水、叶绿素‑a、有色溶解有机物和非藻类颗粒，为波长。

[0041] 水分子吸收 。

[0042] 叶绿素‑a吸收 ，其中， 和 是相关参数。

[0043] 有色溶解有机物和非藻类颗粒吸收 ，其中，是 的光谱斜率。

[0044] 水分子反向散射 。

[0045] SPM反向散射 ，其中，是 的指数律。

[0046] 漫衰减系数 其中，是直射阳光的消光系数， 是直射阳光的前向散射，α是漫反射消光系数减去漫反射前向散射系数，是对直射阳光的后向散射做的积分。

[0047] 由上述内容可知，2SeaColor模型的反演结果由13个波段组成，分别如下：

[0048] （443）【Band 1】，（490）【Band 2】，（560）【Band 3】， （665）【Band 4】，（705）【Band 5】， （740）【Band 6】， （783）【Band 7】， （842）【Band 8】， （440）【Band 9】， （440）【Band 10】， （440）【Band 11】，【Band 12】，【Band 13】。

[0049] 在本发明实施例中，步骤S104包括如下步骤：

[0050] 基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数，包括：

[0051] 基于所述漫衰减系数，计算出多光谱遥感影像数据的漫衰减系数的离散和值；

[0052] 基于主成分分析法对目标数据进行处理，得到多光谱遥感影像数据中每个栅格的特征向量，其中，所述目标数据包括：所述叶绿素‑a、所述有色溶解有机物和非藻类颗粒的吸收系数、所述悬浮颗粒物的后向散射系数和漫所述衰减系数的离散和值；

[0053] 基于所述特征向量和谱角映射器算法，计算出多光谱遥感影像数据的废水污染指数。

[0054] 具体的，基于所述特征向量和谱角映射器算法，计算出多光谱遥感影像数据的废水污染指数，包括如下步骤：

[0055] 确定出多光谱遥感影像数据中每个栅格的特征向量，确定出每个栅格对应的目标栅格；

[0056] 基于每个栅格的特征向量和每个栅格对应的目标栅格的特征向量，计算出多光谱遥感影像数据的废水污染指数。

[0057] 其中，计算初始栅格的特征向量与其它栅格的特征向量之间的角度差，得到角度差集合，其中，所述初始栅格为所述多光谱遥感影像数据中的任意一个栅格，所述其它栅格为所述多光谱遥感影像数据中除所述目标栅格以外的栅格；

[0058] 将所述角度差集合中的最小值对应的其它栅格，确定为所述初始栅格对应的目标栅格。

[0059] 在本发明实施例中，步骤S102中已然得到了 （440）， （440）和（440），接着计算Sentinel‑2影像可见光波段的Kd值的离散和Kd‑PAR。

[0060]

[0061] 其中， 为波段 覆盖的光谱区间，为波段总数（对Sentinel‑2 MSI来说为5）。

[0062] 接着，使用主成分分析法（PCA）减少数据维度，找到能够最佳捕捉废水物质和背景差异的固有光学特性的线性组合。取一个中心矩阵（m×n），计算协方差矩阵：

[0063]

[0064] 其中，COV是一个n×n的对称矩阵。

[0065] 根据特征模量分析，对称矩阵可以被分解为：

[0066]

[0067] 其中，V是一个n×n的矩阵，它的列是表示转换反向的正交特征向量，表示如下：

[0068]

[0069] L是一个n×n的特征值对角矩阵，表示如下：

[0070]

[0071] 主成分矩阵可有中心矩阵A和正交矩阵V相乘得到：

[0072]

[0073] 对单幅影像的不同固有光学特性（IOPs）进行PCA分析，结果可以表示如下：

[0074]

[0075] 其中，等于4， 表示 行 列的二维图像（ = ×）。

[0076] 所以，每个栅格的第一主成分PC1可以表示为：

[0077] PC1 = v11 ·  d‑PAR + v12 ·  chla + v13 ·  dg + v14 ·  b‑spm

[0078] 其中，（v11，v12，v13，v14）为其特征向量。

[0079] 为每个栅格获得的第一个特征向量，采用谱角映射器(SAM)方法，计算n维测试向量和参考向量之间的角度差，也就是将污水污染物对应栅格的特征向量作为参考向量，其他栅格的特征向量作为测试向量进行比较，角度越小，相似性越高:

[0080]

[0081] 其中，是测试向量，是参考向量，代表向量维数，这里n=4。

[0082] 最后，选择相似度最高的特征向量（即角度最小的），定义固有光学特性的最佳线性组合LCI，从而，构建废水污染指数WCI并计算标准化异常：

[0083]

[0084] 其中， 是中心处理和进行线性组合后的像素值，μ是 像素值随时间变化的平均值，σ是随时间变化的标准差。

[0085] 在本发明实施例中，步骤S108包括如下步骤：

[0086] 基于主成分分析算法，对所述废水污染指数时间序列进行处理，得到所述废水污染指数时间序列的第一主成分；

[0087] 基于所述第一主成分和预设置信区间，对所述废水污染指数时间序列进行最小二乘拟合，得到拟合结果；

[0088] 基于所述拟合结果的线性趋势，确定出所述分析结果。

[0089] 在本发明实施例中，使用PCA‑t模式对WCI时间序列数据进行时空模式识别，PCA‑t模式将该数据分解为正交的空间模式和特征向量，在该模式下，时间序列上的每个图像可视为统计变量（列），空间上的像素视为观测值（行），矩阵表示如下：

[0090]

[0091] 其中，表示时间步长（即每个图像）， 表示 行 列的二维图像（ = ×））。

[0092] 具体过程如下：

[0093] 将WCI应用到时序固有光学特性数据上，生成WCI时间序列数据；

[0094] 图像做中心化处理，每个像素减去图像的均值。

[0095] 使用主成分分析法（PCA），其中，得到的协方差矩阵特征值最大的就是第一主成分模式（PC1），它可以解释数据集中大部分可变性。

[0096] 取95%的置信区间，通过对WCI时间序列的最小二乘拟合计算线性趋势，用来评价水质随时间的改善或恶化趋势。

[0097] 在本发明实施例中，根据固有光学特性指定的废水污染指数，本发明实施例能捕捉到废水团和周围水之间的空间差异，并检测到跨时间的异常情况，同时，首次引入PCA和时间分解应用于水色变量的时间序列，从而监测废水污染的时空分布，这样可以监测区域内是否存在重复的污染模式以及潜在的驱动因素。

[0098] 实施例二：

[0099] 本发明实施例还提供了一种基于多光谱影像的废水污染监测装置，该基于多光谱影像的废水污染监测装置用于执行本发明实施例上述内容所提供的基于多光谱影像的废水污染监测方法，以下是本发明实施例提供的装置的具体介绍。

[0100] 如图2所示，图2为上述基于多光谱影像的废水污染监测装置的示意图，该基于多光谱影像的废水污染监测装置包括：

[0101] 获取单元10，用于获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；

[0102] 计算单元20，用于基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；

[0103] 构建单元30，用于基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；

[0104] 分析单元40，用于对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势。

[0105] 在本发明实施例中，通过获取待监测水域的时序多光谱遥感影像数据，并利用2SeaColor解析模型分别对所述时序多光谱遥感影像数据中的各个多光谱遥感影像数据进行反演，得到各个多光谱遥感影像数据的光学特性数据；基于各个多光谱遥感影像数据的光学特性数据，计算出各个多光谱遥感影像数据的废水污染指数；基于各个多光谱遥感影像数据的废水污染指数，构建废水污染指数时间序列；对所述废水污染指数时间序列进行时空模式分析，得到分析结果，其中，所述分析结果用于表征所述待监测水域的水质变化趋势，达到了无需机器学习模型就可以对废水污染进行监测的目的，进而解决了现有的废水污染遥感监测算法的适用性较差的技术问题，从而实现了提高废水污染遥感监测算法适用性的技术效果。

[0106] 实施例三：

[0107] 本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述实施例一中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。

[0108] 参见图3，本发明实施例还提供一种电子设备100，包括：处理器50，存储器51，总线52和通信接口53，所述处理器50、通信接口53和存储器51通过总线52连接；处理器50用于执行存储器51中存储的可执行模块，例如计算机程序。

[0109] 其中，存储器51可能包含高速随机存取存储器（RAM，Random Access Memory），也可能还包括非不稳定的存储器（non‑volatile memory），例如至少一个磁盘存储器。通过至少一个通信接口53（可以是有线或者无线）实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。

[0110] 总线52可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图3中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。

[0111] 其中，存储器51用于存储程序，所述处理器50在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的装置所执行的方法可以应用于处理器50中，或者由处理器50实现。

[0112] 处理器50可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器50中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器50可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(Digital Signal Processing，简称DSP)、专用集成电路(Application Specific Integrated Circuit，简称ASIC)、现成可编程门阵列(Field‑Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器51，处理器50读取存储器51中的信息，结合其硬件完成上述方法的步骤。

[0113] 实施例四：

[0114] 本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器运行时执行上述实施例一中所述方法的步骤。

[0115] 另外，在本发明实施例的描述中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。

[0116] 在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

[0117] 所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

[0118] 另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。

[0119] 最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。