第七届全国大学生GIS应用技能大赛——生态敏感性分析

目录

1. 案例背景

2. 出图展示

3. 数据说明

4. 要求

5. 步骤

5.1 加载显示数据

5.2 DEM数据预处理

5.2.1 镶嵌至新栅格

5.2.2 按掩膜提取

5.3 提取坡度以及重分类

5.3.1 提取坡度

5.3.2 坡度重分类

5.4 高程重分类

5.5 提取坡向以及重分类

5.5.1 提取坡向

5.5.2 坡向重分类

5.6 计算NDVI以及重分类

5.6.1 波段合成

5.6.2 计算NDVI

5.6.3 NDVI重分类

5.7 提取水体以及重分类

5.7.1 计算NDWI

5.7.2 提取水体

5.7.3 计算欧氏距离

5.7.4 欧氏距离重分类

5.8 计算加权总和以及重分类

5.8.1 计算加权总和

5.8.2 加权总和重分类

6. 整饰出图


1. 案例背景

生态环境敏感性是指生态系统对区域内自然和人类活动干扰的敏感程度,它反映区域生态系统在遇到干扰时,发生生态环境问题的难易程度和可能性的大小,并用来表征外界干扰可能造成的后果。即在同样干扰强度或外力作用下,各类生态系统出现区域生态环境问题可能性的大小,研究生态敏感性对保护生态系统及其重要。按照下列要求,进行生态敏感性评价。

2. 出图展示

第七届全国大学生GIS应用技能大赛——生态敏感性分析

3. 数据说明

①DEM数据文件夹:存放所需的DEM数据。

②遥感影像数据文件夹:存放所需的影像数据。

③研究区范围.shp:本次研究的研究范围。

4. 要求

①地形、植被、水体方面的生态因子及其对该地区的敏感性等级见表,请根据表中各因子权重值,加权计算植被覆盖区域的生态敏感性信息,生成该地区的生态敏感性等级分布专题图。(50分)

生态因子及其影响范围所赋属性值

生态因子

二级因子

分类

敏感性等级

赋值

权重

地形因子

坡度

(单位:度)

>60

极高敏感

5

0.2

45-60

高敏感

4

25-45

中敏感

3

10-25

低敏感

2

0-10

非敏感

1

高程

(单位:米)

>1300

极高敏感

5

0.1

1000-1300

高敏感

4

800-1000

中敏感

3

500-800

低敏感

2

<500

非敏感

1

坡向

正北

极高敏感

5

0.1

东北、西北

高敏感

4

正东、正西

中敏感

3

东南、西南

低敏感

2

平地、正南

非敏感

1

用地类型

植被

NDVI>=0.5

极高敏感

5

0.3

NDVI<=0

高敏感

4

0.3<=NDVI<0.5

低敏感

2

0<NDVI<0.3

非敏感

1

水系

水域区域

极高敏感

5

0.3

水域周边300m范围内

高敏感

4

水域周边300—800m范围内

中敏感

3

其他区域

低敏感

2

②对生态敏感性结果进行制图,分为5个等级(敏感性1-5),地图添加文字信息标注每个行政区名称,添加必要的地图元素并排版出图(将结果保存为jpg格式)。(10分)

5. 步骤

5.1 加载显示数据

首先,把整个工程的处理范围设置为与研究区图层的范围相同,如下图所示。

设置环境的处理范围与研究区范围相同

把研究区范围数据导入到GeoScene Pro软件中,并查看其数据信息,如下图所示。

查看研究区范围数据空间参考信息

接着导入DEM数据,并查看其栅格信息,如下图所示。

加载显示DEM栅格数据

查看DEM数据的栅格信息

5.2 DEM数据预处理

5.2.1 镶嵌至新栅格

首先要把以上的四个DEM数据镶嵌成一个数据,这里使用到的是【镶嵌至新栅格】工具,输入栅格设置为原始的四个DEM数据,在栅格数据的空间参考处导入原始DEM数据的空间参考坐标系,像素类型选择16位有符号,像元大小设置为30,波段数设置为1,这些信息都可以在原始的DEM数据属性中查看到,需要与原始DEM数据保持一致,不能想当然设置,如下图所示。

镶嵌至新栅格参数设置

镶嵌结果

然后把镶嵌好的DEM数据和研究区范围矢量数据叠加显示,如下图所示。

同时显示DEM数据和研究区范围数据

5.2.2 按掩膜提取

接下来需要裁剪出研究区范围内的DEM数据,使用到的是【按掩膜提取】工具,输入栅格设置为镶嵌好的DEM数据,输入栅格数据或要素掩膜数据设置为研究区范围的矢量数据,如下图所示。

按掩膜提取研究区DEM数据

研究区DEM数据提取结果

至此,原始的DEM数据预处理就算结束了,接下来进行其他的数据处理。

5.3 提取坡度以及重分类

5.3.1 提取坡度

由上述处理好的研究区DEM数据就可以提取坡度、坡向等数据了。提取坡度用到的是【Slope】工具,输入栅格设置为处理好的研究区DEM数据,其他参数默认即可,如下图所示。

坡度提取

坡度提取结果

5.3.2 坡度重分类

然后根据题目要求对坡度数据重分类,使用到的是【重分类】工具,输入栅格设置为上一步得到的坡度数据,重分类字段选择VALUE字段,按照上述要求中的生态因子及其影响范围所赋属性值表格的内容对坡度数据进行重分类,如下图所示。

坡度重分类

坡度重分类结果

5.4 高程重分类

接下来对高程数据进行重分类,分类依据也是参照题目要求中的生态因子及其影响范围所赋属性值表格,输入栅格设置为研究区的DEM数据,重分类字段选择VALUE字段,如下图所示。

研究区DEM重分类

研究区DEM重分类结果

5.5 提取坡向以及重分类

5.5.1 提取坡向

提取坡向使用的数据是处理好的研究区DEM数据,使用到的工具是【坡向】工具,输入栅格设置为研究区DEM数据,如下图所示。

坡向提取

坡向提取结果

5.5.2 坡向重分类

然后根据题目要求中的生态因子及其影响范围所赋属性值表格对坡向数据进行重分类,使用到的工具是【重分类】工具,输入栅格设置为坡向数据,重分类字段选择VALUE字段,如下图所示。

坡向重分类

坡向重分类结果

5.6 计算NDVI以及重分类

5.6.1 波段合成

首先把原始数据提供的Landsat5数据的1至7波段加载进来,如下图所示。

添加Landsat5数据的1至7波段

Landsat5数据的1至7波段添加结果

然后使用【波段合成】工具进行波段合成操作,输入栅格设置为上一步加载进来的1至7波段遥感影像数据,如下图所示。

遥感影像波段合成

遥感影像波段合成结果

5.6.2 计算NDVI

计算NDVI使用到的是指数菜单中的NDVI工具,近红外波段指数设置为第4波段,红色波段指数设置为第3波段,如下图所示。

NDVI指数工具

NDVI计算的波段参数设置

NDVI计算结果

然后需要提取出研究区范围内的NDVI数据,使用到的是【按掩膜提取】工具,输入栅格设置为上一步得到的NDVI数据,输入栅格数据或要素掩膜数据设置为研究区范围矢量数据,如下图所示。

提取研究区NDVI数据

研究区NDVI提取结果

5.6.3 NDVI重分类

接下来需要按照题目要求中的生态因子及其影响范围所赋属性值表格对研究区的NDVI进行重分类,使用到的是【重分类】工具,输入栅格设置为上一步得到的研究区NDVI数据,重分类字段选择VALUE字段,如下图所示。

研究区NDVI重分类

研究区NDVI重分类结果

5.7 提取水体以及重分类

5.7.1 计算NDWI

计算NDWI使用到的工具是【栅格计算器】工具,在栅格计算器中输入以下公式:(Float("L5119039_03920100524_B20.TIF")- Float("L5119039_03920100524_B40.TIF")) / (Float("L5119039_03920100524_B20.TIF")+ Float("L5119039_03920100524_B40.TIF")),这里所说的输入并非逐个字符地输入,而是通过双击数据来写入需要使用到的数据名称,然后通过基本的运算符连接起来形成一个公式,如下图所示。

计算总范围的NDWI数据

总范围的NDWI数据计算结果

然后需要提取出研究区范围内的NDWI数据,使用到的工具是【按掩膜提取】工具,输入栅格设置为上一步得到的总范围的NDWI数据,输入栅格数据或要素掩膜数据设置为研究区范围矢量数据,如下图所示。

提取研究区NDWI数据

研究区NDWI提取结果

5.7.2 提取水体

根据NDWI指数大于等于0的区域理论上是水体来提取水体数据,使用到的是【栅格计算器】工具,在栅格计算器中输入以下公式:Con( "研究区NDWI">= 0,1),这里所谓的输入也是通过双击研究区NDWI数据来写入数据名称的,并非逐个字符地手动输入,如下图所示。

提取水体

水体提取结果

5.7.3 计算欧氏距离

接下来计算水体的欧氏距离,这里所谓的欧氏距离可以理解为距离水体的远近。计算欧氏距离使用到的是【欧氏距离】工具,输入栅格或要素源数据设置为上一步得到的水体数据,输出像元大小设置为30,这里的像元大小可以通过查看本题中其他的栅格数据的属性信息来得知,需要与其他的栅格数据的像元大小保持一致,其它参数默认即可,如下图所示。

计算欧氏距离

欧氏距离计算结果

然后需要提取出研究区范围内的欧氏距离数据,使用到的工具是【按掩膜提取】工具,输入栅格设置为上一步得到的欧氏距离数据,输入栅格数据或要素掩膜数据设置为研究区范围矢量数据,如下图所示。

提取研究区欧氏距离

研究区欧氏距离提取结果

5.7.4 欧氏距离重分类

接下来需要按照题目要求中的生态因子及其影响范围所赋属性值表格来对欧氏距离进行重分类,使用到的是【重分类】工具,输入栅格设置为上一步得到的研究区欧氏距离数据,重分类字段选择VALUE字段,如下图所示。

研究区欧氏距离重分类

研究区欧氏距离重分类结果

5.8 计算加权总和以及重分类

5.8.1 计算加权总和

最后一步是根据上述得到的各生态因子数据来计算它们对生态环境敏感性的综合影响,也要参照题目要求中的生态因子及其影响范围所赋属性值表格来进行该操作,其中坡度的权重是0.2,高程的权重是0.1,坡向的权重是0.1,NDVI的权重是0.3,欧氏距离的权重是0.3,这里利用【栅格计算器】工具来计算它们的加权结果,如下图所示。

计算各种影响因素的加权总和

加权总和计算结果

5.8.2 加权总和重分类

然后再根据题目要求对加权总和的结果数据进行重分类,需要分为5类,使用到的是【重分类】工具,输入栅格设置为加权总和结果数据,重分类字段设置为VALUE字段,如下图所示。

加权总和重分类

加权总和重分类结果

最后设置一下加权总和重分类结果数据的符号系统就可以出图了,主符号系统选择唯一值,字段选择Value字段,然后添加所有值,在标注列分别把1至5等级设置为不敏感区域,轻度敏感区域,中度敏感区域,高度敏感区域,极其敏感区域,然后再调整一下各等级的填充颜色,如下图所示。

设置加权总和重分类数据的符号系统

6. 整饰出图

第七届全国大学生GIS应用技能大赛——生态敏感性分析

  • 9
    点赞
  • 52
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值