资源环境遥感课件

资源环境遥感实习资源环境遥感实习1实习一植被覆盖度计算实习一植被覆盖度计算2实习目的根据LandsatTM数据计算植被覆盖度实习目的根据LandsatTM数据计算植被覆盖度3原理与方法NDVI为像元NDVI值，fv为像元的植被覆盖度，NDVIv和NDVI0分别为植被覆盖部分和非植被覆盖部分的NDVI值。原理与方法4原理与方法NDVIv和NDVI0这两个参数值的确定是关键，将直接影响到植被覆盖度估算结果。对NDVI统计直方图给定置信区间，取5%和95%频率的NDVI值作为NDVI0和NDVIv值。原理与方法NDVIv和NDVI0这两个参数值的确定是关键，将5实习数据TM_20060504南京地区2006年5月4日的LandsatTM数据TM_20060504.hdrTM_20060504的头文件实习数据TM_200605046计算NDVI主菜单Transform>NDVI，InputFileType选择LandsatTM，计算得到NDVI空间分布图。计算NDVI主菜单Transform>NDVI，Inpu7确定NDVIv和NDVI0值在AvailabelBandsList窗口中右键点击NDVI数据，在弹出的右键菜单中选择QuickStats，对NDVI数据进行统计分析，在统计结果窗口中可以看到NDVI的直方图和累计直方图。在累计直方图中查找5%和95%频率的NDVI值，分别约为-0.194和0.425。因此NDVI0=-0.194，NDVIv=0.425。确定NDVIv和NDVI0值在AvailabelBands8计算植被覆盖度主菜单BasicTools>BandMath，在波段运算窗口输入植被覆盖度计算公式：(b1+0.194)/(0.425+0.194)，计算得到初步结果。计算植被覆盖度主菜单BasicTools>BandM9计算植被覆盖度因为有5%像元的NDVI值小于-0.194，这部分像元的植被覆盖度计算结果为负值，需要通过掩膜操作将这部分像元的植被覆盖度值修改为0，同样需要将5%NDVI值大于0.425的像元植被覆盖度值修改为1。计算植被覆盖度因为有5%像元的NDVI值小于-0.194，这10计算植被覆盖度主菜单BasicTools>Masking>BuildMask，在掩膜定义窗口中设置掩膜生成条件为NDVI>-0.194，生成掩膜。ApplyMask，将掩膜应用于计算得到的植被覆盖度数据，注意在ApplyMaskParameters窗口中将MaskValue设为0，这样可将所有掩掉的像元值修改为0。计算植被覆盖度主菜单BasicTools>Maskin11计算植被覆盖度类似操作将NDVI值大于0.425的像元植被覆盖度值修改为1，注意的是ApplyMaskParameters窗口中将MaskValue需要设为1。计算植被覆盖度类似操作将NDVI值大于0.425的像元植被覆12计算植被覆盖度最终得到植被覆盖度空间分布图计算植被覆盖度最终得到植被覆盖度空间分布图13实习二城市热环境分析实习二城市热环境分析14实习目的基于LandsatTM数据反演地表温度，在此基础上对城市热环境进行分析。实习目的基于LandsatTM数据反演地表温度，在此基础上15原理与方法本实习采用覃志豪的单窗算法反演地表温度，单窗算法的具体形式为：式中，Ts为地表温度，a和b为经验系数（a=-67.35535，b=0.458608），T6为TM6的亮度温度(K)，Ta为大气平均作用温度(K)，C和D为中间参数原理与方法本实习采用覃志豪的单窗算法反演地表温度，单窗算法的16原理与方法中间参数C和D计算公式如下：式中，τ为大气透过率，ε为地表比辐射率原理与方法中间参数C和D计算公式如下：式中，τ为大气透过率，17原理与方法比辐射率计算公式原理与方法比辐射率计算公式18原理与方法大气透过率计算公式大气平均作用温度计算公式式中，w为水汽含量式中，Ta和T0单位均为K原理与方法大气透过率计算公式式中，w为水汽含量式中，Ta和19原理与方法LandsatTM亮度温度的计算公式式中，K1、K2为系数（K1=607.76，K2=1260.56），L为辐射亮度，计算公式如下L=0.055158*DN+1.2378式中，DN值为TM6波段的灰度值原理与方法LandsatTM亮度温度的计算公式式中，K120实习数据njtm：南京地区TM数据njtm.hdr ：njtm的头文件卫星过境时南京地区水汽含量为2.49g/cm2，气温21℃实习数据njtm：南京地区TM数据21辐射定标与亮温计算

主菜单BasicTools>BandMath，在波段运算窗口输入辐射定标公式0.055158*b1+1.2378，得到辐亮度数据。然后输入公式：1260.56/alog(1+607.76/b1)，计算得到亮度温度。辐射定标与亮温计算主菜单BasicTools>Ban22地表比辐射率计算

计算出NDVI，然后在波段运算窗口输入公式：1.0094+0.047*alog(b1)分别运用掩膜操作将NDVI小于0.157以及大于0.727的像元的地表比辐射率修改为0.923和0.994，具体操作步骤可参考试验一中的类似部分。地表比辐射率计算计算出NDVI，然后在波段运算窗口输入公式23地表比辐射率计算取NDVI阈值为-0.07生成水体掩膜，该掩膜保留所有非水体像元，即NDVI值大于等于-0.07的像元。应用该掩膜，将所有水体像元的比辐射率值修改为0.995。地表比辐射率计算取NDVI阈值为-0.07生成水体掩膜，该掩24大气透过率和大气平均作用温度计算

将卫星过境时南京地区水汽含量及气温代入计算方程，计算得到大气透过率和大气平均作用温度分别为0.7442和287.294K大气透过率和大气平均作用温度计算将卫星过境时南京地区水汽含25地表温度反演

首先根据大气透过率、地表比辐射率等通过波段运算计算出中间参数C、D，然后将单窗算法公式输入波段运算工具，分别代入各个参数，计算得到南京地区的地表温度。地表温度反演首先根据大气透过率、地表比辐射率等通过波段运算26地表温度反演地表温度反演27热环境分析

Image窗口菜单Tools>Colormappin