用 Python 搭建可复现的地理空间风险流水线

背景

在气象灾害算法项目中，经常需要将 GRIB/NetCDF 格点数据裁剪、重采样并聚合到行政边界。本文记录一条可复现的流水线思路。

核心步骤

1. 统一读取：用 cfgrib / xarray 打开多源文件，对齐坐标命名
2. 空间裁剪：按 bbox 或矢量掩膜裁剪到业务区域
3. 时间对齐：多时效文件合并、填隙与 max 合成
4. 聚合输出：点面关系预计算（.npz）加速乡镇/县级统计

小结

配置驱动 + 中间结果落盘，能显著降低联调成本，也方便业务侧复跑与验收。