第一次使用spot进行水体提取，这里提供一个简单的方法，就是通过NDWI指数进行水体提取，但是spot影像和其他遥感数据有些区别，所以在提取过程中，需要注意一些问题。

首先进入spot影像下载地址网站进行spot影像下载，模糊查询感兴趣区域的范围，进行影像筛选，因为后续需要通过波段进行指数计算，在这里我下载了HRVIR1和HRVIR2传感器的多光谱影像。

在下载完影像之后，需要对影像进行镶嵌，将多个单幅影像拼接成一个完整的影像。 首先将所有需要镶嵌的影像加载到ENVI中，选择Map – Mosaicking – Georeferenced工具，将所有的影像导入，选择输出，导出镶嵌完成的影像，将影像加载到ENVI中，看下效果。

镶嵌成一幅影像后，但不是我们的研究区域，所以需要对影像进行裁剪。 在ENVI中打开镶嵌完的影像，选择File – Save as – Save as tiff 工具，进入下图界面，选择自己需要裁剪区域的shp文件，点击OK，即可将感兴趣区域进行裁剪。

对裁剪完的影像进行NDWI值的计算，打开Band Math波段运算工具，输入NDWI归一化差异水体指数的公式，其中b1我们选择Green波段，b3选择NIR波段，即可进行NDWI指数的计算。 在这需要注意一个问题，SPOT5影像的波段列表与下载数据的波段排序不太一样按波段本应该为: B3——绿色（0.50-0.59 微米） B2——红色（0.61 - 0.68 微米） B1——近红外（0.78- 0.89 微米） B4——短波红外（1.58 - 1.75微米） 全色——（0.48 - 0.71微米） 但是在下载的影像中: B1——近红外（0.78- 0.89 微米） B2——红色（0.61 - 0.68 微米） B3——绿色（0.50-0.59 微米） B4——短波红外（1.58 - 1.75微米） 全色——（0.48 - 0.71微米） 所以在波段运算选择时要注意。

计算出NDWI值后，通过像元值分布直方图，可以大致判断出水体的值要小于-0.25，所以我们用这个基准进行水体提取。 本来准备使用阈值法进行提取，但是输出的影像总是全为灰色，所以换种方式。 右键上一步输出的NDWI影像，选择New Raster Color Slices，按顺序点击下图的选项，在右边的弹出框中设置最大最小像元值，完成后点击OK，在生成的Slices文件夹右键，选择保存到shp文件，即可将水体范围进行导出。 将输出的水体范围在Arcgis中打开，整体提取效果良好。

本文只是简单的通过NDWI指数进行了水体的提取，在阈值的设置方面不太标准，具体操作中，可以通过DN值进行判断阈值范围，进行更精准的提取。