Landsat温度反演.docx

定量遥感实验2 基于Landsat TM数据和大气校正法的地表温度反演一实验目的学习并掌握基于Landsat 5TM遥感影像，运用辐射传输方程法（大气校正法）对地表温度进行反演的方法。二研究方法实验数据：2005年9月8日长春市Landsat 5TM 图像（景号：118-29和118-30）。实验流程：完整流程涉及Landsat TM 的数据读取、辐射定标、反演模型、遥感反演过程、反演结果验证等内容。软件功能：主要采用ENVI 主模块中的Landsat TM 数据读取、辐射定标、波段运算、结果统计分析等。该方法需要进行大气校正，消除大气层对地表辐射能量的影响，这就需要从卫星观测得到的热辐射能量中扣除大气层的辐射分量，并利用热红外波段（Band 6）范围内的地表发射率作为参数，反演出地表的真实温度。技术路线为：（1）LandsatTM数据预处理：数据读取、辐射定标。（2）相关辅助数据的确定与查找：大气上行辐射以及下行辐射，采用数据当天的大气透过率信息等；（3）采用大气校正法（辐射传输方程法）利用Landsat TM band 6 进行地表温度反演：首先获取地表比辐射率值，其次，计算黑体在热红外波段的辐射亮度；最后，利用普朗克公式的反函数反演出整个研究区域的地表温度分布情况。图1 流程图三实验步骤1、Landsat TM 数据读取和定标（1）数据读取在ENVI 主菜单中，选择File Open External File LandsatGeo TIFF with Metadata ，打开L5118029_02920050908_MTL.txt。其中包含两种数据：可见光波段数据（HRF）和热红外波段数据（HTM ）图2 Landsat TM 数据读取列表（2）数据辐射定标主菜单Basic ToolsPreprocessingCalibration UtilitiesLandsat Calibration ，选择需要定标的数据，包括可见光波段数据和热红外波段数据，该过程功能可以从读取文件中直接获取实验数据的元数据（成像时间，定标参数等等）。图3 辐射定标2、裁剪工程区（1）镶嵌影像主菜单MapMosaickingGeoreferenced，在Map Based Mosaic窗口中Import File，输入两副已经辐射定标的热红外波段遥感影像。影像导入完成后，可修改影像的投影信息，以便提高影像精度。镶嵌模块中，主菜单上选择OptionsChange Base Projection。打开转变投影功能后，根据影像信息设置投影。选择影像，单击右键选择Edit Entry，可以通过此功能，改变影像镶嵌参数。完成点击应用FileApply即可进入镶嵌过程，选择重采样方法等。图4遥感影像镶嵌（2）裁剪工程区①打开要裁剪的区域范围的矢量数据：主菜单FileOpen Vector选择文件类型Shapefile(*.shp)changchun_Disssolve.shp，注意在弹出的Import Vector Files Parameters 对话框选择正确的投影类型（与遥感影像的投影坐标一致）。点击OK之后ENVI自动将矢量文件转为EVF格式。图5 矢量数据读取②将矢量数据转为ROI：在Availabel Vectors List选择数据，通过 FileExport Layers to ROI，在Select Data File to Associate with new ROIS 中选择镶嵌好的热红外影像，再在Export EVF Layers toROI中选择Convert all records of an EVF layer to one ROI，点击OK，从而实现矢量数据转为ROI。图6 Export EVF Layers to ROI对话框③裁剪栅格数据：在ENVI主菜单Basic Tools Subset Data via ROIs，在Select Input File to Subset via ROI 中选择要裁剪的热红外波段影像（镶嵌好的），OK。出现Spatial Subset via ROI Parameters对话框，在Select Input ROIs 中选择建立的 ROI，选项Mask pixels outside of ROI ?选择Yes，在 Mask Background Value默认为0，这样得到长春市的热红外波段遥感影像。图7 Spatial Subset via ROI Parameters对话框对于可见光波段数据的裁剪方法与上述方法相同，由于可见光波段数据和热红外波段数据的分辨率不同，在利用该矢量数据裁剪可见光波段数据的时候需要重复②步骤，完成裁剪操作即可。3、可见光和近红外波段大气校正（快速校正）主菜单