GEE | 1901-2022年全球0.5°标准化降水蒸散指数SPEI数据集

图片

Product: Beguería, Santiago; Vicente Serrano, Sergio M.; Reig-Gracia, Fergus; Latorre Garcés, Borja; 2023; SPEIbase v.2.9 [Dataset]; DIGITAL.CSIC; Version 2.9. doi:10.20350/digitalCSIC/15470

一、数据说明

SPEIbase(全球标准化降水-蒸散指数数据库)提供了全球尺度上干旱状况的长期可靠信息,空间分辨率为0.5度,时间分辨率为月度,可支持从1至48个月的多时间尺度分析。标准化降水-蒸散指数(SPEI)通过标准化变量量化当前气候平衡相对于长期平衡的偏差,确保了时空上的可比性。

二、GEE代码

以SPEI_01_month为例,print月度SPEI变化并导出栅格数据。

图片

图片

<
### 如何在 Google Earth Engine (GEE) 上计算并下载某地区的标准化降水蒸散指数 (SPEI) #### 背景说明 标准化降水蒸散指数 (Standardized Precipitation Evapotranspiration Index, SPEI) 是一种用于评估干旱状况的时间序列指标。它综合考虑了降水量和潜在蒸发量的影响,能够更全面地反映水分平衡状态。通过 GEE 平台可以高效处理大规模地理数据集来实现 SPEI 的计算。 为了完成这一目标,通常需要以下几个关键步骤:定义感兴趣区域、获取气象数据集合(如降水和气温)、利用模型估算实际或潜在蒸散发量、基于这些输入变量构建 SPEI 时间序列以及导出最终结果。 --- #### 定义兴趣区与加载基础影像 假设我们已经有一个名为 `wuhan` 的矢量边界文件作为研究范围,则可以通过如下方式初始化环境: ```javascript // 加载武汉地区特征集合 var wuhan = ee.FeatureCollection("users/yp7454982/wuhan"); Map.centerObject(wuhan, 8); // 设置日期区间 var startDate = ee.Date('2000-01-01'); var endDate = ee.Date('2020-01-01'); // 打印确认对象已成功导入 print("Region of Interest:", wuhan); ``` 上述代码片段展示了如何指定地理位置及其显示级别[^2]。 --- #### 获取必要的气候参数 对于 SPEI 计算而言,至少需取得两个主要因子——月累计降雨量(Precipitation Sum) 和参考作物蒸发速率(Reference Crop Evaporation Rate),即 PET 值。这里推荐采用 CHIRPS 数据源代表前者;而后者可借助 MODIS 或 ERA5-Land 提供的地表辐射通量反演得到。 以下是读取 CHIRPS 雨量资料的一个例子: ```javascript // 导入CHIRPS全球每日降水产品 var chirpsDaily = ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY"); // 过滤时间跨度及空间掩膜操作 var precipitationMonthlySum = chirpsDaily.filterDate(startDate, endDate) .filterBounds(wuhan.geometry()) .reduce(ee.Reducer.sum()); // 将总和转换成每月平均值形式呈现出来 precipitationMonthlySum = precipitationMonthlySum.divide(precipitationMonthlySum.bandNames().size()); ``` 关于 ET₀ 的提取方法,请参阅相关文献[^4]中的具体实践案例描述。 --- #### 实现 SPEI 数学建模过程 一旦拥有了 P 和 ETO 后续就可以着手建立 SPI/SPEI 统计分布拟合流程了。一般情况下我们会选用 Gamma 分布函数去逼近历史记录曲线形状特性从而得出对应的标准分位数映射关系表达式。然而由于篇幅所限在此仅给出伪代码示意框架结构仅供参考学习之用: ```pseudo function computeSpei(precipitationSeries, evapotranspirationSeries){ let combinedDeficit = []; foreach(monthlyDataPoint in range(startYear:endYear)){ deficitValue = monthlyPrecipitation - monthlyEvapotranspirtion; append(deficitValue ,combinedDeficit ); } fittedDistributionParameters = fitGammaDist(combinedDeficit ); zScoresForEachMonth = transformToZScore(fittedDistributionParameters, combinedDeficit ); return zScoresForEachMonth ; } ``` 注意这里的逻辑并未直接体现出来而是留给读者自行补充完善细节部分[^1]。 --- #### 结果输出至本地存储设备 最后一步就是把生成好的栅格图层保存下来以便后续分析应用。下面给出了简单的示例命令行语句帮助理解整个工作流结束后的成果交付环节: ```javascript Export.image.toDrive({ image: speiImage, description: 'spei_wuhan_region', scale: 5000, region: wuhan.geometry(), maxPixels: 1e13 }); ``` 此段脚本负责将处理完毕的地图切片上传到用户的谷歌云端硬盘账户下形成独立项目文件夹管理起来方便调用分享给其他同事协作查看。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

GIS遥感数据处理应用

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

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

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

打赏作者

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

抵扣说明:

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

余额充值