典型区域水体提取指数对比研究

一、预处理

对所得Landsat 8 OLI_TIRS影像进行辐射定标和大气校正。通过辐射定标和大气校正后，随后利用济南市界限图进行矢量数据裁剪，得出研究区域。

二、指数计算

本实验采用NDWI、MNDWI两指数提取水体并进行对比研究。

进行波段运算，分别得到NDWI、MNDWI计算结果图。

NDWI=(green-nir)/(green+nir)

MNDWI=(green-swir)/(green-swir1)

              

                                    (MNDWI)                                                      (NDWI)

三、二值化处理

NDWI图像进行二值化处理，大于等于0的为水体(白色)

        

MNDWI图像进行二值化处理，大于等于0.85的为水体(白色)

        

四、决策树分类

结合两幅图像以及真彩色影像来看，能够明显看出，NDWI图像中存在大量城区、村落、裸地等地物被误判为水体(白色)的情况。相比之下，MNDWI图像更加贴合实际，误判程度相比较小。

为详细比较研究两指数的水体提取情况，采用决策树分类的方式，提取出NDWI图像错判的区域(蓝色)。即NDWI>0且MNDWI<=0.85的区域。

        

五、典型区域比较

在此，以济南市老城区(泉城路大明湖区域)为典型区域，进行分析。

         

(真彩色影像)                                            (NDWI图像)

            

(NDWI图像)                                               (决策树图像)

结合实际，如上图所见，在典型区域中NDWI指数错误将城区及周边建筑分类为水体。根据实际情况来看，可能是由于建筑物高度以及阳光照射后的光线反射等原因，以及NDWI指数本身的缺陷，导致误判。

相比之下，MNDWI指数的精确度较高，但本实验中所取的“0.85”这一界线仅仅是基于本图像的一个结果，不具有普遍性。