科学网

已有 22298 次阅读 2021-3-13 12:37 |个人分类:遥感|系统分类:科研笔记

Sen斜率估计用于计算趋势值，通常与MK非参数检验结合使用，即先计算Sen趋势值，然后使用MK方法判断趋势显著性。

Theil-Sen Median方法又被称为Sen斜率估计，是一种稳健的非参数统计的趋势计算方法。该方法计算效率高，对于测量误差和离群数据不敏感，常被用于长时间序列数据的趋势分析中。

β=mean(xj−xij−i),∀j>i

式中：xj 和 xi 为时间序列数据。β大于0表示时间序列呈现上升趋势；β小于0表示时间序列呈现下降趋势。

Mann-Kendall属于非参数检验方法，与其他参数检验的方法相比，不需要样本遵从一定的分布，受异常值干扰小，更适合顺序变量。Mann-Kendall检验已经在水文、气象趋势变化相关研究中得到了大量的成功应用，用于判断径流、降水、气候等的变化趋势的显著性。

R语言使用Raster包进行栅格计算，利用trend包sen.slope函数进行sen+mk的计算。对于下面代码的阅读，一定要看代码帮助！代码如下：