Gaussian 平滑滤波是 GRACE Level-2 球谐系数处理中最常用的基础空间滤波方法。该方法通过对不同阶数的球谐系数施加随阶数递减的平滑权重，削弱高阶噪声向等效水高结果中的传播。其优点是实现简单、参数清晰、结果稳定；缺点是滤波半径较大时会同时压制真实地球物理信号，导致振幅衰减和空间泄漏。
背景
GRACE Level-2 球谐解在高阶项中通常包含较强噪声。当球谐系数被转换为等效水高时，高阶噪声会在空间域中表现为斑块状扰动或条带背景。Gaussian 平滑的基本思路是降低高阶项权重，使反演结果在空间上更加连续。
该方法适合作为基准滤波方案，也适合用于与 DDK、HSAF 等方法进行对照。但需要注意，Gaussian 平滑并不专门针对南北向相关条带，而是对所有方向进行相同平滑，因此属于各向同性空间滤波。
基本原理
Gaussian 滤波假设空间平滑核只与地表两点之间的角距有关，因此其平滑权重仅依赖球谐阶数 l，与次数 m 无关。半径越大，高阶项衰减越强，滤波结果越平滑。
\[\Delta \hat{C}_{lm}=W_l(r)\Delta C_{lm},\quad \Delta \hat{S}_{lm}=W_l(r)\Delta S_{lm}\]
处理流程
- 读取经过预处理的月度 GSM 球谐系数，包括一阶项替换、C20 替换、GIA 改正和基准期扣除后的异常场。
- 设定平滑半径，例如 300 km、500 km 或更大半径。
- 计算每一阶 l 对应的 Gaussian 权重，并将权重作用于所有 m 次球谐系数。
- 将滤波后的球谐系数综合为格网等效水高 EWH 或 TWSA 产品。
- 保留产品标签、半径参数和上游输入信息。
产品标签与配置示例
该页面中的方法不应被统一保存为 filtered。配置文件至少应保留方法名、产品标签、处理域、上游输入和输出对象。
method: Gaussian
tag: GAUSS_300
radius_km: 300
domain: spherical_harmonic
upstream: PREPROCESSED_GSM
output: EWH_grid效果与局限
能够有效降低高阶随机噪声，使全球 TWSA 图像更加平滑。
滤波半径较大时信号损失不可避免，尤其在流域边界、湖泊、冰盖边缘和海岸带容易造成信号外扩散与振幅低估。
页面图件建议
不同半径 Gaussian 权重随阶数变化曲线。
未滤波与 Gaussian 300 km 滤波后的全球 TWSA 对比图。
典型流域滤波前后时间序列对比。
优点—限制双列卡片。
Gaussian 平滑适合作为基础去噪方案和对照基准，但不宜作为精细区域水储量恢复的唯一依据。