|
& w9 r# n T+ [5 o) }- x
InVEST（Integrated Valuation of Ecosystem Servicesand Tradeoffs）是美国自然资本项目组开发的、用于评估生态系统服务功能量及其经济价值、支持生态系统管理和决策的一套模型系统，它包括陆地、淡水和海洋三类生态系统服务评估模型。目前已经在生态系统服务领域得到广泛应用。
" N' A4 ~5 z2 ?9 t( l7 b. a0 V 生境质量（Habitat Quality）# H/ L' r' `; O2 b" z* b
生物多样性与生态系统服务的生产有着密切的联系。生物多样性本质具有空间化特征，因此，可以通过分析土地利用和土地覆盖图(LULC)及其对生物多样性威胁程度计算得到。我们可以利用InVEST模型，通过输入LULC地图、LULC对每种威胁的敏感性、每种威胁分布和密度的空间数据以及保护区域空间位置等数据，从而评估该地区景观的生境质量、退化水平和生境稀有性等。 7 k3 s- U7 y6 x% O3 n7 d* T+ c
运行原理
5 m/ P6 \ f+ o+ g' v) ]* L0 B, g 输入LULC数据后，首先要考虑哪种土地利用类型应该被列为生境（即生境的适宜性），用来表示类型j的生境适宜性。
6 J- i; E3 M, m9 G# w 模型还需要生境威胁密度数据和其对生境质量的影响，评估这些威胁源导致土地利用类型的退化程度，这里需要四个因素：
7 c; n/ j/ I. f8 n- l, o 1.每一种威胁的相对影响。例如，与农业用地相比，城市建设用地对周围生境质量的影响更强。因此，城市建设用地对生境破坏性的权重更大，这里权重代表威胁r对所有生境的相对破坏性，取值为[0,1]。
: k' J6 t( C9 P7 T$ p7 C 2.生境栅格与威胁的距离。通常，威胁随距离的增加而减小，这里我们可以选择线性或指数距离衰减函数来描述威胁在距离上的衰减： 8 {: T/ O8 K( E* T4 L2 F! X
1 u' \7 F8 d$ f
其中，
" x: Q5 A1 c3 B( j
' T [4 j) ?3 s1 q 表示在栅格y的威胁r对位置x的生境栅格的影响， ; K( H/ \& C H# \5 K
% h! `7 Q4 \4 X/ ^+ P% y 在是栅格x与y的线性距离，
& _4 S8 U# t8 X 2 M2 b# b9 D/ e% ?' z# a
是威胁r的最大作用距离。 . v4 E, X2 l! I7 k; h0 g( B
3.栅格的可接近性。如自然保护区、人类难以到达的高海拔区域等，在这些区域，由于难以接近，威胁对生境的影响非常有限（或无影响）。这里用
5 Q# z/ F2 S$ i+ v4 I
% f1 E9 y1 f9 j' c% O0 | 表示栅格的可接近水平，1表示完全可接近。
: m G9 s7 l, R8 \9 E 4.每种生境类型对每种威胁的相对敏感性。由于每种生境类型对威胁的响应都可能不同，因此需要计算敏感性，修正之前的计算，计算公式如下：
* ] {4 f7 |& [2 a , m/ K, Z0 D, O8 R1 n
其中， 1 B: ? P! q% r8 F+ g C- {' P
$ D4 U0 r% G2 f* `
为第j种生境中栅格x的总威胁水平，y是威胁r的栅格图上的所有栅格， , \0 j- K- e8 [$ ~$ M" [" Z
* j+ r4 e3 ]% a
指威胁r的栅格图中位于网格单元内的栅格， # l0 s: P6 K L) G
% p4 S N/ M6 T* K
表示LULC对威胁r的敏感性，值越接近1越敏感。
1 d; w4 U5 @7 Q! R 最后，根据上述值，得到土地利用类型j中的生境栅格x的生境质量
; S0 Z& }+ e' w* C
. O- N1 y% k- B# N0 q4 R# m ：
$ E9 z& I/ R4 y# A1 E+ M% M
2 l8 P0 P. y$ g* p1 t 其中，z=2.5、k为比例因子（常数）。
8 M: [* Q# d8 I1 O& k 运行界面及输入数据说明
5 p2 p" |( B& B5 ^& E8 F# O% {( Q' V' K
: j9 ?6 `- {5 @: a ^1 X8 f 这是生境质量模型的运行界面，其中至少需要输入4个数据。
0 _ t. z, V `0 w 1.当前LULC数据（必需），这个不多说。不过需要注意的是，输入的栅格数据必须含有地类代码，并在value字段下，而且数据需要定义投影，并具有线性单位且单位为米（可以在arcgis图层属性里查看，这个很重要，不然会报错）
' w+ n$ m- M( Z: i 2.未来LULC数据（可选）：表示未来的土地利用类型，要求同1。
9 E* W$ j+ y4 u( l- g. `6 a( @5 Y& y& R8 R 3.基准LULC数据（可选）：基础景观地图上的LULC类型GIS栅格数据，要求同1（注意：只有1、2、3都输入时，HQ模型才会计算生境稀缺性）。
$ | i0 C7 B6 Z8 s 4.威胁数据（必需）：需要输入一个.csv文件，具有严格的格式要求。表头必须为：THREAT（威胁的名称）、MAX_DIST（每种威胁的最大影响距离，对应原理中的 $ E; K4 W- f! l Z" m2 R
5 Q" S4 y" }. S3 @
）、WEIGHT（威胁的权重，对应原理中的
! H9 `( b; l* P0 ^1 H 5 A% F2 L, v& q5 \8 @
）、DECAY（距离衰减类型，线性还是指数）和CUR_PATH（当前LULC数据对应的威胁因子栅格地图的路径）。这里需要准备好所有威胁因子的栅格地图（威胁从LULC中提取即可。以建设用地为例：从LULC数据中重分类，将建设用地赋值为1，其他地类赋值为0，这里NODATA也要赋值为0） , m0 N! v) I8 c( ]: y! {! _, ]
1 P# @0 h' X% t" j4 \+ Q! l. z 如果输入了多个LULC数据，则要添加FUT_PATH、BASE_PATH两列，输入未来和基准LULC数据对应的威胁因子栅格地图路径，以_c（当前）、_f（未来）、_b（基准）区分。 ' k8 n3 \" R* J4 y: @# Z7 _
5.威胁的可接近性（可选）：对应前面原理的 / D) O" F! M* h6 i$ D, T
$ t) T- }. e8 y: K% _
，需要输入矢量数据。 ( r5 R+ E2 R( x% u0 ], n
6.生境类型及生境类型对威胁的敏感性（必需）：需要输入一个.csv表格，表格需要含有每种土地利用类型对不同威胁的敏感性数据。具有严格的格式要求。表头必须为：LULC（地类代码，必须和LULC栅格地图中的value一一对应）、NAME（地类名称）、HABITAT（生境适宜性，对应原理中的 * a$ x1 b! y7 [- b5 J
$ m) [0 [# x( ^1 _% L9 }
，数值为0到1）、L_THREAT1, L_THREAT2……（每种LULC对于每种威胁的相对敏感性，对应原理中的
& Y0 C3 u# o! L) J 2 O6 N8 G* b7 K7 V( s
。威胁因子有几个，这个就得有几列，且斜体部分必须与威胁因子命名一致）。
, W+ U3 g# B) ^7 o2 i 4 M2 A# V: N7 |( F, ~
（这里在运行时，有可能会提示缺少某个value，可能是NODATA值没有算上，一般在表格后把缺少的值补上，参数全为0就行了） . E" x. ?- _* ^& s7 r
7.半饱和参数（必需）：由于网格单元的退化水平使用半饱和度函数转换为生境质量，因此必须确定半饱和参数。对应原理中的参数k，默认值为0.5。k值的选择仅仅取决于生境质量得分值的扩散和居中趋势。确保对涉及相同的景观的所有运行使用相同的k值。
L3 N& h8 R j1 d 输出结果 - |0 Q4 s% D9 B) [2 q. ~, g
输出结果可以用arcgis打开。主要有这么几个输出结果： & T# t" U0 }: {; f
quality_代表生境质量，根据公式4计算得到。较高数值表示较好的生境质量。非生境的景观区域生境得分值为0。
) J, Q: Q' x3 g3 j deg_表示生境退化水平，根据公式3得到。1代表退化度高，0代表退化度低。非生境的栅格退化得分值为0。 ) O& G5 K; n' x# ]4 j7 y2 Z3 [
rarity_表示生境稀有性，越接近1表示越稀有（这个只有三个LULC数据都输入了才能生成）。 8 b/ g+ l0 R) F, G6 V# I- ~
% Z; w9 I8 i1 B, O/ K( r8 @
5 @" q2 R3 n( d
( p+ e6 d1 l8 w; E v& g( L 声明： 本号旨在传播、传递、交流，对相关文章内容观点保持中立态度。涉及内容如有侵权或其他问题，请与本号联系，第一时间做出撤回。
; N8 N- [, B* G; s 想获取更多数据资源：关注微信公众号：技术科研吧 # m5 L' V0 C. ^3 i V# z
* N0 h U9 y7 v
! G% e8 m% H7 ^: o: r" B- [& A
$ L, e% I! s5 D# |+ o
9 d7 t$ z1 M& q% V5 u | 
