技术领域

本发明涉及水体提取领域，尤其涉及一种水体快速提取方法。

背景技术

水体指数法利用遥感影像多波段信息构建区分水体和非水体的波段运算函数，以达到提取水体的目的。

但采用的主流水体指数法存在以下问题：①高分辨率卫星影像包括的地物信息丰富，传统提取方法极易受暗色建筑、建筑阴影、裸土等信息的干扰，导致细小的河道、小鱼塘、小坑塘难以提取；②高分辨率卫星影像只有红波段(R)、绿波段(G)、蓝波段(B)、近红波段(NIR)，但现有大多水体指数计算都必须采用短波红外(SWIR)。

发明内容

针对现有技术中存在的问题，本发明提供一种水体快速提取方法，应用于遥感影像，所述遥感影像的波段范围包括第一波段、第二波段和第三波段；

所述水体快速提取方法包括：

步骤S1，根据所述第二波段和所述第三波段，获取归一化水体指数并生成第一影像；

步骤S2，根据所述第一波段和所述第三波段，获取归一化蓝色水体指数并生成第二影像；

步骤S3，选取多个水体特征点和非水体特征点，并根据所述水体特征点和所述非水体特征点分别于所述第一影像上的像元灰度值和所述第二影像上的像元灰度值，构建最佳分类线；

步骤S4，根据所述最佳分类线获取融合垂直水体指数；

步骤S5，预设一阈值，并根据所述阈值和所述融合垂直水体指数提取水体。

优选的，所述步骤S1，采用下述公式表示所述归一化水体指数：

其中，

NDWI用于表示所述归一化水体指数；

优选的，所述步骤S2，采用下述公式表示所述归一化蓝色水体指数：

其中，

NBWI用于表示所述归一化蓝色水体指数；

优选的，所述步骤S3包括：

步骤S31，选取多个所述水体特征点和多个所述非水体特征点；

步骤S32，根据所述水体特征点于所述第一影像上的像元灰度值以及所述水体特征点于所述第二影像上的像元灰度值进行线性拟合，获取水体线斜距和水体线截距，构建水体线；

步骤S33，根据所述非水体特征点于所述第一影像上的像元灰度值以及所述非水体特征点于所述第二影像上的像元灰度值进行线性拟合，获取非水体线斜距和非水体线截距，构建非水体线；

步骤S34，根据所述水体线和所述非水体线获取最佳分类线斜距，并根据所述最佳分类线斜距、所述水体线和所述非水体线的交点获取最佳分类线截距，根据所述最佳分类线斜距和所述最佳分类线截距构建所述最佳分类线。

优选的，所述水体线斜距和所述水体线截距采用下述公式表示：

其中，

a

b

x

y

优选的，所述非水体线斜距和所述非水体线截距采用下述公式表示：

其中，

a

b

x

y

优选的，所述最佳分类线斜距采用下述公式表示：

其中，

a

a

a

优选的，所述融合垂直水体指数采用下述公式表示：

其中，

IPWI用于表示所述融合垂直水体指数；

a

NBWI用于表示所述归一化蓝色水体指数；

NDWI用于表示所述归一化水体指数。

上述技术方案具有如下优点或有益效果：本发明提供一种水体快速提取方法，通过获取归一化蓝色水体指数能够增强建筑与水体之间的区别，通过获取融合垂直水体指数能够进一步增强水体与非水体之间的像元灰度值的区别，以快速识别并提取细小的水体，提高水体提取准确度，加快水体提取速度。

附图说明

图1为本发明的一种较优实施方式中，一种水体快速提取方法的流程图；

图2为本发明的一种较优实施方式中，步骤S3的流程图；

图3为本发明的一种较优实施方式中，最佳分类线的示意图；

图4为本发明的一种较优实施方式中，最佳分类线提取水体的示意图；

图5为本发明的一具体实施方式中，一种高景图像；

图6为本发明的一具体实施方式中，针对高景图像的水体提取示意图；

图7为本发明的另一具体实施方式中，一种最宽河流为5m的高景图像；

图8为本发明的另一具体实施方式中，针对最宽河流为5m的高景图像的水体提取示意图；

图9为本发明的另一具体实施方式中，具有鱼塘、坑塘、细小河流的高景图像；

图10为本发明的另一具体实施方式中，针对具有鱼塘、坑塘、细小河流的高景图像的水体提取示意图。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细说明。本发明并不限定于该实施方式，只要符合本发明的主旨，则其他实施方式也可以属于本发明的范畴。

本发明的主旨是解决现有技术中细小河道、小鱼塘、小坑塘难以提取的问题，通过本发明提供的提取方法，构建一种融合垂直水体指数，实现了细小水体的高效提取。

以下提供的具体技术手段均为实现本发明主旨的举例说明，可以理解的是，在不冲突的情况下，以下所举的实施方式，及实施方式中的技术特征均可相互组合。并且，不应当以用于说明本发明可行性的实施方式来限定本发明的保护范围。

下面结合附图和具体实施例对本发明作进一步说明，但不作为本发明的限定。

一种水体快速提取方法，应用于遥感影像，所述遥感影像的波段范围包括第一波段、第二波段和第三波段；

如图1所示，水体快速提取方法包括：

步骤S1，根据第二波段和第三波段，获取归一化水体指数并生成第一影像；

步骤S2，根据第一波段和第三波段，获取归一化蓝色水体指数并生成第二影像；

步骤S3，选取多个水体特征点和非水体特征点，并根据水体特征点和非水体特征点分别于第一影像上的像元灰度值和第二影像上的像元灰度值，构建最佳分类线；

步骤S4，根据最佳分类线获取融合垂直水体指数；

步骤S5，预设一阈值，并根据阈值和融合垂直水体指数提取水体。

具体地，此处的第一波段为波长范围在520mm与590mm之间的蓝波段，第二波段为波长范围在450mm与520mm之间的绿波段，第三波段为波长范围在770mm与890mm之间的近红波段。

本发明为避免暗色建筑、建筑阴影、裸土等信息干扰到细小河道、小池塘等的提取结果，因此通过根据绿波段和近红波段获取归一化水体指数(Normalized DifferenceWater Index，简称NDWI)，以及根据蓝波段和近红波段获取归一化蓝色水体指数(Normalized Blue Water Index，简称NBWI)，并选取多个水体特征点和非水体特征点以构建最佳分类线，根据最佳分类线获取融合垂直水体指数(Integrated PerpendicularWater Index，简称IPWI)，从而提取出细小水体。

本发明优选实施方式中，步骤S1，采用下述公式表示归一化水体指数NDWI：

其中，

NDWI用于表示归一化水体指数；

具体地，步骤S1中可通过波长范围在520nm与590nm的绿波段的像元灰度值，以及波长范围在770nm与890nm的近红波段的像元灰度值获取归一化水体指数NDWI。

本发明优选实施方式中，步骤S2，采用下述公式表示归一化蓝色水体指数：

其中，

NBWI用于表示归一化蓝色水体指数；

具体地，由于利用归一化水体指数NDWI进行水体提取时，很难区分出水体和暗色建筑，而暗色建筑在蓝波段区间具有强的反射率，其他地物又在蓝波段区间的反射率较低，因此，于步骤S2中，通过波长范围在450nm与520nm的蓝波段的像元灰度值，以及波长范围在770nm与890nm的近红波段的像元灰度值获取归一化蓝色水体指数NBWI。

本发明优选实施方式中，如图2所示，步骤S3包括：

步骤S31，选取多个水体特征点和多个非水体特征点；

步骤S32，根据水体特征点于第一影像上的像元灰度值以及水体特征点于第二影像上的像元灰度值进行线性拟合，获取水体线斜距和水体线截距，构建水体线；

步骤S33，根据非水体特征点于第一影像上的像元灰度值以及非水体特征点于第二影像上的像元灰度值进行线性拟合，获取非水体线斜距和非水体线截距，构建非水体线；

步骤S34，根据水体线和非水体线获取最佳分类线斜距，并根据最佳分类线斜距、水体线和非水体线的交点获取最佳分类线截距，根据最佳分类线斜距和最佳分类线截距构建最佳分类线。

具体地，此处基于距离型指数，构建融合垂直水体指数，即判断像元与最佳分类线之间的距离判断像元是水体或是非水体。

因此，于步骤S31，选取多个水体特征点和非水体特征点，具体地，地物可通常划分为四大类：水体、植被、裸土和不透水层。而在遥感影像中，植被、裸土与水体的光谱特性差别明显，不透水层中的暗色建筑光谱特性和水体光谱特性最为相似，因此于此处，非水体特征点主要为不透水层的暗色建筑。

于步骤S32-步骤S33，采用最小二乘法对所有水体特征点的像元灰度值和所有非水体特征点的像元灰度值进行线性拟合。于步骤S32，所有水体特征点可表示为：A(x

其中，x

y

由此，对水体特征点进行建模生成水体线，可采用下述公式表示：

y

其中，水体线斜距可采用下述公式获取：

其中，水体线截距可采用下述公式获取：

由此，将所有水体特征点带入上述公式(4)和公式(5)，即可获取水体线斜距和水体线截距，以确定构建水体线。

于步骤S33，所有非水体特征点可表示为：B(x

其中，x

由此，对非水体特征点进行建模生成非水体线，可采用下述公式表示：

y

其中，非水体线斜距可采用下述公式获取：

其中，非水体线截距可采用下述公式获取：

由此，将所有水体特征点带入上述公式(7)和公式(8)，即可获取非水体线斜距和非水体线截距，以确定构建非水体线。

最后，于步骤S34，根据水体线y

同时联立水体线和非水体线：

即可获取水体线和非水体线的交点(x

将最佳分类线斜距a

本发明优选实施方式中，可根据上述确定的最佳分类线y

IPWI＝η×NBWI+δ×NDWI (10)

其中，NBWI用于表示所述归一化蓝色水体指数，NDWI用于表示所述归一化水体指数，η和δ分别为最优系数，可用于加大水体和暗色建筑像元之间的区别，以进一步增强水体识别效果。

相应的，η和δ可采用下述公式表示：

也就是说，获取的融合垂直水体指数IPWI可通过下述公式确定：

如图4所示，相应的，于步骤S5，预设一阈值θ，并根据阈值θ和融合垂直水体指数IPWI提取水体，具体地，当融合垂直水体指数IPWI大于阈值θ时，可视为水体像元，而当融合垂直水体指数IPWI不大于阈值θ时，可视为非水体像元。

运用本技术方案，针对如图5所示的高景图像，可获取如图6所示的提取结果，针对如图7所示的最宽河流为5m的高景图像，可获取如图8所示的水体提取结果，针对如图9所示的具有鱼塘、坑塘、细小河流的高景图像，可获得如图10所示的水体提取结果。

由此可见，本发明提供一种水体快速提取方法，通过获取归一化蓝色水体指数能够增强建筑与水体之间的区别，通过获取融合垂直水体指数能够进一步增强水体与非水体之间的像元灰度值的区别，以快速识别并提取细小的水体，提高水体提取准确度，加快水体提取速度。

以上仅为本发明较佳的实施方式，并非因此限制本发明的实施方式及保护范围，对于本领域技术人员而言，应当能够意识到凡运用本说明书及图示内容所作出的等同替换和显而易见的变化所得到的方案，均应当包含在本发明的保护范围内。