METHOD NOTE

DDK 正则化滤波

基于误差协方差与正则化约束的各向异性滤波产品。

正则化协方差约束各向异性DDK1–DDK8
处理域球谐系数域
典型输入预处理后的 GRACE 球谐系数
典型输出DDK1–DDK8 滤波 EWH/TWSA 产品
页面摘要

DDK 滤波是一类基于误差协方差、信号协方差和 Tikhonov 正则化思想构建的各向异性滤波方法。与 Gaussian、Fan 等经验平滑滤波不同，DDK 方法显式引入误差统计信息，能够在去条带与信号保持之间取得相对稳定的折中。

背景

GRACE 球谐解中的误差具有明显的各向异性和相关性。传统平滑滤波主要依赖经验半径，难以充分利用误差协方差信息。DDK 滤波通过正则化框架构造滤波算子，使滤波强度与误差统计和先验信号约束相关，因此在 GRACE 数据处理中被广泛采用。

基本原理

DDK 滤波可理解为对球谐系数最小二乘解的正则化改造。其核心思想是在保留观测解的同时，引入先验约束，抑制不稳定或噪声主导的系数分量。不同参数组合形成 DDK1 至 DDK8 等不同强度产品。

\[\hat{x}=(N+\alpha R)^{-1}Nx\]

处理流程

  1. 准备预处理后的 GRACE 球谐系数。
  2. 选择 DDK 类型，例如 DDK1、DDK2、DDK3、DDK4 等。
  3. 调用对应 DDK 滤波矩阵或产品算子。
  4. 将滤波后的球谐系数转换为 EWH/TWSA。
  5. 在产品标签中保留 DDK 编号，不能只写 DDK。

产品标签与配置示例

该页面中的方法不应被统一保存为 filtered。配置文件至少应保留方法名、产品标签、处理域、上游输入和输出对象。

method: DDK
tag: DDK4
domain: spherical_harmonic
regularization: anisotropic
upstream: PREPROCESSED_GSM
output: EWH_grid

效果与局限

效果

DDK1 和 DDK2 滤波强度较大，条带抑制效果明显，但可能导致真实信号损失。

局限

DDK5 至 DDK8 滤波较弱，可能保留较多残余条带；DDK3 和 DDK4 通常在去条带和信号保持之间更均衡。

页面图件建议

DDK1–DDK8 全球 TWSA 对比图。

DDK4 与 Gaussian/Fan/HSAF 的频谱对比图。

DDK 强度梯度说明条。

推荐使用场景卡片。

使用判断

DDK4 适合作为稳健对照产品，但仍需注意其正则化先验可能造成局部信号衰减。