长篇幅详解辐射定标、大气校正、监督分类、掩膜统计、植被覆盖度操作

长篇幅详解辐射定标、大气校正、监督分类、掩膜统计、植被覆盖度操作

一、序言

Landsat卫星是遥感初学者甚至遥感从业者、相关研究学者接触最多的一个遥感数据源。今天我们通过一个示例来学习ENVI一些相关的操作。

1.1 实验区简介

大家在浏览谷歌地球或者谷歌地图的时候，有没有发现陕西北部的卫星图会有上图这样的一个明显的轮廓特征。这个轮廓鲜明的地区就是陕西的吴起县，如图下所示：

1998年，对吴起来说是一个至关重要的年份：这一年，该县率先在全国开启退耕还林探索，在全县范围内推行“封山禁牧、舍饲养畜”政策，一场旨在改变生存环境的生态环境建设在吴起县轰轰烈烈展开。20年，吴起累计退耕还林243.79万亩，因封得最早、退得最快、面积最大，成为全国退耕还林第一县。这个卫星图上清晰可见的轮廓正是人们改造自然的见证。

1.2 数据查询选取

首先打开ArcMAP，加载全国的县级区划矢量，使用属性查找功能选择吴起县：

这样吴起县便被选中，然后右键可以导出选中的吴起县的区划为单独一个SHP文件：


然后加载Landsat的PATH和ROW的矢量文件，标注出PATH和ROW：

我们可以知道覆盖吴起县需要128034和128035两景Landsat数据：

接下来我们可以将当前视图的图框坐标系设置为UTM-48N：

打开编辑，画一个比吴起县大的矩形，作为我们的研究区范围（因为我们需要吴起县和县域外一部分数据做比对，所以就没有直接用吴起县的区划矢量），然后导出这个矩形为一个新的SHP文件：

接下来可以进入USGS查询数据，USGS可以导入研究区的SHP文件，将SHP文件压缩为ZIP文件即可被导入：


接下来选择卫星查询即可，我们本次选用的数据是1997年5月19日的Landsat5和2018年6月14日Landsat8。USGS查询数据的时候，可以先通过贴图的方式来看看选取的数据是否理想，比如云量这些：

二、预处理

下面我们以2018年的两景数据为例来介绍一般的预处理过程，涉及数据裁切、辐射定标、掩膜的使用、大气校正等操作。首先打开ENVI软件打开数据，常规的方法是使用OpenAs里面找到Landsat卫星：

偷懒的方法是将解压后文件夹里的MTL.txt文件直接拖拽进ENVI视窗即可：

切忌自己使用Laystack工具私下将这些波段组合在一起，不仅费事费力，而且自己组合后的文件不带定标参数，做不了后续操作。

2.1辐射定标

常规的方法一般是对整景数据进行定标，但是这样的话，假如研究区比较小，反而会增加数据量和运行时间，所以我们可以先裁剪，再进行辐射定标。
裁剪的时候需要将MASK选项设为YES，这样就是根据实际形状裁剪，默认的NO则是裁剪外接矩形
<