官网:https://naturalcapitalproject.stanford.edu/software/invest
软件【版本:3.90】:
网盘链接:https://pan.baidu.com/s/1T9EWK8yuMSy0uACViDjylQ
提取码:49s9
生境质量是指生态环境能够提供适合自然生态条件的能力,具有较强的地域性。生境质量的高低能反映区域生境的破碎程度以及对生境退化的抗干扰能力。计算结果范围为0 - 1,值越大代表研究区对生境退化的抗干扰能力越强。计算公式为:
式中,Qxj表示土地利用类型j中的栅格x的生境质量;Hj表示土地利用类型J的生境适合性;Dxj表示土地利用类型j中的栅格x所受胁迫水平;Z表示归一化常量;k是缩放常数。式中Dxj的计算公式如下:
式中,R为胁迫因子数量;r为胁迫因子;y为胁迫因子r的栅格数;Yr为胁迫因子所占栅格数;Wr为胁迫因子权重,取值范围为0 - 1;irxy为胁迫因子r对生境的每个栅格产生的影响(指数或线性);βx为生境抗干扰水平;Sjr为不同生境对每种胁迫因子的相对敏感性。
1 数据准备
1.1 威胁源
打开ArcGIS,导入土地利用数据
[spatial analyst tools]-[reclass]-[reclassify]
1.2胁迫因子属性表、敏感性表
通过查找研究区相关文献获得参数,且保存为“.csv”格式
2 路径设置
Workspace:工作空间,成果输出在选择的文件夹中
Curent land cover:土地利用数据
Threats Data:胁迫因子属性表
Accessibility to threats(可选):例如生态红线等矢量数据,保护区内的地区不受胁迫影响,生境质量不下降。
Senxitivity:敏感性表
Half-saturation Constant:默认值0.5
3 学习笔记
3.1胁迫源命名后缀原则
_c :现在胁迫源数据
_f :未来胁迫源数据
_b :基准胁迫源数据
3.2 threats 表应和胁迫源放在同一路径
3.3 threats 表的 Threat列应和威胁源的名称一一对应
3.4 senxitivity 表的LULC列应和土地利用数据的分类ID一一对应
3.5 土地利用数据下载自中国科学院资源环境科学数据中心(http://www.resdc.cn),分辨率1km * 1km。
3.6 若需加入道路等数据,则先把道路数据转为栅格数据,再按照胁迫源所用方法进行重分类