ENVI水体提取流程实践与问题

卫片解译中，水体属于比较好判别的一类。有水体指数、阈值法、监督分类等一系列方法进行提取，今天介绍利用水体指数和ENVI提取水体的过程。

1.所用数据

Landsat8卫片1景、ENVI5.3。

2.辐射定标与水陆分离

由于要对卫片进行水体指数的计算，因此需要对数据进行辐射定标以及大气校正的预处理。
辐射定标就是将图像的数字量化值(DN)转化为辐射亮度值或者反射率或者表面温度等物理量的处理过程。在工具栏中搜索Radiometric Correction，点开进行辐射定标。

点击下方的apply flaash settings即可匹配大气校正需要的参数。

最终得到辐射定标的结果，如下：

我们可以看一下水体的光谱曲线，相对较符合。

我们试着在定标的结果上做水体指数。在工具栏中搜索band math，添加(b3-b5)/(b3+b5)方程，其中b3为绿波段，b5为近红外。


结果较好，在-1-1之间，可以进行水体的提取。

绘制水体兴趣区，对兴趣区进行统计分析，得到阈值。



我们选择0.7204作为阈值进行水体的提取（需要多次尝试，找到最好的感觉）。

利用Band Threshold to ROI进行陆水分离。


效果还不错。

3.大气校正与水体提取

大气校正可以去除了空气中水汽颗粒等因子的影响，利用辐射定标后的文件可以进行大气校正。工具栏选择FLAASH Atmospheric Correction。具体参数设置等为题可以参考http://blog.sina.com.cn/s/blog_764b1e9d0101cb2l.html进行设置，有点儿复杂，不再赘述。

结果

我们可以看到，大气校正后，像元值在几百到几千不等，这是由于大气校正的过程中，给各像元值扩大了10000倍，我们可以通过band math对像元值进行还原。但不管还原后还是直接做水体指数，结果都不在-1-1之间，我也不知道为啥。。。。。不过即使这样，我们也可以用上述方法对水体进行阈值划分，提取水体，只不过阈值不在-1-1中取。

4.利用envi自带公式提取水体

envi中有mndwi的自带公式，我们可以利用它对大气校正后的图像做水体指数的计算。
工具栏搜索Spectral Indices，选择大气校正结果。

选择计算mndwi。

利用自带公式算完后，结果依然不在-1-1之间。。。。。。不清楚是应该怎么搞了，这大气校正的结果到底应该砸用呢。不过，也不影响提取水体，直接利用前面的过程提取就好啦。

有大神可以解答一下疑问，大气校正结果到底咋用啊？高程参数咋看啊？膜拜