结果包含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