基于GIS缓冲区分析的荒漠河岸林空间分布获取方法与流程

本发明涉及干旱区植被遥感监测技术领域，具体为一种基于gis（地理信息系统）缓冲区分析的荒漠河岸林空间梯度分布获取方法。

背景技术：

荒漠是西部干旱区的一种重要地理景观，它与山地及绿洲形成了干旱区特殊的耦合关系，构成了西部干旱区山盆体系的核心内容。其中，依水而生的荒漠河岸林生态系统则是绿洲系统的天然防护体系，对实现荒漠、山地与绿洲的协调发展具有重要意义。荒漠河岸林是干旱内陆河流域的独特景观，又称杜加依林，集中分布在干旱区内陆河流域沿岸，主要依靠地下水生存，是复杂而重要的生态系统，支撑着很多临界环境的变化过程，成为抵抗风沙、抑制沙漠化、维护区域生态平衡、保护生物多样性以及保障绿洲农牧业生产的重要屏障。依水而生，无水而亡，河流的洪枯变化、河道的摆动迁移以及河流泛滥所补给的地下潜水均对荒漠河岸林的形成、发展和衰亡具有重大影响。因此有必要对其空间分布进行实时的监测与分析，获取其保护范围，为干旱荒漠区的生态治理以及可持续发展提供有益参考。一些学者（陈永金，李卫红，刘加珍，等.输水堤防工程对塔里木河中游荒漠河岸林生态系统的影响[j].自然科学进展，2009，19(5):505-512.）主要通过样方法以及定点观测资料，在小尺度上分析了荒漠河岸林年际变化以及空间分布对河流各种转变的响应。该监测方法结果真实，但却有着一定的局限性。对于植被稀疏，受河水影响范围较广的荒漠河岸林来说，该方法监测范围较小，数据获取不便，效率较低，精度不高，且耗费人力物力，无法对植被的变化进行实时动态观察以及全面客观的比较分析。现如今，随着卫星遥感处理技术与gis地理信息系统的快速发展，一些学者采用遥感数据，在gis基础平台下，对荒漠区内陆河流域荒漠河岸林的植被覆盖度变化以及覆被变化进行了研究（万红梅，李霞，董道瑞，等.塔里木河下游植被覆盖变化遥感定量分析[j].自然资源学报，2013.）。但研究都局限于部分监测点数据，在大尺度上，却较少有对荒漠河岸林的连续梯度空间分布以及河水影响范围确定的分析（白元，徐海量，刘新华，等.塔里木河干流荒漠河岸林的空间分布与生态保护[j].自然资源学报，2013，28(5):776-785.），同时也缺乏相应高效便捷的梯度空间分布信息提取分析方法。

技术实现要素：

本发明提供了一种基于gis缓冲区分析的荒漠河岸林空间分布获取方法，目的是科学、高效、动态以及准确地获取多时段大尺度干旱区河岸林的空间梯度分布结构，以实时监测荒漠河岸林对河流变化的响应，对其保护范围确定提供依据，为干旱荒漠区的生态治理以及可持续发展提供有益参考。

为实现上述目的，本发明是通过以下技术方案实现的：

遥感是一种以物理手段、数学方法和地学分析为基础的综合应用技术，具有宏观、综合、动态和快速的特点。特别是在解决宏观尺度的环境问题时，卫星遥感可重复获取多种空间、不同时相和不同波谱分辨率的地球信息。现如今，对地观测技术的迅速发展和高分辨率商业卫星的投入使用，为生态环境监测提供了便捷。而gis地理空间信息系统又为遥感数据的提取分析提供了一个有利的平台，其特有的空间定位、缓冲区分析、叠置分析等空间分析功能为传统地物分布空间分析方法注入了信息与科技元素，使得这些方法更具科学性、计量性与信息性。因此，gis空间分析方法与遥感技术的结合是植被获取空间分布信息的有效途径。

本发明一种基于gis缓冲区分析的荒漠河岸林空间分布获取方法，其具体包括以下步骤：

步骤一：下载干旱区目标河流所在地的tm遥感影像，借用envi软件对遥感数据进行图像预处理、遥感解译、植被指数（ndvi）计算，提取荒漠河岸林ndvi遥感数据；

步骤二：在googleearth软件3d视窗中勾画河道边界矢量线，并导入gis基础平台，与遥感影像相应河道位置进行统一的配准与校正；

步骤三：在gis平台下，以河道边界线为基础，按具体研究要求生成其外围一定缓冲精度以及范围的多重缓冲区矢量图；

步骤四：在gis平台下，按照特定的坐标系统和参考信息，将多重缓冲区矢量图与荒漠河岸林ndvi遥感数据进行叠置，沿着距河道距离顺序提取各个缓冲带内ndvi均值，并对垂直河道距离自变量和ndvi均值因变量数据，采用曲线拟合进行荒漠河岸林水平空间梯度分布结构分析。

在上述步骤中，所述tm遥感影像由landsat系列卫星拍摄而成，空间分辨率为30m。

在上述步骤中，所述荒漠河岸林垂直河道水平分布结构的获取，是通过距河道距离不同距离下，ndvi这一环境因子的梯度分布来反映。

在上述步骤中，步骤一所述图像预处理进一步包括：

辐射定标、大气校正、几何校正、波段融合（红波、近红外波段）、影像拼接、研究区剪裁；

所述遥感解译进一步包括：地物分类、精度评价以及荒漠河岸林遥感影像掩膜提取；

进一步地，所述荒漠河岸林遥感影像掩膜提取进一步包括：

基于地物分类结果，在envi软件中建立荒漠河岸林掩膜文件，将其导入gis平台，与遥感影像进行叠置，再利用gis栅格剪切工具提取荒漠河岸林遥感影像。

所述植被指数（ndvi）计算包括：

利用envi软件的bandmath工具计算提取出的荒漠河岸林遥感影像ndvi值，计算模型为：

上式中：b为近红外波与红波的地表反射率；nir为近红外波段；red为红波段；

随后再用bandmath工具除去ndvi异常值，使得ndvi范围在[0，1]，计算公式为：，其中b1为ndvi值，此公式是将大于等于1的数值归为1，将小于等于0的数值归为0。

在上述步骤中，步骤二所述河道边界矢量线是在googleearth的3d视窗中，用添加多边形工具按最大比例尺，沿着河道边界画出的闭合矢量线，该矢量线自带wgs-84地理坐标；

所述配准与校正是在gis基础平台下，将河道边界矢量线与遥感栅格数据按照一定标准赋予的统一空间参考投影信息和坐标系统，对其进行空间配准，使得河道边界矢量线与遥感图像中的河道边界位置具有一一对应的空间关系。

在上述步骤中，步骤三所述多重缓冲区建立在河道边界线外围，由数条宽度（缓冲精度）一定又连续的缓冲带组成，其中根据研究需要的不同，缓冲范围灵活，缓冲精度也可任意设定。

与现有技术相比，本发明的有益技术效果在于：

本发明提供了一种结合了遥感技术与gis空间分析功能的河岸林空间梯度分布的获取方法，借此可了解距离河道各距离的植被生长情况，确定河岸林水平梯度分布结构及其保护范围，适用于干旱区大尺度区域河流生态环境研究。该方法效率高、数据获取便捷准确、人力物力耗费少，可实时动态的对环境进行监测，为干旱荒漠区的生态治理以及可持续发展提供有益参考。

附图说明

图1为本发明一种基于gis缓冲区分析的荒漠河岸林空间分布获取方法的技术路线图；

图2为2000年及2014年地物分类结果；

图3为googleearth3d视窗中的河道边界线；

图4为研究区2014年荒漠河岸林ndvi遥感影像与1km缓冲精度的20条缓冲带矢量图的叠置分析；

图5为2014年20km范围内荒漠河岸林ndvi均值垂直河道空间梯度分布趋势。

具体实施方式

下面结合附图对本发明的具体实施方式进行详细说明。

塔里木河两岸拥有世界上面积最大、分布最密集的以胡杨林为主体的荒漠河岸林植物群落。其丰富的产出，为生物多样性提供适宜生境，为水生生物提供有机食物，有着强大的生态功能，为荒凉的内陆盆地带来了一片生机。然而，源流区人口的增长以及对水资源的不合理开发与浪费，使得塔里木河干流荒漠河岸林生态系统受到了一定影响。随着时间的增长，为研究塔里木河干流荒漠河岸林垂直河道的空间梯度分布结构变化，下面以塔里木河中游河段为研究对象，对本发明进行实例分析，技术路线参见图1。

步骤一：数据下载与处理

1）从地理空间数据云（geospatialdatacloud）网站下载覆盖研究区域的行列号分别为143/31、143/32、144/31、144/32的landsattm遥感影像，所述tm遥感影像由landsat系列卫星拍摄而成，空间分辨率为30m，易于获得，且适用于荒漠植被稀疏地区。为达到可比性，影像采集时间分别为2000年和2014年的8-9月份，为塔里木河丰水期，云覆盖率<5%，植被覆盖度高；

2）借用envi软件对遥感数据进行图像预处理，包括辐射定标、大气校正、几何校正、波段融合（红波、近红外波段）、影像拼接和研究区裁剪；

3）借用envi软件对遥感数据进行遥感解译，包括地物分类、精度评价以及荒漠河岸林遥感影像掩膜提取，具体地，采用最大似然比(maximumlikelihood)监督分类算法，结合googleearth同期遥感影像以及研究区特征地物点位信息，根据目视判读方法对遥感影像进行地物分类，由于塔里木河中游靠近河道地区人烟较少，因此土地利用主要分为天然植被、耕地、未利用地以及其他地物四类，2000年及2014年地物分类结果参见图2；

4）采用随机抽样的方法进行精度评价，参考googleearth同期高分辨率遥感影像以及1：10万土地利用现状数据选择样本点，通过图像判读建立混淆矩阵，计算相关指标进行评价。结果表明：2000年、2014年的总精度和kappa指数分别为84.4%、85.4%和0.77、0.79，达到最低判别精度0.7的要求，能够满足大区域宏观土地覆被分析的需要；

5）基于遥感影像地物分类结果，在envi软件中建立荒漠河岸林掩膜文件，将其导入gis平台，与遥感影像进行叠置，再利用gis栅格剪切工具掩膜提取荒漠河岸林遥感影像；

6）在envi软件中，利用bandmath工具计算研究区荒漠河岸林遥感影像ndvi（植被指数）值，计算模型为：

上式中：b为近红外波与红波的地表反射率；nir为近红外波段；red为红波段；

随后再用bandmath工具除去ndvi异常值，使得ndvi范围在[0，1]，计算公式为：，其中b1为ndvi值，此公式是将大于等于1的数值归为1，将小于等于0的数值归为0。

步骤二：河道边界矢量线提取以及图像配准与校正

在googleearth的3d视窗中，用添加多边形工具按最大比例尺，沿着河道边界画出闭合矢量线（参见图3），该矢量线自带wgs-84地理坐标，再以kml的格式导出。

配准与校正，是在gis基础平台下，利用kml转矢量图层工具导入河道边界线矢量文件，与遥感栅格数据按照一定标准赋予的统一空间参考投影信息和坐标系统，对其进行空间配准，使得河道边界线与遥感图像中的河道边界位置具有一一对应的空间关系。

步骤三：基于gis的缓冲区梯度分析

在gis平台下，以河道边界矢量线为基础，按具体研究要求生成其外围一定缓冲精度（1km）以及范围的多重缓冲区矢量图，多重缓冲区建立在河道边界线外围，由数条宽度（缓冲精度）一定又连续的缓冲带组成，其中根据研究需要的不同，缓冲范围灵活，缓冲精度也可任意设定。具体地，生成其外围1km缓冲精度、范围为20km的20条缓冲区矢量图，参见图4。

步骤四：提取缓冲带的ndvi均值信息，进行荒漠河岸林水平空间梯度分布结构分析

在gis平台下，利用转栅格工具对缓冲区进行栅格化，按照特定的坐标系统和参考信息，将河道外多重缓冲区与荒漠河岸林ndvi遥感数据进行叠置分析，沿着距河道距离顺序提取各个缓冲带内ndvi均值，然后以垂直河道距离为自变量，ndvi均值数据为因变量，进行曲线拟合，以此来反映荒漠河岸林的水平空间梯度分布结构。

具体地，将2000年与2014年荒漠河岸林空间梯度分布进行比较，分析塔河中游荒漠河岸林在14年间的变化。参见图5。

本发明一种基于gis缓冲区分析的荒漠河岸林空间分布获取方法，该方法结合了遥感技术与gis空间分析功能，对荒漠河岸林的连续梯度空间分布以及河水影响范围确定的分析，效率高，数据获取便捷准确，大大减少了人力物力耗费，并且可以实时动态的对环境进行监测。

尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。