对于Landsat 5 TM6波段，根据Nichol等（Nichol，2005）的计算方法提取其地表温度，计算公式如下：

式中，L6为传感器处的辐射亮度值；DN为像元灰度值；gain、bias为TM6波段的增益与偏置，分别取0.056 322、1.238；Tb为像元亮度温度；Ts为地表温度；K1和K2为定标参数，取值分别为607.76W/（m2·sr·μm）、1 260.56 K；中心波长λ取11.475μm；α取1.438×10-2 mK；ε6为基于TM6波段的地表比辐射率。

Landsat8则基于波段10利用辐射传输方程反演地表温度（Barsi，2014），其计算公式为：

式中，C1、C2为常数，分别取值14 387.7μm·K和1.191 04×109 W·μm4/（m2·sr）；λ10取10.9μm；ε10为TIRS10波段的地表比辐射率；τ10为大气在TIRS10波段的透过率；I10↑和I10↓分别为大气向上、向下辐射亮度。参考中纬度夏季标准大气剖面，依据影像成像时间和中心经纬度，采用插值大气剖面方法，获取τ10为0.91，I10↑为0.68W/（m2·sr·μm），I10↓为1.19W/（m2·sr·μm），本书首先假定水体的地表比辐射率为0.995，分别采用Yu et al.发表于2014年和Sobrino et al.发表于2004年文献的NDVI阈值法获取地表比辐射率ε10和ε6。根据LST算法公式4-3和4-4，在ERDAS中利用空间建模工具建立模型，分别计算TM band 6和Landsat 8 band 10的地表温度。模型如图4-1、4-2。

图4-1　TM band 6地表温度ERDAS计算模型

图4-2　Landsat 8 band 10地表温度ERDAS计算模型(www.daowen.com)