[发明专利]一种森林动态变化模式自动提取方法

权利要求书

1.一种森林动态变化模式自动提取方法，其特征在于：包括如下步骤，

S01：逐像元建立研究区多年逐日增强型植被指数以及积雪指数时序数据集；

S02：逐像元逐年分别提取植被覆盖强度指标、离散度指标以及持续度指标；

S03：逐像元逐年提取亮度指标；

S04：逐像元检测植被覆盖强度指标、离散度指标、持续度指标以及亮度指标的多年变化趋势；

S05：建立森林动态变化模式识别流程图；

S06：实现森林变化模式自动提取，获得研究区森林动态变化分布图。

2.根据权利要求1所述的一种森林动态变化模式自动提取方法，其特征在于：在所述步骤S01中，利用MODIS波段反射率数据，通过如下方式建立研究区每个栅格单元的多年逐日增强型植被指数MODISEVI以及归一化差分积雪指数MODISNDSI时序数据，计算MODISEVI2以及NDSI时序数据，分别为：

EVI2=2.5*(Nir-Red)(Nir+Red+1),]]>

NDSI=Green-SWIR6Green+SWIR6,]]>

其中，Nir，Red，Green，SWIR6分别为MODIS的近红外、红光、绿光以及中红外6波段的反射率；基于无云日计算得到MODISEVI以及NDSI时序数据，利用线性插值方法获得研究区多年的逐日MODISEVI以及NDSI时序数据，采用Whittakersmoother数据平滑方法，逐像元构建多年时空连续的逐日MODISEVI和NDSI时序数据集。

3.根据权利要求1所述的一种森林动态变化模式自动提取方法，其特征在于：在所述步骤S02中，逐像元逐年基于植被指数时序数据的高值区域的均值，构建该像元每年的植被覆盖强度指标；逐像元逐年基于植被指数时序数据的中高值区域的极差以及标准差，构建该像元每年的植被覆盖离散度指标，其中，DM＝(Max-Quan₂)(SD_m)，DM为植被指数中高值区域内数据的极差与标准差的乘积，Max为对应年份增强型植被指数时序数据的最大值，Quan₂为对应年份增强型植被指数时序数据的二分位数；逐像元逐年基于植被指数时序曲线上大于等于预设阈值的持续时间，构建该像元每年的植被覆盖持续度指标。

4.根据权利要求1所述的一种森林动态变化模式自动提取方法，其特征在于：在所述步骤S03中，逐像元逐年基于归一化差分积雪指数的年均值，构建该像元每年的亮度指标。