地表温度反演单通道算法

地表温度反演的单通道⽅法辩异利⽤遥感数据反演地表温度（LST）的物理基础是基于普朗克定律（Planck）量化所构成的热辐射传输⽅程。

根据卫星传感器光谱分辨率设置，反演⽅法分为单波段算法、双波段法（劈窗算法）和多波段算法。

graph TB A[LST反演] -->B(单波段算法) A --> C[双波段算法也称劈窗算法] A --> D(多波段算法) B --> E[辐射传输⽅程] B --> F[单通道算法] B --> G[单窗算法]对于但波段算法中常⽤的辐射传输⽅程法、单通道算法、单窗算法，通过查阅资料他们的主要区别如下。

1、辐射传输⽅程的⼤⽓校正法基本原理是：⾸先估计⼤⽓对地表热辐射的影响, 然后把这部分⼤⽓影响从卫星⾼度上传感器所观测到的热辐射总量中减去, 从⽽得到地表热辐射强度, 再把这⼀热辐射强度转化为相应的地表温度。

2、单通道算法单通道算法（Single-Channel Method, SC）是Jimenez-Munoz和Sobrino在对Planck函数在某个温度值附近作⼀阶Taylor级数展开⽽得出的⼀种普适性单通道算法，该算法可以针对任何⼀种热红外数据反演地表温度。

Jimenez-Munoz, JC & Sobrino, JA. A generalized single-channel method for retrieving land surface temperature from remote sensing data[J]. JOURNAL OF GEOPHYSICAL RESEARCH-ATMOSPHERES, 108(D22):46883、单窗算法单窗算法(Mono-window Algorithm)是覃志豪等根据地表热辐射传导⽅程, 推导出的适⽤于从⼀个热波段遥感数据中推演地表温度的算法，是覃志豪等使⽤中值定理，根据热辐射传输⽅程对Planck函数进⾏线性化⼀阶Taylor级数展开，通过简化⼤⽓向上辐射亮度和⼤⽓向下辐射亮度的计算模型。

基于Landsat 8数据单通道算法反演地表温度——以济南市为例夏安全;齐建国;姜振飞;马津【摘要】当前地表温度反演的遥感数据源多为Landsat TM/ETM+、MODIS数据,Landsat 8热红外数据的使用还不是很多,许多针对Landsat 8热红外数据的地表温度反演算法虽然被提出,但是否能满足不同试验区的精度要求还有待考究,同时,经验模型的使用可能会对求解地表比辐射率和大气透过率等地表温度反演参数造成不同程度的影响.因此,本研究以山东省济南市为研究区,Landsat 8数据为数据源,分别利用大气校正法、JM_SC10算法、TIRS10_SC算法,结合分类回归树(classification and regression tree,简称CART)算法与中等光谱分辨率大气透过率算法计算机模型(moderate spectral resolution atmospheric transmittance algorithm and computer model,简称MODTRAN),构建适合济南地区的温度反演参数,实现地表温度的反演,并以济南市16个气象站的温度数据为基准进行精度验证.结果表明,3种算法反演温度平均误差为1.78℃,TIRS10_SC算法反演精度最高,其次是大气校正法、JM_SC10算法.在一定误差要求下,3种算法均可应用于济南地区的地表温度反演.【期刊名称】《江苏农业科学》【年(卷),期】2017(045)020【总页数】5页(P254-258)【关键词】Landsat8;地表温度反演;单通道算法;大气透过率;地表比辐射率【作者】夏安全;齐建国;姜振飞;马津【作者单位】山东农业大学信息科学与工程学院,山东泰安271018;山东农业大学信息科学与工程学院,山东泰安271018;山东农业大学信息科学与工程学院,山东泰安271018;山东农业大学信息科学与工程学院,山东泰安271018【正文语种】中文【中图分类】S127陆地表面温度(land surface temperature，简称LST)是地表能量平衡和资源环境变化的重要参数。

1、单通道算法模型为：Ts=r*[(y1*L6+y2)/e6+y3]+t2、大气平均作用温度Ta的近似估计温度换算：T＝t＋273.15本图为7月份拍摄，对于中纬度夏季平均大气Ta=16.0110+0.92621T0取乌鲁木齐市平均气温为25摄氏度时Ta = 312.157533、大气透射率t6的估计t6=0.974290-0.08007w，0.4≤w≤1.6。

w为水分含量，单位(g/cm2)，这里，取w=1.0，计算得到t6=0.894224、地表比辐射率的估计典型地表类型的比辐射率ew=0.995ev=0.986em=0.970Pv=[(NDVI- NDVIs)/(NDVIv- NDVIs)]2NDVI 为归一化植被指数, 取NDVIv=0.70 和NDVIs=0（分别取自5%及95%数据）e surface=0.9625+0.0614Pv-0.0461Pv^25、像元亮度温度计算T6=K2/ln(1+K1/L6)其中：K1=607.76，K2=1260.56（覃志豪，用陆地卫星TM6数据演算地表温度的单窗算法）；L6为遥感器接收的辐射强度。

6、遥感器接收的辐射强度计算L6=(15.303-1.238)*b1/255.0+1.238 （覃志豪，用陆地卫星TM6数据演算地表温度的单窗算法）；b1为第六波段像元灰度值（DN值）。

LMAX_BAND6 = 15.303LMIN_BAND6 = 1.238单通道算法模型表达式中, Ts是陆地表面温度, L6是卫星高度上遥感传感器测得的辐射强度(W·m- 2·sr- 1·m- 1),e6是地表发射率; r, t,y1,y2,y3是中间变量, 分别由以下公式计算:R=1/[(c2*L6/T6^2)*(λ^4*L6/c1+1/λ)]y1=0.14714w^2- 0.15583w+1.1234=1.11471y2=-1.1836w^2- 0.37607w- 0.52894=0.27859y3=-0.04554w^2+1.8719w- 0.39071=1.43565t=T6-r*L6其中, C1和C2是Plank 函数的常量, C1=1.19104×10^8W·μm4·m- 2·sr- 1, C2=14387.7 μm·K; T6是卫星高度上传感器所探测到的像元亮度温度, 单位为K; λ是有效作用波长, TM6 有效作用波长为11.457μm;w为大气水分含量在此我们也取w=取w=1.0操作步骤。

热红外遥感（Infrared Remote Sensing）是指传感器工作波段限于红外波段范围之内的遥感。

即利用星载或机载传感器收集、记录地物的热红外信息，并利用这种热红外信息来识别地物和反演地表参数如温度、湿度和热惯量等。

目前有很多的卫星携带了热红外传感器，包括ASTER、AVHRR、MODIS、TM/ETM+/ TIRS等。

目前，地表温度反演算法主要有以下三种：大气校正法（也称为辐射传输方程：Radiative Transfer Equation——RTE）、单通道算法和分裂窗算法。

本实例是基于大气校正法，利用Landsat8 TIRS反演地表温度。

基本原理：首先估计大气对地表热辐射的影响, 然后把这部分大气影响从卫星传感器所观测到的热辐射总量中减去, 从而得到地表热辐射强度, 再把这一热辐射强度转化为相应的地表温度。

具体实现为：卫星传感器接收到的热红外辐射亮度值Lλ由三部分组成：大气向上辐射亮度L↑，地面的真实辐射亮度经过大气层之后到达卫星传感器的能量；大气向下辐射到达地面后反射的能量。

卫星传感器接收到的热红外辐射亮度值Lλ的表达式可写为（辐射传输方程）：Lλ = [εB(T S) + (1-ε)L↓]τ+ L↑（1.1）式中，ε为地表比辐射率，T S为地表真实温度(K)，B(T S)为黑体热辐射亮度，τ为大气在热红外波段的透过率。

则温度为T的黑体在热红外波段的辐射亮度B(T S)为：B(T S) = [Lλ - L↑- τ(1-ε)L↓]/τε（1.2）T s可以用普朗克公式的函数获取。

T S = K2/ln(K1/ B(T S)+ 1) （1.3）对于TM，K1 =607.76 W/(m2*µm*sr)，K2 =1260.56K。

对于ETM+，K1=666.09 W/(m2*µm*sr)，K2 =1282.71K。

对于TIRS Band10，K1= 774.89 W/(m2*µm*sr)，K2 = 1321.08K。

地表温度反演的算法综述作者：张微程武学倪静来源：《绿色科技》2014年第12期摘要：指出了在一些关于地理环境的研究中，地表温度是一个不可或缺的因子，而传统的测量方法无法获知地表温度，随着遥感的应用，对地表温度进行反演可以迅速得知地表温度。

国内外学者提出了多种地表温度反演的算法，主要有单通道法、单通道多角度法、多通道多角度法、昼夜法、分裂窗法等。

其中应用最广泛的是分裂窗法。

关键词：地表温度;反演;算法中图分类号：P237文献标识码：A文章编号：1674-9944（2014）12-0220-031引言地表温度是一个复杂的概念，难以准确测量得出具体数值。

而它在地理学中有很重要的应用，关于地球表面的很多研究中都少不了地表温度这一参数。

如何尽可能准确地获得地表温度这一参数，成为一个很重要的课题。

随着遥感技术的广泛应用，可以大范围监测地球表面，使地表温度的获取有了可能。

许多学者在遥感技术的基础上，提出了一些算法，来反演地表温度。

各种算法各有特点，各有所长。

本文主要对目前各学者研发的地表温度反演算法进行了总结与归纳，方便后续的研究者根据不同研究对象选择合适的算法来解决问题。

2地表温度2.1地表温度概况地表温度（Land Surface Temperature， LST），即地面的温度。

太阳的短波辐射到达地面，一部分被反射到宇宙空间，一部分被地面吸收，被吸收的这部分太阳辐射使得地面增温，经过对地面温度的测量，得到的数值即地表温度。

地表温度在地理学各个圈层中起着十分重要的作用，促进大气圈、生物圈、土壤圈、水圈、岩石圈等的能量交换和物质循环，因而地表温度在地理学的研究中有着广泛的应用。

如在对土壤的研究中可以通过地表温度推算出干旱程度;在对石油、铀矿的探测中可以利用地表温度来寻找;在对城市地表温度的研究可以得出城市的热岛效应概况等。

2.2地表温度反演地表温度的测量可以用热动力学的方法测量（用温度计测量）或辐射学测量（用遥感方法测量）。

《遥感数字图像处理》第六讲实习2 Landsat TM6 地表温度反演EX2 Retrieving Earth Surface Temperature from Landsat TMBand Sixth Imagery一目标1、TM Level 1 数据的热红外波段辐射定标：学会阅读头文件，找出所需定标参数；利用定标参数将TM图像热红外波段DN值转换为辐射亮度；2、运用单通道法，反演地表温度反演。

二要求1、提交实习报告，作为第2次平时成绩；2、按5组分组，组长负责本组同学，组织沟通和交流，并督促完成实习报告；科代表负责全班的实习报告收集。

三地表温度的反演——单窗算法技术流程1、如图所示，卫星传感器接收到的热红外辐射亮度 L由三部分组成：λλλλλλλλτετεo o o o s L L ↓↑-++=)1()T (B L其中，s T 为地表真实温度，)T (B s λ表示温度为s T 的黑体在热红外波段的辐射亮度，Lo ↓λ表示大气向下辐射亮度，L o↑λ表示大气向上辐射亮度，λτo为大气在热红外波段的透过率，λε为地表发射率。

上式移项得到：L L o o o s ↓↑---=λλλλλλλλεετε1L )T (B ）（2、大气参数的确定方法通过模拟大气对辐射传输的影响，可以为计算大气效应提供了一种有效的方法。

由于没有卫星过境时的同步气象数据，不能很好的模拟当时的大气状况。

这里我们参考中纬度夏季标准大气剖面，采用MODTRAN 模拟得到各个大气参数：Lo ↓λ表示大气向下辐射亮度，模拟结果为1.68 Wm -2um -1Sr -1，L o↑λ表示大气向上辐射亮度，模拟结果为1.74 Wm -2um -1Sr -1，λτo 为大气在热红外波段的透过率，模拟结果为0.77。

3、地表发射率λε的确定方法（1）可以根据不同地物类型，赋值给出：（2）使用国外研究者的经验公式（利用λε和NDVI 之间的拟合公式给出） 根据Van 的经验公式：)ln(047.00094.1NDVI +=λεVan 经验公式是在自然地表上总结出来的，在应用于非自然地表地区（如城市地表）时必须进行订正。

利用单通道算法对MERSI数据进行地表温度的反演研究摘要选取2008年8月16日MERSI 1km遥感数据,以新民为中心截取170×170像元区域作为研究区,选用具有普适性的单通道算法对其进行地表温度的反演。

结果表明:MERSI LST数值低于观测值12k左右,低于MODIS LST18k左右。

虽然MERSI LST反演结果未能达到理想偏差范围,但是反演温度的趋势与实际观测值和MODIS LST趋势基本一致,最高值均出现在康平,最低值则都出现在本溪。

关键词单通道算法;MERSI数据;地表比辐射率;大气水汽含量;地表温度1云检测利用VC++对基于HDF5格式的MERSI数据进行提取、数据预处理、云检测,判断所选用数据是否可用。

在有云的情况下,热红外波段得到的不是真实的地表温度。

为了有效地消除云的干扰,从遥感图象上获得真实的地表温度信息,首先要进行云像元识别,即进行云像元的检测,以达到去除云像元的目的。

云在3波段(0.65μm)具有高反射率,而且0.65μm遥感图象对区分陆地和云的边界较好;云在近红外波段18波段(0.94μm)的波谱特征主要与大气中的含水量有关,它主要反映大气中的水汽特征,即水汽吸收谷。

由于在0.65μm与0.94μm处,云与各种地物波谱特征形成明显反差,因此将其归一化处理,不仅可以突出云的信息,而且可以部分消除太阳高度角、卫星扫描角及大气程辐射的影响,归一化云检测指数CDI(Cloud Detection Index)表示为:选用FY3A_MERSI_GBAL_L1_20080903_0145_1000M_ MS.HDF和FY3A_MERSI_GBAL _L1_20080816_0225_1000M_ MS.HDF数据,以沈阳(东经123.4°,北纬41.8°)为中心截取100×100像元。

云的云检测指数为正(CDI>O),土壤的云检测指数接近0(CDI≈0),植被的云检测指数为负(CDI0的像元是9513个。

地表温度反演的算法综述张微;程武学;倪静【期刊名称】《绿色科技》【年(卷),期】2014(000)012【摘要】指出了在一些关于地理环境的研究中，地表温度是一个不可或缺的因子，而传统的测量方法无法获知地表温度，随着遥感的应用，对地表温度进行反演可以迅速得知地表温度。

国内外学者提出了多种地表温度反演的算法，主要有单通道法、单通道多角度法、多通道多角度法、昼夜法、分裂窗法等。

其中应用最广泛的是分裂窗法。

%In some researches that related to geographical environment , land surface temperature is an essential factor .However ,land surface temperature cannot be measured through traditional method .Atpresent ,because of the application of remote sensing ,through land surface temperature retrieval ,land surface temperature can quickly be measured .Domestic and abroad scholars have proposed a variety of algorithms of land surface temperature retrieval , mainly including single-channel algorithm ,single-channel multi-angle algorithm ,multi-channels multi-angle algorithm ,round-the-clock algorithm and split-window algorithm ,etc .One of the most widely used algorithms is split-window algorithm .【总页数】3页(P220-222)【作者】张微;程武学;倪静【作者单位】四川师范大学西南土地资源评价与监测教育部重点实验室，四川成都610068;四川师范大学西南土地资源评价与监测教育部重点实验室，四川成都610068; 中国科学院成都生物研究所，四川成都610081;四川师范大学西南土地资源评价与监测教育部重点实验室，四川成都610068【正文语种】中文【中图分类】P237【相关文献】1.地表温度反演的算法综述 [J], 朱贞榕;程朋根;桂新;腾月;童成卓2.气象卫星在地表温度反演及应用方面的综述 [J], 朝鲁门3.矿区地表温度反演算法比较 [J], 吴蓉蓉4.基于Landsat8数据的地表温度反演算法研究——以乌鲁木齐市为例 [J], 马术;蒲智ndsat-8影像地表温度反演算法对比分析 [J], 朱江山;杨丽萍;白宇兴;苏志强;孔金玲因版权原因，仅展示原文概要，查看原文内容请购买。