一种快速评估人类活动对植被覆盖变化影响的方法技术

本发明专利技术公开了一种快速评估人类活动对植被覆盖变化影响的方法，以长时间序列MOD13Q1 NDVI为数据源，采用MRT软件对数据进行分析处理形成NDVI数据集，采用线性回归模型、最小二乘法趋势分析及残差分析法，评估人类活动对植被覆盖变化影响情况，本发明专利技术能够基于长时间序列MODIS NDVI数据及气象数据，快速、准确、客观开展大尺度的人类活动对植被覆盖变化影响研究，充分挖掘了遥感卫星数据与气象数据进行大区域尺度生态环境影响评估的潜力。

A rapid assessment of the effects of human activity on vegetation cover change

The invention discloses a fast evaluation method of human activities impact on the vegetation cover change, with long time sequence MOD13Q1 NDVI as data source, analyze the formation of the NDVI data set of data by MRT software, using linear regression model, least squares analysis trend and residual analysis, evaluation of human activities on vegetation cover change and the invention can long time series MODIS NDVI data and meteorological data based on fast, accurate and objective to carry out large-scale human activities on vegetation cover change research, fully excavate the evaluation of remote sensing satellite data and meteorological data of regional ecological environmental impact potential.

【技术实现步骤摘要】
一种快速评估人类活动对植被覆盖变化影响的方法
本专利技术属于环境质量评估领域，涉及一种环境质量评估方法，具体涉及应用长时间序列MODIS卫星遥感数据以及气象数据，评估人类活动对植被覆盖变化影响的方法。
技术介绍
植被是陆地生态系统的主体，是生态环境监测的指示器。植被覆盖变化同时受到气候和人类活动的双重驱动。随着城市化进程的不断加速，在局部区域人类活动对植被覆盖变化的影响程度甚至超过了气候变化，量化人类活动对植被覆盖变化的影响可以有针对性地为人类活动对生态环境质量评估提供数据支撑，同时也能够为区域经济可持续发展和生态环境建设提供技术支持。传统的研究方法大多数为实地考察，定性评估，这种方法耗时、费力，不适宜进行大区域尺度评估，且不能直观反映出人类活动对植被覆盖变化的影响区域。遥感技术的不断发展，为大尺度范围内的动态监测定性评估提供了可能。归一化植被指数(NormalDifferenceVegetationIndex,NDVI)是监测植被覆盖变化最常用的指标之一，同时也是衡量区域生态环境的重要信息。MODISNDVI数据具有时间和空间分辨率高、覆盖范围广、数据易获取等特点，这也为动态评估人类活动对植被覆盖变化的影响提供了坚实的数据保障。
技术实现思路
本专利技术的目的是提供了一种基于长时间序列MODISNDVI数据以及气象数据，经专业软件处理分析，快速准确评估人类活动对植被覆盖变化影响的方法。为了实现上述目的，本专利技术的技术方案为：提供一种快速评估人类活动对植被覆盖变化影响的方法，以长时间序列MODISNDVI为数据源，采用MRT软件对数据进行分析处理获得区域时间序列NDVI数据集，采用线性回归模型、最小二乘法趋势分析及残差分析法，评估人类活动对植被覆盖变化影响情况，具体步骤如下:(1)NDVI数据下载与预处理：a、下载研究区域时间序列的MOD13Q1NDVI，数据空间分辨率为250m，时间分辨率为16d，使用MRT软件对数据进行镶嵌拼接，将数据由HDF转换为Geotiff格式，投影坐标由Sinusoidal重投影到WGS84/Albers坐标系；b、采用ArcGIS中掩膜提取对NDVI数据进行裁剪，通过最大值合成法计算得到研究区域年最大NDVI值，形成研究区域NDVI时间序列数据集；(2)研究区域气象数据处理：a、搜集研究区域内所有气象站点的日均气温与日降水量数据，将日均气温进行算数平均计算获取气象站点的年均气温，对日降水量求和获取各气象站点的年降水量；b、利用ArcGIS软件中的克里格插值法对研究区域各气象站点的年均气温和年降水量进行空间插值，生成空间分辨率为250m×250m的时间序列年均气温和年降水量栅格数据集，投影到WGS84/Albers坐标系下；(3)计算时间序列预测NDVI值：建立NDVI和年均气温、年降水量的二元一次回归关系模型，在ArcGIS中利用公式(1)模拟计算每个栅格单元上的NDVI值；公式(1)如下所示：NDVI预测值＝β0+β1T+β2P(1)式中，NDVI预测值为由年降水和年均气温因子预测的NDVI值，β0为常数项，β1、β2为回归方程的待定系数，待定系数可以采用最小二乘法计算得到；T、P分别为年均气温、年降水量；(4)计算研究区域的人类活动对植被覆盖变化的影响：a、将步骤(1)中得到的年最大NDVI值，即通过遥感观测得到的真实NDVI值减去步骤(3)计算的预测NDVI值，得到时间序列NDVI残差数据集，若残差大于0，表明人类活动对植被生长有促进作用；若残差小于0，表明人类活动不利于植被生长；若残差等于0，表明人类活动对植被变化影响微弱；b、采用非参数化趋势度与Mann-Kendall检验方法分析人类活动对植被覆盖变化影响的变化趋势，在ArcGIS中利用公式(2)计算各像元的残差NDVI变化趋势，非参数化趋势度计算公式(2)如下所示：式中，β为残差NDVI变化趋势，i、j为时间序列，NDVIi、NDVIj分别表示第i、j时间的NDVI值，若β>0，表明残差NDVI呈上升趋势，即人类活动有利于植被生长；若β<0，表明残差NDVI呈下降趋势，即人类活动不利于植被生长；若β＝0，表明残差NDVI序列没有变化趋势；所述Mann-Kendall是一种非参数统计检验方法，用来判断变化趋势的显著性；假设时间序列(x1，x2，…，xn)是n个独立的、随机变量同分布的样本，定义检验统计量为S，S计算公式如公式(3)所示：式中，sign（)为符号函数，当Xi-Xj小于、等于或大于零时，sign(Xi-Xj)分别为-1、0或1；Z统计量的计算公式如公式(4)所示：在给定显著水平α下，当|Z|＞μ1-α/2时，表示序列在α水平下存在显著的变化，α一般取值为0.05。c、基于ArcGIS空间分析功能，统计各变化趋势空间分布及面积比例。本专利技术快速评估人类活动对植被覆盖变化影响的方法具有以下的有益效果：1、本专利技术能够基于长时间序列MODISNDVI数据及气象数据，快速、准确、客观开展大尺度的人类活动对植被覆盖变化影响研究，充分挖掘了遥感卫星数据与气象数据进行大区域尺度生态环境影响评估的潜力。2、本专利技术应用计算机图形处理技术以及空间化模拟技术，能够定量评估人类活动对植被覆盖变化影响的空间分布特征和影响程度。3、本专利技术克服传统方法耗时耗力定性研究人类活动对植被覆盖变化的局限性，实现了大尺度人类活动对植被覆盖变化的低廉、高效、定量化。附图说明图1为海南岛2001-2015年均预测NDVI空间分布图；图2为海南岛2001-2015年均NDVI残差空间分布图；图3为海南岛2001-2015年NDVI残差变化趋势图；图4为海南岛NDVI残差变化趋势Mann-Kendall显著性检验图；具体实施方式为了详细说明本专利技术的
技术实现思路
、构造特征、以下结合实施方式并配合附图作进一步说明。海南岛是我国重要的生态涵养功能区，良好的生态环境是国际旅游岛建设的突出优势。随着国际旅游岛建设的推进，新型工业化、城镇化建设步伐加快提速，人类活动对海南岛生态环境的变化产生了极大的影响，准确、客观地分析人类活动对海南岛植被覆盖变化的影响能够为相关部门合理规划布局以及海南岛可持续发展提供科学数据支撑。快速评估人类活动对海南岛植被覆盖变化影响的方法，包括以下步骤：(1)NDVI数据下载与预处理：a、从美国国家航空航天局NASA网站下载MOD13Q1NDVI数据，研究区域在全球正弦投影中的编号为h28v06和h28v07，时间跨度为2001年到2015年共15年共计690景数据。使用MRT软件对数据进行镶嵌拼接，将数据由HDF转换为Geotiff格式，投影坐标由Sinusoidal重投影到WGS84/Albers坐标系；b、采用ArcGIS中掩膜提取功能，基于研究区矢量文件对MODISNDVI数据进行裁剪，通过最大值合成法计算得到研究区域年最大NDVI值，形成海南岛年NDVI时间序列数据集；(2)研究区域气象数据处理：a、在中国气象数据共享服务网中下载海南岛所有气象站点的日均气温与日降水量数据，将日均气温进行算数平均计算获取气象站点的年均气温，对日降水量求和获取各气象站点的年降水量；b、利用ArcGIS软件中的克里格插值法对本文档来自技高网...

【技术保护点】
一种快速评估人类活动对植被覆盖变化影响的方法，包括以下步骤：(1)NDVI数据下载与预处理：a、下载研究区域时间序列的MOD13Q1NDVI，数据空间分辨率为250m，时间分辨率为16d，使用MRT软件对数据进行镶嵌拼接，将数据由HDF转换为Geotiff格式，投影坐标由Sinusoidal重投影到WGS84/Albers坐标系；b、采用ArcGIS中掩膜提取对NDVI数据进行裁剪，通过最大值合成法计算得到研究区域年最大NDVI值，形成研究区域NDVI时间序列数据集；(2)研究区域气象数据处理：a、搜集研究区域内所有气象站点的日均气温与日降水量数据，将日均气温进行算数平均计算获取气象站点的年均气温，对日降水量求和获取各气象站点的年降水量；b、利用ArcGIS软件中的克里格插值法对研究区域各气象站点的年均气温和年降水量进行空间插值，生成空间分辨率为250m×250m的时间序列年均气温和年降水量栅格数据集，投影到WGS84/Albers坐标系下；(3)计算时间序列预测NDVI值：建立NDVI和年均气温、年降水量的二元一次回归关系模型，在ArcGIS中利用公式(1)模拟计算每个栅格单元上的NDVI值；公式(1)如下所示：NDVI预测值＝β0+β1T+β2P  (1)式中，NDVI预测值为由年降水量和年均气温预测的NDVI值，β0为常数项，β1、β2为回归方程的待定系数，待定系数可以采用最小二乘法计算得到；T、P分别为年均气温、年降水量；(4)计算研究区域的人类活动对植被覆盖变化的影响：a、将步骤(1)中得到的年最大NDVI值，即通过遥感观测得到的真实NDVI值减去步骤(3)计算的预测NDVI值，得到时间序列NDVI残差数据集，若残差大于0，表明人类活动对植被生长有促进作用；若残差小于0，表明人类活动不利于植被生长；若残差等于0，表明人类活动对植被变化影响微弱；b、采用非参数化趋势度与Mann‑Kendall检验方法分析人类活动对植被覆盖变化影响的变化趋势，在ArcGIS中利用公式(2)计算各像元的残差NDVI变化趋势，非参数化趋势度计算公式(2)如下所示：...

【技术特征摘要】
1.一种快速评估人类活动对植被覆盖变化影响的方法，包括以下步骤：(1)NDVI数据下载与预处理：a、下载研究区域时间序列的MOD13Q1NDVI，数据空间分辨率为250m，时间分辨率为16d，使用MRT软件对数据进行镶嵌拼接，将数据由HDF转换为Geotiff格式，投影坐标由Sinusoidal重投影到WGS84/Albers坐标系；b、采用ArcGIS中掩膜提取对NDVI数据进行裁剪，通过最大值合成法计算得到研究区域年最大NDVI值，形成研究区域NDVI时间序列数据集；(2)研究区域气象数据处理：a、搜集研究区域内所有气象站点的日均气温与日降水量数据，将日均气温进行算数平均计算获取气象站点的年均气温，对日降水量求和获取各气象站点的年降水量；b、利用ArcGIS软件中的克里格插值法对研究区域各气象站点的年均气温和年降水量进行空间插值，生成空间分辨率为250m×250m的时间序列年均气温和年降水量栅格数据集，投影到WGS84/Albers坐标系下；(3)计算时间序列预测NDVI值：建立NDVI和年均气温、年降水量的二元一次回归关系模型，在ArcGIS中利用公式(1)模拟计算每个栅格单元上的NDVI值；公式(1)如下所示：NDVI预测值＝β0+β1T+β2P(1)式中，NDVI预测值为由年降水量和年均气温预测的NDVI值，β0为常数项，β1、β2为回归方程的待定系数，待定系数可以采用最小二乘法计算得到；T、P分别为年均气温、年降水量；(4)计算研究区域的人类活动对植被覆盖变化的影响：a、将步骤(1)中得到的年最大NDVI值，即通过遥感观测得到的真实NDVI值减去步骤(3)计算的预测NDVI值，得到时间序列NDVI残差数据集，若残差大于0，表明人类活动对植被生长有促进作用；若残差小于0，表明人类活动...

【专利技术属性】
技术研发人员：罗红霞，戴声佩，刘恩平，谢铮辉，方纪华，李茂芬，
申请(专利权)人：中国热带农业科学院科技信息研究所，
类型：发明
国别省市：海南,46