一种基于Landsat OLI遥感影像的区域性水体制图方法技术

本发明专利技术公开了一种基于Landsat OLI遥感影像的区域性水体制图方法，包括以下内容：1）普适性水体指数模型构建；2）最佳水体指数阈值自动选取；3）区域性水体提取；4）水体提取结果验证；5）区域性水体制图。本发明专利技术在现有水体指数模型基础上，发展了一种适合于Landsat OLI遥感影像的区域水体提取指数模型，采用该模型波段运算后，可以将水体和非水体类型的值域进行有效分离，其中水体表现较为集中的值域范围，并采用了图像熵算法克服时空差异带来的节地物反射率异质性，该算法自动选取最佳水体提取阈值，减少人工确定阈值的不确定性和主观性，提高了水体制图精度和效率。

A Regional Water System Mapping Method Based on Landsat OLI Remote Sensing Image

【技术实现步骤摘要】
一种基于LandsatOLI遥感影像的区域性水体制图方法
本专利技术涉及的是一种区域性水体遥感制图方法，特别是一种基于LandsatOLI遥感影像的区域性水体制图方法，属于遥感制图

技术介绍
随着现代水体制图方法研究的深入和进步，水体提取指数模型的建构方法也在不断提升和改进：1996年Mcfeeters第一次利用绿光(Green，0.52-0.60um)波段和近红外(NIR，0.76-0.96um)波段构建了归一化水体指数(NDWI)，2004年Roger和Kearney利用红光(Red，0.62-0.69um)波段和短波红外(SWIR1，1.55～1.75μm)波段提出了另一种水体提取的指数模型，该模型可以有效鉴别水体和裸土；2005年徐涵秋教授利用绿光(Green，0.52～0.60μm)和短波红外(SWIR1，1.55～1.75μm)构建了归一化差异水体指数模型(MNDWI)，该模型可以有效鉴别水体和低反射率的人工建筑(如老旧房屋和柏油马路等)；2014年Feyisa等利用绿光(Green，0.525-0.600μm)波段、短波红外(SWIR1：1.560-1.660μm，SWIR2：2.100-2.300μm)波段、近红外(0.845-0.885μm)波段和蓝光(Blue，0.450-0.515μm)波段构建了自动水体提取指数模型，该模型可以有效鉴别水体和阴影，适合在城市地区和山地地区的水体提取。现有主要水体提取指数模型如表1所示：表1现有主要水体提取指数模型然而，上述水体提取指数模型的建构存在如下缺陷：(1)现有水体指数模型较少考虑到不同地理空间的遥感影像对水体反射值的影响，大部分水体指数模型阈值选择采用人工方式，导致水体提取结果具有较强主观性，且阈值试探过程费时费工；(2)现有水体指数模型大多没有考虑到季节变化对水体反射值的影响，由于不同季节相同地物的反射率存在一定的差异，采用不变的水体指数阈值提取水体时，将导致提取的水体精度不一致，波动性大；(3)现有指数模型阈值选取时主观性强，水体提取精度随选取阈值的变化而变化。LandsatOLI遥感影像是Landsat系列卫星中Landsat8搭载的传感器--运营性陆地成像仪(OperationalLandImager：OLI)采集的影像，是目前Landsat系列卫星中最新的在轨运行多光谱遥感成像仪。利用Landsat遥感影像提取水体方法较多，其中运算复杂度较低而效率较高的方法发展多种水体指数然后阈值化操作提取水体。因此，将LandsatOLI遥感影像作为水体模型建构和水体图像绘制的基础，能够获得更加精确的水体模型，特别是在区域性水体的建模过程中能够体现其显著的优势。
技术实现思路
本专利技术的目的在于克服现有水体指数模型的基础上，在其基础上提出一种基于LandsatOLI遥感影像的普适性水体指数模型。该模型充分考虑到地理空间位置变化和季节变化对地物反射率的影响，采用了图像熵方法自动调节水体指数阈值，减少人为干预，提高水体提取的精度、效率和稳定性。本专利技术的技术解决方案：基于LandsatOLI遥感影像的区域性水体制图方法，其具体步骤如下：1)普适性水体指数模型构建；定义普适性水体指数模型UWI的指数方程为：其中分别对应遥感影像选取的波段的地表反射值，Mean(bi+bj)表示对于选取波段地表反射值的均值和；bi选取LandsatOLI影像的波段七：band7，SWIR2：2.100-2.300μm或波段六：band6，SWIR1：1.560-1.660μm，bj选取LandsatOLI影像的波段四：band4，Red：0.630-0.680μm或波段三：band3，Green：0.525-0.600μm或波段二：band2，Blue：0.450-0.515μm；公式(1)可以组合为六种指数，分别为：2)最佳水体指数阈值自动选取；针对遥感影像存在的时空差异性，采用图像熵的图像分割方法，对步骤1)中水体指数模型运算后不同时空的影像进行阈值的自动选取，通过迭代选取适合该影像的最佳阈值，实现较高精度水体的提取，具体操作为将图像分割为前景为水体，逻辑值1，背景为非水体，逻辑值为0，其具体步骤如下：(1)将图像划分为k个灰度级别，给定一个初始阈值T0，0≤T0≤k-1；(2)用阈值T0将图像分成两个区域，即前景区域R1和背景区域R0，小于T0为背景，大于等于T0为前景，其每个区域中各灰度级别的概率密度可表示为：(3)计算背景中各灰度级别像素的累计概率：则前景各灰度级别的像素的累计为：P1(T0)＝1-P(T0)；(6)(4)计算背景和前景各对应的图像熵：(5)当阈值为T0时，计算图像总熵：H(T0)＝H0(T0)+H1(T0)；(8)(6)随着阈值变化，分别计算各阈值时对应的图像总熵；(7)比较各阈值时图像总熵，求最大图像总熵为H(Tj)(0≤j≤k-1)，并作为最终的最佳阈值Tj：图像中灰度大于Tj的像元分割为前景，小于等于Tj的像元为背景。3)区域性水体提取；根据不同地理空间区域和不同季节的地物反射率存在的差异，通过步骤2)中得到的图像熵自动选取最佳阈值Tj，对不同区域不同季节的影像自动调节阈值，提取水体的条件为：4)水体提取结果验证；(1)验证设计采用2.5m高分辨率GoogleEarth遥感影像对提取的水体进行验证：在LandsatOLI原始影像区域中随机选取300个点，将这些随机生成点的坐标导入对应年份的GoogleEarth影像中，通过目视解译判别对于该位置的真实地物类型为水体或非水体，并进行计数统计。(2)精度分析对不同时空的遥感影像，随机生成300个验证点，提取各点坐标导入到对应年份的GoogleEarth影像中，当验证点落在云上或阴影上，验证时将这些点作为非有效点排除统计之外；从GoogleEarth影像上识别所有效验证点的类别为水体或非水体，分别统计落在水体和非水体中有效验证点总数，将提取后的影像与GoogleEarth影像配准后逐点验证，统计出水体提取影像中包含落在水体上的正确点数和落在非水体上的正确点数，分别计算并比较各区域水体提取精度，总体精度和Kappa系数。5)区域性水体制图：将对象区域的维度变化和季节变化对水体反射率的影响作为重要考虑因素，进行区域水体提取及制图：提取及制图过程采用Python编程语言，逐步实现原始LandsatOLI遥感图像的读取，水体指数模型波段运算，图像熵算法实现最佳水体提取阈值选取，水体提取和图像输出；其中，水体提取的效果验证采用人工解译计数方式进行，利用ArcGIS软件对原始遥感影像自动生成300个验证点，将各点的坐标导入到GoogleEarth中，将300个点地物类别在GoogleEarth中通过人工解译并计数各地物类别的点数，将点的坐标导入到提取的水体图像中，识别提取的水体中点的个数，并建立混淆矩阵，计算提取精度和Kappa系数。本专利技术的优点：(1)提出一种普适性区域水体指数模型，该模型可以根据选择的波段不同，组合六种指数模型，且每种指数模型均可以有效进行区域性水体提取；(2)有效解决了地理空间变化和季节变化引起水体反射率的时空差异，解决了采用同一阈值对不同影像提取水体时存在精度差异的问本文档来自技高网...

【技术保护点】
1.基于Landsat OLI遥感影像的区域性水体制图方法，其特征是具体包括如下步骤：1)普适性水体指数模型构建；2)最佳水体指数阈值自动选取；3)区域性水体提取；4)水体提取结果验证；5)区域性水体制图。

【技术特征摘要】
1.基于LandsatOLI遥感影像的区域性水体制图方法，其特征是具体包括如下步骤：1)普适性水体指数模型构建；2)最佳水体指数阈值自动选取；3)区域性水体提取；4)水体提取结果验证；5)区域性水体制图。2.根据权利要求1所述的基于LandsatOLI遥感影像的区域性水体制图方法，其特征是所述步骤1)具体包括如下步骤：定义普适性水体指数模型UWI的指数方程为：其中分别对应遥感影像选取的波段的地表反射值，Mean(bi+bj)表示对于选取波段地表反射值的均值和；bi选取LandsatOLI影像的波段七：band7，SWIR2：2.100–2.300μm或波段六：band6，SWIR1：1.560–1.660μm，bj选取LandsatOLI影像的波段四：band4，Red：0.630–0.680μm或波段三：band3，Green：0.525–0.600μm或波段二：band2，Blue：0.450–0.515μm；公式(1)可以组合为六种指数，分别为：3.根据权利要求1所述的基于LandsatOLI遥感影像的区域性水体制图方法，其特征是所述步骤2)具体包括如下步骤：针对遥感影像存在的时空差异性，采用图像熵的图像分割方法，对步骤1)中水体指数模型运算后不同时空的影像进行阈值的自动选取，通过迭代选取适合该影像的最佳阈值，实现较高精度水体的提取，具体操作为将图像分割为前景为水体，逻辑值1，背景为非水体，逻辑值为0，其具体步骤如下：(1)将图像划分为k个灰度级别，给定一个初始阈值T0，0≤T0≤k-1；(2)用阈值T0将图像分成两个区域，即前景区域R1和背景区域R0，小于T0为背景，大于等于T0为前景，其每个区域中各灰度级别的概率密度可表示为：(3)计算背景中各灰度级别像素的累计概率：则前景各灰度级别的像素的累计为：P1(T0)＝1-P(T0)；(6)(4)计算背景和前景各对应的图像熵：(5)当阈值为T0时，计算图像总熵：H(T0)＝H0(T0)+H1(T0)；(8)(6)随着阈值变化，分别计算各阈值时对应的图像总熵；(7)比较各阈值时图像总熵，求最大图像总熵为H(Tj)(0≤j≤k-1)，并作为最终的最佳阈值Tj：图像中灰度大于Tj的像元分割为前景，小于等于T...

【专利技术属性】
技术研发人员：李明诗，王玉亮，
申请(专利权)人：南京林业大学，
类型：发明
国别省市：江苏,32