使用GEE计算KGE

问题背景

随机森林回归需要增加一个KGE的指标,但是网上这方面的资料比较少,于是自己按照公式在GEE上实现了一下。

具体代码

// 计算均值
var meanObs = observationValidation.reduce(ee.Reducer.mean(), [0]);
var meanPred = predictionValidation.reduce(ee.Reducer.mean(), [0]);

var meanObsRepeat = meanObs.repeat(0, observationValidation.length().get([0]))
var meanPredRepeat = meanPred.repeat(0, predictionValidation.length().get([0]))


// 计算标准差
var stdObs = observationValidation.reduce(ee.Reducer.stdDev(), [0]);
var stdPred = predictionValidation.reduce(ee.Reducer.stdDev(), [0]);

// 计算相关系数 r
var covariance = observationValidation.subtract(meanObsRepeat)
                  .multiply(predictionValidation.subtract(meanPredRepeat))
                  .reduce(ee.Reducer.mean(), [0]);
var r = covariance.divide(stdObs.multiply(stdPred));

// 计算相对标准差 β (Beta)
var beta = meanPred.divide(meanObs);

// 计算变异系数 γ (Gamma)
var gamma = stdPred.divide(stdObs);

// 计算KGE
var kge = ((r.subtract(1).pow(2))
           .add(beta.subtract(1).pow(2))
           .add(gamma.subtract(1).pow(2))
           .sqrt())
           .multiply(-1).add(1);

// 输出结果
print('KGE:', kge);

代码中,observationValidation和predictionValidation分别是实际值和预测值,按照如下公式进行计算。
在这里插入图片描述
在我的数据上进行计算,KGE可达0.92.

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值