下载地址:Downloads (fragstats.org)
参考景观格局指数计算方法及代表的生态学意义(待补充)
11fragstats 常用景观指数 - 知乎 (zhihu.com)
Fragstats指标
SHDI: Shannon多样性指数,在景观 Diverstiy 下,景观
Pi = 斑块类型(类)i 占据的景观比例。
(1)SHDI 等于负的所有的斑块类型的香浓总和,即每种斑块类型的比例丰度乘以该斑块比例。请注意,Pi 基于总景观面积 (A),不包括存在的任何内部背景。
(2)当景观仅包含 1 个斑块(即无多样性)时 SHDI = 0。SHDI 随着不同斑块类型数量的增加或斑块类型之间的面积比例分布变得更加公平而增加。
(3)香农多样性指数是社区生态多样性的一种流行衡量标准,适用于景观。香农指数比辛普森多样性指数对稀有斑块类型更敏感。
AWMPFD(FRAC_AM): 面积加权的平均斑块分形指标,Area-weighted mean patch
fractal dimension,在 Shape 下第三行,类型/景观
取值范围:1 ≤ FRAC ≤ 2。 对于周长非常简单的形状(例如正方形),FRAC 接近 1,对于具有高度复杂的平面填充周长的形状,FRAC 接近 2。
在一定程度上反映了人类活动对景观格局的影响,受人类干扰程度小的景观分维数值高,受人为干扰程度越大的景观分维数值越低【文献1】。
CONTAG: 蔓延度指数,Contagion,景观,% ,Aggregation 下
取值范围:0% ≤ Contag ≤ 100%。
高蔓延度值说明景观中的某种优势斑块类型形成了良好的连接性;反之则表明 景观是具有多种要素的密集格局,景观的破碎化程度较高【文献1】。
表示景观不同斑块类型的团聚程度或延展趋势【文献2】。
COHESION: 斑块内聚力指数,Patch Cohesion Index,类型/景观,Aggregation 下
官网没有详细描述,但范围极有可能也是0 % -100%。
描述景观中生物体的生境连接情况【文献3】。
凝聚度可衡量相应景观类型自然连接性的程度,关键类型占景观的比例减少并分割成不连接的斑块,则其值接 近于0,比例增加其值增加【文献4,p13】。
FN: 景观破碎化指数【文献5】
A 为景观总面积;Ni 为第 i 类斑块 的数量;Ni 也可作 N,即所有类型的斑块总数。该指数表示单位面积上某类斑块或所有类型斑块的个数。前者可比较不同时期同类斑块的景观破碎度,后者可比较不同时期的景观总体破碎度。值越大,景观破碎度越高。
N 即 Number of Patches(NP或NA),在Aggregation选项卡下
参考文献:
【1】青弋江流域土地利用/景观格局对水质的影响
【2】长江流域景观格局与生态系统水质净化服务的关系
【3】基于 Google Earth Engine 平台的湿地景观空间格局演变分析:以白洋淀为例
【4】云蒙湖流域土地利用变化的尺度效应研究
【5】基于景观破碎度分析的温州三垟湿地乡土景观保护研究
Arcgis处理
我的理解:
内部背景就是属于研究区的范围,但是你不计算这些 class(类型),那你就把这些class赋值为999,忽略它们不输出。
外部背景就是我们在arcgis中常见的 Nodata,给它赋值为一个负值即可。
(1)导出土地利用类型
将 nodata 设置为 -9999,注意不要采用在栅格计算器把nodata赋值为-9999的方式
(2)创建渔网格
工具箱 —— 数据管理 —— 采样 —— 创建渔网 (10.8版本)
选择输出位置,渔网大小选择土地利用范围
学习的视频中给出的单元格高宽单位是m,但我实际是度,根据研究区大小选择了0.075°
注:实际和数据框的坐标系有关,如果是投影坐标系,单位就是m,如果是地理坐标系,单位就是度,具体可以在数据框中查看。
(2)选择出与研究区相交的渔网格
选择 —— 按位置选择 —— 与原图层相交
给格网设置一个从 1开始的递增字段
(3)打散格网
分析工具 —— 提取分析 ——按属性分割
(4)批量 按掩膜提取
全部设置好之后在需要输数据的地方右键,勾选模型参数,然后重命名为自己认识的样子,保存之后就可以用了!arcgis中模型构建器的保存转为模型小工具
(5)fragstats批量导入tif
先导入一个,然后export batch导出 fbt 文件,用excel打开,数据自文本,按逗号分割,填充好之后另存为 csv 文件,import batch导入
(6)创建类型描述fcd文件,用txt写完把后缀名改为fcd,save run as导出,在excel中使用自文本加载,然后链接到shp文件的属性表,把N/A改为0,链接需要导出一下,才能把N/A变成0
不计算class用不到这个
(7)run,运行,save run as保存运行结果,
导出后逗号分隔,另存为 xlxs 文件,然后使用数据分割提取出 gridID
以上Arcgis处理攻略转载自B站老襄:BV1N54y1f71g