ENVI中使用水体指数法NDWI提取水体

        在进行地类解译时，为了提高工作效率，经常会采用自动解译和人工修订相结合的方式进行。因此，经常会用到ENVI或者易康进行地类的自动解译，以NDWI提取水体为例，在ENVI中提取水体流程如下：

（1）数据预处理，理论上来说在使用遥感影像进行科学研究或者工程实践中时首先需要进行遥感影像的预处理工作。但仅从自动解译的角度来讲，影像是否进行预处理对自动解译几乎没有太大影响。

（2）使用Band math工具进行水体指数NDWI的计算。

Band math工具在工具箱中的位置如下：Band algebra>Band math

 （3）下面是对NDWI指数的使用，其中B2我们选择Green波段，B4选择NIR波段，之后即可生成NDWI的指数数据。

（4）经过水体指数NDWI计算之后的结果如下左图，图中高亮的白色部分即为经过NDWI方法增强后的水体部分，其他的植被、裸地等地类均被抑制呈现灰黑色。右图为432波段合成的经过大气校正后的Landsat TM数据。

 

（5）在得到NDWI指数结果之后，通过设定相应的阈值，再次使用Band math工具，进行提取水体结果。这里我们设定的阈值为0.2，即NDWI的值大于0.2的部分即为水体。在Band math中设置B1 gt 0.2,其中gt即代表大于的意思。设置阈值提取结果如下图，可以明显看出影像中的水体几乎被完整提取出，但由于卫星影像质量的问题，左侧的白色为影像原始的条带影响所致。

（6）在得到水体结果之后我们需要将栅格结果转出为矢量成果。栅格转换为矢量采用Classification>post Classification>Classfication to vector工具进行转换。在转换时需要设置DN值，由于我们经过指数法处理之后，将整景影像转换为了二值图，水体的DN值显示为1。

 经过转换之后的结果如下图所示。图中蓝色部分即为我们使用NDWI指数提取到的水体结果。图中可以看到水体部分基本被完全提出，边界效果也较为良好。但仍存在一些小图斑或者误提的现象。而针对小图斑的问题，我们也可以继续采取一些去除小图斑的操作，并且根据自己的需要进行设定最小上图面积从而改善这一结果。

 

 （7）在转换为EVF的结果之后，如果需要转换为GIS中所常用的SHP数据，我们还需要使用Vector>Classic EVF to Shapefile工具进行将其转换为shp数据结果。

        上面即为在ENVI中采用指数法进行提取水体的全部步骤，如果要使用其他的指数方法进行提取如植被、建设用地等其他地类，其原理和步骤是一样的，采用上述的方法均可进行有效提取。下篇将给大家展示如何使用Ecognition软件采用面向对象的方法进行提取水体。