一种遥感影像内河洲滩智能化提取方法与流程

本发明属于遥感技术领域，涉及一种遥感影像内河洲滩智能化提取方法。

背景技术：

内河洲滩是维系河流系统生态健康的关键区域，随着社会经济的不断发展，对于内河洲滩资源的有效利用尤其重要，一些内河洲滩上有人居住，在河流汛期时，由于水位上涨，部分内河洲滩的区域可能会被洪水淹没，这对内河洲滩居民的生命财产安全会造成重大威胁，随着遥感技术的发展，从遥感影像提取内河洲滩成为可能，并将成为内河洲滩监测的重要技术手段，遥感影像内河洲滩提取，对于内河洲滩资源的合理利用、防洪安全、生态环境建设与保护具有重要意义。

技术实现要素：

本发明针对内河洲滩提取研究领域存在的问题，提供一种遥感影像内河洲滩智能化提取方法。

本发明的技术方案为一种遥感影像内河洲滩智能化提取方法，基于遥感水体指数与神经网络，包含以下步骤：

步骤1，在得到多光谱遥感影像后，利用贝叶斯正则化反馈神经网络对影像进行水陆识别，获得水陆二值化图像；

步骤2，在得到由步骤1处理后的水陆二值化图像后，对二值化图像进行非目标去除，所述非目标包括池塘、水田和部分船舶；

步骤3，进行内河洲滩提取，实现方式如下，

(3.1)对于洲滩大面积部分，采用区域生长方法进行内河洲滩提取，首先人工随机选取生长点，再执行区域生长算法进行八邻域扩张，从而得到目标内河洲滩；

(3.2)对于数量多、面积小的洲滩部分区域，采用手动框型选取方式获取框内内河洲滩像素；

步骤4，内河洲滩提取结果优化与渲染，实现方式如下,

(4.1)内河洲滩提取结果优化，设置为船舶大小阈值为t1，设置岛屿陆地大小阈值为t2，遍历图像中的像素，使用船舶大小阈值t1获取可能船舶像素，并将该像素作为种子点，然后使用区域生长算法来获取可能船舶对象集，然后进行判断，若船舶对象集像素数量<t2，则判定为船舶对象，予以去除，完成内河洲滩提取结果的优化；

(4.2)内河洲滩提取结果渲染，将影像中的水体部分、内河洲滩部分、岸上陆地及其他地物分别设置为不同颜色，以此渲染并显示内河洲滩提取结果。

进一步的，步骤1中利用贝叶斯正则化反馈神经网络对影像进行水陆识别的具体实现方式如下，

以遥感影像的蓝色波段、红色波段、mndwi水体指数作为输入特征向量；

在原图像上选取训练样本与测试样本，然后进行贝叶斯正则化反馈神经网络的训练，所述贝叶斯正则化反馈神经网络包括3层，其中输入层与输入向量特征值数量相等有3个神经元，隐含层包含10个神经元，输出层与分类目标数目相等由2个神经元构成，分类目标数目为水体和非水体两类；

最后利用训练好的贝叶斯正则化反馈神经网络对图像进行预测，获得水陆二值化图像。

进一步的，所述贝叶斯正则化反馈神经网络采用贝叶斯正则化算法作为训练函数，从输入层到隐含层的激活函数设置为tansig函数，从隐含层到输出层的激活函数设置为logsig函数。

进一步的，步骤2的具体实现方式如下，

(2.1)首先人工随机选取生长点，再执行区域生长算法进行八邻域扩张，从而得到大面积目标水体；另外采用手动框型选取方式作为补充，实现小面积目标水体的获取，手动框型选取为用户指定矩形框并获取框内水体；

(2.2)通过用户手动选取矩形框，并去除框内船舶点像素。

进一步的，还包括步骤5，进行内河洲滩提取结果的精度评价，实现方式如下，

首先获得内河洲滩岸线，再随机生成p个样本点，并目视解译这p个样本点，分为水体点及非水体点两类，而后建立混淆矩阵，以目视解译结果作为参考值，以内河洲滩提取效果作为预测值，进行基于混淆矩阵精度评价指标的内河洲滩提取结果精度评价，最终获得内河洲滩提取精度评价结果。

进一步的，p的取值为200，样本点的选取原则为：在内河洲滩岸线八邻域范围内随机选取100个样本点，并在整幅影像范围内随机选取100个样本点。

本发明利用内河洲滩及水体的空间形态特性及光谱特性，通过联合遥感水体指数与反馈神经网络并经过多个步骤对内河洲滩进行了有效提取，并进行了基于混淆矩阵的内河洲滩提取结果精度评价，验证了方法的有效性。

附图说明

图1是本发明实施例的一种遥感影像内河洲滩智能化提取方法流程图。

图2是本发明实施例的反馈神经网络结构示意图。

图3是本发明实施例的神经网络水陆识别步骤流程图。

图4是本发明实施例的非目标去除步骤流程图。

图5是本发明实施例的内河洲滩提取步骤流程图。

图6是本发明实施例的内河洲滩提取结果优化与渲染步骤流程图。

图7是本发明实施例的内河洲滩提取结果精度评价步骤流程图。

具体实施方式

本发明主要基于遥感技术，考虑到遥感影像内河洲滩形态及光谱特性，提出的一种遥感影像内河洲滩智能化提取实验方法及系统。本方法充分考虑了内河洲滩提取过程中的各种问题，通过mndwi水体指数与反馈神经网络方法以及相关影像处理技术进行内河洲滩的提取。通过本发明可以获得更为精确完整的内河洲滩提取结果。

本发明提供的方法能利用计算机软件技术进行编程实现。参见图1，实施例以长江内河洲滩提取为例对本发明流程进行具体说明，如下：

步骤1，在获取多光谱遥感影像后，对影像进行基于遥感水体指数和神经网络的水陆识别，将遥感水体指数与神经网络进行联合，设计了相应的贝叶斯正则化反馈神经网络结构，与标准反馈神经网络相比，该网络具有更好的泛化和识别能力，与深度网络相比，该网络训练时间短。其识别方法为，首先计算多光谱遥感影像的mndwi水体指数，再构建针对水陆识别的贝叶斯正则化反馈神经网络结构，此神经网络以多光谱遥感影像的蓝色波段、红色波段、mndwi水体指数作为输入特征向量，在原影像上选取200个训练样本，基于levenberg-marquardt优化进行神经网络训练，最后利用训练好的神经网络对影像进行水陆识别，获得水陆二值化影像；

实施例的具体过程说明如下：

在读取多光谱遥感影像后，利用影像的绿色波段和中红外波段进行mndwi水体指数的计算，mndwi利用建筑物与土壤及其他地物与水体之间的波谱差异，采用创造波段差及比值的方式对水体信息进行强调，经实验验证具有较高的水体提取精度，能够有效的对水体目标进行提取，因此选用mndwi水体指数作为神经网络的输入特征向量分量，mndwi水体指数计算公式为：mndwi＝(green-mir)/(green+mir)，此外由于遥感影像蓝色波段及红色波段对水体比较敏感，能够有效的将水体与其他地物区分开来，因此本发明选用影像中的蓝色波段、红色波段及mndwi水体指数作为神经网络的输入特征向量，其表达式为：x＝(blue,red,mndwi)^t，本文中神经网络的分类目标是将影像中的地物分为水体及非水体两类，水体目标矩阵表达形式为：ywater＝[10]^t，非水体目标矩阵表达形式为：ynon-water＝[01]^t，综上理由，本发明中反馈神经网络结构被分为3层，其中输入层与输入向量特征值数量相等有3个神经元，隐含层设置为一层共10个神经元时具有较好的分类效果，因此隐含层包含一层共10个神经元，输出层与分类目标数目相等由2个神经元构成，本发明反馈神经网络结构如图2所示。

在训练样本选取方面，以长江内河洲滩遥感影像为例，首先在影像上随机生成200个点，再使用envi5.1软件对这200个点进行目视解译，目视解译结果分为水体与非水体两类，然后将这些样本点作为训练样本进行反馈神经网络的训练，由于贝叶斯正则化算法具有泛化能力强，分类精度高的优势，本发明神经网络采用贝叶斯正则化算法(trainbr)作为训练函数，从输入层到隐含层的激活函数设置为tansig函数，从隐含层到输出层的激活函数设置为logsig函数，训练的最大次数为1000次，训练全局最小误差被设置为0.005，显示频率为1000，动量因子设置为0.85，需要注意的是训练样本要进行数据归一化操作，训练结束后，再利用训练好的反馈神经网络对整幅影像进行预测，得到未经过非目标去除步骤的二值化影像。

如图3所示，具体流程为，计算影像mndwi指数，再选取训练样本和测试样本，之后利用训练样本与测试样本进行反馈神经网络训练，然后再进行神经网络水陆识别最终获得水陆二值化影像。

步骤2，在得到由步骤1处理后的水陆二值化影像后，对二值化影像进行非目标去除，非目标去除主要包括两部分内容，实现方式如下，

一是获取目标水体，排除非目标水体干扰，在目标水体周围可能存在着一些小型水体，如池塘、水田等，而且在某些内河洲滩内部也存在着需要保留的目标水体，这些水体主要特点是数量多、面积小，针对这个问题，本发明结合自动与手动非目标去除方法对目标水体进行提取，对于大面积水体，本发明采用区域生长方法进行水体提取，首先人工随机选取生长点，再执行区域生长算法进行八邻域扩张，从而得到目标水体；对于内河洲滩内部的具有数量多、面积小特征的目标水体，采用手动框型选取方式获取目标水体，水体框型选取为用户指定矩形框并获取框内水体。

二是部分船舶去除，为避免对内河洲滩提取造成干扰，需要对影像中的部分船舶进行去除，考虑到不同分辨率遥感影像中船舶尺寸差别较大，本发明采用部分船舶手动框型去除方法进行非目标去除，该方法主要思想是用户手动选取矩形框，并去除框内船舶点像素。

实施例的具体过程说明如下：

以长江内河洲滩提取为例，在第一步获取目标水体的过程中，对于长江大面积水域由人工随机选取生长点并执行区域生长算法，从而获取长江主干水域，其中人工随机选取生长点以文件形式存储并读取，区域生长算法具体思想是由水体生长点出发获取其空间邻域内其他水体点，并将其他水体点也作为生长点，并重复这个过程，直至获得完整的水体连通区域，对于长江内河洲滩中由于水位上涨而产生的大量小型水体，采用框型选取方式获取目标水体，用户指定一矩形区域并获取该矩形区域内的所有水体点，其中矩形区域的选择需要用户交互完成。

由于长江内河洲滩周围可能有一些船舶，这些船舶与洲滩部分小型区域大小相似，且位置较近，为避免对后续内河洲滩提取造成影响，通过目视解译方式进行部分船舶手动框型去除，用户指定一矩形区域并去除该矩形区域内的所有船舶点，其中矩形区域的选择需要用户交互完成。

如图4所示，具体流程为，首先加载水体生长种子点，即获取长江水体主干部分，再进行水体手动框型选取，即获取小型目标水体，然后进行部分船舶手动框型去除，最终获得非目标去除后的二值化影像。

步骤3，在获取步骤2经非目标去除后的水陆二值化影像后，内河洲滩提取步骤得以展开，实现方式如下，

内河洲滩提取分为两步进行，这是由内河洲滩的特点所决定的，对于某一内河洲滩来说，其原本可能为一整体，然而若发生洪水后，该洲滩的地理形态可能发生较大变化，影像上可能呈现出零散稀疏，本发明结合自动与手动方法对目标内河洲滩进行提取，对于洲滩大面积部分，本发明采用区域生长方法进行内河洲滩提取，首先人工随机选取生长点，再执行区域生长算法进行八邻域扩张，从而得到目标内河洲滩；对于数量多、面积小特征的洲滩部分区域，采用手动框型选取方式获取这一部分区域，内河洲滩框型选取主要思想是用户手动选取矩形框，并获取框内内河洲滩像素。

实施例的具体过程说明如下：

以长江内河洲滩提取为例，第一步采用区域生长算法获取长江内河洲滩主体部分，具体做法是人工随机选取生长点，再将这些生长点作为区域生长算法的种子点进行长江内河洲滩主体部分的生长获取，其中这些生长点以文件形式存储并读取；第二步对于数量大、面积小的长江内河洲滩部分区域，为获取更为准确的内河洲滩提取效果，通过目视解译方式进行内河洲滩手动框型选取，用户指定一矩形区域并获取该矩形区域内的所有内河洲滩点，其中矩形区域的选择需要用户交互完成，通过该步骤即可获得完整的长江内河洲滩。

如图5所示，具体流程为，首先加载内河洲滩生长种子点，即获取长江内河洲滩主体部分，再进行内河洲滩手动框型选取，即获取长江内河洲滩部分区域，最终获得长江内河洲滩提取结果。

步骤4，为合理显示内河洲滩提取结果，本发明还需进行内河洲滩提取结果优化与渲染步骤，该步骤主要包括两部分内容，实现方式如下,

一是内河洲滩提取结果优化，在步骤2中，可能对步骤3内河洲滩提取造成影响的部分船舶被去除，然而影像中仍然存在着大量距离内河洲滩较远的船舶未被去除，因此为合理显示内河洲滩提取结果，需要进行船舶自动去除，该方法利用船舶和目标水体之间的空间关系，即船舶被水体所包围，船舶大小阈值设置为t1，岛屿陆地大小的阈值设置为t2。遍历影像中的像素，并使用船舶大小阈值t1获取可能的船舶像素，并将该像素作为种子点，然后使用区域生长算法来获取完整对象，如果获得的完整对象的大小小于t2，它被判断为船舶并被去除，在阈值选取方面，考虑到不同分辨率遥感影像中船舶尺寸差别较大，阈值选取需由操作人员设定。

二是内河洲滩提取结果渲染，利用之前步骤获得的相关结果进行内河洲滩提取的结果渲染，其中影像中的水体部分设置为蓝色、内河洲滩部分被设置为黄色、岸上陆地及其他地物被设置为黑色，以此渲染并显示内河洲滩提取结果。

实施例的具体过程说明如下：

以长江内河洲滩提取为例，首先遍历长江内河洲滩影像，设置船舶大小阈值为t1，若影像中某一像素的左右两方t1像素范围内都存在水体像素，则该像素可判定为可能船舶像素，之后对可能船舶像素集进行区域生长获得可能船舶对象集，这个过程中若可能船舶对象集中有待进行区域生长的像素在之前的区域生长中已经进行过遍历则不会进行重复生长，设置阈值t2，若可能船舶对象包含像素数量小于t2，则判断为船舶进行去除，考虑到实际应用场景的复杂性，阈值t1及t2由实际操作人员设置，通过该步骤，影像中的船舶将会被自动去除，内河洲滩提取结果的优化得以完成。

在内河洲滩提取结果的渲染中，影像中的水体部分设置为蓝色、内河洲滩部分被设置为黄色、岸上陆地及其他地物被设置为黑色，以此渲染并显示长江内河洲滩提取结果。

如图6所示，具体流程为，首先遍历影像像素，再设置阈值t1获取可能船舶像素，执行区域生长获取可能船舶对象集，然后进行判断，若船舶对象集像素数量<t2，则判定为船舶对象，予以去除，完成内河洲滩提取结果的优化，最后进行内河洲滩提取结果的渲染。

步骤5，获取内河洲滩提取结果后，最后进行内河洲滩提取结果的精度评价，实现方式如下，

首先进行内河洲滩的岸线提取，并在岸线八邻域范围以及整幅影像范围内随机生成样本点，并进行目视解译，从而建立混淆矩阵，并基于混淆矩阵相关评价指标进行内河洲滩提取结果的精度评价。

实施例的具体过程说明如下：

以长江内河洲滩提取为例，在获取步骤3获得的长江内河洲滩提取结果后，进行岸线提取，其主要思想是：岸线属于非水体像素，在岸线像素的八个邻域中必须有一个水体像素，根据岸线的上述空间邻域关系遍历影像中的每个像素，即可获得内河洲滩岸线。然后再随机生成200个样本点，需要注意的是样本点的选取策略是在内河洲滩岸线八邻域范围内随机选取100个样本点，并在整幅影像范围内随机选取100个样本点，这样充分利用了内河洲滩的线状与面状特征，并目视解译这200个样本点，分为水体点及非水体点两类，而后建立混淆矩阵，以目视解译结果作为参考值，以内河洲滩提取效果作为预测值，进行基于混淆矩阵精度评价指标的内河洲滩提取结果精度评价，最终获得内河洲滩提取精度评价结果。

如图7所示，具体流程为，进行内河洲滩岸线提取，再选取样本点，建立混淆矩阵，最后得到内河洲滩提取精度评价结果。

本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。