GEE | 2000年至今MODIS蒸散发ET变化分析 (日/月/季节/年)

一、说明

MOD16A2.061是基于Penman-Monteith方程的8天复合蒸散量产品,空间分辨率为500米,时间范围自2001年起至今。基于GEE对ET进行时间变化分析,包括日、月、季节和年,并导出对应的栅格数据。

各位同学可以根据需求更改研究区、时间段和分辨率,当然也可以分析其他产品数据。

二、日变化特征

图片

图片

 

三、月尺度

以下是使用 Google Earth Engine (GEE) 计算 MODIS 散发的示例代码: ```javascript var startDate = '2019-01-01'; var endDate = '2019-12-31'; var region = ee.Geometry.Polygon([[73.5, 18.5], [135, 18.5], [135, 54], [73.5, 54], [73.5, 18.5]]); // 加载 MODIS 数据集 var modis = ee.ImageCollection('MODIS/006/MOD16A2') .filterDate(startDate, endDate) .filterBounds(region); // 定义计算散发的函数 var et = function(image) { var lw = image.select('LW_incoming_1000m'); var sw = image.select('SW_incoming_1000m'); var et = image.select('ET'); var t = image.select('LST_Day_1km').multiply(0.02).subtract(273.15); var esat = t.expression('0.6108 * exp(17.27 * b / (b + 237.3))', {'b': t}); var delta = t.expression('4098 * esat / pow(b + 237.3, 2)', {'b': t}); var ra = sw.expression('(solar_constant * 24 * 3600 / pi) * (d * sin(phi) * sin(delta) + cos(phi) * cos(delta) * sin(d))', { 'solar_constant': 1367, 'pi': 3.14159, 'phi': image.select('cos_phi'), 'delta': delta, 'd': lw.expression('acos(-1 * tan(phi) * tan(delta))', {'phi': image.select('cos_phi'), 'delta': delta}) }); var g = lw.subtract(ra).multiply(0.408); var et_image = g.divide(g.add(delta.multiply(900 / t.add(273.15).add(17.27).pow(2)))).multiply(24 * 3600 / 1000); return image.addBands(et_image.rename('ET')); }; // 计算散发 var etCollection = modis.map(et); // 取出散发数据 var etData = etCollection.select('ET'); // 显示散发数据 Map.addLayer(etData, {min: 0, max: 10, palette: ['blue', 'green', 'yellow', 'red']}, 'ET'); ``` 注:此代码计算的是 MODIS 散发数据集(MOD16A2)中的每散发数据,并在地图上显示。在代码中,我们首先加载了 MODIS 数据集,然后定义了一个计算散发的函数 et,最后对 MODIS 数据集应用该函数并取出散发数据。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

GIS遥感数据处理应用

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值