多尺度时间序列遥感影像趋势和断点检测方法与流程

本发明涉及遥感影像处理技术领域，尤其涉及一种多尺度时间序列遥感影像趋势和断点检测方法。

背景技术：

时间序列趋势和断点检测方法是目前进行遥感影像时间序列分析的重要方法之一，该技术方法体系为植被物候和陆地表层特征参数长时间序列分析提供了可行的方案，有助于揭示全球气候变化和人为活动下的陆地表层变化特征和生态环境变化过程，为自然资源可持续管理提供依据。由于受季节、趋势和噪声影响，植被指数时间序列数据为非平稳序列，主要包含了季节性变化、长期和短期波动等要素。

目前，已有许多方法可以从时间序列数据中分离出季节项、噪声项和趋势成分，并以此为基础进行时间序列趋势和断点检测以及分析生态系统的状态转换和路径。这些方法已经在遥感影像处理领域得到较为广泛的应用，为植被物候变化分析提供了重要基础。

但目前现有的时间序列趋势和断点检测方法在进行季节项、噪声项和趋势成分分离时候，忽视了时间序列数据存在的尺度效应，导致获得的时间序列趋势和断点是唯一尺度的结果，针对的情况比较单一。

技术实现要素：

针对现有技术中存在的问题，本发明实施例提供一种多尺度时间序列遥感影像趋势和断点检测方法及装置。

本发明实施例提供一种多尺度时间序列遥感影像趋势和断点检测方法，包括：

接收时间序列遥感影像，分析所述遥感影像中的时间序列信号，通过预设的离散小波对所述时间序列信号进行离散小波变化，得到分解信号；

获取预设的目标尺度，根据所述目标尺度提取所述分解信号，得到对应的尺度信号，并根据所述离散小波变化得到所述尺度信号的细节成分和近似成分；

根据所述近似成分和所述时间序列信号的长度计算所述时间序列信号的时间序列趋势；

根据所述细节成分计算并筛选得到所述时间序列信号的断点、断点类型及断点出现的时间。

在其中一个实施例中，所述方法还包括：

获取预设的线性回归模型，将所述近似成分和所述时间序列信号的长度代入所述线性回归模型，计算所述时间序列信号的时间序列趋势。

在其中一个实施例中，所述方法还包括：

通过如下线性回归模型的公式计算所述时间序列趋势：

其中，ai，是近似成分，xj是交叉小波变换近似成分ai的序号(j＝1,2,3,...,n)，n是时间序列信号的长度，ai,j是在j处的交叉小波变换近似成分ai值，s是近似成分ai的变化率；

当s大于0时，输出所述时间序列趋势增长的结果；

当s小于0时，输出所述时间序列趋势降低的结果。

在其中一个实施例中，所述方法还包括：

通过m-k检验法计算所述时间序列趋势的显著性，输出所述时间序列趋势的显著性计算结果。

在其中一个实施例中，所述方法还包括：

获取预设的峰/谷检测算法，通过所述细节成分及所述峰/谷检测算法检测所述时间序列信号的断点。

在其中一个实施例中，所述预设的峰/谷检测算法，包括：

其中，所述djⁱ是位置j的di，dj+1ⁱ、dj-1ⁱ是位置j+1、j-1的di，di是细节成分。

在其中一个实施例中，所述方法还包括：

通过预设的判别函数检测所述时间序列信号的断点的断点类型；

所述断点类型包括：峰型、谷型、突然减低型、突然增加型。

在其中一个实施例中，所述方法还包括：

获取预设的趋势局部变化函数及预设的阈值，将所述趋势局部变化函数的计算结果与所述阈值进行对比；

当所述结算结果大于所述阈值时，保存所述时间序列的断点的检测结果；

当所述结算结果大于所述阈值时，剔除所述时间序列的断点的检测结果。

本发明实施例提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述多尺度时间序列遥感影像趋势和断点检测方法的步骤。

本发明实施例提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述多尺度时间序列遥感影像趋势和断点检测方法的步骤。

本发明实施例提供的时间序列遥感影像趋势和断点检测方法及装置，提出一种基于离散小波变换的遥感影像多尺度时间序列趋势和断点检测方法，可以在不同尺度上进行时间序列趋势、断点类型以及断点出现时间的检测和分析，以适应复杂影响环境演化过程的分析。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例中多尺度时间序列遥感影像趋势和断点检测方法的流程图；

图2为本发明实施例中多尺度时间序列遥感影像趋势和断点检测装置的结构图；

图3为本发明实施例中电子设备结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图1为本发明实施例提供的多尺度时间序列遥感影像趋势和断点检测方法的流程示意图，如图1所示，本发明实施例提供了一种多尺度时间序列遥感影像趋势和断点检测方法，包括：

步骤s101，接收时间序列遥感影像，分析所述遥感影像中的时间序列信号，通过预设的离散小波对所述时间序列信号进行离散小波变化，得到分解信号。

具体地，服务器接收时间序列遥感影像，分析得到遥感影像中的时间序列信号，遥感影像可以为植物覆盖检测得到的图像，从遥感影像中得到的时间序列信号可以为跟随时间变化的植被的覆盖率、温度的变化率等等因素随时间的变化关系，然后通过离散小波对时间序列信号中的遥感影像逐像元进行分解，得到分解信号，其中，可以选择dmey等小波作为母小波进行离散小波变化。

步骤s102，获取预设的目标尺度，根据所述目标尺度提取所述分解信号，得到对应的尺度信号，并根据所述离散小波变化得到所述尺度信号的细节成分和近似成分。

具体地，获取预设的目标尺度，目标尺度可以根据生态环境的要求自由变化，比如尺度可以为一个月内温度的变化率，也可以为一百年内植被的覆盖率等等，根据目标尺度可以提取分解信号，将处在不同尺度的分解信号合并在一起得到尺度信号，并根据离散小波变化得到尺度包含的细节成分和近似成分，总的来说，细节成分可以大致表示信号的波动情况，而近似成分可以大致表示信号的波动趋势。

步骤s103，根据所述近似成分和所述时间序列信号的长度计算所述时间序列信号的时间序列趋势。

具体地，在尺度包含的近似成分的基础上，结合实践序列信号的长度可以计算得到时间序列信号的时间序列趋势，其中，时间序列趋势也包括趋势变化的显著性，计算方法可以为线性回归统计等，其中，具体的线性回归计算方法可以为获取时间序列信号的长度，以及预设的线性回归模型，通过线性回归模型，将近似成分与时间序列信号的长度代入其中的计算公式，计算公式可以为：

其中，ai，是近似成分，xj是交叉小波变换近似成分ai的序号(j＝1,2,3,...,n)，n是时间序列信号的长度，ai,j是在j处的交叉小波变换近似成分ai值，s是近似成分ai的变化率。

根据上述公式中s的计算结果，可以判断时间序列信号的时间序列趋势，即当s大于0时，可以输出时间序列趋势增长的结果；当s小于0时，输出时间序列趋势降低的结果。

步骤s104，根据所述细节成分计算并筛选得到所述时间序列信号的断点及断点类型。

具体地，根据细节成分可以计算得到时间序列信号的断点，因为细节成分可以代表时间序列的波动情况，包括波动方向的转变、波动幅度的转变等等，通过细节成分可以进一步判断时间序列信号的断点、断点类型及断点出现的时间，具体的判断方法可以通过获取预设的峰/谷检测算法，通过相对应的细节成分及峰/谷检测算法可以筛选时间序列信号的断点，峰/谷检测算法可以为：

其中，所述djⁱ是位置j的di，同理，dj+1ⁱ、dj-1ⁱ是位置j+1、j-1的di，di是细节成分。

上述检测算法中，di检测到的峰和谷则对应了原始时间序列信号的断点。

本发明实施例提供的一种多尺度时间序列遥感影像趋势和断点检测方法及装置，提出一种基于离散小波变换的遥感影像多尺度时间序列趋势和断点检测方法，可以在不同尺度上进行时间序列趋势、断点类型以及断点出现时间的检测和分析，以适应复杂影响环境演化过程的分析。

在上述实施例的基础上，所述多尺度时间序列遥感影像趋势和断点检测方法，还包括：

通过m-k检验法计算所述时间序列趋势的显著性，输出所述时间序列趋势的显著性计算结果。

在本发明实施例中，由于m-k检验法对样本的分布没有限制，且对于异常值的干扰不明显，所以在时间序列趋势检测基础上采用了m-k检验来验证时间序列趋势的显著性，具体的检测方法可以包括以下步骤：

假设该时间序列是n个相互独立的样本，对于时间序列信号ai的任意两个点k，j≤n，且如果k≠j，则可以定义检验统计量s

由于s服从正态分布，且均值为0，所以s的方差可以定义为：

当n>10，则定义标准的正态系统变量为：

通过m-k检验法，也可以检测时间序列趋势，即当z大于0时，可以认为时间序列的趋势是上升的，当z小于0时，时间序列的趋势是下降的。此外，当z的绝对值大于1.96，则可以认识时间序列趋势变化在95％的置信度上是显著的。

本发明实施例通过m-k检验法计算时间序列趋势的显著性，以适应复杂影响环境演化过程的分析。

在上述实施例的基础上，所述多尺度时间序列遥感影像趋势和断点检测方法，还包括：

通过预设的判别函数检测所述时间序列信号的断点的断点类型；

所述断点类型包括：峰型、谷型、突然减低型、突然增加型。

在本发明实施例中，由于离散小波变换近似成分ai-1可以由反映时间序列趋势的ai和反映时间序列断点的di重构，所以选择ai-1序列进行原始时间序列断点检测。本实施例构建了一个判别函数g(j)来分离ai-1序列的峰/谷以及突变：

其中，0是没有突变的序列，p是峰，v是谷，ad是突然下降，ai是突然上升，f(j)是运用di进行的峰/谷检测算法获得结果，f(j)＝1意味着在di时间序列i处存在断点。

本发明实施例通过判别函数分离时间序列趋势中的各个类型断点，以适应复杂影响环境演化过程的分析。

在上述实施例的基础上，所述多尺度时间序列遥感影像趋势和断点检测方法，还包括：

获取预设的趋势局部变化函数及预设的阈值，将所述趋势局部变化函数的计算结果与所述阈值进行对比；

当所述结算结果大于所述阈值时，保存所述时间序列的断点的检测结果；

当所述结算结果大于所述阈值时，剔除所述时间序列的断点的检测结果。

在本实施例中，在获取到时间序列信号的断点及断点类型之后，为了判断检测到的断点的有效性和可用性，构建趋势局部变化函数评价断点变化的趋势和程度：

对于g(j)不为0的时间序列断点，本研究基于断点j和下一个断点z之间ai-1的变化进行时间序列断点的程度评估。

另外，为了验证断点检测的准确性，在上述通过趋势局部变化函数评价断点变化的趋势和程度之后，还可以采用分段线性建模方法进行时间序列的模拟，并运用模拟的残差均分根(rmse)进行断点的评估。

具体的评估步骤为：通过预设的趋势局部变化函数h(i)以及预设的阈值(k)进行断点的评价，当h(i)大于k时，保存时间序列断点检测结果；当h(i)小于k时，剔除所述序列断点检测结果。通过时间序列原始信号和断点检测结果，基于预设的分段线性回归方法进行原始信号模拟，模拟的残差均分根(rmse)和预设的评估阈值(r，推荐设置为0.05)进行断点的评估。当rmse小于r时，输出时间序列断点检测结果，包括时间序列断点类型和对应时间，即峰/谷检测算法中的j；当h(i)大于r时，时间序列断点检测结果不通过检验，通过降低趋势局部变化函数检验中的阈值k进行循环，直到rmse达到小于r的要求。

本发明实施例通过趋势局部变化函数完成后续的断点筛选及模型评估。

图2为本发明实施例提供的一种时间序列遥感影像趋势和断点检测装置，包括：接收模块201、获取模块202、第一计算模块203、第二计算模块204，其中：

接收模块201，用于接收时间序列遥感影像，分析遥感影像中的时间序列信号，通过预设的离散小波对时间序列信号进行离散小波变化，得到分解信号。

获取模块202，用于获取预设的目标尺度，根据目标尺度提取分解信号，得到对应的尺度信号，并根据离散小波变化得到尺度信号的细节成分和近似成分。

第一计算模块203，用于根据近似成分和时间序列信号的长度计算时间序列信号的时间序列趋势。

第二计算模块204，用于根据细节成分计算并筛选得到时间序列信号的断点、断点类型及断点出现的时间。

在一个实施例中，装置还可以包括：

第二获取模块，用于获取预设的线性回归模型，将近似成分和时间序列信号的长度代入线性回归模型，计算时间序列信号的时间序列趋势。

在一个实施例中，装置还可以包括

公式模块，用于存储如下线性回归模型的公式计算所述时间序列趋势：

其中，ai，是近似成分，xj是交叉小波变换近似成分ai的序号(j＝1,2,3,...,n)，n是时间序列信号的长度，ai,j是在j处的交叉小波变换近似成分ai值，s是近似成分ai的变化率。当s大于0时，输出时间序列趋势增长的结果；当s小于0时，输出时间序列趋势降低的结果。

在一个实施例中，装置还可以包括：

第三计算模块，用于通过m-k检验法计算时间序列趋势的显著性，输出时间序列趋势的显著性计算结果。

在一个实施例中，装置还可以包括：

第三获取模块，用于获取预设的峰/谷检测算法，通过细节成分及峰/谷检测算法检测时间序列信号的断点。

在一个实施例中，装置还可以包括：

算法模块，用于存储预设的峰/谷检测算法，包括：

其中，djⁱ是位置j的di，同理，dj+1ⁱ、dj-1ⁱ是位置j+1、j-1的di，di是细节成分。

在一个实施例中，装置还可以包括：

检测模块，用于通过预设的判别函数检测时间序列信号的断点的断点类型；断点类型包括：峰型、谷型、突然减低型、突然增加型。

在一个实施例中，装置还可以包括：

第四获取模块，用于获取预设的趋势局部变化函数及预设的阈值，将所述趋势局部变化函数的计算结果与所述阈值进行对比，当所述结算结果大于所述阈值时，保存所述时间序列的断点的检测结果，当所述结算结果大于所述阈值时，剔除所述时间序列的断点的检测结果。

图3示例了一种电子设备的实体结构示意图，如图3所示，该电子设备可以包括：处理器(processor)301、存储器(memory)302、通信接口(communicationsinterface)303和通信总线304，其中，处理器301，存储器302，通信接口303通过通信总线304完成相互间的通信。处理器301可以调用存储器302中的逻辑指令，以执行如下方法：接收时间序列遥感影像，分析所述遥感影像中的时间序列信号，通过预设的离散小波对所述时间序列信号进行离散小波变化，得到分解信号；获取预设的目标尺度，根据所述目标尺度提取所述分解信号，得到对应的尺度信号，并根据所述离散小波变化得到所述尺度信号的细节成分和近似成分；根据所述近似成分和所述时间序列信号的长度计算所述时间序列信号的时间序列趋势；根据所述细节成分计算并筛选得到所述时间序列信号的断点、断点类型及断点出现的时间。

此外，上述的存储器302中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。

另一方面，本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各实施例提供的传输方法，例如包括：接收时间序列遥感影像，分析所述遥感影像中的时间序列信号，通过预设的离散小波对所述时间序列信号进行离散小波变化，得到分解信号；获取预设的目标尺度，根据所述目标尺度提取所述分解信号，得到对应的尺度信号，并根据所述离散小波变化得到所述尺度信号的细节成分和近似成分；根据所述近似成分和所述时间序列信号的长度计算所述时间序列信号的时间序列趋势；根据所述细节成分计算并筛选得到所述时间序列信号的断点、断点类型及断点出现的时间。

以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。