资源环境遥感实习实习一植被覆盖度计算实习目的•根据LandsatTM数据计算植被覆盖度原理与方法•NDVI为像元NDVI值,fv为像元的植被覆盖度,NDVIv和NDVI0分别为植被覆盖部分和非植被覆盖部分的NDVI值。00NDVINDVINDVINDVIfvv原理与方法•NDVIv和NDVI0这两个参数值的确定是关键,将直接影响到植被覆盖度估算结果。•对NDVI统计直方图给定置信区间,取5%和95%频率的NDVI值作为NDVI0和NDVIv值。实习数据•TM_20060504–南京地区2006年5月4日的LandsatTM数据•TM_20060504.hdr–TM_20060504的头文件计算NDVI•主菜单TransformNDVI,InputFileType选择LandsatTM,计算得到NDVI空间分布图。确定NDVIv和NDVI0值•在AvailabelBandsList窗口中右键点击NDVI数据,在弹出的右键菜单中选择QuickStats,对NDVI数据进行统计分析,在统计结果窗口中可以看到NDVI的直方图和累计直方图。•在累计直方图中查找5%和95%频率的NDVI值,分别约为-0.194和0.425。•因此NDVI0=-0.194,NDVIv=0.425。计算植被覆盖度•主菜单BasicToolsBandMath,在波段运算窗口输入植被覆盖度计算公式:(b1+0.194)/(0.425+0.194),计算得到初步结果。计算植被覆盖度•因为有5%像元的NDVI值小于-0.194,这部分像元的植被覆盖度计算结果为负值,需要通过掩膜操作将这部分像元的植被覆盖度值修改为0,同样需要将5%NDVI值大于0.425的像元植被覆盖度值修改为1。计算植被覆盖度•主菜单BasicToolsMaskingBuildMask,在掩膜定义窗口中设置掩膜生成条件为NDVI-0.194,生成掩膜。•ApplyMask,将掩膜应用于计算得到的植被覆盖度数据,注意在ApplyMaskParameters窗口中将MaskValue设为0,这样可将所有掩掉的像元值修改为0。计算植被覆盖度•类似操作将NDVI值大于0.425的像元植被覆盖度值修改为1,注意的是ApplyMaskParameters窗口中将MaskValue需要设为1。计算植被覆盖度•最终得到植被覆盖度空间分布图实习二城市热环境分析实习目的•基于LandsatTM数据反演地表温度,在此基础上对城市热环境进行分析。原理与方法•本实习采用覃志豪的单窗算法反演地表温度,单窗算法的具体形式为:s6aT[a(1CD)(b(1CD)CD)TDT]/C式中,Ts为地表温度,a和b为经验系数(a=-67.35535,b=0.458608),T6为TM6的亮度温度(K),Ta为大气平均作用温度(K),C和D为中间参数原理与方法•中间参数C和D计算公式如下:CD(1)[1(1)]式中,τ为大气透过率,ε为地表比辐射率原理与方法•比辐射率计算公式0.9230.1571.00940.047ln0.1570.7270.9940.727NDVINDVINDVINDVI0.995水体原理与方法•大气透过率计算公式•大气平均作用温度计算公式1.0314120.11536w019.27040.91118aTT式中,w为水汽含量式中,Ta和T0单位均为K原理与方法•LandsatTM亮度温度的计算公式K2T6ln(1K1/L)式中,K1、K2为系数(K1=607.76,K2=1260.56),L为辐射亮度,计算公式如下L=0.055158*DN+1.2378式中,DN值为TM6波段的灰度值实习数据•njtm:南京地区TM数据•njtm.hdr:njtm的头文件•卫星过境时南京地区水汽含量为2.49g/cm2,气温21℃辐射定标与亮温计算•主菜单BasicToolsBandMath,在波段运算窗口输入辐射定标公式0.055158*b1+1.2378,得到辐亮度数据。•然后输入公式:1260.56/alog(1+607.76/b1),计算得到亮度温度。地表比辐射率计算•计算出NDVI,然后在波段运算窗口输入公式:1.0094+0.047*alog(b1)•分别运用掩膜操作将NDVI小于0.157以及大于0.727的像元的地表比辐射率修改为0.923和0.994,具体操作步骤可参考试验一中的类似部分。地表比辐射率计算•取NDVI阈值为-0.07生成水体掩膜,该掩膜保留所有非水体像元,即NDVI值大于等于-0.07的像元。•应用该掩膜,将所有水体像元的比辐射率值修改为0.995。大气透过率和大气平均作用温度计算•将卫星过境时南京地区水汽含量及气温代入计算方程,计算得到大气透过率和大气平均作用温度分别为0.7442和287.294K地表温度反演•首先根据大气透过率、地表比辐射率等通过波段运算计算出中间参数C、D,然后将单窗算法公式输入波段运算工具,分别代入各个参数,计算得到南京地区的地表温度。地表温度反演热环境分析•Image窗口菜单ToolsColormappingDensitySlice,运用密度分割将地表温度分为若干个不同的热岛等级,使其从冷色到暖色调变化,最终得到总共包含若干个热岛强度级别的南京市热岛强度空间分布图热环境分析
本文标题:资源环境遥感
链接地址:https://www.777doc.com/doc-4117006 .html