17 - grace数据处理 - 补充 - 地下水储量计算过程分解 -- 陆地水储量变化

0 引言

  最近关于GRACE地下水储量计算的有关过程多被问起，接下来几篇将介绍关于地下水储量计算的有关内容。由水量平衡方程可以将地下水储量的计算过程分解为4个部分，第一部分计算陆地水储量变化、第二部分计算地表水储量变化、第三部分计算冰后回弹（GIA）效应、第四部分计算地下水储量变化。本篇简单介绍下第一部分的内容，主要是GRACE有关数据处理的过程，实际是对前面方法的一个整合 。有需要对GRACE数据处理过程进行了解的可以比较详细的看前面的内容。

1 Grace陆地水储量过程整合

  下面程序包含了GRACE数据解析(以CSR数据为例)、1阶项替换、C20项替换、C30项替换、去相关滤波和高斯滤波等过程，反演了青藏高原地区的陆地质量变化。下面为部分代码，最后有获取方式：

clc;clear;
%% 数据解析
file0