The invention discloses a vegetation monotone change trend detection method, which uses the EMD decomposition method to decompose the NDVI time series reconstructed by each pixel to obtain the trend component, and then carries out the monotonicity test on the trend component to obtain the trend result of each pixel. It does not depend on the technique of first order regression, which must be solved by least squares, and no matter what the seasonal trend or local trend is, the monotonic increasing or decreasing trend of the whole sequence can be obtained. The invention also provides a vegetation monotone change trend detection system, a device and a computer readable storage medium, which can also realize the above technical effect.
【技术实现步骤摘要】
一种植被单调变化趋势检测方法及相关装置
本专利技术涉及遥感数据数据分析领域,更具体地说,涉及一种植被单调变化趋势检测方法、系统、装置及计算机可读存储介质。
技术介绍
遥感长时间序列数据提供了关于全球地表动态变化以及变化趋势等重要信息,这些信息的有效提取可以为生产、社会发展、环境保护等各项重要决策提供依据。特别是植被的长期变化信息是气候变化中的关键因素,植被作为陆地生物圈的主要组成部分,植被变化趋势被认为与土地退化关联,植被状态常用于评估自然和农用土地的生产率和退化。在几十年的时间框架内,植被的突然减少,一般认为是由一些短期过程导致的,如火灾、农作物收获或灾害等,植被的突然增加则被认为可能由降雨事件或雪盖的减少引起;植被的渐变则认为是体现了植被对全球变化的适应过程,如大洋振荡、持续的气候变化、年际降雨减少或大气中二氧化碳浓度增加等。植被的渐变,单调地变得更绿或更不绿,这种单调变化的趋势也称为“绿化”或“褐化”。遥感数据中表征植被信息的特征是植被指数(NDVI,NormalizedDifferenceVegetationIndex)。植被“绿化”或“褐化”的趋势可以通过NDVI的变化趋势表达。因此植被“绿化”或“褐化”信息可通过检测NDVI单调变化的趋势获得。一般的时间序列趋势检测常用的方法有一元回归变化斜率法、Sen趋势度估计方法、曼肯德尔(Mann-Kendall)检验方法及季节性Mann-Kendall方法等。NDVI序列不同于一般时间序列的特性有两点:强烈的季节性特点和数据在时间的有相关性特点。数据在时间上的相关性特点会使得由一元回归变化斜率法得到的N ...
【技术保护点】
1.一种植被单调变化趋势检测方法,其特征在于,包括:提取NDVI时间序列图像中每个像元的时间序列;对每个像元的所述时间序列进行重构得到每个像元的重构NDVI时间序列;对每个像元的所述重构NDVI时间序列进行EMD分解,得到对应每个像元的所述重构NDVI时间序列的趋势分量;对每个所述趋势分量进行单调性检验,得到每个像元的趋势结果。
【技术特征摘要】
1.一种植被单调变化趋势检测方法,其特征在于,包括:提取NDVI时间序列图像中每个像元的时间序列;对每个像元的所述时间序列进行重构得到每个像元的重构NDVI时间序列;对每个像元的所述重构NDVI时间序列进行EMD分解,得到对应每个像元的所述重构NDVI时间序列的趋势分量;对每个所述趋势分量进行单调性检验,得到每个像元的趋势结果。2.根据权利要求1所述的方法,其特征在于,所述对每个所述趋势分量进行单调性检验,得到每个像元的趋势结果之后,还包括:将所述每个像元的趋势结果进行可视化处理得到每个像元的单调性趋势图。3.根据权利要求2所述的方法,其特征在于,所述将所述每个像元的趋势结果进行可视化处理得到每个像元的单调性趋势图,包括:将每个像元的上升趋势使用绿色调表示绿化,下降趋势使用黄色调表示褐化,利用色调饱和度表示趋势结果的趋势显著程度。4.根据权利要求1所述的方法,其特征在于,所述对每个像元的所述时间序列进行重构得到每个像元的重构NDVI时间序列,包括:S301,确定每个像元的所述时间序列中的异常点;S302,将每个像元的时间序列中的每个所述异常点更换为正常点,得到每个像元更新后的NDVI时间序列,作为第一时间序列;S303,利用第一S-G滤波获得的缓变曲线与每个所述第一时间序列进行对比,确定每个第一时间序列中在相同时间点低于所述缓变曲线值的点作为异常值,并将每个所述异常值更换为所述缓变曲线中与所述异常值对应的时间点的值,得到每个更新后的NDVI时间序列,作为第二时间序列;S304,利用第二S-G滤波对每个所述第二时间序列进行滤波得到滤波后的第二时间序列;S305,利用每个所述第一时间序列与对应的每个所述滤波后的第二时间序列计算对应每个第一时间序列与第二时间序列的拟合残差指数,将本次迭代后的残差指数小于上次迭代后的残差指数的时间序列作为第一时间序列,返回S303;将本次迭代后的残差指数不小于上次迭代后的残差指数的时间序列作为重构NDVI时间序列。5.根据权利要求4所述的方法,其特征在于,所述将每个像元的时间...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。