基于gis与遥感的福州市植被生态环境质量综合分析技术流程与具体实验步骤.docVIP

PAGE PAGE 1 基于GIS与遥感的福州市植被生态环境质量综合分析 ——技术流程与具体实验步骤 一、技术流程 具体的实验流程如图1所示。 图1.技术流程图 二、实验步骤 （一）数据预处理 （1）波段叠加 1、2006年11月5日的Landsat5 TM原始数据如图2所示。 图2.原始数据 2.在ENVI软件中点击菜单Basic Tools，选择子菜单Layer Stacking，如图3所示。 图3.波段叠加菜单（ENVI4.5） 3.在弹出的窗口中设置参数：选择要进行叠加的波段数据以及结果数据名称和存放路径，如图4所示。生成的多波段影像如图5所示。 图4.设置波段叠加的参数 图5.生成的多波段影像（RGB为5，4，3组合） （2）投影转换 利用Arcmap软件中的投影转换工具Project，对福州市行政界线图层进行投影转换。如图6所示。 图6.投影转换 （3）影像裁剪 1.在ENVI软件中将福州市的矢量边界图层与遥感影像叠加显示。首先打开遥感影像图，然后在视图中点击菜单Overlay下的子菜单Vectors…，在弹出的窗口中点击菜单File，选择Open Vector File…，如图7所示。 图7.叠加矢量图的菜单 2.在弹出的窗口中设置参数如下：选择输入的矢量边界图层，选择输出的.evf数据文件名及其存放路径，选择矢量图曾的投影坐标系统，然后点击OK，在弹出的窗口中设置矢量图层的显示颜色为红色，点击Apply。如图8所示。遥感影像图与福州市矢量边界图层叠加显示如图9所示。 图8.设置要叠加矢量图的参数 图9. 遥感影像图与福州市矢量边界图层叠加显示 3.利用掩膜工具将福州市矢量边界图层作为掩膜，裁剪出福州市的遥感影像图：点击菜单Basic Tools，选择子菜单Masking中的Apply Mask。在弹出的窗口中选择要进行裁剪的遥感影像图，并点击Mask Option按钮的下拉箭头，选择Build Mask..，建立掩膜，如图10所示。 图10. 建立掩膜菜单 4.在弹出来的Mask Definition窗口中点击Option按钮的下拉箭头，选择菜单Import EVFs…,导入矢量边界图层。如图11所示。 图11.导入矢量边界图层 5.点击OK，在弹出的窗口中选择与矢量边界图层关联的影像图，点击OK，回到Mask Definition窗口中选择输出掩膜的文件名及其存储路径。如图12所示。 图12.设置建立掩膜所需的参数 6.点击OK,回到Apply Mask Input File窗口，点击Spatial Subset按钮，在弹出的窗口中点击ROI/EVF按钮，选择矢量图层，这样裁剪的空间范围就是矢量边界的范围。如图13所示。 图13.选择裁剪的空间范围 7.点击OK,可以看到裁剪的空间范围值也相应改变了。点击OK，回到Apply Mask Input File窗口，点击OK。 图14. 空间范围值的改变 8.在弹出的窗口中设置掩膜的值为255，这样生成的影像的背景色是白色的。选择输出的文件名及其存放路径。裁剪出的福州市遥感影像图如图15所示。 图15. 福州市遥感影像图（RGB543） （二）提取遥感信息 （1）植被指数提取 1. 植被指数的计算公式为 NDVI =(NIR-R)/(NIR+R)。打开福州市遥感影像，利用波段运算工具，输入公式：float(b4)-float(b3))/(float(b4)+float(b3))。 图16.植被指数的波段运算 2.点击OK，在弹出的窗口中进行变量与波段之间的配对。然后选择输出文件名及其存放路径。结果如图17所示。 图17. 变量与波段之间的配对以及生成的植被指数影像图 （2）热度指数提取 1.热红外波段的遥感数据TM6 (10.4～12.5um) ,对温度敏感,经辐射定标后,直接得到辐射温度,可以根据地物辐射响应的差异,反映不同的地类,因而TM6图像数据可以作为温度的相对指标。故在本实验中,以辐射校正处理后的TM6作为热度指数。辐射校正的公式以及参数值如图18所示。 图18.辐射校正的公式以及参数值 2.在Band Math窗口中输入公式：((15.30-1.24)*b6)/255+1.24。点击OK，在弹出的窗口中进行变量匹配。生成的热度指数影像图如图19所示。 图19.热度指数的波段运算与热度指数影像图 （3）湿度指数提取 在Band Math窗口中输入公式：0.1509*(b1)+0.1973*(b2)+0.3279*(b3)+0.3406*(b4)-0.7112*(b5)-0.4572*(b7)。 点击OK，在弹出的窗口中进