地表温度反演算法.doc

地表温度作为地球环境分析的重要指标,而遥感技术作为现代重要的对地观测手段,使得基于遥感图像的地表温度反演的研究越来越多。主要的地表温度反演方法有:大气校正法,单窗算法,单通道法等等。本文介绍用辐射传输方程法对地表温度进行反演。技术流程:例子数据为2002年9月2日的襄樊市LandsatETM+数据。根据数据的特点以及地表温度反演研究的技术要求,采用的技术路线为:先对LandsatETM+数据进行预处理:数据读取、辐射定标、大气校正、襄樊区域裁剪,利用大气校正,即:辐射传输方程法对其影像热红外波段数据进行操作反演,实现襄樊市地区的地表真实温度的反演研究。具体的处理流程如下:具体的实现步骤如下:第一步:准备数据热红外数据使用的是Landsat的第六波段,已经做了传感器定标、几何校正、工程区裁剪,详细流程参考上面的流程图。文件为TM6-rad-subset-jz-。由TM影像(已经过大气校正)生成的NDVI数据,已经利用主菜单->BasicTools->ResizeData(SFatial/SFectral)重采样为60米分辨率,与TMi6数据保持一致,文件名为:TM-NDVI-。第二步:地表比辐射率计算物体的比辐射率是物体向外辐射电磁波的能力表征。它不仅依赖于地表物体的组成,而且与物体的表面状态(表面粗糙度等)及物理性质(介电常数、含水量等)有关,并随着所测定的波长和观测角度等因素有关。在大尺度上对比辐射率精确测量的难度很大,目前只是基于某些假设获得比辐射率的相对值,本文主要根据可见光和近红外光谱信息来估计比辐射率。(一)植被覆盖度计算计算植被覆盖度Fv采用的是混合像元分解法,将整景影像的地类大致分为水体、植被和建筑,具体的计算公式如下:FV=(NDVI-NDVIS)/(NDVIV-NDVIS)(2)其中,NDVI为归一化差异植被指数,取NDVIV==,且有,,FV取值为1;,FV取值为0。利用ENVI主菜单->BasicTools->BandMath,在公式输入栏中输入:()*1+(b1lt0.)*0+()*((b1-)/(-))b1:选择NDVI图像得到植被盖度图像。(二)地表比辐射率计算根据前人的研究,将遥感影像分为水体、城镇和自然表面3种类型。本专题采取以下方法计算研究区地表比辐射率:,自然表面和城镇像元的比辐射率估算则分别根据下式(3)(4)进行计算:εsurface=+-(3)εbuilding=+-(4)式中,εsurface和εbuilding分别代表自然表面像元和城镇像元的比辐射率。利用ENVI主菜单->BasicTools->BandMath,在公式输入栏中输入:(b1le0)*+()*(+*b2-*b2^2)+()*(+*b2-*b2^2)b1:NDVI值;b2:植被覆盖度值。得到地表比辐射率数据。第三步:计