1.软件下载地址:https://github.com/fengweiigg/GRACE_Matlab_Toolbox
2.由于目前使用的数据大多都是RL06产品,而测地所冯伟老师的工具箱处理的是RL05产品,所以需要进行一些修改:
(1)gmt_replace_C20函数中RL05替换成RL06
(2)gmt_replace_C21_S21_C22_S22函数中的get_days改为gmt_get_days, RL05替换成RL06
(3)Degree 1和Degree 2 都选用最新的RL06版本
3.控制文件(GMT_Control_File)的书写
每一行代表的含义可参考博主 扎不下村村长的 “GRAMAT控制文件”一文
4.运行(以下每一步都需要一定的时间生成结果,只有弹出success的对话框才说明生成结果成功)
(1)点击 STEP3中的open,打开控制文件,工具箱会自动读取控制文件的信息,点击运行直接出结果
(2)泄露改正
(3)接下来是将生成的球谐分析,将生成的球谐数据文件转换为grid_data.mat
(4)将上一步生成的格网数据转换为时间序列Grid2Series,以亚马逊区域为例:
(5)继续分析,可以得到全球趋势项变化。
5.常见错误
(1)在进行第一步低阶项替换和高斯平滑的时候,如果出现程序运行错误,可以根据进度条卡住的位置以及相关错误提示进行修改。目前,常见的错误是GRACE数据无法识别,主要原因是存储GRACE数据的文件夹的文件名缺失某些关键字(例如GRACE、ICGEM)导致程序无法识别
(2)对于一些没有错误提示的问题,多半是C21_S21、C22_S22、C20_RL06等文件中的说明文件没有删除或者调整正确(调整为如下图所示)