一种浅水湖泊面积提取的方法

技术领域

本发明属于遥感技术与信息科学技术领域，特别是涉及一种浅水湖泊面积提取的方法。

背景技术

遥感技术是大范围、实时动态监测湖泊面积的主要技术之一。基于遥感影像提取湖泊面积的方法主要有指数(水体指数等)阈值、决策树分类、深度学习等。指数阈值方法主要根据遥感影像波段特点，利用中波红外、短波红外、红、绿等波段的算术运算构建水体指数图像，进而设定阈值来提取湖泊水域面积，这类方法提取精度跟水体指数计算效果、阈值设定好坏等直接相关，并且容易受到影像噪声、不纯净水体、水生植物等的影响。决策树分类主要是对不同波段图像进行阈值分类，并对所有分类结果进行决策级融合，实现湖泊面积的提取，这类方法提取精度与每个阈值分类结果、决策级融合规则等直接相关，而且该方法通常不具有泛化能力，一个湖泊面积的提取模型基本不能用于其他湖泊，工程化能力较差。深度学习方法通过大量先验湖泊样本，训练语义分割深度学习网络参数，并对湖泊区域影像进行处理，从而实现湖泊面积的提取，这类方法虽然具有较好的泛化能力，但是受到样本数量质量、类型全面性、分布均匀性等因素限制，需要消耗大量人力进行样本制作和大量时间进行样本迭代，工程化成本较高。

此外，由于很多浅水湖泊中常常伴生有大量水生植物，成片的、高于水面的水生植物往往造成遥感影像上该区域地表覆被表现为植被类型，以上三种方法基本无法成功提取水生植物覆盖区域的面积，极大影响了湖泊面积提取的精度，显著降低了湖泊面积数据的准确性。

发明内容

本发明的目的是提供一种浅水湖泊面积提取的方法，以解决上述现有技术存在的问题。

为实现上述目的，本发明提供了一种浅水湖泊面积提取的方法，包括以下步骤：

获取水体和水生植物的样本影像和标签数据，并对所述标签数据进行腐蚀和编码处理，获得正负波段图像；

基于所述正负波段图像和样本影像对交互式语义分割模型进行训练，获得训练后的交互式语义分割模型；

基于遥感影像获取待提取湖泊的水体指数和植被指数；

对所述水体指数和植被指数分别进行腐蚀处理，基于处理后的数据构建待提取湖泊的正负波段图像；

基于所述训练后的交互式语义分割模型对待提取湖泊的遥感影像和正负波段图像进行处理，获得湖泊面积提取结果。

可选地，获取水体和水生植物的样本影像和标签数据的过程包括：获取水体和水生植物相关监测数据以及对应的遥感影像，基于所述水体和水生植物相关监测数据裁切同空间位置区域的遥感影像，获得水体及水生植物图斑区域影像块；基于交互式语义分割模型的预训练权重对所述水体及水生植物图斑区域影像块进行交互式分割处理，获得水体和水生植物的样本影像和标签数据。

可选地，对所述标签数据进行腐蚀处理的过程包括：对所述标签数据基于网格标签进行正负像元的全向递归腐蚀运算，获得标签数据的正、负点。

可选地，对所述标签数据进行编码处理的过程包括：基于加权距离编码算法将所述标签数据的正、负点转换为正负波段图像。

可选地，获取待提取湖泊的水体指数和植被指数的过程包括：获取待提取湖泊水面面积最大范围以及对应的遥感影像，基于所述待提取湖泊水面面积最大范围裁切同空间位置区域的遥感影像，获得待提取湖泊图斑区域影像块；基于所述待提取湖泊图斑区域影像块获取所述待提取湖泊影像的水体指数和植被指数。

可选地，对所述水体指数和植被指数分别进行腐蚀处理的过程包括：对待提取湖泊内的水体指数进行全向腐蚀运算，提取待提取湖泊内的水体种子点；对待提取湖泊内的植被指数进行定向腐蚀运算，定向方向为最近水体种子点的同方向，提取待提取湖泊内的植被种子点；对待提取湖泊外的植被指数进行定向形态学运算，定向方向为最近水体种子点与最近植被种子点最小值的反方向，提取待提取湖泊外的植被种子点。

可选地，构建待提取湖泊的正负波段图像的过程包括：将待提取湖泊图斑区域影像块等间隔分为3*3区域；在每个区域的待提取湖泊内选择一个离待提取湖泊中心点相对近的水体种子点和植被种子点作为正点；在每个区域的待提取湖泊外选择一个离待提取湖泊中心点相对近的植被种子点作为负点；基于所述正点和负点构建待提取湖泊的正负波段图像。

可选地，获得湖泊面积提取结果的过程包括：基于训练后的交互式语义分割模型对待提取湖泊的正负波段图像以及待提取湖泊影像的合并图像进行交互式语义分割处理；对处理结果进行二值化、矢量化，获得湖泊面积提取初步结果；采用待提取湖泊水面面积最大范围对湖泊面积提取初步结果进行矢量裁切，得到湖泊面积提取结果。

本发明的技术效果为：

本发明生成的浅水湖泊面积，不仅无需人工设定阈值，而且无需遥感解译专业知识，通过智能设定交互式深度学习输入数据的正负极点波段，对于包含水生植物的浅水湖泊的提取精度更高，可实现全自动、高精度的湖泊面积提取。

附图说明

构成本申请的一部分的附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：

图1为本发明实施例中的浅水湖泊面积提取的方法流程示意图；

图2为本发明实施例中的待提取湖泊图斑区域影像块示意图；

图3为本发明实施例中的浅水湖泊面积提取效果；其中，(a)为湖泊区域遥感影像(R：4，G：3，B：2)，(b)为水体指数图像，(c)为植被指数图像；(d)为待提取湖泊内植被种子点图像，(e)为待提取湖泊湖泊外植被种子点图像，(f)为待提取湖泊内水体种子点图像，(g)为待提取湖泊的正点和负点图像，(h)为浅水湖泊面积提取结果图像。

具体实施方式

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。

需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。

实施例一

如图1所示，本实施例中提供了一种浅水湖泊面积提取的方法，其中，该湖泊面积提取的方法包括：利用自然资源调查监测数据，基于交互式语义分割网络和预训练权重参数，制作水体和水生植物的样本影像与标签数据；利用全向腐蚀和加权距离编码由标签数据构建正、负点及正负波段图像，与样本影像一并用于交互式语义分割模型参数更新，获得交互式语义分割模型；依据自然资源调查监测数据获得待提取浅水湖泊区域最大范围；利用Landsat7/8等包含中红外、红、绿波段的遥感影像，计算待提取浅水湖泊区域影像的水体指数和植被指数，分别经全向和定向腐蚀运算提取浅水湖泊区域内、区域外的水体和植被种子点；以湖泊区域内、外的水体种子点和植被种子点为正、负点，构建正、负波段图像；利用交互式语义分割模型，对提取结果错误的像元进行交互式的语义分割修正，获得湖泊面积提取结果。

具体的实施方法包括以下步骤：

收集第三次全国国土调查数据、地理国情监测数据、海南岛陆域水体遥感提取数据集、全国地表类型遥感影像样本数据集等自然资源调查监测数据，处理形成水体相关监测数据。

自然资源调查监测数据可表示为：

R＝{P

其中，R为自然资源调查监测数据集，P

针对每类调查监测数据P

P

其中，P

针对全部调查监测数据集R，对P

R’＝{UNION

其中，R’为水体及水生植物相关监测数据，P

利用水体及水生植物相关监测数据R’裁切同空间位置区域的光学遥感影像，基于轻量级交互式语义分割网络预训练权重对裁切影像块进行交互式分割处理，制作浅水水面的样本影像与标签数据。

作为具体的实施例，如图2所示，以水体及水生植物相关监测数据R’的每个图斑最小外截矩形向外缓冲，用缓冲后矩形裁切同空间位置区域的光学遥感影像，获得水体及水生植物图斑区域影像块。按照图斑最小外截矩形长宽的30％等比例向外部四个方向分别缓冲，缓冲后矩形长宽均为原影像的160％，考虑到语义分割模型需要对影像进行多次池化以获取该影像大步长的上下文语义信息，经3-4次池化后影像大小为原影像尺寸的1/5-1/10，最外层像元能包含水体外的背景信息，以利于水体及水生植物的正负样本的平衡。

基于HRNet18s+OCR轻量级交互式语义分割网络模型，参考遥感影像，采用目视解译的方式，在遥感影像上勾绘包含水生植物的水面面积图斑，获得水体的样本影像和标签数据。按照不大于50％的比例，对样本数据进行平移、旋转、翻转、缩放、增加噪音、亮度拉伸、色彩变换、拼接、填充、叠置等增广处理，经抽样和欠抽样，形成类别分布均衡的水体样本数据。

sample＝{Sample

其中，sample为水体样本数据集，Sample

对标签数据进行网格控制的腐蚀运算，构建正、负点波段图像，基于与样本影像一并输入HRNet18s+OCR交互语义分割网络进行训练，获得交互式语义分割模型获得交互式语义分割网络模型。

作为具体的实施例，对标签数据按网格进行正负像元的全向递归腐蚀运算，获得标签数据的正、负点。

Point_positive＝{Q

Point_nagetive＝{Q

B

其中，Point_positive为正点像元，Point_nagetive为负点像元，Sample

对标签数据的正、负点，采用加权距离编码算法将正负点转换为正、负波段图像。对全部正点、全部负点分别进行编码。

D

D

Point_mid＝(Σx

Gray

其中，D

利用HRNet18s+OCR交互式语义分割预训练参数，对水体样本数据、标签数据和正负点波段图像进行交互式语义分割模型更新训练。

由自然资源调查监测数据获得浅水湖泊水面面积最大范围，利用HRNet18s+OCR交互式语义分割模型对Landsat 7/8等包含中红外、红、绿波段的遥感影像进行处理，提取浅水湖泊水域面积。

自然资源调查监测数据获得浅水湖泊水面面积最大范围；

P

对Landsat 7/8等包含中红外、红、绿波段的遥感影像，同步骤2的方式获得待提取浅水湖泊区域影像，并计算待提取影像的水体指数和植被指数；

根据浅水湖泊水面面积最大范围，对湖泊区域内的水体指数进行全向腐蚀运算，提取待提取影像的水体种子点；

B＝B

其中，NDWI

根据浅水湖泊水面面积最大范围，对湖泊区域内的植被指数进行定向腐蚀运算，定向方向为最近水体种子点的同方向，提取湖泊区域内的植被种子点；

B＝{B

其中，NDVI

对湖泊区域外的植被指数进行定向形态学运算，定向方向为最近水体种子点与最近植被种子点最小值的反方向，提取湖泊区域外的植被种子点；

B＝{Bi＝1 if i in最近水体/植被种子点最小值的反向,Bi＝0else}

其中，NDVIout为湖泊区域外的植被指数图像，B为3*3卷积核。

湖泊面积提取：

将待提取区域等间隔分为3*3区域；对湖泊区域内的水体种子点和植被种子点，每个区域内分别选择一个离湖泊区域中心点相对近的水体种子点和植被种子点作为正点；对湖泊区域外的植被种子点，每个区域内分别选择一个离湖泊区域中心点相对近的植被种子点作为负点。

以正、负点构建正点距离波段和负点距离波段图像；利用训练好的交互式语义分割网络模型，对正点波段图像、负点波段图像、原始影像等合并图像进行交互式语义分割处理。

对深度学习结果进行二值化、矢量化，获得湖泊面积提取初步结果；采用湖泊区域范围对湖泊面积提取初步结果进行矢量裁切，得到湖泊面积提取结果。如图3所示，为本发明实施例中的浅水湖泊面积提取效果；其中，(a)为湖泊区域遥感影像(R：4，G：3，B：2)，(b)为水体指数图像，(c)为植被指数图像；(d)为待提取湖泊内植被种子点图像，(e)为待提取湖泊湖泊外植被种子点图像，(f)为待提取湖泊内水体种子点图像，(g)为待提取湖泊的正点和负点图像，(h)为浅水湖泊面积提取结果图像。

以上所述，仅为本申请较佳的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应该以权利要求的保护范围为准。