GRACE专题--提供一个转换ICGEM网站的gfc文件为mat的代码
果然程序还是得多看,多编,然后融会贯通。
目前主要有CSR、GFZ、JPL三大机构提供GRACE球谐系数,ICGEM网站提供的都是以gfc为后缀的文件。
GSM-2_2002095-2002120_GRAC_UTCSR_BA01_0600.gfc
下面的代码可以读取gfc文件,并转成mat文件。
点击生成的mat文件,比如Cnm系数

其中的C00为1,使用时由于全球质量守恒,因此需要将其置0。此外,C10为0,由于GRACE卫星对地心运动不敏感,需要补充一阶项。另外C20(表中的第三行、第一个)不准确,也需要替换。下面是经过处理的数据。

后期还需要进行滤波处理。希望对大家有所帮助!!