文档库 最新最全的文档下载
当前位置:文档库 › MODIS数据地表温度反演分裂窗算法的IDL实现_姜立鹏

MODIS数据地表温度反演分裂窗算法的IDL实现_姜立鹏

MODIS数据地表温度反演分裂窗算法的IDL实现_姜立鹏
MODIS数据地表温度反演分裂窗算法的IDL实现_姜立鹏

(完整word版)MODIS数据反演地表温度

表1 MODIS 部分波段及其参数[14] 波段 光谱范围 信噪比 主要用途 分辨率 1 620~670nm 128 陆地、云边界 250m 2 841~876nm 201 陆地、云边界 250m 19 915~965nm 250 大气水汽 1000m 31 10.780~11.280μm 0.05 地球表面和 云顶温度 1000m 32 11.770~12.270μm 0.05 1000m 劈窗算法介绍 McMillin (1975年)最早提出了劈窗算法,最先是用于海面温度的反演,这种方法是利用2个相邻的热红外窗口大气水汽吸收特性的差异,把海面温度表达成2个热红外窗口亮度温度的线性组合。Price (1984年)最先把劈窗算法推广到陆面温度的反演,通过引入比辐射率改正项来减小因陆地表面比辐射率变化而引起的误差。Becker 从理论上证明了用分裂窗技术反演地表温度的可行性,并且第一次从理论上给出了使用分裂窗技术时大气和比辐射率对地表温度反演的影响。Becker 和Li 根据热辐射传导的地方性特征,提出了著名的局地劈窗算法,已得到了较广泛的应用。Wan 和Dozier 在Becker 和Li 的研究基础上,于1996年提出了一种广义的地表温度反演劈窗算法。Sobrino 和Becker 用Lowtran 7对不同的大气、观测角度以及地表参数进行模拟,得出了各参数的表达式。在这些表达式里,大气和比辐射率的作用是耦合在一起的。而Sobrino 等则通过某些近似把这2种作用分开了,通过对大气向下热辐射的近似解和对Planck 辐射函数的线性化。覃志豪等推导了劈窗算法,该算法仅需要2个因素来进行地表温度的演算,即大气透过率和地表比辐射率 [15][ 16] 。在众多的劈窗算法中,覃志豪等提出的算法由于需要参数少、计算简单且精 度较高,被认为是较好的算法之一。本文主要针对这一算法进行介绍。 覃志豪 [15] 等提出的针对MODIS 数据反演地表温度的劈窗算法使用的公式如下: 0131232Ts A AT A T =+- (1) 其中:Ts 是地表温度,31T 、32T 分别是MODIS 第31、32通道的亮温。0A 、1A 、2A 是系数,分别定义如下: 01264.6036368.72575A E E =-+ 1110.440817A A E =++ 220.473453A A E =+

地表温度反演实验报告

遥感原理与及应用 地表温度反演实验报告 专业:地理信息系统 班级: XXXXXXXX 姓名: XXX 学号: XXXXXX 成绩: 指导教师: XXX 2014年12月17日 一. 实验目的 1. 根据实际需要,学会在网上(如中国科学院遥

感与数字地球研究所数据共享网)下载研究区内的遥感数据; 2. 掌握在ENVI中实现简单的地表温度反演的原理与步骤。 二. 实验任务 1. 在中国科学院遥感与数字地球研究所数据共享网上订购并下载覆盖郫县的TM影像; 2. 在ENVI中实现简单的地表温度反演算法。 三. 实验数据 在中国科学院遥感与数字地球研究所数据共享网上下载的覆盖郫县地区的TM影像。

四. 实验原理 图1 TM 影像地表温度反演流程 1. 地表温度(Land Surface Temperature)反演公 式为: 2 1(1)K LST K In R ε=+, 其中,R m DN d =?+,2111607.76K W m sr m μ---=???,21260.56K K =。 2. 根据TM 辐射定标原理,热红外波段表观辐亮 度可以进一步写作: max min 6min 255L L R DN L -=?+, 其中LmaxBand6=15.303 , LminBand6=1.238。 3. 地表比辐射率ε为同温度下地表辐射能与黑体 辐射能的比率,其可以表示为: 1.0090.047(In )(0)NDVI NDVI ε=+>,

其中,4343 TM TM NDVI TM TM -=+,当0NDVI <=时(如水体)地表比辐射率取常数1。 五. 实验步骤 1. TM 数据下载 数 据查询和下载网址https://www.wendangku.net/doc/432076196.html,/query.html ,界面如图2 所示。 图2 中国科学院遥感与数字地球研究所数据共享 网址界面

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

单窗算法反演地表温度教程 1.1 算法原理 1.1.1 单窗算法 单窗算法(M W算法)是覃志豪于2001年提出的针对TM 数据只有一个热红外波段的地面温度反演算法。经过众多学者验证,单窗算法具有很高的反演精度,且同样适用于ETM +和land sat 8数据。公式如下: 式中,LS T为地表温度(K),T sens or 是传感器上的亮度温度(K),T a 是大气平均温度(K );a 、b 为参考系数,当地表温度为0-70℃时,a = -67.355351,b = 0.458606;C 、D 为中间变量,计算公式为: 式中,为地表比辐射率,为地面到传感器的大气总透射率。 因此单窗算法反演地表温度的关键是计算得到亮度温度T senso 、地表比辐射率、大气透射率和大气平均作用温度Ta 。 1.1.2 参数计算 1.1. 2.1 辐射亮温计算 利用Pla nck 公式将图像像元对应传感器辐射强度值转换为对应的亮度温度值。公式如下 6666666666/)))1(()1((C T D T D C D C b D C a T a sensor s -++--+-- =

式中,T senso 为亮度温度值;λL 影像预处理后得到的光谱辐射值,单 位为 )/(2m sr m w μ??,K 1 、K2为常量,可由数据头文件获取。 计算图像辐射亮温之前,需采用辐射定标参数将像元灰度值DN 转换为热辐射强度值,公式如下: 式中,ML 为增益参数,A L为偏移参数,该参数可直接在影像通文件数据中获取,且ENV I软件中已经集成,不需要自己在查找。 1.1. 2.2 地表比辐射率计算 根据覃志豪针对TM 影像提出的混合像元分解法来确定区域地表福辐射率。对于城市区域,我们简单的将其分为水体、自然表面和建筑表面三种,因此针对混合像元尺度上的地表比辐射率通过下式来估算: 式中,为混合像元的地表比辐射率;P V 为植被覆盖率;R V为植被的温度比率;R M 为建筑表面的温度比率; V 表示植被法地表比辐射率,m表示建筑表面的地表比辐射率;d 表示辐射校正项。 根据覃志豪经验公式, V =0.986;m =0.972。

地表温度反演

地表温度反演

目录 一:单窗算法 (3) 1.1实验原理 (3) 1.1.1TM/ETM波段的热辐射传导方程: (3) 1.1.2化简后最终的单窗体算法模型为: (3) 1.1.3大气平均作用温度Ta的近似估计 (3) 1.1.4大气透射率t6的估计 (3) 1.1.5地表比辐射率的估计 (4) 1.1.6像元亮度温度计算 (4) 1.1.7遥感器接收的辐射强度计算 (4) 1.2操作步骤 (5) 1.2.1研究区示意图 (5) 1.3实验结果 (7) 1.3.1灰度图像 (7) 1.3.2密度分割后图像 (7) 二:单通道算法 (8) 2.1实验原理 (8) 2.1.1单通道算法模型为: (8) 2.1.2大气平均作用温度Ta的近似估计 (8) 2.1.3大气透射率t6的估计 (8) 2.1.5像元亮度温度计算 (8) 2.1.6遥感器接收的辐射强度计算 (9) 2.2操作步骤 (9) 研究区示意图 (9) 2.2.1计算L6 (10) 2.2.2T6e6的求算 (10) 2.2.3计算R (10) 2.2.4计算t (10) 2.3实验结果 (11) 2.3.1温度反演灰度图像 (11) 2.3.2密度分割后的图像 (11) 三:辐射方程 (12) 3.1实验过程 (12) 3.1.1数据准备 (12) 3.1.2地表比辐射率的估计 (12) 3.1.3计算同温度下黑体的辐射亮度值 (12) 3.1.4反演地表温度 (13) 3.2温度反演结果 (13)

一:单窗算法 1.1实验原理 1.1.1TM/ETM波段的热辐射传导方程: B6(T6)=t6(q)[ ε6B6(Ts)+(1-ε6)I6~]+I6_ Ts是地表温度; T6是TM6的亮度温度; t6是大气透射率; ε6是地表辐射率。 B6(T6)表示TM6遥感器所接收到的热辐射强度; B6(Ts)是地表在TM6波段区间内的实际热辐射强度,直接决取于地表温度; I6~和I6_分别是大气在TM6波段区间内的向上和向下热辐射强度。 1.1.2化简后最终的单窗体算法模型计算Ts(地表温度): Ts={a(1-C-D)+[b(1-C-D)+C+D]T6-DTa}/C 式中 C6=τ6ε6(ε6为比辐射率,τ6为透射率) D6=(1-τ6)[1+t6(1-ε6)] a =-67.355351,b=0.458606 1.1.3大气平均作用温度Ta的近似估计 温度换算:T=t+273.15 本图为9月份拍摄,对于中纬度夏季平均大气Ta=16.0110+0.92621T0 取平均气温为25摄氏度时Ta = 312.15753 1.1.4大气透射率τ6的估计 τ6=0.974290-0.08007w,0.4≤w≤1.6。 w为水分含量,单位(g/cm2),这里,取w=1.0,计算得到τ6=0.89422

ENVI支持下地表温度反演

[转载]ENVI下利用ETM+数据反演地表温度 (2012-05-15 08:31:18) 转载▼ 标签: 转载 原文地址:ENVI下利用ETM+数据反演地表温度作者:ENVIIDL 地表温度作为地球环境分析的重要指标,而遥感技术作为现代重要的对地观测手段,使得基于遥感图像的地表温度反演的研究越来越多。主要的地表温度反演方法有:大气校正法,单窗算法,单通道法等等。本文介绍用辐射传输方程法对地表温度进行反演。 技术流程: 例子数据为2002年9月2日的襄樊市Landsat ETM+数据。根据数据的特点以及地表温度反演研究的技术要求,采用的技术路线为:先对Landsat ETM+数据进行预处理:数据读取、辐射定标、大气校正、襄樊区域裁剪,利用大气校正,即:辐射传输方程法对其影像热红外波段数据进行操作反演,实现襄樊市地区的地表真实温度的反演研究。具体的处理流程如下:

具体的实现步骤如下: 第一步:准备数据 热红外数据使用的是Landsat的第六波段,已经做了传感器定标、几何校正、工程区裁剪,详细流程参考上面的流程图。文件为TM6-rad-subset-jz-xiangfan.img。

由TM影像(已经过大气校正)生成的NDVI数据,已经利用主菜单->Basic Tools- >Resize Data(SFatial/SFectral)重采样为60米分辨率,与TMi6数据保持一致,文件名为:TM-NDVI-60m.img。 第二步:地表比辐射率计算 物体的比辐射率是物体向外辐射电磁波的能力表征。它不仅依赖于地表物体的组成, 而且与物体的表面状态(表面粗糙度等)及物理性质(介电常数、含水量等)有关,并随着所测 定的波长和观测角度等因素有关。在大尺度上对比辐射率精确测量的难度很大,目前只是 基于某些假设获得比辐射率的相对值,本文主要根据可见光和近红外光谱信息来估计比辐 射率。 (一)植被覆盖度计算 计算植被覆盖度Fv采用的是混合像元分解法,将整景影像的地类大致分为水体、植被和建筑,具体的计算公式如下: F V = (NDVI- NDVI S)/(NDVI V - NDVI S) (2) 其中,NDVI为归一化差异植被指数,取NDVI V = 0.70和NDVI S = 0.00,且有,当某 个像元的NDVI大于0.70时,F V取值为1;当NDVI小于0.00,F V取值为0。 利用ENVI主菜单->Basic Tools->Band Math,在公式输入栏中输入: (b1 gt 0.7)*1+(b1 lt 0.)*0+(b1 ge 0 and b1 le 0.7)*((b1-0.0)/(0.7-0.0))b1:选择NDVI图像 得到植被盖度图像。 (二)地表比辐射率计算

基于劈窗算法的青岛地区地表温度反演研究

龙源期刊网 https://www.wendangku.net/doc/432076196.html, 基于劈窗算法的青岛地区地表温度反演研究作者:孙乐乐金宝轩 来源:《安徽农学通报》2017年第23期 摘要:地表温度对于地球资源环境监测的意义十分重要,在农业、环境生态相关领域研 究以及地表分析研究中都是重要的参数指标。气象测量难以大面积提供面状地温数据,而遥感监测成本低、时相性好、可提供面状数据等特点使其在地表温度监测中具有独特优势。该研究根据Landsat-8遥感数据,通过调整地表比辐射率估算方法后的劈窗地表温度反演算法实现了青岛地区地表温度的反演。并对结果与青岛市气象局网站提供的气象数据进行对比验证,反演结果总体符合实际情况,为沿海城市热能空间分布研究和城区规划提供参考依据。 关键词:地表温度;劈窗算法;Landsat-8 中图分类号 S511 文献标识码 A 文章编号 1007-7731(2017)23-0012-3 Abstract: As the key factor in the study of environment, ecology and land, land surface temperature(LST) plays a vital role in resource and environment monitoring. Comparing with traditional ground meteorological monitoring, remote sensing has more advantages in LST retrieving for it’s spatial continuity and lower cost. With the Landsat-8 data, we adjust the calculation of land emissivity and apply it to retrieve the LSD by a Split-Window algorithm. The analysis result evaluates the effectiveness of the Split-Window algorithm in peninsula regions and the result can provide a reference for the study of urban heat distribution and urban planning. Key words: Land surface temperature;Split-window algorithm;Landsat-8 地表温度LST(Land Surface Temperature)是地球资源环境监测评价中的一项重要指标,它直观反映了区域热能分布,对环境、生态相关领域研究以及地表分析研究都有重要意义。目前,部分学者已进行了相关研究并提出了基于遥感热红外波段的地表温度反演算法,按照所用波段的数量总体上分为单窗算法[1-3]与劈窗算法[4-6]2种。其中,覃志豪单窗算法[1]、Offer Rozenstein劈窗算法[6]是较为代表性的地表温度反演算法。劈窗算法早期多基于MODIS、NOAA/AVHRR等具有2个热红外波段的影像数据实现,其中MODIS数据在地表温度反演中效果较为理想,但由于其分辨率较低的缺陷使得其只能应用于大范围地区的地表温度反演。中小区域的地表温度反演多采用单通道遥感影像进行单窗算法地表温度反演[2,4]。 Landsat-8卫星于2013年发射升空,其携带OLI和TIRS 2个传感器。其中OLI (Operational Land Image)是陆地成像仪,TIRS(Thermal Infrared Sensor)是热红外传感器。TIRS传感器具有2个热红外波段10、11波段,其分辨率为100m,并且其两个热红外波段的波宽及中心波长与MODIS数据相近。因此Landsat-8数据在进行利用劈窗算法进行地表温度反演具有独特优势。本文利用Landsat-8数据基于劈窗算法对青岛地区地表温度进行反演,并结合气象数据对于结果进行评价分析。

landsat 遥感影像地表温度反演 教程(大气校正法)

基于辐射传输方程的Landsat数据地表温度反演教 程 一、数据准备 Landsa 8遥感影像数据一景,本教程以重庆市2015年7月26日的=行列号为(128,049)影像(LC81280402016208LGN00)为例。 同时需提前查询影像的基本信息(详见下表) 标识日期采集时 间 中心经度中心纬度 LC81280402016208LGN002016/7/263:26:56106.1128830.30647…………………………注:基本信息在影像头文件中均可查询到,采集时间为格林尼治时间。 二、地表温度反演的总体流程

三、具体步骤 1、辐射定标 地表温度反演主要包括两部分,一是对热红外数据,二是多光谱数据进行辐射定标。 (1)热红外数据辐射定标 选择Radiometric Correction/Radiometric Calibration。在File Selection对话框中,选择数据LC81230322013132LGN02_MTL_Thermal,单击Spectral Subset选择Thermal Infrared1(10.9),打开Radiometric Calibration面板。

Scale factor 不能改变,否则后续计算会报错。保持默认1即可。 Scale factor 不能改变,否则后续计算会报错。保持默认1即可。 (2)多光谱数据辐射定标 选择要校正的多光谱数 据“LC81230322013132LGN02_MTL_MultiSpectral” 进行辐射定标。 因为后续需要对多光谱数据进行大气校正,可直接单击Apply Flaash Settings,如下图。

landsat遥感影像地表温度反演教程(大气校正法)

基于辐射传输方程的Landsat数据地表温度反演教程一、数据准备 Landsa 8遥感影像数据一景,本教程以市2015年7月26日的=行列号为(128,049)影像(LC81280402016208LGN00)为例。 同时需提前查询影像的基本信息(详见下表) 注:基本信息在影像头文件中均可查询到,采集时间为格林尼治时间。 二、地表温度反演的总体流程 三、具体步骤 1、辐射定标

地表温度反演主要包括两部分,一是对热红外数据,二是多光谱数据进行辐射 定标。 (1)热红外数据辐射定标 选择Radiometric Correction/Radiometric Calibration。在File Selection对话框 中,选择数据LC81230322013132LGN02_MTL_Thermal,单击Spectral Subset 选择Thermal Infrared1(10.9),打开Radiometric Calibration面板。 (2)多光谱数据辐射定标 选择要校正的多光谱数据“LC81230322013132LGN02_MTL_MultiSpectral” 进行辐射定标。 Scale factor 不能改变,否则后续 计算会报错。保持默认1即可。

因为后续需要对多光谱数据进行大气校正,可直接单击Apply Flaash Settings,如下图。 注意与热红外数据辐射定标是的差 别,设置后Scale factor值为0.1。 2、大气校正 本教程选择Flaash 校正法。FLAASH Atmospheric Correction,双击此工具,打开辐射定标的数据,进行相关的参数设置进行大气校正。 注意:如果在多光谱数据辐射定标时Scale factor值忘记设置,可在本步骤中打开辐射定标数时设置single scale faceor 值为0.1,若已设置,则默认值为1即可。 1)Input Radiance Image:打开辐射定标结果数据; 2)设置输出反射率的路径,由于定标时候;

地表温度反演

《地表温度反演》实验报告院系:资源与环境科学学院 专业及班级:地信08-1 学号:20081207019 姓名:李荣立 指导教师:丁建丽 2011年12月

目录 一:单窗算法 (3) 1.1实验原理 (3) 1.1.1TM/ETM波段的热辐射传导方程: (3) 1.1.2化简后最终的单窗体算法模型为: (3) 1.1.3大气平均作用温度Ta的近似估计 (3) 1.1.4大气透射率t6的估计 (3) 1.1.5地表比辐射率的估计 (4) 1.1.6像元亮度温度计算 (4) 1.1.7遥感器接收的辐射强度计算 (4) 1.2操作步骤 (5) 1.2.1研究区示意图 (5) 1.3实验结果 (7) 1.3.1灰度图像 (7) 1.3.2密度分割后图像 (7) 二:单通道算法 (8) 2.1实验原理 (8) 2.1.1单通道算法模型为: (8) 2.1.2大气平均作用温度Ta的近似估计 (8) 2.1.3大气透射率t6的估计 (8) 2.1.5像元亮度温度计算 (8) 2.1.6遥感器接收的辐射强度计算 (9) 2.2操作步骤 (9) 研究区示意图 (9) 2.2.1计算L6 (10) 2.2.2T6e6的求算 (10) 2.2.3计算R (10) 2.2.4计算t (10) 2.3实验结果 (11) 2.3.1温度反演灰度图像 (11) 2.3.2密度分割后的图像 (11) 三:辐射方程 (12) 3.1实验过程 (12) 3.1.1数据准备 (12) 3.1.2地表比辐射率的估计 (12) 3.1.3计算同温度下黑体的辐射亮度值 (12) 3.1.4反演地表温度 (13) 3.2温度反演结果 (13)

Landsat8 TIRS 地表温度反演

热红外遥感(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。 从上可知此类算法需要2个参数:大气剖面参数和地表比辐射率。大气剖面参数在NASA提供的网站(https://www.wendangku.net/doc/432076196.html,/)中,输入成影时间以及中心经纬度可以获取大气剖面参数。适用于只有一个热红外波段的数据,如Landsat TM /ETM+/TIRS数据。 主要内容就是使用BandMath工具计算公式(1.2)和公式(1.3),处理流程如下图所示。

地表温度反演单通道算法

1、单通道算法模型为: Ts=r*[(y1*L6+y2)/e6+y3]+t 2、大气平均作用温度Ta的近似估计 温度换算:T=t+273.15 本图为7月份拍摄,对于中纬度夏季平均大气Ta=16.0110+0.92621T0 取乌鲁木齐市平均气温为25摄氏度时Ta = 312.15753 3、大气透射率t6的估计 t6=0.974290-0.08007w,0.4≤w≤1.6。 w为水分含量,单位(g/cm2),这里,取w=1.0,计算得到t6=0.89422 4、地表比辐射率的估计 典型地表类型的比辐射率 ew=0.995 ev=0.986 em=0.970 Pv=[(NDVI- NDVIs)/(NDVIv- NDVIs)]2 NDVI 为归一化植被指数, 取NDVIv=0.70 和NDVIs=0(分别取自5%及95%数据) e surface=0.9625+0.0614Pv-0.0461Pv^2 5、像元亮度温度计算 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.303 LMIN_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.11471 y2=-1.1836w^2- 0.37607w- 0.52894=0.27859 y3=-0.04554w^2+1.8719w- 0.39071=1.43565 t=T6-r*L6 其中, C1和C2是Plank 函数的常量, C1=1.19104×10^8 W·μm4·m- 2·sr- 1, C2=14387.7 μm·K; T6是卫星高 度上传感器所探测到的像元亮度温度, 单位为K; λ

定量遥感_地表温度反演

遥感数字影像处理 作品名称:黄河三角洲地表温度反演 +学号: 小组成绩:

一、概述 1、作业背景: 地表温度是很多环境模型的一个重要参数,在大气与地表的能量与物质交换,天气预报,全球洋流循环,气候变化等研究领域有重要的应用。利用热红外遥感可以得到大围的地表温度面状信息,与传统的地表温度测量方式相比,具有快速、便捷、测量围大、信息连续等特点,因此利用热红外遥感数据反演地面温度得到了广泛的应用 2、作业意义: 黄河三角洲是黄河携带大量泥沙在渤海凹陷处沉积形成的冲积平原,位处黄河入处的黄河三角洲自然保护区正是以保护河口湿地生态系统和珍稀、濒危鸟类为主的湿地类型保护区。以利津为顶点,北到徒骇河口,南到小清河口,呈扇状三角形,面积5,450平方公里。地面平坦,在海拔10公尺以下。向东撒开的扇状地形,海拔高程低于15米,面积达5450平方公里。三角洲属,温带季风性气候。四季分明,光照充足,区自然资源丰富。 黄河口湿地生态旅游区占地23万亩,都处在黄河三角洲之,地貌以芦苇沼泽,湿地为主,其次为河口滩地,带翅碱蓬盐滩湿地,灌丛疏林湿地以及人工槐林湿地等。集自然景观与人文景观为一体,既有沧海桑田的神奇与壮阔,又有黄龙入海的壮观和长河落日的静美,是人们休闲、度假、观光科普的最佳场所。 二、数据介绍 数据来自地理空间数据云,Landsat 4-5 TM(陆地卫星4、5号,1982年发射后运行至今,携带有TM传感器)的相关遥感影像作为研究数据,研究黄河三角洲温度分布状况。 实验数据:2010年9月11号黄河三角洲图像(中心经度:118.8878w,中心纬度:37.4815n) 三、基本概念及技术流程图 3.1、基本概念:

landsat 遥感影像地表温度反演教程

基于辐射传输方程的Landsat 数据地表温度反演教程 一、数据准备 Landsa 8遥感影像数据一景,本教程以重庆市2015年7月26日的=行列号为(128,049)影像(LC8LGN00)为例。 同时需提前查询影像的基本信息(详见下表) 二、地表温度反演的总体流程 三、具体步骤 1、辐射定标 地表温度反演主要包括两部分,一是对热红外数据,二是多光谱数据进行辐射定标。 (1)热红外数据辐射定标 选择Radiometric Correction/Radiometric Calibration 。在File Selection 对话框中,选择数据LC8LGN02_MTL_Thermal ,单击Spectral Subset 选择Thermal Infrared1(),打开Radiometric Calibration 面板。 (2)多光谱数据辐射定标 选择要校正的多光谱数据“进行辐射定标。 Settings ,如下图。 2、大气校正

本教程选择Flaash 校正法。FLAASH Atmospheric Correction,双击此工具,打开辐射定标的数据,进行相关的参数设置进行大气校正。 注意:如果在多光谱数据辐射定标时Scale factor值忘记设置,可在本步骤中打开辐射定标数时设置single scale faceor 值为,若已设置,则默认值为1即可。 1)Input Radiance Image:打开辐射定标结果数据; 2)设置输出反射率的路径,由于定标时候; 3)设置输出FLAASH校正文件的路径,最优状态:路径所在磁盘空间足够大; 4)中心点经纬度Scene Center Location:自动获取; 5)选择传感器类型:Landsat-8 OLI;其对应的传感器高度以及影像数据的分辨率自动读取; 6) 设置研究区域的地面高程数据; 7)影像生成时的飞行过境时间:在layer manager中的Lc8数据图层右键选择View Metadata,浏览time字段获取成像时间; 注:也可以从元文件“”中找到,具体名称:DATE_ACQUIRED = 2013-05-12;SCENE_CENTER_TIME = 02:55:; 8) 大气模型参数选择:Sub-Arctic Summer(根据成像时间和纬度信息选择); 9) 气溶胶模型Aerosol Model:Urban,气溶胶反演方法Aerosol Retrieval:2-band(K-T); 10) 其他参数按照默认设置即可。 11) 多光谱参数设置中, K-T反演选择默认模式:Defaults->Over-Land Retrieval standard (600:2100) 波谱响应函数:默认指向.. \Program Files\Exelis\ENVI51\classic\filt_func\ 把它重新指向:..\Program Files\Exelis\ENVI51\resource\filterfuncs\ 注:这是因为版本的一个小bug,即Classic中的L8的波谱响应函数不正确,另外一个一劳永逸的方法是:将

定量遥感:地表温度反演

作品名称:黄河三角洲地表温度反演 姓名+学号: 小组成绩:

一、概述 1、作业背景: 地表温度是很多环境模型的一个重要参数,在大气与地表的能量与物质交换,天气预报,全球洋流循环,气候变化等研究领域有重要的应用。利用热红外遥感可以得到大范围的地表温度面状信息,与传统的地表温度测量方式相比,具有快速、便捷、测量范围大、信息连续等特点,因此利用热红外遥感数据反演地面温度得到了广泛的应用 2、作业意义: 黄河三角洲是黄河携带大量泥沙在渤海凹陷处沉积形成的冲积平原,位处黄河入海口处的黄河三角洲自然保护区正是以保护河口湿地生态系统和珍稀、濒危鸟类为主的湿地类型保护区。以利津为顶点,北到徒骇河口,南到小清河口,呈扇状三角形,面积5,450平方公里。地面平坦,在海拔10公尺以下。向东撒开的扇状地形,海拔高程低于15米,面积达5450平方公里。三角洲属,温带季风性气候。四季分明,光照充足,区内自然资源丰富。 黄河口湿地生态旅游区占地23万亩,都处在黄河三角洲之内,地貌以芦苇沼泽,湿地为主,其次为河口滩地,带翅碱蓬盐滩湿地,灌丛疏林湿地以及人工槐林湿地等。集自然景观与人文景观为一体,既有沧海桑田的神奇与壮阔,又有黄龙入海的壮观和长河落日的静美,是人们休闲、度假、观光科普的最佳场所。 二、数据介绍 数据来自地理空间数据云,Landsat 4-5 TM(陆地卫星4、5号,1982年发射后运行至今,携带有TM传感器)的相关遥感影像作为研究数据,研究黄河三角洲温度分布状况。 实验数据:2010年9月11号黄河三角洲图像(中心经度:118.8878w,中心纬度:37.4815n) 三、基本概念及技术流程图 3.1、基本概念:

基于单窗算法反演地表温度的ENVI操作教程(word文档良心出品)

单窗算法反演地表温度教程 1.1算法原理 1.1.1单窗算法 单窗算法(MW算法)是覃志豪于2001年提出的针对TM数据只有一个热红外波段的地面温度反演算法。经过众多学者验证,单窗算法具有很高的反演精度,且同样适用于ETM+和Iandsat 8数据。公式如下: T s (a6(1C6 D6)(b6(1C6 D6 )C6 D6)T s ensor D6T a)/C6 式中,LST为地表温度(K),T sensor是传感器上的亮度温度(K),T a 是大气平均温度(K); a、b为参考系数,当地表温度为0-70C时,a =-67.355351, b = 0.458606 C、D为中间变量,计算公式为: C = ez D = (l- r)[l-F (1- r)r] 式中,&为地表比辐射率,<为地面到传感器的大气总透射率。 因此单窗算法反演地表温度的关键是计算得到亮度温度T senso、地表比辐射率匚、大气透射率?和大气平均作用温度T a。 1.1.2参数计算1.1. 2.1辐射亮温计算 利用Planek公式将图像像元对应传感器辐射强度值转换为对应的亮度温度值。公式如下

£ hid + 41) 式中,T senso为亮度温度值;L影像预处理后得到的光谱辐射值, 2 单位为w/(m sr m),K1、K2为常量,可由数据头文件获取。 计算图像辐射亮温之前,需采用辐射定标参数将像元灰度值DN 转换为热辐射强度值,公式如下: 厶I二x DN十虫丄 式中,M L为增益参数,A L为偏移参数,该参数可直接在影像通文件数据中获取,且ENVI软件中已经集成,不需要自己在查找。 1.1. 2.2地表比辐射率计算 根据覃志豪针对TM影像提出的混合像元分解法来确定区域地表福辐射率。对于城市区域,我们简单的将其分为水体、自然表面和建筑表面三种,因此针对混合像元尺度上的地表比辐射率通过下式来估算: €= PR £+Q —F)/? e +/ r r v% r 7mm

地表温度完整算法

地表温度完整演算公式 一、DN值反演热辐射强度: L(λ)= 0.1238 + 0.005632156Qdn 二、热辐射强度反演地表亮温: T6=K2/ln(1 +K1/L(λ)) 其中,T6为TM6的象元亮度温度(K),K1和K2为发射前预设的常量,对于Landsat 5的TM数据,K1=60.776 mW cm-2sr-1μm-1,K2=1 260.56K。 三、地表亮温推算地表温度:(需确定三个参数:地表辐射率、大气平均温度、大气投射率) (1)地表辐射率推算:(主要根据NDVI推算地表辐射率) ①自然表面地表辐射率: ε=1.0094+0.047ln(NDVI) ②城镇用地地表辐射率: ε= 0.9589+0.086Pv-0.0671Pv2 Pv= [(NDVI- NDVI s)/(NDVI v- NDVI s)]2 其中,NDVI为归一化植被指数,取NDVIv=0.70和NDVIs=0.05,且有,当某个像元的NDVI大于0.70时,Pv取值为1;当NDVI小于0.05,Pv取值为0。 ③水体地表辐射率: 水体辐射率一般取值 0.995. (2)大气平均温度推算:

①热带平均大气( 北纬15°, 年平均) Ta= 17.9769+0.91715T0 ②中纬度夏季平均大气( 北纬45°, 7 月) Ta= 16.0110+0.92621T0 ③中纬度冬季平均大气( 北纬45°, 1 月) Ta= 19.2704+0.91118T0 根据万州所经纬度坐标:北纬38度48分,东经102度22分,其遥感影像两幅为2010年8月一幅为10月份,因此选取公式: Ta= 16.0110+0.92621T0 来推算大气平均温度。 其中:T0为距地表2米左右的温度,根据当时所处月份和卫星运行周期,我选取的T0 为 305K。 (3)大气透射率推算: 一般情况下,大气水分含量在0.4~3.0 g/cm2 变动区间,根据专家的研究结果大气透射率可以通过与水分含量建立线性函数关系来推算。 由于万州正处于8月到10月之间,重庆夏天较热,其温度应该在35度左右,所以根据当时实际情况,推测万州早上9-10点钟的

地表温度反演单窗算法

一:单窗算法 1、TM/ETM波段的热辐射传导方程: B6(T6)=t6(q)[e6B6(Ts)+(1-e6)I6~]+I6_ Ts是地表温度; T6是TM6的亮度温度; t6是大气透射率; e6是地表辐射率。 B6(T6)表示TM6遥感器所接收到的热辐射强度; B6(Ts)是地表在TM6波段区间内的实际热辐射强度,直接决取于地表温度; I6~和I6_分别是大气在TM6波段区间内的向上和向下热辐射强度。 2、化简后最终的单窗体算法模型为: Ts={a(1-C-D)+[b(1-C+D)+C+D]T6-DTa}/C 式中 C6=t6e6(e6为比辐射率,t6为透射率) D6=(1-t6)[1+t6(1-e6)] a =-67.355351,b=0.458606 3、大气平均作用温度Ta的近似估计 温度换算:T=t+273.15 本图为7月份拍摄,对于中纬度夏季平均大气Ta=16.0110+0.92621T0 取乌鲁木齐市平均气温为25摄氏度时Ta = 312.15753 4、大气透射率t6的估计 t6=0.974290-0.08007w,0.4≤w≤1.6。 w为水分含量,单位(g/cm2),这里,取w=1.0,计算得到t6=0.89422 5、地表比辐射率的估计 典型地表类型的比辐射率 ew=0.995 ev=0.986 em=0.970 Pv=[(NDVI- NDVIs)/(NDVIv- NDVIs)]2 NDVI 为归一化植被指数, 取NDVIv=0.70 和NDVIs=0(分别取自5%及95%数据) e surface=0.9625+0.0614Pv-0.0461Pv^2 6、像元亮度温度计算 T6=K2/ln(1+K1/L6) 其中:K1=607.76,K2=1260.56(覃志豪,用陆地卫星TM6数据演算地表温度的单窗算法); L6为遥感器接收的辐射强度。 7、遥感器接收的辐射强度计算 L6=(15.303-1.238)*b1/255.0+1.238 (覃志豪,用陆地卫星TM6数据演算地表温度的单窗算法);b1为第六波段像元灰度值(DN值)。 LMAX_BAND6 = 15.303 LMIN_BAND6 = 1.238

遥感反演地表温度

1、 裁剪出出济南市区 2、 分别利用ENVI 、ERDAS 反演地表温度(LST )、NDVI ,对LST 进行彩色显示。 3、 分析LST 、NDVI 的关系。 反演公式 具体流程: 图像的DN 值 辐射亮度 辐射亮温 地表温度。 反演时从图像数值(DN )转换成绝对辐射亮度值时的公式、从辐射亮度值转成辐射亮温时的公式、从亮温转换成地表温度时的公式分别是: min min max 6255)(L L L DN L tm +-?=、 )1/ln(/12+=λL K K T 、 ε ρλl n )/(1T T T s += 其中:6tm L 为TM 传感器所接收到的辐射亮度(mW .cm -2s r-1.um -1),max L 、min L 分别是传感器所接收到的最大和最小的辐射强度,即对应于DN =255和DN =0时的最大和最小辐射强度。对于Landsat5的TM 6波段, 1K =60.77mW .cm -2s r-1.um -1,2K =1260.56K 。S T 为地表温度(K ) ;T 为辐射温度(K );λ为有效波谱范围内的最大灵敏度值,λ=11.5um ,ρ=/hc δ=1.438×10-2mk ,其中δ=1.38 ×10-23/J k ,为玻尔兹曼常数,h =6.626×10-34Js ,为Plank′s 常数,c =2.998 ×108/m s ,为光速。一般地,有植被覆盖的地表取ε=0.95,没有植被覆盖的地表取ε=0.92(Weng ,2004[16])。 min L =0.1238 255 )(min max L L - =0.005632156 )1/ln(/12+=λL K K T 1260.56 / LOG ( 1 + 60.766 / $n8_fu ) $n1_12736l / (1 + (0.0000115 * $n1_12736l /0.01438) * LOG (0.95 ) )

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

单窗算法反演地表温度教程 1.1 算法原理 1.1.1 单窗算法 单窗算法(MW 算法)是覃志豪于2001年提出的针对TM 数据只有一个热红外波段的地面温度反演算法。经过众多学者验证,单窗算法具有很高的反演精度,且同样适用于ETM+和landsat 8数据。公式如下: 6666666666/)))1(()1((C T D T D C D C b D C a T a sensor s -++--+--= 式中,LST 为地表温度(K ),T sensor 是传感器上的亮度温度(K ),T a 是大气平均温度(K );a 、b 为参考系数,当地表温度为0-70℃时,a = -67.355351,b = 0.458606;C 、D 为中间变量,计算公式为: 式中,为地表比辐射率,为地面到传感器的大气总透射率。 因此单窗算法反演地表温度的关键是计算得到亮度温度T senso 、地表比辐射率、大气透射率和大气平均作用温度T a 。 1.1.2 参数计算 1.1. 2.1 辐射亮温计算 利用Planck 公式将图像像元对应传感器辐射强度值转换为对应的亮度温度值。公式如下

式中,T senso 为亮度温度值;λL 影像预处理后得到的光谱辐射值, 单位为 )/(2m sr m w μ??,K1 、K2为常量,可由数据头文件获取。 计算图像辐射亮温之前,需采用辐射定标参数将像元灰度值DN 转换为热辐射强度值,公式如下: 式中,M L 为增益参数,A L 为偏移参数,该参数可直接在影像通文件数据中获取,且ENVI 软件中已经集成,不需要自己在查找。 1.1. 2.2 地表比辐射率计算 根据覃志豪针对TM 影像提出的混合像元分解法来确定区域地表福辐射率。对于城市区域,我们简单的将其分为水体、自然表面和建筑表面三种,因此针对混合像元尺度上的地表比辐射率通过下式来估算: 式中,为混合像元的地表比辐射率;P V 为植被覆盖率;R V 为植被的温度比率;R M 为建筑表面的温度比率; V 表示植被法地表比辐射率,m 表示建筑表面的地表比辐射率;d 表示辐射校正项。 根据覃志豪经验公式, V =0.986;m =0.972。

相关文档
相关文档 最新文档