基于单窗算法反演地表温度的ENVI操作教程word文档良心出品.docx
- 文档编号:1224828
- 上传时间:2022-10-19
- 格式:DOCX
- 页数:10
- 大小:739.82KB
基于单窗算法反演地表温度的ENVI操作教程word文档良心出品.docx
《基于单窗算法反演地表温度的ENVI操作教程word文档良心出品.docx》由会员分享,可在线阅读,更多相关《基于单窗算法反演地表温度的ENVI操作教程word文档良心出品.docx(10页珍藏版)》请在冰豆网上搜索。
基于单窗算法反演地表温度的ENVI操作教程word文档良心出品
单窗算法反演地表温度教程
1.1算法原理
1.1.1单窗算法
单窗算法(MW算法)是覃志豪于2001年提出的针对TM数据只有一个热红外波段的地面温度反演算法。
经过众多学者验证,单窗算法具有很高的反演精度,且同样适用于ETM+和Iandsat8数据。
公式如下:
Ts(a6(1C6D6)(b6(1C6D6)C6D6)TsensorD6Ta)/C6
式中,LST为地表温度(K),Tsensor是传感器上的亮度温度(K),Ta是大气平均温度(K);a、b为参考系数,当地表温度为0-70C时,a=-
67.355351,b=0.458606C、D为中间变量,计算公式为:
C=ez
D=(l-r)[l-F(1-r)r]
式中,&为地表比辐射率,<为地面到传感器的大气总透射率。
因此单窗算法反演地表温度的关键是计算得到亮度温度Tsenso、地表比辐射率匚、大气透射率•和大气平均作用温度Ta。
1.1.2参数计算1.1.2.1辐射亮温计算
利用Planek公式将图像像元对应传感器辐射强度值转换为对应的亮度温度值。
公式如下
£
hid+41)
式中,Tsenso为亮度温度值;L影像预处理后得到的光谱辐射值,
2
单位为w/(msrm),K1、K2为常量,可由数据头文件获取。
计算图像辐射亮温之前,需采用辐射定标参数将像元灰度值DN
转换为热辐射强度值,公式如下:
厶I二xDN十虫丄
式中,Ml为增益参数,Al为偏移参数,该参数可直接在影像通文件数据中获取,且ENVI软件中已经集成,不需要自己在查找。
1.1.2.2地表比辐射率计算
根据覃志豪针对TM影像提出的混合像元分解法来确定区域地表福辐射率。
对于城市区域,我们简单的将其分为水体、自然表面和建筑表面三种,因此针对混合像元尺度上的地表比辐射率通过下式来估算:
€=PR£+Q—F)/?
e+/
rrv%r7mm 式中,三为混合像元的地表比辐射率;Pv为植被覆盖率;Rv为植被的温度比率;Rm为建筑表面的温度比率;=v表示植被法地表比辐射率,fm表示建筑表面的地表比辐射率;*表示辐射校正项。 根据覃志豪经验公式,匸v=0.986;・讣=0.972 根据Aster提供的常用地物比辐射率光谱库Landsat8数据特点,对于B10波段,匸v=0.98672;m=0.96767。 对于B11波段,匸v=0.98990;Tm=0.977515° 在地表相对较平整的区域,一般可取de=在地表起伏较大的区域,de可根据植被的构成比例估计。 研究表明热辐射相互作用在植被与裸土分别占一半时达到最大,所以可以根据如下经验公式来估计de: 当Pv<0.5时,de=0.0038Pv 当Pv>0.5时,de=0.0038(1-Pv) 当Pv=0.5时,de=0.0019 RV和Rm的确定采用覃志豪等提出的经验公式计算: RV=0.9332+0.0585*PV Rm=0.9886+0.1287*PV PV为植被覆盖率可用下式计算: PV=(NDVI-NDVIs)/(NDVIv-NDVIs) 其中NDVI是归一化植被指数,NDVIV、NDVIS分别是植被和裸土的NDVI值,由于没有详细的区域植被和土壤光谱资料,采用覃志豪经验值以NDVIV=0.70和NDVIS=0.05来进行植被覆盖度的近似估计。 1.123大气平均作用温度的计算 对于大气平均温度,覃志豪等也提供了一个经验公式,该公式是根据Modtran所提供的标准大气推导出来的,如下表所示,其中To为近地表温度(K)。 3-4大气平均作用凰度与地面附近气温舵关系腳 人气的剖而类醴 7;芍珀拘关垂 哭国1976均大气 7,^2505Ta 热帶屮均丸气 T.=17.9765+0.51715r0 屮纬度旻李半均人气 ^=1ClO110+0.9? «21r0 屮纬度冬季平均大气 I;-19.2704+0.91118rD 般研究区在重庆的,选择中纬度夏季平均大气 1.124大气透射率的计算 由于单窗算法对输入的参数较敏感,因此需要较精确的大气透射率参数来反演地表温度。 影响大气透射率的因素较多,如气压、气温、气溶胶含量、大气水分含量、03、C02、CO、NH4等对热辐射传导均有不同程度的作用,从而使地表的热辐射在大气中的传导产生衰减。 因此,准确的大气透射率求算比较复杂,需要较详细的大气剖面数据,但在实际应用中很难获取详细的大气剖面数据。 研究表明大气水分含量对大气透射率的变化影响较大,而其它因素由于很难在短时间内发生较大变化而对大气透射率的影响较小。 因此,水分含量就成为大气 透射率估计的主要考虑因素。 覃志豪等根据大气中的水汽含量动态变化对大气透射率的影响 程度进行了研究,并对大气透射率进行了估计,如表3.1。 其中: w是 大气水汽含量。 以重庆为研究区的研究一般选用是气温较高时的估算公式。 大气透射率估计方程 大气剖面水分含量大气透射率估计方相关系数标准误 r口w/(g•cm-2)程平方(R2)差(SEE) t6=0.974290- 0.08007W t6=1.031412- 0.11536W 041.6 0.99611|0.002368 高气温 1.6-3.0 0.9982710.002539 低气温 041.6 1.6-3.0 t6=0.982007- 0.09611W t6=1.053710- 0.14142W 0.994630.003340 0.998990.002375 接下来只需要计算出大气水分含量便可计算得到大气透射率。 学 者杨景梅根据地面、高空气象要输资料,提出利用地面湿度计算整层 大气可降水量,由此计算得出整层大气水汽含量的方法(杨景梅邱金恒,用地面湿度参量计算我国整层大气可降水量及有效水汽含量的方法的研究,大气科学,2002.01)公式如下: &)=cQ+q职 式中,3为大气水汽含量(g/cm2),W为整层大气可降水量(cm); COC1为经验系数,以重庆为研究对象,一般可选取成都的经验值来 计算。 CO二0.1274;C仁0.6878.计算大气可降水量公式如下: W-坷总 式中,W为整层大气可降水量(cm),e为地面水汽压(hPa),a0al为经验系数,可从以下公式计算得出。 0_03+2.74J7+0.15)(卩>33°) 0.04+£(^<33°) =\ 0.20—£了<3y) J0.05 d.; (p-25.0)J+0.25 fo9>2(y) 乩-\ ■-0.9(^<204) 0.066 ((p-33)2十4.41 式中,©为地理纬度(°,H为海拔高度。 重庆市主城区纬度取29.5,平均高程取0.35km,则有。 3=0.19604e+0.04691 其中: e是绝对水汽压(hPa)。 采用WMO在1996年公布的《气象仪器和观测方法指南》中公布的 饱和水汽压E公式。 对于纯水面用公式计算,适用于温度-45~60C。 E=6.112exp[17.62*To/(243.12+To)] RH=e/E*100 推算出: e=6.112exp[17.62*To/(243.12+To)]*RH 式中,RH为相对湿度,可以通过气象观测数据获得,To为气温「C)。 1.2地表温度反演步骤 ▼ 1.2.1影像预处理 (1)辐射定标 选择RadiometricCorrection/RadiometricCalibration。 在FileSelection对话框中,选择数据LC81230322013132LGN02_MTL_Thermal,单击SpectralSubset选择ThermalInfrared1(10.9),打开RadiometricCalibration面板。 (2)大气校正 本教程选择Flaash校正法。 FLAASHAtmosphericCorrection,双 击此工具,打开辐射定标的数据,进行相关的参数设置进行大气校正。 注意: 如果在多光谱数据辐射定标时Scalefactor值忘记设置,可在本步骤中打 开辐射定标数时设置singlescalefaceor值为0.1,若已设置,则默认值为1即 1)InputRadianeeImage打开辐射定标结果数据; 2)设置输出反射率的路径,由于定标时候; 3)设置输出FLAASH校正文件的路径,最优状态: 路径所在磁盘空间足够大; 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-ArcticSummer(根据成像时间和纬度信息选择); 9)气溶胶模型AerosolModel: Urban,气溶胶反演方法Aerosol Retrieval: 2-band(K-T); 10)其他参数按照默认设置即可。 11)多光谱参数设置中, K-T反演选择默认模式: Defaults->Over-LandRetrievalstandard (600: 2100) 波谱响应函数: 默认指向.. 'ProgramFiles\Exelis\ENVI51\classic\filt_func\landsat8_oli.sli 把它重新指向: ..'Program Files\Exelis\ENVI51\resource\filterfuncs\landsat8_oli.sli 注: 这是因为ENVI5.1版本的一个小bug,即Classic中的L8的波 谱响应函数不正确,另外一个一劳永逸的方法是: 将 “\Exelis\ENVI51\resource\filterfuncs中的”Landsat8_oli.sli和 Landsat8_oli.hdr两个文件拷贝覆盖: “...\ENVI51\classic\filt_func中 的两个文件。 否则SWIR1波段大气校正后的结果全为0。 见下图 12)高级参数设置: 根据内存大小设置TileSize(Mb): 100(8g物理 内存),其他参数默认即可,详细见下图 经过FLAASH校正的影像基本去除了空气中水汽颗粒等因子的影响,植被的波谱曲线趋于正常。 1.2.1.1地表比辐射率计算 (1)植被覆盖度计算 第一步,利用多光谱遥感影像,计算研究区域归一化植被指数NDVI,ENVI软件中有工具集成,使用即可。 第二步,利用波段运算工具,根据覃志豪经验公式,计算研究区植被覆盖度。 (b1gt0.7)*1+(b
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 算法 反演 地表 温度 ENVI 操作 教程 word 文档 良心 出品