根据Landset 8计算植被覆盖度

    之前处理了山东省的Modis合成产品,这次用landset8的原始数据算下植被覆盖度练练手。

简要说下步骤:

  1. 数据集介绍
  2. 辐射定标
  3. 大气校正
  4. 计算NDVI
  5. 导出Excel
  6. 计算植被覆盖度

1.数据集介绍

先看一下landset的参数和各波段的信息,本想用landset7的数据,众所周知03年以后坏掉了,看了下各位大神开发的修复程序,还是决定先不搞它了,8和7的区别如下
根据Landset 8计算植被覆盖度_第1张图片

2.辐射定标

然后预处理呗,辐射定标,大气校正,辐射定标还是很简单的
根据Landset 8计算植被覆盖度_第2张图片

3.大气校正

但是但是!到了大气校正,flaash的参数设置还挺多的,而且导入文件的时候还蹦出来个选项卡,粗略查了一下资料问了下老师结论如下:根据Landset 8计算植被覆盖度_第3张图片
根据Landset 8计算植被覆盖度_第4张图片
这两个分别代表:1.从文件中读取,不同的波段设置不同的参数。
2.对所有波段用一个单一的参数。
老师告诉我根据数据的类型不同,选项不同,比如地理位置和植被类型
根据Landset 8计算植被覆盖度_第5张图片
一开始还忘了设置K-T变换参数报错,在Defaults里面选择OLI(600…),其余K-T设置如图
根据Landset 8计算植被覆盖度_第6张图片
另外网上说要修改波谱库的指向,我修改完发现两次大气校正结果一样
根据Landset 8计算植被覆盖度_第7张图片
处理吧,时间还是蛮久的!
根据Landset 8计算植被覆盖度_第8张图片
第一次校正,发现结果的DN值被扩大了,又是各种找资料,结论是:都变大,算植被覆盖度的话并不会受影响。
根据Landset 8计算植被覆盖度_第9张图片

4.计算NDVI

接下来band math,归一化算NDVI呗,常规操作,另外大家要是有一步计算出来的公式欢迎分享给我,band math是最浪费我时间的操作了,想公式有点困难。。根据Landset 8计算植被覆盖度_第10张图片

5.导出到Excel

根据Excel计算植被覆盖度,我取5%置信区间
根据Landset 8计算植被覆盖度_第11张图片

6.计算植被覆盖度并查看结果

VFC计算公式

输入公式计算,基本完成,这一景图是选的青藏高原一角,时间是18年2月份,所以肯定冰雪覆盖,水体我改成了白色。
根据Landset 8计算植被覆盖度_第12张图片

你可能感兴趣的:(攻略,笔记)