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

基于GIS和遥感福州市植被生态环境质量综合分析

——技术步骤和具体试验步骤

一、技术步骤

具体试验步骤图1所表示。

图1.技术步骤图

二、试验步骤

（一）数据预处理

（1）波段叠加

1、11月5日Landsat5TM原始数据图2所表示。

图2.原始数据

2.在ENVI软件中点击菜单BasicTools，选择子菜单LayerStacking，图3所表示。

图3.波段叠加菜单（ENVI4.5）

3.在弹出窗口中设置参数：选择要进行叠加波段数据和结果数据名称和存放路径，图4所表示。生成多波段影像图5所表示。

图4.设置波段叠加参数图5.生成多波段影像（RGB为5，4，3组合）

（2）投影转换

利用Arcmap软件中投影转换工具Project，对福州市行政界线图层进行投影转换。图6所表示。

图6.投影转换

（3）影像裁剪

1.在ENVI软件中将福州市矢量边界图层和遥感影像叠加显示。首先打开遥感影像图，然后在视图中点击菜单Overlay下子菜单Vectors…，在弹出窗口中点击菜单File，选择OpenVectorFile…，图7所表示。

图7.叠加矢量图菜单

2.在弹出窗口中设置参数以下：选择输入矢量边界图层，选择输出.evf数据文件名及其存放路径，选择矢量图曾投影坐标系统，然后点击OK，在弹出窗口中设置矢量图层显示颜色为红色，点击Apply。图8所表示。遥感影像图和福州市矢量边界图层叠加显示图9所表示。

图8.设置要叠加矢量图参数

图9.遥感影像图和福州市矢量边界图层叠加显示

3.利用掩膜工具将福州市矢量边界图层作为掩膜，裁剪出福州市遥感影像图：点击菜单BasicTools，选择子菜单Masking中ApplyMask。在弹出窗口中选择要进行裁剪遥感影像图，并点击MaskOption按钮下拉箭头，选择BuildMask..，建立掩膜，图10所表示。

图10.建立掩膜菜单

4.在弹出来MaskDefinition窗口中点击Option按钮下拉箭头，选择菜单ImportEVFs…,导入矢量边界图层。图11所表示。

图11.导入矢量边界图层

5.点击OK，在弹出窗口中选择和矢量边界图层关联影像图，点击OK，回到MaskDefinition窗口中选择输出掩膜文件名及其存放路径。图12所表示。

图12.设置建立掩膜所需参数

6.点击OK,回到ApplyMaskInputFile窗口，点击SpatialSubset按钮，在弹出窗口中点击ROI/EVF按钮，选择矢量图层，这么裁剪空间范围就是矢量边界范围。图13所表示。

图13.选择裁剪空间范围

7.点击OK,能够看到裁剪空间范围值也对应改变了。点击OK，回到ApplyMaskInputFile窗口，点击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.在BandMath窗口中输入公式：((15.30-1.24)*b6)/255+1.24。点击OK，在弹出窗口中进行变量匹配。生成热度指数影像图图19所表示。

图19.热度指数波段运算和热度指数影像图

（3）湿度指数提取

在BandMath窗口中输入公式：0.1509*(b1)+0.1973*(b2)+0.3279*(b3)+0.3406*(b4)-0.7112*(b5)-0.4572*(b7)。

点击OK，在弹出窗口中进行变量匹配。生成湿度指数影像图图20所表示。

图20.湿度指数波段运算和湿度指数影像图

（3）土壤亮度指数提取

在BandMath窗口中输入公式：0.3037*