Google Earth Engine（GEE）计算降水量偏差

今天来简单分享下如何在GEE计算某一地区降水量的偏差

目标：

在GEE计算某一地区降水量的偏差，以百分比显示，并制作统计图

以河南省为研究区

GEE实现代码：

首先确定研究区和使用的数据集

var chirps = ee.ImageCollection("UCSB-CHG/CHIRPS/PENTAD")
// 选择的数据集为"UCSB-CHG/CHIRPS/PENTAD"
//CHIRPS的数据是以Pentad形式出现的
var roi = ee.FeatureCollection("users/lilei655123/HeNan_shi");
Map.centerObject(roi,7)
var styling = {color:"red",fillColor:"00000000"};
Map.addLayer(roi.style(styling),{},"geometry")

CHIRPS Pentad提供了1981年以来的数据,为全球降雨量数据集

确定计算月平均值的时期（10年）

// 我们可以计算出每月的平均值
var lpaYears = ee.List.sequence(2010, 2020);
var months = ee.List.sequence(1, 12);

创建一个每月的总数集合

var monthlyImages = lpaYears.map(function(year) {return months.map(function(month) {var filtered = chirps.filter(ee.Filter.calendarRange(year, year, 'year')).filter(ee.Filter.calendarRange(month, month, 'month'))var monthly = filtered.sum()return monthly.set({'month': month, 'year': year})})
}).flatten()
print(monthlyImages);

我们可以计算出所有年份中每个月的平均值

var shortTermMeans10y = months.map(function(month) {var filtered = monthlyCol.filter(ee.Filter.eq('month', month))var monthlyMean = filtered.mean()return monthlyMean.set('month', month)
})
var monthlyRainfall = ee.ImageCollection.fromImages(shortTermMeans10y)

现在我们采取2021年的数据，并计算出月平均值

var filtered = chirps.filter(ee.Filter.date('2021-01-01', '2021-12-31')).filter(ee.Filter.bounds(roi));
// 计算月平均降雨量
var monthlyTotals = months.map(function(month) {return filtered.filter(ee.Filter.calendarRange(month, month, 'month')).sum().set('month', month);
});

计算偏差，单位为%

var deviation = months.map(function(month) {var shortTermMean10y = monthlyRainfall.filter(ee.Filter.eq('month', month)).first()var monthlyObserved = observedRainfall.filter(ee.Filter.eq('month', month)).first()var deviation = (monthlyObserved.subtract(shortTermMean10y).divide(shortTermMean10y)).multiply(100).set('month', month)return deviation
})

创建统计图

var chart = ui.Chart.image.series({imageCollection: deviation, region: roi, reducer: ee.Reducer.mean(), scale: 10000,xProperty: 'month'
}).setOptions({interpolateNulls: true,lineWidth: 1,pointSize: 3,title: '降雨量均值偏差',vAxis: {title: '偏差 %'},hAxis: {title: 'Month', gridlines: {count: 12}}
});
print(chart);

运行结果如下:

完整代码请在公众号后台回复“0730降水量偏差”

感谢关注，欢迎转发！

声明：仅供学习使用！

希望关注的朋友们转发，如果对你有帮助的话记得给小编点个赞或者在看！

## ****更多内容请关注微信公众号“生态遥感监测笔记”**