1.实验目的和内容
1.1 实验目的
基于Landsat8/TIRS数据反演地表温度,在此基础上研究北京城市热岛效应。掌握Landsat8数据的亮度温度计算、比辐射率计算、地表温度反演以及城市地表热岛分析等内容。
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,使用波段