根据Landset 8计算植被覆盖度

之前处理了山东省的Modis合成产品，这次用landset8的原始数据算下植被覆盖度练练手。
先看一下landset的参数和各波段的信息，本想用landset7的数据，众所周知03年以后坏掉了，看了下各位大神开发的修复程序，还是决定先不搞它了，8和7的区别如下

然后预处理呗，辐射定标，大气校正，辐射定标还是很简单的

但是但是！到了大气校正，flaash的参数设置还挺多的，而且导入文件的时候还蹦出来个选项卡，粗略查了一下资料问了下老师结论如下：

这两个分别代表:1.从文件中读取，不同的波段设置不同的参数。
2.对所有波段用一个单一的参数。
老师告诉我根据数据的类型不同，选项不同，比如地理位置和植被类型

一开始还忘了设置K-T变换参数报错，在Defaults里面选择OLI（600…）,其余K-T设置如图

另外网上说要修改波谱库的指向，我修改完发现两次大气校正结果一样

处理吧，时间还是蛮久的！

第一次校正，发现结果的DN值被扩大了，又是各种找资料，结论是：都变大，算植被覆盖度的话并不会受影响。

接下来band math，归一化算NDVI呗，常规操作，另外大家要是有一步计算出来的公式欢迎分享给我，band math是最浪费我时间的操作了，想公式有点困难。。
根据Excel计算植被覆盖度，我取5%置信区间

输入公式计算，基本完成，这一景图是选的青藏高原一角，时间是18年2月份，所以肯定冰雪覆盖，水体我改成了白色。