01
处理流程
点击任一处理节点，可展开查看该步骤对应的科学处理约定。
01 获取 获取月尺度 Level-2 GSM Stokes 系数及其产品元数据。
原理
输入应被明确视作 Level-2 重力场解，通常为月尺度 GSM 球谐系数。产品名称、处理中心、版本、最高阶次、时间标识、缺失月份和参考框架必须完整保留。
公式
配置示例
- 产品
- HUST-Grace2024 / CSR RL06 / GFZ RL06 / JPL RL06
- 最高阶次
- 60
- 月度解
- true
质量检查
- 核对产品版本和最高阶次。
- 任何插值前保留原始月份标识。
- 原始输入文件不应被覆盖修改。
02 预处理 完成低阶项替换、GIA 改正、背景场统一和异常基准扣除。
原理
预处理用于将月尺度系数统一为可比较的异常场。典型操作包括一阶项地心改正、C20/C30 SLR 替换、GIA 改正、必要的去相关前置处理，以及扣除明确基准期的平均场。
公式
配置示例
- 一阶项替换
- true
- C20/C30
- SLR
- GIA 模型
- ICE-6G_D
- 基准期
- 2004-01..2009-12
质量检查
- 明确基准期。
- 记录所有替换产品来源。
- 不要把已改正和未改正结果混用同一产品标识。
03 滤波 将去相关、空间平滑、DDK 与 HSAF 作为独立产品路由。
原理
滤波用于削弱南北条带相关噪声和高阶随机噪声，同时尽可能保持真实地表质量信号。Gaussian、Fan、PnMm 去相关、DDK 和 HSAF 应保存为不同产品，而不是互相覆盖。
公式
配置示例
- GAUSS_300
- 300 km Gaussian 平滑
- FAN_300_P3M6
- Fan 滤波结合 P3M6 去相关
- DDK4
- 各向异性正则化产品
- HSAF
- Hankel 谱分析滤波
质量检查
- 每个滤波结果必须有唯一产品标签。
- 记录滤波半径、多项式阶数和起始阶次。
- 比较信号损失、残余条带和泄漏误差。
04 合成 将改正后的球谐系数转换为格网等效水高异常。
原理
球谐合成将改正后的系数转换为等效水高。负荷 Love 数、地球半径、水密度、格网分辨率、掩膜和单位必须在不同产品之间保持一致。
公式
配置示例
- 输出单位
- cmEWH
- 格网
- 1deg 或 0.5deg
- Love 数
- Wahr / Farrell
- 最高阶次
- 60
质量检查
- 确认 cmEWH/mmEWH 单位。
- 对比实验使用同一格网约定。
- 掩膜和海岸缓冲区设置需要版本化。
05 验证 检查图像预览、流域平均、泄漏校正、趋势和独立产品一致性。
原理
验证应同时包含图像质量控制和定量诊断。流域平均时间序列、年振幅、趋势、泄漏校正、Mascon 对比、水文模型和测高资料是不同层级的证据。
公式
配置示例
- 诊断项
- map_preview, basin_mean, annual_amp, trend, rmse
- 参考数据
- CSR_Mascon, GLDAS, Hydroweb
质量检查
- 先检查图像，再进行统计。
- 报告泄漏校正和尺度因子假设。
- 尽量引入至少一种独立产品进行对照。
02
数据产品
数据产品应具有明确命名、版本和可审计路径。
03
方法记录
记录不同滤波方法及其产品路由约束。
软件发布
安装包建议通过 GitHub Releases 分发，网站仅保留文档说明、版本信息和可靠下载入口。
可靠的数据处理必须透明、可测试、可复现。
先规范流程，再讨论结果。