CN113433075A - 一种简单的黑臭水体遥感快速监测方法 - Google Patents

Abstract

本发明涉及一种黑臭水体遥感监测方法、装置、电子设备及介质，方法包括：获取遥感影像数据；使用水体提取模型对所述遥感影像数据进行处理得到水体指数和植被指数，其中所述水体提取模型包括改进的归一化差分水体指数MNDWI算法、归一化差分水体指数NDWI算法和归一化差分植被指数NDVI算法；基于所述水体指数和植被指数，判断是否存在水体；使用水体类型识别模型对所述水体的水质进行识别。本发明的实施例通过黑臭水体提取模型，有效的获得了的黑臭水体信息遥感影像数据，能够较完善地利用光谱差异特征。

Description

一种简单的黑臭水体遥感快速监测方法

技术领域

本发明涉及水体质量监测技术领域，尤其涉及一种黑臭水体遥感 监测方法、装置、电子设备及存储介质。

背景技术

现有的黑臭水体识别模型方法往往仅使用单一的模型进行分析， 然而黑臭水体由于污染物成分的差异，河道下垫面构成的差异，其光 谱特征具有不确定性，黑臭水体面积小、河段破碎，其遥感信息容易 受到河岸绿化树木、建筑物等周边环境的影响，因此现有的方法常常 具有一定的局限性，在快速监测不同地理条件和不同水质的黑臭水体 时存在亟待解决的问题。

现有的黑臭水体识别模型方法往往仅使用单一的模型进行分析， 然而黑臭水体由于污染物成分的差异，河道下垫面构成的差异，其光 谱特征具有不确定性，黑臭水体面积小、河段破碎，其遥感信息容易 受到河岸绿化树木、建筑物等周边环境的影响，因此现有的方法常常 具有一定的局限性，在快速监测不同地理条件和不同水质的黑臭水体 时存在亟待解决的问题。

发明内容

本发明的目的是提供一种黑臭水体遥感监测方法、装置、电子设 备及存储介质，用以解决现有技术中的问题。具体地，本发明实施例 提供了以下技术方案：

第一方面，本发明实施例提供了一种黑臭水体遥感监测方法，包 括：

获取遥感影像数据；

使用水体提取模型对所述遥感影像数据进行处理得到水体指数 和植被指数，其中所述水体提取模型包括改进的归一化差分水体指数 MNDWI算法、归一化差分水体指数NDWI算法和归一化差分植被指 数NDVI算法；

基于所述水体指数和植被指数，判断是否存在水体；

使用水体类型识别模型对所述水体的水质进行识别。

可选地，所述基于所述水体指数和植被指数，判断是否存在水体 包括：

在植被较为稀疏的地区，当MNDWI_t0≥a₁成立时，判断存在水体， 其中MNDWI_t0是t₀时刻的MNDWI值，a₁是预设阈值；

或，在植被较为茂盛的地区，当MNDWI_ti＜a₁，NDVI_ti≥a₂， MNDWI_ti-j≥a₁并且NDVI_ti-j＜a₂成立时，判断存在水体，其中MNDWI_ti， MNDWI_ti-j，NDVI_ti和NDVI_ti-j分别是t_i、t_i-j时刻的MNDWI值和NDVI 值，j<i，a₂是预设阈值。

可选地，所述使用水体类型识别模型对所述水体的水质进行识别 包括：

对于非细小水体，建立包括黑臭水体子模型1、黑臭水体子模型 2和黑臭水体变异系数的水体类型识别模型，其中黑臭水体子模型1 的算法为：

其中R_g是绿波段反射率，R_r是红波段反射率；

黑臭水体子模型2的算法为：

其中R_g是绿波段反射率，R_b是蓝波段反射率；

黑臭水体变异系数的算法为：

其中

R_i为由R_g、R_r和R_b组成的集合中的元素。

可选地，当HCS₁≥a₃、HCS₂≥a₄并且HCC≥a₅成立时，判断所述 水体是正常水体，其中a₃、a₄和a₅是预设阈值；

或，当a₆≤HCS₁＜a₃、a₇≤HCS₂＜a₄并且a₈≤HCC＜a₅成立时， 判断所述水体是轻度黑臭水体，其中a₆、a₇和a₈是预设阈值；

或，当HCS₁＜a₉、HCS₂＜a₁₀并且HCC＜a₁₁成立时，判断所述水体 是严重黑臭水体，其中a₉、a₁₀和a₁₁是预设阈值。

可选地，所述使用水体类型识别模型对所述水体的水质进行识别 包括：

对于细小水体，建立包括高程变异系数的水体类型识别模型，其 中高程变异系数的算法为：

其中H_i为预定规格的网格高程、

为 研究区预定规格的网格高程均值，n为预定规格的矢量网格覆盖的数 字高程模型栅格数量。

可选地，当L≤a₁₂，|S_i-S_j|≥a₁₃并且DXQF≤a₁₄成立时，判断 所述水体是黑臭水体，其中L为细小水体的边界长度，a₁₂为预设阈 值；S_i和S_j为距离最近的两个细小斑块的中心点，a₁₃为相邻两个斑 块的最短距离，a₁₄为预设阈值。

第二方面，本发明实施例还提供了一种黑臭水体遥感监测装置， 包括：

数据获取模块，用于获取遥感影像数据；

水体提取模块，用于使用水体提取模型对所述遥感影像数据进行 处理得到水体指数和植被指数，其中所述水体提取模型包括改进的归 一化差分水体指数MNDWI算法、归一化差分水体指数NDWI算法和 归一化差分植被指数NDVI算法；

水体判断模块，用于基于所述水体指数和植被指数，判断是否存 在水体；

水质识别模块，用于使用水体类型识别模型对所述水体的水质进 行识别。

第三方面，本发明实施例还提供了一种电子设备，包括存储器、 处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征 在于，所述处理器执行所述程序时实现如第一方面所述黑臭水体遥感 监测方法的步骤。

第四方面，本发明实施例还提供了一种非暂态计算机可读存储介 质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执 行时实现如第一方面所述黑臭水体遥感监测方法的步骤。

第五方面，本发明实施例还提供了一种计算机程序产品，所计算 机程序产品包括有计算机程序，该计算机程序被处理器执行时实现如 第一方面所述黑臭水体遥感监测方法的步骤。

附图说明

图1为本发明一实施例提供的黑臭水体遥感监测方法的流程图；

图2为本发明一实施例提供的黑臭水体遥感监测装置的结构示意 图；以及

图3为本发明一实施例提供的电子设备的结构示意图。

具体实施方式

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

在本发明提供的实施例中，可以通过各种渠道获取各种格式的遥 感数据，如从国家卫星数据中心，商用卫星数据中心等处，获取如高 分2号、资源3号和哨兵2号等等的卫星采集到的影像数据，获取到的 影像数据可以包括多个波段通道，可以包括数字高程模型(Digital Elevation Model，DEM)数据。本发明实施例提供的方法可以用于水 系调查，生态恢复，土地资源利用调查等等。

在本发明提供的实施例中，可以对获取到的影像数据进行各种预 处理，主要是对数字图像处理，包括但不限于用于校正在成像、记录、 传输或回放过程中引入的数据错误、噪声与畸变的图像恢复手段，如 辐射校正、几何校正等；用于方便传输、储存和提高数据处理效率的 数据压缩手段；用于去模糊、提高影像质量的影像增强手段等。

对获取到的影像数据的预处理还可以包括大气校正，遥感所利用 的各种辐射能均要与地球大气层发生相互作用：或散射、或吸收，从 而使能量衰减，并使光谱分布发生变化。有时必须对遥感数据进行大 气校正，从水体中提取生物物理变量(如水体中的叶绿素、悬浮泥沙、 温度)时，就必须对遥感数据进行大气校正。如果数据未经校正，就 可能会丢失这些重要成分的反射率的微小差别信息。此外，如果需要 将某影像中提取的物理量与另一不同时相影像中提取的同一物理量 相比较，就必须对遥感数据进行大气校正。例如，可以使用完整的遥 感图像处理平台(The Environment for Visualizing Images，ENVI)中 的FLAASH工具进行辐射定标和大气校正，从而得到较精准的反射率 数据。

参考图1，图1为本发明一实施例提供的黑臭水体遥感监测方法的 流程图，该黑臭水体遥感监测方法包括如下步骤：

步骤110：获取遥感影像数据；

步骤120：使用水体提取模型对所述遥感影像数据进行处理得到 水体指数和植被指数，其中所述水体提取模型包括改进的归一化差分 水体指数MNDWI算法、归一化差分水体指数NDWI算法和归一化差 分植被指数NDVI算法；

步骤130：基于所述水体指数和植被指数，判断是否存在水体；

步骤140：使用水体类型识别模型对所述水体的水质进行识别。

根据本发明的实施例，构建水体提取模型。目前使用遥感影像提 取水体的算法很多，在本发明提供的实施例中，采用水体指数算法， 具体地，首先计算水体指数和植被指数，对于超过4个波段的遥感影 像，可以选择改进的归一化差分水体指数(ModifiedNormalized Difference Water Index，MNDWI)算法或者归一化差分水体指数(Normalized Difference Water Index，NDWI)以得到水体指数，并采 用归一化差分植被指数(Normalized Difference Vegetation Index, NDVI)以得到植被指数。具体地，MNDWI参照下式计算：

其中R_g是绿波段反射率，R_swir是短波红外波段反射率。

NDWI参照下式计算：

其中R_g是绿波段反射率，R_nir为近红外波段反射率。NDVI参照下式 计算：

其中R_r是红波段反射率。对于只有四个波段的遥感影像，选择NDWI 算法，使用NDWI值代替MNDWI值。在本发明的实施例中，通过 采用包括两种水体指数算法的水体提取模型，较好地适应了获取到 的遥感影像数据，能够较完善地利用光谱特征数据。

接下来在本实施例中，基于上面计算得到的水体指数和植被指 数，对河道或者水系中是否有水进行判断。在植被较为稀疏的地区， 有下式

MNDWI_t0≥a₁ (4)当式(4)为真时，判断该水系中有水，其中MNDWI_t0是t₀时刻的MNDWI 值，a₁是预设阈值，可以通过多种方法确定合适的阈值，如根据具 体数据通过统计学方法判断出合适阈值，阈值也可以是根据具体情 况而变动的。

在植被较为茂盛的地区，需要考虑植被的影响，因此有

MNDWI_ti＜a₁，NDVI_ti≥a₂，MNDWI_ti-j≥a₁，NDVI_ti-j＜a₂ (5)

当式(5)中各条件均为真时，判断该水系中有水，其中MNDWI_ti， MNDWI_ti-j，NDVI_ti和NDVI_ti-j分别是t_i、t_i-j时刻的MNDWI值和NDVI 值，其中j<i，t_i-j为t_i时刻之前的时刻，其时间间隔根据具体情况而 设置，a₂是预设阈值。通过上述判断，本实施例监测到了遥感影像 中存在水的河道或水系，接下来对其中水的质量进行判断，即监测 是否为黑臭水体。

在当前步骤中，首先需要建立水体类型识别模型，其中涉及到 三个相对重要的水体识别算法，分别为

其中，

R_g是绿波段反射率，R_r是红波段反射率，R_b是蓝波段反射率，R_i为由R_g、R_r和R_b组成的集合中的元素。HCS₁表 示黑臭水体子模型1，HCS₂表示黑臭水体子模型2，HCC表示黑臭 水体变异系数。通过构建包括两个黑臭水体子模型和黑臭水体变异 系数的水体类型识别模型，在本发明的实施例中，考虑到了由于污 染物成分的差异和河道下垫面构成的差异所造成的光谱特征的不确 定性，能够更精准地判断具有不同污染物的水体的水质状况。

接下来使用上述水体类型识别模型，对水系或河道中水体的水 质进行识别，在本发明实施例中，可以将水质级别分为三类，正常 水体、轻度黑臭水体和严重黑臭水体。其中，当满足下式时，

HCS₁≥a₃

HCS₂≥a₄

HCC≥a₅ (9)

该水体是正常水体，其中a₃、a₄和a₅是预设阈值，可以通过多种方 法确定合适的阈值，如根据具体数据通过统计学方法判断出合适阈 值，阈值也可以是根据具体情况而变动的。

当满足下式时，

a₆≤HCS₁＜a₃

a₇≤HCS₂＜a₄

a₈≤HCC＜a₅ (10)

该水体是轻度黑臭水体，其中a₆、a₇和a₈是预设阈值。

当满足下式时，

HCS₁＜a₉

HCS₂＜a₁₀

HCC＜a₁₁ (11)

该水体是严重黑臭水体，其中a₉、a₁₀和a₁₁是预设阈值。

至此，可以根据本发明的实施例，完成对一般尺寸的河道和水 系的水质监测。根据本发明的实施例，还提供了基于遥感技术的对 细小水体的水质监测方法。在本文中提到的细小水体是相对于遥感 图像的像元(即影像单元，是遥感成像过程中的采样点)而言的，当监测到的水体为混合像元时，即一个像元中还包括了除水体特征 以外的其他地表特征时，则判定该水体为细小水体。

当根据遥感数据获得细小水体的空间分布信息时，通过上述的 水体提取模型提取水体分布信息，并矢量化，水体由于具有水体光 谱特征的共性，河流与地表背景之间光谱反差明显，边界清晰，局 部区域内的河流水体光谱相似性强，易聚集成斑块状，因此会得到 矢量化的斑块，然后开展如下判定：

(a)L≤a₁₂

(b)|S_i-S_j|≥a₁₃

(c)

DXQF≤a₁₄

式中，L为细小水体的边界长度，a₁₂为预设阈值；S_i和S_j为距 离最近的两个细小斑块的中心点，a₁₃为相邻两个斑块的最短距离； DXQF为高程变异系数，H_i为预定规格的网格高程、H为研究区预定 规格的网格高程均值，n为预定规格的矢量网格覆盖的数字高程模型栅格数量，a₁₄为预设阈值。当满足上述判定时，该细小水体为黑 臭水体。根据本发明的实施例，通过细小水体建立水质识别模型， 解决了黑臭水体面积小、河段破碎，其遥感信息容易受到河岸绿化 树木、建筑物等周边环境的影响的问题，能够较快速地实现对复杂 地理条件下细小水体的水质监测。

图2为本发明一实施例提供的黑臭水体遥感监测装置的结构示 意图，参照图2，本实施例提供的黑臭水体遥感监测装置包括：

数据获取模块210，用于获取遥感影像数据；

水体提取模块220，用于使用水体提取模型对所述遥感影像数据 进行处理得到水体指数和植被指数，其中所述水体提取模型包括改进 的归一化差分水体指数MNDWI算法、归一化差分水体指数NDWI算 法和归一化差分植被指数NDVI算法；

水体判断模块230，用于基于所述水体指数和植被指数，判断是 否存在水体；

水质识别模块240，用于使用水体类型识别模型对所述水体的水 质进行识别。

由于本发明实施例提供的黑臭水体遥感监测装置，可以用于执 行上述实施例所述的黑臭水体遥感监测方法，其工作原理和有益效 果类似，故此处不再详述，具体内容可参见上述实施例的介绍。

在本实施例中，需要说明的是，本发明实施例的装置中的各个 模块可以集成于一体，也可以分离部署。上述模块可以合并为一个 模块，也可以进一步拆分成多个子模块。

图3示例了一种电子设备的实体结构示意图，如图3所示，该 电子设备可以包括：处理器(processor)310、通信接口(Communications Interface)320、存储器(memory)330和通信总线340，其中，处理器 310，通信接口320，存储器330通过通信总线340完成相互间的通 信。处理器310可以调用存储器330中的逻辑指令，以执行黑臭水 体遥感监测方法，该方法包括：获取遥感影像数据；使用水体提取 模型对所述遥感影像数据进行处理得到水体指数和植被指数，其中 所述水体提取模型包括改进的归一化差分水体指数MNDWI算法、 归一化差分水体指数NDWI算法和归一化差分植被指数NDVI算法； 基于所述水体指数和植被指数，判断是否存在水体；使用水体类型 识别模型对所述水体的水质进行识别。

此外，上述的存储器330中的逻辑指令可以通过软件功能单元 的形式实现并作为独立的产品销售或使用时，可以存储在一个计算 机可读取存储介质中。基于这样的理解，本发明的技术方案本质上 或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软 件产品的形式体现出来，该计算机软件产品存储在一个存储介质中， 包括若干指令用以使得一台计算机设备(可以是个人计算机，服务 器，或者网络设备等)执行本发明各个实施例所述方法的全部或部 分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

另一方面，本发明还提供一种计算机程序产品，所述计算机程 序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所 述计算机程序包括程序指令，当所述程序指令被计算机执行时，计 算机能够执行上述各方法所提供的黑臭水体遥感监测方法，该方法 包括：获取遥感影像数据；使用水体提取模型对所述遥感影像数据 进行处理得到水体指数和植被指数，其中所述水体提取模型包括改 进的归一化差分水体指数MNDWI算法、归一化差分水体指数NDWI 算法和归一化差分植被指数NDVI算法；基于所述水体指数和植被指数，判断是否存在水体；使用水体类型识别模型对所述水体的水 质进行识别。

又一方面，本发明还提供一种非暂态计算机可读存储介质，其 上存储有计算机程序，该计算机程序被处理器执行时实现以执行上 述各提供的黑臭水体遥感监测方法，该方法包括：获取遥感影像数 据；使用水体提取模型对所述遥感影像数据进行处理得到水体指数 和植被指数，其中所述水体提取模型包括改进的归一化差分水体指 数MNDWI算法、归一化差分水体指数NDWI算法和归一化差分植 被指数NDVI算法；基于所述水体指数和植被指数，判断是否存在 水体；使用水体类型识别模型对所述水体的水质进行识别。

以上所描述的装置实施例仅仅是示意性的，其中所述作为分离 部件说明的单元可以是或者也可以不是物理上分开的，作为单元显 示的部件可以是或者也可以不是物理单元，即可以位于一个地方， 或者也可以分布到多个网络单元上。可以根据实际的需要选择其中 的部分或者全部模块来实现本实施例方案的目的。本领域普通技术 人员在不付出创造性的劳动的情况下，即可以理解并实施。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了 解到各实施方式可借助软件加必需的通用硬件平台的方式来实现， 当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者 说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该 计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、 磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个 人计算机，服务器，或者网络设备等)执行各个实施例或者实施例 的某些部分所述的方法。

最后应说明的是：以上实施例仅用以说明本发明的技术方案， 而非对其限制；尽管参照前述实施例对本发明进行了详细的说明， 本领域的普通技术人员应当理解：其依然可以对前述各实施例所记 载的技术方案进行修改，或者对其中部分技术特征进行等同替换； 而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实 施例技术方案的精神和范围。

Claims (10)

1.一种黑臭水体遥感监测方法，其特征在于，包括：

获取遥感影像数据；

使用水体提取模型对所述遥感影像数据进行处理得到水体指数和植被指数，其中所述水体提取模型包括改进的归一化差分水体指数MNDWI算法、归一化差分水体指数NDWI算法和归一化差分植被指数NDVI算法；

基于所述水体指数和植被指数，判断是否存在水体；

使用水体类型识别模型对所述水体的水质进行识别。

2.根据权利要求1所述的黑臭水体遥感监测方法，其特征在于，所述基于所述水体指数和植被指数，判断是否存在水体包括：

在植被较为稀疏的地区，当MNDWI_t0≥a₁成立时，判断存在水体，其中MNDWI_t0是t₀时刻的MNDWI值，a₁是预设阈值；

或，在植被较为茂盛的地区，当MNDWI_ti＜a₁，NDVI_ti≥a₂，MNDWI_ti-j≥a₁并且NDVI_ti-j＜a₂成立时，判断存在水体，其中MNDWI_ti，MNDWI_ti-j，NDVI_ti和NDVI_ti-j分别是t_i、t_i-j时刻的MNDWI值和NDVI值，j<i，a₂是预设阈值。

3.根据权利要求1所述的黑臭水体遥感监测方法，其特征在于，所述使用水体类型识别模型对所述水体的水质进行识别包括：

对于非细小水体，建立包括黑臭水体子模型1、黑臭水体子模型2和黑臭水体变异系数的水体类型识别模型，其中黑臭水体子模型1的算法为：

其中R_g是绿波段反射率，R_r是红波段反射率；黑臭水体子模型2的算法为：

其中R_g是绿波段反射率，R_b是蓝波段反射率；黑臭水体变异系数的算法为：

其中

R_i为由R_g、R_r和R_b组成的集合中的元素。

4.根据权利要求3所述的黑臭水体遥感监测方法，其特征在于，当HCS₁≥a₃、HCS₂≥a₄并且HCC≥a₅成立时，判断所述水体是正常水体，其中a₃、a₄和a₅是预设阈值；

或，当a₆≤HCS₁＜a₃、a₇≤HCS₂＜a₄并且a₈≤HCC＜a₅成立时，判断所述水体是轻度黑臭水体，其中a₆、a₇和a₈是预设阈值；

或，当HCS₁＜a₉、HCS₂＜a₁₀并且HCC＜a₁₁成立时，判断所述水体是严重黑臭水体，其中a₉、a₁₀和a₁₁是预设阈值。

5.根据权利要求1所述的黑臭水体遥感监测方法，其特征在于，所述使用水体类型识别模型对所述水体的水质进行识别包括：

对于细小水体，建立包括高程变异系数的水体类型识别模型，其中高程变异系数的算法为：

其中H_i为预定规格的网格高程、

为研究区预定规格的网格高程均值，n为预定规格的矢量网格覆盖的数字高程模型栅格数量。

6.根据权利要求5所述的黑臭水体遥感监测方法，其特征在于，当L≤a₁₂，|S_i-S_j|≥a₁₃并且DXQF≤a₁₄成立时，判断所述水体是黑臭水体，其中L为细小水体的边界长度，a₁₂为预设阈值；S_i和S_j为距离最近的两个细小斑块的中心点，a₁₃为相邻两个斑块的最短距离，a₁₄为预设阈值。

7.一种黑臭水体遥感监测装置，其特征在于，包括：

数据获取模块，用于获取遥感影像数据；

水体提取模块，用于使用水体提取模型对所述遥感影像数据进行处理得到水体指数和植被指数，其中所述水体提取模型包括改进的归一化差分水体指数MNDWI算法、归一化差分水体指数NDWI算法和归一化差分植被指数NDVI算法；

水体判断模块，用于基于所述水体指数和植被指数，判断是否存在水体；

水质识别模块，用于使用水体类型识别模型对所述水体的水质进行识别。

8.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至6任一项所述黑臭水体遥感监测方法的步骤。

9.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1至6任一项所述黑臭水体遥感监测方法的步骤。

10.一种计算机程序产品，所计算机程序产品包括有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1至6任一项所述黑臭水体遥感监测方法的步骤。

Images