GEE实战第二期：西南八省的极端降水分析

最新推荐文章于 2025-09-07 22:29:54 发布

全域智图

最新推荐文章于 2025-09-07 22:29:54 发布

阅读量580

点赞数 22

CC 4.0 BY-SA版权

文章标签： javascript

本文链接：https://blog.csdn.net/2301_76972702/article/details/150111904

在气候变化背景下，极端降水事件的频率与强度呈现变化趋势。为了评估风险、分析空间分布和时间变化，我们需要一种高效、可复用的方法来处理全球或区域降水数据。

今天分享的这段 Google Earth Engine (GEE) 脚本，结合 NASA GPM IMERG V07 降水数据，可以对你自定义的研究区域（ROI）进行极端降水分析、可视化和导出，整个过程自动化完成。

🧠 用到的知识点

1. GEE 矢量与栅格处理

从资产加载 ROI（FeatureCollection）
按属性筛选（可选）并 dissolve 成一个分析几何

2. 时间序列影像集合处理

使用 filterDate 筛选时间范围
map 和 sum() 聚合半小时/小时降水数据成日累计

3. 指标计算

平均日降水（mean()）
最大日降水（max()）
极端降水天数（阈值判断 + 掩膜 + sum()）
3 日滚动累计最大值（滑动时间窗计算）
最大日降水发生日期（qualityMosaic() 提取）

4. 可视化与图表

Map.addLayer() 显示多指标空间分布
ui.Chart.image.series() 绘制日序列图
ui.Chart.feature.byFeature() 绘制月累计柱状图

5. 数据导出

Export.image.toDrive() 输出多波段 GeoTIFF
Export.table.toDrive() 导出点位日降水 CSV

🔍 核心函数解析

① 日降水构造
将 IMERG 小时级数据聚合成每日总降水量，生成包含降水值和时间戳的影像集合，作为后续分析基础。

② 阈值天数统计
判断每天的降水量是否 ≥ 设定阈值（如 50 mm/day），累加得到极端降水天数影像，用于评估极端事件频率。

③ rolling3(IC)
以三天为滑动窗口，对日降水集合进行累计，取研究期内每个像元的最大累计值，更好反映持续性强降水的影响。

④ 峰值日期提取
利用 qualityMosaic() 在最大降水发生时保留对应的日期波段，实现极端事件时空定位。

⑤ 月累计计算
按月对日降水影像进行累计，输出多波段影像（每波段为一个月），适合季节性与气候特征分析。

📊 运行结果

运行脚本后，你将得到：

地图图层

平均日降水分布
最大日降水分布
≥ 阈值天数分布
3 日累计最大值分布
最大降水发生日期分布

图表输出

ROI 质心点位的日降水时间序列
ROI 平均日降水时间序列
ROI 平均月累计降水柱状图

可选导出（将 DO_EXPORT 设置为 true）

多波段 GeoTIFF（平均值、最大值、阈值天数、3 日累计最大值、峰值日期）
月累计多波段 GeoTIFF
点位日降水 CSV

🛠 使用指南

准备 ROI 资产
- 将你的研究区域（行政区、流域、多边形等）作为 FeatureCollection 上传到 GEE 资产
- 记录资产路径，如 users/username/roi 或 projects/project-id/assets/roi

修改参数区
在脚本顶部设置：

完整代码关注公众号全域智图


// ===================== 参数区（只改这里） =====================
var ASSET_FC_ID   = 'projects/land-111/assets/8province'; // ← 必填：你的 FeatureCollection 资产ID
var USE_PROPERTY  = false;           // 是否按属性筛选资产里的要素
var PROPERTY_NAME = 'NAME';          // 属性名（当 USE_PROPERTY=true 时生效）
var PROPERTY_VAL  = 'Beijing';       // 属性值（当 USE_PROPERTY=true 时生效）

var START_DATE = '2024-01-01';       // 起始（含）
var END_DATE   = '2025-01-01';       // 结束（不含）
var THRESH_MM  = 50;                 // 极端降水阈值（mm/day）
var SCALE_M    = 10000;              // 统计/导出的标称分辨率（米）
var CRS        = 'EPSG:4326';        // 导出坐标系

var DO_EXPORT  = false;              // 是否导出结果到 Google Drive
var EXPORT_FOLDER = 'GEE_exports';   // 导出文件夹名
// ============================================================