一种乡村生境质量评估及预测方法

技术领域

本发明属于城乡规划领域，具体涉及一种乡村生境质量评估及预测方法。

背景技术

生境质量是指生态环境的优劣程度，它以生态学理论为基础，在特定的时间和空间范围内，从生态系统层次上反映生态环境对人类生存及社会经济持续发展的适宜程度，根据人类的具体要求对生态环境的性质及变化状态的结果进行评定。乡村相对于城市往往拥有更优良的生态环境，但近年来，由于经济社会的快速发展，人类活动的过度干扰使乡村区域生境质量也面临威胁。生境质量有其自身内在的规律，掌握了这些规律，必然可正确的预测生境质量状况，因此，对乡村区域的生境质量进行量化评估，预测其未来状况，在此基础上制定针对性的生境优化策略不失为一种提升乡村地区生境质量的有效途径。

当前生境评估常常利用单一模型评估现状，而无法科学预测未来，因此无法根据未来发展合理安排对应政策；或是仅根据土地利用情况得出生境质量，而不考虑物种丰富度格局，得出的结果较为片面。本发明将综合模拟预测未来不同发展情景下的生境质量。

合理预测生境质量的未来发展状况，才能安排对应政策。而生境质量的未来发展状况与未来气候变化息息相关，其中排放情景是对潜在的辐射活跃排放物(如温室气体、气溶胶等)未来发展的一种可能的表述，也是预估未来气候变化的基础。基于未来大气辐射强度变化，联合国政府间气候变化专门委员会(IPCC)第五次气候变化评估报告中提出温室气体的典型浓度路径排放情景(RCPs)，可用于气候模式和土地利用变化等各项预测模型。RCPs以稳定浓度为特征，综合考虑未来气候变化、温室气体排放、社会经济变化和土地利用变化，提供了几种不同的发展情景，包括RCP2.6、RCP4.5、RCP6.0和RCP8.0四种；其中RCP4.5和RCP6.0均属于在政府干预下的中间发展情形，但由于RCP4.5情境下耕地面积减少，有可能超过最低红线，不符合实际，故不考虑此类情景。选择了其中RCP 2.6,RCP 6.0,RCP 8.5三种作为本专利研究的未来情景。

RCP2.6是温室气体浓度非常低的情景模式；全球范围内能源利用类型的改变，使温室气体排放显著减少，是全球作物面积增加最大的排放情景。RCP6.0是政府干预下的气候情景，在这种情境下至2100年，人口数目增至100亿；模拟各种政策的制定减少了温室气体的排放，但与RCP2.6相比，排放量的缓解程度较低，耕地的面积增长对森林面积的影响程度小。RCP8.5是在无气候变化政策干预时的基线情景，特点是温室气体排放和浓度不断增加，在此情景下，随着全球人口大幅增长、收入缓慢增长以及技术变革和能源效率改变导致的化石燃料消耗变大。

目前，国内外对生物多样性的评估和测算的方法主要集中于两种方式：一种是基于地面调查方法，重点是关注生物本底调查的评估研究；另一种是基于生态模型模拟的评价方法，如生态风险分析的EVR模型(Ecological Value at Risk)、生态系统服务和权衡的综合评估模型(InVEST)等。但单个模型预测的精度及时间跨度有所局限，例如InVEST模型无法预测未来生境状况，无法较好地用于评估土地利用变化及其导致的生态环境影响，使得InVEST模型用于预测未来一些特定情境下的生境质量变化成为了难题。

发明内容

针对现有技术的不足，本发明的目的在于提供一种乡村生境质量评估及预测方法。

本发明的目的可以通过以下技术方案实现：

一种乡村生境质量评估及预测方法，包括以下步骤：

S1，对研究区域进行等距网格划分；

S2，建立乡村生境质量评估数据库，数据库包括：乡村历史土地利用数据、当前社会经济数据、RCPs情境下的未来预测数据以及环境变量数据；

S3，建立Flus模型，将乡村历史土地利用数据作为影响数据，利用Flus模型筛选出训练样本，确定土地利用变化因子X，并得出各类用地的出现概率；

S4，将RCPs情境下未来气候预测数据和历史土地利用数据作为影响因素，在Flus模型内多次迭代模拟出，未来不同RCPs情境下的土地利用分布的栅格数据；

S5，将Flus模型得出的土地利用的栅格数据输入InVEST模型，来得出生境质量数据及威胁元分布数据；

S6，选取代表物种，并建立研究区域的物种地理分布数据集；

S7，利用S6得到的物种地理分布数据集及S1中的环境变量驱动MaxEnt模型，在RCPs情境下，分别对野生保护物种的潜在分布区进行模拟预测；

S8，提取模型预测结果阈值，绘制野生保护植物丰富度图，得到野生保护物种丰富度的空间分布的等级；

S9，将生境质量分布现状及未来预测和物种丰富度分布图进行叠加，得到生境质量综合评价图表。

进一步地，当前社会经济数据包括：人口数据、国内生产总值、距市中心距离以及路网；

未来预测数据包括各类用地百分比、气候、经济以及人口；

环境变量数据包括：19个气候因子数据、DEM高程数据、坡度数据、坡向数据以及NDVI数据。

进一步地，所述S3中，将历史土地利用、区位、自然环境、社会经济、气候因素作为影响数据，利用模型中人工神经网络的适宜性概率估算模块，筛选出训练样本，确定土地利用变化因子X，输入神经网络得出各类用地的出现概率。

进一步地，所述S3中，确定土地利用变化因子X的公式为：

X＝(x

式中X

进一步地，所述S4中，模拟出RCPs未来土地利用分布的步骤包括：

S41，将RCPs情境下未来气候预测数据、历史土地利用数据作为影响因素，初步得出RCPs情境下未来土地利用需求，

S42，基于自适应惯性机制的元胞自动机模块，定义自适应系数自动调整每类用地的惯性，该系数根据未来不同情景下的土地利用需求和当前各类用地的实际数量进行多次迭代，最终模拟出RCPs未来土地利用分布。

进一步地，所述S42中的迭代公式为：

式中：I

进一步地，所述S5中，生境质量数据的计算步骤为：

1)根据土地利用数据和表中数据进行生境退化度计算，表达式为：

式中，D

2)生境质量得分表达式为：

式中，Q

进一步地，所述S7中，模拟预测野生保护物种的潜在分布区步骤为：

S71，将收集到的研究区当前的物种的地理分布数据集及环境变量输入MaxEnt模型，通过已知物种出现点的分布数据和环境特征变量，算出目标物种分布的约束条件，探索生态需求下最大熵的可能分布；

S72，通过构建模型将模拟结果投射到研究区域，据此来预测目标物种在研究地的潜在生境分布及适宜性。

进一步地，所述S7中，在包含已知信息的前提下，假设随机变量α，包含A1，A2，A3，…，An共n种可能结果，则其熵值：

式中，H(α)为熵值，P1，P2，P3，…，Pn为每种出现的概率。

进一步地，于，所述S8中，绘制野生保护植物丰富度图的步骤为：

S81，将每个野生保护物种的最大熵运行阈值结果导入Arc GIS中，通过Arc GIS软件图层属性中自带的阈值划分功能对每个保护物种的丰富度格局进行区间划分，处理为10个区间，将每个区间可视化为不同灰度，灰度与物种丰富度呈负相关；

S82，对所选保护物种的潜在分布图层进行叠加，最后得出研究区野生保护植物物种丰富度图。

本发明的有益效果：能够对研究区域的未来生境质量作出较为准确的预测，并能将其与研究区现有生境质量与物种丰富度格局结合进行综合分析，从而为该区域生境优化策略得出较为科学全面的参考依据。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明的整体方法流程示意图；

图2是本发明实施例研究区域划分示意图；

图3是本发明实施例研究区域气候数据可视化示意图；

图4是本发明实施例研究区域人口密度数据可视化示意图；

图5是本发明实施例研究区域交通区位数据可视化示意图；

图6是本发明实施例研究区域地形数据示意图；

图7是本发明实施例研究与区NDVI数据示意图；

图8是本发明实施例未来RCP2.6、6.0、8.5情景下的土地利用示意图；

图9是本发明实施例当前和未来在RCP2.6情景下的生境质量数值示意图；

图10是本发明实施例研究区域野生保护植物分布点示意图；

图11是本发明实施例MaxEnt模型运行结果示意图；

图12是本发明实施例生境质量综合评价图表。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

如图1所示，一种乡村生境质量评估及预测方法包括以下步骤：

S1，参照国家统计局印发的《统计用区划代码和城乡划分代码编制规则》(国统字〔2009〕91号)，根据国家统计局官网上查询的城乡分类三位数代码，首位代码为1的区域划定为城镇，首位代码为2的区域划定为乡村；在Arcgis10.7中对城市和乡村区域进行划分；

明确研究对象与范围之后，对研究区域进行等距网格划分，选取1km*1km的网格，等距分布于研究区域各部分，并对网格进行编号，以西北角第一个网格为编号1，从左向右、从上向下依次进行编号，1、2、3、...、n。

S2，建立乡村生境质量评估数据库，进行数据预处理；数据库包括：乡村历史土地利用数据、当前社会经济数据(人口、国内生产总值、距市中心距离、路网)、未来预测数据(RCPs情境下各类用地百分比、气候、经济、人口)以及研究区域的23个环境变量数据(主要为19个气候因子数据、地形数据(DEM高程、坡度、坡向)和NDVI数据)；

其中，乡村历史土地利用数据从地理国情监测云平台获取，下载研究区过去20年历史土地利用数据，精度为30m，在Arc GIS中进行栅格化处理。根据不同省份土地利用情况进行土地利用类型简化；

当前人口数据从欧空局人口分布数据获取，下载研究区目前人口密度利用数据，精度为30m，在Arc GIS中进行栅格化处理；

国内生产总值从各省统计年鉴或全球变化科学研究数据出版系统获取，GDP以区为单位，在Arc GIS中进行栅格化处理；

距市中心距离从世界城市化前景网站(World UrbanizationProspects)获取，下载研究区目前距市中心距离数据，精度为30m，在Arc GIS中进行欧氏距离运算和栅格化处理；

路网从全球道路数据集网站(Global roads open assess data set)获取，下载研究区目前路网数据，精度为30m，在Arc GIS中进行栅格化处理。

未来气候模式采用RCP2.6、RCP6.0和RCP8.5三种排放情景，RCPs情境下温度、降水、日照强度均从WorldClim网站获取；例如，其中RCP2.6情景假设21世纪全球平均温度上升0.2-1.8℃，全球平均降水量增加2.1％；在Arc GIS中对当前温度、降水数据进行更改；

研究区所覆盖的生物气候因子数据从世界气候数据库(http：//www.WordClim.org)获取，共19个气候因子：bio1～bio19，不同研究区覆盖的气候因子数量不同；当前气候模式采用1950-2000年19个生物气候因子，如表1：

表1生物气候因子数据

DEM高程数据从地理空间数据云网站获取，下载研究区DEM影像，精度为30m；在ArcGIS中通过Arctoolbox中的Spatital Analyst tools-Surface-Slope和Aspect命令提取坡度、坡向，最终获得研究区域的海拔、坡度、坡向，分辨率均为30m×30m对下载的气候变量数据掩膜提取得到研究区范围内的图层；坐标系统一为WGS1984，栅格大小统一设置，最后转化成ASCII格式的文件用于模型运行；

NDVI(植被指数)数据是根据植被的光谱特征利用卫星可见光和近红外波段探测数据组合而来的，有助于识别植被存在度；首先从地理空间数据云下载Landsat 8影像数据，选取2020年同月份云量少的影像，在ENVI中通过band math计算NDVI，得到NDVI指数栅格图，裁剪出研究区范围，进而得到研究区的NDVI栅格数据，作为环境因子用于模型预测；坐标系统一为WGS1984，栅格大小统一设置，最后转化成ASCII格式的文件用于模型运行。

S3，根据Flus(未来土地利用模拟模型)模型要求，将历史土地利用、区位、自然环境、社会经济、气候因素作为影响数据，利用模型中人工神经网络的适宜性概率估算模块，筛选出训练样本，确定土地利用变化因子X，输入神经网络得出各类用地的出现概率；

数据输入Flus模型后，在基于人工神经网络的适宜性概率估算模块，采用随机采样的方法筛选出训练样本，确定驱动土地利用变化的因子X，作为神经网络中的输入层神经元(x

X＝(x

式中Xi为第1个采样点抽取的第i个驱动银子的变量，T为转置矩阵；

输出层中每个神经元都将生成一个0～1之间的值，该数值表示该像元发展为某类用地类型的概率，其值越高表示未来发展成该类用地的可能性越大；得出关于未来RCP2.6、4.5、8.5情境下各类用地出现概率的图像，该结果用于下一步骤。

S4，将RCPs情境下未来气候预测数据、历史土地利用数据作为影响因素，初步得出RCPs情境下未来土地利用需求，在基于自适应惯性机制的元胞自动机模块，定义自适应系数自动调整每类用地的惯性，该系数根据未来不同情景下的土地利用需求和当前各类用地的实际数量进行多次迭代，如式(2)所示，最终模拟出RCPs未来土地利用分布：

式中：I

S5，将Flus模型得出的土地利用的栅格数据输入InVEST(生态系统服务和权衡的综合评估模型)模型，根据乡村各类土地的特点，定义InVEST模型所需的威胁因子类型、胁迫距离、相对敏感度三个影响因素，运行InVEST得出生境质量数据及威胁源分布数据；具体步骤为：

参照国土资源部修订的《土地利用现状分类》(GB/T 21010—2017)，我国土地现状分类类型包括耕地、园地、林地、草地、商服用地、工矿仓储用地、住宅用地、公共管理与公共服务用地、特殊用地、交通运输用地、水域及水利设施用地、其他土地。因耕地、商服用地、工矿仓储用地、住宅用地、公共管理与公共服务用地、特殊用地、交通运输用地、其他土地(其他农用地)8种地类人类活动剧烈，将其作为非生境地类；而园地、林地、草地和水域及水利设施用地受人类活动影响较小的地类则作为生境地类；非生境地类在InVEST模型中列为威胁因子，对其量化分析：

第一个因子是每一种威胁源的对所有生境的相对破坏性；根据不同威胁类型对生境的破坏性给予退化源权重wr，能够选取0到1的任何一个数值，如表2所示：

表2威胁源权重

第二个因子是威胁源对各类生境的影响距离；威胁的程度随栅格与威胁源距离的增加而减小，因此距离威胁最近的那些栅格单元将受到较高的影响；根据InVEST规定，有线性和指数距离衰减两个函数来描述威胁在空间的衰减，如表3所示：

表3威胁源对各类生境的影响距离

第三个因子是每一种生境类型对每一种威胁源的相对敏感性，用于修正上一步计算的总影响；如表4所示：

表4威胁源的相对敏感性

根据研究区域存在的威胁源和生境类型制作表格(表5、表6)，并将各类威胁源的历史分布数据和未来分布数据在Arc GIS中进行栅格化处理，导出tif格式数据，将表格和tif数据一起输入InVEST模型；

表5威胁源数据

表6威胁源对各类用地影响

根据土地利用数据和表中数据进行生境退化度计算，表达式为：

式中，D

生境质量得分表达式为：

式中，Q

S6，以选取的代表物种(野生保护植物)建立研究区域的物种地理分布数据集，最终用于模型模拟；

从中国植物主题数据库(http://www.plant.csdb.cn/)获取野生保护植物名录；通过全球生物多样性信息网络(Global Biodiversity Information Facility，https://www.gbif.org/)获得各个物种的地理分布数据，将记录数少于五条的物种剔除，最终得到达到模型运行要求的物种，导入Excel中进行整理，去除重复点；将整理好的该省野生保护植物物种数据导入Arc GIS中，通过叠加分析得到研究区域内的野生保护植物分布点数据；将该数据导出，输入到Excel表中，并另存为*.csv格式，形成研究区域的物种地理分布数据集，每组数据包括该物种的学名及其分布点，具体到经纬度。

S7，利用收集到的研究区当前的物种的地理分布数据集及上述23个环境变量驱动MaxEnt(最大熵模型)模型，在RCPs情境下，分别对野生保护物种的潜在分布区进行模拟预测；

模型通过已知物种“出现点”的分布数据和环境特征变量，依据相应算法运算出目标物种分布的约束条件，探索生态需求下最大熵的可能分布，熵最大时物种的概率分布即满足物种生境条件的空间范围；通过构建模型将模拟结果投射到研究区域，据此来预测目标物种在研究地的潜在生境分布及适宜性；在包含已知信息的前提下，熵值最大时，冗余信息被排除，假设随机变量α，包含A1，A2，A3，…，An共n种可能结果，则其熵值：

式中，H(α)为熵值，P1，P2，P3，…，Pn为每种出现的概率。

S8，提取模型预测结果阈值，绘制野生保护植物丰富度格局，利用自然断点法分为10个等级，得到野生保护物种丰富度的空间分布的等级；

首先，将每个野生保护物种的最大熵运行阈值结果导入Arc GIS中，通过Arc GIS软件图层属性中自带的阈值划分功能对每个保护物种的丰富度格局进行区间划分，处理为10个区间，将每个区间可视化为不同灰度，灰度与物种丰富度呈负相关；其次，对所选保护物种的潜在分布图层进行叠加，最后得出研究区野生保护植物物种丰富度图。

S9，综合分析InVEST模型和MaxEnt模型评价结果，将两种模型模拟的生境质量分布现状及未来预测和物种丰富度分布格局进行叠加，按照叠加所得的生境质量综合评价图表，将生境质量的综合评价等级由高到低共划分为10个等级，灰度越小，数值越大，等级越高，则说明研究区生物多样性综合状况越好；

识别在未来不同的发展情景下，研究区最佳生境范围和生境优化重点区域，分析生境空间优先保护区域，最终为未来研究区生境优化策略提供较为科学的参考依据。

实施例：

本发明利用基于Maxent和InVEST模型的生境质量评估，在三种RCP情景下，通过模型的组合分析以得出研究区域生境质量的综合评估、预测研究区域生境质量在不同RCP情景下的未来状况，并以此为依据提出研究区具体的生境优化策略；

下面以某省为例，RCPs情景下基于Maxent和InVEST模型的生境质量评估的方法包括：

S1：参照国家统计局印发的《统计用区划代码和城乡划分代码编制规则》(国统字〔2009〕91号)，根据国家统计局官网上查询的城乡分类三位数代码，代码第一位为1的划分为城镇，第一位为2的划分为乡村；在Arcgis10.7中对该省城市和乡村区域进行划分；得到某省乡村用地区域，如图2所示。

S2：建立模型数据库，进行数据预处理；数据分为三类，分别应用于FLUS、InVEST、Maxent三个模型。

应用于FLUS模型的数据有：

1)该省历史土地利用数据(J):J12000年该省历史土地利用数据、J2 2005年该省历史土地利用数据、J32010年该省历史土地利用数据；

2)气候数据(Q)，作为FLUS模型驱动力因子：Q1月平均气温、Q2年平均温差、Q3年均降水量、Q4太阳辐射强度；如图3所示；

3)社会经济数据(S)，作为FLUS模型驱动力因子：S1该省分市地区生产总值、S2该省人口密度；如图4所示

4)交通区位数据(T):T1该省路网、T2到市中心距离，T3到城镇中心距离、T4到高速公路距离、T5到主干道距离、T6到铁路距离；如图5所示

未来预测数据(W)：RCP4.5、RCP6.0和RCP8.5情境下W1温度变化、W2降水变化、W3日照强度变化。

应用于InVEST模型的数据有：

威胁数据(X)：X1威胁因子，X2胁迫距离,X3各类土地对威胁因子的相对敏感度。

应用于MAXENT模型的数据有：

1)研究区国家重点野生保护植物地理分布数据；

2)19个气候因子数据；

3)研究区DEM高程数据；如图6所示；

4)研究区NDVI数据；如图7所示。

S3：利用上述数据库，按照FLUS使用手册说明定义气候、社会经济、交通区位定义为驱动力因子，利用该省2000、2005、2010年历土地利用数据，在Flus模型中人工神经网络的适宜性概率估算模块，采用随机取样的方法筛选出训练样本，确定驱动土地利用变化的因子X，作为神经网络的输入层神经元(xi，i＝1,2,…,n)：

X＝(x1(1),x2(1),…,xn(1))T(1)

式中Xi为第1个采样点抽取的第i个驱动银子的变量，T为转置矩阵。

输出层中每个神经元都将生成一个0～1之间的值，该数值表示该像元发展为某类用地类型的概率，其值越高表示未来发展成该类用地的可能性越大。得出关于未来RCP2.6、6.0、8.5情境下各类用地出现概率的图像，该结果用于S4。

S4：按照FLUS使用手册说明定义RCPs情境下各类用地百分比、未来人口经济数据，未来气候变化数据、历史土地利用数据作为影响因素，初步得出RCPs情境下未来土地利用需求。在基于自适应惯性机制的元胞自动机模块，输入步骤3得出的适宜性概率数据，输入RCPs情境下未来土地利用需求作为土地利用类型变化数量的目标，运行元胞自动机。元胞自动机将根据未来不同情景下的土地利用需求和当前各类用地的实际数量进行多次迭代(式2)，最终模拟出未来RCP 2.6、6.0、8.5情景下的土地利用分布；如图8所示；

式中：I

S5：参照国家标准化管理委员会和中华人民共和国国家质量监督检验检疫局引发的《土地利用现状分类》(GB/T 21010—2017)，我国土地现状分类类型包括耕地、园地、林地、草地、商服用地、工矿仓储用地、住宅用地、公共管理与公共服务用地、特殊用地、交通运输用地、水域及水利设施用地、其他土地。因耕地、商服用地、工矿仓储用地、住宅用地、公共管理与公共服务用地、特殊用地、交通运输用地、其他土地(其他农用地)8种地类人类活动剧烈，将其作为非生境地类；而园地、林地、草地和水域及水利设施用地受人类活动影响较小的地类则作为生境。选择该省广泛拥有的用地类型作为数据，其中非生境地类在InVEST模型中列为威胁因子，将各类威胁源的历史分布数据和未来分布数据在Arc GIS中进行栅格化处理，导出tif格式文件，将表格和tif数据一起输入InVEST模型。(表7、表8)

表7输入InVEST的threat文件

表8输入InVEST的sensitivitu文件

根据土地利用数据和表中数据进行生境退化度计算，表达式为：

式中，D

生境质量得分表达式为：

式中，Q

S6：从中国植物主题数据库(http://www.plant.csdb.cn/)获取该省野生保护植物名录，共13种。通过全球生物多样性信息网络(Global Biodiversity InformationFacility，https://www.gbif.org/)获得该13个物种的地理分布数据，将记录数少于五条的物种剔除，最终得到达到模型运行要求的物种共11种，导入Excel中进行整理。将整理好的该省野生保护植物物种数据导入Arc GIS中，通过叠加分析得到研究区域内的野生保护植物分布点数据，如图10所示。将该数据导出，输入到Excel表中，并另存为*.csv格式，形成研究区域的物种地理分布数据集，每组数据包括该物种的学名及其分布点，具体到经纬度。

S7：将所收集的研究区域当前的生物气候变量数据、高程数据和NDVI数据，输入MaxEnt模型，模型通过已知物种“出现点”的分布数据和环境特征变量，依据相应算法运算出目标物种分布的约束条件，探索生态需求下最大熵的可能分布，熵最大时物种的概率分布即满足物种生境条件的空间范围。通过构建模型将模拟结果投射到研究区域，据此来预测目标物种在研究地的潜在生境分布及适宜性。在包含已知信息的前提下，熵值最大时，冗余信息被排除，假设随机变量α，包含A1，A2，A3，…，An共n种可能结果，则其熵值：

式中，H(α)为熵值，P1，P2，P3，…，Pn为每种出现的概率。

S8：首先，将每个野生保护物种的最大熵运行阈值结果(如图11所示)，导入ArcGIS中，通过Arc GIS软件图层属性中自带的阈值划分功能对每个保护物种的丰富度格局进行区间划分，处理为10个区间，将每个区间可视化为不同灰度，灰度与物种丰富度呈负相关。其次，对所选保护物种的潜在分布图层进行叠加，最后得出研究区野生保护植物物种丰富度图。

S9：对InVEST模型进行与S8相同的处理，与研究区野生保护植物物种丰富度图进行灰度叠加，作为生物多样性的代表。按照叠加所得的生境质量综合评价图表，如图12所示；将生境质量的综合评价等级由高到低共划分为10个等级，灰度越小，数值越大，等级越高，则说明研究区生物多样性综合状况越好。

在此图表基础上，结合现有乡村区域筛选得到生物多样性保护优先地区，分析生境空间优先保护区域，最终为乡村区域生境优化提供直接有效的参考依据。

在本说明书的描述中，参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。