一种基于遥感影像反演的水质污染溯源系统和方法与流程

1.本发明涉及环境遥感技术领域，尤其涉及一种基于遥感影像反演的水质污染溯源系统和方法。


背景技术：

2.近年来，水污染事故频发不断。而现有的水质污染溯源主要是通过全方位式溯源法，即在河流、湖泊沿岸设置常规监测设备，监测获取实测数据，同时人工进行污染源走访排查，并针对不同的排污企业和源头建立污染样本数据库，通过实地取样与数据库进行比对进行水质污染溯源。此外，现有技术中还存在通过应用卫星的遥感数据，反演得到常规元素含量和重金属含量，对水质污染情况进行溯源。
3.但是全方位式溯源法前期布置时间长，导致进度缓慢，且难以有效投入使用；而现有卫星溯源方法完全依赖经验公式，或依赖只进行了简单线性组合后的经验公式，导致可靠性低、误差较大且不能够整合分析所有水质的污染情况，适应性差。


技术实现要素：

4.基于此，有必要针对上述技术问题，提供一种基于遥感影像反演的水质污染溯源系统和方法。
5.一种基于遥感影像反演的水质污染溯源系统，包括：遥感反演层、数据统合层和分析溯源层；所述遥感反演层与所述数据统合层连接，所述数据统合层与所述分析溯源层连接；所述遥感反演层，用于获取遥感影像数据，根据所述遥感影像数据反演地面数据，所述地面数据包括有常规水质指标和水质重金属指标；所述数据统合层，用于对所述遥感影像数据和地面数据进行哈希表映射，获取显像结果；所述分析溯源层，用于根据所述显像结果进行分项分析溯源和综合分析溯源。
6.在其中一个实施例中，所述遥感反演层包括有遥感影像获取模块、常规水质指标反演模块和水质重金属指标反演模块，所述常规水质指标反演模块和水质重金属指标反演模块并行。
7.在其中一个实施例中，所述常规水质指标反演模块包括有遥感影像获取单元、遥感影像降差单元、指标反演单元和第一显像单元；所述遥感影像获取单元用于通过有线或无线协议获取遥感影像数据；所述遥感影像降差单元用于通过遥感预处理和人机交互的方式，降低所述遥感影像数据的像差，第一遥感影像数据；所述指标反演单元用于根据第一遥感影像数据对常规水质指标进行反演；所述第一显像单元用于进行格式转换和显示反演结果；所述水质重金属指标反演模块包括有云图再验单元、公式生成单元、数据核算单元和第二显像单元；所述云图再验单元用于降低所述遥感影像数据的像差，获取第二遥感影像数据；所述公式生成单元用于生成和呈递反演公式；所述数据核算单元用于根据所述第二遥感影像数据和反演公式，对水质重金属指标进行反演，并验证反演结果的置信度；所述第二显像单元用于显示反演结果和置信度。
8.在其中一个实施例中，所述数据统合层包括遥感影像统合模块、地面数据统合模块和数据映射模块；所述遥感影像统合模块和地面数据统合模块并行，分别用于对所述遥感影像数据和地面数据进行统合；所述据映射模块用于根据哈希表对统合后的遥感影像数据和地面数据进行映射。
9.在其中一个实施例中，所述遥感影像统合模块和地面数据统合模块设置有相同结构的电路模型，分别为遥感电路模型和地面电路模型，所述遥感电路模型的输入行宽为常规水质监测指标的个数，所述地面电路模型的输入行宽为水质重金属指标的个数。
10.在其中一个实施例中，所述电路模型包括有a/d数模转换器、第一延时器、高速离散序列读写器、高速格式转换器、第二延时器、寄存器、d/a转换器和显像设备；所述a/d数模转换器接收时序脉冲，将连续信号采样为离散数列，所述连续信号为遥感影像数据或地面数据；所述第一延时器用于延迟所述离散数据进入所述高速离散序列读写器的时间；所述高速离散序列读写器用于对所述离散数据进行统合，获取统合数据；所述高速格式转换器用于将所述统合数据转换为二进制格式；所述第二延时器用于延迟将统合数据发送至所述寄存器的时间；所述寄存器用于寄存二进制格式的统合数据；所述d/a转换器用于将所述统合数据的数字量转换为模拟量；所述显像设备用于显示转换为模拟量后的统合数据。
11.在其中一个实施例中，所述数据映射模块包括有哈希表映射器和映射检查器，所述哈希表映射器与所述映射检查器连接；所述哈希表映射器用于将所述遥感影像和地面数据作为输入信号，根据哈希表对所述输入信号进行映射，获取显像结果作为输出信号；所述映射检查器用于是否存在一个输入信号对映了多个输出信号；在存在一个输入信号对映了多个输出信号时，将除第一个输出信号以外的所有输出信号反馈回哈希表进行重新映射。
12.在其中一个实施例中，所述分析溯源层包括有分项溯源分析模块和综合溯源分析模块；所述分项分析模块包括有坐标转换单元、物质分布标记单元和相关性分析单元；所述坐标转换单元用于根据欧式距离将经纬度信息转换为二维空间坐标系；所述物质浓度标记单元用于根据所述地面数据，在二维空间坐标系中标记对应水质指标的物质浓度分布；所述相关性分析单元用于根据地面数据，获取各水质指标间的相关关系；所述综合溯源分析模块用于根据所述相关关系和物质浓度分布，获取各个水质指标的溯源结果。
13.一种基于遥感影像反演的水质污染溯源方法，包括以下步骤：采用遥感反演层获取遥感影像数据，根据所述遥感影像数据反演获取地面数据，所述地面数据包括有常规水质指标和水质重金属指标；根据数据统合层对所述遥感影像数据和地面数据进行统合，并通过哈希表对统合后的遥感影像和地面数据进行映射，获取显像结果；根据分析溯源层对所述显像结果进行分项分析溯源和综合分析溯源，获取各个水质指标的溯源结果。
14.在其中一个实施例中，所述根据数据统合层对所述遥感影像数据和地面数据进行统合，并通过哈希表对统合后的遥感影像和地面数据进行映射，获取显像结果，还包括：将所述统合后的遥感影像和地面数据作为输入信号f(x)，输入电路模型中；在电路模型检测到异常时序脉冲或驱动电平不稳时向第一延时器发送中断信号，并向冲激响应h(x)发送一段测试序列f
test
(x)，根据所述测试序列经过冲激响应后的返回结果y
test
(x)判断异常种类；其中，若返回结果y
test
(x)中，x间隔不一，则为时序异常，清除当前时序脉冲信号，向所述第一延时器请求新的脉冲信号，并自动同步到当前处理的时序；若返回结果y
test
(x)与预存结果y
truth
(x)不符，则为数据异常，丢弃整组异常数据，请求从输入信号f(x)中重新获取数据；
若接收到的返回结果y
test
(x)不完整或未接受到y
test
(x)信号，则为硬件异常，追溯f
test
(x)丢失发生的位置并标记，在发送标记后中止系统，显示带有标记信息的异常。
15.相比于现有技术，本发明的优点及有益效果在于：
16.1、本发明能够在保证可靠性和准确度的同时，直接通过遥感影像数据反演获取地面数据，提升了工作效率，并通过数据统合层对遥感数据进行映射处理，获取对应的显像结果，确保数据可靠性，还能够进行分项分析溯源和综合分析溯源，整合分析所有水质的污染情况，实现全面溯源。
17.2、本发明能够通过电路模型对遥感影像数据和地面数据进行多次同步验证，确保了映射显像结果的准确性。
18.3、本发明能够广泛应用于环保检测、城市规划、应急监测、农业监管与产业布局规划等方面。
附图说明
19.图1为一个实施例中一种基于遥感影像反演的水质污染溯源系统的结构示意图；
20.图2为图1中常规水质指标反演模块的结构示意图；
21.图3为图1中水质重金属指标反演模块的结构示意图；
22.图4为一个实施例中电路模型的结构示意图；
23.图5为图1中数据映射模块的结构示意图；
24.图6为图1中分项溯源分析模块的结构示意图；
25.图7为一个实施例中一种基于遥感影像反演的水质污染溯源方法的流程示意图；
26.图8为图7中步骤s102的流程示意图。
27.附图中，遥感反演层10、遥感影像获取模块11、常规水质指标反演模块12、遥感影像获取单元121、遥感影像降差单元122、指标反演单元123、第一显像单元124、水质重金属指标反演模块13、云图再验单元131、公式生成单元132、数据核算单元133、第二显像单元134、数据统合层20、遥感影像统合模块21、地面数据统合模块22、数据映射模块23、哈希表映射器231、映射检查器232、a/d数模转换器241、第一延时器242、高速离散序列读写器243、高速格式转换器244、第二延时器245、寄存器246、d/a转换器247、显像设备248、分析溯源层30、分项溯源分析模块31、坐标转换单元311、物质分布标记单元312和相关性分析单元313、综合溯源分析模块32。
具体实施方式
28.为了使本发明的目的、技术方案及优点更加清楚明白，下面通过具体实施方式结合附图对本发明做进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
29.在一个实施例中，如图1所示，提供了一种基于遥感影像反演的水质污染溯源系统，包括：遥感反演层10、数据统合层20和分析溯源层30；遥感反演层10与数据统合层20连接，数据统合层20与分析溯源层30连接；遥感反演层10，用于获取遥感影像数据，根据遥感影像数据反演地面数据，地面数据包括有常规水质指标和水质重金属指标；数据统合层20，用于对遥感影像数据和地面数据进行哈希表映射，获取显像结果；分析溯源层30，用于根据
显像结果进行分项分析溯源和综合分析溯源。
30.在本实施例中，通过遥感反演层10获取遥感影像数据，并反演获取地面数据，包括常规水质指标和水质重金属指标；通过数据统合层20对遥感影像数据和地面数据进行哈希表映射，获取显像结果，并根据分析溯源层30对显像结果进行分项分析溯源和综合分析溯源，能够在保证可靠性和准确度的同时，直接通过遥感影像数据反演获取地面数据，提升了工作效率，并通过数据统合层20对遥感数据进行映射处理，获取对应的显像结果，确保数据可靠性，还能够进行分项分析溯源和综合分析溯源，整合分析所有水质的污染情况，实现全面溯源。
31.其中，遥感反演层10包括有遥感影像获取模块11、常规水质指标反演模块12和水质重金属指标反演模块13，常规水质指标反演模块12和水质重金属指标反演模块13并行。
32.具体地，遥感影像获取模块11可以根据现有高分系列卫星、资源系列卫星、哨兵系列卫星或探测波段位于可见光波段及红外波段的在营卫星获取，也可以从互联网上下载公开的遥感影像数据。常规水质指标反演模块12和水质重金属指标反演模块13为并行。
33.其中，常规水质指标为溶解氧、电导率、浑浊度、高锰酸盐、氨氮、总磷和总氮；水质重金属指标包括有铜离子、锌离子、六价铬离子和氰化物。常规水质指标和水质重金属指标均为水质指标。
34.如图2和图3所示，常规水质指标反演模块12包括遥感影像获取单元121、遥感影像降差单元122、指标反演单元123和第一显像单元124；遥感影像获取单元121用于通过有线或无线协议获取遥感影像数据；遥感影像降差单元122用于通过遥感预处理和人机交互的方式，降低遥感影像数据的像差，获取第一遥感影像数据；指标反演单元123用于根据第一遥感影像数据对常规水质指标进行反演；第一显像单元124用于进行格式转换和显示反演结果；水质重金属指标反演模块13包括有云图再验单元131、公式生成单元132、数据核算单元133和第二显像单元134；云图再验单元131用于降低遥感影像数据的像差，获取第二遥感影像数据；公式生成单元132用于生成和呈递反演公式；数据核算单元133用于根据第二遥感影像数据和反演公式，对水质重金属指标进行反演，并验证反演结果的置信度；第二显像单元134用于显示反演结果和置信度。
35.具体地，常规水质指标反演模块12根据遥感影像获取单元121的有线或无线协议获取遥感影像数据，遥感影像降差单元12用于通过遥感预处理，例如图像定标、图像增强、图像纠正、正射校正、图像镶嵌、数据融合、图像变换、信息提取和内容分类等方式，以及人机交互的方式降低遥感影像数据的像差，获取第一遥感影像数据，便于后续计算；指标反演单元123根据第一遥感影像数据反演获取常规水质指标含量，通过第一显像单元124进行格式转换并显示。
36.具体地，水质重金属指标反演模块13获取遥感影像数据，根据云图再验单元131，采用人机交互的方法降低遥感影像的像差，获取第二遥感影像数据，通过公式生成单元132生成反演所需的公式，数据核算单元133根据生成的最优公式反演水质重金属指标含量，并检验其置信度，最后通过第二显像单元134显示水质重金属指标和对应的置信度。
37.具体地，对遥感影像数据进行遥感预处理和遥感影像降差的过程，主要发生在遥感影像降差单元122和云图再验单元131；根据处理后的遥感影像数据反演出地面数据，即常规水质指标和水质重金属指标，主要通过指标反演单元123和数据核算单元133完成。
38.其中，数据统合层20包括遥感影像统合模块21、地面数据统合模块22和数据映射模块23；遥感影像统合模块21和地面数据统合模块22并行，分别用于对遥感影像数据和地面数据进行统合；数据映射模块23用于根据哈希表对统合后的遥感影像数据和地面数据进行映射。
39.具体地，遥感影像统合模块21和地面数据统合模块22并行，分别对遥感影像数据和地面数据进行统合，并通过数据映射模块23采用哈希表对统合后的遥感影像数据和地面数据进行映射，获取的显像结果发送至分析溯源层30。
40.其中，遥感影像统合模块21和地面数据统合模块22设置有相同结构的电路模型，分别为遥感电路模型和地面电路模型，遥感电路模型的输入行宽为常规水质监测指标的个数，地面电路模型的输入行宽为水质重金属指标的个数。
41.具体地，遥感电路模型和地面电路模型的输入行宽不同，遥感电路模型的输入行宽为常规水质监测指标的个数，为7；地面电路模型的输入行宽为水质重金属指标的个数，为4。
42.如图4所示，电路模型包括有a/d数模转换器241、第一延时器242、高速离散序列读写器243、高速格式转换器244、第二延时器245、寄存器246、d/a转换器247和显像设备248；a/d数模转换器241接收时序脉冲，将连续信号采样为离散数列，连续信号为遥感影像数据或地面数据；第一延时器242用于延迟离散数据进入高速离散序列读写器243的时间；高速离散序列读写器243用于对离散数据进行统合，获取统合数据；高速格式转换器244用于将统合数据转换为二进制格式；第二延时器245用于延迟将统合数据发送至寄存器246的时间；寄存器246用于寄存二进制格式的统合数据；d/a转换器247用于将统合数据的数字量转换为模拟量；显像设备248用于显示转换为模拟量后的统合数据。
43.具体地，电路模型接受一个时序脉冲t(x)，在该时序脉冲的作用下，a/d数据转换模块将连续输入信号f(t)，按t采样为离散数据f(x)，在第一延时器242的作用下，电路模型能够有足够的时间处理f(x
‑
τ)的数据，而不用担心数据丢失和时序冲突的情况；离散输入信号传入高速离散序列读写器243，在第一延时器242的作用下，高速离散序列读写器243能够在τ
–
a的时间内，将离散输入信号传递给高速格式转换器244，在第二延时器245的作用下，高速格式转换器246在τ
–
b的时间内将离散输入信号转换为寄存器247能够存储的二进制格式，为了满足实时性和避免冲突，有：
44.τ
‑
a τ
‑
b＝τ；
ꢀꢀꢀ
(1)
45.τ＝a b；
ꢀꢀꢀ
(2)
46.其中，a为高速离散序列读写器243的弛豫时间，b为高速格式转换器244的弛豫时间。
47.根据上式(2)，选择弛豫时间和小于τ的两种元器件，搭建遥感影像统合模块21和地面数据统合模块22。
48.在第二延时器245接受一个高电平信号后，将时间n内的信号序列呈递给寄存器247，之后自动触发信号调整为低电平触发模型，一般地，需满足n≤τ，且最佳工作点为n≈τ；然后通过d/a转换器247将数字信号转换为模拟信号，并通过显像设备248进行显示。
49.如图5所示，数据映射模块23包括有哈希表映射器231和映射检查器232，哈希表映射器231与映射检查器232连接；哈希表映射器231用于将遥感影像和地面数据作为输入信
号，根据哈希表对输入信号进行映射，获取显像结果作为输出信号；映射检查器232用于是否存在一个输入信号对映了多个输出信号；在存在一个输入信号对映了多个输出信号时，将除第一个输出信号以外的所有输出信号反馈回哈希表进行重新映射。
50.具体地，由于哈希表映射的方法可能出现重映射，因此设置映射检查器232。在通过哈希表映射器231对输入信号映射完成后，通过映射检查器232检查是否存在一个输入信号对映了多个输出信号，若存在，则将第一输出信号以外的所有输出信号反馈回哈希表映射器231，进行重新映射。由于采用的是反馈形式，可能会出现向一个输入查询映射关系得到深度较大的情况，提供的映射函数悦合理，这种情况出现的概率越小。
51.如图6所示，分析溯源层30包括有分项溯源分析模块31和综合溯源分析模块32；分项分析模块31包括有坐标转换单元311、物质分布标记单元312和相关性分析单元313；坐标转换单元311用于根据欧式距离将经纬度信息转换为二维空间坐标系；物质浓度标记单元312用于根据地面数据，在二维空间坐标系中标记对应水质指标的物质浓度分布；相关性分析单元313用于根据地面数据，获取各水质指标间的相关关系；综合溯源分析模块32用于根据相关关系和物质浓度分布，获取各个水质指标的溯源结果。
52.具体地，分项溯源分析模块31受理论分析法和影像结合法的支持，综合溯源分析模块32受综合分析法的支持。
53.例如，在选取共2973组遥感影像数据后，通过上述特征对遥感影像数据进行地面数据反演，根据坐标转换单元311采用欧式距离将经纬度信息转换为二维空间坐标系，在遥感地图上获得各水质指标的物质浓度分布，并通过物质浓度标记单元312进行对应的标记，采用相关性分析单元313对地面数据进行皮尔逊相关性分析，获得各个水质指标间的相关关系，如表1所示，为11项水质指标(包括7中常规水质指标和4中水质重金属指标)的皮尔逊相关性表。
54.表1 水质指标的皮尔逊相关性表
[0055][0056]
在本实施例中，分项溯源分析模块31的溯源分析具体包括：
[0057]
其中，关于溶解氧：水中溶解氧水温、空气含氧量、水生植物光合作用有关，一般地，在20℃、大气压强100kpa条件下，纯水中溶解氧大约为9mg/l，当水中溶解氧小于5mg/l时，水质质量变差。在表1中，与溶解氧相关性显著的指标共九种(按强弱排序)：氰化物、氨氮、总氮、总磷、铜离子、锌离子、高锰酸盐、浑浊度、六价铬离子。可知，除电解率外，其他水
质指标均对溶解氧有较大影响，其中影响最大的为氰化物。
[0058]
因此，造成水质溶解氧下降的溯源方向为：(1)水温变化大；(2)空气含氧量低；(3)水生植物光合作用差；(4)与氰化物排放相关的因素。
[0059]
其中，关于电解率：水中导电性越好，电解率越大，它与温度和掺杂浓度有很大的关系，一般地，纯水的导电率≤2μs/cm，电解率变高，水质质量变差。在表1中，与电解率相关性显著的指标共十种(按强弱排序)：浑浊度、六价铬离子、铜离子、锌离子、总磷、高锰酸盐、氰化物、总氮、氨氮。可知，其他水质指标对溶解氧有较大影响，其中影响最大的为浑浊度。
[0060]
因此，造成水质电解率上升的溯源方向为：(1)水温变化大；(2)掺杂浓度高；(3)与水体浑浊程度相关的因素。
[0061]
其中，关于浑浊度：水中浑浊度表示水中含有悬浮物和胶体物质而呈浑浊状态，造成通过水的光线被散射或光线透过时受到阻碍的程度。水的浑浊度主要由掺杂引起，包括泥沙、粘土、水生生物、细菌、病毒、高分子有机物等。一般地，以sio2计，蒸馏水中浑浊度不高于1mg/l，浑浊度变高，水质质量变差。在表1中，与浑浊度相关性显著的指标共九种(按强弱排序)：六价铬离子、铜离子、锌离子、总氮、氰化物、电解率、氨氮、高锰酸盐、总磷。可知，其他水质指标对溶解氧有较大影响，其中影响最大的为六价铬离子。
[0062]
因此，造成水质浑浊度上升的溯源方向：(1)掺杂浓度高；(2)与六价铬离子排放相关的因素。
[0063]
其中，关于高锰酸盐：水中高锰酸盐指数不能作为理论需氧量或总有机物含量的指标，标准水质中，高锰酸盐含量应总小于溶解氧量。高锰酸盐变高，水质质量变差。在表1中，与高锰酸盐相关性显著的指标共九种(按强弱排序)：铜离子、锌离子、六价铬离子、电解率、浑浊度、总磷、总氮、氰化物、溶解氧。可知，除氨氮外，其他水质指标均对溶解氧有较大影响，其中影响最大的为铜离子。
[0064]
因此，造成水质高锰酸盐上升的溯源方向：(1)掺杂浓度高；(2)与铜离子排放相关的因素。
[0065]
其中，关于氨氮：水中氨氮是水体中以游离氨和铵离子形式存在的氮，是水体的营养素，可导致水富营养化现象产生，与水中藻类、浮游生物等含量有关，水中氨氮变高，水质质量变差。在表1中，与氨氮相关性显著的指标共八种(按强弱排序)：总氮、浑浊度、溶解氧、氰化物、锌离子、铜离子、电解率、六价铬离子。可知，除高锰酸盐和总磷外，其他水质指标均对溶解氧有较大影响，其中影响最大的为总氮。
[0066]
因此，造成水质氨氮上升的溯源方向为：(1)水中藻类、浮游生物等含量升高；(2)与总氨相关的因素。
[0067]
其中，关于总磷：水中总磷是水样经消解后将各种形态的磷转变为正磷酸盐后的测定结果，其主要来源为生活污水、化肥、有机磷农药和近代洗涤剂排放，水中藻类生物依赖总磷生长，过量磷会导致水体富营养化。一般地，水质中总磷含量不高于0.05mg/l，总磷变高，水质质量变差。在表1中，与总磷相关性显著的指标共九种(按强弱排序)：铜离子、六价铬离子、锌离子、电解率、浑浊度、高锰酸盐、溶解氧、氰化物、总氮。可知，除氨氮外，其他水质指标均对溶解氧有较大影响，其中影响最大的为铜离子。
[0068]
因此，造成水质总磷上升的溯源方向为：(1)生活污水排放；(2)农业化肥排放；(3)农业农药排放；(4)与铜离子排放相关的因素。
[0069]
其中，关于总氮：水中总磷是水中各种形态无机和有机氮的总和，常被用于表示水体受营养物质污染的程度。一般地，水质中总氮含量不高于0.05mg/l，总氮变高，水质质量变差。在表1中，与总氮相关性显著的指标共十种(按强弱排序)：氨氮、锌离子、铜离子、六价铬离子、浑浊度、溶解氧、高锰酸盐、电解率、总磷、氰化物。可知，其他水质指标对溶解氧有较大影响，其中影响最大的为氨氮。
[0070]
因此，造成水质总氮上升的溯源方向：(1)水体营养物质变多；(2)与氨氮排放相关的因素。
[0071]
其中，关于水中重金属离子：与各种离子相关性显著的指标(按强弱排序)如下：
[0072]
铜离子：锌离子、六价铬离子、浑浊度、总氮、高锰酸盐、总磷、电解率、氰化物、溶解氧、氨氮。
[0073]
锌离子：铜离子、六价铬离子、浑浊度、总氮、高锰酸盐、电解率、总磷、氰化物、溶解氧、氨氮。
[0074]
六价铬离子：锌离子、铜离子、浑浊度、氰化物、总氮、高锰酸盐、电解率、总磷、溶解氧、氨氮。
[0075]
氰化物：六价铬离子、浑浊度、溶解氧、铜离子、电解率、氨氮、锌离子、总磷、高锰酸盐、总氮。
[0076]
因此，造成水质重金属元素上升的溯源方向：(1)化工污水排放；(2)土壤污染严重；(3)与浑浊度相关的因素；(4)与溶解氧相关的因素。
[0077]
在本实施例中，综合溯源分析模块32根据上述分项溯源分析模块31获取的溯源结果，获取的溯源方向如下：
[0078]
溶解氧低：水温变化大、空气含氧量低、水生植物光合作用差、与氰化物排放相关的因素。
[0079]
电解率高：水温变化大、掺杂浓度高、与水体浑浊程度相关的因素。
[0080]
浑浊度高：掺杂浓度高、与水六价铬离子排放相关的因素。
[0081]
高锰酸盐高：掺杂浓度高、与铜离子排放相关的因素。
[0082]
氨氮高：水中藻类、浮游生物等含量升高、与总氮相关的因素。
[0083]
总磷高：生活污水排放、农业化肥排放、农业农药排放、与铜离子排放相关的因素。
[0084]
总氮高：水体营养物质变多、与氨氮排放相关的因素。
[0085]
水中重金属离子高：化工污水排放、土壤污染严重、与浑浊度相关的因素、与溶解氧相关的因素。
[0086]
采用综合分析法替换上述“与*相关的因素”，过程如下：
[0087]
对氰化物受污染原因溯源，追溯到水中重金属排放异常；对水体浑浊程度加深原因溯源，追溯到掺杂浓度变高、水中重金属排放异常；对铜离子浓度增大原因进行溯源，追溯到水中重金属排放异常；对总氨浓度增大原因进行溯源，追溯到水体富营养化、水中藻类生物滋生、水中浮游生物含量升高。
[0088]
在一个实施例中，如图7所示，提供了一种基于遥感影像反演的水质污染溯源方法，包括以下步骤：
[0089]
步骤s701，采用遥感反演层获取遥感影像数据，根据遥感影像数据反演获取地面数据，地面数据包括有常规水质指标和水质重金属指标。
[0090]
步骤s702，根据数据统合层对遥感影像数据和地面数据进行统合，并通过哈希表对统合后的遥感影像和地面数据进行映射，获取显像结果。
[0091]
步骤s703，根据分析溯源层对显像结果进行分项分析溯源和综合分析溯源，获取各个水质指标的溯源结果。
[0092]
在本实施例中，采用遥感影像反演层获取遥感影像数据，并根据遥感影像数据反演获取地面数据的常规水质指标和水质重金属指标，根据数据统合层对遥感影像数据和地面数据进行统合，并通过哈希表对统合后的遥感影像和地面数据进行映射，获取显像结果，最后采用分析溯源层对显像结果进行分项分析溯源和综合分析溯源，获取各个水质指标的溯源结果，能够在保证可靠性和准确度的同时，直接通过遥感影像数据反演获取地面数据，提升了工作效率，并对遥感数据进行映射处理，获取对应的显像结果，确保数据可靠性，还能够进行分项分析溯源和综合分析溯源，整合分析所有水质的污染情况，实现全面溯源。
[0093]
如图8所示，步骤s102还包括：将统合后的遥感影像和地面数据作为输入信号f(x)，输入电路模型中；在电路模型检测到异常时序脉冲或驱动电平不稳时向第一延时器发送中断信号，并向冲激响应h(x)发送一段测试序列f
test
(x)，根据测试序列经过冲激响应后的返回结果y
test
(x)判断异常种类；其中，若返回结果y
test
(x)中，x间隔不一，则为时序异常，清除当前时序脉冲信号，向第一延时器请求新的脉冲信号，并自动同步到当前处理的时序；若返回结果y
test
(x)与预存结果y
truth
(x)不符，则为数据异常，丢弃整组异常数据，请求从输入信号f(x)中重新获取数据；若接收到的返回结果y
test
(x)不完整或未接受到y
test
(x)信号，则为硬件异常，追溯f
test
(x)丢失发生的位置并标记，在发送标记后中止系统，显示带有标记信息的异常。
[0094]
具体地，通过上述时序同步算法，增强电路模型的稳定性；且能够保障两种异常的完全恢复：时序异常和数据异常，保障一种异常的精确定位：硬件异常，在遇到未预料的异常时，将提示并停止统合。
[0095]
在一个实施例中，以2021年7月27日的某市某区的遥感断面数据(之后称为实验数据)为例，分析7种常规水质指标和4种水质重金属指标的污染源。
[0096]
具体实施步骤如下：
[0097]
步骤s701，获取到哨兵二号的遥感波段数据如下：
[0098]
表2 下载的断面数据(部分)
[0099]
band1band2band3band4band5band6band7band8band8aband9band11band121167130016841874174777780660446921424019012221374174219361786795800590450227227181115913461706192017447517595554121302191731308144017541926181588490972053842825518213731480183019841867940100273163857030120713881566186820381899102211148727386313582631128136617581982187088992368155726331525510851294167818861812100110718667473214553711027118612961186128111671252121210702460954733
[0100]
在步骤s701中，获取到的反演数据如下表3所示：
[0101]
表3 反演得到的11种水质指标(单位mg/l)(部分)
[0102]
o2elentumnntptncuzncr6cn5.41960.78388.185.720.320.036.020.000.060.000.016.011211.83370.485.980.330.026.060.000.060.000.016.531078.04403.046.110.410.028.560.000.060.000.013.191458.18225.075.490.010.056.070.000.070.000.012.27800.02141.333.690.020.025.880.000.080.000.010.714598.2899.593.230.100.024.890.000.080.000.010.03201.78376.125.040.320.016.760.000.070.000.014.39263.43321.964.310.030.094.670.000.090.000.0112.53381.6984.423.710.110.123.580.000.110.000.00
[0103]
注：反演出的11种水质指标数据的小数点位数并非只有2位，此处旨在方便展示。
[0104]
步骤s702，对遥感影像数据和地面数据进行统合，并通过哈希表对统合后的遥感影像和地面数据进行映射，获取显像结果，在遥感地图上显示对应水质指标的物质浓度分布。
[0105]
步骤s703，通过上述两个步骤的分析，发现某市某区清溪河水域的重金属和常规水质指标(主要是tp、ntu)异常，结论为：
[0106]
(1)水质中浑浊度增加，可能是泥沙流入、水位增加、其他河流汇入造成。
[0107]
(2)水质中总磷增加，可能是水中藻类生物和浮游生物增长。
[0108]
(3)需勘察流域附近是否有异常排水、堤岸降低等情况。
[0109]
进一步地，通过工程人员实地考察发现，清溪河部分水域正处于退潮期，岸上泥沙被卷覆到水面上，引起水面浑浊度增加和水体富营养化，同时，在勘察点附近坐落有发电站，可能会有异常排水现象，其余水域正常，如治理得当，该现象应在几周到一月左右会逐渐消失。
[0110]
本领域普通技术人员可以理解实现上述实施例所述方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(read
‑
only memory，rom)或随机存储记忆体(random access memory，ram)等。
[0111]
显然，本领域的技术人员应该明白，上述本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在计算机存储介质(rom/ram、磁碟、光盘)中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。所以，本发明不限制于任何特定的硬件和软件结合。
[0112]
以上内容是结合具体的实施方式对本发明所做的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护
范围。