[发明专利]草原植被覆盖时空动态变化的遥感定量监测方法

说明书

技术领域

本发明涉及植被覆盖动态遥感监测技术领域，特别是涉及一种基于植被NDVI数据的草原植被覆盖时空动态变化的遥感定量监测方法。

背景技术

植被是陆地表面最突出的土地覆盖类型，具有明显的季节变化与年际变化特征，其时空变化在全球能量循环及物质的生物化学循环中具有重要的影响，植被覆盖的变化研究在生态环境调查、水土保持研究、蒸散量研究等领域都有广泛的应用，同时，开展植被时空演变动态监测与评价，是评估生态保护措施、环境管理政策和全球变化研究的重要基础。近年来，随着全球气候变化及人为因素干扰，我国北方的草原荒漠化、草地退化现象严重，利用遥感手段监测草原植被时空动态变化特征，认识草原植被变化规律，对于草原植被生态系统保护及治理具有重要意义。

遥感技术作为一种技术手段可以长时间序列、大范围的获取地面植被覆盖信息。其中，归一化植被指数(Normalized DifferenceVegetation Index，NDVI)是目前应用较广泛的一种植被覆盖状况监测指标，它消除了大部分大气、仪器定标和地形等对指数的影响，增强了对植被的响应能力。经研究，植被NDVI数据与植被覆盖具有很高的正相关性，NDVI值越大，植被覆盖越高，因此，植被NDVI数据在一定程度上能够反映草原植被的覆盖状况，通过研究植被NDVI数据的动态变化，进而实现植被覆盖的时空动态研究。

而目前国内利用植被NDVI数据监测草原植被覆盖时空动态变化主要分析NDVI平均值变化情况，而未建立多指标的定量评价方法，这方面研究仍还有限，不能全面、深入、定量地支撑草原植被动态变化特征研究；而国外开展草原植被覆盖时空动态变化指标监测较少，研究缺少系统性，未能得到很好地应用。

发明内容

(一)要解决的技术问题

本发明要解决的技术问题是如何提高利用植被NDVI数据进行草原植被覆盖时空动态变化研究的有效性，克服现有技术中的定量研究的不足，形成基于植被NDVI数据的草原植被覆盖动态变化遥感定量监测方法。

(二)技术方案

为了解决上述技术问题，本专利提供一种草原植被覆盖时空动态变化的遥感定量监测方法，所述方法包括如下步骤：

步骤S1：预先设定统计年份的起始年份作为定量监测草原植被覆盖时空动态变化的时间范围；

步骤S2：利用植被的归一化植被指数NDVI数据制作预设定统计年份内每年的NDVI时序文件，然后消除NDVI时序文件的噪声干扰；

步骤S3：在所述消除噪声干扰后的时序NDVI文件中，计算得出植被NDVI数据中关于草原植被覆盖时空动态变化的若干个指标；

步骤S4：分别计算所述若干个指标各自在预设定统计年份内的年际动态变化率；

步骤S5：在预设定统计年份内，根据Mann-Kendall方法对所述若干个指标进行趋势分析，获取各指标在预设定统计年份内的统计检验值，获得各指标在预设定统计年份内的动态变化趋势；

步骤S6：基于所述若干个指标在预设定统计年份内的年际动态变化率及统计检验值，来定量研究草原植被覆盖的动态变化趋势。

所述步骤S2中，利用遥感影像处理工具处理植被NDVI数据，生成预设定统计年份内每年的NDVI时序文件；然后在NDVI时序文件噪声处，利用均值法，即计算其前后植被NDVI数据的平均值并以该平均值替换噪声值，来消除噪声干扰。

所述步骤S3中，

所述若干个指标包括：植被NDVI数据的年平均值AVE、植被NDVI数据的年最大值MAX、植被NDVI数据的年最小值MIN、植被NDVI数据的年最大值出现日期DMAX、植被NDVI数据的年最小值出现日期DMIN以及植被NDVI数据的季节动态值RREL；

其中，

所述植被NDVI数据的年平均值AVE表示每年草原植被覆盖的平均状况；

所述植被NDVI数据的年最大值MAX与植被NDVI数据的年最小值MIN代表草原植被覆盖最高与最低状况；

所述植被NDVI数据的年最大值出现日期DMAX以及植被NDVI数据的年最小值出现日期DMIN用旬数表示，分别代表植被覆盖最高值与最低值所对应的日期；

所述NDVI季节动态值RREL表示植被的年际变化状况，其计算公式为：

RREL＝[DMAX-DMIN]/[AVE]。

所述步骤S4中，