景观干扰度（arcgis和fragstats）

参考的up主后续视频说这个计算干扰度的方法是错的，麻烦大家看一下https://b23.tv/JiBtKZ8

软件：arcgis和fragstats

数据：分类后矢量数据

fragstats软件下载地址：https://www.umass.edu/landeco/research/fragstats/downloads/fragstats_downloads.html

安装即用

景观干扰度计算操作步骤来自

​​​​​​https://www.bilibili.com/video/BV1k54y1V7mS/spm_id_from=333.788.recommend_more_video.4

一、shp矢量数据转为tif栅格数据

1、根据要求重新分类：

选中矢量数据open attribute table，add field ”newclass”，select by attribute，field calculator

根据实际情况给 ”newclass”赋值

2、转换为投影坐标系：

ctrl+f 搜索 project(data management)

将坐标系转换到projected coordinate systems，

选择WGS 1984 UTM xxN

XX=经度/6+31

转自http://blog.sina.com.cn/s/blog_456959200102wrog.html

3、转为栅格数据

ctrl+f 搜索 feature to raster

field选择 newclass，out put cell size根据需要设置栅格文件像元值大小

如图输出的栅格数据一个像元是10000m2

4、创建栅格属性表

ctrl+f 搜索 build raster attribute table

input raster选择刚刚导出的栅格数据

5、导出栅格数据

export data，据说路径要全英文

二、fragstats计算景观指数

不要关arcmap，打开景观软件fragstats

1、导入数据

点击new，点击add layer

先选中geo tiff grid ，再选中刚刚导出的tif数据

2、计算景观指数

点击alysis parameter，选中moving window，勾选ndscape metrics景观指数，窗口大小根据需要调整

这里也要选中

pd代表破碎度，split代表分离度，division代表分维度。

点击run，点击proceed

三、arcgis计算景观干扰度

1、arcgis加载fragstats输出的破碎度、分离度和分维度tif图像

删除黑色背景（-999值）

ctrl+f 搜索 raster calculator

分别对3幅图像输入如图所示公式

2、归一化

还是打开raster calculator

以分维度为例输入：(“division”-valuemin)/(valuemax-valuemin)

3、计算景观干扰度

景观干扰度=0.5*破碎度+0.3*分离度+0.2*分维度

依旧打开raster calculator

输入：“pdnd”*0.5+”splitnd”*0.3+””division”*0.2

结果如图

使用渔网工具归一化操作略