Landsat8TIRS地表温度反演.docxVIP

热红外遥感( Infrared Remote Sensing )是指传感器工作波段限于红外波段范 围之内的遥感。 即利用星载或机载传感器收集、 记录地物的热红外信息， 并利用 这种热红外信息来识别地物和反演地表参数如温度、 湿度和热惯量等。 目前有很 多的卫星携带了热红外传感器，包括 ASTER、AVHRR 、MODIS 、TM/ETM+/ TIRS 等。 目前，地表温度反演算法主要有以下三种： 大气校正法 (也称为辐射传输方 程： Radiative Transfer Equation —— RTE)、单通道算法和分裂窗算法。 本实例是基于大气校正法，利用 Landsat8 TIRS 反演地表温度。 基本原理：首先估计大气对地表热辐射的影响 , 然后把这部分大气影响从卫 星传感器所观测到的热辐射总量中减去 , 从而得到地表热辐射强度 , 再把这一热 辐射强度转化为相应的地表温度。 具体实现为：卫星传感器接收到的热红外辐射亮度值 Lλ由三部分组成：大 气向上辐射亮度 L↑，地面的真实辐射亮度经过大气层之后到达卫星传感器的能 量；大气向下辐射到达地面后反射的能量。 卫星传感器接收到的热红外辐射亮度 值 L λ的表达式可写为(辐射传输方程)： Lλ = [ εB(TS) + (1- ε)L↓]τ+ L ↑(1.1) 式中， ε为地表比辐射率， TS 为地表真实温度 (K)，B(TS)为黑体热辐射亮度， τ为 大气在热红外波段的透过率。则温度为 T 的黑体在热红外波段的辐射亮度 B(TS) 为： B(TS) = [L λ - L↑- τ(1- ε)L↓]/ τε(1.2) Ts可以用普朗克公式的函数获取。 TS = K 2/ln(K 1/ B(T S)+ 1) (1.3) 对于 TM，K1 =607.76 W/(m 2* μ m*sr) K，2 =1260.56K 。 对于 ETM+ ，K1=666.09 W/(m 2* μ m*sr) K，2 =1282.71K 。 对于 TIRS Band10 ，K1= 774.89 W/(m 2* μ m*sr) K，2 = 1321.08K 。 从上可知此类算法需要 2 个参数：大气剖面参数和地表比辐射率。大气剖 面参数 在 NASA 提供的网站 ( /) 中，输入成影时 间以及中心经纬度可以获取大气剖面参数。适用于只有一个热红外波段的数据， 如 Landsat TM /ETM+/TIRS 数据。 主要内容就是使用 BandMath 工具计算公式( 1.2 )和公式( 1.3 )，处理 流程如下图所示。 注：OLI 大气校正是可选项，根据前人研究成果，大气校正处理对结果影响不是 很大。 图 1.1 基于大气校正法的 TIRS 反演流程图 1、图像辐射定标 （1） 在 主 界 面 中 ， 选 择 File → Open ， 在 文 件 选 择 对 话 框 中 选 择 LC81230322013276LGN00_MTL.txt 文件，ENVI 自动按照波长分为五个 数据集：多光谱数据（ 1-7 波段），全色波段数据（ 8 波段），卷云波段数 据（9 波段），热红外数据（ 10 ，11 波段）和质量波段数据（ 12 波段）。 2 2 ）在 Toolbox 工具箱中，选择 Radiometric Correction/Radiometric Calibration 。 在 File Selection 对 话 框 中 ， 选 择 数 据 LC81230322013132LGN02_MTL_Thermal ，单击 Spectral Subset 选择 Thermal Infrared1 （10.9 ），打开 Radiometric Calibration 面板。 （3 ）在 Radiometric Calibration 面板中，设置以下参数： 定标类型（ Calibration Type ）：辐射亮度值（ radiance ）。 其他选择默认参数。 （4 ）选择输出路径和文件名 1-Lband10_rad.dat ，单击 OK 按钮执行定标处理。得到 Band10 辐射亮度图像。 2 、 地表比辐射率计算 TIRS 的 Band10 热红外波段与 TM/ETM+6 热红外波段具有近似的波谱范 围，本例采用 TM/ETM+6 相同的地表比辐射率计算方法。使用 Sobrino 提出 的 NDVI 阈值法计算地表比辐射率。 ε=0.004Pv+0.986 （1.4 ） 其中， Pv 是植被覆盖度，用以下公式计算： Pv = ［（NDVI- NDVI Soil ）／（NDVI Veg - NDVI Soil ）］ （1.5） 其中， NDVI 为归一化植被指数， NDVI Soi