专利名称：草原植被覆盖时空动态变化的遥感定量监测方法
技术领域：
本发明涉及植被覆盖动态遥感监测技术领域，特别是涉及一种基于植被NDVI数据的草原植被覆盖时空动态变化的遥感定量监测方法。
背景技术：
植被是陆地表面最突出的土地覆盖类型，具有明显的季节变化与年际变化特征， 其时空变化在全球能量循环及物质的生物化学循环中具有重要的影响，植被覆盖的变化研究在生态环境调查、水土保持研究、蒸散量研究等领域都有广泛的应用，同时，开展植被时空演变动态监测与评价，是评估生态保护措施、环境管理政策和全球变化研究的重要基础。 近年来，随着全球气候变化及人为因素干扰，我国北方的草原荒漠化、草地退化现象严重， 利用遥感手段监测草原植被时空动态变化特征，认识草原植被变化规律，对于草原植被生态系统保护及治理具有重要意义。遥感技术作为一种技术手段可以长时间序列、大范围的获取地面植被覆盖信息。 其中，归一化植被指数(Normalized DifferenceVegetation Index, NDVI)是目前应用较广泛的一种植被覆盖状况监测指标，它消除了大部分大气、仪器定标和地形等对指数的影响， 增强了对植被的响应能力。经研究，植被NDVI数据与植被覆盖具有很高的正相关性，NDVI 值越大，植被覆盖越高，因此，植被NDVI数据在一定程度上能够反映草原植被的覆盖状况， 通过研究植被NDVI数据的动态变化，进而实现植被覆盖的时空动态研究。而目前国内利用植被NDVI数据监测草原植被覆盖时空动态变化主要分析NDVI平均值变化情况，而未建立多指标的定量评价方法，这方面研究仍还有限，不能全面、深入、定量地支撑草原植被动态变化特征研究；而国外开展草原植被覆盖时空动态变化指标监测较少，研究缺少系统性，未能得到很好地应用。
发明内容
(一)要解决的技术问题本发明要解决的技术问题是如何提高利用植被NDVI数据进行草原植被覆盖时空动态变化研究的有效性，克服现有技术中的定量研究的不足，形成基于植被NDVI数据的草原植被覆盖动态变化遥感定量监测方法。( 二 )技术方案为了解决上述技术问题，本专利提供一种草原植被覆盖时空动态变化的遥感定量监测方法，所述方法包括如下步骤步骤Sl 预先设定统计年份的起始年份作为定量监测草原植被覆盖时空动态变化的时间范围；步骤S2 利用植被的归一化植被指数NDVI数据制作预设定统计年份内每年的 NDVI时序文件，然后消除NDVI时序文件的噪声干扰；步骤S3 在所述消除噪声干扰后的时序NDVI文件中，计算得出植被NDVI数据中关于草原植被覆盖时空动态变化的若干个指标；步骤S4 分别计算所述若干个指标各自在预设定统计年份内的年际动态变化率；步骤S5 在预设定统计年份内，根据Marm-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中，所述植被NDVI数据的年平均值AVE、所述植被NDVI数据的年最大值MAX、所述植被NDVI数据的年最小值MIN以及所述植被NDVI数据的季节动态值RREL的年际动态变化率以变异系数CVt来表示；所述植被NDVI数据的年最大值出现日期DMAX以及所述植被NDVI数据的年最小值出现日期DMIN的年际动态变化率以标准差SDt来表示；
其中，CVt与SDt的计算公式分别为
权利要求
1.一种草原植被覆盖时空动态变化的遥感定量监测方法，其特征在于，所述方法包括如下步骤步骤Si 预先设定统计年份的起始年份作为定量监测草原植被覆盖时空动态变化的时间范围；步骤S2 利用植被的归一化植被指数NDVI数据制作预设定统计年份内每年的NDVI时序文件，然后消除NDVI时序文件的噪声干扰；步骤S3 在所述消除噪声干扰后的时序NDVI文件中，计算得出植被NDVI数据中关于草原植被覆盖时空动态变化的若干个指标；步骤S4 分别计算所述若干个指标各自在预设定统计年份内的年际动态变化率； 步骤S5 在预设定统计年份内，根据Marm-Kendall方法对所述若干个指标进行趋势分析，获取各指标在预设定统计年份内的统计检验值，获得各指标在预设定统计年份内的动态变化趋势；步骤S6 基于所述若干个指标在预设定统计年份内的年际动态变化率及统计检验值， 来定量研究草原植被覆盖的动态变化趋势。
2.如权利要求1所述的草原植被覆盖时空动态变化的遥感定量监测方法，其特征在于，所述步骤S2中，利用遥感影像处理工具处理植被NDVI数据，生成预设定统计年份内每年的NDVI时序文件；然后在NDVI时序文件噪声处，利用均值法，即计算其前后植被NDVI数据的平均值并以该平均值替换噪声值，来消除噪声干扰。
3.如权利要求1所述的草原植被覆盖时空动态变化的遥感定量监测方法，其特征在于，所述步骤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]。
4.如权利要求3所述的草原植被覆盖时空动态变化的遥感定量监测方法，其特征在于，所述步骤S4中，所述植被NDVI数据的年平均值AVE、所述植被NDVI数据的年最大值MAX、所述植被 NDVI数据的年最小值MIN以及所述植被NDVI数据的季节动态值RREL的年际动态变化率以变异系数CVt来表示；所述植被NDVI数据的年最大值出现日期DMAX以及所述植被NDVI数据的年最小值出现日期DMIN的年际动态变化率以标准差SDt来表示； 其中，CVt与SDt的计算公式分别为
5.如权利要求1所述的草原植被覆盖时空动态变化的遥感定量监测方法，其特征在于，所述步骤S5中，具体包括如下步骤步骤501 预设定置信度为95%，趋势检验的显著水平α = 0. 05 ； 步骤502 利用下式计算各个指标在预设定统计年份内的统计数据S
6.如权利要求5所述的草原植被覆盖时空动态变化的遥感定量监测方法，其特征在于，所述步骤S6中，定量研究草原植被覆盖的时空动态变化趋势的过程，具体包括如下步骤步骤S601 判断所述统计检验值&的趋势变化显著性；步骤S602 如果&为正数，且& > ，则表明草原植被覆盖有显著的上升趋势，如果&为负数，且，则表明草原植被覆盖有显著的下降趋势，如果|ZS|，则表明草原植被覆盖无变化趋势。
全文摘要
本发明涉及一种基于植被NDVI数据的草原植被覆盖时空动态变化的遥感定量监测方法，属于植被覆盖动态遥感监测技术领域。为了解决草原植被覆盖动态变化遥感监测定量研究较少的现状，本发明所提供的定量监测方法具体为利用NDVI时序文件计算获得反映草原植被覆盖状况的NDVI年均值、年最大值、年最小值、年最大值出现日期、年最小值出现日期和季节动态性等六个NDVI指标，通过对一定时期内上述NDVI指标进行年际动态变化率计算和Mann-Kendall趋势分析，定量监测草原植被覆盖的时空动态变化过程、物候特征及动态变化趋势。本发明适用于所有的植被NDVI数据。
文档编号G06F17/18GK102184162SQ20111003440
公开日2011年9月14日 申请日期2011年2月1日 优先权日2011年2月1日
发明者李营, 王桥 申请人:环境保护部卫星环境应用中心