维纳滤波是GRACE数据处理的一种空间滤波方法，它是一种各项同性滤波器，通过设计滤波器，对信号进行线性卷积得到的实际输出信号，使其与期望输出信号满足最小二乘，从而得到维纳滤波函数。通过matlab代码结合网上资源写了计算阶方差的方法，并实现了维纳滤波计算到平滑函数的过程。该程序包包含测试数据、主调函数和相关子函数。

结合GRACE卫星数据和全球陆面数据同化系统GLDAS数据,反演了2004―2009年连续72个月的海河流域地下水储量变化。 在此基础上,结合2004―2009 年海河流域水资源公报的降水量、地下水开采量数据,建立了地下水年开采量与GRACE地下水储量年变化、年降水量的二元回归模型。 利用GRACE卫星数据和GLDAS数据反演的地下水储量年变化与由地下水位观测数据计算出的地下水储量年变化相关性较强,其R2 为0.804;基于GRACE地下水储量年变化数据与年降水量数据,对地下水年开采量的估算结果良好,建立

该资源包含地下水储量计算和泄露误差改正的Matlab代码，有测试数据，可以直接运行得到青藏高原地区的地下水储量变化结果。该程序输入数据为陆地水储量变化、地表水储量变化、gai改正数据，这3个数据用其它程序计算，都在资源中。有关的内容或理论可以查看系列文章【https://blog.csdn.net/weixin_43339605/cat】。如有问题可以留言或私信讨论。

本压缩文件包含Gldas数据处理的Malab代码和测试数据，程序可直接运行，结果输出为文件，需要出图的可以用Gmt进行绘图。本程序简单介绍：由水量平衡方程可以将地下水储量的计算过程分解为以下部分，`第一部分计算陆地水储量变化`、`第二部分计算地表水储量变化`、`第三部分计算冰后回弹改正`、`第四部分计算地下水储量变化`。本篇简单介绍下第二部分的内容，主要是GLDAS水文模型数据的有关处理过程，同样也是对前面几篇博文方法的一个整合或总结 。详细理论和介绍可以参考[https://blog.csdn.net/weixin_43339605/category_12556003.html]系列博文，希望有所帮助，同时遇到问题也可以留言交流。

GRACE数据处理：根据水平衡方程，计算地下水储量变化，要知道陆地质量变化和地表水储量变化，本程序为地下水储量变化计算的一步，用于处理GRACE数据，反演得到陆地质量变化（陆地水储量变化），该程序包含测试数据，可直接运行，如运行出错可更换matlab版本到2019。具体理论及过程可以查看系列文章（https://blog.csdn.net/weixin_43339605/category_12556003.html），如有问题可以留言讨论。

关于Grace模型数据的介绍可以《00 - GRACE卫星及数据简述》，数据由3家机构发布，这里做一个关于数据读取的补充，源码来自(https://github.com/csdms-contrib/slepian_delta)，直接运行slepian_delta中的程序会出现意想不到的错误，下面分享的程序是在此基础上的改进和优化，仅提取了读GRACE数据的相关方法，能够进行3家机构数据的读取和低阶项的替换

尺度因子计算

该资源包含Matlab程序和测试数据，以长江中下游平原为测试区，程序简单符合GRACE数据处理理论，设置好输入就可以输出得到需要的结果。可以参照博文《02 - GRACE数据处理步骤简叙》进行理解。如有问题可以留言交流。

。。。

GRACE卫星监测地表质量迁移