基于Landsat8/TIRS的北京城市热岛研究（附练习数据下载）

1．实验目的和内容

1.1 实验目的
基于Landsat8/TIRS数据反演地表温度，在此基础上研究北京城市热岛效应。掌握Landsat8数据的亮度温度计算、比辐射率计算、地表温度反演以及城市地表热岛分析等内容。

1.2实验数据

1. 北京市行政区划数据
2. 2017年7月北京市Landsat8数据
   文末有练习数据下载

1.3 实验内容
1.数据预处理
2.土地覆盖分类
3.辐射定标与亮温计算
4.地表比辐射率计算
5.大气透过率计算
6.地表温度计算
7.城市热岛效应分析

2．实验内容

2.1数据预处理
在ENVI中将数据转换为.tif格式后在arcmap里面进行进行镶嵌、裁剪、掩膜处理后如下图所示：
图 1：北京市镶嵌掩膜后OLI数据

2.2土地覆盖分类

2.2.1地物特征的选取
由于提供的OLI数据分辨率低、地物模糊、可识别度低，为此我在arcmap软件中使用了在线地图小插件对比地物进行建筑、裸地、水体、林地和农田这5大类地物训练区的划分，如下图所示：
图 2：地物特征选取过程图

2.2.2监督分类
将采集的地物特征训练数据输出为特征文件*.GSG，并且采用最大拟然分类进行对北京土地覆盖的监督分类，最后进行符号化后制图输出，如下图所示：

图 3：北京市土地覆盖分类结果图

2.3辐射定标与亮温计算

使用波段运算公式“3.34420E-04*b1+0.1”先计算得到TIRS第1、2波段的辐射亮度然后继续使用波段运算公式分别运算得到TIRS第一波段和第二波段的亮度温度。

2.4地表比辐射率计算

2.4.1计算植被覆盖度
使用波段运算公式“(float(b5)-b4)/(float(b5)+b4)”计算北京市NDVI指数。接着计算NDVI数据统计直方图，掩膜数据为北京市范围，读取出各年份累积频率为5%和95%的NDVI值分别作为裸植被像元NDVIv和纯植被像元NDVI0，使用波段