基于单窗算法反演地表温度的ENVI操作教程

1、单窗算法反演地表温度教程1.1 算法原理1.1.1 单窗算法单窗算法（MW算法）是覃志豪于2001年提出的针对TM数据只有一个热红外波段的地面温度反演算法。经过众多学者验证，单窗算法具有很高的反演精度，且同样适用于ETM+和landsat8数据。公式如下：Ts=（a6（1-C6-D6）也（146-口6）C6D6）Tsensor-D6Ta）/C6式中，LST为地表温度（K）,Tsensor是传感器上白亮度温度（K）,Ta是大气平均温度（K）;a、b为参考系数，当地表温度为0-70C时，a=-67.355351b=0.458606C、D为中间变量，计算公式为：C=2J=（1-r）l+（l-f）r式

2、中，E为地表比辐射率，雷为地面到传感器的大气总透射率。因此单窗算法反演地表温度的关键是计算得到亮度温度Tsenso、地表比辐射率反大气透射率下和大气平均作用温度Ta。1.1.2 参数计算 辐射亮温计算利用Planck公式将图像像元对应传感器辐射强度值转换为对应的亮度温度值。公式如下式中，Tsenso为亮度温度值；J影像预处理后得到的光谱辐射值,，2、单位为w/(m,sjNm),K1、K2为常量，可由数据头文件获取。计算图像辐射亮温之前，需采用辐射定标参数将像元灰度值DN转换为热辐射强度值，公式如下：.AX-式中，Ml为增益参数，Al为偏移参数，该参数可直接在影像通文件数据中获取，

3、且ENVI软件中已经集成，不需要自己在查找。 地表比辐射率计算根据覃志豪针对TM影像提出的混合像元分解法来确定区域地表福辐射率。对于城市区域，我们简单的将其分为水体、自然表面和建筑表面三种，因此针对混合像元尺度上的地表比辐射率通过下式来估算：£=PR£+(1尸)R+/rV1'、1'7nimt式中，色为混合像元的地表比辐射率；Pv为植被覆盖率；Rv为植被的温度比率；Rm为建筑表面的温度比率；£v表示植被法地表比辐射率，¥m表示建筑表面的地表比辐射率；*表示辐射校正项。根据覃志豪经验公式，v=0.986;%=0.972根据Ast

4、er提供的常用地物比辐射率光谱库Landsat8数据特点，对于B10波段，fv=0.98672;后m=0.96767。对于B11波段，v=0.98990;%=0.977515。在地表相对较平整的区域，一般可取de=维地表起伏较大的区域,d£可根据植被的构成比例估计。研究表明热辐射相互作用在植被与裸土分别占一半时达到最大，所以可以根据如下经验公式来估计d£：当Pvw0.5日Id&=0.0038Pv当Pv>0.5时,d&=0.0038(Pv)当Pv=0.5时,d£=0.0019RV和Rm的确定采用覃志豪等提出的经验公式计算：RV=0.9332+0

5、.0585*PVRm=0.9886+0.1287*PVPV为植被覆盖率可用下式计算：PV=(NDVI-NDVIs)/(NDVIv-NDVIs)其中NDVI是归一化植被指数，NDVIV、NDVIS分别是植被和裸土的NDVI值,由于没有详细的区域植被和土壤光谱资料，采用覃志豪经验值以NDVIV=0.70和NDVIS=0.05来进行植被覆盖度的近似估计。 大气平均作用温度的计算对于大气平均温度，覃志豪等也提供了一个经验公式，该公式是根据Modtran所提供的标准大气推导出来的，如下表所示，其中T0为近地表温度（K）。袤3T大气平均作用温度与地面附近气温的关系叫大气的剖间类型兀与”的关系

6、美国1网6年平均大气兀二二5.为96+0W的45To拉带平均大气丁17976£+0号1715T口中纬度夏季平均人气Ta=160110+0.92621To中纬度冬季平均大气T,-1$.2704+0J1118Z；般研究区在重庆的，选择中纬度夏季平均大气 大气透射率的计算由于单窗算法对输入的参数较敏感，因此需要较精确的大气透射率参数来反演地表温度。影响大气透射率的因素较多，如气压、气温、气溶胶含量、大气水分含量、03、CO2、CO、NH4等对热辐射传导均有不同程度的作用，从而使地表的热辐射在大气中的传导产生衰减。因此，准确的大气透射率求算比较复杂，需要较详细的大气剖面数据，但

7、在实际应用中很难获取详细的大气剖面数据。研究表明大气水分含量对大气透射率的变化影响较大，而其它因素由于很难在短时间内发生较大变化而对大气透射率的影响较小。因此，水分含量就成为大气透射率估计的主要考虑因素。覃志豪等根据大气中的水汽含量动态变化对大气透射率的影响程度进行了研究，并对大气透射率进行了估计，如表3.1。其中:w是大气水汽含量。以重庆为研究区的研究一般选用是气温较高时的估算公式大气透射率估计方程大气剖面水分含里w/（gcm-2）高气温低气温0.4-1.61.6-3.00.4-1.61.6-3.0相关系大气透射率估计方程数平方(R2)p6=0.97429O.08007w0.99611r6=

8、1.03141-0.11536w0.99827r6=0.9820O-0.09611w0.99463p6=1.0537914142w0.99899标准误差(SEE)0.0023680.0025390.0033400.002375接下来只需要计算出大气水分含量便可计算得到大气透射率。学者杨景梅根据地面、高空气象要输资料，提出利用地面湿度计算整层大气可降水量，由此计算得出整层大气水汽含量的方法（杨景梅邱金恒，用地面湿度参量计算我国整层大气可降水量及有效水汽含量的方法的研究，大气科学，2002.01）公式如下：=q+q犷式中，为大气水汽含量（g/cm2）,W为整层大气可降水量（cm）；coC1为经验系

9、数，以重庆为研究对象，一般可选取成都的经验值来计算。C0=0.1274;C1=0.6878计算大气可降水量公式如下：式中，W为整层大气可降水量（cm）,e为地面水汽压（hPa）,a0al为经验系数，可从以下公式计算得出（9之33。）（中3号）b03exp(-L?沂+2.74耳40.15)0.04exp(0.6H)-&+也0.17+&(<P>Ha)40.204(p<33.)0.051 (-25画+0.25,0">20")丸-i-1-090020。)0.066""l(p-33)"+4.41式中，。为地理纬度（

10、°）,H为海拔高度。重庆市主城区纬度取29.5,平均高程取0.35km,则有。=0.19604e+0.04691其中：e是绝对水汽压（hPa）。采用WMO在1996年公布的气象仪器和观测方法指南中公布的饱和水汽压E公式。对于纯水面用公式计算，适用于温度-45-60CE=6.112exp17.6270/(243.12+T0)RH=e/E*100推算出:e=6.112exp17.62.1大气丽而畚敝/*I同褊度F黑体福肺臣度计靠I刖表温度“.1.2.1影像预处理(1)辐射定标选择RadiometricCorrection/RadiometricCalibration。在FileSelec

11、tion对话框中，选择数据LC81230322013132LGN02_MTL_Thermal,单击SpectralSubset选择ThermalInfrared1(10.9),打开RadiometricCalibration面板。(2)大气校正本教程选择Flaash校正法。FLAASHAtmosphericCorrection,双击此工具，打开辐射定标的数据，进行相关的参数设置进行大气校正。注意：如果在多光谱数据辐射定标时Scalefactor值忘记设置，可在本步骤中打开辐射定标数时设置singlescalefaceor值为0.1,若已设置，则默认值为1即T0/(243.12+T0)*RH式中

12、，RH为相对湿度，可以通过气象观测数据获得，To为气温(C)。1.2地表温度反演步骤轴射屋尿川-TQu大，校正-!NDviTrn/亮联图窿/*地表比附射串计端一XgtputDatsTyp?ScaleFcctcrL00一RedometricCi1i*，二krScalefactor不能改变，报错。保持默认1即可oJipnlvTIaASMSrilingsOutputFiLcmuc:II值DisplayresultOOKCancel1）InputRadianceImage打开辐射定标结果数据；2）设置输出反射率的路径，由于定标时候；3）设置输出FLAASH校正文件的路径，最优状态：路径所在磁盘空间足够

13、大；4）中心点经纬度SceneCenterLocation自动获取；5）选择传感器类型：Landsat-8OLI;其对应的传感器高度以及影像数据的分辨率自动读取；6）设置研究区域的地面高程数据；7）影像生成时的飞行过境时间：在layermanager中的Lc8数据图层右键选择ViewMetadata,浏览time字段获取成像时间；注：也可以从元文件“LC81230322013132LGN02_MTL.tXt中找到，具体名称：DATE_ACQUIRED=2013-05-12;SCENE_CENTER_TIME=02:55:26.6336980Z;8）大气模型参数选择：Sub-ArcticSumm

14、er（根据成像时间和纬度信息选择）；9）气溶胶模型AerosolModel:Urban,气溶胶反演方法AerosolRetrieval:2-band(K-T);10)其他参数按照默认设置即可。11)多光谱参数设置中,K-T反演选择默认模式：Defaults->Over-LandRetrievalstandard(600:2100)波谱响应函数：默认指向ProgramFilesExelisENVI51classicfilt_funclandsat8_oli.sli把它重新指向：.ProgramFilesExelisENVI51resourcefilterfuncslandsat8_oli.

15、sli注：这是因为ENVI5.1版本的一个小bug,即Classic中的L8的波谱响应函数不正确，另外一个一劳永逸的方法是：将“filterfuncs中的Landsat8_oli.sli和Landsat8_oli.hdr两个文件拷贝覆盖：“filt_func中的两个文件。否则SWIR1波段大气校正后的结果全为0。见下图.MultispectraISettinqs12)高级参数设置：根据内存大小设置TileSize(Mb):100(8g物理内存)，其他参数默认即可，详细见下图经过FLAASH校正的影像基本去除了空气中水汽颗粒等因子的影响,植被的波谱曲线趋于正常。1.2,1,1地表比辐射率计算(1

16、)植被覆盖度计算第一步，利用多光谱遥感影像，计算研究区域归一化植被指数NDVI,ENVI软件中有工具集成，使用即可。第二步，利用波段运算工具，根据覃志豪经验公式，计算研究区植被覆盖度。(b1gt0.7)*1+(b1lt0.05)*0+(b1ge0.05andb1le0.7)*(b1-0.05)/(0.7-0.05)b1:NDVI(2)地表比辐射率计算第一步，计算de(b1le0.5)*0.0038*b1+(b1eq0.5)*0.0038*(1-b1)+(b1gt0.5)*0.0019b1:植被覆盖度Pv第二步，计算RV与Rm,在波段运算器中分别输入下列式中：Rv:0.9332+0.0585*b1Rm:0.9886+0.1287*b1b1:植被覆盖度Pv第三步，计算地表比辐射率利用B10波段计算地表比辐射率：b1*b2*0.98672+(1-b1)*b3*0.96767+b4利用B11波段计算地表比辐射率：b1*b2*0.9899+(1-b1)*b3*0.97515b1:植被覆盖度Pv;b2:植被表面温度比率RV;b3:建筑表面温度比率Rm;b4:指的是de1.2.2 辐射亮温计算K1K2需要根据不同影像查表得到：(132