【GRACE学习_2】冯伟老师工具箱主要函数的功能

1.gmt_readgsm
该功能可用于读取GRACE level-2 GSM文件。GSM文件有两种格式，这两种格式都可以被该函数读取。一种是官方的GRACE格式;另一种是ICGEM格式。它们都包含相同的GRACE GSM数据，但格式不同。

2.gmt_replace_degree_1:
由于GRACE所使用的参考系原点是地球的质心，在这个参考系中，一阶斯托克斯系数为零。然而，当我们解释中心的质量变化时，需要考虑一阶项的变化。在这个函数中，用Swenson et al.(2008)估计的一阶项来代替原来的GRACE 一阶项。

3.gmt_replace_C20:
由于原始的C20系数具有较大的不确定性，因此使用卫星激光测距(SLR)的独立估计值来替代GRACE解的原始估计值。在这个函数中，使用了Cheng等人(2013)基于SLR的C20值。

4.gmt_destriping:
该函数提供了几种去条带化方法，以减少GRACE SH系数中的相关误差。

5.gmt_gaussian_filter:
在该函数中，可以对SH系数使用自定义半径的高斯滤波器来抑制SH系数中的高频噪声。

6.gmt_gc2mc:
该函数可将GRACE GSM产品的大地水准面SH系数转换为等效水高质量变化SH系数。

7.gmt_gc2lc:
该函数可以将GRACE GSM产品的大地水准面SH系数转换为载荷变形的SH系数

8.gmt_cs2grid:
在该函数中，谱域的SH系数可以转换为空间域中指定分辨率的网格值，即0.25度、0.5度或1度。

9.gmt_grid2cs:
在该函数中，可以将全局网格化的值转换为指定最大度的SH系数。

10.gmt_cs_error:
基于Wahr et al.(2006)提出的方法，该函数可用于计算GRACE在SH系数中的测量误差。

11.gmt_grid2map:
该函数基于m_map映射工具可视化网格值的全局空间模式

12.gmt_grid2series:
该函数可用于从网格中检索特定区域质量变化的时间序列。除了网格化值之外，边界文件还应作为输入数据给出，其中包括边界点的总数及其经纬度。

13.gmt_harmonic:
该函数可用于对时间序列进行谐波分析。该函数通过最小二乘拟合可以估计出年周期、半年周期和趋势。去除这些估计分量后，可以得到时间序列的年际变化。

相关文献：https://doi.org/10.1007/s12145-018-0368-0