地表温度反演——以Landsat8数据为例-20210324.docxVIP

PAGE 19 地表温度反演——以Landsat8数据为例 地表温度是表征陆地表面能量平衡和气候变化的重要指标，是区域和全球尺度地表物理过程的一个关键参数，能提供地表能量平衡状态的时空变化信息，在气象、水文、地质、生态、城市环境和灾害监测等众多研究领域都有广泛应用需求。本章基于Landsat8卫星遥感数据，利用PIE-Basic软件对河南省孟津县地区地表温度进行反演。 主要内容 本次课程将利用PIE-Basic软件，使用Landsat8卫星影像数据进行地表温度反演，包括图像预处理（辐射定标、大气校正、图像裁剪）、NDVI计算、地表比辐射率计算、同温度下黑体辐射亮度计算、地表温度计算等操作。 学习目标 （1）了解基于Landsat8数据的地表温度反演流程； （2）熟练掌握PIE-Basic软件图像预处理、波段运算等操作，了解基于Landsat数据进行地表温度反演的方法。 实验区域与数据 本案例选择的Landsat8卫星于2013年10月17日所拍摄的河南省孟津县区域影像数据，详细信息如下表所示： 表 STYLEREF 1 \s 3 SEQ 表 \* ARABIC \s 1 1实验影像数据列表 数据名称 分辨率 波段数量 坐标系 Landsat8影像1景 30m 11 WGS84 Landsat8卫星图像文件内容如下： 操作流程 基于Landsat8影像的地表温度反演流程如下图所示： 图 STYLEREF 1 \s 4 SEQ 图 \* ARABIC \s 1 1基于Landsat8卫星影像的地表温度反演流程 详细操作步骤 图像预处理 波段合成 波段合成将Landsat8数据中的1、2、3、4、5、6、7波段进行波段合成，得到一幅多波段合成图像（多光谱影像）。 在“常用功能”标签下的“图像运算”组，单击【波段合成】工具，弹出“波段合成”对话框，如下图所示： 图 STYLEREF 1 \s 5 SEQ 图 \* ARABIC \s 1 1波段合成界面 文件选择：点击对话框右上角“…”图标，打开“文件选择”对话框，添加待合成波段影像，即添加Landsat8数据中的1-7波段影像，数据选择完成后，单击【确定】返回“波段合成”对话框； 在波段合成对话框，通过单击文件选择显示框右侧的“∧”“∨”按钮可以调整波段顺序，将波段顺序调整为“1、2、3、4、5、6、7”，单击“x”按钮可以删除选中的波段； 输出文件：设置输出文件及其保存路径。 所有参数设置完成后，单击【确定】，软件进行波段合成处理。 辐射定标 辐射定标是将记录的原始DN值转换为大气外层表面反射率（或称为辐射亮度值）。用户需要计算地物的光谱反射率或光谱辐射亮度时，或者需要对不同时间、不同传感器获取的图像进行比较时，都必须将图像的亮度灰度值转换为绝对的辐射亮度，这个过程就是辐射定标。目的是为了消除传感器本身的误差，确定传感器入口处的准确辐射值。 在“图像预处理”标签下的“辐射校正”组，单击【辐射定标】工具，弹出“辐射定标”对话框，如下图所示： 图 STYLEREF 1 \s 5 SEQ 图 \* ARABIC \s 1 2辐射定标设置界面 输入文件：在此选择波段合成后的影像文件； 元数据文件：输入原始数据文件夹LC81250362013290LGN01下的_MTL.txt文件； 定标类型：选择表观反射率/亮温； 输出文件：设置输出文件名及文件保存路径。 所有参数设置完成后，单击【确定】，软件进行辐射定标处理。 大气校正 大气校正的目的是消除大气对太阳辐射、大气对目标辐射产生的吸收和散射作用的影响，从而获得目标反射率、辐射率、地表温度等真实物理模型参数。大多数情况下，大气校正同时也是反演地物真实反射率的过程。 PIE-Basic软件的大气校正功能是基于6S大气辐射传输模型。6S模型假定在无云大气的情况下，考虑了水汽、CO2、O3和O2的吸收、分子和气溶胶的散射以及非均一地面和双向反射率的问题。光谱积分的步长为2.5nm，可以模拟机载观测、设置目标高程、解释BRDF作用和临近效应，增加了两种吸收气体的计算（CO、N2O）。采用SOS (successive order of scattering) 方法计算散射作用以提高精度。在PIE-Basic软件中，对经过辐射定标后的多光谱数据进行大气校正处理，得到多光谱表观反射率图像。 在“图像预处理”标签下的“辐射校正”组，单击【大气校正】工具，打开“大气校正”设置对话框。 图 STYLEREF 1 \s 5 SEQ 图 \* ARABIC \s 1 4增强型Frost滤波参数设置界面 数据类型：此步骤选择的数据类型需要与输入数据类型一致，在此输入影像文件为上一步辐射定标后得到的表观反射率影像，在此选择表观反射率； 输入文件：在