GEE：非参数检验（Non-Parametric Trend Analysis，Sen+Mann-Kendall）

追风少年6626

已于 2023-11-17 21:55:00 修改

阅读量437

点赞数 1

分类专栏： gee 文章标签：云计算

于 2023-06-28 11:28:49 首次发布

本文链接：https://blog.csdn.net/sdjkefsda/article/details/131433025

版权

gee 专栏收录该内容

4 篇文章 3 订阅

订阅专栏

结果包含MK检验、Sen斜率、方差Var、显著性（z、P-value）

遇到的一些问题，需要留意：

求取senslope时，x轴对应的时间长度要注意，官方给的是’day'，一般分析年际变化

var slope = function(i, j) { // i and j are images
  return ee.Image(j).subtract(i)
      .divide(ee.Image(j).date().difference(ee.Image(i).date(), 'days'))
      .rename('slope')
      .float();
};

MK中求符号时，官方给的代码如下：其含义是当两个时段内差值的绝对值大于1时，才会进行赋值。根本原因是，官方用的数据集有比例因子，其两个时段内差值基本在几百，很容易得到-1，0，或1。但是当你在对如NDVI、EVI进行计算时，由于其差值在0.001左右，sign返回值会一值为0。而真实的sign函数，是<0时，赋值-1，>0时，赋值1。故该函数需要修改，可以修改为
```
var sign = function(i, j) { // i and j are images
  return ee.Image(j).neq(i) // Zero case
      .multiply(ee.Image(j).subtract(i).clamp(-1, 1)).int();
};
```
```
var sign = function (i, j) {
    var my_sign = ee.Image(j).neq(i)
        .multiply(ee.Image(j).subtract(i));
    return my_sign.where(my_sign.gt(0),1).where(my_sign.lt(0), -1)
};
```

代码出处：GEE社区

源代码地址：https://code.earthengine.google.com/89e83a99e7c0bc90c5c66ee03746becb