一种通过数据融合提取水体信息的方法

技术领域

本发明涉及水资源管理领域，尤其涉及一种通过数据融合提取水体信息的方法。

背景技术

水是生命之源，水资源对一个国家和社会的发展起着尤为重要的作用。传统的水体信息获取方法主要是对人工野外监测和水文监测站获取的数据进行分析。虽然野外监测数据精度最高，但是野外监测周期长、成本高、监测范围有限。此外，水体的季节性变化大，年际变化明显，这些都不利于获取长时序、大范围的水体信息。

遥感技术在监测水体方面的应用越来越广泛。遥感影像覆盖范围广、周期短，能够获取高精度、高时空分辨率的水体信息，具有不受研究区地形的限制且不会对研究地点造成破坏等优点。目前，提取水体的技术按照影像类型可以分为利用光学影像提取、SAR影像提取以及将两种影像融合进行提取等。光学影像具有丰富的波段信息，但是影像受天气状况的影响，在多云雷雨等天气状况下，影像质量差，无法提取水体信息；SAR具有全天时、全天候监测的能力，并且能够穿云透雾，不受天气的影响，但是SAR在提取水体信息的时候易受地形以及相干斑噪声的影响，造成提取结果精度较低。

针对上述问题，本发明提出一种通过数据融合提取水体信息的方法。该方法通过数据融合，建立光学影像计算得到的水体指数与SAR影像预处理得到的后向散射强度系数之间的关系，克服光学影像受天气影响的限制和SAR影像噪声较多的问题，能够获取高精度、高时空分辨率的水体信息。通过该方法，能够实现水体的识别和定位，获取水体的空间位置、水体面积的时空变化规律，能够为水资源的监测和管理提供依据，对水利、资源环境等部门提供决策服务。

发明内容

本发明的目的是为了将两种遥感影像进行融合获取高时空分辨率的水体信息，获取水体的时空变化特征，从而提出了一种通过数据融合提取水体信息的方法，能够有效的将两种数据融合，获取高精度的水体信息。

为了实现上述目的，本发明提出的提通过数据融合提取水体信息的方法步骤如下：

S1、获取覆盖研究区范围的光学影像数据和SAR观测影像数据，为了获得较好的光学影像数据，光学影像当天的天气状况需为晴天；

S2、对所述S1中的光学影像数据进行辐射定标、大气校正和几何校正的数据预处理。辐射定标是建立地物光谱反射率与传感器获取的DN值之间的定量关系。主要依据的公式：

L＝M×Q+A

其中，L-转换后的辐亮度值，M-图像增益，Q-卫星载荷观测值(DN值)，A-图像偏置。大气校正的目的是消除大气层中的各种气体和杂质的影响，获取地物真实的反射率和辐射率。

S3、基于所述S2中的预处理结果，获取每个像元上的地表反射率值，通过波段计算，获取水体指数，选取提取效果最佳的水体指数图；

S4、将所述S1中的SAR观测影像数据进行辐射定标、地形校正、地理编码和滤波预处理，获取地物的后向散射系数图。本步骤中辐射定标的目的是将影像上的强度值转换为后向散射系数。辐射定标的基本公式为：

其中σ

S5、利用滑动窗口回归拟合的方法，将所述S3中确定的水体指数图与所述S4中获得的后向散射系数图建立关系并进行数据融合，获得水体信息。

优选的，所述S2中的大气校正预处理采用FLAASH方法。

优选的，所述S3中的水体指数采用Otsu阈值法对选取的归一化差异水体指数NDWI、改进的归一化差异水体指数MNDWI、自动水体提取指数AWEInsh和AWEIsh进行二分类，基于GoogleEarth影像对提取结果进行验证，选取提取效果最好的水体指数图。

优选的，所述S4中选用RD模型进行地形校正，选用SRTMDEM数据进行地理编码，不仅可将SAR数据从斜距坐标系转换到地理坐标系，而且可消除透视、收缩、顶底位移等现象。SAR影像固有的相干斑噪声对水体信息的提取也会造成很大的干扰，本发明采用的Lee滤波来消除相干斑噪声。为了更好地区分水体和非水体，最后将经过预处理的两幅影像上的像元值由线性刻度转为对数刻度(dB)。

优选的，所述S2中对光学影像数据进行预处理可选用ENVI或SNAP软件，所述S4中对SAR观测影像数据预处理可选用SNAP、ENVI/Sarscape、PIE-SAR、PoLSAR软件中的一种。

优选的，还包括有对所述S5中数据融合结果进行精度验证，验证的指标分别为整体精度、用户精度、产品精度和Kappa系数。

与现有技术相比，本发明的有益效果为：

(1)本发明基于光学影像和SAR影像数据，利用滑动窗口回归拟合的方法，建立水体指数与后向散射系数之间的相关关系，提取水体信息。这种方法能够有效克服在水体提取时光学影像受天气状况的限制和SAR影像斑点噪声较多的缺点，使两种遥感数据在水体提取中得到更充分的利用。

(2)本发明能够精确提取高精度、高时空分辨率的水体信息，掌握水体在长时序中的时空变化规律。此外，遥感影像的图幅大，覆盖区域广，在低成本的条件下，可以高效快速的获取大范围的水体信息。掌握水体的时空分布和变化规律，建立水资源数据库，为水利、资源管理等部门合理规划水资源的有效利用提供理论依据。

附图说明

图1为本发明提出的一种通过数据融合提取水体信息的方法的流程示意图；

图2为利用本发明的方法提取的水体信息精度验证图；

图3为利用本发明的方法提取的水体结果图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。

本发明提出一种通过数据融合提取水体信息的方法，具体的方法步骤如下：

S1、获取覆盖研究区范围的两景光学影像数据，其编号分别为L

S2、将两景光学影像L

S3、计算水体指数。基于光学数据预处理结果，最终获取每个像元上的地表反射率值，通过波段计算，获取水体指数。本发明选取4种运用比较广泛水体指数，主要包括：归一化差异水体指数NDWI、改进的归一化差异水体指数MNDWI、自动水体提取指数AWEI

AWEI

AWEI

其中G-绿光，B-蓝光，NIR-近红外，MIR-中红外，SWIR1、SWIR2-短波近红外。

选用Otsu阈值法进行二分类，将研究区的地物分为水体和非水体。将提取结果叠加GoogleEarth影像进行验证，选取提取效果最好的水体指数图。

S4、将SAR影像S

其中σ

S5、将光学影像的水体指数图与SAR影像的后向散射系数图进行影像融合。首先对S

在第j个滑动窗口中，自变量Yj为：

在自变量图层X1，X2，…Xn第j个滑动窗口中，对应的自变量Xij(1≤i≤n)为：

将Yj，Xij分别变成列向量yj，xij，

假设k为系数图层，C为常数项图层，根据拟合方程：

最终获得n个中心像元的系数矩阵图层，如第i个自变量第j个窗口对应的系数矩阵kij：

某个常数项图层，在第j个窗口对应的常数矩阵Cj:

最终得到n个系数图层k(k

S

Y′＝kX′+C

计算得到拟合值Y'。

S6、将融合结果进行精度验证。将拟合得到的Y'通过Otsu阈值法进行二分类，将地物分为水体和非水体。基于GoogleEarth影像，对光学影像L

需要注意的是，本发明中光学影像的处理软件可以选用ENVI、SNAP其中的一款软件。而雷达影像的预处理软件可以为常用的SNAP、ENVI/Sarscape、PIE-SAR、PoLSAR等其中的一款软件。而上述步骤S2、S3、S4中的计算都可以通过ENVI、SNAP、PIE-SAR等软件计算所得，S5中的计算可以通过ArcGIS、Matlab、Python、C++等计算获得。

下面通过利用本发明提出的方法应用于实际案例中进一步的对本发明的可行性进一步补充证明。

本方案先获取张家口地区的两景光学影像数据Landsat8OIL(数据时间分别为：2019-07-30和2019-08-15)和两景SAR影像数据Sentinel-1(数据时间分别为：2019-08-02和2019-08-14)。利用本发明提出的数据融合方法，将两种数据进行融合，获取高精度的水体信息。主要步骤如下：

步骤一：利用ENVI软件对Landsat8OIL数据进行预处理，预处理的过程包括辐射定标、大气校正和几何校正。利用Landsat8自带的定标文件进行辐射定标；大气校正选取FLAASH工具；选取控制点进行几何校正，控制点均匀分布在图像上，校正误差控制在0.5个像元内。通过数据预处理，获得Landsat8地表反射率数据。

步骤二：对Landsat8地表反射率数据进行波段计算得到水体指数图像(NDWI、MNDWI、AWEI

步骤二：利用SNAP软件对Sentinel-1数据进行预处理，预处理的过程包括辐射定标、地形校正、地理编码和滤波处理。其中辐射定标参数从Sentinel-1数据中获取，地形校正的选用RD(RangeDoppler)模型，选用SRTMDEM作为地形纠正数据，选用Lee滤波算法(滤波窗口：7×7)进行滤波处理，最终将后向散射系数值由线性刻度值转为dB值。通过数据预处理，最终获得VV、VH两种极化方式下的后向散射系数图(编号VV、VH)，空间分辨率为20m。

步骤三：通过Arcgis软件将VV、VH极化方下的后向散射系数图利用最近邻法进行重采样，重采样后的分辨率为30m。通过栅格计算，计算VV×VV、VH×VH和VV×VH，获得编号为VVVV、VHVH和VVVH三幅影像。

将2019年8月15日的MNDWI水体指数图作为因变量Y，将2019年8月14日的VV、VH、VVVV、VHVH和VVVH作为自变量X。根据本发明提出的滑动窗口回归拟合的方法，在不同的滑动窗口下对Y、X进行回归拟合，最终获得系数图层、常数项图层、决定系数R

以上，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。