不透水面提取1 BCI

1. 对预处理影像进行水体mask

参考这两篇文章尝试做了研究区域的水体掩膜

http://blog.sina.com.cn/s/blog_764b1e9d01019ch8.html

http://blog.sina.com.cn/s/blog_aaeaa2f80102vz1f.html

使用的是MNDWI指数，计算完通过目视解译选择水体的阈值（display / histogram）

ENVI 5.3 ：Toolbox/Raster Management/Masking/Build Mask，mask definition中的options 可以选择各种制作掩膜的方式，此处选择的是import data range, 即阈值范围，输出，水体掩膜制作完成。

此时还需要注意，水体掩膜部分为1，其余部分为0，这是在apply mask中mask value 输入1。

或者进行掩膜反转，band math  是一个很好的工具。（b1 eq 0）*1+(b1 eq 1)*0   然后进行应用。

但是有个很重要的问题，做完掩膜后，其他区域的值好像变了？？（还未找到答案）

2.  进行穗帽变换并计算BCI数值

参考该文章 http://blog.sina.com.cn/s/blog_764b1e9d0101efk9.html#cmt_55D7DA7E-7F000001-9874000F-8B5-8A0

Landsat8 TC变换  与landsat 7 ETM类似，只是8多了第一波段，于是将做完水体掩膜后的文件save as, 单击spectral subset，选择后6个波段即可。 Transform- Tassled Cap,  file type 为 landsat 7 ETM, 输出结果即可。

Brightness TC1 H

Greenness TC2 V

Wetness TC3 L

统计最值，分别计算TC分量，然后计算BCI

3. 统计不透水面的面积

统计栅格面积有很多种方法，可以参考 ENVI下的统计分析功能 http://blog.sina.com.cn/s/blog_764b1e9d0100rdk7.html

这里为了方便，依然选择是统计掩膜文件。通过目视解译，选择出阈值，输出掩膜文件，computer statistics，其中不透水面区域count有1040971个像元，该影像分辨率为30*30， 即该区域的面积为93687.39 hm2。