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

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系数。

这个函数的作用是将GRACE卫星重力场数据中的大地水准面（geoid）的球谐系数转换为等效水位高（mass changes in equivalent water height）的
球谐系数。具体来说，该函数根据Wahr等人1998年的论文中的方程式（14）中的系数和与度数相关的因子，将球谐系数乘以一个系数，以将负荷
（在kg/m2中）转换为相当于水位高度的质量系数（以m为单位）。

该函数首先会检查输入的系数是不是CS格式（经纬格点上的球谐系数），如果是，则提取最大阶数并继续运算。如果输入的系数是SC格式
（普通球谐系数），则转换为CS格式，并提取最大阶数。然后，该函数计算密度因子dens，然后使用来自Han和Wahr的Love数，
以及从Wahr等人1998年论文方程式（13）中计算出来的因子factor，将每行系数与factor相乘。最后，输出转换后的等效水位高的球谐系数，
格式为CS。

如果输入的文件格式不正确，则函数会输出错误信息

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

这段代码是一个 MATLAB 函数，它的功能是将一个网格数据转换为一个一维时间序列，用于后续分析和可视化。

函数的输入参数包括：

grid：网格数据，可以是二维的经纬度网格，也可以是三维的经纬度和时间网格。
dir_msk：边界线数据的文件路径。
type：边界线数据的类型，可以是 line 或 mask。
bound：用于边界线数据的范围，当 type 为 mask 时使用。
函数的输出参数是一个一维数组 plot_region，其中包含了网格数据的时间序列。

该函数的主要功能包括：

判断输入的网格数据是二维的还是三维的。
根据输入的边界线数据，创建一个网格掩模，将网格中边界线以外的部分置为 0，边界线以内的部分置为 1。
遍历网格数据，将不在掩模范围内的网格点的值设为 NaN。
将经过掩模处理的网格数据转换为一维时间序列 plot_region，用于后续分析和可视化。
需要注意的是，这段代码中的很多数值都是针对特定的网格数据的，比如经纬度范围、网格分辨率等等，所以在使用时需要根据具体情况进行修改

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

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