一种基于数学形态学的遥感图像提取连续河流骨架线方法

技术领域

本发明属于遥感图像处理领域，特别涉及一种基于数学形态学的遥感 图像提取连续河流骨架线方法。

背景技术

目前，水利行业对河流骨架线的提取主要依靠实地测绘，人工勾画， 工作量大，工作周期长。特别是随着影像分辨率的提高，用户对河流骨架 线的精度提出了更高的需求。从遥感影像中提取河流骨架线，能有效的降 低人工勾画的工作量，且更能保证骨架线位于河流中心，这也成为遥感技 术与水利业务结合的必然结果。但此方法提取出的骨架线受水体提取漏提 影响，会发生断裂，从而严重与自然语义不符，也对骨架线的分析、应用 造成了阻碍。

发生断裂是线提取面临的重大问题，目前连接断裂线的主流方法有：1) 最小点对法，根据断线的形态特点，连接距离最短，或方向差异性最小的 点，该方法的缺点是发生错误的概率较高；2)膨胀细化法，通过膨胀将断 裂部分连接，再进行细化，该方法的缺点是对断裂较大的图像会造成骨架 线严重变形。针对上述方法的不足，学者针对等高线、道路、边界等不同 对象继续进行了探索，并取得了一定成果。但由于河流所在图像通常范围 较大，干扰性断点很多，且断裂长短不一，河流走向复杂，目前相关研究 较少，上述的断线连接方法也难以适用。

发明内容

针对现有技术中河流骨架线提取断裂的问题，本发明提供了一种基于 数学形态学的遥感图像提取连续河流骨架线方法，该在获取初始断裂骨架 线基础上，对断裂处进行自适应连接，并对连接线进行去伪，从而获取连 续的河流骨架线。

为解决上述技术问题，本发明采用如下的技术方案：

一种基于数学形态学的遥感图像提取连续河流骨架线方法，其特征包 括以下步骤：

步骤1，基础河流骨架线提取，包括：

步骤1.1，对影像进行水体提取：计算影像的NDWI指数，其中，R_green为影像绿波段，R_NIR为影像的近红外波段；然后 对影像的NDWI指数进行直方图统计，采用谷值法计算分割阈值，提取水 体；

步骤1.2，对提取出的水体进行细化，提取初始骨架线：基于Rosenfeld细 化算法提

取水体骨架线，骨架像元值为1，背景像元值为0；

步骤2，断裂骨架线自适应连接，包括：

步骤2.1，统计骨架像元八邻域情况，骨架像元的值DN_i＝∑_j∈IDN_j，I 为像元i八邻域集合，从而将骨架像元分为值为端点、普通点、节点三类点， 所述三类点包括端点(DN_i＝1)、普通点(DN_i＝2)、节点(DN_i＝3)；

步骤2.2，根据图像断裂情况确定膨胀尺度，具体是：对所有端点计算 与其距离最短的端点之间的距离的1/2，并进行统计端点在各距离值的聚集 情况，形成曲线，由于河流骨架线的端点主要由断裂造成，用曲线表征骨 架线断裂情况，曲线陡降结束点的距离值为膨胀尺度Scale；

步骤2.3，对图像端点进行局部膨胀，膨胀尺度为Scale；

步骤2.4，对图像端点类型进行分类，对膨胀后影像进行腐蚀-膨胀操作， 腐蚀尺度为Scale+k，再次膨胀尺度为Scale+m，根据再次膨胀情况选择执 行：

选择步骤一：若端点的膨胀区域消失，则为断裂形式单一情况，直接 对首次膨胀所得区域进行Rosenfeld细化获取连接线；

选择步骤二：若膨胀区域未消失，则所膨胀区域包含的端点为河流骨 架有分支时的多断点聚集情况，提取膨胀区域边界上的骨架像元，河流对 于该区域而言，为一头进、一头出的形态，以膨胀区域边界上骨架点作为 起始骨架像元，向膨胀区域内进行跟踪，寻找并连接断点或节点，获得连 接线；

步骤3，连接线去伪：

步骤3.1，利用初始骨架将穿过河流内部的连接线切断，条件如下：

$\left(\begin{array}{c}{\mathrm{DN}}_i=2or\mathrm{><mn>3</mn>}\\ x_j=1,j\in I^{,}\\ \forall {\mathrm{DN}}_k\ne 1,k\in J^{,}\end{array}\right)$

式中，DN_i为骨架像元i的值，I’为i的八邻域余i本身的集合，j为I’中 的像元，如果为连接线像元，则值x_j＝1，J’为j的八邻域和j本身集合，即 骨架像元i的八邻域中，如果有连接线像元j，则将该连接线像元删除，但 由于有的连接线的一端距离端点附近的骨架线过近，因此还需加最后一个 限定条件，J’中的所有像元k都不为骨架线端点；

步骤3.2，删除伪连接线，首先对连接线像元进行骨架像元八邻域情况 的统计，再进行删除，具体是从连接线切断点开始，跟踪连接线，并将其 值赋为背景像元值，其中，

像元为连接线切断点的判断条件如下：

$\left(\begin{array}{c}{x}_i=1\\ \forall {\mathrm{DN}}_j=0,j\in I^{,}\end{array}\right)$

跟踪结束的判断条件为：

${\mathrm{DN}}_i=1or\forall x_j\ne 1,j\in I.$

在上述的一种基于数学形态学的遥感图像提取连续河流骨架线方法， 所述的背景像元值为0；初始骨架点值为1，骨架点经统计过后，端点值为 1，普通点值为2，节点值为3。

在上述的一种基于数学形态学的遥感图像提取连续河流骨架线方法， 所述的膨胀尺度统计以整十数为分界，求出的所有端点的最小距离的1/2 用进一法化为整十数，从而体现统计特性；膨胀尺度通常不超过50，尺度 为统计数陡降的结束点。

在上述的一种基于数学形态学的遥感图像提取连续河流骨架线方法， 所述的腐蚀尺度为Scale+k，再次膨胀尺度为Scale+m，k∈(0,5)， m∈(k,10)，以保证腐蚀时能将断裂形式单一的端点全部去除，膨胀时能 将所有河流骨架有分支时的多端点聚集断点包含。

本发明具有如下特点和有益效果：1、本发明采用了自适应的断线连接 方式连接骨架，利用最小点对法思想获取膨胀尺度，从端点聚集情况入手， 自适应连接不同情况断点，断裂形式单一时，用端点局部膨胀细化法获取 的连接线，河流骨架有分支造成多端点聚集时，通过分析河流一头入、一 头出的形态，连接断裂骨架。2、本发明引入了连接线去伪的机制，以往的 断线连接技术中，人们都着重考虑如何将断线连接起来，但是连接过程中， 往往会出现错连现象，特别是对于河流而言，不确定因素大，干扰众多， 必须引入去伪的机制，以保证连接的正确性，通过分析河流不能流入其本 身的特点，利用初始骨架线切断所有与其内部相交的连接线，并对被切断 的连接线进行消除，从而去除伪连接线。

附图说明

图1为膨胀尺度统计选择示例图。

图2a为断裂形式单一断点情况膨胀示意图。

图2b为断裂形式单一断点情况膨胀区域腐蚀示意图。

图2c为断裂形式单一断点情况膨胀区域细化示意图。

图3a为河流骨架有分支时的多端点聚集情况膨胀示意图。

图3b为河流骨架有分支时的多端点聚集情况膨胀区域腐蚀示意图。

图3c为河流骨架有分支时的多端点聚集情况再次膨胀示意图。

图3d为河流骨架有分支时的多端点聚集情况断裂骨架线跟踪连接示意 图。

图4a为普通端点膨胀示意图。

图4b为普通端点错误连接示意图。

图4c为错误连接线打断并跟踪删除示意图。

图4d为错误连接线去除结果示意图。

图5为本发明的方法流程示意图。

具体实施方式

下面将详细描述本发明的具体实施例。应当注意，这里描述的实施例只 用于举例说明，并不用于限制本发明。

1、提取初始河流骨架线。

计算影像的NDWI指数，其中，R_green为影像绿波 段，R_NIR为影像的近红外波段；

对影像的NDWI指数进行直方图统计，采用估值法计算分割阈值，提 取水体；

利用Rosenfeld细化算法提取水体骨架线，骨架像元值为1，背景像元值为 0；

2、自适应连接断裂骨架线。

统计骨架像元八邻域情况，骨架像元的值DN_i＝∑_j∈IDN_j，I为像元i 八邻域集合，从而将骨架像元分为值为端点(DN_i＝1)、普通点(DN_i＝2)、 节点(DN_i＝3)三类点；

根据图像断裂情况确定膨胀尺度，对所有端点计算与其距离最短的端 点之间的欧式距离，将该值的1/2用进一法保留为整十数，并进行统计，形 成曲线，如图1所示，数量陡降处结束点的距离值为膨胀尺度Scale，即图 一中的红点处，此图应取30；

对图像端点进行膨胀，膨胀尺度为30；

对图像端点类型进行分类，对膨胀后影像进行腐蚀后膨胀，腐蚀尺度 为30+3，再次膨胀尺度为30+5，以保证腐蚀时能将断裂形式单一的端点全 部去除，膨胀时能将所有河流骨架有分支时的多端点聚集断点包含。膨胀 区域消失的端点为断裂形式单一情况，直接对膨胀后区域进行细化，膨胀 区域相交的端点将被相连，获得连接线如图2所示；膨胀区域未消失的端 点为河流骨架有分支时的多端点聚集情况膨胀区域内的端点，提取膨胀区 域边界上的骨架像元，河流对于该区域而言，为一头进、一头出的形态， 通常，相对边界上的骨架像元所在的骨架线方向差异性更小，通过膨胀区 域边界位置关系，确定目标骨架像元，从骨架像元开始向膨胀区域内进行 跟踪，寻找并连接断点，获得连接线如图3所示；

3、去除伪连接线。

利用初始骨架将穿过河流内部的连接线切断，条件如下：

$\left(\begin{array}{c}{\mathrm{DN}}_i=2or\mathrm{><mn>3</mn>}\\ x_j=1,j\in I^{,}\\ \forall {\mathrm{DN}}_k\ne 1,k\in J^{,}\end{array}\right)$

式中，DN_i为骨架像元i的值，I’为i的八邻域余i本身的集合，j为I’中 的像元，如果为连接线像元，则值x_j＝1，J’为j的八邻域和j本身集合，即 骨架像元i的八邻域中，如果有连接线像元j，则将该连接线像元删除，但 由于有的连接线的一端距离端点附近的骨架线过近，因此还需加最后一个 限定条件，J’中的所有像元k都不为骨架线端点。

删除伪连接线，即从连接线切断点开始，扫描当前点的8领域，如果存 在连接线像元，将当前点赋值为背景像元，并将下一个连接线像元作为当 前像元继续扫描，像元为连接线切断点的判断条件如下：

$\left(\begin{array}{c}{x}_i=1\\ \forall D{N}_j=0,j\in I^{,}\end{array}\right)$

跟踪结束的判断条件为：

${\mathrm{DN}}_i=1or\forall x_j\ne 1,j\in I.$

伪连接线的形成与去除如图4所示。